The movement of fluid in the physical domain is driven by various properties. For the purpose of bringing the behavior of fluid flow to light and developing a mathematical model, those properties have to be defined precisely as to provide transition between the physical and the numerical domain. Velocity, pressure, temperature, density, and viscosity are the main properties that should be considered simultaneously when conducting a fluid flow examination. In accordance with the physical incidents such as combustion, multiphase flow, turbulent, mass transport, etc., those properties diversify enormously, which can be categorized into kinematic, transport, thermodynamic, and other miscellaneous properties
Thermo-fluid incidents directed by governing equations are based on the laws of conservation. The Navier-Stokes equations are the broadly applied mathematical model to examine changes on those properties during dynamic and/or thermal interactions. The equations are adjustable regarding the content of the problem and are expressed based on the principles of conservation of mass, momentum, and energy
Although some sources specify the expression of Navier-Stokes equations merely for conservation of momentum, some of them also use all equations of conservation of the physical properties. Regarding the flow conditions, the Navier-Stokes equations are rearranged to provide affirmative solutions in which the complexity of the problem either increases or decreases. For instance, having a numerical case of turbulence according to the pre-calculated Reynolds number, an appropriate turbulent model has to be applied to obtain credible results.
Despite the fact that motion of fluid is an exploratory topic for human beings, the evolution of mathematical models emerged at the end of 19
century after the industrial revolution. The initial appropriate description of the viscous fluid motion had been indicated in the paper “Principia” by Sir Isaac Newton (1687) in which dynamic behavior of fluids under constant viscosity was investigated
. Later, Daniel Bernoulli (1738) and Leonhard Euler (1755) subsequently derived the equation of inviscid flow which is now expressed as Euler’s inviscid equations. Even though Claude-Louis Navier (1827), Augustin-Louis Cauchy (1828), Siméon Denis Poisson (1829), and Adhémar St.Venant (1843) had carried out studies to explore the mathematical model of fluid flow, they had overlooked the viscous (frictional) force. In 1845, Sir George Stokes had derived the equation of motion of a viscous flow by adding Newtonian viscous terms, thereby the Navier-Stokes Equations had been brought to their final form which has been used to generate numerical solutions for fluid flow ever since
, and Stokes
, the mathematicians behind the Navier-Stokes equations
The observation method of fluid flow based on kinematic properties is a fundamental issue for generating a convenient mathematical model. Movement of fluid can be investigated with either Lagrangian or Eulerian methods
. Lagrangian description of fluid motion is based on monitoring a fluid particle which is large enough to detect properties. Between the initial coordinates at time
and coordinates of the same particle at time
millions of separate particles have to be examined through the path that is almost impossible to follow. In the Eulerian method, any specific particle across the path is not followed; instead, the velocity field as a function of time and position is examined. The missile example (Figure 3) precisely fits to emphasize these methods.
The Lagrangian formulation of motion is always time dependent. As
are the initial coordinates of a particle;
are coordinates of the same particle at time
. Description of motion for Lagrangian:
In the Eulerian method,
are the components of velocity at the point
is the time. The velocity components
are the unknowns which are functions of the independent variables
. The description of motion with the Eulerian method for any particular value of
The equations of conservation in the Eulerian system in which fluid motion is described are expressed as Continuity Equation for mass, Navier-Stokes Equations for momentum and Energy Equation for the first law of Thermodynamics. The equations are all considered simultaneously to examine fluid and flow fields.
The mass in the control volume can be neither created nor destroyed in accordance with physical laws. The conservation of mass, also expressed as Continuity Equation, states that the mass flow difference throughout system between inlet- and outlet-section is zero:
is velocity and gradient operator
While the density is constant, the flow is assumed incompressible and then continuity is simplified as below, which indicates a steady-state process:
The momentum in a control volume is kept constant, which implies conservation of momentum that we call ‘The Navier-Stokes Equations’. The description is set up in accordance with the expression of Newton’s Second Law of Motion:
is the net force applied to any particle,
is the acceleration, and
is the mass. In case the particle is a fluid, it is convenient to divide the equation to volume of particle to generate a derivation in terms of density as follows:
is the force exerted on the fluid particle per unit volume, and
is the applied force on the whole mass of fluid particles as below:
is the density of fluid and
is the gravitational acceleration. External forces which are deployed through the surface of fluid particles,
is expressed by pressure and viscous forces as shown below:
is expressed as stress tensor. According to the general deformation law of Newtonian viscous fluid given by Stokes,
is expressed as
Hence, Newton’s equation of motion can be specified in the form as follows:
Substitution of equation (10) into equation (11) results in the Navier-Stokes equations of Newtonian viscous fluid in one equation:
: Momentum convection
: Mass force
: Surface force
: Viscous force
where static pressure
and gravitational force
. The equation (12) is convenient for fluid and flow fields both transient and compressible.
indicates the substantial derivative as follows:
If the density of fluid is accepted to be constant, the equations are greatly simplified in which the viscosity coefficient
is assumed constant and
in equation (12). Thus, the Navier-Stokes equations for an incompressible three-dimensional flow can be expressed as follows:
For each dimension when the velocity is
are unknowns where a solution is sought by application of both continuity equation and boundary conditions. Besides, the energy equation has to be considered if any thermal interaction is available in the problem.
The Conservation of Energy is the first law of thermodynamics which states that the sum of the work and heat added to the system will result in the increase of energy of the system:
is the heat added to the system,
is the work done on the system, and
is the increment in the total energy of the system. One of the common types of energy equation is:
: Local change with time
: Convective term
: Pressure work
: Heat flux
: Heat dissipation term
The Navier-Stokes equations are extremely thorough in order to simulate physical incidents. Nonetheless, having a non-linear structure and various complexities, it is hardly possible to conduct an exact solution of those equations. Thus, with regard to the physical domain, both approaches and assumptions are partially applied to grind the equations. Beyond this, though, some assumptions have to be applied as to provide a reliable model in which the equation is carried out, up to a further step in terms of complexity such as turbulence.
The mathematical model merely gives ties among parameters which are part of the whole process. Hence, the solution of the Navier-Stokes equations can be realized with either analytical or numerical methods. The analytical method is the process that only compensates solutions in which non-linear and complex structures in the Navier-Stokes equations are ignored within several assumptions. It is only valid for simple / fundamental cases such as Couette flow, Poisellie flow, etc
. On the other hand, almost every case in fluid dynamics comprises non-linear and complex structures in the mathematical model which cannot be ignored to sustain reliability. Hence, the solution of the Navier-Stokes equations are carried out within several numerical methods, the omnipresence of Ordinary Differential Equations (ODEs) and Partial Differential Equations (PDEs). A computational analysis of fluid flow can be described as shown in Figure 4.
The analysis of fluid flow can be conducted in either steady (time-independent) or unsteady (time-dependent) condition depending on the physical incident. In case the fluid flow is steady, it means the motion of fluid and parameters do not rely on change in time, the term
where the continuity and momentum equations are re-derived as follows:
The Navier-Stokes equation in
While the steady flow assumption negates the effect of some non-linear terms and provides a convenient solution, variation of density is a hurdle that keeps the equation in a complex formation.
Due to the malleable structure of fluids, the compressibility of particles is a significant issue. Despite the fact that all types of fluid flow are compressible in a various range regarding molecular structure, most of them can be assumed to be incompressible in which the density changes are negligible. Thus, the term
is thrown away regardless of whether the flow is steady or not, as below:
The Navier-Stokes equation in
As incompressible flow assumption provides reasonable equations, the application of steady flow assumption concurrently enables us to ignore non-linear terms where
. Moreover, the density of fluid in high speed cannot be accepted as incompressible in which the density changes are important. “The Mach Number” is a dimensionless number that is convenient to investigate fluid flow, whether incompressible or compressible
is the Mach number,
is the velocity of flow, and
is the speed of sound at
at sea level. As in equation (24), when the Mach number is lower than 0.3, the assumption of incompressibility is acceptable. On the contrary, the change in density cannot be negligible in which density should be considered as a significant parameter. For instance, if the velocity of a car is higher than 100
, the suitable approach to conduct credible numerical analysis is the compressible flow. Apart from velocity, the effect of thermal properties on the density changes has to be considered in geophysical flows
The Reynolds number, the ratio of inertial and viscous effects, is also effective on Navier-Stokes equations to truncate the mathematical model. While
, the viscous effects are presumed negligible where viscous terms in Navier-Stokes equations are thrown away. The simplified form of Navier-Stokes equation, described as Euler equation, can be specified as follows
The Navier-Stokes equation in
Even though viscous effects are relatively important for fluids, the inviscid flow model partially provides a reliable mathematical model as to predict real process for some specific cases. For instance, high-speed external flow over bodies is a broadly used approximation where inviscid approach reasonably fits. While
, the inertial effects are assumed negligible where related terms in Navier-Stokes equations drop out. The simplified form of Navier-Stokes equations is called either creeping flow or Stokes flow
The Navier-Stokes equation in
Having tangible viscous effects, creeping flow is a suitable approach to investigate the flow of lava, swimming of microorganisms, flow of polymers, lubrication, etc.
The behavior of the fluid under dynamic conditions is a challenging issue that is compartmentalized as laminar and turbulent. The laminar flow is orderly at which motion of fluid can be predicted precisely. Except that, the turbulent flow has various hindrances, therefore it is hard to predict the fluid flow which shows a chaotic behavior. The Reynolds number, the ratio of inertial forces to viscous forces, predicts the behavior of fluid flow whether laminar or turbulent regarding several properties such as velocity, length, viscosity, and also type of flow. Whilst the flow is turbulent, a proper mathematical model is selected to carry out numerical solutions. Various turbulent models are available in literature and each of them has a slightly different structure to examine chaotic fluid flow.
Turbulent flow can be applied to the Navier-Stokes equations in order to conduct solutions to chaotic behavior of fluid flow. Apart from the laminar, transport quantities of the turbulent flow, it is driven by instantaneous values. Direct numerical simulation (DNS) is the approach to solving the Navier-Stokes equation with instantaneous values. Having district fluctuations varies in a broad range, DNS needs enormous effort and expensive computational facilities. To avoid those hurdles, the instantaneous quantities are reinstated by the sum of their mean and fluctuating parts as follows:
are velocity components and
is temperature. The differences among values are shown in Figure 5 both for steady and unsteady conditions:
Instead of instantaneous values which cause non-linearity, carrying out a numerical solution with mean values provides an appropriate mathematical model which is named “The Reynolds-averaged Navier-Stokes (RANS) equation”
. The fluctuations can be negligible for most engineering cases which cause a complex mathematical model. Thus, RANS turbulence model is a procedure to close the system of mean flow equations. The general form of The Reynolds-averaged Navier-Stokes (RANS) equation can be specified as follows:
The Navier-Stokes equation in
The turbulence model of RANS can also vary regarding methods such as k-omega, k-epsilon, k-omega-SST, and Spalart-Allmaras which have been used to seek a solution for different types of turbulent flow.
Likewise, large eddy simulation (LES) is another mathematical method for turbulent flow which is also comprehensively applied for several cases. Tough LES ensures more accurate results than RANS, it requires much more time and computer memory. As in DNS, LES considers to solve the instantaneous Navier-Stokes equations in time and three-dimensional space
The simple form of the Navier-Stokes equations only encompasses the change in properties such as velocity, pressure, and density under dynamic conditions for one phase laminar flow. Most engineering applications require further mathematical models to simulate physical incidents with the aim of obtaining affirmative results in the numerical domain. Some of the most common engineering problems and relevant mathematical models as to carry out numerical simulations can be given as below:
The Navier-Stokes equations cannot compensate the physical model of the flow at very small scales such as the motion of single bacteria — also called microfluidics. Thereby, it is convenient to either change or reinstate the Navier-Stokes equations with a suitable mathematical model. The Knudsen number (Kn) is a dimensionless number that is the ratio of the mean free path of molecular structure to observation scale. The preferred model in accordance with the Knudsen number is shown in Figure 6:
The Navier-Stokes equations have been considered an important issue by animation / video game companies
. For the purpose of generating a real-time animation, the application of Navier-Stokes equations provides stunning results when it comes to enhancing reality
. They are also part of the reason that modern video games appear to be realistic in many more ways than several years ago. Just think about the movement of flags in the wind: doesn’t this look absolutely realistic in modern games?!
: White, Frank (1991).Viscous Fluid Flow. 3rd Edition. McGraw-Hill Mechanical Engineering. ISBN-10: 0072402318.
: Stokes, George (1851). “On the Effect of the Internal Friction of Fluids on the Motion of Pendulums”. Transactions of the Cambridge Philosophical Society. 9: 8–106.
: White, Frank (2002). Fluid Mechanics. 4th edition. McGraw-Hill Higher Education. ISBN: 0-07-228192-8.
: Cebeci, T., Shao, J.P., Kafyeke, F., Laurendeau, E (2005). Computational Fluid Dynamics for Engineers. Horizon Publishing Inc. ISBN: 0-9766545-0-4.
: Bird, R.B., Stewart, W.E. and Lightfoot, E.N. (2001). Transport Phenomena, 2th edition. John Wiley & Sons. ISBN 0-471-41077-2.
Last updated: January 29th, 2019
We appreciate and value your feedback.