• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.53-59
• /
• 2007

An axisymmetric supersonic screeching jet is numerically simulated to examine the length scales of screech frequency as well as screech tone generation mechanism. The axisymmetric Reynolds-averaged Navier-Stokes equations in conjuction with a modified Spalart-Allmaras turbulence model are employed. It is demonstrated that the axisymmetric jet screech tones can be simulated correctly and the numerical results are in good agreement with the experimental data. Instability waves, shock-cell structures and the phenomena of shock motion are investigated in detail to identify the screech tone generation mechanism. Shock spacings and standing wave length are analyzed to determine the dominent length scale crucial to the screech frequency formulation.

• Advances in aircraft and spacecraft science
• /
• v.5 no.4
• /
• pp.445-458
• /
• 2018

Turbulent airflow in channels of rectangular cross section with symmetric centerbodies is studied numerically. Shock wave configurations formed in the channel and in front of the entrance are examined. Solutions of the unsteady Reynolds-averaged Navier-Stokes equations are obtained with finite-volume solvers of second-order accuracy. The solutions demonstrate an expulsion/swallowing of the shocks with variations of the free-stream Mach number or angle of attack. Effects of the centerbody length and thickness on the shock wave stability and flow bifurcation are examined. Bands of the Mach number and angle of attack, in which there exist non-unique flow fields, are identified.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• Journal of Ocean Engineering and Technology
• /
• v.11 no.2
• /
• pp.70-76
• /
• 1997

수치격자의 크기제한에 의한 자유표면 유동해석 문제를 효율적으로 다루기 위하여 자유표면의 모든 격자를 x 방향으로 4,8,12개로 등분할 하고, y 방향으로는 4개로 잘라서 계산하였다. 이중격자 또는 삼중격자로 Navier-Stokes 방정식의 각항에 크기가 다른 격자를 사용해 효율을 향상시키는 계산방법의 연장으로, 본 논문에서는 자유표면 방정식에 보다 세분화된 격자를 적용해, Marker Particle 이동 및 자유수면 형성에 효율향상을 줄수 있는 수치방법을 도입하였다. 계산결과에 의하면 초기사용 격자가 coarse한 경우가 본방법의 효과가 커짐을 알 수 있고 대상물로는 층류유동에서 Wigley모형과 낭류유동의 S103 모형이다.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2003.10a
• /
• pp.226-231
• /
• 2003

The viscous interaction of stern wave is studied by simulating the free-surface flows, including sub-breaking phenomena around a high speed catamaran hull advancing on calm water. The Navier-Stokes equation is solved by a finite difference method where the body-fitted coordinate system, the wall function and the triple-grid system are invoked. The numerical appearance of the sub-breaking waves is qualitatively supported by the experimental observation They are also applied to study precisely on the stern flow of S-103 as to which extensive experimental data are available. For the catamaran, computations are carried out for the mono ana twin hulls.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.7
• /
• pp.2347-2355
• /
• 1996

Purpose of the present study is to develop a computational design program for shape optimization, combining the numerical optimization technique with the flow analysis code. The present methodology is then validated in three cases of aerodynamic shape optimization. In the numerical optimization, a feasible direction optimization algorithm and shape functions are considered. In the flow analysis, the Navier-Stokes equations are discretized by a cell-centered finite volume method, and Roe's flux difference splitting TVD scheme and ADI method are used. The developed design code is applied to a transonic channel flow over a bump, and an external flow over a NACA0012 airfoil to minimize the wave drag induced by shock waves. Also a separated subsonic flow over a NACA0024 airfoil is considered to determine a maximum allowable thickness of the airfoil without separation.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 1996.11a
• /
• pp.107-116
• /
• 1996

Numerical simulations have been conducted to characterize unsteady flow structures in an axisymmetric supersonic inlet diffuser with sinusoidal pressure oscillations at the diffuser exit. The formulation is based on the unsteady Navier-Stokes equations and turbulence closure is achieved using a two-layer model with a too-Reynolds-number scheme for the near-wall treatment. The governing equations are formulated in an integral form, and are discretized by the four-stage Runge-Kutta scheme for temporal terms and the Harten-Yee upwind TVD scheme for convective terms. Results indicated that the inlet shock characteristics are significantly modified by acoustic oscillations originating from the combustor. The characteristics of shock/boundarv-layer interactions (such as the size of separation bubble, terminal shock shape, and vorticity intensity) are also greatly iufluenced by the shock oscillation due to acoustic waves.

• Ocean Systems Engineering
• /
• v.7 no.3
• /
• pp.195-223
• /
• 2017

This paper presents 3D numerical simulations of a Free Standing Hybrid Riser under Vortex Induced Vibration, with prescribed motion on the top to replace the motion of the buoyancy can. The model is calculated using a fully implicit discretization scheme. The flow field around the riser is computed by solving the Navier-Stokes equations numerically. The fluid domain is discretized using the overset grid approach. Grid points in near-wall regions of riser are of high resolution, while far field flow is in relatively coarse grid. Fluid-structure interaction is accomplished by communication between fluid solver and riser motion solver. Simulation is based on previous experimental data. Two cases are studied with different current speeds, where the motion of the buoyancy can is approximated to a 'banana' shape. A fully three-dimensional CFD approach for VIV simulation for a top side moving Riser has been presented. This paper also presents a simulation of a riser connected to a platform under harmonic regular waves.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.783-786
• /
• 2002

The unsteady turbulent flow characteristics of NACA4406 airfoil at close proximity to the pound are numerically investigated, especially focused on the noise generation mechanism near the blunt trailing edge. The unsteady two-dimensional compressible Wavier-Stokes equations with a Spalart-Allmaras turbulence closure model are solved by the 6th-order compact scheme and the 4th-order Runge-Kutta scheme. The computation shows a noise generation by a feedback mechanism at the blunt tailing edge, where the acoustic-fluidic coupling occurs between the wall-reflected sound waves and the periodically disturbed turbulent shear layer.