• The Transactions of the Korea Information Processing Society
• /
• v.7 no.1
• /
• pp.57-64
• /
• 2000

We describe a parallel implementation for the large-eddy simulation(LES) of stratified and rotating turbulence based on MPI. The parallelization strategy is specified by eliminating the tridiagonal solver with explicit method and by domain decompositions for solving the poisson equation. In this simulation we have run on CRAY-T3E under the message passing platform MPI with a various domain decomposition and the scalability of this parallel code of LES are also presented. The result shows that we can gain up to 16 times faster speed up on 64 processors with xyz-directional domain decomposition and scalable up to $128{\times}128{\times}$ which processing time is almost similar to that of $40{\times}40{\times}40$ on a single processor machine with a sequential code.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.17 no.11
• /
• pp.1060-1071
• /
• 2005

Passive gas dispersions over a 1/1000 scale terrain model at Eiffel type wind tunnel were reproduced by numerical simulation. Large eddy simulation was used to treat the sub-grid scale turbulences. The terrain features were represented by millions of point forces densely distributed over the solid surface using the virtual boundary method. The model simulations agreed very well with the experiments in a consistent fashion for all wind directions. The measured profiles of the wind speeds as well as the tracer gas concentrations were nicely simulated by the CFD model at most locations scattered over the model terrain. With scale factor adjusted and the thermal stratification effects incorporated, the CFD model was expected to provide reliable information on pollutant dispersions over the real complex terrains.

• Ecology and Resilient Infrastructure
• /
• v.5 no.4
• /
• pp.199-209
• /
• 2018

In river design, consideration of bed shear stresses is necessary to secure stability of levee and floodplain. In this study distributions of bed shear stresses in compound open channels are analyzed through numerical simulation for various width and depth. LES solver in OpenFOAM is applied to 12 cases of compound channel shapes considering secondary flow which effects distributions of bed shear stresses. By the results time averaged velocity distributions, secondary currents, and distributions of bed shear stresses are analyzed. Overall distributions of bed shears in floodplain show that higher shear stresses are seen in left of floodplain and the shears decrease toward right of floodplain. However, high local variations in shear stresses are shown due to the secondary flow effects. In shallow floodplain, bed shear stresses show low value below 0.8 times of averaged bed shear. In deep floodplain, bed shear stresses show high value over 1.2 - 1.4 times of averaged bed shear.

• Journal of the Korean Society of Safety
• /
• v.28 no.1
• /
• pp.57-62
• /
• 2013

Aerodynamic characteristics of cross section shape is an important parameter for the wind response and structural stability of long span bridges. Numerical simulation methods have been introduced to estimate the aerodynamic characteristics for more detailed flow analysis and cost saving in place of existing wind tunnel experiment. In this study, the computational fluid dynamics(CFD) simulation and large eddy simulation( LES) technique were used to estimate lift, drag and moment coefficients of four cross sections. The Strouhal numbers were also determined by the fast Fourier transform of time series of the lift coefficient. The values from simulations and references were in a good agreement with average difference of 16.7% in coefficients and 8.5% in the Strouhal numbers. The success of the simulations is expected to attribute to the practical use of numerical estimation in construction engineering and wind load analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.161-162
• /
• 2006

• 한국연소학회:학술대회논문집
• /
• 2015.12a
• /
• pp.327-328
• /
• 2015

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.628-637
• /
• 2004

The present study investigates the performance of dynamic mixed model (DMM; Zang et ai.) in a channel with strong wall injection through a Large eddy simulation (LES) technique. The DMM results are compared with those of DNS and the results obtained with popular dynamic Smagorinsky model (DSM). Better agreement is achieved when using the DMM with box filter than DSM and coarse DNS in predicting the first and second order statistics as well as large-scale structures of velocity and temperature fields. Such favorable features of DMM are attributed to the fact that it explicitly calculates the modified Leonard stress term and only models the remaining cross and the SGS Reynolds stress terms and, thus, it reduces the excessive burden put on the model coefficient of DSM. Also it is demonstrated that the DMM can be successfully extended to the prediction of temperature (passive scalar) field where strong streamwise inhomogeneity exists.

• Wind and Structures
• /
• v.23 no.4
• /
• pp.313-350
• /
• 2016

Taipei 101 Tower, which has 101 stories with height of 508 m, is located in Taipei where typhoons and earthquakes commonly occur. It is currently the second tallest building in the world. Therefore, the dynamic performance of the super-tall building under strong wind actions requires particular attentions. In this study, Large Eddy Simulation (LES) integrated with a new inflow turbulence generator and a new sub-grid scale (SGS) model was conducted to simulate the wind loads on the super-tall building. Three-dimensional finite element model of Taipei 101 Tower was established and used to evaluate the wind-induced responses of the high-rise structure based on the simulated wind forces. The numerical results were found to be consistent with those measured from a vibration monitoring system installed in the building. Furthermore, the equivalent static wind loads on the building, which were computed by the time-domain and frequency-domain analysis, respectively, were in satisfactory agreement with available wind tunnel testing results. It has been demonstrated through the validation studies that the numerical framework presented in this paper, including the recommended SGS model, the inflow turbulence generation technique and associated numerical treatments, is a useful tool for evaluation of the wind loads and wind-induced responses of tall buildings.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.34 no.7
• /
• pp.79-88
• /
• 2006

In the present paper, the swirl flow structure and flame characteristics of turbulent premixed combustion in a model gas turbine combustor are investigated using large eddy simulation(LES). A G-equation flamelet model is employed to simulate the unsteady flame behavior. When inlet swirl number is increased, the distinct flow structures, such as the shapes of corner recirculation and center toroidal recirculation zone, are observed and the flame length is shorted gradually. Also, the phenomena of flashback are identified at strong swirl intensity. In order to get the accurate description of unsteady flame behavior, the predictive ability of the acoustic wave in a combustor is primarily evaluated. It is found that the vortex generated near the edge of step plays an important role in the flame fluctuation. Finally it is examined systematically that the flame and heat release fluctuation are coupled strongly to the vortex shedding generated by swirl flow and acoustic wave propagation from the analysis of flame-vortex interaction.

• Journal of the Society of Naval Architects of Korea
• /
• v.49 no.5
• /
• pp.440-446
• /
• 2012

A twisted cylinder has been newly designed by rotating the elliptic cross section along the spanwise direction in order to reduce the drag and vorticies in wake region. The flow around the twisted cylinder at a subcritical Reynolds number (Re) of 3000 is investigated to analyze the effect of twisted spiral pattern on the drag reduction and vortex suppression using large eddy simulation (LES). The instantaneous wake structures of the twisted cylinder are compared with those of a circular and a wavy cylinder at the same Re. The shear layer of the twisted cylinder covering the recirculation region is more elongated than that of the circular and the wavy cylinder. Successively, vortex shedding of the twisted cylinder is considerably suppressed, compared with those of the circular and the wavy cylinder. Consequently, the mean drag coefficient and the fluctuating lift of the twisted cylinder are less than those of the circular and the wavy cylinder.