4.8.2. Transport Equations for the Intermittency Transition Model

The transport equation for intermittency is the following:

The formulation of the source terms is similar in nature to the corresponding terms in the - transition model. The transition source term is defined as:

where is the strain rate magnitude and

The destruction/relaminarization source is defined as follows:

where is the magnitude of the absolute vorticity rate, and . The transition onset is controlled by the following functions:

(4–196)

where is the wall distance and is a correlation used to trigger the transition model:

In this correlation, is the critical momentum thickness Reynolds number, and and are locally defined variables that approximate the freestream turbulence intensity and the pressure gradient parameter, respectively. A local formulation that provides levels in the middle of the boundary layer similar to the freestream turbulence intensity, , is expressed as:

The formula for reads as:

(4–197)

For numerical robustness, was bounded as follows:

(4–198)

The function accounts for the influence of the pressure gradient on transition, and is defined as:

The transition model has been calibrated against a wide range of generic as well as turbomachinery and external aeronautical test cases. While you should not have to adjust the model coefficients in most situations, the model provides access to the coefficients so that you can fine-tune the model. This should only be done based on detailed experimental data. You can adjust the following constants in the correlations: , , , , , , , , , .

The constant defines the minimal value of critical number, whereas the sum of and defines the maximal value of . The controls how fast the decreases as the turbulence intensity () increases. The and constants adjust the value of critical number in areas with favorable and adverse pressure gradient, respectively. The constant becomes active in regions with separation, correcting the value if necessary. Additionally, the size of the separation bubble can be tuned by another constant, , which enters the model through the source terms directly rather than the correlations. The constants and define the limits of the function; while they typically do not need to be adjusted, it is possible to do so (for details, contact your support engineer).

For crossflow transition, a widely used experimental correlation is the C1 correlation proposed by Arnal (1984) [223]. According to that correlation, transition occurs when the following condition is met:

where

The crossflow Reynolds number is defined as:

In the previous equations, is the crossflow Reynolds number, the index refers to a coordinate system that is aligned with the freestream velocity (where refers to the boundary layer edge), is the crossflow velocity normal to the direction of the freestream velocity. The value of 150 is the calibration constant of the Arnal C1 criterion corresponding to a critical Reynolds number.

The transition model accounts for the crossflow transition through a separate correlation, which approximates the C1 correlation through a local formulation. The formulation uses only local information and has the following functional form:

The function accounts for the influence of pressure gradient (shape factor); the function accounts for the ratio of the crossflow strength relative to the streamwise strength; the Reynolds number effect is included through . The correlation can be adjusted by the constant . The transition location due to the crossflow instability will move upstream as increases.

The function is computed as follows:

where is a local variable used to approximate the pressure gradient parameter and is defined similarly to (see Equation 4–197 and Equation 4–198):

The crossflow indicator function is a non-dimensional measure of the local crossflow strength relative to the streamwise strength. It is computed based on the wall normal change of the normalized vorticity vector , which serves as a measure of the three-dimensionality of the velocity profile and which becomes zero for 2D flow.

Here, is the vorticity vector and is the wall normal vector.

The indicator function is proportional to the local change of the flow angle and therefore a measure of the overall crossflow strength.

The triggering function that accounts for the crossflow effects is defined as follows:

The source term in the equation is then triggered, either by the primary triggering function (see Equation 4–196) of the Intermittency transition model, or by : whichever has a maximal value.

The factor of 100 is used in in order to speed up the transition process.

The default values for the model constants are the following:

= 100.0, = 1000.0, = 1.0, = 1.0, = -0.5, = 0.0, = 1.5, = 3.0, = 1.0, = 1.0