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

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• 한국전산유체공학회지
• /
• 제15권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.

• 한국전산유체공학회지
• /
• 제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.

• 대한조선학회논문집
• /
• 제59권3호
• /
• pp.141-148
• /
• 2022

A hybrid turbulence model has developed by combining a sub-grid scale model using dynamic k equation in LES with k-𝜔 SST model of RANS equation. To ascertain potential applicability of the hybrid turbulence model, fully developed turbulent channel flows at Re_𝜏=180 have been simulated of which computational domain has a top wall with coarse cells and a bottom wall with fine cells. The streamwise mean velocity and turbulent intensity profiles showed a good agreement with DNS data when using the hybrid model rather than using a single model in k-𝜔 SST or dynamic k equation models. Computational simulations of turbulent flows around KVLCC2 with a pre-swirl duct have been mainly performed using the hybrid turbulence model. Compared to the results obtained from RANS simulation with k-𝜔 SST model as well as LES with dynamic k equation SGS model, turbulent wakes of the duct in the present simulation using the hybrid turbulence model were very similar to that of LES. Also, the resistances acting on hull, rudder and duct in hybrid turbulence model were similar to those in RANS simulation whereas the viscous forces acting on the hull in LES had a significant error due to coarse cells inappropriate to the sub-grid scale model.

• 한국전산유체공학회지
• /
• 제16권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.

• 대한기계학회논문집B
• /
• 제28권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 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).