Research Papers: Multiphase Flows

Mathematical Model for Viscoplastic Fluid Hammer

[+] Author and Article Information
Gabriel M. Oliveira, Admilson T. Franco

Research Center for Rheology
and Non-Newtonian Fluids (CERNN),
Post-Graduate Program in Mechanical
and Materials Engineering (PPGEM),
Federal University of
Technology - Paraná (UTFPR),
Av. Sete de Setembro,
Curitiba, PR 3165, Brazil

Cezar O. R. Negrão

Research Center for Rheology
and Non-Newtonian Fluids (CERNN),
Post-Graduate Program in Mechanical
and Materials Engineering (PPGEM),
Federal University of
Technology - Paraná (UTFPR),
Av. Sete de Setembro,
Curitiba, PR 3165, Brazil
e-mail: negrao@utfpr.edu.br

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received January 22, 2015; final manuscript received June 29, 2015; published online August 10, 2015. Assoc. Editor: Satoshi Watanabe.

J. Fluids Eng 138(1), 011301 (Aug 10, 2015) (8 pages) Paper No: FE-15-1054; doi: 10.1115/1.4031001 History: Received January 22, 2015

The current work presents a mathematical model to simulate “viscoplastic fluid hammer”-overpressure caused by sudden viscoplastic fluid deceleration in pipelines. The flow is considered one-dimensional, isothermal, laminar, and weakly compressible and the fluid is assumed to behave as a Bingham plastic. The model is based on the mass and momentum balance equations and solved by the method of characteristics (MOC). The results show that the overpressures taking place in viscoplastic fluids are smaller than those occurring in Newtonian fluids and also that two pressure gradients-one negative and one positive-are possibly noted after pressure stabilization. The pressure stabilizes nonuniformly on the pipeline because viscoplastic fluids present yield stresses. Overpressure magnitudes depend not only on the ratio of pressure wave inertia to viscous effect but also on the Bingham number. The pipeline designer should take into account the viscoplastic fluid behavior reported in this paper when engineering a new pipeline system.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.


Streeter, V. L. , and Wylie, E. B. , 1974, “Waterhammer and Surge Control,” Annu. Rev. Fluid Mech., 6, pp. 57–73. [CrossRef]
Mitsoulis, E. , 2007, “Flows of Viscoplastic Materials: Models and Computations,” Rheol. Rev., pp. 135–178.
Ghidaoui, M. S. , Zhao, M. , McInnis, D. A. , and Axworthy, D. H. , 2005, “A Review of Water Hammer Theory and Practice,” ASME Appl. Mech. Rev., 58(1), pp. 49–76. [CrossRef]
Bergant, A. , Simpson, A. R. , and Tijsseling, A. S. , 2006, “Water Hammer With Column Separation: A Historical Review,” J. Fluids Struct., 22(2), pp. 135–171. [CrossRef]
Brunone, B. , Karney, B. W. , Mecarelli, M. , and Ferrante, M. , 2000, “Velocity Profiles and Unsteady Pipe Friction in Transient Flow,” ASCE J. Water Resour. Plann. Manage., 126(4), pp. 236–244. [CrossRef]
Brunone, B. , Ferrante, M. , and Cacciamani, M. , 2004, “Decay of Pressure and Energy Dissipation in Laminar Transient Flow,” ASME J. Fluid. Eng., 126(6), pp. 928–934. [CrossRef]
Wahba, E. M. , 2013, “Non-Newtonian Fluid Hammer in Elastic Circular Pipes: Shear-Thinning and Shear-Thickening Effects,” J. Non-Newtonian Fluid Mech., 198, pp. 24–30. [CrossRef]
Bird, R. B. , Armstrong, R. C. , and Hassager, O. , 1987, Dynamics of Polymeric Liquids, 2nd ed., Vol. 1, Wiley, New York, Chap. 2.
Oliveira, G. M. , Negrão, C. O. R. , and Franco, A. T. , 2012, “Pressure Transmission in Bingham Fluids Compressed Within a Closed Pipe,” J. Non-Newtonian Fluid Mech., 169–170, pp. 121–125. [CrossRef]
Sestak, J. , Cawkwell, M. G. , Charles, M. E. , and Houskas, M. , 1987, “Start-Up of Gelled Crude Oil Pipelines,” J. Pipelines, 6(1), pp. 15–24.
Cawkwell, M. G. , and Charles, M. E. , 1987, “An Improved Model for Start-Up of Pipelines Containing Gelled Crude Oil,” J. Pipelines, 7(1), pp. 41–52.
Chang, C. , Rønningsen, H. P. , and Nguyen, Q. D. , 1999, “Isothermal Start-Up of Pipeline Transporting Waxy Crude Oil,” J. Non-Newtonian Fluid Mech., 87(2–3), pp. 127–154. [CrossRef]
Davidson, M. R. , Nguyen, Q. D. , Chang, C. , and Rønningsten, H. P. , 2004, “A Model for Restart of a Pipeline With Compressible Gelled Waxy Crude Oil,” J. Non-Newtonian Fluid Mech., 123(2–3), pp. 269–280. [CrossRef]
Oliveira, G. M. , Rocha, L. L. V. , Franco, A. T. , and Negrão, C. O. R. , 2010, “Numerical Simulation of the Start-Up of Bingham Fluid Flows in Pipelines,” J. Non-Newtonian Fluid Mech., 165(19–20), pp. 1114–1128. [CrossRef]
Vinay, G. , Wachs, A. , and Agassant, J. F. , 2006, “Numerical Simulation of Weakly Compressible Bingham Flows: The Restart of Pipeline Flows of Waxy Crude Oils,” J. Non-Newtonian Fluid Mech., 136(2–3), pp. 93–105. [CrossRef]
Vinay, G. , Wachs, A. , and Frigaard, I. , 2007, “Start-Up Transients and Efficient Computation of Isothermal Waxy Crude Oil Flows,” J. Non-Newtonian Fluid Mech., 143(2–3), pp. 141–156. [CrossRef]
Wachs, A. , Vinay, G. , and Frigaard, I. , 2009, “A 1.5d Numerical Model for the Start Up of Weakly Compressible Flow of a Viscoplastic and Thixotropic Fluid in Pipelines,” J. Non-Newtonian Fluid Mech., 159(1–3), pp. 81–94. [CrossRef]
Negrão, C. O. R. , Franco, A. T. , and Rocha, L. L. V. , 2011, “A Weakly Compressible Flow Model for the Restart of Thixotropic Drilling Fluids,” J. Non-Newtonian Fluid Mech., 166(23–24), pp. 1369–1381. [CrossRef]
Anderson, J. D. , 1990, Modern Compressible Flow: With Historical Perspective, 2nd ed., McGraw-Hill, New York, Chaps. 1 and 3.
Swamee, P. K. , and Aggarwal, N. , 2011, “Explicit Equations for Laminar Flow of Bingham Plastic Fluids,” J. Pet. Sci. Eng., 76(3–4), pp. 178–184. [CrossRef]
Melrose, J. C. , Savins, J. G. , Foster, W. R. , and Parisha, E. R. , 1958, “A Practical Utilization of the Theory of Bingham Plastic Flow in Stationary Pipes and Annuli,” Pet. Trans. AIME, 213, pp. 316–324.
Wylie, E. B. , Streeter, V. L. , and Suo, L. , 1993, Fluid Transients in Systems, Prentice Hall, NJ, Chap. 3.
Wahba, E. M. , 2008, “Modeling the Attenuation of Laminar Fluid Transients in Piping Systems,” Appl. Math. Modell., 32(12), pp. 2863–2871. [CrossRef]
Oliveira, G. M. , Franco, A. T. , Negrão, C. O. R. , Martins, A. L. , and Silva, R. A. , 2013, “Modeling and Validation of Pressure Propagation in Drilling Fluids Pumped Into a Closed Well,” J. Pet. Sci. Eng., 103, pp. 61–71. [CrossRef]
Holmboe, E. L. , and Rouleau, W. T. , 1967, “The Effect of Viscous Shear on Transients in Liquid Lines,” ASME J. Basic Eng., 89(1), pp. 174–180. [CrossRef]
Sobey, R. J. , 2004, “Analytical Solutions for Unsteady Pipe Flow,” J. Hydroinf., 6, pp. 187–207.


Grahic Jump Location
Fig. 1

Illustration of the physical problem

Grahic Jump Location
Fig. 2

Spatial and temporal grid

Grahic Jump Location
Fig. 3

Effect of mesh size for λ = 1.0 and BR = 0.5: (a) time evolution of pressure at Z = 0.6 and (b) axial pressure distribution at T = 0.4

Grahic Jump Location
Fig. 4

Comparisons of the numerical results with experimental data [25] and with an analytical solution. Time evolution of pressure: (a) at the valve position and (b) at pipe midpoint position.

Grahic Jump Location
Fig. 7

Effect of the Bingham number on the time evolution of (a) pressure and (b) velocity at pipe midpoint position for λ = 1.0

Grahic Jump Location
Fig. 8

Time evolution of the axial pressure distribution for λ = 0.1: (a) BR = 0 and (b) BR = 0.5

Grahic Jump Location
Fig. 9

Effect of the Bingham number on the time evolution of pressure at the valve position for (a) λ = 0.1, (b) λ = 1.0, (c) λ = 5.0, and (d) λ = 10.0

Grahic Jump Location
Fig. 10

Effect of λ on the final pressure distribution for (a) BR = 0.25 and (b) BR = 0.75

Grahic Jump Location
Fig. 11

Final pressure distribution as a function of the Bingham number for λ = 1.0

Grahic Jump Location
Fig. 5

Time evolution of the axial pressure distribution for λ = 1.0: (a) and (b) BR = 0 and (c) and (d) BR = 0.5

Grahic Jump Location
Fig. 6

Time evolution of the axial velocity distribution for λ = 1.0: (a) BR = 0 and (b) BR = 0.5



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In