• Journal of computational fluids engineering
• /
• v.15 no.2
• /
• pp.21-27
• /
• 2010

Immersed boundary lattice Boltzmann method (IBLBM) has been applied to simulate a turbulent flow over circular cylinder in a flow field effectively. Although IBLBM is very effective method to simulate the flow over a complex shape of obstacle in the flow field regardless of the constructed grids in the calculation domain, the results, however, become numerically unstable in high reynolds number flow. The most effective suggestion to archive the numerical stability in high Reynolds number flow is applying the multiple relaxation time (MRT) model instead of single relaxation time(SRT) model in the collision term of lattice Boltzmann equation. In the research MRT model for IBLBM was introduced and comparing the numerical results obtained by applying SRT and MRT. The hydraulic characteristic of cylinder in a flow field between two parallel plate at the range of $Re{\leqq}2000$represented and it is also compared the drag and lifting coefficients of the cylinder calculated by IBLBM with SRT and MRT model.

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

• Journal of Power System Engineering
• /
• v.4 no.3
• /
• pp.40-47
• /
• 2000

The turbulent flow with wake, reattachment and recirculation flow is very important from the viewpoint of engineering. But that is still difficult because of especially the unsteady problems which are related with the vehicle dynamics and the aerodynamics noise. This paper evaluate LES that can analyze about all fluid flow region including the laminar, transition and turbulent. So we compare the results of LES with those of PIV measurement and Reynolds averaging models. In conclusion, LES predicts flow behavior better than Reynolds averaging models.

• Journal of computational fluids engineering
• /
• v.1 no.1
• /
• pp.64-70
• /
• 1996

Three dimensional turbulent flow fields around ships are simulated by a numerical method. Reynolds Averaged Navier-Stokes equations are used where Reynolds stresses are approximated by Baldwin-Lomax and Sub-Grid Scale(SGS) turbulence models. Body-fitted coordinate system is introduced to conform three dimensional ship geometries. The governing equations are discretized by a finite volume method. Temporal derivatives are approximated by the forward differencing and the convection terms are approximated by the QUICK or Kawamura scheme. The 2nd-order centered differencing is used for other spatial derivatives. Pressure and velocity fields are simultaneously iterated by the Highly Simplified Marker-And-Cell method. To verify the numerical method and turbulence models, flow fields around ships are simulated and compared to the experiments.

• Journal of computational fluids engineering
• /
• v.20 no.3
• /
• pp.94-103
• /
• 2015

The turbulent flow characteristics in the channel flow are investigated using large eddy simulation(LES) of FDS code, built in NIST(USA), in which the near-wall flow is solved by Werner-Wengle wall function. The periodic flow condition is applied in streamwise direction to get the fully developed turbulent flow and symmetric condition is applied in lateral direction. The height of the channel is H=1m, and the length of the channel is 6H, and the lateral length is H. The total grid is $32{\times}32{\times}32$ and $y^+$ is kept above 11 to fulfill the near-wall flow requirement. The Smagorinsky model is used to solve the sub-grid scale stress. Smagorinsky constant $C_s$ is 0.2(default in FDS). Three cases of Reynolds number(10,700, 26,000, 49,000.), based on the channel height, are analyzed. The simulated results are compared with direct numerical simulation(DNS) and particle image velocimetry(PIV) experimental data. The linear low-Re eddy viscosity model of Launder & Sharma and non-linear low-Re eddy viscosity model of Abe-Jang-Leschziner are utilized to compare the results with LES of FDS. Reynolds normal stresses, Reynolds shear stresses, turbulent kinetic energys and mean velocity flows are well compared with DNS and PIV data.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.24-29
• /
• 2011

Immersed boundary lattice Boltzmann method has been applied to analyze the characteristics of the self-propelled fish motion swimming robot. The airfoil NACA0012 with caudal fin stroke model was considered to examine the characteristics. The foil in steady forward motion and a combination of steady-state harmonic deformation produces thrust through the formation of a flow downstream from the trailing edge. The harmonic motion of the foil causes unsteady shedding of vorticity from the trailing edge, while forming the vortices at the leading edge as well. The resultant thrust is developed by the pressure difference formed on the upper and lower surface of the airfoil. and the time averaged thrust coefficient increases as Re increase in the region of $Re{\leqq}700$. The suggested numerical method is suitable to develop the fish-motion model to control the swimming robot, however It would need to extend in 3D analysis to examine the higher Re and to determine the more detail mechanism of thrust production.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.265-273
• /
• 2018

Direct numerical and large eddy simulations of transitional flows around studs installed on flat plate and bulbous bow have been performed to investigate an effectiveness of turbulence stimulators on laminar-to-turbulence transition at a very low speed. The flow velocity was determined to be 0.366m/s corresponding to 4 knots of full-scale ship speed when the objective ship was Kriso container ship. The spatial evolution of skin friction coefficient disclosed that a fully development of turbulence was observed behind the second stud installed on flat plate while a rapid transition from laminar to turbulence gave rise to the fully development of turbulence behind the first stud installed on bulbous bow. A comparison of streamwise mean velocity profiles showed that the viscous sublayer and log-layer were in good agreement with previous results although the friction velocity of Smagrosinsky sub-grid scale model was about 10% larger than that of direct numerical simulation. While the turbulence intensities of bulbous bow was similar to those of flat plate in inner region, larger intensities of turbulence were observed in outer region of bulbous bow than those of flat plate.

• Journal of Korean Society for Atmospheric Environment
• /
• v.19 no.5
• /
• pp.595-610
• /
• 2003

The statistical indexes such as RMSE (Root Mean Square Error), Mean Bias error, and IOA (Index of agreement) are used to evaluate 3 Dimensional wind and temperature fields predicted by operational meteorological model RAMS (Regional Atmospheric Meteorological System) implemented in CARIS (Chemical Accident Response Information System) for the dispersion forecast of hazardous chemicals in case of the chemical accidents in Korea. The operational atmospheric model, RAMS in CARIS are designed to use GDAPS, GTS, and AWS meteorological data obtained from KMA (Korean Meteorological Administration) for the generation of 3-dimensional initial meteorological fields. The predicted meteorological variables such as wind speed, wind direction, temperature, and precipitation amount, during 19 ∼ 23, August 2002, are extracted at the nearest grid point to the meteorological monitoring sites, and validated against the observations located over the Korean peninsula. The results show that Mean bias and Root Mean Square Error are 0.9 (m/s), 1.85 (m/s) for wind speed at 10 m above the ground, respectively, and 1.45 ($^{\circ}C$), 2.82 ($^{\circ}C$) for surface temperature. Of particular interest is the distribution of forecasting error predicted by RAMS with respect to the altitude; relatively smaller error is found in the near-surface atmosphere for wind and temperature fields, while it grows larger as the altitude increases. Overall, some of the overpredictions in comparisons with the observations are detected for wind and temperature fields, whereas relatively small errors are found in the near-surface atmosphere. This discrepancies are partly attributed to the oversimplified spacing of soil, soil contents and initial temperature fields, suggesting some improvement could probably be gained if the sub-grid scale nature of moisture and temperature fields was taken into account. However, IOA values for the wind field (0.62) as well as temperature field (0.78) is greater than the 'good' value criteria (> 0.5) implied by other studies. The good value of IOA along with relatively small wind field error in the near surface atmosphere implies that, on the basis of current meteorological data for initial fields, RAMS has good potentials to be used as a operational meteorological model in predicting the urban or local scale 3-dimensional wind fields for the dispersion forecast in association with hazardous chemical releases in Korea.

• Journal of Korean Society for Atmospheric Environment
• /
• v.26 no.2
• /
• pp.137-150
• /
• 2010

A meso-urban meteorological model (Urbanized MM5; uMM5) with urban canopy parameterization (UCP) was applied to the high-resolution simulation of meteorological fields in a complex coastal urban area and the assessment of urban impacts. Multi-scale simulations with the uMM5 in the innermost domain (1-km resolution) covering the Busan metropolitan region were performed during a typical sea breeze episode (4~8 August 2006) with detailed fine-resolution inputs (urban morphology, land-use/land-cover sub-grid distribution, and high-quality digital elevation model data sets). An additional simulation using the standard MM5 was also conducted to identify the effects of urban surface properties under urban meteorological conditions. Results showed that the uMM5 reproduced well the urban thermal and dynamic environment and captured well the observed feature of sea breeze. When comparison with simulations of the standard MM5, it was found that the uMM5 better reproduced urban impacts on temperature (especially at nighttime) and urban wind flows: roughness-induced deceleration and UHI (Urban Heat Island)-induced convergence.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.1
• /
• pp.40-48
• /
• 2003

Large eddy simulation is used to investigate the compressible flow over a open cavity, 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 results show a typical flow pattern of the shear layer mode of 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. Predicted acoustic resonant frequency is in good agreement with that of the empirical formula. The mean flow streamlines are nearly horizontal along the mouth of the cavity. The pressure has its minimum value in the vortex core inside the cavity.