The Low Reynolds Number k-ε Turbulence Model
When the accuracy provided by wall functions in the k-ε model is not enough, a so called low Reynolds number model can be used. “Low Reynolds number” refers to the region close to the wall where viscous effects dominate.
Most low Reynolds number k-ε models adapt the turbulence transport equations by introducing damping functions. This module includes the AKN model (after the inventors Abe, Kondoh, and Nagano; Ref. 12). The AKN k-ε model for compressible flows reads (Ref. 8 and Ref. 12):
(3-121)
where
(3-122)
and
(3-123)
lw is the distance to the closest wall.
Realizability Constraints are applied to the low Reynolds number k-ε model.
Wall Distance
The wall distance variable, lw, is provided by a mathematical Wall Distance interface that is included when using the low Reynolds number k-ε model. The solution to the wall distance equation is controlled using the parameter lref. The distance to objects larger than lref is represented accurately, while objects smaller than lref are effectively diminished by appearing to be farther away than they actually are. This is a desirable feature in turbulence modeling since small objects would get too large an impact on the solution if the wall distance were measured exactly.
The most convenient way to handle the wall distance variable is to solve for it in a separate study step. A Wall Distance Initialization study type is provided for this purpose and should be added before the actual Stationary or Time Dependent study step.
Wall Boundary Conditions
Automatic Wall Treatment
The automatic wall treatment is a way to obtain an accurate low-Reynolds-number formulation when the mesh allows it, and to fall back on a wall function formulation when the mesh is coarse. It is a blending between the solutions in the linear sublayer and the logarithmic layer respectively. A formulation for ω-based methods was described by Menter et al. in Ref. 9. The Low-Reynolds number model uses a similar formulation and defines
(3-124)
and
(3-125)
δw is the distance to the closest wall. The boundary conditions for the momentum equations are a no-penetration condition u n = 0 and a shear stress condition
(3-126)
In Equation 3-126, u+ = U||/uτ with
(3-127)
where in turn
(3-128)
Here, κv, is the von Kárman constant (default value 0.41), U|| is the velocity parallel to the wall and B is a constant that by default is set to 5.2.
These expression can be combined with the lift-off concept shown in Figure 3-7 which gives δw = hw/2. The k-equation formally fulfills both at the wall and in the log-layer, so this condition is applied for all δw+.
The conditions for the turbulent dissipation, ε, is given by the Wolfshtein model which is commonly employed in two-layer k-ε implementations (Ref. 13):
(3-129)
The resulting wall resolution, δw+, is available as the postprocessing variable. Delta_wPlus.
Low Reynolds Number Wall Treatment
The damping terms in the equations for k and ε allow for a no slip condition to be applied to the velocity, that is u=0 which is the case when Wall Treatment is set to Low Re.
Since all velocities must disappear on the wall, so must k. Hence, k=0 on the wall.
The correct wall boundary condition for ε is
where n is the wall normal direction. This condition is however numerically very unstable. Therefore, ε is not solved for in the cells adjacent to a solid wall and instead the analytical relation
(3-130)
is prescribed in those cells (using the variable εw, which only exists in those cells). Equation 3-130 can be derived as the first term in a series expansion of
For the expansion to be valid, it is required that
is the distance, measured in viscous units, from the wall to the center of the wall adjacent cell. The boundary variable Distance to cell center in viscous units, lplus_cc, is available to ensure that the mesh is fine enough. Observe that it is unlikely that a solution is obtained at all if
Inlet Values for the Turbulence Length Scale and Intensity
The guidelines given in Inlet Values for the Turbulence Length Scale and Turbulent Intensity for selecting the turbulence length scale, LT, and the turbulence intensity, IT, apply also to the low-Reynolds number k-ε model.
Initial Values
The low-Reynolds number k-ε model has the same default initial guess as the standard k-ε model (see Initial Values) but with replaced by lref.
The default initial value for the wall distance equation (which solves for the reciprocal wall distance) is 2/lref.
In some cases, especially for stationary solutions, a fast way to convergence is to first solve the model using the ordinary k-ε model and then to use that solution as an initial guess for the low-Reynolds number k-ε model. The procedure is then as follows:
1
2
3
Add a new Stationary with Initialization study.
4
In the Wall Distance Initialization study step, set Values of variables not solved for to Solution from the first study. This is to propagate the solution from the first study down to the second step in the new study.
5
Scaling for Time-Dependent Simulations
The low-Reynolds number k-ε model applies absolute scales of the same type as the k-ε model (see Scaling for Time-Dependent Simulations).
In the COMSOL Multiphysics Reference Manual: