Physics:Dean number

From HandWiki
Short description: Dimensionless group in fluid mechanics

The Dean number (De) is a dimensionless group in fluid mechanics, which occurs in the study of flow in curved pipes and channels. It is named after the British scientist W. R. Dean, who was the first to provide a theoretical solution of the fluid motion through curved pipes for laminar flow by using a perturbation procedure from a Poiseuille flow in a straight pipe to a flow in a pipe with very small curvature.[1][2]

Physical Context

Schematic of a pair of Dean vortices that form in curved pipes.

If a fluid is moving along a straight pipe that after some point becomes curved, the centripetal forces at the bend will cause the fluid particles to change their main direction of motion. There will be an adverse pressure gradient generated from the curvature with an increase in pressure, therefore a decrease in velocity close to the convex wall, and the contrary will occur towards the outer side of the pipe. This gives rise to a secondary motion superposed on the primary flow, with the fluid in the centre of the pipe being swept towards the outer side of the bend and the fluid near the pipe wall will return towards the inside of the bend. This secondary motion is expected to appear as a pair of counter-rotating cells, which are called Dean vortices.

Definition

The Dean number is typically denoted by De (or Dn). For a flow in a pipe or tube it is defined as:

[math]\displaystyle{ \mathit{De} = \frac{\sqrt{\frac{1}{2}\,(\text{inertial forces})(\text{centripetal forces})}}{\text{viscous forces}} = \frac{\sqrt{\frac{1}{2}\,(\rho\,D^2\,R_c\,\frac{v^2}{D}) (\rho\,D^2\,R_c\,\frac{v^2}{R_c})}}{\mu \frac{v}{D} D\,R_c} = \frac{\rho\,D\,v}{\mu} \sqrt{ \frac{D}{2\,R_c} } = \textit{Re} \, \sqrt{ \frac{D}{2\,R_c} } }[/math]

where

  • [math]\displaystyle{ \rho }[/math] is the density of the fluid
  • [math]\displaystyle{ \mu }[/math] is the dynamic viscosity
  • [math]\displaystyle{ v }[/math] is the axial velocity scale
  • [math]\displaystyle{ D }[/math] is the diameter (for non-circular geometry, an equivalent diameter is used; see Reynolds number)
  • [math]\displaystyle{ R_c }[/math] is the radius of curvature of the path of the channel.
  • [math]\displaystyle{ \textit{Re} }[/math] is the Reynolds number.

The Dean number is therefore the product of the Reynolds number (based on axial flow [math]\displaystyle{ v }[/math] through a pipe of diameter [math]\displaystyle{ D }[/math]) and the square root of the curvature ratio.

Turbulence transition

The flow is completely unidirectional for low Dean numbers (De < 40~60). As the Dean number increases between 40~60 to 64~75, some wavy perturbations can be observed in the cross-section, which evidences some secondary flow. At higher Dean numbers than that (De > 64~75) the pair of Dean vortices becomes stable, indicating a primary dynamic instability. A secondary instability appears for De > 75~200, where the vortices present undulations, twisting, and eventually merging and pair splitting. Fully turbulent flow forms for De > 400.[3] Transition from laminar to turbulent flow has also been examined in a number of studies, even though no universal solution exists since the parameter is highly dependent on the curvature ratio.[4] Somewhat unexpectedly, laminar flow can be maintained for larger Reynolds numbers (even by a factor of two for the highest curvature ratios studied) than for straight pipes, even though curvature is known to cause instability.[5]

The Dean equations

The Dean number appears in the so-called Dean equations.[6] These are an approximation to the full Navier–Stokes equations for the steady axially uniform flow of a Newtonian fluid in a toroidal pipe, obtained by retaining just the leading order curvature effects (i.e. the leading-order equations for [math]\displaystyle{ a/r \ll 1 }[/math]).

We use orthogonal coordinates [math]\displaystyle{ (x,y,z) }[/math] with corresponding unit vectors [math]\displaystyle{ (\hat{\boldsymbol{x}},\hat{\boldsymbol{y}},\hat{\boldsymbol{z}}) }[/math] aligned with the centre-line of the pipe at each point. The axial direction is [math]\displaystyle{ \hat{\boldsymbol{z}} }[/math], with [math]\displaystyle{ \hat{\boldsymbol{x}} }[/math] being the normal in the plane of the centre-line, and [math]\displaystyle{ \hat{\boldsymbol{y}} }[/math] the binormal. For an axial flow driven by a pressure gradient [math]\displaystyle{ G }[/math], the axial velocity [math]\displaystyle{ u_z }[/math] is scaled with [math]\displaystyle{ U=Ga^2/\mu }[/math]. The cross-stream velocities [math]\displaystyle{ u_x, u_y }[/math] are scaled with [math]\displaystyle{ (a/R)^{1/2} U }[/math], and cross-stream pressures with [math]\displaystyle{ \rho a U^2/L }[/math]. Lengths are scaled with the tube radius [math]\displaystyle{ a }[/math].

In terms of these non-dimensional variables and coordinates, the Dean equations are then

[math]\displaystyle{ De \left( \frac{\mathrm{D} u_x}{\mathrm{D} t} - u_z^2 \right) = -De \frac{\partial p}{\partial x} + \nabla^2 u_x }[/math]
[math]\displaystyle{ De \frac{\mathrm{D} u_y}{\mathrm{D} t} = -De\frac{\partial p}{\partial y} + \nabla^2 u_y }[/math]
[math]\displaystyle{ De \frac{\mathrm{D} u_z}{\mathrm{D} t} = 1 + \nabla^2 u_z }[/math]
[math]\displaystyle{ \frac{\partial u_x}{\partial x} + \frac{\partial u_y}{\partial y} = 0 }[/math]

where

[math]\displaystyle{ \frac{\mathrm{D}}{\mathrm{D} t} = u_x \frac{\partial}{\partial x} + u_y \frac{\partial}{\partial y} }[/math]

is the convective derivative.

The Dean number De is the only parameter left in the system, and encapsulates the leading order curvature effects. Higher-order approximations will involve additional parameters.

For weak curvature effects (small De), the Dean equations can be solved as a series expansion in De. The first correction to the leading-order axial Poiseuille flow is a pair of vortices in the cross-section carrying flow from the inside to the outside of the bend across the centre and back around the edges. This solution is stable up to a critical Dean number [math]\displaystyle{ De_c \approx 956 }[/math].[7] For larger De, there are multiple solutions, many of which are unstable.

References

  1. Dean, W. R. (1927). "Note on the motion of fluid in a curved pipe". Phil. Mag. 4 (20): 208–223. doi:10.1080/14786440708564324. http://www.informaworld.com/openurl?genre=article&issn=1941-5982&volume=4&issue=20&spage=208. 
  2. Dean, W. R. (1928). "The streamline motion of fluid in a curved pipe". Phil. Mag.. Series 7 5 (30): 673–695. doi:10.1080/14786440408564513. http://www.informaworld.com/openurl?genre=article&issn=1941-5982&volume=5&issue=30&spage=673. 
  3. Ligrani, Phillip M. "A Study of Dean Vortex Development and Structure in a Curved Rectangular Channel With Aspect Ratio of 40 at Dean Numbers up to 430", U.S. Army Research Laboratory (Contractor Report ARL-CR-l44) and Lewis Research Center (NASA Contractor Report 4607), July 1994. Retrieved on 11 July 2017.
  4. Kalpakli, Athanasia (2012). Experimental study of turbulent flows through pipe bends (Thesis). Stockholm, Sweden: Royal Institute of Technology KTH Mechanics. pp. 461–512.
  5. Taylor, G. I. (1929). "The criterion for turbulence in curved pipes". Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences 124 (794): 243–249. doi:10.1098/rspa.1929.0111. Bibcode1929RSPSA.124..243T. 
  6. Mestel, J. Flow in curved pipes: The Dean equations, Lecture Handout for Course M4A33, Imperial College.
  7. Dennis, C. R.; Ng, M. (1982). "Dual solutions for steady laminar-flow through a curved tube". Q. J. Mech. Appl. Math. 35 (3): 305. doi:10.1093/qjmam/35.3.305. 

Further reading