• Wind and Structures
• /
• 제1권2호
• /
• pp.127-144
• /
• 1998

In this paper, the wake structures behind a square cylinder at the Reynolds number of 22,000 are simulated using the large eddy simulation, and the main features of the wake structure associated with unsteady vortex-shedding are investigated. The Smagorinsky model is used for parametrization of the subgrid scales. The finite element method with isoparametric linear elements is employed in the computations. Unsteady computations are performed using the explicit method with streamline upwind scheme for the advection term. The time integration incorporates a subcycling strategy. No-slip condition is enforced on the wall surface. A comparative study between two-and three-dimensional computations puts a stress on the three-dimensional effects in turbulent flow simulations. Simulated three-dimensional wake structures are compared with numerical and experimental results reported by other researchers. The results include time-averaged, phase-averaged flow fields and numerically visualized vortex-shedding pattern using streaklines. The results show that dynamics of the vortex-shedding phenomenon are numerically well reproduced using the present method of finite element implementation of large eddy simulation.

• Journal of Advanced Marine Engineering and Technology
• /
• 제28권1호
• /
• pp.98-108
• /
• 2004

This study concerns a large eddy simulation (LES) of turbulent flow around a bluff body inside a sudden expansion cylindrical chamber, a configuration which resembles a premixed gas turbine combustor The simulation code is constructed by using the general coordinate system based on the physical contravariant velocity components. The Smagorinsky model is employed and the calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of the inlet pipe. The combined grid technique and cylindrical grid are tested in the numerical simulation with complex geometry. The predicted turbulent statistics are evaluated by comparing with LDV measurement data. The numerical flow visualizations depict the behavior of turbulent mixing process behind the flame holder.

• 한국화재소방학회논문지
• /
• 제32권3호
• /
• pp.42-50
• /
• 2018

Large Eddy Simulation (LES) 기반의 Fire Dynamics Simulator (FDS) Version(5와 6)에 따른 아격자 난류 및 연소모델의 변화가 헵탄 풀화재의 평균 화염높이에 미치는 영향이 검토되었다. 화재시뮬레이션을 위한 열발생률은 동일한 조건에서 수행된 실험을 통해 제공되었으며, FDS Version에 따른 평균 화염높이의 예측성능은 기존 상관식과의 비교를 통해 평가되었다. FDS 5와 6에 각각 적용된 Smagorinsky 및 Deardorff 난류 모델은 평균 유동장, 화염형상 및 화염높이에 큰 영향을 주지 않는다. 그러나 평균 화염높이를 포함한 풀화재 특성 차이는 주로 FDS 5와 6에 각각 적용된 혼합분율 및 Eddy Dissipation Concept (EDC) 연소모델 차이로 인하여 발생됨을 확인하였다. 마지막으로 FDS 6에 비해 FDS 5는 상당히 긴 화염높이의 예측결과를 제공하며, 기존 상관식과 보다 일치하는 평균 화염높이 예측이 가능함을 확인하였다.

• Wind and Structures
• /
• 제14권5호
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.

• 한국방재학회 논문집
• /
• 제10권1호
• /
• pp.111-117
• /
• 2010

수중구조물에 의한 파랑의 변형을 예측하기 위해 3차원 수치모형을 도입하여 수치모형 실험을 수행하였다. 본 수치모형은 Navier-Stokes 방정식을 유한차분법을 이용하여 계산하는 동수압 모형으로서, 난류의 해석을 위해서 상대적으로 큰 에디(eddy)만을 고려하는 SANS(Spatially Averaged Navier-Stokes) 방정식의 해를 구하는 LES(large-eddy-simulation) 기반의 수치모형이다. 엇갈림 격자체계에서 유한차분법을 사용하여 지배방정식을 해석하는 모형으로서 수치기법으로 Two-step projection 기법을 사용하여 SANS 방정식을 계산하였으며, Bi-CGSTAB 기법을 이용하여 Poisson 방정식의 해를 구하고 압력장을 계산하였다. 또한, 자유수면의 추적을 위하여 2차 정확도의 VOF(volume-of-fluid) 기법을 사용하였다. 먼저 선형파를 일정 수심상에서 조파시켜 해석해와 비교한 후 수중구조물이 설치된 지형에 적용하여 파랑의 변형을 수치모의하여 수리모형 실험결과와 비교 및 분석하였다.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2003년도 추계학술대회 논문집
• /
• pp.204-212
• /
• 2003

As the first step to investigate the nonlinear interactions between turbulence and marine structures inside a viscous NWT, a LES technique was applied to the turbulent channel flow for $Re_{T}=150$, in this paper. The employed models were 4 types, such as the Smagorinsky model, the Dynamic SGS model, the Structure Function model and the Generalized Normal Stress model. The simulated data in time-series for the LESs were averaged in both time and space and performed statistical analysis. And results of the LESs were compared with those of a DNS developed in the present study and two spectral methods by Yoon et al.(2003) & Kim et al.(1987). It seems to be quite difficult to evaluate their performances to the present problem, but is seen that the accuracy of LESs are still related to the number of grids(or fine grid size).

• Wind and Structures
• /
• 제22권2호
• /
• pp.161-184
• /
• 2016

Long-span suspension bridges have evolved through the years and with them, the bridge girder decks improved as well, changing their shapes from standard box-deck girders to twin box and multi-box decks sections. The aerodynamic characteristics of the new generation of twin and multiple-decks are investigated nowadays, to provide the best design wind speeds and the optimum dimensions such bridges could achieve. The multi-box Megane bridge deck is one of the new generation bridge decks, consisting of two side decks for traffic lanes and two middle decks for railways, linked between them with connecting beams. Three-dimensional CFD simulations were performed by employing the Large Eddy Simulation (LES) algorithm with a standard Smagorinsky subgrid-scale model, for $Re=9.3{\times}10^7$ and angles of attack ${\alpha}=-4^{\circ}$, $-2^{\circ}$, $0^{\circ}$, $2^{\circ}$ and $4^{\circ}$. Also, a wind tunnel experiment was performed for a scaled model, 1:80 of the Megane bridge deck section, for $Re=5.1{\times}10^5$ and the aerodynamic static coefficients were found to be in good agreement with the results obtained from the CFD-LES model. However the aerodynamic coefficients determined individually, from the CFD-LES model, for each of the traffic and railway decks of the Megane bridge, varied significantly, especially for the downstream traffic deck. Also the pressure distribution and the effect of the spacing between the connecting beams, on the wind speed profiles showed a slight increase in turbulence above the downstream traffic and railway decks.

• 대한기계학회논문집B
• /
• 제25권7호
• /
• pp.905-914
• /
• 2001

A large eddy simulation(LES) is performed for turbulent flow around a bluff body inside a sudden expansion cylinder chamber, a configuration which resembles a premixed gas turbine combustor. To promote turbulent mixing and to accommodate flame stability, a flame holder is installed inside the combustion chamber. The Smagorinsky model is employed and the calculated Reynolds number is 5,000 based on the bulk velocity and the diameter of the inlet pipe. The simulation code is constructed by using a general coordinate system based on the physical contravariant velocity components. The predicted turbulent statistics are evaluated by comparing them with the laser-doppler velocimetry (LDV) measurement data. The agreement of LES with the experimental data is shown to be satisfactory. Emphasis is placed on the time-dependent evolutions of turbulent vortical structure behind the flame holder. The numerical flow visualizations depict the behavior of large-scale vortices. The turbulent mixing process behind the flame holder is analyzed by visualizing the sectional views of vortical structure.

• 대한기계학회논문집B
• /
• 제32권9호
• /
• pp.720-727
• /
• 2008

Time delay effects on near-wall turbulent structures are investigated by performing a large-eddy simulation of a transient turbulent flow in a pipe. To elucidate the time delay effects on the near-wall turbulence, we selected the dimensionless acceleration parameter which was used in the previous study. Various turbulent statistics revealed the distinctive features of the delay. It was shown that the dynamic Smagorinsky model is valid to capture the alterations of the turbulence physics well. A dimensionless time for the responses of the flow quantities was introduced to give the detailed information on the delay of the nearwall turbulence. The conditionally-averaged flow fields associated with Reynolds shear stress producing events show that sweep and ejections are closely related to the delays of the turbulence production and the turbulence propagation toward the pipe center. The present study suggested that the enhanced anisotropy of the turbulence in the initial and transient stages would be a challenging problem to standard turbulence models.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 춘계학술대회논문집E
• /
• pp.386-391
• /
• 2001

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.