Continuum Mechanics
This article is about ‘continuum mechanics’ and we first consider what is meant by this.
Continuum mechanics is a branch of mechanics that deals with the mechanical behaviour of materials modelled as a continuous mass rather than as discrete particles. Augustin-Louis Cauchy was the first to formulate such models.
It may sound somewhat insignificant at first, but in the formulations of physics it makes a huge difference whether one understands the structure of matter as a continuum or as a lattice of discrete particles. Ultimately, both formulations must merge in the limiting case of infinitesimal lattice distances.
Modelling an object as a continuum assumes that the substance of the object completely fills the space it occupies. Modelling objects in this way ignores the fact that matter is made of atoms, and so is not continuous. On scales much greater than that of inter-atomic distances, such models are highly accurate.
Fundamental physical laws such as the conservation of mass, the conservation of momentum, and the conservation of energy may be applied to such models to derive differential equations describing the behaviour of such objects.
Continuum mechanics deals with physical properties of eg. fluids which are independent of any particular coordinate system in which they are observed. These physical properties are then represented by tensors, which are mathematical objects that have the required property of being independent of a specific coordinate system.
We also need to introduce the concept of classical field theories. What is this?
A classical field theory is a physical theory that predicts how physical fields interact with matter through field equations. The term ‘classical field theory’ is reserved for describing those physical theories that describe electromagnetism and gravitation, two of the fundamental forces of nature.
A physical field can be thought of as the assignment of a physical quantity at each point of space and time.
Fluid dynamics has fields of pressure, density, and flow rate that are connected by conservation laws for energy and momentum. The mass continuity equation is a continuity equation, representing the conservation of mass.
In physics and engineering, fluid dynamics is a subdiscipline of fluid mechanics that describes the flow of fluids—liquids and gases. 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 the First Law of Thermodynamics).
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.
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 well-known equation of state of an ideal gas:
p · V = n · R · T
where p is pressure, ρ is density, T the absolute temperature, while R is the gas constant and M is the molar mass for a particular gas.
- changes due to different field strengths,
- any time dependences of the field.
- the partial time derivative at a fixed location (local part)
- the change due to the movement of the fluid element because of which the temperature T must be evaluated at a “shifted” location (convective part)
Gauss’ Theorem
Let us return to fluid dynamics and look at these physical observables:
The actual amount of gas inside a closed control volume V with the closed surface S is:
M = ∫ ρ(x, y, z) dxdydz where the integral extends over the entire control volume.
At this point we introduce the mass flow vector as follows:
Thus we find the relation for the mass conservation law in integral representation:
Gauss’ divergence theorem will not be derived here but applied to the above relation:
This yields
or
Recall the definition of the substantial derivative
Applying this now to the density yields
Combining this with the relation for mass conservation in differential form is
This relation is now the Lagrange version of the mass conservation law.
Euler Equation
We now need a relationship that describes the effects of the internal forces in the fluid in a control volume.
We make the assumption of an ideal liquid which means that with the exception of the pressure, we do not allow any other physical effects which may alter the momentum or the energy of the fluid i.e., no friction and no heat conduction.
Our fluid element is defined as a small volume e.g. in the form of a cube which feels pressure from all sides.
For the pressure change in x-direction we find:
The resulting force (note that pressure is defined a force per area) is then
The components of the force in all three spatial dimensions is
In a flow, the velocity is generally neither constant in time nor constant in space. The fluid therefore flows either slower or faster past the lateral surfaces of the fluid volume under consideration. Frictional forces occur, which are all the greater the more viscous the fluid. However, since we assume a frictionless flow and thus a non-viscous fluid, such shear forces do not occur in the flow.
Other forces that act on the fluid, but generally cannot be neglected, are field forces such as those caused by gravity. Air flow on earth or fluid flow in pipes, for example, are significantly influenced by gravity.
In addition to the pressure force Fp, field forces Fg also act on a fluid element (the weight force is only meant to be representative of other possible field forces at this point). The sum of both forces then ultimately corresponds to the resulting accelerating force Facc, which influences the movement of the fluid element.
According to Newton’s second axiom, the accelerating force Facc leads to a corresponding change in velocity (acceleration), which depends on the mass dm of the fluid element under consideration. This material acceleration is also called total or substantial acceleration asub. The control volume is abbreviated as dV:
The substantial or material acceleration of a fluid element has two reasons:
- A flow and thus the speed of a fluid element changes not only in time but also in space.
- The temporal change therefore causes a “local acceleration” and the local change, which is nothing other than a flow transport from one location to another one, causes a “convective acceleration”.
One could now conclude that in a stationary flow no acceleration acts on the fluid elements. But this is not the case!
In a stationary flow, the velocity no longer changes in time at a fixed location. In general, a fluid element must permanently change its velocity when flowing, for example at narrow points in a pipe. Despite the lack of temporal change at a fixed location, the fluid element is nevertheless accelerated when it changes its location. This occurs through flow transport and is therefore referred to as “convective acceleration”.
Let us summarize at this point:
- The local acceleration is due to the time-varying flow velocity of a transient flow at a fixed location,
- the convective acceleration is due to the flow velocity changing from place to place,
- both components together result in the observable acceleration and are therefore referred to as substantial or material acceleration,
- the local component is omitted in the case of a stationary flow, as the velocity at a fixed location then no longer changes in time.
Energy Conservation
Let us consider a sample fluid element on a streamline:
In order to allow flowing in the indicated direction we postulate:
The total derivative of the velocity v of the fluid element is then
Division by dt
and eventually
An accelerating tangential force acts on the fluid element of mass dm in direstion s:
This equation only describes the effect of a temporal or local change of the velocity, but not the cause for the acceleration which only stems from the external forces acting on the fluid element.
Note that
This is the effective force due to the pressure difference between the left and the right side of the fluid element.
Let us consider the different types of acceleration i.e., the local and convective ones:
The formulation of the convective acceleration in three dimensions is a bit more difficult beacuse each velocity component depends on three spatial coordinates:
which translates into the three-dimensional vector notation:
or in a compact notation
and now using the Nabla-operator:
Thus the substantial derivative for the velocity reads:
Recall
We set the last two equations equal to each other
and eventually find the so-called Euler equation:
We now apply this equation to just streamline in one dimension:
In one dimension we work with
Multiplication with ds and ρ yields
ρ dv v + dp = -g ρ dz
or
dp + ρ dv v + g ρ dz = 0
Integration then gives
and eventually
This describes the energy conservation and is no more than the Bernoulli equation.
The Navier-Stokes equation (for incompressible fluids)
Before we start to think about some generalizations of what we learned so far, it is worth while to take a look at the history of all these things.
Isaac Newton published his three-volume Principia with the laws of motion in 1686 and also defined the viscosity of a linear viscous fluid in the second book. This type of fluid is nowadays referred to as a Newtonian fluid.
In 1755, Euler derived the Euler equations from the laws of motion, with which the behaviour of viscosity-free fluids (liquids and gases) can be calculated.
French mathematician d’Alembert introduced the Eulerian approach, derived the local mass balance and formulated the d’Alembert paradox, according to which no force is exerted on a body in the direction of flow by the flow of viscosity-free fluids – which Euler had already proved earlier.
Because of this and other paradoxes of viscosity-free flows, it was clear that Euler’s equations of motion had to be supplemented.
Four researchers of that time ie. Navier, Poisson, de Saint-Venant and Stokes independently from each other formulated the momentum theorem for Newtonian fluids in differential form in the first half of the 19th century.
Navier (1827) and Poisson (1831) established the momentum equations after considering the effect of intermolecular forces.
In 1843, de Saint-Venant published a derivation of the momentum equations from Newton’s linear viscosity approach, two years before Stokes did so in 1845. However, the name Navier-Stokes equations became established for the momentum equations.
A significant advance in the understanding of viscous fluids was provided by Prandtl in 1904 with his boundary layer theory. From the middle of the 20th century, numerical fluid mechanics developed to such an extent that with its help solutions to the Navier-Stokes equations could be found for practical problems, which – as it turns out – correspond well with the real flow processes.
The Navier Stokes equation is a mathematical model of the flow of linear-viscous Newtonian liquids and gases referred to as fluids. We will see that the equations are an extension of the Euler equations of fluid mechanics by terms describing viscosity.
The Navier-Stokes equation hence is a universal formula to describe the flow of fluids.
Just to not discourage the reader it is worth to mention that the Navier-Stokes equation could not be solved in the three-dimensional case until today and is one of the millennium problems of mathematics.
Analytically, the Navier-Stokes equation can so far only be solved under simplifying assumptions which we will also use in our article.
In the two-dimensional case, the existence of such solutions has already been successfully investigated. However, the Navier-Stokes equation can be solved numerically, which is used for the simulation of fluids, for example.
Especially in physics, the Navier-Stokes equations means the momentum equation for flows. In a broader sense i.e., in fluid mechanics, this momentum equation is extended by the continuity equation and the energy equation and then forms a system of non-linear partial differential equations of second order.
This is the fundamental mathematical model of fluid mechanics. In particular, the equations also represent turbulence and boundary layers.
If the Navier-Stokes equations are formulated by using dimensionless quantities, various dimensionless ratios such as the Reynolds number or the Prandtl number can be derived.
In what follows we will derive the Navier-Stokes equation for incompressible fluids, which just means at this point that the density is constant and doesn’t change across the three-dimensional movement of the fluid through a pipe. We will also look at the difference between compressible and incompressible fluids and how to simplify the Navier-Stokes equation for certain systems.
We want to understand how the Navier-Stokes equation is derived from Newton’s law. To be able to describe the flow of a volume element ΔV with a mass of Δm = ρ · ΔV we start from the well-known 2nd law of Newton:
As we derived earlier, the expression for the material derivative of velocity is:
To repeat the meaning we consider the temporal change of at the same location and the temporal change when the fluid flows from one location to another.
Therefore, the term is also called substantial acceleration.
We know the right hand side of Newton’s 2nd law and next we need to know the forces acting on a volume element ΔV of the fluid. As we will see there are various forces which need to taken into consideration. We recall the formula for the force due to the pressure differences between two different locations.
We now also consider friction forces which also contribute to the right hand side of Newton’s 2nd law. So, if the friction within a liquid is not negligible, a force or pressure is necessary to move a liquid uniformly in relation e.g. to a pipe system. How great the necessary thrust force is depends on the viscosity of the fluid.
Considering the fluid as a stack of layers, the movement of the layers will cause friction between them if the velocity of the layers is different. If, for example, a cover glass is placed on a drop of liquid and moved slowly and evenly along the base on the thin layer of liquid, a force F is required to maintain the movement. This force is proportional to the area A of the glass, the speed v of the movement and the viscosity; moreover, the force is inversely proportional to the thickness Δx of the liquid layer. Altogether, therefore, the force F necessary to overcome friction is
More generally speaking we expose the liquid to two plates which are moving with different velocities v1 and v2
Thus we find for the force F to overcome the friction
A volume element dv experiences frictional forces at the top and the bottom. The area A becomes dy · dz
F1 is the force acting on top of dV at position x:
F2 is the force acting on the bottom of dV at position x+dx
The resulting force F1 + F2 is then
Note that
is the second derivative of vz, with respect to dx
The resulting force for F2 with regard to variations on x-direction now becomes
However, this is not yet the entire force component in the z-direction because we also need to include changes of velocity in y- and z-directions:
Similar expressions can be derived for Fx and Fy and we introduce at this point the Laplace operator. This will allow us to write the formula for friction forces in a compact manner.
The expression for the friction force thus becomes:
There are many other forces which may contribute and to include them in our calculation we recall
As an example we consider the gravitational acceleration which may play a role in many of our applications. Then our expression for external forces takes the form
We now consider the total expression for any type of forces in accordance with Newton’s 2nd Law:
The left hand side of Newton’s 2nd Law corresponds to the total force:
which implies this complicated looking expression:
We divide this by dV and arrive at the Navier-Stokes equation for for incompressible fluids then
To better understand the meaning of all these terms it therefore makes sense to illustrate the individual summands with an example i.e., we apply the Navier-Stokes logic to a crowd of people.
What do all the summands mean in this case?
We can easily see here the analogies to the movement of particles or infinitesimally small volume elements in a liquid or gas flow, just replace ‘people’ by ‘fluid elements’.
If we set herein η = 0 we are back to frictionless movements of the liquid and the Navier-Stokes equation becomes the Euler equation which we derived before and which is prevailingly used for calculations of flows in pipe systems:
We can further simplify this equation if we consider very viscous fluids where the time dependence of the velocity field and the inertia term can be neglected i.e.,
which then leads us to the simpler form:
This is the so called Stokes Equation which is used in the case of very viscous fluids.
Before we continue the discussion of the Navier-Stokes equation, we need to take a closer look at another partial differential equation which is of relevance for our considerations, the continuity equation.
The Continuity Equation
In physics, the ‘law of conservation’ is the formulation of the observed fact that the value of a quantity, called a conservation quantity, does not change in certain physical process. In a closed system, conservation variables do not change. Typical conservation quantities are mass, electrical charge, energy, momentum, angular momentum and others.
the corresponding current density is defined as
A continuity equation is a particular partial differential equation associated with a conservation variable (see above). It links the temporal change of the density belonging to this conservation variable with the spatial change of its current density.
Recall the form of the substantial derivation which represents the total differential
This means in our case that if there are no sources in the flow then we don’t see any changes in density. That’s the case of incompressible fluids. The situation is different if we look at compressible fluids.
In fluid mechanics, the continuity law for (incompressible) fluids follows from the continuity equation.
The continuity law states (in integral form) that the mass flow of a fluid in a pipe is independent of where it is measured.
The differential form is the continuity equation.
It applies in both frictionless and frictional cases for steady-state ie. time-independent and for transient flows of incompressible fluids, but not for transient flows of compressible fluids.
For incompressible (non-compressible) fluids, continuity also applies to the volume flow.
Again, for illustration purposes, for incompressible fluids we note that the same volume emerges from a pipe section that simultaneously enters at the other end.
The entering volume is V1 = A1 · v1, the leaving volume is V2 = A2 · v2. Because of V1 = V2 we in the narrower part of the pipe the displacement x2 is greater than x1 by the same factor as the cross-section A1 is greater than A2. The same applies to the flow velocities averaged over the cross-section
On the other hand, for compressible and non-stationary fluids that can change their density, the following applies to the mass flow rate:
The mass flow rate is conservation variable i.e., the mass which enters from one side has to leave on the other side, but compared with the incompressible case, we now have to consider also the density which may be subject to change, i.e.,
For example, the density of the fluid can change if the temperature of the fluid changes between the beginning and the end of the pipe. If the density decreases, a larger volume must come out in the same time.
As an illustration we look at cars in a traffic jam due to narrowing lanes. They also behave according to the law of continuity when the lane narrows. The distance between the cars may be considered as density. If the cross-section is large, the density is low, the speed is high and the traffic flows freely. In the congestion before the narrowing, the density is high and the speed is low. In the constriction, the cross-section is small, the speed and density are medium, and the vehicle throughput is the same in all cases – provided that no car leaves or joins the road.
Navier-Stokes equation for compressible fluids and nondimensionalisation
This section is the most difficult one to understand and is added for the sake of a complete presentation of the topic.
It remains to consider the Navier-Stokes equation for compressible fluids. Essentially, the extension of the equation derived so far consists in the addition of a further term, the derivation of which we will spare ourselves here:
To the left of the equals sign is the substantial acceleration of the fluid elements and the term formed with the Nabla operator represents
Please note that this is the already known Navier-Stokes equation plus a term which caters for the effects of being now a compressible fluid.
As it was mentioned before so far nobody was able to solve this equation analytically; only numerical procedures could be developed following proper approximation schemes by truncating polynomials, but it remains a numerical solution only and not a thorough one.
What can we do to make our life simpler now?
In order to make this very confusing equation a little more “digestible”, we resort to a common method in mathematics, the nondimensionalisation of the Navier-Stokes equation. We do not deal with this in detail, but apply it without questioning or even deriving the individual steps.
But a little bit of background will be helpful:
Nondimensionalisation is the partial or complete removal of dimensioned quantities such as units of measurement from a physical equation by a suitable substitution of all variables. With the help of dimensionless variables and dimensioned constants, effects resulting from the choice of the system of units can be eliminated and intrinsic constants of the system, such as characteristic lengths, times or frequencies, can be found. The technique is therefore suitable for simplifying systems of differential equations. That’s exactly what we need now!
The Navier-Stokes equations can be nondimensionalised with characteristic measures of the entire flow area for the length L, the velocity v∞ and the density ρ∞. This gives rise to the dimensionless quantities ie. we replace the variables (parameters without an asterisk) by these ones:
This eventually leads to the dimensionless Navier-Stokes equation which looks very similar to the previous version but exclusively works with dimensionless (or scaled) variables:
with
This is the Reynolds number, a dimensionless parameter, which characterizes the flow in terms of the ratio of inertial to shear forces. In the following we omit the indices which were just introduced to outline the nondimensionalization procedure.
The Reynolds number is thus defined as
The characteristic length L, or reference length, is to be defined for the respective problem.
For flow bodies, the length of the body in the direction of flow is usually chosen. For resistance bodies, the width or height transverse to the direction of flow is usually taken as the characteristic length, for pipe flows the radius or diameter of the pipe and for flumes the depth or width at the flume surface.
The kinematic viscosity ν (Greek: nu) of the fluid differs from the dynamic viscosity η = ν ρ by the factor ρ. Hence the Reynolds number can also be expressed in terms of the kinematic viscosity ν:
If the Reynolds number exceeds a (problem-dependent) critical value, a previously laminar flow becomes susceptible to the smallest disturbances. Accordingly, for Re > Recritical a change from laminar to turbulent flow is to be expected. In ideal fluids (incompressible, no friction) there is no viscosity and the Reynolds number is infinite because the quotient of velocity and viscosity becomes infinite.