# Fluid dynamics

Last updated

In physics and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluidsliquids and gases. It has several subdisciplines, including aerodynamics (the study of air and other gases in motion) and hydrodynamics (the study of liquids in motion). Fluid dynamics has a wide range of applications, including calculating forces and moments on aircraft, determining the mass flow rate of petroleum through pipelines, predicting weather patterns, understanding nebulae in interstellar space and modelling fission weapon detonation.

## Contents

Fluid dynamics offers a systematic structure—which underlies these practical disciplines—that embraces empirical and semi-empirical laws derived from flow measurement and used to solve practical problems. The solution to a fluid dynamics problem typically involves the calculation of various properties of the fluid, such as flow velocity, pressure, density, and temperature, as functions of space and time.

Before the twentieth century, hydrodynamics was synonymous with fluid dynamics. This is still reflected in names of some fluid dynamics topics, like magnetohydrodynamics and hydrodynamic stability, both of which can also be applied to gases. [1]

## Equations

The foundational axioms of fluid dynamics are the conservation laws, specifically, conservation of mass, conservation of linear momentum, and conservation of energy (also known as First Law of Thermodynamics). These are based on classical mechanics and are modified in quantum mechanics and general relativity. They are expressed using the Reynolds transport theorem.

In addition to the above, fluids are assumed to obey the continuum assumption. Fluids are composed of molecules that collide with one another and solid objects. However, the continuum assumption assumes that fluids are continuous, rather than discrete. Consequently, it is assumed that properties such as density, pressure, temperature, and flow velocity are well-defined at infinitesimally small points in space and vary continuously from one point to another. The fact that the fluid is made up of discrete molecules is ignored.

For fluids that are sufficiently dense to be a continuum, do not contain ionized species, and have flow velocities small in relation to the speed of light, the momentum equations for Newtonian fluids are the Navier–Stokes equations—which is a non-linear set of differential equations that describes the flow of a fluid whose stress depends linearly on flow velocity gradients and pressure. The unsimplified equations do not have a general closed-form solution, so they are primarily of use in Computational Fluid Dynamics. The equations can be simplified in a number of ways, all of which make them easier to solve. Some of the simplifications allow some simple fluid dynamics problems to be solved in closed form.[ citation needed ]

In addition to the mass, momentum, and energy conservation equations, a thermodynamic equation of state that gives the pressure as a function of other thermodynamic variables is required to completely describe the problem. An example of this would be the perfect gas equation of state:

${\displaystyle p={\frac {\rho R_{u}T}{M}}}$

where p is pressure, ρ is density, T the absolute temperature, while Ru is the gas constant and M is molar mass for a particular gas.

### Conservation laws

Three conservation laws are used to solve fluid dynamics problems, and may be written in integral or differential form. The conservation laws may be applied to a region of the flow called a control volume. A control volume is a discrete volume in space through which fluid is assumed to flow. The integral formulations of the conservation laws are used to describe the change of mass, momentum, or energy within the control volume. Differential formulations of the conservation laws apply Stokes' theorem to yield an expression which may be interpreted as the integral form of the law applied to an infinitesimally small volume (at a point) within the flow.

• Mass continuity (conservation of mass): The rate of change of fluid mass inside a control volume must be equal to the net rate of fluid flow into the volume. Physically, this statement requires that mass is neither created nor destroyed in the control volume, [2] and can be translated into the integral form of the continuity equation:
${\displaystyle {\partial \over \partial t}\iiint _{V}\rho \,dV=-\,{}}$${\displaystyle {\scriptstyle S}}$${\displaystyle {}\,\rho \mathbf {u} \cdot d\mathbf {S} }$
Above, ${\displaystyle \rho }$ is the fluid density, u is the flow velocity vector, and t is time. The left-hand side of the above expression is the rate of increase of mass within the volume and contains a triple integral over the control volume, whereas the right-hand side contains an integration over the surface of the control volume of mass convected into the system. Mass flow into the system is accounted as positive, and since the normal vector to the surface is opposite the sense of flow into the system the term is negated. The differential form of the continuity equation is, by the divergence theorem:
${\displaystyle \ {\partial \rho \over \partial t}+\nabla \cdot (\rho \mathbf {u} )=0}$
• Conservation of momentum: Newton's second law of motion applied to a control volume, is a statement that any change in momentum of the fluid within that control volume will be due to the net flow of momentum into the volume and the action of external forces acting on the fluid within the volume.
${\displaystyle {\frac {\partial }{\partial t}}\iiint _{\scriptstyle V}\rho \mathbf {u} \,dV=-\,{}}$${\displaystyle _{\scriptstyle S}}$${\displaystyle (\rho \mathbf {u} \cdot d\mathbf {S} )\mathbf {u} -{}}$${\displaystyle {\scriptstyle S}}$${\displaystyle {}\,p\,d\mathbf {S} }$${\displaystyle \displaystyle {}+\iiint _{\scriptstyle V}\rho \mathbf {f} _{\text{body}}\,dV+\mathbf {F} _{\text{surf}}}$
In the above integral formulation of this equation, the term on the left is the net change of momentum within the volume. The first term on the right is the net rate at which momentum is convected into the volume. The second term on the right is the force due to pressure on the volume's surfaces. The first two terms on the right are negated since momentum entering the system is accounted as positive, and the normal is opposite the direction of the velocity ${\displaystyle \mathbf {u} }$ and pressure forces. The third term on the right is the net acceleration of the mass within the volume due to any body forces (here represented by fbody). Surface forces, such as viscous forces, are represented by ${\displaystyle \mathbf {F} _{\text{surf}}}$, the net force due to shear forces acting on the volume surface. The momentum balance can also be written for a moving control volume. [3]
The following is the differential form of the momentum conservation equation. Here, the volume is reduced to an infinitesimally small point, and both surface and body forces are accounted for in one total force, F. For example, F may be expanded into an expression for the frictional and gravitational forces acting at a point in a flow.
${\displaystyle \ {D\mathbf {u} \over Dt}=\mathbf {F} -{\nabla p \over \rho }}$
In aerodynamics, air is assumed to be a Newtonian fluid, which posits a linear relationship between the shear stress (due to internal friction forces) and the rate of strain of the fluid. The equation above is a vector equation in a three-dimensional flow, but it can be expressed as three scalar equations in three coordinate directions. The conservation of momentum equations for the compressible, viscous flow case are called the Navier–Stokes equations. [2]
${\displaystyle \ \rho {Dh \over Dt}={Dp \over Dt}+\nabla \cdot \left(k\nabla T\right)+\Phi }$
Above, h is the specific enthalpy, k is the thermal conductivity of the fluid, T is temperature, and ${\displaystyle \Phi }$ is the viscous dissipation function. The viscous dissipation function governs the rate at which mechanical energy of the flow is converted to heat. The second law of thermodynamics requires that the dissipation term is always positive: viscosity cannot create energy within the control volume. [4] The expression on the left side is a material derivative.

### Compressible vs incompressible flow

All fluids are compressible to an extent; that is, changes in pressure or temperature cause changes in density. However, in many situations the changes in pressure and temperature are sufficiently small that the changes in density are negligible. In this case the flow can be modelled as an incompressible flow. Otherwise the more general compressible flow equations must be used.

Mathematically, incompressibility is expressed by saying that the density ρ of a fluid parcel does not change as it moves in the flow field, i.e.,

${\displaystyle {\frac {\mathrm {D} \rho }{\mathrm {D} t}}=0\,,}$

where D/Dt is the material derivative, which is the sum of local and convective derivatives. This additional constraint simplifies the governing equations, especially in the case when the fluid has a uniform density.

For flow of gases, to determine whether to use compressible or incompressible fluid dynamics, the Mach number of the flow is evaluated. As a rough guide, compressible effects can be ignored at Mach numbers below approximately 0.3. For liquids, whether the incompressible assumption is valid depends on the fluid properties (specifically the critical pressure and temperature of the fluid) and the flow conditions (how close to the critical pressure the actual flow pressure becomes). Acoustic problems always require allowing compressibility, since sound waves are compression waves involving changes in pressure and density of the medium through which they propagate.

### Newtonian vs non-Newtonian fluids

All fluids are viscous, meaning that they exert some resistance to deformation: neighbouring parcels of fluid moving at different velocities exert viscous forces on each other. The velocity gradient is referred to as a strain rate; it has dimensions ${\displaystyle T^{-1}}$. Isaac Newton showed that for many familiar fluids such as water and air, the stress due to these viscous forces is linearly related to the strain rate. Such fluids are called Newtonian fluids. The coefficient of proportionality is called the fluid's viscosity; for Newtonian fluids, it is a fluid property that is independent of the strain rate.

Non-Newtonian fluids have a more complicated, non-linear stress-strain behaviour. The sub-discipline of rheology describes the stress-strain behaviours of such fluids, which include emulsions and slurries, some viscoelastic materials such as blood and some polymers, and sticky liquids such as latex, honey and lubricants.[ citation needed ]

### Inviscid vs viscous vs Stokes flow

The dynamic of fluid parcels is described with the help of Newton's second law. An accelerating parcel of fluid is subject to inertial effects.

The Reynolds number is a dimensionless quantity which characterises the magnitude of inertial effects compared to the magnitude of viscous effects. A low Reynolds number (Re<<1) indicates that viscous forces are very strong compared to inertial forces. In such cases, inertial forces are sometimes neglected; this flow regime is called Stokes or creeping flow.

In contrast, high Reynolds numbers (Re>>1) indicate that the inertial effects have more effect on the velocity field than the viscous (friction) effects. In high Reynolds number flows, the flow is often modeled as an inviscid flow, an approximation in which viscosity is completely neglected. Eliminating viscosity allows the Navier–Stokes equations to be simplified into the Euler equations. The integration of the Euler equations along a streamline in an inviscid flow yields Bernoulli's equation. When, in addition to being inviscid, the flow is irrotational everywhere, Bernoulli's equation can completely describe the flow everywhere. Such flows are called potential flows, because the velocity field may be expressed as the gradient of a potential energy expression.

This idea can work fairly well when the Reynolds number is high. However, problems such as those involving solid boundaries may require that the viscosity be included. Viscosity cannot be neglected near solid boundaries because the no-slip condition generates a thin region of large strain rate, the boundary layer, in which viscosity effects dominate and which thus generates vorticity. Therefore, to calculate net forces on bodies (such as wings), viscous flow equations must be used: inviscid flow theory fails to predict drag forces, a limitation known as the d'Alembert's paradox.

A commonly used[ citation needed ] model, especially in computational fluid dynamics, is to use two flow models: the Euler equations away from the body, and boundary layer equations in a region close to the body. The two solutions can then be matched with each other, using the method of matched asymptotic expansions.

A flow that is not a function of time is called steady flow. Steady-state flow refers to the condition where the fluid properties at a point in the system do not change over time. Time dependent flow is known as unsteady (also called transient [6] ). Whether a particular flow is steady or unsteady, can depend on the chosen frame of reference. For instance, laminar flow over a sphere is steady in the frame of reference that is stationary with respect to the sphere. In a frame of reference that is stationary with respect to a background flow, the flow is unsteady.

Turbulent flows are unsteady by definition. A turbulent flow can, however, be statistically stationary. According to Pope: [7]

The random field ${\displaystyle U(x,t)}$ is statistically stationary if all statistics are invariant under a shift in time.

This roughly means that all statistical properties are constant in time. Often, the mean field is the object of interest, and this is constant too in a statistically stationary flow.

Steady flows are often more tractable than otherwise similar unsteady flows. The governing equations of a steady problem have one dimension fewer (time) than the governing equations of the same problem without taking advantage of the steadiness of the flow field.

### Laminar vs turbulent flow

Turbulence is flow characterized by recirculation, eddies, and apparent randomness. Flow in which turbulence is not exhibited is called laminar. The presence of eddies or recirculation alone does not necessarily indicate turbulent flow—these phenomena may be present in laminar flow as well. Mathematically, turbulent flow is often represented via a Reynolds decomposition, in which the flow is broken down into the sum of an average component and a perturbation component.

It is believed that turbulent flows can be described well through the use of the Navier–Stokes equations. Direct numerical simulation (DNS), based on the Navier–Stokes equations, makes it possible to simulate turbulent flows at moderate Reynolds numbers. Restrictions depend on the power of the computer used and the efficiency of the solution algorithm. The results of DNS have been found to agree well with experimental data for some flows. [8]

Most flows of interest have Reynolds numbers much too high for DNS to be a viable option, [9] given the state of computational power for the next few decades. Any flight vehicle large enough to carry a human (L > 3 m), moving faster than 20 m/s (72 km/h) is well beyond the limit of DNS simulation (Re = 4 million). Transport aircraft wings (such as on an Airbus A300 or Boeing 747) have Reynolds numbers of 40 million (based on the wing chord dimension). Solving these real-life flow problems requires turbulence models for the foreseeable future. Reynolds-averaged Navier–Stokes equations (RANS) combined with turbulence modelling provides a model of the effects of the turbulent flow. Such a modelling mainly provides the additional momentum transfer by the Reynolds stresses, although the turbulence also enhances the heat and mass transfer. Another promising methodology is large eddy simulation (LES), especially in the guise of detached eddy simulation (DES)—which is a combination of RANS turbulence modelling and large eddy simulation.

### Subsonic vs transonic, supersonic and hypersonic flows

While many flows (e.g. flow of water through a pipe) occur at low Mach numbers, many flows of practical interest in aerodynamics or in turbomachines occur at high fractions of M=1 (transonic flows) or in excess of it (supersonic or even hypersonic flows). New phenomena occur at these regimes such as instabilities in transonic flow, shock waves for supersonic flow, or non-equilibrium chemical behaviour due to ionization in hypersonic flows. In practice, each of those flow regimes is treated separately.

### Reactive vs non-reactive flows

Reactive flows are flows that are chemically reactive, which finds its applications in many areas such as combustion(IC engine), propulsion devices (Rockets, jet engines etc.), detonations, fire and safety hazards, astrophysics etc. In addition to conservation of mass, momentum and energy, conservation of individual species (for example, mass fraction of methane in methane combustion) need to be derived, where the production/depletion rate of any species are obtained by simultaneously solving the equations of chemical kinetics.

### Magnetohydrodynamics

Magnetohydrodynamics is the multi-disciplinary study of the flow of electrically conducting fluids in electromagnetic fields. Examples of such fluids include plasmas, liquid metals, and salt water. The fluid flow equations are solved simultaneously with Maxwell's equations of electromagnetism.

### Relativistic fluid dynamics

Relativistic fluid dynamics studies the macroscopic and microscopic fluid motion at large velocities comparable to the velocity of light. [10] This branch of fluid dynamics accounts the relativistic effects both from the special theory of relativity and the general theory of relativity. The governing equations are derived in Riemannian geometry for Minkowski spacetime.

### Other approximations

There are a large number of other possible approximations to fluid dynamic problems. Some of the more commonly used are listed below.

## Terminology

The concept of pressure is central to the study of both fluid statics and fluid dynamics. A pressure can be identified for every point in a body of fluid, regardless of whether the fluid is in motion or not. Pressure can be measured using an aneroid, Bourdon tube, mercury column, or various other methods.

Some of the terminology that is necessary in the study of fluid dynamics is not found in other similar areas of study. In particular, some of the terminology used in fluid dynamics is not used in fluid statics.

### Terminology in incompressible fluid dynamics

The concepts of total pressure and dynamic pressure arise from Bernoulli's equation and are significant in the study of all fluid flows. (These two pressures are not pressures in the usual sense—they cannot be measured using an aneroid, Bourdon tube or mercury column.) To avoid potential ambiguity when referring to pressure in fluid dynamics, many authors use the term static pressure to distinguish it from total pressure and dynamic pressure. Static pressure is identical to pressure and can be identified for every point in a fluid flow field.

A point in a fluid flow where the flow has come to rest (i.e. speed is equal to zero adjacent to some solid body immersed in the fluid flow) is of special significance. It is of such importance that it is given a special name—a stagnation point. The static pressure at the stagnation point is of special significance and is given its own name—stagnation pressure. In incompressible flows, the stagnation pressure at a stagnation point is equal to the total pressure throughout the flow field.

### Terminology in compressible fluid dynamics

In a compressible fluid, it is convenient to define the total conditions (also called stagnation conditions) for all thermodynamic state properties (e.g. total temperature, total enthalpy, total speed of sound). These total flow conditions are a function of the fluid velocity and have different values in frames of reference with different motion.

To avoid potential ambiguity when referring to the properties of the fluid associated with the state of the fluid rather than its motion, the prefix "static" is commonly used (e.g. static temperature, static enthalpy). Where there is no prefix, the fluid property is the static condition (i.e. "density" and "static density" mean the same thing). The static conditions are independent of the frame of reference.

Because the total flow conditions are defined by isentropically bringing the fluid to rest, there is no need to distinguish between total entropy and static entropy as they are always equal by definition. As such, entropy is most commonly referred to as simply "entropy".

## Related Research Articles

In physics, the Navier–Stokes equations, named after French engineer and physicist Claude-Louis Navier and Anglo-Irish physicist and mathematician George Gabriel Stokes, describe the motion of viscous fluid substances.

The vorticity equation of fluid dynamics describes evolution of the vorticity ω of a particle of a fluid as it moves with its flow, that is, the local rotation of the fluid . The equation is:

In fluid dynamics, the drag coefficient is a dimensionless quantity that is used to quantify the drag or resistance of an object in a fluid environment, such as air or water. It is used in the drag equation in which a lower drag coefficient indicates the object will have less aerodynamic or hydrodynamic drag. The drag coefficient is always associated with a particular surface area.

In physics and fluid mechanics, a boundary layer is the layer of fluid in the immediate vicinity of a bounding surface where the effects of viscosity are significant.

In fluid dynamics, the Euler equations are a set of quasilinear hyperbolic equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. The equations represent Cauchy equations of conservation of mass (continuity), and balance of momentum and energy, and can be seen as particular Navier–Stokes equations with zero viscosity and zero thermal conductivity. In fact, Euler equations can be obtained by linearization of some more precise continuity equations like Navier–Stokes equations in a local equilibrium state given by a Maxwellian. The Euler equations can be applied to incompressible and to compressible flow – assuming the flow velocity is a solenoidal field, or using another appropriate energy equation respectively. Historically, only the incompressible equations have been derived by Euler. However, fluid dynamics literature often refers to the full set – including the energy equation – of the more general compressible equations together as "the Euler equations".

A continuity equation in physics is an equation that describes the transport of some quantity. It is particularly simple and powerful when applied to a conserved quantity, but it can be generalized to apply to any extensive quantity. Since mass, energy, momentum, electric charge and other natural quantities are conserved under their respective appropriate conditions, a variety of physical phenomena may be described using continuity equations.

In fluid mechanics or more generally continuum mechanics, incompressible flow refers to a flow in which the material density is constant within a fluid parcel—an infinitesimal volume that moves with the flow velocity. An equivalent statement that implies incompressibility is that the divergence of the flow velocity is zero.

In fluid dynamics, d'Alembert's paradox is a contradiction reached in 1752 by French mathematician Jean le Rond d'Alembert. D'Alembert proved that – for incompressible and inviscid potential flow – the drag force is zero on a body moving with constant velocity relative to the fluid. Zero drag is in direct contradiction to the observation of substantial drag on bodies moving relative to fluids, such as air and water; especially at high velocities corresponding with high Reynolds numbers. It is a particular example of the reversibility paradox.

Stokes flow, also named creeping flow or creeping motion, is a type of fluid flow where advective inertial forces are small compared with viscous forces. The Reynolds number is low, i.e. . This is a typical situation in flows where the fluid velocities are very slow, the viscosities are very large, or the length-scales of the flow are very small. Creeping flow was first studied to understand lubrication. In nature this type of flow occurs in the swimming of microorganisms and sperm and the flow of lava. In technology, it occurs in paint, MEMS devices, and in the flow of viscous polymers generally.

In fluid dynamics, drag is a force acting opposite to the relative motion of any object moving with respect to a surrounding fluid. This can exist between two fluid layers or a fluid and a solid surface. Unlike other resistive forces, such as dry friction, which are nearly independent of velocity, drag forces depend on velocity. Drag force is proportional to the velocity for a laminar flow and the squared velocity for a turbulent flow. Even though the ultimate cause of a drag is viscous friction, the turbulent drag is independent of viscosity.

Aeroacoustics is a branch of acoustics that studies noise generation via either turbulent fluid motion or aerodynamic forces interacting with surfaces. Noise generation can also be associated with periodically varying flows. A notable example of this phenomenon is the Aeolian tones produced by wind blowing over fixed objects.

Fluid mechanics is the branch of physics concerned with the mechanics of fluids and the forces on them. It has applications in a wide range of disciplines, including mechanical, civil, chemical and biomedical engineering, geophysics, oceanography, meteorology, astrophysics, and biology.

In fluid mechanics, added mass or virtual mass is the inertia added to a system because an accelerating or decelerating body must move some volume of surrounding fluid as it moves through it. Added mass is a common issue because the object and surrounding fluid cannot occupy the same physical space simultaneously. For simplicity this can be modeled as some volume of fluid moving with the object, though in reality "all" the fluid will be accelerated, to various degrees.

The intent of this article is to highlight the important points of the derivation of the Navier–Stokes equations as well as its application and formulation for different families of fluids.

Pressure-correction method is a class of methods used in computational fluid dynamics for numerically solving the Navier-Stokes equations normally for incompressible flows.

Volume viscosity is a material property relevant for characterizing fluid flow. Common symbols are or . It has dimensions, and the corresponding SI unit is the pascal-second (Pa·s).

In fluid dynamics, hydrodynamic stability is the field which analyses the stability and the onset of instability of fluid flows. The study of hydrodynamic stability aims to find out if a given flow is stable or unstable, and if so, how these instabilities will cause the development of turbulence. The foundations of hydrodynamic stability, both theoretical and experimental, were laid most notably by Helmholtz, Kelvin, Rayleigh and Reynolds during the nineteenth century. These foundations have given many useful tools to study hydrodynamic stability. These include Reynolds number, the Euler equations, and the Navier–Stokes equations. When studying flow stability it is useful to understand more simplistic systems, e.g. incompressible and inviscid fluids which can then be developed further onto more complex flows. Since the 1980s, more computational methods are being used to model and analyse the more complex flows.

The Reynolds number is an important dimensionless quantity in fluid mechanics used to help predict flow patterns in different fluid flow situations. At low Reynolds numbers, flows tend to be dominated by laminar (sheet-like) flow, while at high Reynolds numbers turbulence results from differences in the fluid's speed and direction, which may sometimes intersect or even move counter to the overall direction of the flow. These eddy currents begin to churn the flow, using up energy in the process, which for liquids increases the chances of cavitation. The Reynolds number has wide applications, ranging from liquid flow in a pipe to the passage of air over an aircraft wing. It is used to predict the transition from laminar to turbulent flow, and is used in the scaling of similar but different-sized flow situations, such as between an aircraft model in a wind tunnel and the full size version. The predictions of the onset of turbulence and the ability to calculate scaling effects can be used to help predict fluid behaviour on a larger scale, such as in local or global air or water movement and thereby the associated meteorological and climatological effects.

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 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.

## References

1. Eckert, Michael (2006). The Dawn of Fluid Dynamics: A Discipline Between Science and Technology. Wiley. p. ix. ISBN   3-527-40513-5.
2. Anderson, J. D. (2007). Fundamentals of Aerodynamics (4th ed.). London: McGraw–Hill. ISBN   0-07-125408-0.
3. Nangia, Nishant; Johansen, Hans; Patankar, Neelesh A.; Bhalla, Amneet Pal S. (2017). "A moving control volume approach to computing hydrodynamic forces and torques on immersed bodies". Journal of Computational Physics. 347: 437–462. arXiv:. Bibcode:2017JCoPh.347..437N. doi:10.1016/j.jcp.2017.06.047.
4. White, F. M. (1974). Viscous Fluid Flow. New York: McGraw–Hill. ISBN   0-07-069710-8.
5. Shengtai Li, Hui Li "Parallel AMR Code for Compressible MHD or HD Equations" (Los Alamos National Laboratory) Archived 2016-03-03 at the Wayback Machine
6. "Transient state or unsteady state? -- CFD Online Discussion Forums". www.cfd-online.com.
7. See Pope (2000), p. 75.
8. See, for example, Schlatter et al, Phys. Fluids 21, 051702 (2009); doi : 10.1063/1.3139294
9. See Pope (2000), p. 344.
10. Landau, Lev Davidovich; Lifshitz, Evgenii Mikhailovich (1987). Fluid Mechanics. London: Pergamon. ISBN   0-08-033933-6.