scholarly journals Deconvolution of induced spatial discretization filters subgrid modeling in LES: application to two-dimensional turbulence

2021 ◽  
Vol 2090 (1) ◽  
pp. 012064
Author(s):  
A Boguslawski ◽  
K Wawrzak ◽  
A Paluszewska ◽  
B J Geurts

Abstract The paper presents a new approximate deconvolution subgrid model for Large Eddy Simulation in which corrections to implicit filtering due to spatial discretization are integrated explicitly. The top-hat filter implied by second-order central finite differencing is a key example, which is discretised using the discrete Fourier transform involving all the mesh points in the computational domain. This discrete filter kernel is inverted by inverse Wiener filtering. The inverse filter obtained in this way is used to deconvolve the resolved scales of the implicitly filtered velocity field on the computational grid. Subgrid stresses are subsequently calculated directly from the deconvolved velocity field. The model was applied to study decaying two-dimensional turbulence. Results were compared with predictions based on the Smagorinsky model and the dynamic Germano model. A posteriori testing in which Large Eddy Simulation is compared with filtered Direct Numerical Simulation obtained with a Fourier spectral method is included. The new model presented strictly speaking applies to periodic problems. The idea of recovering a high-order inversion of the numerically induced filter kernel can be extended to more general non-periodic problems, also in three spatial dimensions.

Author(s):  
Xiaofeng Yang ◽  
Saurabh Gupta ◽  
Tang-Wei Kuo ◽  
Venkatesh Gopalakrishnan

A comparative cold flow analysis between Reynolds-averaged Navier–Stokes (RANS) and large eddy simulation (LES) cycle-averaged velocity and turbulence predictions is carried out for a single cylinder engine with a transparent combustion chamber (TCC) under motored conditions using high-speed particle image velocimetry (PIV) measurements as the reference data. Simulations are done using a commercial computationally fluid dynamics (CFD) code CONVERGE with the implementation of standard k-ε and RNG k-ε turbulent models for RANS and a one-equation eddy viscosity model for LES. The following aspects are analyzed in this study: The effects of computational domain geometry (with or without intake and exhaust plenums) on mean flow and turbulence predictions for both LES and RANS simulations. And comparison of LES versus RANS simulations in terms of their capability to predict mean flow and turbulence. Both RANS and LES full and partial geometry simulations are able to capture the overall mean flow trends qualitatively; but the intake jet structure, velocity magnitudes, turbulence magnitudes, and its distribution are more accurately predicted by LES full geometry simulations. The guideline therefore for CFD engineers is that RANS partial geometry simulations (computationally least expensive) with a RNG k-ε turbulent model and one cycle or more are good enough for capturing overall qualitative flow trends for the engineering applications. However, if one is interested in getting reasonably accurate estimates of velocity magnitudes, flow structures, turbulence magnitudes, and its distribution, they must resort to LES simulations. Furthermore, to get the most accurate turbulence distributions, one must consider running LES full geometry simulations.


2018 ◽  
Vol 180 ◽  
pp. 02054
Author(s):  
Martin Lasota ◽  
Petr Šidlof

The phonatory process occurs when air is expelled from the lungs through the glottis and the pressure drop causes flow-induced oscillations of the vocal folds. The flow fields created in phonation are highly unsteady and the coherent vortex structures are also generated. For accuracy it is essential to compute on humanlike computational domain and appropriate mathematical model. The work deals with numerical simulation of air flow within the space between plicae vocales and plicae vestibulares. In addition to the dynamic width of the rima glottidis, where the sound is generated, there are lateral ventriculus laryngis and sacculus laryngis included in the computational domain as well. The paper presents the results from OpenFOAM which are obtained with a large-eddy simulation using second-order finite volume discretization of incompressible Navier-Stokes equations. Large-eddy simulations with different subgrid scale models are executed on structured mesh. In these cases are used only the subgrid scale models which model turbulence via turbulent viscosity and Boussinesq approximation in subglottal and supraglottal area in larynx.


2011 ◽  
Vol 133 (2) ◽  
Author(s):  
Hassan Raiesi ◽  
Ugo Piomelli ◽  
Andrew Pollard

The performance of some commonly used eddy-viscosity turbulence models has been evaluated using direct numerical simulation (DNS) and large-eddy simulation (LES) data. Two configurations have been tested, a two-dimensional boundary layer undergoing pressure-driven separation, and a square duct. The DNS and LES were used to assess the k−ε, ζ−f, k−ω, and Spalart–Allmaras models. For the two-dimensional separated boundary layer, anisotropic effects are not significant and the eddy-viscosity assumption works well. However, the near-wall treatment used in k−ε models was found to have a critical effect on the predictive accuracy of the model (and, in particular, of separation and reattachment points). None of the wall treatments tested resulted in accurate prediction of the flow field. Better results were obtained with models that do not require special treatment in the inner layer (ζ−f, k−ω, and Spalart–Allmaras models). For the square duct calculation, only a nonlinear constitutive relation was found to be able to capture the secondary flow, giving results in agreement with the data. Linear models had significant error.


2019 ◽  
Vol 213 ◽  
pp. 02076
Author(s):  
Jan Sip ◽  
Frantisek Lizal ◽  
Jakub Elcner ◽  
Jan Pokorny ◽  
Miroslav Jicha

The velocity field in the area behind the automotive vent was measured by hot-wire anenemometry in detail and intensity of turbulence was calculated. Numerical simulation of the same flow field was performed using Computational fluid dynamics in commecial software STAR-CCM+. Several turbulence models were tested and compared with Large Eddy Simulation. The influence of turbulence model on the results of air flow from the vent was investigated. The comparison of simulations and experimental results showed that most precise prediction of flow field was provided by Spalart-Allmaras model. Large eddy simulation did not provide results in quality that would compensate for the increased computing cost.


2020 ◽  
Vol 19 (3-5) ◽  
pp. 207-239
Author(s):  
Saman Salehian ◽  
Reda R Mankbadi

The focus of this work is on understanding the effect of water injection from the launch pad on the noise generated during rocket’s lift-off. To simplify the problem, we consider a supersonic jet impinging on a flat plate with water injection from the impingement plate. The Volume of Fluid model is adopted in this work to simulate the two-phase flow. A Hybrid Large Eddy Simulation – Unsteady Reynolds Averaged Simulation approach is employed to model turbulence, wherein Unsteady Reynolds Averaged Simulation is used near the walls, and Large Eddy Simulation is used elsewhere in the computational domain. The numerical issues associated with simulating the noise of two-phase supersonic flow are addressed. The pressure fluctuations on the impingement plate obtained from numerical simulations agree well with the experimental data. Furthermore, the predicted effect of water injection on the far-field broadband noise is consistent with that of the experiment. The possible mechanisms for noise reduction by water injection are discussed.


Author(s):  
Sunil Patil ◽  
Danesh Tafti

Large eddy simulations of swirling flow and the associated convective heat transfer in a gas turbine can combustor under cold flow conditions for Reynolds numbers of 50,000 and 80,000 with a characteristic Swirl number of 0.7 are carried out. A precursor Reynolds averaged Navier-Stokes (RANS) simulation is used to provide the inlet boundary conditions to the large-eddy simulation (LES) computational domain, which includes only the can combustor. A stochastic procedure based on the classical view of turbulence as a superposition of the coherent structures is used to simulate the turbulence at the inlet plane of the computational domain using the mean flow velocity and Reynolds stress data from the precursor RANS simulation. To further reduce the overall computational resource requirement and the total computational time, the near wall region is modeled using a zonal two layer model (WMLES). A novel formulation in the generalized co-ordinate system is used for the solution of effective tangential velocity and temperature in the inner layer virtual mesh. The WMLES predictions are compared with the experimental data of Patil et al. (2011, “Experimental and Numerical Investigation of Convective Heat Transfer in Gas Turbine Can Combustor,” ASME J. Turbomach., 133(1), p. 011028) for the local heat transfer distribution on the combustor liner wall obtained using robust infrared thermography technique. The heat transfer coefficient distribution on the liner wall predicted from the WMLES is in good agreement with experimental values. The location and the magnitude of the peak heat transfer are predicted in very close agreement with the experiments.


2020 ◽  
Vol 6 (4) ◽  
Author(s):  
N. Satish ◽  
K. Venkatasubbaiah

Abstract Pulsation and acceleration of liquid metal turbulent flow through a horizontal channel has been numerically studied using a large eddy simulation (LES) technique. The effect of inlet pulsation and acceleration on flow and heat transfer characteristics of low Prandtl number liquid metal flow have been investigated and reported here. Results have been presented for different Reynolds numbers, different amplitudes, and frequency with constant bottom wall thickness. The flow field is modeled as unsteady-state two-dimensional incompressible turbulent-forced convection flow. Turbulence is modeled using a LES technique. Two-dimensional unsteady-state heat conduction equation is solved to know the temperature distribution in the solid region. Finite difference method solver is developed for solving the governing equations using sixth-order accuracy of compact schemes. The average Nusselt number shows cyclic variation with respect to time in pulsation flows. The enhancement of heat transfer with pulsation at amplitude 0.4 and frequency 100 Hz is 6.51%. The rate of heat transfer increases in pulsation flow compared to quasi-steady flow. The inlet acceleration shows a significant effect on flow characteristics. The present results are compared with direct numerical simulation (DNS) results available in the literature and matching well with DNS data.


Author(s):  
Ying Huai ◽  
Amsini Sadiki

In this work, Large Eddy Simulation (LES) has been carried out to analyze the turbulent mixing processes in an impinging jet configuration. To characterize and quantify turbulent mixing processes, in terms of scalar structures and degree of mixing, three parameters have been basically introduced. They are “mixedness parameter”, which represents the probability of mixed fluids in computational domain, the Spatial Mixing Deficiency (SMD) and the Temporal Mixing Deficiency (TMD) parameters for characterizing the mixing at different scalar scale degrees. With help of these parameters, a general mixing optimization procedure has then been suggested and achieved in an impinging jet configuration. An optimal jet angle was estimated and the overall mixing degree with this jet angle reached around six times more than the original design. It turns out that the proposed idea and methodology can be helpful for practical engineering design processes.


Sign in / Sign up

Export Citation Format

Share Document