See Theory for the Single-Phase Flow Interfaces and
Theory for the Turbulent Flow Interfaces in the
CFD Module User’s Guide for a description of the theory related to laminar and turbulent single-phase flow interfaces.
where Qted is the thermoelastic damping heat source (SI unit: W/(m
3)). This term is not included by default but must be added by selecting the corresponding check box.
where H is the enthalpy. The vector
where k is the turbulent kinetic energy, which in turn is defined by
The correlation between uj″ and
H″ in
Equation 4-140 is the turbulent transport of heat. It is modeled analogously to the laminar conductive heat flux
Taking the inner product between and Equation 4-147 results in an equation for the resolved kinetic energy, which can be subtracted from
Equation 4-146 with the following result:
According to Wilcox (Ref. 27), it is usually a good approximation to neglect the contributions of
k for flows with Mach numbers up to the supersonic range. This gives the following approximation of
Equation 4-148:
Larsson (Ref. 28) suggests to make the split
Equation 4-150 is completely analogous to the laminar energy equation of
Equation 4-15 and can be expanded using the same theory to get the temperature equation similar to
Equation 4-18 (see for example
Ref. 28):
This is a relatively exact model for PrT, while still quite simple. In
Ref. 29, it is compared to other models for
PrT and found to be a good approximation for most kinds of turbulent wall bounded flows except for turbulent flow of liquid metals. The model is given by
Weigand and others (Ref. 30) suggested an extension of
Equation 4-151 to liquid metals by introducing
where Re∞, the Reynolds number at infinity must be provided either as a constant or as a function of the flow field. This is entered in the Model Inputs section of the Fluid feature.
Analogous to the single-phase flow wall functions (see Wall Functions described for the Wall boundary condition), there is a theoretical gap between the solid wall and the computational domain for the fluid and temperature fields. This gap is often ignored when the computational geometry is drawn.
where ρ is the fluid density,
Cp is the fluid heat capacity, and
uτ is the friction velocity.
T+ is the dimensionless temperature and is given by (
Ref. 31):
λ is the thermal conductivity, and
κ is the von Kármán constant equal to 0.41.
The velocity at δw+c is then calculated from a velocity wall function, according to the turbulence model and wall treatment chosen to keep consistency. The velocity wall function used for the determination of
Uc when the L-VEL or a LES turbulence model is chosen is based on a blending of the linear and logarithmic laws. It gives very similar results to the wall function actually used for the calculation of the velocity at the first cell which is based on the wall Reynolds number.
The distance between the computational fluid domain and the wall, δw, is always
hw/2 for automatic wall treatment where
hw is the height of the mesh cell adjacent to the wall.
hw/2 is almost always very small compared to any geometrical quantity of interest, at least if a boundary layer mesh is used. For wall function,
δw is at least
hw/2 and can be bigger if necessary to keep
δw+ higher than
11.06. The computational results should be checked so that the distance between the computational fluid domain and the wall,
δw, is everywhere small compared to any geometrical quantity of interest. The distance
δw is available for evaluation on boundaries.