Explicit algebraic stress model explained

The algebraic stress model arises in computational fluid dynamics. Two main approaches can be undertaken. In the first, the transport of the turbulent stresses is assumed proportional to the turbulent kinetic energy; while in the second, convective and diffusive effects are assumed to be negligible. Algebraic stress models can only be used where convective and diffusive fluxes are negligible, i.e. source dominated flows. In order to simplify the existing EASM and to achieve an efficient numerical implementation the underlying tensor basis plays an important role. The five-term tensor basis that is introduced here tries to combine an optimum of accuracy of the complete basis with the advantages of a pure 2d concept. Therefore a suitable five-term basis is identified. Based on that the new model is designed and validated in combination with different eddy-viscosity type background models.

Integrity basis

In the frame work of single-point closures (Reynolds-stress transport models = RSTM) still provide the best representation of flow physics. Due to numeric requirements an explicit formulation based on a low number of tensors is desirable and was already introduced originally most explicit algebraic stress models are formulated using a 10-term basis:

b _{ij} = \sum ^{10} _{λ=1} G ^{(λ)} T ^{( λ )} _{ij} The reduction of the tensor basis however requires an enormous mathematical effort, to transform the algebraic stress formulation for a given linear algebraic RSTM into a given tensor basis by keeping all important properties of the underlying model. This transformation can be applied to an arbitrary tensor basis. In the present investigations an optimum set of basis tensors and the corresponding coefficients is to be found.

Projection method

The projection method was introduced to enable an approximate solution of the algebraic transport equation of the Reynolds-stresses. In contrast to the approach of the tensor basis is not inserted in the algebraic equation, instead the algebraic equation is projected. Therefore, the chosen basis tensors does not need to form a complete integrity basis. However, the projection will fail if the basis tensor are linear dependent. In the case of a complete basis the projection leads to the same solution as the direct insertion, otherwise an approximate solution in the sense is obtained.

An example

In order to prove, that the projection method will lead to the same solution as the direct insertion, the EASM for two-dimensional flows is derived. In two-dimensional flows only the tensors are independent.

T ^{( 1 )} _{ij} = s _{ij} T ^{( 2 )} _{ij} = s _{ik} w _{kj} - w _{ik} s _{kj} T ^{( 3 )} _{ij} = s _{ik} s _{kj} - s _{mk} s _{km} (1)/(3) \delta _{ij} The projection leads then to the same coefficients. This two-dimensional EASM is used as starting point for an optimized EASM which includes three-dimensional effects. For example the shear stress variation in a rotating pipe cannot be predicted with quadratic tensors. Hence, the EASM was extended with a cubic tensor. In order to do not affect the performance in 2D flows, a tensor was chosen that vanish in 2d flows. This offers the concentration of the coefficient determination in 3d flows. A cubic tensor, which vanishes in 3d flow is: T ^{( 5 )} _{ij} = w _{ik} s _{kl} s _{lj} - s _{ik} s _{kl} w _{lj} The projection with tensors T⁽¹⁾, T⁽²⁾, T⁽³⁾ and T⁽⁵⁾ yields then the coefficients of the EASM.

Limitation of C_μ

A direct result of the EASM derivation is a variable formulation of C_μ.As the generators of the extended EASM where chosen to preserve the existing 2D formulation the expression of C_μ remains unchanged:

C \mu = (- A _{1} g)/(g ^{2- (2)/(3)} A ^{2} _{3} η _{1} -2 A ^{2} η _{2}) A_i are the constants of the underlying pressure-strain model.Since η₁ is always positive it might be possible that C_μ becomes singular. Therefore in the first EASM derivation of a regularization was introduced, which prevent a singular by cutting the range of η₁. However, Wallin et al. pointed out that the regularization deteriorated the performance of the EASM. In their model the methodology was refined to account for the coefficient.g = C _{1} - 2 b _{ij}

This leads to a weak non-linear conditional equation for the EASM coefficients and an additional equation for g must be solved. In 3D the equation of g is of 6th order, wherefore a closed solution is only possible in 2D flows, where the equation reduces to 3rd order. In order to circumvent the root finding of a polynomial equation quasi self-consistent approach. He showed that by using a C_μ expression of a realizable linear model instead of the EASM-C_μ expression in the equation of g the same properties of g follows. For a wide range of and the quasi self-consistent approach is almost identical to the fully self-consistent solution. Thus the quality of the EASM is not affected with the advantage of no additional non-linear equation. Since in the projections to determine the EASM coefficients the complexity is reduced by neglecting higher order invariants.

References

1. Gatski, T.B. and Speziale, C.G., "On explicit algebraic stress models for complex turbulent flows". J. Fluid Mech.
2. Rung, T., "Entwicklung anisotroper Wirbelzähigkeitsbeziehungen mit Hilfe von Projektionstechniken", PHD-thesis, Technische Universität Berlin, 2000
3. Taulbee, D.B., "An improved algebraic Reynolds stress model and corresponding nonlinaer stress model", Phys. Fluids, 28, pp 2555–2561, 1992
4. Lübcke, H., Rung, T. and Thiele, F. "Prediction of the Spreading Mechanism of 3D Turbulent Wall Jets with Explicit Reynolds-Stress Closures", Eng. Turbulence Modelling and Experiments 5, Mallorca, 2002
5. Wallin, S. and Johansson, A.V., "A new explicit algebraic Reynolds stress turbulence model including an improved near-wall treatment", Flow Modelling and Turbulence Measurements IV
6. Taulbee, D.B., "An improved algebraic Reynolds stress model and corresponding nonlinear stress model"
7. Jongen, T. and Gatski, T.B., "General explicit algebraic stress relations and best approximations for three-dimensional flows", Int. J. Engineering Science