• 대한기계학회논문집
• /
• 제19권7호
• /
• pp.1691-1701
• /
• 1995

Two aspects of Large-Eddy Simulation(LES) are investigated in order to improve its performance. The first one is on how to determine the model coefficient in conjunction with a dynamic subgrid-scale model, and the second one is on a wall-layer model(WLM) which allows one to skip near-wall regions to save a large number of grid points otherwise required. Especially, a WLM suitable for a separated flow is considered. Firstly, an averaging technique to calculate the model coefficient of dynamic subgrid-scale modeling(DSGSM) is introduced. The technique is based on the concept of local averaging, and useful to stabilize numerical solution in conjunction with LES of complex turbulent flows using DSGSM. It is relatively simple to implement, and takes very low overhead in CPU time. It is also able to detect the region of negative model coefficient where the "backscattering" of turbulence energy occurs. Secondly, a wall-layer model based on a local turbulence intensity is considered. It locally determines wall-shear stresses depending on the local flow situations including separation, and yields better predictions in separated regions than the conventional WLM. The two techniques are tested for a turbulent obstacle flow, and show the direction of further improvements.rovements.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• pp.819-822
• /
• 2002

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2011년 춘계학술대회논문집
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 추계 학술대회논문집
• /
• pp.221-225
• /
• 2004

Evaluation of turbulence models is performed for a better prediction of thermal stratification in an upper plenum of a liquid metal reactor by applying them to the experiment conducted at JNC. The turbulence models tested in the present study are the two-layer model, the $\kappa-\omega$ model, the v2-f model and the low-Reynolds number differential stress-flux model. When the algebraic flux model or differential flux model are used for treating the turbulent heat flux, there exist little differences between turbulence models in predicting the temporal variation of temperature. However, the v2-f model and the low-Reynolds number differential stress-flux model better predict the steep gradient o( temperature at the interface of thermal stratification, and only the v2-f model predicts properly the oscillation of temperature. The LES Is needed for a better prediction of the amplitude and frequency of the temperature fluctuation.

• 한국전산유체공학회지
• /
• 제13권3호
• /
• pp.28-34
• /
• 2008

In engineering application of large eddy simulation, there are still questions as follows grid dependency on numerical results, the effect of upwind scheme against a calculation instability, appropriate boundary conditions dealing with turbulence fluctuation and the performance of SGS models. In this study, in order to develop the LES to the engineering application, large eddy simulation was carried out to investigate the effect of upwind scheme, turbulent subgrid model and the grid dependancy of the flow around a wall-mounted cube in a channel at Re=40,000 based on cubic height and bulk mean velocity. The computed velocities, turbulence quantities, separation and reattachment length were evaluated compared with the experimental results of R. Matinuzzi and C. Tropea.

• International Journal of Fluid Machinery and Systems
• /
• 제2권4호
• /
• pp.346-352
• /
• 2009

The objective of the present work is to improve numerical predictions of unsteady turbulent swirling flows in the draft tubes of hydraulic power plants. We present Large Eddy Simulation (LES) results on a simplified draft tube consisting of a straight conical diffuser. The basis of LES is to solve the large scales of motion, which contain most of the energy, while the small scales are modeled. LES strategy is here preferred to the average equations strategies (RANS models) because it resolves directly the most energetic part of the turbulent flow. LES is now recognized as a powerful tool to simulate real applications in several engineering fields which are more and more frequently found. However, the cost of large-eddy simulations of wall bounded flows is still expensive. Bypass methods are investigated to perform high-Reynolds-number LES at a reasonable cost. In this study, computations at a Reynolds number about 2 $10^5$ are presented. This study presents the result of a new near-wall model for turbulent boundary layer taking into account the streamwise pressure gradient (adverse or favorable). Validations are made based on simple channel flow, without any pressure gradient and on the data base ERCOFTAC. The experiments carried out by Clausen et al. [1] reproduce the essential features of the complex flow and are used to develop and test closure models for such flows.

• 한국수자원학회논문집
• /
• 제52권10호
• /
• pp.697-706
• /
• 2019

다상흐름 모델링 기법과 하이브리드 난류 모델링 기법을 결합한 수치모형을 이용하여 사각형 수로에서의 중력류를 수치모의 하였다. 이 연구에서 적용한 다상흐름 해석기법은 밀도가 큰 중력류 유체, 상대적으로 밀도가 작은 주변류 유체 그리고 자유수면 위에서 흐르는 공기를 3개의 상으로 처리하며, 각 상에 대해서 분리된 흐름 지배방정식을 적용한다. 난류흐름은 벽경계 근처에서는 RANS 모드로 모의하고 벽에서 떨어진 영역에서는 LES 모드로 해석하는 하이브리드 RANS/LES 방법의 일종인 IDDES 기법을 이용하여 해석한다. 이 연구에서 적용한 모델링 기법은 중력류의 머리의 전파속도를 실험값과 일치하게 잘 예측하는 것으로 나타났다. 수치해석 결과는 아울러 낮은 레이놀즈수 난류모형을 이용한 RANS 수치모의에서 이용되는 정도의 격자해상도에서도 큰 규모의 Kelvin-Helmholtz 형식의 경계면 와의 발달과 이들 와가 지속적으로 3차원 형식의 붕괴를 거쳐 작은 난류구조로 분해되면서 난류에너지가 소산되는 현상을 성공적으로 예측함을 보여준다. 적용한 수치모의 기법은 공학적으로 접근 가능한 격자해상도에서 돌출-쪼개짐 흐름 불안정을 동반한 중력류 머리부분의 3차원 거동 특성을 잘 재현하며, 이 결과는 보다 높은 격자해상도에서 구해진 LES 결과에 상응하는 것으로 나타났다. 이 연구결과는 하이브리드 난류모델링 기법과 다상흐름 해석기법을 병합한 수치모형이 자연상태에서 복잡한 중력류의 물리적 거동을 예측하는데 공학적으로 유망한 방법임을 보여준다.

• 한국전산유체공학회지
• /
• 제9권2호
• /
• pp.23-29
• /
• 2004

A numerical simulation of an incompressible cavity flow is conducted using the hybrid turbulence model. The model adopted is a modified type of DES using k- ε two-equation model. Cavity geometry and flow condition are based on Cattafesta's experiment. Computational results are compared with the results of Cattafesta's experiment. The simulation successfully predicts the oscillatory features and the Strouhal number of the oscillation compares very favorably with that of the dominant mode of experimental data. Vorticity contours obtained from the simulation data are consistent with the smoke visualization of the Cattafesta's experiment. The coherent structures of cavity flow are also investigated using Q criterion.

• Wind and Structures
• /
• 제13권6호
• /
• pp.557-580
• /
• 2010

A new inflow turbulence generation method and a combined dynamic SGS model recently developed by the authors were applied to evaluate the wind effects on 508 m high Taipei 101 Tower. Unlike the majority of the past studies on large eddy simulation (LES) of wind effects on tall buildings, the present numerical simulations were conducted for the full-scale tall building with Reynolds number greater than $10^8$. The inflow turbulent flow field was generated based on the new method called discretizing and synthesizing of random flow generation technique (DSRFG) with a prominent feature that the generated wind velocity fluctuations satisfy any target spectrum and target profiles of turbulence intensity and turbulence integral length scale. The new dynamic SGS model takes both advantages of one-equation SGS model and a dynamic production term without test-filtering operation, which is particular suitable to relative coarse grid situations and high Reynolds number flows. The results of comparative investigations with and without generation of inflow turbulence show that: (1) proper simulation of an inflow turbulent field is essential in accurate evaluation of dynamic wind loads on a tall building and the prescribed inflow turbulence characteristics can be adequately imposed on the inflow boundary by the DSRFG method; (2) the DSRFG can generate a large number of random vortex-like patterns in oncoming flow, leading to good agreements of both mean and dynamic forces with wind tunnel test results; (3) The dynamic mechanism of the adopted SGS model behaves adequately in the present LES and its integration with the DSRFG technique can provide satisfactory predictions of the wind effects on the super-tall building.

• 한국전산유체공학회지
• /
• 제16권2호
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.