0
Research Papers: Fundamental Issues and Canonical Flows

Developing Region Solution for High Reynolds Number Laminar Flows of Pseudoplastic and Dilatant Fluids in Circular Ducts

[+] Author and Article Information
Massimo Capobianchi

Professor
Mem. ASME
Department of Mechanical Engineering,
Gonzaga University,
502 E. Boone Avenue,
Spokane, WA 99258-0026
e-mail: capobianchi@gonzaga.edu

Patrick McGah

Mem. ASME
Department of Mechanical Engineering,
University of Washington,
Box 352600,
Seattle, WA 98195
e-mail: pmcgah@u.washington.edu

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received August 8, 2015; final manuscript received November 7, 2016; published online February 6, 2017. Assoc. Editor: Kausik Sarkar.

J. Fluids Eng 139(4), 041202 (Feb 06, 2017) (11 pages) Paper No: FE-15-1547; doi: 10.1115/1.4035242 History: Received August 08, 2015; Revised November 07, 2016

This article reports the results of a numerical computation of the length and total pressure drop in the entrance region of a circular tube with laminar flows of pseudoplastic and dilatant fluids at high Reynolds numbers (i.e., approximately 400 or higher). The analysis utilizes equations for the apparent viscosity that span the entire shear rate regime, from the zero to the infinite shear rate Newtonian regions, including the power law and the two transition regions. Solutions are thus reported for all shear rates that may exist in the flow field, and a shear rate parameter is identified that quantifies the shear rate region where the system is operating. The entrance lengths and total pressure drops were found to be bound by the Newtonian and power law values, the former being approached when the system is operating in either the zero or the infinite shear rate Newtonian regions. The latter are approached when the shear rates are predominantly in the power law region but only if, in addition, the zero and infinite shear rate Newtonian viscosities differ sufficiently, by approximately four orders of magnitude or more. For all other cases, namely, when more modest differences in the limiting Newtonian viscosities exist, or when the system is operating in the low- or high-shear rate transition regions, then intermediate results are obtained. Entrance length and total pressure drop values are provided in both graphical form, and in tabular and correlation equation form, for convenient access.

FIGURES IN THIS ARTICLE
<>
Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.

References

Skelland, A. H. P. , 1967, Non-Newtonian Flow and Heat Transfer, Wiley, New York, pp. 5–7.
Dunleavy, J. E., Jr. , and Middleman, S. , 1966, “ Correlation of Shear Behavior of Solutions of Polyisobutylene,” Trans. Soc. Rheol., 10(1), pp. 157–168. [CrossRef]
Capobianchi, M. , and Irvine, T. F., Jr ., 1992, “ Predictions of Pressure Drop and Heat Transfer in Concentric Annular Ducts With Modified Power Law Fluids,” Heat Mass Transfer, 27(4), pp. 209–215.
Capobianchi, M. , 2008, “ Pressure Drop Predictions for Laminar Flows of Extended Modified Power Law Fluids in Rectangular Ducts,” Int. J. Heat Mass Transfer, 51(5–6), pp. 1393–1401. [CrossRef]
Cross, M. M. , 1965, “ Rheology of Non-Newtonian Fluids: A New Flow Equation for Pseudoplastic Systems,” J. Colloid Sci., 20(5), pp. 417–437. [CrossRef]
Shah, R. K. , and London, A. L. , 1978, Laminar Flow Forced Convection in Ducts, A Source Book for Compact Heat Exchanger Analytical Data (Advances in Heat Transfer (Supplement 1)), T. F. Irvine, Jr. and J. P. Hartnett , eds., Academic Press, New York, pp. 85–99.
Collins, M. , and Schowalter, W. R. , 1963, “ Behavior of Non-Newtonian Fluids in the Entry Region of a Pipe,” AIChE J., 9(6), pp. 804–809. [CrossRef]
Fargie, D. , and Martin, B. W. , 1971, “ Developing Laminar Flow in a Pipe of Circular Cross-Section,” Proc. R. Soc. London A, 321(1547), pp. 461–476. [CrossRef]
Bogue, D. C. , 1959, “ Entrance Effects and Prediction of Turbulence in Non-Newtonian Flow,” Ind. Eng. Chem., 51(7), pp. 874–878. [CrossRef]
Chebbi, R. , 2002, “ Laminar Flow of Power-Law Fluids in the Entrance Region of a Pipe,” Chem. Eng. Sci., 57(21), pp. 4435–4443. [CrossRef]
Mashelkar, R. A. , 1974, “ Hydrodynamic Entrance-Region Flow of Pseudoplastic Fluids,” Proc. Inst. Mech. Eng., 188(1974), pp. 683–689. [CrossRef]
Matras, Z. , and Nowak, Z. , 1983, “ Laminar Entry Length Problem for Power Law Fluids,” Acta Mech., 48(1), pp. 81–90. [CrossRef]
Sparrow, E. M. , Lin, S. H. , and Lundgren, T. S. , 1964, “ Flow Development in the Hydrodynamic Entrance Region of Tubes and Ducts,” Phys. Fluids, 7(3), pp. 338–347. [CrossRef]
Langhaar, H. L. , 1942, “ Steady Flow in the Transition Length of a Straight Tube,” ASME J. Appl. Mech., 9(2), pp. A55–A58.
Durst, F. , Ray, S. , Ünsal, B. , and Bayoumi, O. A. , 2005, “ The Development Lengths of Laminar Pipe and Channel Flows,” ASME J. Fluids Eng., 127(6), pp. 1154–1160. [CrossRef]
Christiansen, E. B. , and Lemmon, H. E. , 1965, “ Entrance Region Flow,” AIChE J., 11(6), pp. 995–999. [CrossRef]
Vrentas, J. S. , Duda, J. L. , and Bargeron, K. G. , 1966, “ Effect of Axial Diffusion of Vorticity on Flow Development in Circular Conduits: Part I. Numerical Solutions,” AIChE J., 12(5), pp. 837–844. [CrossRef]
Hornbeck, R. W. , 1964, “ Laminar Flow in the Entrance Region of a Pipe,” Appl. Sci. Res., Sect. A, 13(1), pp. 224–232. [CrossRef]
Poole, R. J. , and Ridley, B. S. , 2007, “ Development-Length Requirements for Fully Developed Laminar Pipe Flow of Inelastic Non-Newtonian Liquids,” ASME J. Fluids Eng., 129(10), pp. 1281–1287. [CrossRef]
Mehrotra, A. K. , and Patience, G. S. , 1990, “ Unified Entry Length for Newtonian and Power-Law Fluids in Laminar Pipe Flow,” Can. J. Chem. Eng., 68(4), pp. 529–533. [CrossRef]
Soto, R. J. , and Shah, V . L. , 1976, “ Entrance Flow of a Yield-Power Law Fluid,” Appl. Sci. Res., 32(1), pp. 73–85. [CrossRef]
Bejan, A. , 1984, Convection Heat Transfer, Wiley, New York, pp. 28–39.
Brewster, R. A. , and Irvine, T. F., Jr ., 1987, “ Similitude Consideration in Laminar Flow of Modified Power Law Fluids in Circular Ducts,” Heat Mass Transfer, 21(2–3), pp. 83–86.
Patankar, S. V. , 1980, Numerical Heat Transfer and Fluid Flow, Hemisphere Publishing Company, New York, pp. 11–153.
Brewster, R. A. , 2013, “ Pressure Drop Predictions for Laminar Fully-Developed Flows of Purely-Viscous Non-Newtonian Fluids in Circular Ducts,” ASME J. Fluids Eng., 135(10), p. 101106. [CrossRef]

Figures

Grahic Jump Location
Fig. 1

Schematic flow curves for pseudoplastic and dilatant fluids and definition of the five flow regions: region I—zero shear rate Newtonian region; region II—low shear rate transition region; region III—power law region; region IV—high-shear rate transition region; and region V—infinite shear rate Newtonian region

Grahic Jump Location
Fig. 2

Schematic of the entrance region problem with definition of coordinate axes andcorresponding velocity components. HBL—hydrodynamic boundary layer; HER—hydrodynamic entrance region; and HFDR—hydrodynamic fully developed region.

Grahic Jump Location
Fig. 3

u+ velocity profiles for n = 0.50, Log10(β*) = 2.00, and Log10(R*) = −4.00 at x+= 0.0001, 0.005, 0.02284, 0.05534, 0.075, and 0.1228 (from left to right along the horizontal axis). The profile drawn in heavy-weight line, at x+= 0.1228, is the x+ location where the 99% centerline velocity criterion is met.

Grahic Jump Location
Fig. 4

u+ velocity profiles for n = 1.00 (Newtonian fluid) at x+= 0.0001, 0.005, 0.02284, and 0.05534 (from left to right along the horizontal axis). The profile drawn in heavy-weight line, at x+= 0.05534, is the x+ location where the 99% centerline velocity criterion is met.

Grahic Jump Location
Fig. 5

u+ velocity profiles for n = 1.5, Log10(β*) = −2.00, and Log10(R*) = 4.00 at x+= 0.0001, 0.005, and 0.02284 (from left to right along the horizontal axis). The profile drawn in heavy-weight line, at x+= 0.02284, is the x+ location where the 99% centerline velocity criterion is met.

Grahic Jump Location
Fig. 6

xe+ results using the 99% centerline velocity criterion. From inner to outer curve in each family, n = 0.90, 0.80, 0.70 0.60, and 0.50 for pseudoplastic fluids, and n = 1.10, 1.20, 1.30 1.40, and 1.50 for dilatant fluids. Also shown are horizontal linesat the Newtonian and power law asymptotes, the latter computed at | Log10(β*)|=10 with | Log10(R*)|=20.

Grahic Jump Location
Fig. 7

xe+ results using the 101% f·ReM criterion. From inner to outer curve in each family, n = 0.90, 0.80, 0.70 0.60, and 0.50 for pseudoplastic fluids, and n = 1.10, 1.20, 1.30 1.40, and 1.50 for dilatant fluids. Also shown are horizontal lines at the Newtonian and power law asymptotes, the latter computed at | Log10(β*)|=10 with | Log10(R*)|=20. Numerical values for the points indicated by “” are tabulated in Table 1.

Grahic Jump Location
Fig. 8

fx·ReM¯ results using the 101% f·ReM criterion. From inner to outer curve in each family, n = 0.90, 0.80, 0.70 0.60, and 0.50 for pseudoplastic fluids, and n = 1.10, 1.20, 1.30 1.40, and 1.50 for dilatant fluids. Also shown are horizontal lines at the Newtonian and power law asymptotes, the latter computed at | Log10(β*)|=10 with | Log10(R*)|=20.

Grahic Jump Location
Fig. 9

f·ReM (i.e., HFDR) results. From inner to outer curve in each family, n = 0.90, 0.80, 0.70 0.60, and 0.50 for pseudoplastic fluids, and n = 1.10, 1.20, 1.30 1.40, and 1.50 for dilatant fluids. Also shown are horizontal lines at the Newtonian and power law asymptotes, the latter computed at | Log10(β*)|=10 with | Log10(R*)|=20.

Tables

Errata

Discussions

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