• 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 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.

• 대한기계학회논문집B
• /
• 제30권5호
• /
• pp.438-445
• /
• 2006

We performed the simulation of the unsteady three dimensional flow over a square cylinder in a wind tunnel in moderate Reynolds number range, $100{\sim}2500$ by using LBM. SGS model was applied for the turbulent flow. Frist of all we compared LBM(Lattice Boltzmann Method) solution of Poiseuille flow applied Farout and bounce back boundary conditions with the analytical and FOAM solutions to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured grids and prescribed uniform velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback and wind tunnel boundary conditions were applied to the cylinder walls and the boundaries of calculation domain respectively. The maximum Strouhal number of the vortex shedding is 0.2025 at Re = 750. and the number maintains the constant value of 0.18 when Re>1000. We also predicted that the critical reynolds number of the turbulent flow is in the range of $250{\sim}500$.

• 한국산학기술학회논문지
• /
• 제10권3호
• /
• pp.507-514
• /
• 2009

본 연구에서는 재순환 유동 현상을 포함하는 후향 계단에서 난류 유동장에 대한 LES의 예측성능을 검토하였다. LES의 난류모델로서 Localized Dynamic ksgs-equation 모델이 적용되었으며, 계산시간의 절감을 위하여 16개의 프로세서를 이용한 병렬계산이 수행되었다 후향 계단의 층류 유동에 대한 직접수치모사(DNS)의 수행 결과, 본 계산 결과는 기존의 실험 및 수치결과를 매우 잘 예측하였다. 또한 중간 및 높은 Re 수에 해당되는 난류 영역의 LES 결과는 평균 재순환 유동특성을 비교적 잘 예측하였다. 위 결과를 통해 본 연구에서 개발된 LES 프로그램은 향후 실용 연소기에서 연소 불안정성 및 연소 소음 등의 해석에 유용할 것으로 기대된다.

• 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.

• 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
• /
• 제24권10호
• /
• pp.1409-1416
• /
• 2000

In this study, turbulent flows in a planar combustor which has a square rib-type flame holder are numerically investigated by Large Eddy Simulation(LES). Firstly, the flow fields with or without jet injection downstream of the flame-holder are examined using uniform inlet velocity. Comparison of the present LES results with experimental one shows a good agreement. Secondly, to investigate mixing of oxidizer(air) and fuel injected behind the flame holder, the scalar-transport equation is introduced and solved. From the instantaneous flow and scalar fields, complex and intense mixing phenomena between fuel and jet are observed. It is shown that the ratio of jet to blocked air velocity is an important factor to determine the flow structure. Especially, when the ratio is large enough, the fuel jet penetrates the main vortices shed from the flame holder, resulting in significant changes in the flow and scalar fields.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2003년도 춘계학술대회
• /
• pp.1806-1813
• /
• 2003

The stirred tank reactor is one of the most commonly used devices in industry for achieving mixing and reaction. Here we report on results obtained from the large eddy simulations of flow inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius. Stereoscopic PIV measurements (Hill et $al.^{(1)}$) along with the theoretical model of the impeller-induced flow (Yoon et $al.^{(2)}$) were used in defining the impeller-induced flow as superposition of circumferential, jet and tip vortex pair components. Large eddy simulation of flow in a stirred tank was carried out for the three different Reynolds numbers of 4000, 16000 and 64000. The effect of different Reynolds numbers is well observed in both instantaneous and time averaged flow fields. The instantaneous and mean vortex structures are identified by plotting an isosurfaces of swirling strength for all Reynolds numbers. The Reynolds number dependency of the nondimeansional eddy viscosity, resolve scale and subgrid scale dissipations is clearly shown in this study.

• 대한기계학회논문집B
• /
• 제27권9호
• /
• pp.1290-1298
• /
• 2003

The stirred tank reactor is one of the most commonly used devices in industry for achieving mixing and reaction. Here we report on results obtained from the large eddy simulations of flow inside the tank performed using a spectral multi-domain technique. The computations were driven by specifying the impeller-induced flow at the blade tip radius. Stereoscopic PlY measurements (Hill et al. $^{(1)}$) along with the theoretical model of the impeller-induced flow (Yoon et al. $^{(2)}$) were used in defining the impeller-induced flow as superposition of circumferential, jet and tip vortex pair components. Large eddy simulation of flow in a stirred tank was carried out for the three different Reynolds numbers of 4000, 16000 and 64000. The effect of different Reynolds numbers is well observed in both instantaneous and time averaged flow fields. The instantaneous and mean vortex structures are identified by plotting an isosurfaces of swirling strength for all Reynolds numbers. The Reynolds number dependency of the non-dimensional eddy viscosity, resolved scale and subgrid scale dissipations is clearly shown in this study.

• 동력기계공학회지
• /
• 제16권4호
• /
• pp.51-59
• /
• 2012

The computational flow analysis using LES technique was carried out to investigate the flow characteristics of a sirocco fan under the maximum flowrate condition. The commercial SC/Tetra software was used for this unsteady and three-dimensional numerical analysis. In consequence, because a flow is unstable within the range of about 50% of a housing depth from a bellmouth around the cutoff region, the passing flow through the blade cascade occurred on the X-Y plane is a slow or a reverse with approaching to the housing inlet. Also, the secondary flow shows on the radial plane of a housing, and its vortex center exists within about 33% of a housing depth from a bellmouth except the cutoff region. Moreover, the flow occurring on the exit plane of a sirocco fan shows a complex secondary flow.