Large-Eddy Simulation of Shock-Turbulence Interaction Using the Approximate Deconvolution Model in a Finite Volume Scheme

Author(s):  
Roland von Kaenel ◽  
Leonhard Kleiser ◽  
Nikolaus Adams ◽  
Jan Vos
2003 ◽  
Vol 125 (2) ◽  
pp. 375-381 ◽  
Author(s):  
R. von Kaenel ◽  
N. A. Adams ◽  
L. Kleiser ◽  
J. B. Vos

The approximate deconvolution model for large-eddy simulation is formulated for a second-order finite volume scheme. With the approximate deconvolution model, an approximation of the unfiltered solution is obtained by repeated filtering, and given a good approximation of the unfiltered solution, the nonlinear terms of the Navier-Stokes equations are computed directly. The effect of scales not represented on the numerical grid is modeled by a relaxation regularization involving a secondary filter operation. A turbulent channel flow at a Mach number of M=1.5 and a Reynolds number based on bulk quantities of Re=3000 is selected for validation of the approximate deconvolution model implementation in a finite volume code. A direct numerical simulation of this configuration has been computed by Coleman et al. Overall, our large-eddy simulation results show good agreement with our filtered direct numerical simulation data. For this rather simple configuration and the low-order spatial discretization, differences between approximate deconvolution model and a no-model computation are found to be small.


2002 ◽  
Vol 124 (4) ◽  
pp. 829-835 ◽  
Author(s):  
R. von Kaenel ◽  
N. A. Adams ◽  
L. Kleiser ◽  
J. B. Vos

The approximate deconvolution model for large-eddy simulation is formulated for a second-order finite volume scheme. With the approximate deconvolution model, an approximation of the unfiltered solution is obtained by repeated filtering, and given a good approximation of the unfiltered solution, the nonlinear terms of the Navier-Stokes equations are computed directly. The effect of scales not represented on the numerical grid is modeled by a relaxation regularization involving a secondary filter operation. A turbulent channel flow at a Mach number of M=1.5 and a Reynolds number based on bulk quantities of Re=3000 is selected for validation of the approximate deconvolution model implementation in a finite volume code. A direct numerical simulation of this configuration has been computed by Coleman et al. Overall, our large-eddy simulation results show good agreement with our filtered direct numerical simulation data. For this rather simple configuration and the low-order spatial discretization, differences between approximate deconvolution model and a no-model computation are found to be small.


Author(s):  
Junli Bai ◽  
Ning Ma ◽  
Xiechong Gu

Freak wave is an unexpectedly large wave in ocean with extreme height and abnormal shape. The viscous effect is important in prediction of the flow patterns of the freak wave due to its strong nonlinearity. Thus, compared with the potential theory, more accurate information of flow field of the freak wave can be obtained by using the computational fluid dynamics (CFD) method. In laboratory test and numerical study, the focused wave is often adopted to substitute the freak wave in real sea. In this paper, we present a high accurate numerical model for large eddy simulation of the focused wave. In this model, the space filtered Navier-Stokes equations are solved on non-staggered grids by the finite volume. The fourth order compact scheme is adopted for discretization of both convection and diffusion terms of the governing equations. The standard fourth-order Runge-Kutta method is used for time advancement. The velocity-pressure coupling is ensured at each stage and the discretized equations are solved by strongly implicit procedure (SIP) method. The turbulence is simulated by the Smagorinsky model while the free surface is captured by using of the volume of fluid (VOF) method. The model is firstly validated by simulation of the cavity flow and linear wave. The simulation results are compared with theoretical values and published results, respectively. Finally, large eddy simulation of focused wave is presented. The comparison of the numerical results and measured data reveals that the proposed model is capable of reproducing the propagation and evolution of the focused wave.


2018 ◽  
Vol 35 (5) ◽  
pp. 2025-2049 ◽  
Author(s):  
Hongbin Liu ◽  
Xinrong Su ◽  
Xin Yuan

Purpose Adopting large eddy simulation (LES) to simulate the complex flow in turbomachinery is appropriate to overcome the limitation of current Reynolds-Averaged Navier–Stokes modelling and it provides a deeper understanding of the complicated transitional and turbulent flow mechanism; however, the large computational cost limits its application in high Reynolds number flow. This study aims to develop a three-dimensional GPU-enabled parallel-unstructured solver to speed up the high-fidelity LES simulation. Design/methodology/approach Compared to the central processing units (CPUs), graphics processing units (GPUs) can provide higher computational speed. This work aims to develop a three-dimensional GPU-enabled parallel-unstructured solver to speed up the high-fidelity LES simulation. A set of low-dissipation schemes designed for unstructured mesh is implemented with compute unified device architecture programming model. Several key parameters affecting the performance of the GPU code are discussed and further speed-up can be obtained by analysing the underlying finite volume-based numerical scheme. Findings The results show that an acceleration ratio of approximately 84 (on a single GPU) for double precision algorithm can be achieved with this unstructured GPU code. The transitional flow inside a compressor is simulated and the computational efficiency has been improved greatly. The transition process is discussed and the role of K-H instability playing in the transition mechanism is verified. Practical/implications The speed-up gained from GPU-enabled solver reaches 84 compared to original code running on CPU and the vast speed-up enables the fast-turnaround high-fidelity LES simulation. Originality/value The GPU-enabled flow solver is implemented and optimized according to the feature of finite volume scheme. The solving time is reduced remarkably and the detail structures including vortices are captured.


2010 ◽  
Vol 2 (2) ◽  
pp. 237-249 ◽  
Author(s):  
M. A. Uddin ◽  
C. Kato ◽  
N. Oshima ◽  
M. Tanahashi ◽  
T. Miyauchi

Large eddy simulation (LES) in homogeneous isotropic turbulence is performed by using the Finite element method (FEM) and Finite volume vethod (FVM) and the results are compared to show the performance of FEM and FVM numerical solvers. The validation tests are done by using the standard Smagorinsky model (SSM) and dynamic Smagorinsky model (DSM) for subgrid-scale modeling. LES is performed on a uniformly distributed 643 grids and the Reynolds number is low enough that the computational grid is capable of resolving all the turbulence scales. The LES results are compared with those from direct numerical simulation (DNS) which is calculated by a spectral method in order to assess its spectral accuracy. It is shown that the performance of FEM results is better than FVM results in this simulation. It is also shown that DSM performs better than SSM for both FEM and FVM simulations and it gives good agreement with DNS results in terms of both spatial spectra and decay of the turbulence statistics. Visualization of second invariant, Q, in LES data for both FEM and FVM reveals the existence of distinct, coherent, and tube-like vortical structures somewhat similar to those found in instantaneous flow field computed by the DNS. Keywords: Large eddy simulation; Validation; Smagorinsky model; Dynamic Smagorinsky model; Tube-like vortical structure; Homogeneous isotropic turbulence. © 2010 JSR Publications. ISSN: 2070-0237 (Print); 2070-0245 (Online). All rights reserved.DOI: 10.3329/jsr.v2i2.2582              J. Sci. Res. 2 (2), 237-249 (2010) 


Sign in / Sign up

Export Citation Format

Share Document