High Fidelity Simulation of Turbulent Jet and Identification of Acoustic Sources

A turbulent round jet at the Reynolds number of 7,290 is simulated by direct numerical simulation to study jet noise source and its propagation. A large domain of 60 diameters in the axial direction, 30 diameters in the radial direction, and full 360 degrees is chosen to minimize the effect of boundary conditions. The spatial statistics of the simulated flow and the peak amplitude of acoustic radiation at 45-degree angle agree well with existing literature, confirming the quality of the data. Motivated by the Ffowcs-Williams and Hawkings equation, acoustic spectra at three representative locations of the potential core, the shear layer, and the ambient indicate that the fluctuation magnitudes are highest in the shear layer, followed by the potential core and the ambient. The turbulent characteristics are isotropic in the potential core and the shear layer, and the ambient exhibits an intermittent behaviour

Numerically stable formulations of convective terms for turbulent compressible flows

A systematic analysis of the discrete conservation properties of non-dissipative, central-difference approximations of the compressible Navier-Stokes equations is reported. A general triple splitting of the nonlinear convective terms is considered, and energy-preserving formulations are fully characterized by deriving a two-parameter family of split forms. Previously developed formulations reported in literature are shown to be particular members of this family; novel splittings are introduced and discussed as well. Furthermore, the conservation properties yielded by different choices for the energy equation (i.e. total and internal energy, entropy) are analyzed thoroughly. It is shown that additional preserved quantities can be obtained through a suitable adaptive selection of the split form within the derived family. Local conservation of primary invariants, which is a fundamental property to build high-fidelity shock-capturing methods, is also discussed in the paper. Numerical tests performed for the Taylor-Green Vortex at zero viscosity fully confirm the theoretical findings, and show that a careful choice of both the splitting and the energy formulation can provide remarkably robust and accurate results

Numerically stable formulations of convective terms for turbulent compressible flows

A systematic analysis of the discrete conservation properties of non-dissipative, central-difference approximations of the compressible Navier–Stokes equations is reported. A generalized splitting of the nonlinear convective terms is considered, and energy-preserving formulations are fully characterized by deriving a two-parameter family of split forms. Previously developed formulations reported in literature are shown to be particular members of this family; novel splittings are introduced and discussed as well. Furthermore, the conservation properties yielded by different choices for the energy equation (i.e. total and internal energy, entropy) are analyzed thoroughly. It is shown that additional preserved quantities can be obtained through a suitable adaptive selection of the split form within the derived family. Local conservation of primary invariants, which is a fundamental property to build high-fidelity shock-capturing methods, is also discussed in the paper. Numerical tests performed for the Taylor–Green Vortex at zero viscosity fully confirm the theoretical findings, and show that a careful choice of both the splitting and the energy formulation can provide remarkably robust and accurate results.Peer ReviewedPostprint (author's final draft

Investigation of low-dissipation monotonicity-preserving scheme for direct numerical simulation of compressible turbulent flows

© 2014 Elsevier Ltd. The influence of numerical dissipation on direct numerical simulation (DNS) of decaying isotropic turbulence and turbulent channel flow is investigated respectively by using the seventh-order low-dissipation monotonicity-preserving (MP7-LD) scheme with different levels of bandwidth dissipation. It is found that for both benchmark test cases, small-scale turbulence fluctuations can be largely suppressed by high level of scheme dissipation, while the appearance of numerical errors in terms of high-frequency oscillations could destabilize the computation if the dissipation is reduced to a very low level. Numerical studies show that reducing the bandwidth dissipation to 70% of the conventional seventh-order upwind scheme can maximize the efficiency of the MP7-LD scheme in resolving small-scale turbulence fluctuations and, in the meantime preventing the accumulation of non-physical numerical errors. By using the optimized MP7-LD scheme, DNS of an impinging oblique shock-wave interacting with a spatially-developing turbulent boundary layer is conducted at an incoming free-stream Mach number of 2.25 and the shock angle of 33.2°. Simulation results of mean velocity profiles, wall surface pressure, skin friction and Reynolds stresses are validated against available experimental data and other DNS predictions in both the undisturbed equilibrium boundary layer region and the interaction zone, and good agreements are achieved. The turbulence kinetic energy transport equation is also analyzed and the balance of the equation is well preserved in the interaction region. This study demonstrates the capability of the optimized MP7-LD scheme for DNS of complex flow problems of wall-bounded turbulent flow interacting with shock-waves

An Entropy Stable Nodal Discontinuous Galerkin Method for the Two Dimensional Shallow Water Equations on Unstructured Curvilinear Meshes with Discontinuous Bathymetry

We design an arbitrary high-order accurate nodal discontinuous Galerkin spectral element approximation for the nonlinear two dimensional shallow water equations with non-constant, possibly discontinuous, bathymetry on unstructured, possibly curved, quadrilateral meshes. The scheme is derived from an equivalent flux differencing formulation of the split form of the equations. We prove that this discretisation exactly preserves the local mass and momentum. Furthermore, combined with a special numerical interface flux function, the method exactly preserves the mathematical entropy, which is the total energy for the shallow water equations. By adding a specific form of interface dissipation to the baseline entropy conserving scheme we create a provably entropy stable scheme. That is, the numerical scheme discretely satisfies the second law of thermodynamics. Finally, with a particular discretisation of the bathymetry source term we prove that the numerical approximation is well-balanced. We provide numerical examples that verify the theoretical findings and furthermore provide an application of the scheme for a partial break of a curved dam test problem