×

Improved discontinuity-capturing finite element techniques for reaction effects in turbulence computation. (English) Zbl 1177.76192

Summary: Recent advances in turbulence modeling brought more and more sophisticated turbulence closures (e.g., \(k-\varepsilon , k-\varepsilon -v ^{2}-f\), second moment closures), where the governing equations for the model parameters involve advection, diffusion and reaction terms. Numerical instabilities can be generated by the dominant advection or reaction terms. Classical stabilized formulations such as the Streamline-Upwind/Petrov-Galerkin (SUPG) formulation [A. N. Brook and T. J. R. Hughes, Comput. Methods Appl. Mech. Eng. 32, 199–259 (1982; Zbl 0497.76041); T. J. R. Hughes and T. E. Tezduyar, Comput. Methods Appl. Mech. Eng. 45, 217–284 (1984; Zbl 0542.76093)] are very well suited for preventing the numerical instabilities generated by the dominant advection terms. A different stabilization however is needed for instabilities due to the dominant reaction terms. An additional stabilization term, called the diffusion for reaction-dominated (DRD) term, was introduced by T. E. Tezduyar and Y. J. Park [Comput. Methods Appl. Mech. Eng. 59, 307–325 (1986; Zbl 0593.76096)] for that purpose and improves the SUPG performance. In recent years, a new class of variational multi-scale (VMS) stabilization [T. J. R. Hughes, Comput. Methods Appl. Mech. Eng. 127, No. 1–4, 387–401 (1995; Zbl 0866.76044)] has been introduced, and this approach, in principle, can deal with advection-diffusion-reaction equations. However, it was pointed out in [G. Hauke, Comput. Methods Appl. Mech. Eng. 191, No. 27–28, 2925–2947 (2002; Zbl 1005.76057)] that this class of methods also need some improvement in the presence of high reaction rates.
In this work, we show the benefits of using the DRD operator to enhance the core stabilization techniques such as the SUPG and VMS formulations. We also propose a new operator called the DRDJ (DRD with local variation jump) term, targeting the reduction of numerical oscillations in the presence of both high reaction rates and sharp solution gradients. The methods are evaluated in the context of two stabilized methods: the classical SUPG formulation and a recently-developed VMS formulation called the V-SGS [A. Corsini et al. Comput. Methods Appl. Mech. Eng. 194, No. 45–47, 4797–4823 (2005; Zbl 1093.76032)]. Model problems and industrial test cases are computed to show the potential of the proposed methods in the simulation of turbulent flows.

MSC:

76M10 Finite element methods applied to problems in fluid mechanics
76V05 Reaction effects in flows
76F99 Turbulence
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Hauke G (2002) A simple subgrid scale stabilized method for the advection-diffusion-reaction equation. Comput Methods Appl Mech Eng 191:2925–2947 · Zbl 1005.76057 · doi:10.1016/S0045-7825(02)00217-7
[2] Tezduyar TE, Park YJ (1986) Discontinuity capturing finite element formulations for nonlinear convection-diffusion-reaction equations. Comput Methods Appl Mech Eng 59:307–325 · Zbl 0593.76096 · doi:10.1016/0045-7825(86)90003-4
[3] Brooks AN, Hughes TJR (1982) Streamline upwind/Petrov–Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier–Stokes equations. Comput Methods Appl Mech Eng 32:199–259 · Zbl 0497.76041 · doi:10.1016/0045-7825(82)90071-8
[4] Corsini A, Rispoli F, Santoriello A (2005) A variational multiScale higher-order finite element formulation for turbomachinery flow computations. Comput Methods Appl Mech Eng 194:4797–4823 · Zbl 1093.76032 · doi:10.1016/j.cma.2004.11.013
[5] Lien FS, Kalitzin G, Durbin PA (1998) RANS modeling for compressible and transitional flows. Proceedings of the summer program 1998, Center for Turbulence Research, Stanford University, California pp 267–286
[6] Corsini A, Rispoli F, Santoriello A (2004) A new stabilized finite element method for advection-diffusion-reaction equations using quadratic elements. MFF – modelling fluid flow – the state of the art. Springer, Berlin Heidelberg New York, pp 247–266 · Zbl 1112.76385
[7] Tezduyar TE (2003) Computation of moving boundaries and interfaces and stabilization parameters. Int J Numer Methods Fluids 43:555–575 · Zbl 1032.76605 · doi:10.1002/fld.505
[8] Hughes TJR (1995) Multiscale phenomena: Green’s functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origin of stabilized methods. Comput Methods Appl Mech Eng 127:387–401 · Zbl 0866.76044 · doi:10.1016/0045-7825(95)00844-9
[9] Hughes TJR, Feijòo GR, Mazzei L, Quincy J-B (1998) The variational multiscale method-a paradigm for computational mechanics. Comput Methods Appl Mech Eng 166:3–24 · Zbl 1017.65525 · doi:10.1016/S0045-7825(98)00079-6
[10] Durbin PA (1995) Separated flow computations with the k- -v 2 model. AIAA J 33:659–664 · doi:10.2514/3.12628
[11] Elazar Y, Shreeve RP (1990) Viscous flow in a controlled diffusion compressor cascade with increasing incidence. ASME J Turbomach 112:256–266 · doi:10.1115/1.2927642
[12] Chen WL, Lien FS, Leschziner MA (1998) Computational prediction of flow around higly loaded compressor cascade blades with non-linear eddy-viscosity models. Int J Heat Fluid Flow 19:307–319 · doi:10.1016/S0142-727X(98)10001-2
[13] Tezduyar TE (1992) Stabilized finite element formulations for incompressible flow computations. Adv App Mech 28:1–44 · Zbl 0747.76069 · doi:10.1016/S0065-2156(08)70153-4
[14] Saad Y (1993) A flexible inner-outer preconditioned GMRES algorithm. SIAM J Sci Comput 14:461–469 · Zbl 0780.65022 · doi:10.1137/0914028
[15] Corsini A, Rispoli F, Santoriello A (2005) Quadratic Petrov-Galerkin finite elements for advective–reactive features in turbomachinery CFD. Int J Numer Methods Heat Fluid Flow 15(8):894–925 · doi:10.1108/09615530510625147
[16] Hughes TJR, Tezduyar TE (1984) Finite element methods for first-order hyperbolic systems with particular emphasis on the compressible Euler equations. Comput Methods Appl Mech Eng 45:217–284 · Zbl 0542.76093 · doi:10.1016/0045-7825(84)90157-9
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.