General Single-Phase Flow Theory
The Single-Phase Fluid Flow interfaces are based on the Navier-Stokes equations, which in their most general form read
(3-9)
(3-10)
(3-11)
where
•
ρ is the density (SI unit: kg/m3)
•
u is the velocity vector (SI unit: m/s)
•
p is pressure (SI unit: Pa)
•
τ is the viscous stress tensor (SI unit: Pa)
•
F is the volume force vector (SI unit: N/m3)
•
Cp is the specific heat capacity at constant pressure (SI unit: J/(kg·K))
•
T is the absolute temperature (SI unit: K)
•
q is the heat flux vector (SI unit: W/m2)
•
Q contains the heat sources (SI unit: W/m3)
•
S is the strain-rate tensor:
The operation “:” denotes a contraction between tensors defined by
(3-12)
This is sometimes referred to as the double dot product.
Equation 3-9 is the continuity equation and represents conservation of mass. Equation 3-10 is a vector equation which represents conservation of momentum. Equation 3-11 describes the conservation of energy, formulated in terms of temperature. This is an intuitive formulation that facilitates boundary condition specifications.
To close the equation system, Equation 3-9 through Equation 3-11, constitutive relations are needed. For a Newtonian fluid, which has a linear relationship between stress and strain, Stokes (Ref. 1) deduced the following expression:
(3-13)
The dynamic viscosity, μ (SI unit: Pa·s), for a Newtonian fluid is allowed to depend on the thermodynamic state but not on the velocity field. All gases and many liquids can be considered Newtonian. Examples of non-Newtonian fluids are honey, mud, blood, liquid metals, and most polymer solutions. With the CFD Module, you can model flows of non-Newtonian fluids using the predefined models. The following models which describe the dynamic viscosity for non-Newtonian fluids are avalible: power law, Carreau, Bingham-Papanastasiou, Herschel-Bukley-Papanastasiou, and Casson-Papanastasiou models. The Heat Transfer Module treats all fluids as Newtonian according to Equation 3-13.Other commonly used constitutive relations are Fourier’s law of heat conduction and the ideal gas law.
In theory, the same equations describe both laminar and turbulent flows. In practice, however, the mesh resolution required to simulate turbulence with the Laminar Flow interface makes such an approach impractical.
There are several books where derivations of the Navier-Stokes equations and detailed explanations of concepts such as Newtonian fluids can be found. See, for example, the classical text by Batchelor (Ref. 3) and the more recent work by Panton (Ref. 4).
Many applications describe isothermal flows for which Equation 3-11 is decoupled from Equation 3-9 and Equation 3-10. Nonisothermal flow and the temperature equation are described in the Heat Transfer and Nonisothermal Flow Interfaces chapter.
2D Axisymmetric Formulations
A 2D axisymmetric formulation of Equation 3-9 and Equation 3-10 requires to be zero. That is, there must be no gradients in the azimuthal direction. A common additional assumption is, however, that . In such cases, the -equation can be removed from Equation 3-10. The resulting system of equations is both easier to converge and computationally less expensive compared to retaining the -equation. The default 2D axisymmetric formulation of Equation 3-9 and Equation 3-10 therefore assumes that
You can activate the Swirl Flow property which reduces the above assumptions to and reintroduces the -equation into Equation 3-10.