This paper briefly introduces the work done up to 1998 during the past twenty years for numerical modeling of ocean process focussing on the neighbouring seas of Korean Peninsula. Modeling of global ocean dynamics has also been performed as a pathway to understand the regional ocean dynamics. The ocean simulation produces a vast amount of multidimensional multivariate dataset therefore adoption of scientific and technical visualization techniques were essential to properly understand the physics involved.

A numerical prediction method has been proposed to predict non-linear free surface oscillation in a three-dimensional container. The fluid motions are numerically predicted with Navier-Stokes equations discretized in a Lagrangian scheme with sufficient numerical accuracy. The profile of a free surface is precisely represented with three-dimensional body-fitted coordinates (BFC), which are regenerated in each computational step on the basis of the arbitrary Lagrangian-Eulerian (ALE) formulation. In order to confirm the reliability of the computational method, it was firstly applied to three-dimensional flows within complicated-shaped rigid boundaries, such as curved pipes and ducts. Than it was applied to benchmark computations related to free surface oscillations. Following these basic verifications, non-linear sloshings in a cylindrical tank and transitions from sloshing to swirling motions were numerically predicted. Throughout these computations, the applicability of the present computational method has been confirmed and some of the predicted free surface motions were visualized as sequential images and animations to understand their dynamic futures.

A mutigrid convergence acceleration technique is presented for computing hypersonic inviscid and viscous flows in equilibrium state. The governing equations are solved using an explicit Runge-Kutta method. Curve fitting data in NASA Reference Publication 1181, 1260 are used to calculate equilibrium properties. In order to ensure stability, damped prolongation and modified implicit residual smoothing are proposed. Blunt body test cases are presented to demonstrate the robustness and the efficiency in performance of the proposed methods

A two-dimensional Navier-Stokes code based on AUSMPW+ scheme has been developed to simulate the hypersonic flowfield of hypersonic shock-shock interaction. AUSMPW+ scheme is a new hybrid flux splitting scheme, which is improved by introducing pressure-based weight functions to eliminate the typical drawbacks of AUSM-type schemes, such as non-monotone pressure solutions. To study the real gas effects, three different gas models are taken into account in this paper: perfect gas, equilibrium flow and nonequilibrium flow. It has been investigated how each gas model influences on the peak surface loading, such as wall pressure and wall heat transfer, and unsteady flowfield structure in the region of shock-shock interaction. With the results, the value of peak pressure is not sensitive to the real gas effects nor to the wall catalyticity. However, the value of peak heat transfer rates is affected by the real gas effects and the wall catalyticity. The structure of the flowfield also changes drastically in the presence of real gas effects.

The development of a parallelized aerodynamic simulation process involving moving bodies is presented. The implementation of this process is demonstrated using a fully systemized Chimera methodology for steady and unsteady problems. This methodology consist of a Chimera hole-cutting, a new cut-paste algorithm for optimal mesh. interface generation and a two-step search method for donor cell identification. It is fully automated and requires minimal user input. All procedures of the Chimera technique are parallelized on the Cray T3E using the MPI library. Two and three-dimensional examples are chosen to demonstate the effectiveness and parallel performance of this procedure.

This paper deals with the numerical simulation of the behavior of single bubble rising near the free surface. Volume fraction of fluid (VOF) method with continuum surface force (CSF) model, the well known method for two phase flow simulation is adopted. A bubble of spherical shape positioned beneath the free surface is assumed at the initial stage. The difference according to the fluid properties of surrounding medium is examined. Simulation results are depicted and explained with the time history of bubble shape, velocity field and vorticity distribution.

The two different approaches for solving the nonlinear ship wave problem are disscussed in the present paper. The first one is based on a panel method which neglects the viscous effects. Another one is based on a finite volume method which solves RANS equations. The present paper has been focused on the advantages/disadvantages of the above two methods. The developed methods are applied to calculating the flow around Wigley hull & Series 60 hull to validate the performance of the present nonlinear method, Although the two methods employ quite different numerical approaches, the calculated wave patterns from the both methods show good agreements with the experiments.

A numerical simulation method has been under development for solving turbulent flows around a ship model in maneuvering motion using the Reynolds Averaged Navier-Stokes equations. The method used second-order finite differences, collocated grids, pressure-Poisson equation and four-stage Runge-Kutta scheme as key components of the solution method. A modified Baldwin-Lomax model is used for the turbulence closure. This paper presents a preliminary result of the computational study on turbulent flows past a ship model in drift motion. Calculations are carried out for a Series 60 $C_B=0.6$ ship model, for which detailed experimental data are available. The results of the present calculations are compared with the experimental data for hydrodynamic forces acting on the model as well as velocity distributions at longitudinal sections. Only fair agreements has been achieved. The computational results show the complex asymmetrical shear flow patterns including three-dimensional separations followed by formation of bilge vortices both in bow and stern regions.

극초음속 유동장의 정확한 해석을 위해 AUSMPW+ 수치기법과 충격파 포착시 생기는 수치오차를 제거하기 위해 충격파 정렬 기법(Shock-Aligned Grid Technique)을 개발하였다. AUSMPW+ 수치기법은 자체 수치점성이 적은 수치기법으로 점성 경계층 계산시 정확한 계산결과를 보여주며 기존의 AUSM 계열이 가지는 문제점인 물성치의 진동 현상을 제거한 수치기법이다. 원통형과 무딘 물체 주위의 극초음속 유동장 해석을 통해 공력이 진동현상 없이 정확하게 계산됨을 확인하였다. 그리고 충격파 정렬 기법의 특성을 파악하기 위해 충격파 반사문제와 충격파-충격파 상호작용 문제를 해석하여 수치오차 없이 충격파를 포착할 수 있음을 보였다. 또한 화학적 평형 비평형 유동 영역까지 충격파 정렬 격자 기법을 확장하였다.

The shock wave diffracted behind a finite wedge is partially scattered after interacting with a starting vortex originated from the sharp vertex of the wedge. The shock is divided into the accelerated and decelerated shocks. The decelerated shock then interacts with the small vortexlets brought about by the vortex instability, producing weak compression waves. The shock-shock interaction produces Mach stems. Through this successive process, the shock attenuated. In this study, these complicated shock phenomena are computed using Euler equations and compared with experimental results obtained by the authors.

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 at the center axis and on the sonic surface in the shear layer, and the pressure difference is resolved across these waves interacted with the turbulence mixing layer. In this paper, the axi-symmetric Navier-Stokes equation has been used with two equation $k-{\varepsilon}$ turbulence closure model. The second order TVD scheme with flux limiters, based on the flux vector split by the smooth eigenvalue split, has been used to capture internal shocks and other discontinuities. The correction term for the compressible flow and the damping function are used in the turbulence model. 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.

An implicit turbulent flow solver is developed for 2-D unstructured hybrid meshes. Spatial discretization is accomplished by a cell-centered finite volume formulation using an upwind flux differencing. Time is advanced by an implicit backward Euler time stepping scheme. Flow turbulence effects are modeled by the Spalart-Allmaras one equation model, which is coupled with wall function. The numerical method is applied for flows on a flat plate, the NACA 0012 airfoil, and the Douglas 3 element airfoil. The results are compared with experimental data.

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

In this article, combination of the FAS-FMG multi-grid method and the Krylov subspace method was presented in solving two dimensional driven-cavity flows. Three algorithms of the Krylov subspace method, CG, CGSTAB(Bi-CG Stabilized) and GMRES method were tested with MILU preconditioner. As a smoother of the pressure correction equation, the MILU-CG is recommended rather than MILU-GMRES(k) or MILU-CGSTAB, since the MILU-GMRES(k) preconditioner has too much computation on the coarse grid compared to the MILU-CG one. As for the momentum equation, relatively cheap smoother like SIP solver may be sufficient.

본 연구에서는 3차원, 비압축성 유동에 적용 가능한 순수 벡터 포텐셜 방법 (Pure Vector Potential Method, PVPM for short)의 정식화를 수행하였고, 2차원의 경우에 적용 가능한 수치적 방법과 경계조건의 부여법을 상세히 설명하였다. 이 방법의 타당성을 검증하기 위하여 급작스럽게 출발하는 벽구동 웅덩이 내부 유동과 엄밀해가 존재하는 Stokes의 제 1문제를 수치해석하였다. 계산 결과들은 이 두 유동의 중요한 거동을 매우 잘 묘사하였다.

A flowfields around a NACA0012 airfoil pitching about a 1/4 chord and plunging in vertical displacement are analyzed by solving two-dimensional compressible Navier-Stokes equations. A steady solution was solved first as a validation of the code used and the results were compared with experimental data. Then as a unsteady case, the oscillatory airfoil was solved to compare the results with experimental data. Oscillating rate of pitching and plunging motion was set to have analogy and the magnitude of plunging was set using the magnitude of pitching angle of attack. Finally combined pitching and plunging motion was solved to show the effect of 2 different types of oscillating motion of the airfoil.

박리를 수반하는 진동하는 이차원 익형 주위의 비정상 유동장에 대해 SST 난류 모델을 이용하여 해석을 수행하였다. SST 모델은 정상 유동장 해석에서 기존의 난류 모델에 비해 우수한 성능을 보인다고 알려져 있으나 큰 박리영역에서 공력 계수의 진동 현상을 보이는 등 비정상 유동장 해석에 문제점을 보이고 있다. 본 논문에서는 이러한 공력 계수의 진동 현상이 SST 모델을 이차원으로 확장하는 과정에서 발생한 것임을 밝히고 이에 대한 보완을 통하여 수정된 SST 모델을 제시하고자 한다. SST 모델의 기본이 되었던 BSL 모델 및 SST 모델, 수정된 SST 모델을 사용하여 정상 유동장과 비정상 유동장 해석을 수행하여 각 모델의 난류 유동장 해석 특성을 비교하고 이를 통하여 수정된 SST 모델이 박리를 수반하는 비정상 유동장 해석에서 원래의 SST 모델에 비해 향상된 결과를 나타남을 알 수 있었다.

The unsteady transonic viscous flow has been analyzed over an oscillatory airfoil. The CSCM(Conservative Supra Characteristic Method) upwind flux difference splitting method and the iterative time marching scheme having first order accuracy in time and second to third order accuracy in space was applied on dynamic meshes. A steady flow field of Mach number 0.7 has been calculated for the verification of unsteady algorithm. The time-accurate unsteady calculations have been done for NACA 0012 airfoil oscillating around quarter chord about freestream Mach number 0.6 on dynamic meshes. The results have been compared with the AGARD Case 3 experimental data. The periodic characteristics have been compared with the experimental results.

This paper describes the flow analysis of the automobile with crosswind effects of $15^{\circ},\;30^{\circ}\;and\;45^{\circ}$ of yaw angle. The governing equations of the 3-D unsteady incompressible Navier-Stokes equations are solved by the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror. The A- and C-pillar vortex and other flow phenomena around the ground vehicle are evidently shown.

The three-dimensional turbulent cascade flows with and without endwall fences are numerically investigated by solving the incompressible Navier-Stokes equations with a high-Reynolds number $k-{\varepsilon}$ turbulence closure model. A projection method based algorithm is used in the finite-volume formulation, with the second order upwind-differencing scheme for the convective terms. First, assessments on accuracy of the present method are made by comparing the static pressure distributions at the mid-span of the cascade with measured data, and also by confirming the experimental observations on the choice of an optimal fence height for the secondary flow control. In understanding the three-dimensional nature of the secondary flow in turbine cascade, the limiting streamline patterns and the static pressure contours at the suction surface of the blade as well as on the cascade endwall are employed to visualize the effectiveness of the endwall fence for the secondary flow control. Analysis on the streamwise vorticity contour maps along the cascade with the three-dimensional representation of their iso-surfaces reveals the strucuture of the complicated vortical flow in the turbine cascade with endwall fence, and also leads to an understanding on formation of the counter-rotating streamwise vortex over the endwall fence, in explaining the mechanisms of controlling the secondary flow and also for the proper selection of an optimal fence height.

Numerical study has been peformed to develop a computer code for the design & optimal operating condition of SNCR(Selective Noncatalytic Reactor) for NOx reduction together with the analysis of the performance of vortex tube. Especially for the SNCR of the scale of industrial boiler, the required mixing and residence time of $NH_3$ solution was successfully tested numerically by the implementation of some baffle setups in a combustor. The introduction of interesting phenomena of vortex tube and similar system is made together with a theoretical hypothesis and simple cold flow simulation for the flow field analysis.

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

The airfoil design optimization procedures based on the Navier-Stokes equations were developed, This procedure enables more realistic and practical transonic airfoil designs. The modified Hicks-Henne functions were used to generate the shape of airfoils. Five Hick-Henne functions were used to design upper surface of airfoil only. To enhance the ability of Hick-Henne function to generate various airfoil shape with limited number of functions, the positions of control points were adjusted through optimization procedure. The design procedure was applied to the single-point design for the drag minimization problem with lift and area constraints. The result shows the capability of the procedure to generate much realistic airfoils with very small drag-creep in the low transonic regime. This is mainly due to the viscosity effect of Navier-Stokes flow analysis. However, in the higher transonic range tile drag-creep appears. The multi-point design is shown to be an effective way to avoid the drag-creep and improve off-design performance which is very similar in the Euler design.

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

The present study proposes modified temperature-dependent non-Newtonian viscosity model and investigates flow characters and heat transfer enhancement of the viscoelastic non-Newtonian fluid in a 2:1 rectangular duct. The proposed modified temperature dependent viscosity model has non-zero value near the high temperature and high shear rate region while on the existing viscosity models have zero value. Two versions of thermal boundary conditions involving difference combination of heated walls and adiabatic walls are analyzed in this study. The combined effect of temperature dependent viscosity, buoyancy, and secondary flow caused by second normal stress difference are ail considered. The Reiner-Rivlin model is adopted as a viscoelastic fluid model to simulate the secondary flow caused by second normal stress difference. Calculated Nusselt numbers by the modified temperature-dependent viscosity model gives under prediction than the existing temperature-dependent viscosity model in the regions of thermally developed with same secondary normal stress difference coefficients with experimental results in the regions of thermally developed. The heat transfer enhancement of the viscoelastic fluid in a 2:1 rectangular duct is highly dependent on the secondary flow caused by the magnitude of second normal stress difference.

Three dimensional numerical analysis for the turbulent reactive flow and radiative heat transfer in the walking beam type of a reheat furnace in POSCO has been carried out by the industrial code FLUENT. Computations an based on the conservation equations of mass, momentum, energy and species with the $k-{\varepsilon}$ turbulence model and mixture fraction/PDF(Probability Density Function) approach for the combustion rate. Radiative heat transfer is computed by the discrete ordinates radiation model in combination with the weighted-sum-of-gray-gas model for the absorption coefficient of gas medium. The predicted temperture distribution in the reheat furnace and energy flow fractions are in reasonable agreement with the measurement data.

Safety Depressurization System of the Korean Next Generation Reactor prevents the Reactor Coolant System from over-pressurization by discharging the coolant with high pressure and temperature into the In-containment Refueling Water Storage Tank(IRWST) during an accident. If temperature in the IRWST rises above the temperature limit of $200\;^{\circ}F$ due to the discharged coolant, an unstable steam condensation may occur and cause large load on the IRWST wall. To investigate whether this condition can be reached or not for the design basis accident, the flow and temperature distributions of water in the IRWST wire calculated by using CFX 4.2 computer code. The results show that the local water temperature does not exceeds the temperature limit within the transient time of 5 seconds.

The $LiBr-H_{2}O$ absorption process on a horizontal tube has been analyzed using the numerical method which incorporates the fully elliptic Navier-Stokes equations for the momentum equations, the energy and mass-diffusion equations. On a staggered grid, the SIMPLER algorithm with the QUICK scheme is used to solve these equations along with the MAC method for the free surface tracking. With the assumption that the absorbent is linear, calculations have been made for various inlet temperature and flow-rate conditions. The detailed results of the parametric study, such as the temperature, concentration, absorption mass flux and wall heat flux distributions are presented. The self-sustained feature of the absorption process is clearly elaborated. The analyses have also been carried out for multiple tube arrangement and the results show that the absorption rate converges after a few tube rows.