• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• 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.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.25 no.3
• /
• pp.182-188
• /
• 2016

This study investigates cavity flows through a guide grill above a resonator. Vortex distributions and intake flows are simulated for various shapes of the guide grill. The flows are assumed to be compressible, unsteady, and turbulent. Numerical simulations are conducted using a large eddy simulation (LES) model. To analyze the effect of the guide grill shape, three cavity lengths (0.2H, 0.6H, and 1.0H) and cavity angles ($30^{\circ}$, $45^{\circ}$ and $60^{\circ}$) are considered based on resonator height (H). The results show that the vortex generated in the resonator by cavity flow increases with cavity length. Thus, the intake flow is minimum at the smallest cavity length and angle. However, when cavity length is equal to resonator height, the intake flow decreases. The maximum intake flow occurs at a cavity angle $45^{\circ}$ at higher cavity lengths owing to the interaction between the vortex in the resonator and intake flow.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.414-419
• /
• 2011

Missile am fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 60 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure distribution, etc. are analyzed according to the angle of attack. The results at a low angle of attack are compared with other results before a stall condition. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. Unsteady velocity field, periodic vortex shedding, the unsteady pressure distribution on the airfoil surface, and the acoustic fields are analyzed. The effects of these unsteady characteristics in the aerodynamic coefficients are analyzed.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.1-6
• /
• 2013

Missile and fighter aircraft have been challenged by low restoring nose-down pitching moment at high angle of attach. The consequence of weak nose-down pitching moment can be resulting in a deep stall condition. Especially, the pressure oscillation has a huge effect on noise generation, structure damage, aerodynamic performance and safety, because the flow has strong unsteadiness at high angle of attack. In this paper, the unsteady aerodynamics coefficients were analyzed at high angle of attack up to 50 degrees around two dimensional NACA0012 airfoil. The two dimensional unsteady compressible Navier-Stokes equation with a LES turbulent model was calculated by OHOC (Optimized High-Order Compact) scheme. The flow conditions are Mach number of 0.3 and Reynolds number of $10^5$. The lift, drag, pressure, entropy distribution, etc. are analyzed according to the angle of attack. The results of average lift coefficients are compared with other results according to the angle of attack. From a certain high angle of attack, the strong vortex formed by the leading edge are flowing downstream as like Karman vortex around a circular cylinder. The primary and secondary oscillating frequencies are analyzed by the effects of these unsteady aerodynamic characteristics.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• pp.1984-1989
• /
• 2008

In this study, fire simulations were performed to analyze the characteristics of the fire driven flow and the effects of the platform screen door on the smoke flow in the station, when the fire occurred in the center of the platform. Soongsil Univ. station (line number 7, 47m in depth underground) was chosen which was the one of the deepest underground subway stations in the Seoul metro, SMRT. The parallel computational method was employed to compute the heat and mass transfer eqn's with 6 CPUs of the linux clustering machine. The fire driven flow was simulated with using FDS code in which LES method was applied. The Heat release rate was 10MW and The Ultrafast model was applied for the growing model of the fire source. The 10,000,000 structured grids were used.

• Ocean Systems Engineering
• /
• v.9 no.1
• /
• pp.21-42
• /
• 2019

In this paper, the Finite-Analytic Navier-Stokes (FANS) code is coupled with an in-house finite-element code to study the dynamic interaction between a floating buoy and its mooring system. Hydrodynamic loads on the buoy are predicted with the FANS module, in which Large Eddy Simulation (LES) is used as the turbulence model. The mooring lines are modeled based on a slender body theory. Their dynamic responses are simulated with a nonlinear finite element module, MOORING3D. The two modules are coupled by transferring the forces and displacements of the buoy and its mooring system at their connections through an interface module. A free-decay model test was used to calibrate the coupled method. In addition, to investigate the capability of the present coupled method, numerical simulations of two degree-of-freedom vortex-induced motion of a CALM buoy in uniform currents were performed. With the study it can be verified that accurate predictions of the motion responses and tension responses of the CALM buoy system can be made with the coupling CFD-FEM method.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.25 no.2
• /
• pp.24-33
• /
• 2021

Numerical study is conducted to analyze combustion instability in the partially premixed combustor. The simulations are performed according to fuel conditions, and Large Eddy Simulation(LES) model and PaSR combustion model are implemented in the solver. Comparison with the experimental result is conducted to confirm the validity of simulation, and quantitative and qualitative agreement is confirmed. The flame characteristics in the combustor are subsequently investigated, and the association with the occurrence of combustion instability is clarified. According to the simulation results, the flame length varies greatly depending on the fuel conditions. When the flame length becomes sufficiently long, flame-vortex interactions occurred around the wall sections, which works as the main cause of combustion instability.

• Journal of Coastal Disaster Prevention
• /
• v.1 no.1
• /
• pp.1-9
• /
• 2014

In order to examine the characteristics of tsunami run-up heights on impermeable/permeable slope, a numerical wave tank by upgrading LES-WASS-3D was used in this study. Then, the model were compared with existing hydraulic model test for its verification. The numerical results well reproduced experimental results of solitary wave deformation, propagation and run-up height under various conditions. Also, the numerical simulation with a slope boundary condition has been carried out to understand solitary wave run-up on impermeable/permeable slope. It is shown that the run-up heights on permeable slope is 52.64-63.2% smaller than those on the impermeable slope because of wave energy dissipation inside the porous media. In addition, it is revealed that the numerical results with slope boundary condition agreed well with experimental results in comparison with the results by using stair type boundary condition.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.10 no.1
• /
• pp.111-117
• /
• 2010

A three-dimensional numerical model is employed to investigate wave deformation due to a submerged structure. The three-dimensional numerical model solves the spatially averaged Navier-Stokes equations for two-phase flows. The LES(large-eddy-simulation) approach is adopted to model the turbulence effect by using the Smagorinsky SGS(sub-grid scale) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate VOF(volume-of-fluid) method is used to track the distorted and broken free surface. A simple linear wave is generated on a constant depth and compared with analytical solutions. The model is then applied to study wave deformation due to a submerged structure and the predicted results are compared with available laboratory measurements.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.44 no.4
• /
• pp.343-351
• /
• 2016

Cryogenic spray characteristics of a nitrogen swirl injector operating in supercritical environment have been numerically investigated. By comparing the equation of states(EOS) used for supercritical condition, SRK EOS was applied to predict the nitrogen thermodynamic property under supercritical environment. A Chung's method was implemented for the calculation of viscosity and conductivity and Takahashi's correlation based on Fuller's Theorem was implemented for the calculation of diffusion coefficient. By injecting the nitrogen with 5 bar differential pressure into 50 bar chamber filled with nitrogen, numerical simulation has been conducted. The dynamic Smagorinsky sub-grid scale (SGS) model has been compared with the algebraic Smagorinsky SGS model using FFT frequency analysis. The instability at the liquid film and gas core inside injector and the propagation of pressure oscillation into the injector has been investigated. The spreading angle of swirl injector obtained by numerical calculation has been validated with experimental result.