Euler equations (fluid dynamics) explained

In fluid dynamics, the Euler equations are a set of partial differential equations governing adiabatic and inviscid flow. They are named after Leonhard Euler. In particular, they correspond to the Navier–Stokes equations with zero viscosity and zero thermal conductivity.

The Euler equations can be applied to incompressible and compressible flows. The incompressible Euler equations consist of Cauchy equations for conservation of mass and balance of momentum, together with the incompressibility condition that the flow velocity is divergence-free. The compressible Euler equations consist of equations for conservation of mass, balance of momentum, and balance of energy, together with a suitable constitutive equation for the specific energy density of the fluid. Historically, only the equations of conservation of mass and balance of momentum were derived by Euler. However, fluid dynamics literature often refers to the full set of the compressible Euler equations – including the energy equation – as "the compressible Euler equations".

The mathematical characters of the incompressible and compressible Euler equations are rather different. For constant fluid density, the incompressible equations can be written as a quasilinear advection equation for the fluid velocity together with an elliptic Poisson's equation for the pressure. On the other hand, the compressible Euler equations form a quasilinear hyperbolic system of conservation equations.

The Euler equations can be formulated in a "convective form" (also called the "Lagrangian form") or a "conservation form" (also called the "Eulerian form"). The convective form emphasizes changes to the state in a frame of reference moving with the fluid. The conservation form emphasizes the mathematical interpretation of the equations as conservation equations for a control volume fixed in space (which is usefulfrom a numerical point of view).

History

The Euler equations first appeared in published form in Euler's article "Principes généraux du mouvement des fluides", published in Mémoires de l'Académie des Sciences de Berlin in 1757 (although Euler had previously presented his work to the Berlin Academy in 1752). Prior work included contributions from the Bernoulli family as well as from Jean le Rond d'Alembert.^[1]

The Euler equations were among the first partial differential equations to be written down, after the wave equation. In Euler's original work, the system of equations consisted of the momentum and continuity equations, and thus was underdetermined except in the case of an incompressible flow. An additional equation, which was called the adiabatic condition, was supplied by Pierre-Simon Laplace in 1816.

During the second half of the 19th century, it was found that the equation related to the balance of energy must at all times be kept for compressible flows, and the adiabatic condition is a consequence of the fundamental laws in the case of smooth solutions. With the discovery of the special theory of relativity, the concepts of energy density, momentum density, and stress were unified into the concept of the stress–energy tensor, and energy and momentum were likewise unified into a single concept, the energy–momentum vector.

Incompressible Euler equations with constant and uniform density

In convective form (i.e., the form with the convective operator made explicit in the momentum equation), the incompressible Euler equations in case of density constant in time and uniform in space are:

where:

is the flow velocity vector, with components in an N-dimensional space

u_1,u_2,...,u_N

	D\boldsymbol\Phi
	Dt

	\partial\boldsymbol\Phi
	\partialt

+v ⋅ \nabla\boldsymbol\Phi

, for a generic function (or field)

\boldsymbol\Phi

denotes its material derivative in time with respect to the advective field

and

\nablaw

is the gradient of the specific (with the sense of per unit mass) thermodynamic work, the internal source term, and

\nabla ⋅ u

is the flow velocity divergence.

represents body accelerations (per unit mass) acting on the continuum, for example gravity, inertial accelerations, electric field acceleration, and so on.

The first equation is the Euler momentum equation with uniform density (for this equation it could also not be constant in time). By expanding the material derivative, the equations become: $\begin + (\mathbf \cdot \nabla) \mathbf &= -\nabla w + \mathbf,\\ \nabla \cdot \mathbf &= 0.\end$

In fact for a flow with uniform density

\rho₀

the following identity holds:

\nabla w \equiv \nabla \left(\frac p \right) = \frac 1 \nabla p,

where

is the mechanic pressure. The second equation is the incompressible constraint, stating the flow velocity is a solenoidal field (the order of the equations is not causal, but underlines the fact that the incompressible constraint is not a degenerate form of the continuity equation, but rather of the energy equation, as it will become clear in the following). Notably, the continuity equation would be required also in this incompressible case as an additional third equation in case of density varying in time or varying in space. For example, with density nonuniform in space but constant in time, the continuity equation to be added to the above set would correspond to:

\frac = 0.

So the case of constant and uniform density is the only one not requiring the continuity equation as additional equation regardless of the presence or absence of the incompressible constraint. In fact, the case of incompressible Euler equations with constant and uniform density discussed here is a toy model featuring only two simplified equations, so it is ideal for didactical purposes even if with limited physical relevance.

The equations above thus represent respectively conservation of mass (1 scalar equation) and momentum (1 vector equation containing

scalar components, where

is the physical dimension of the space of interest). Flow velocity and pressure are the so-called physical variables.

In a coordinate system given by

\left(x_1,...,x_N\right)

the velocity and external force vectors

and

have components

(u_1,...,u_N)

and

\left(g_1,...,g_N\right)

, respectively. Then the equations may be expressed in subscript notation as:

$\begin + \sum_^N &= g_i,\\ \sum_^N &= 0.\end$

where the

and

subscripts label the N-dimensional space components, and

\delta_ij

is the Kroenecker delta. The use of Einstein notation (where the sum is implied by repeated indices instead of sigma notation) is also frequent.

Properties

Although Euler first presented these equations in 1755, many fundamental questions or concepts about them remain unanswered.

In three space dimensions, in certain simplified scenarios, the Euler equations produce singularities.^[2]

Smooth solutions of the free (in the sense of without source term: g=0) equations satisfy the conservation of specific kinetic energy: $\left(\frac u^2 \right) + \nabla \cdot \left(u^2 \mathbf u + w \mathbf u\right) = 0.$

In the one-dimensional case without the source term (both pressure gradient and external force), the momentum equation becomes the inviscid Burgers' equation: $+ u = 0.$

This model equation gives many insights into Euler equations.

Nondimensionalisation

In order to make the equations dimensionless, a characteristic length

r₀

, and a characteristic velocity

u₀

, 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:

\begin u^* & \equiv \frac, & r^* & \equiv \frac, \\[5pt] t^* & \equiv \frac t, & p^* & \equiv \frac, \\[5pt] \nabla^* & \equiv r_0 \nabla.\end

and of the field unit vector:

\hat\equiv \frac g.

Substitution of these inversed relations in Euler equations, defining the Froude number, yields (omitting the * at apix):

Euler equations in the Froude limit (no external field) are named free equations and are conservative. The limit of high Froude numbers (low external field) is thus notable and can be studied with perturbation theory.

Conservation form

See also: Conservation equation. The conservation form emphasizes the mathematical properties of Euler equations, and especially the contracted form is often the most convenient one for computational fluid dynamics simulations. Computationally, there are some advantages in using the conserved variables. This gives rise to a large class of numerical methodscalled conservative methods.

The free Euler equations are conservative, in the sense they are equivalent to a conservation equation: $\frac+ \nabla \cdot \mathbf F =,$ or simply in Einstein notation: $\frac+\frac= 0_i,$ where the conservation quantity

in this case is a vector, and

is a flux matrix. This can be simply proved.

Notes and References

Darrigol . O. . Frisch . U. . 2008 . From Newton's mechanics to Euler's equations . Physica D: Nonlinear Phenomena . en . 237 . 14–17 . 1855–1869 . 10.1016/j.physd.2007.08.003. 2008PhyD..237.1855D .
Elgindi . Tarek M. . 2021-11-01 . Finite-time singularity formation for $C^$ solutions to the incompressible Euler equations on $\mathbb^3$ . Annals of Mathematics . 194 . 1904.04795 . 3 . 10.4007/annals.2021.194.3.2 . 0003-486X.