A low-dissipation DG method for the under-resolved simulation of low Mach number turbulent flows

The low-frequency character of two model problems is exploited in order to illustrate the acoustic consequences of the interactions between chemically reacting (or relaxing) inhomogeneities and flames or constrictions in ducts. The monopole of the former is associated with heat transfer in a fluid which exhibits variations in its specific heats, while in the latter there is an extension of the classical phenomenon associated with the pulsations of an inhomogeneity of the fluid compressibility. This second mechanism is found to be insignificant, but the heat-conduction source is considered to be very powerful at sufficiently low Mach numbers; in fact, to first order it is independent of the flow Mach number for laminar, as well as a certain class of turbulent, flows.

Download Full-text

A Time-Accurate Projection Method for Low-Mach Number Variable Density Flows in Curvilinear Grids

Volume 1: Symposia, Parts A, B and C ◽

10.1115/fedsm2009-78327 ◽

2009 ◽

Author(s):

F. N. Fard ◽

B. Lessani

Keyword(s):

Mach Number ◽

Turbulent Flows ◽

Projection Method ◽

Time Integration ◽

Variable Density ◽

Integration Scheme ◽

Coordinate Systems ◽

Low Mach Number ◽

Curvilinear Coordinate ◽

Time Integration Scheme

A time-accurate numerical algorithm is proposed for low Mach number variable density flows in curvilinear coordinate systems. In order to increase the stability of the method, a predictor-corrector time integration scheme, coupled with the projection method, is employed. The projection method results in a constant-coefficient Poisson equation for the pressure in both the predictor and corrector steps. The continuity equation is fully satisfied at each step. To prevent the pressure odd-even decoupling typically encountered in collocated grids, a flux interpolation technique is developed. The spatial discretization method offers computational simplicity and straightforward extension to 3D curvilinear coordinate systems, which are essential in the simulation of turbulent flows in complex geometries. The accuracy and stability of the algorithm are tested with a series of numerical experiments, and the results are validated against the available data in the literature.

Download Full-text

Large Eddy Simulations of Low-Mach Number Aeroacoustics for Complex Wall-Bounded Turbulent Flows via a Multiblock/Characteristic Interface Approach

12th AIAA/CEAS Aeroacoustics Conference (27th AIAA Aeroacoustics Conference) ◽

10.2514/6.2006-2633 ◽

2006 ◽

Cited By ~ 1

Author(s):

Jungsoo Suh ◽

Steven Frankel

Keyword(s):

Mach Number ◽

Turbulent Flows ◽

Large Eddy Simulations ◽

Low Mach Number ◽

Large Eddy

Download Full-text

Artificial Compressibility Method and Preconditioning Method for Solving Two Dimensional Incompressibile Flow

ASME-JSME-KSME 2011 Joint Fluids Engineering Conference: Volume 1, Symposia – Parts A, B, C, and D ◽

10.1115/ajk2011-01007 ◽

2011 ◽

Cited By ~ 1

Author(s):

Hyungro Lee ◽

Einkeun Kwak ◽

Seungsoo Lee

Keyword(s):

Mach Number ◽

Turbulent Flows ◽

Stokes Equations ◽

Navier Stokes ◽

Two Dimensional ◽

Navier Stokes Equations ◽

Low Mach Number ◽

Artificial Compressibility ◽

Artificial Compressibility Method ◽

Preconditioning Method

In this study, two commonly used numerical methods for the analysis of incompressible flows (or low Mach number flows), Chorins’ artificial compressibility method and Wiess and Smith’s preconditioning method are compared. Also, the convergence characteristics of two methods are numerically investigated for two-dimensional laminar and turbulent flows. Although the two methods have similar governing equations, the eigensystems and other details are very different. The eigensystems of the artificial compressibility method and the preconditioning method are analytically examined. An artificial compressibility code that solves the incompressible RANS (Reynolds Averaged Navier-Stokes) equations is newly developed for the study. An artificial compressibility code and a well-verified existing low Mach number code uses Roe’s approximate Riemann solver in conjunction with a cell centered finite volume method. Using MUSCL extrapolation with nonlinear limiters, 2nd order spatial accuracy is achieved while maintaining TVD (total variation diminishing) property. AF-ADI (approximate factorization-alternate direction implicit) method is used to get the steady solution for both codes. Menter’s k–ω SST turbulence model is used for the analysis of turbulent flows. Navier-Stokes equations and the turbulence model equations are solved in a loosely coupled manner.

Download Full-text

LES and DNS of Turbulent Flows in Weakly Compressible Condition

Volume 1: Symposia, Parts A and B ◽

10.1115/fedsm2007-37300 ◽

2007 ◽

Author(s):

Takeo Kajishima ◽

Takashi Ohta

Keyword(s):

Mach Number ◽

Turbulent Flow ◽

Turbulent Flows ◽

Plane Channel ◽

Scale Model ◽

Turbulent Shear ◽

Quadrant Analysis ◽

Number Range ◽

Low Mach Number ◽

Flow Scheme

Flow field of low Mach number (e.g. M <0.3) is usually simulated by the incompressible flow scheme due to the severe limitation of time-increment in the compressible flow scheme. In this work, we propose a modification to the usual incompressible scheme, based on the elliptic equation for pressure, to improve the accuracy for turbulent flows considering weak compressibility. Two examples will be shown to validate our method. (1) LES (Large-Eddy Simulation) was conducted for turbulent flow around NACA0012 airfoil. Particular attention was focused on the influence of compressibility, despite the low Mach number range. In addition, new subgrid scale model of one-equation type using dynamic procedure was compared with traditional Smagorinsky model. Our method successfully reproduced the separation bubble near the leading edge, resulted in the improvement in the intensity of pressure fluctuation. (2) DNS (Direct Numerical Simulation) of turbulent flow in a plane channel is carried out, taking wall temperature difference into account. As a result of the density fluctuation in near-wall eddies, asymmetric profiles are observed in turbulence statistics. By the 4-quadrant analysis of turbulent shear stress, it is found that the ejection events in the vicinity of the walls are particularly affected by the density variation.

Download Full-text