• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.45 no.3
• /
• pp.233-240
• /
• 2017

The flowfield around transonic wing-body configuration was simulated using in-house CFD code and compared with the experimental data to understand the influence of several features of CFD(Computational Fluid Dynamics) ; grid dependency, turbulence models, spatial discretization, and viscosity. The wing-body configuration consists of a simple planform RAE Wing 'A' with an RAE 101 airfoil section and an axisymmetric body. The in-house CFD code is a compressible Euler/Navier-Stokes solver based on unstructured grid. For the turbulence model, the $k-{\omega}$ model, the Spalart-Allmaras model, and the $k-{\omega}$ SST model were applied. For the spatial discretization method, the central differencing scheme with Jameson's artificial viscosity and Roe's upwind differencing scheme were applied. The results calculated were generally in good agreement with experimental data. However, it was shown that the pressure distribution and shock-wave position were slightly affected by the turbulence models and the spatial discretization methods. It was known that the turbulent viscous effect should be considered in order to predict the accurate shock wave position.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.32 no.1
• /
• pp.55-68
• /
• 2020

In order to investigate the hydraulic characteristics of a boundary layer streaming over the beach cusps appeared in swells prevailing mild seas, we numerically simulated the shoaling process of Edge waves over the beach cusp. Synchronous Edge waves known to sustain the beach cusps could successfully be duplicated by generating two obliquely colliding Edge waves in front of beach cusps. The amplitude A_B and length L_B of Beach Cusp were elected to be 1.25 m and 18 m, respectively based on the measured data along the Mang-Bang beach. Numerical results show that boundary layer streaming was formed at every phase of shoaling process without exception, and the maximum boundary layer streaming was observed to occur at the crest of sand bar. In RUN 1 where the shortest waves were deployed, the maximum boundary layer streaming was observed to be around 0.32 m/s, which far exceeds the amplitude of free stream by two times. It is also noted that the maximum boundary layer streaming mentioned above greatly differs from the analytical solution by Longuet-Higgins (1957) based on wave Reynolds stress. In doing so, we also identify the recovery procedure of natural beaches in swells prevailing mild seas, which can be summarized such as: as the infra-gravity waves formed in swells by the resonance wave-wave interaction arrives near the breaking line, the sediments ascending near the free surface by the Phase II waves orbital motion were carried toward the pinnacle of foreshore by the shoreward flow commenced at the steep front of breaking waves, and were deposited near the pinnacle of foreshore due to the infiltration.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.117-122
• /
• 1998

A numerical study is carried out to investigate the transient process of combustion phenomena associated with hypersonic propulsion devices. Reynolds averaged Navier-Stokes equations for reactive flows are used as governing equations with a detailed chemistry mechanism of hydrogen-air mixture and two-equation SST turbulence modeling. The governing equations are discretized by a high order accurate upwind scheme and solved in a fully coupled manner with a fully implicit time accurate method. At first, oscillating shock-induced combustion is analyzed and the comparison with experimental result gives the validity of present computational modeling. Secondly, the model ram accelerator experiment was simulated and the results show the detailed transient combustion mechanisms. Thirdly, the evolution of oblique detonation wave is simulated and the result shows transient and final steady state behavior at off-stability condition. Finally, shock wave/boundary layer interaction in combustible mixture is studied and the criterion of boundary layer flame and oblique detonation wave is identified.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.3
• /
• pp.375-383
• /
• 2002

The flowfield of transverse jet in a supersonic air stream subjected to shock wave turbulent boundary layer interactions is simulated numerically by Generalized Taylor Galerkin(GTG) finite element methods. Effects of turbulence are taken into account with a two-equation (k-$\varepsilon$) model with a compressibility correction. Injection pressures and slot widths are varied in the present study. Pressure, separation extents, and penetration heights are compared with experimental data. Favorable comparisons with experimental measurements are demonstrated.

• 한국신재생에너지학회:학술대회논문집
• /
• 2011.11a
• /
• pp.38.1-38.1
• /
• 2011

This study investigates the effects of Pierson-Moskowitz, Jonswap spectrum that are typical irregular wave spectrums for wind turbine system with jacket support structure. Also various offshore environmental parameters based on korean local condition were used in our study. The loads acting on the system was considered by referring to the Design Load Case from IEC guide line. And improved von Karman model was used as a turbulence model. As a result, various significant wave height and peak spectral period cause noticeable difference of extreme and fatigue loads prediction.

• Journal of Ship and Ocean Technology
• /
• v.5 no.1
• /
• pp.38-54
• /
• 2001

A finite difference method for calculating turbulent flow and surface wave fields around a ship model is evaluated through the comparison with the experimental data of a Series 60 $C_B$=0.6 ship model. The method solves the Reynolds-averaged Navior-Stokes Equations using the non-staggered grid system, the four-stage Runge-Kutta scheme for the temporal integration of governing equations and the Bladwin-Lomax model for the turbulence closure. The free surface waves are captured by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and free-surface conforming grids are generated at each time step so that one of the grid surfaces coincides always with the free surface. The computational results show an overall close agreement with the experimental data and verify that the present method can simulate well the turbulent boundary layers and wakes as well as the free-surface waves.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.6
• /
• pp.1510-1517
• /
• 1995

A theoretical method is described to simulate the propagation of turbulent premixed flames in a closed vessel. The objective is to develop and test an efficient technique to predict the propagation speed of flame as well as the geometric structure of the flame surfaces. Flame is advected by the statistically generated turbulent flow field and propagates as a wave by solving twodimensional Hamilton-Jacobi equation. In the simulation of the unburned gas flow field, following turbulence properties were satisfied: mean velocity field, turbulence intensities, spatial and temporal correlations of velocity fluctuations. It is assumed that these properties are not affected by the expansion of the burned gas region. Predictions were compared with existing experimental data for flames propagating in a closed vessel charged with hydrogen/air mixture with various turbulence intensities and Reynolds numbers. Comparisons were made in flame radius growth rate, rms flame radius fluctuations, and average perimeter and fractal dimensions of the flame boundaries. Two dimensional time dependent simulation resulted in correct trends of the measured flame data. The reasonable behavior and high efficiency proves the usefulness of this method in difficult problems of flame propagation such as in internal combustion engines.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.57-66
• /
• 1999

Numerical Analysis has been done for the supersonic off-design jet flow due to the pressure difference between the jet and the ambient fluid. The difference of pressure generates an oblique shock or an expansion wave at the nozzle exit. The waves reflect repeatedly on the center axis and the sonic surface in the shear layer. The pressure difference is resolved across these reflected waves. In this paper, the axi-symmetric Navier-Stokes equation has been used with the κ-ε turbulence model. The second order TVD scheme with flux limiters, based on the flux vector split with the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. Numerical calculations have been done to analyze the off-design jet flow due to the pressure difference. The variation of pressure along the flow axis is compared with an experimental result and other numerical result. The characteristics of the interaction between the shock cell and the turbulence mixing layer have been analyzed.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.3
• /
• pp.24-33
• /
• 2008

In this study, flow and creation and dissipation of vorticity around rectangular floating breakwater is investigated both experimentally and numerically. The PIV system(Particle image velocimetry) is employed to obtain the velocity field in the vorticity of rectangular structure. The numerical model, combined with ${\kappa}-{\varepsilon}$ turbulence model and the VOF method based on RANS equation, is used to analyze the turbulence structure. In the results of this study, the vorticity is found around conner of rectangular structure at all time domain, and creation and dissipation of vorticity are closely related to wave period. Separation points of phase of vortex due to flow separation for longer period waves are faster then for shorter period waves.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.6 no.2
• /
• pp.75-84
• /
• 2002

A numerical analysis of unsteady motion in solid rocket motors with a nozzle has been conducted. The numerical formulation including modified $\kappa$-$\varepsilon$ turbulence model treats the complete conservation equation for the gas phase and the one-dimensional equations in the radial direction for the condensed phase. A fully coupled implicit scheme based on a dual time-stepping integration algorithm has been adopted to solve the governing equations. After obtaining a steady state solution, pulse and periodic oscillations of pressure are imposed at the head-end to simulate acoustic oscillations of a travelling-wave motion in the combustion chamber. Various steady and unsteady state features in the combustion chamber of a rocket motor has been analyzed as results of numerical calculations.