• Journal of Advanced Marine Engineering and Technology
• /
• v.34 no.8
• /
• pp.1084-1093
• /
• 2010

Large eddy simulation of turbulent premixed flame stabilized by the bluff body is performed by using sub-grid scale combustion model based on the G-equation describing the flame front propagation. The basic idea of LES modeling is to evaluate the filtered-front speed, which should be enhanced in the grid scale by the scale fluctuations. The dynamic subgrid scale models newly introduced into the G-equation are validated by the premixed combustion flow behind the triangle flame holder. The calculated results can predict the velocity and temperature of the combustion flow in good agreement with the experiment data.

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.3
• /
• pp.369-375
• /
• 2006

Large eddy simulation has been applied to simulate turbulent flow around a cubic obstacle mounted on a channel surface for a Reynolds number of 40000(based on the incoming bulk velocity and the obstacle height) using a Smagorinsky model and a Lagrangian dynamic model. In order to develop the LES to the practical engineering application, the effect of upwind scheme, turbulent sub-grid scale model were investigated. The computed velocities. turbulence quantifies, separation and reattachment length were evaluated by compared with the previous experimental results.

• Journal of Advanced Marine Engineering and Technology
• /
• v.29 no.8
• /
• pp.849-854
• /
• 2005

The laminar flame concept in turbulent reacting flow is considered applicable to many practical combustion systems For turbulent premixed combustion under widely used flamelet concept, the flame surface is described as an infinitely thin propagating surface that such a Propagating front can be represented as a level contour of a continuous function G. In this study, for the Purpose of validating the LES of G-equation combustion model. LES of turbulent Premixed combustion with dynamic SGS model of G-equation in turbulent channel flow are carried out A constant density assumption is used. The Predicted flame propagating speed is goof agreement with the DNS result of G. Bruneaux et al.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.

• Journal of the Society of Naval Architects of Korea
• /
• v.59 no.3
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• International Journal of Aeronautical and Space Sciences
• /
• v.16 no.2
• /
• pp.254-263
• /
• 2015

In this study, numerical simulations are carried out to investigate the dynamic SGS model effects in a Kerosene-LOx coaxial swirl injector under high pressure conditions. The turbulent model is based on large-eddy simulation (LES) with real-fluid transport and thermodynamics. To assess the effect of the dynamic subgrid-scale (SGS) model, the dynamic SGS model is compared with that of the algebraic SGS model. In a swirl injector under supercritical pressure, the characteristics of temporal pressure fluctuation and power spectral density (PSD) present comparable discrepancies dependant on the SGS models, which affect the mixing characteristics. Mixing efficiency and the probability density (PDF) function are conducted for a statistical description of the turbulent flow fields according to the SGS models. The back-scattering of turbulent kinetic energy is estimated in terms of the film thickness of the swirl injector.

• Wind and Structures
• /
• v.5 no.2_3_4
• /
• pp.369-378
• /
• 2002

The application of Large Eddy Simulation (LES) in a curvilinear coordinate system to the flow around a square cylinder is presented. In order to obtain sufficient resolution near the side of the cylinder, we use an O-type grid. Even with a curvilinear coordinate system, it is difficult to avoid the numerical oscillation arising in high-Reynolds-number flows past a bluff body, without using an extremely fine grid used. An upwind scheme has the effect of removing the numerical oscillations, but, it is accompanied by numerical dissipation that is a kind of an additional sub-grid scale effect. Firstly, we investigate the effect of numerical dissipation on the computational results in a case where turbulent dissipation is removed in order to clarify the differences between the effect of numerical dissipation. Next, the applicability and the limitations of the present method, which combine the dynamic SGS model with acceptable numerical dissipation, are discussed.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1240-1247
• /
• 2006

Large eddy simulation of turbulent premixed flame in turbulent channel flow is studied by using G-equation. A flamelet model for the premixed flame is combined with a dynamic subgrid combustion model for the filtered propagation flame speed. The objective of this work is to investigate the validity of the dynamic subgrid G-equation model to a complex turbulent premixed flame. The effect of model parameters of the dynamic sub grid G-equation on the turbulent flame speed is investigated. In order to consider quenching of laminar flames on the wall, wall-quenching damping function is employed in this calculation. In the present study, a constant density turbulent channel flow is used. The calculation results are evaluated by comparing with the DNS results of Bruneaux et al.

• Nuclear Engineering and Technology
• /
• v.41 no.7
• /
• pp.893-906
• /
• 2009

The present work presents the development of a Large Eddy Simulation (LES) methodology viable for complex geometries and suitable for the simulation of rod-bundles. The use of LES and Direct Numerical Simulation (DNS) allows for a deeper analysis of the flow field and the use of stochastical tools in order to obtain additional insight into rod-bundle hydrodynamics. Moreover, traditional steady-state CFD simulations fail to accurately predict distributions of velocity and temperature in rod-bundles when the pitch (P) to diameter (D) ratio P/D is smaller than 1.1 for triangular lattices of cylindrical pins. This deficiency is considered to be due to the failure to predict large-scale coherent structures in the region of the gap. The main features of the code include multi-block capability and the use of the fractional step algorithm. As a Sub-Grid-Scale (SGS) model, a Dynamic Smagorinsky model has been used. The code has been tested on plane channel flow and the flow in annular ducts. The results are in excellent agreement with experiments and previous calculations.

• Journal of computational fluids engineering
• /
• v.9 no.1
• /
• pp.1-9
• /
• 2004

Large eddy simulation is used to investigate the compressible flow over a cavity with high aspect ratio. The sub-grid scale stresses are modeled using the dynamic model. The compressible Navier-Stokes equations are solved with the sixth order accurate compact finite difference scheme in the space and the 4th order Runge-Kutta scheme in the time. The buffer Bone techniques are used for non-reflecting boundary conditions. The results show the shear layer oscillation over the cavity. The votical disturbances, the roll-up of vorticity, and impingement and scattering of vorticity at the downstream cavity edge can be seen in the shear layer. Several peaks for the resonant frequencies are found in the spectra of the vertical velocity at the center-line. The most energetic Peak near the downstream edge is different from that at the center part of the cavity The pressure has its minimum value in the vortex core inside the cavity, and becomes very high at the downstream face of the cavity. The variation of the model coefficient predicted by the dynamic model is quite large between 0 and 0.3. The model coefficient increases in the stream-wise evolution of the shear layer and sharply decreases near the wall due to the wall effect.