Commonly used SI units are given in parentheses although the equations are general in nature and other units can be entered into them or units can be removed at all by nondimensionalization.
Note that only we use column vectors (in the Cartesian coordinate system) above for clarity, but the equation is written using physical components (which are neither covariants ("column") nor contravariants ("row") ).[5] However, if we chose a non-orthogonal curvilinear coordinate system, then we should calculate and write equations in covariant ("row vectors") or contravariant ("column vectors") form.
After an appropriate change of variables, it can also be written in conservation form:
where j is the momentum density at a given space-time point, F is the flux associated to the momentum density, and s contains all of the body forces per unit volume.
Differential derivation
Let us start with the generalized momentum conservation principle which can be written as follows: "The change in system momentum is proportional to the resulting force acting on this system". It is expressed by the formula:[6]
where is momentum at time t, and is force averaged over . After dividing by and passing to the limit we get (derivative):
Let us analyse each side of the equation above.
Right side
The X component of the forces acting on walls of a cubic fluid element (green for top-bottom walls; red for left-right; black for front-back).In the top graph we see approximation of function (blue line) using a finite difference (yellow line). In the bottom graph we see "infinitely many times enlarged neighborhood of point " (purple square from the upper graph). In the bottom graph, the yellow line is completely covered by the blue one, thus not visible. In the bottom figure, two equivalent derivative forms have been used: ], and the designation was used.
Surface forces act on walls of the cubic fluid element. For each wall, the X component of these forces was marked in the figure with a cubic element (in the form of a product of stress and surface area e.g. with units ).
Explanation of the value of forces (approximations and minus signs) acting on the cube walls.
It requires some explanation why stress applied to the walls covering the coordinate axes takes a minus sign (e.g. for the left wall we have ). For simplicity, let us focus on the left wall with tension . The minus sign is due to the fact that a vector normal to this wall is a negative unit vector. Then, we calculated the stress vector by definition , thus the X component of this vector is (we use similar reasoning for stresses acting on the bottom and back walls, i.e.: ).
The second element requiring explanation is the approximation of the values of stress acting on the walls opposite the walls covering the axes. Let us focus on the right wall where the stress is an approximation of stress from the left wall at points with coordinates and it is equal to . This approximation suffices since, as goes to zero, goes to zero as well, by definition of partial derivative.
A more intuitive representation of the value of approximation in point has been shown in the figure below the cube. We proceed with similar reasoning for stress approximations .
Adding forces (their X components) acting on each of the cube walls, we get:
After ordering and performing similar reasoning for components (they have not been shown in the figure, but these would be vectors parallel to the Y and Z axes, respectively) we get:
We can then write it in the symbolic operational form:
There are mass forces acting on the inside of the control volume. We can write them using the acceleration field (e.g. gravitational acceleration):
Left side
Let us calculate momentum of the cube:
Because we assume that tested mass (cube) is constant in time, so
where Ω represents the control volume. Since this equation must hold for any control volume, it must be true that the integrand is zero, from this the Cauchy momentum equation follows. The main step (not done above) in deriving this equation is establishing that the derivative of the stress tensor is one of the forces that constitutes Fi.[1]
The Cauchy momentum equation can also be put in the following form:
Cauchy momentum equation(conservation form)
simply by defining:
where j is the momentum density at the point considered in the continuum (for which the continuity equation holds), F is the flux associated to the momentum density, and s contains all of the body forces per unit volume. u ⊗ u is the dyad of the velocity.
Here j and s have same number of dimensions N as the flow speed and the body acceleration, while F, being a tensor, has N2.[note 1]
In the Eulerian forms it is apparent that the assumption of no deviatoric stress brings Cauchy equations to the Euler equations.
Convective acceleration
An example of convective acceleration. The flow is steady (time-independent), but the fluid decelerates as it moves down the diverging duct (assuming incompressible or subsonic compressible flow).
A significant feature of the Navier–Stokes equations is the presence of convective acceleration: the effect of time-independent acceleration of a flow with respect to space. While individual continuum particles indeed experience time dependent acceleration, the convective acceleration of the flow field is a spatial effect, one example being fluid speeding up in a nozzle.
Regardless of what kind of continuum is being dealt with, convective acceleration is a nonlinear effect. Convective acceleration is present in most flows (exceptions include one-dimensional incompressible flow), but its dynamic effect is disregarded in creeping flow (also called Stokes flow). Convective acceleration is represented by the nonlinear quantity u ⋅ ∇u, which may be interpreted either as (u ⋅ ∇)u or as u ⋅ (∇u), with ∇u the tensor derivative of the velocity vector u. Both interpretations give the same result.[7]
Advection operator vs tensor derivative
The convective acceleration (u ⋅ ∇)u can be thought of as the advection operatoru ⋅ ∇ acting on the velocity field u.[7] This contrasts with the expression in terms of tensor derivative ∇u, which is the component-wise derivative of the velocity vector defined by [∇u]mi = ∂m vi, so that
where the Feynman subscript notation ∇a is used, which means the subscripted gradient operates only on the factor a.
Lamb in his famous classical book Hydrodynamics (1895),[8] used this identity to change the convective term of the flow velocity in rotational form, i.e. without a tensor derivative:[9][10]
where the vector is called the Lamb vector. The Cauchy momentum equation becomes:
Using the identity:
the Cauchy equation becomes:
In fact, in case of an external conservative field, by defining its potential φ:
In case of a steady flow the time derivative of the flow velocity disappears, so the momentum equation becomes:
And by projecting the momentum equation on the flow direction, i.e. along a streamline, the cross product disappears due to a vector calculus identity of the triple scalar product:
If the stress tensor is isotropic, then only the pressure enters: (where I is the identity tensor), and the Euler momentum equation in the steady incompressible case becomes:
In the steady incompressible case the mass equation is simply:
that is, the mass conservation for a steady incompressible flow states that the density along a streamline is constant. This leads to a considerable simplification of the Euler momentum equation:
The convenience of defining the total head for an inviscid liquid flow is now apparent:
in fact, the above equation can be simply written as:
That is, the momentum balance for a steady inviscid and incompressible flow in an external conservative field states that the total head along a streamline is constant.
Irrotational flows
The Lamb form is also useful in irrotational flow, where the curl of the velocity (called vorticity) ω = ∇ × u is equal to zero. In that case, the convection term in reduces to
Stresses
The effect of stress in the continuum flow is represented by the ∇p and ∇ ⋅ τ terms; these are gradients of surface forces, analogous to stresses in a solid. Here ∇p is the pressure gradient and arises from the isotropic part of the Cauchy stress tensor. This part is given by the normal stresses that occur in almost all situations. The anisotropic part of the stress tensor gives rise to ∇ ⋅ τ, which usually describes viscous forces; for incompressible flow, this is only a shear effect. Thus, τ is the deviatoric stress tensor, and the stress tensor is equal to:[11]
where I is the identity matrix in the space considered and τ the shear tensor.
All non-relativistic momentum conservation equations, such as the Navier–Stokes equation, can be derived by beginning with the Cauchy momentum equation and specifying the stress tensor through a constitutive relation. By expressing the shear tensor in terms of viscosity and fluid velocity, and assuming constant density and viscosity, the Cauchy momentum equation will lead to the Navier–Stokes equations. By assuming inviscid flow, the Navier–Stokes equations can further simplify to the Euler equations.
The divergence of the stress tensor can be written as
The effect of the pressure gradient on the flow is to accelerate the flow in the direction from high pressure to low pressure.
As written in the Cauchy momentum equation, the stress terms p and τ are yet unknown, so this equation alone cannot be used to solve problems. Besides the equations of motion—Newton's second law—a force model is needed relating the stresses to the flow motion.[12] For this reason, assumptions based on natural observations are often applied to specify the stresses in terms of the other flow variables, such as velocity and density.
External forces
The vector field f represents body forces per unit mass. Typically, these consist of only gravity acceleration, but may include others, such as electromagnetic forces. In non-inertial coordinate frames, other "inertial accelerations" associated with rotating coordinates may arise.
Often, these forces may be represented as the gradient of some scalar quantity χ, with f = ∇χ in which case they are called conservative forces. Gravity in the z direction, for example, is the gradient of −ρgz. Because pressure from such gravitation arises only as a gradient, we may include it in the pressure term as a body force h = p − χ. The pressure and force terms on the right-hand side of the Navier–Stokes equation become
It is also possible to include external influences into the stress term rather than the body force term. This may even include antisymmetric stresses (inputs of angular momentum), in contrast to the usually symmetrical internal contributions to the stress tensor.[13]
Nondimensionalisation
In order to make the equations dimensionless, a characteristic length r0 and a characteristic velocity u0 need to be defined. These should be chosen such that the dimensionless variables are all of order one. The following dimensionless variables are thus obtained:
Substitution of these inverted relations in the Euler momentum equations yields:
and can be eventually conservation equations. The limit of high Froude numbers (low external field) is thus notable for such equations and is studied with perturbation theory.
↑ In 3D for example, with respect to some coordinate system, the vector j has 3 components, while the tensors σ and F have 9 (3×3), so the explicit forms written as matrices would be:
Note, however, that if symmetrical, F will only contain 6 degrees of freedom. And F's symmetry is equivalent to σ's symmetry (which will be present for the most common Cauchy stress tensors), since dyads of vectors with themselves are always symmetrical.
Related Research Articles
In physics, specifically in electromagnetism, the Lorentz force is the combination of electric and magnetic force on a point charge due to electromagnetic fields. A particle of charge q moving with a velocity v in an electric field E and a magnetic field B experiences a force of
The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).
In fluid dynamics, potential flow or irrotational flow refers to a description of a fluid flow with no vorticity in it. Such a description typically arises in the limit of vanishing viscosity, i.e., for an inviscid fluid and with no vorticity present in the flow.
The vorticity equation of fluid dynamics describes the 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 governing equation is:
Poisson's equation is an elliptic partial differential equation of broad utility in theoretical physics. For example, the solution to Poisson's equation is the potential field caused by a given electric charge or mass density distribution; with the potential field known, one can then calculate electrostatic or gravitational (force) field. It is a generalization of Laplace's equation, which is also frequently seen in physics. The equation is named after French mathematician and physicist Siméon Denis Poisson.
In fluid dynamics, Stokes' law is an empirical law for the frictional force – also called drag force – exerted on spherical objects with very small Reynolds numbers in a viscous fluid. It was derived by George Gabriel Stokes in 1851 by solving the Stokes flow limit for small Reynolds numbers of the Navier–Stokes equations.
Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.
In the calculus of variations, a field of mathematical analysis, the functional derivative relates a change in a functional to a change in a function on which the functional depends.
This is a list of some vector calculus formulae for working with common curvilinear coordinate systems.
In physics, the Hamilton–Jacobi equation, named after William Rowan Hamilton and Carl Gustav Jacob Jacobi, is an alternative formulation of classical mechanics, equivalent to other formulations such as Newton's laws of motion, Lagrangian mechanics and Hamiltonian mechanics.
In differential geometry, the four-gradient is the four-vector analogue of the gradient from vector calculus.
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. In technology, it occurs in paint, MEMS devices, and in the flow of viscous polymers generally.
Toroidal coordinates are a three-dimensional orthogonal coordinate system that results from rotating the two-dimensional bipolar coordinate system about the axis that separates its two foci. Thus, the two foci and in bipolar coordinates become a ring of radius in the plane of the toroidal coordinate system; the -axis is the axis of rotation. The focal ring is also known as the reference circle.
In mathematics, parabolic cylindrical coordinates are a three-dimensional orthogonal coordinate system that results from projecting the two-dimensional parabolic coordinate system in the perpendicular -direction. Hence, the coordinate surfaces are confocal parabolic cylinders. Parabolic cylindrical coordinates have found many applications, e.g., the potential theory of edges.
Oblate spheroidal coordinates are a three-dimensional orthogonal coordinate system that results from rotating the two-dimensional elliptic coordinate system about the non-focal axis of the ellipse, i.e., the symmetry axis that separates the foci. Thus, the two foci are transformed into a ring of radius in the x-y plane. Oblate spheroidal coordinates can also be considered as a limiting case of ellipsoidal coordinates in which the two largest semi-axes are equal in length.
There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.
The derivation of the Navier–Stokes equations as well as their application and formulation for different families of fluids, is an important exercise in fluid dynamics with applications in mechanical engineering, physics, chemistry, heat transfer, and electrical engineering. A proof explaining the properties and bounds of the equations, such as Navier–Stokes existence and smoothness, is one of the important unsolved problems in mathematics.
In continuum mechanics, objective stress rates are time derivatives of stress that do not depend on the frame of reference. Many constitutive equations are designed in the form of a relation between a stress-rate and a strain-rate. The mechanical response of a material should not depend on the frame of reference. In other words, material constitutive equations should be frame-indifferent (objective). If the stress and strain measures are material quantities then objectivity is automatically satisfied. However, if the quantities are spatial, then the objectivity of the stress-rate is not guaranteed even if the strain-rate is objective.
Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.
The streamline upwind Petrov–Galerkin pressure-stabilizing Petrov–Galerkin formulation for incompressible Navier–Stokes equations can be used for finite element computations of high Reynolds number incompressible flow using equal order of finite element space by introducing additional stabilization terms in the Navier–Stokes Galerkin formulation.
This page is based on this Wikipedia article Text is available under the CC BY-SA 4.0 license; additional terms may apply. Images, videos and audio are available under their respective licenses.