Non-dimensionalization and scaling of the Navier–Stokes equations

From HandWiki

In fluid mechanics, non-dimensionalization of the Navier–Stokes equations is the conversion of the Navier–Stokes equation to a nondimensional form. This technique can ease the analysis of the problem at hand, and reduce the number of free parameters. Small or large sizes of certain dimensionless parameters indicate the importance of certain terms in the equations for the studied flow. This may provide possibilities to neglect terms in (certain areas of) the considered flow. Further, non-dimensionalized Navier–Stokes equations can be beneficial if one is posed with similar physical situations – that is problems where the only changes are those of the basic dimensions of the system.

Scaling of Navier–Stokes equation refers to the process of selecting the proper spatial scales – for a certain type of flow – to be used in the non-dimensionalization of the equation. Since the resulting equations need to be dimensionless, a suitable combination of parameters and constants of the equations and flow (domain) characteristics have to be found. As a result of this combination, the number of parameters to be analyzed is reduced and the results may be obtained in terms of the scaled variables.

Need for non-dimensionalization and scaling

In addition to reducing the number of parameters, non-dimensionalized equation helps to gain a greater insight into the relative size of various terms present in the equation.[1][2] Following appropriate selecting of scales for the non-dimensionalization process, this leads to identification of small terms in the equation. Neglecting the smaller terms against the bigger ones allows for the simplification of the situation. For the case of flow without heat transfer, the non-dimensionalized Navier–Stokes equation depend only on the Reynolds Number and hence all physical realizations of the related experiment will have the same value of non-dimensionalized variables for the same Reynolds Number.[3]

Scaling helps provide better understanding of the physical situation, with the variation in dimensions of the parameters involved in the equation. This allows for experiments to be conducted on smaller scale prototypes provided that any physical effects which are not included in the non-dimensionalized equation are unimportant.

Non-dimensionalized Navier–Stokes equation

The incompressible Navier–Stokes momentum equation is written as:

[math]\displaystyle{ \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} = -\frac 1 \rho \nabla p + \nu \nabla^2 \mathbf{u} + \mathbf{g}. }[/math] [4][5]

where ρ is the density, p is the pressure, ν is the kinematic viscosity, u is the flow velocity, and g is the body acceleration field.

The above equation can be non-dimensionalized through selection of appropriate scales as follows:

Scale dimensionless variable
Length L [math]\displaystyle{ \mathbf r^*\ = \frac{\mathbf r}{L} }[/math] and [math]\displaystyle{ \nabla^*\ = L\nabla }[/math]
Flow velocity U [math]\displaystyle{ \mathbf u^*\ = \frac{\mathbf u}{U}\, }[/math]
Time L/U [math]\displaystyle{ t^*\ = \frac{t}{L/U}\, }[/math]
Pressure: there is no natural selection for the pressure scale. Where dynamic effects are dominant i.e. high velocity flows
[math]\displaystyle{ p^* = \frac{p}{\rho U^2} }[/math]

Where viscous effects are dominant i.e. creeping flows

[math]\displaystyle{ p^* = \frac{p L}{\mu U} }[/math]

Substituting the scales the non-dimensionalized equation obtained is:

[math]\displaystyle{ \frac{\partial \mathbf{u^*}}{\partial t^*} + (\mathbf{u^*} \cdot \nabla^*) \mathbf{u^*}\ = -\nabla^* p^* + \frac{1}{Re} \nabla^{*2} \mathbf{u^*} + \frac{1}{Fr^2} \hat{g}. }[/math][4]

 

 

 

 

(1)

where [math]\displaystyle{ Fr }[/math] is the Froude number and [math]\displaystyle{ Re }[/math] is the Reynolds number ([math]\displaystyle{ Re = UL/\nu }[/math]).

Flows with large viscosity

For flows where viscous forces are dominant i.e. slow flows with large viscosity, a viscous pressure scale μU/L is used. In the absence of a free surface, the equation obtained is

[math]\displaystyle{ Re \left( \frac{\partial \mathbf{u^*}}{\partial t^*} + (\mathbf{u^*} \cdot \nabla^* )\mathbf{u^*} \right)\ = -\nabla^* p^* + \nabla^{*2} \mathbf{u^*}. }[/math]

 

 

 

 

(2)

Stokes regime

Scaling of equation (1) can be done, in a flow where inertia term is smaller than the viscous term i.e. when Re → 0 then inertia terms can be neglected, leaving the equation of a creeping motion.

[math]\displaystyle{ Re \frac{\partial \mathbf{u^*}}{\partial t^*} = -\nabla^* p^* + \nabla^{*2} \mathbf{u^*}. }[/math]

Such flows tend to have influence of viscous interaction over large distances from an object.[citation needed] At low Reynolds number the same equation reduces to a diffusion equation, named Stokes equation

[math]\displaystyle{ -\nabla^* p^* + \nabla^{*2} \mathbf{u^*}= \mathbf 0. }[/math]

Euler regime

Similarly if Re → ∞ i.e. when the inertia forces dominates, the viscous contribution can be neglected. The non-dimensionalized Euler equation for an inviscid flow is

[math]\displaystyle{ \frac{\partial \mathbf{u^*}}{\partial t} + (\mathbf{u^*} \cdot \nabla^* )\mathbf{u^*}\ = -\nabla^* p^*. }[/math][6]

When density varies due to both concentration and temperature

Density variation due to both concentration and temperature is an important field of study in double diffusive convection. If density changes due to both temperature and salinity are taken into account, then some more terms add to the Z-Component of momentum as follows:[7] [8]

[math]\displaystyle{ \frac{\partial W}{\partial t} + U \frac{\partial W}{\partial X} + W \frac{\partial W}{\partial Z}\ = -\frac{1}{\rho_o}\frac{\partial p_d}{\partial Z} + v \left(\frac{\partial^2 W}{\partial X^2} + \frac{\partial^2 W}{\partial Z^2}\right)\ - g \left(\beta_{s}\nabla{S} - \beta_{T}\nabla{T}\right) }[/math]

Where S is the salinity of the fluid, βT is the thermal expansion coefficient at constant pressure and βS is the coefficient of saline expansion at constant pressure and temperature.

Non dimensionalizing using the scale:

[math]\displaystyle{ S^* = \frac{S - S_B}{S_T - S_B} }[/math] and [math]\displaystyle{ T^* = \frac{T - T_B}{T_T - T_B} }[/math]

we get

[math]\displaystyle{ \frac{\partial W^*}{\partial t^*} + U^* \frac{\partial W^*}{\partial X^*} + W^* \frac{\partial W^*}{\partial Z^*}\ = -\frac{\partial p_d}{\partial Z^*} + Pr \left(\frac{\partial^2 W^*}{\partial X^{*2}} + \frac{\partial^2 W^*}{\partial Z^{*2}}\right)\ - {Ra_s Pr_s S} + {Ra_T Pr_T T} }[/math]

where ST, TT denote the salinity and temperature at top layer, SB, TB denote the salinity and temperature at bottom layer, Ra is the Rayleigh Number, and Pr is the Prandtl Number. The sign of RaS and RaT will change depending on whether it stabilizes or destabilizes the system.

References

Footnotes

  1. Versteeg H.K, An introduction to computational fluid dynamics: the finite volume method, 2007, prentice hall, 9780131274983
  2. Patankar Suhas V., Numerical heat transfer and fluid flow, 1980, Taylor & Francis, 9780891165224
  3. Salvi Rodolfo, The Navier Stokes equation theory and numerical methods, 2002, M. Dekker, 9780824706722
  4. 4.0 4.1 Fox, Robert W.; Alan T. McDonald; Philip J. Pritchard (2006). Introduction to fluid mechanics (6th ed.). Hoboken, NJ: Wiley. p. 213–215. ISBN 9780471735588. https://archive.org/details/introductiontofl00foxr_614. 
  5. Tritton, D.J. (1988). Physical fluid dynamics (2nd ed.). Oxford [England]: Clarendon Press. pp. 55–58. ISBN 0198544898. 
  6. White, Frank M. (2003). Fluid mechanics (5th ed.). Boston: McGraw-Hill. pp. 188–189. ISBN 9780072402179. https://archive.org/details/fluidmechanics00whit_184. 
  7. On the relationship between finger width, velocity, and fluxes in thermohaline convection, 2009, K. R. Sreenivas, O. P. Singh & J. Srinivasan, Phys. Fluids (American Institute of Physics) 21(2), pp. 026601.
  8. Approximation of the hydrostatic Navier–Stokes system for density stratified flows by a multilayer model. Kinetic interpretation and numerical validation, E. Audusse, M.-O. Bristeau, M. Pelanti, J. Sainte-Marie, Université Paris 13, Institut Galilée, 99 avenue Jean-Baptiste Clément, 93430 Villetaneuse, France. INRIA Rocquencourt, B.P. 105, 78153 Le Chesnay Cedex, France. Saint-Venant Laboratory, 6 quai Watier, 78400 Chatou, France.

Other

  • Y. Cengel and J. Cimbala, FLUID MECHANICS: Fundamentals and Applications, 4th Edition, McGraw-Hill Education, 2018 (see p521, section 10.2. Nondimensionalized Equations of Motion).

Further reading