4.7.2. Transport Equations for the Transition SST Model

The transport equation for the intermittency is defined as:

(4–169)

The transition sources are defined as follows:

(4–170)

where is the strain rate magnitude, is an empirical correlation that controls the length of the transition region, and and hold the values of 2 and 1, respectively. The destruction/relaminarization sources are defined as follows:

(4–171)

where is the vorticity magnitude. The transition onset is controlled by the following functions:

(4–172)

(4–173)

(4–174)

where is the wall distance and is the critical Reynolds number where the intermittency first starts to increase in the boundary layer. This occurs upstream of the transition Reynolds number and the difference between the two must be obtained from an empirical correlation. Both the and correlations are functions of .

The constants for the intermittency equation are:

(4–175)

The transport equation for the transition momentum thickness Reynolds number is

(4–176)

The source term is defined as follows:

(4–177)

(4–178)

(4–179)

(4–180)

The model constants for the equation are:

(4–181)

The boundary condition for at a wall is zero flux. The boundary condition for at an inlet should be calculated from the empirical correlation based on the inlet turbulence intensity.

The model contains three empirical correlations. is the transition onset as observed in experiments. This has been modified from Menter et al.  [433] in order to improve the predictions for natural transition. It is used in Equation 4–176. is the length of the transition zone and is substituted in Equation 4–169. is the point where the model is activated in order to match both and , and is used in Equation 4–173. These empirical correlations are provided by Langtry and Menter [335].

(4–182)

The first empirical correlation is a function of the local turbulence intensity, :

(4–183)

where is the turbulent energy.

The Thwaites’ pressure gradient coefficient is defined as

(4–184)

where is the acceleration in the streamwise direction.

4.7.2.1. Separation-Induced Transition Correction

The modification for separation-induced transition is:

(4–185)

Here, is a constant with a value of 2.

The model constants in Equation 4–185 have been adjusted from those of Menter et al. [433] in order to improve the predictions of separated flow transition. The main difference is that the constant that controls the relation between and was changed from 2.193, its value for a Blasius boundary layer, to 3.235, the value at a separation point where the shape factor is 3.5  [433]. The boundary condition for at a wall is zero normal flux, while for an inlet, is equal to 1.0.

4.7.2.2. Coupling the Transition Model and SST Transport Equations

The transition model interacts with the SST turbulence model by modification of the -equation (Equation 4–101), as follows:

(4–186)

(4–187)

(4–188)

where and are the original production and destruction terms for the SST model. Note that the production term in the -equation is not modified. The rationale behind the above model formulation is given in detail in Menter et al.  [433].

In order to capture the laminar and transitional boundary layers correctly, the mesh must have a of approximately one. If the is too large (that is, > 5), then the transition onset location moves upstream with increasing . It is recommended that you use the bounded second order upwind based discretization for the mean flow, turbulence and transition equations.

4.7.2.3. Transition SST and Rough Walls

When the Transition SST Model is used together with rough walls, the roughness correlation must be enabled in the Viscous Model Dialog Box. This correlation requires the geometric roughness height as an input parameter, since, for the transition process from laminar to turbulent flow, the geometric roughness height is more important than the equivalent sand-grain roughness height . Guidance to determine the appropriate equivalent sand-grain roughness height (based on the geometric roughness height, shape and distribution of the roughness elements) can be obtained, for example, from Schlichting and Gersten [578] and Coleman et al. [117].

The roughness correlation is a modification of the built-in correlation for and is defined as:

(4–189)

The new defined is then used in the correlations for and . represents the transition momentum thickness Reynolds number. The value specified for the geometric roughness will apply to all walls. In case a different value is required for different walls, a user-defined function can be specified. It is important to note that the function for K is used in the volume (not at the wall). The function therefore must cover the region of the boundary layer and beyond where it should be applied. As an example, assume a roughness strip on a flat plate (x-streamwise direction, y wall normal z-spanwise) at the location (m) and a boundary layer thickness in that region of perhaps (spanwise extent (m)). The following pseudo-code will switch between roughness K0 everywhere (could be zero) and K1 at the transition strip. The height in the y-direction does not have to be exactly the boundary layer thickness – it can be much larger – as long as it does not impact other walls in the vicinity.

(4–190)