Underwater explosion properties for TNT, an ideal explosive, and DXD-04, a nonideal explosive, were numerically modeled with a one-dimensional Lagrangian hydrodynamic code. The equation of state parameters for detonation products for TNT and DXD-04 were obtained from the BKW code, assuming complete reaction. Burn of TNT was modeled by using the Chapman-Jouguet(CJ) volume burn technique, a programmed-burn technique, assuming instantaneous detonation reaction. Burn of DXD-04 was modeled by using the same technique and by using the reaction rate calibrated from two-dimensional steady-state detonation experiments. The calculations for TNT reproduced the experimental peak pressure of the shock wave propagating through water with an error of $3.0\%$ and the experimental oscillation period of the bubble formed of detonation products with an error of $2.3\%$. For DXD-04, the CJ volume burn technique could not reproduce the experimental observations. When the reaction rate calibrated from two-dimensional steady-state detonation experimental data, the calculated peak pressure was slightly higher by $7.3\%$ than the experimental data, but the calculated shock profile was in good agreement. The bubble period was reproduced with an error of $1.8\%$. These results demonstrated that underwater explosion properties for an ideal explosive can be predicted by using a programmed burn technique, and that, however, those for a nonideal explosive can be predicted only when a well-calibrated reaction rate is used.

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.

An efficient 3-dimensional compressible solver is developed using the second-order upwind TVD scheme and the multigrid diagonalized ADI method. The multigrid method is improved so that the present DADI algorithm obtains better convergence rates. Results are computed on Cray C90 computer for transonic unsaperated flows past ONERA-M6 wing to demonstrate the accuracy and efficiency. The results show good agreement with experimetal data. A reduction of four orders of residual for 3-dimensional transonic flow is obtained about 99 seconds.

Time marching algorithms applied to compressible Navier-Stokes equation have a convergence problem at low Mach number. It is mainly due to the eigenvalue stiffness and pressure singularity as Mach number approaches to zero. Among the several methods to overcome the shortcomings of time marching scheme, time derivative preconditioning method have been used successfully. In this numerical analysis, we adopted a preconditioner of K.H. Chen and developed a two-dimensional, axisymmetric Navier-Stokes program. The steady state driven cavity flow and backward facing step flow problems were computed to confirm the accuracy and the robustness of preconditioned algorithm for low Mach number flows. And the transonic and supersonic flows insice the JPL axisymmetric nozzle internal flow is exampled to investigate the effects of preconditioning at high Mach number flow regime. Test results showed excellent agreement with the experimental data.

The development of unsteady three-dimensional incompressible viscous solver based on unsteady physical curvilinear coordinate system is presented. A 12-point finite analytic scheme based on local uniform grid spacing is extended for nonuniform grid spacing. The formulation of a condition is suggested to avoid the oscillation of the series summations produced by the application of the method of separation of variables. SIMPLER and pressure Poisson equation techniques are used for solving a velocity-pressure coupled problem. The matrix is solved using the Generalized Minimal RESidual (GMRES) method to enhance the convergence rate of unsteady flow solver and the Kinematic boundary condition of a free surface flow. It is demonstrated that the numerical solutions of these equations are less mesh sensitive.

3-dimensional Euler solver is parallelized. The spatial discretization method is the 2nd order TVD scheme and DADI method with multigrid is used as a time integration. In order to parallelize this solver, the domain decomposition method with overlapped grid and message passing techniques are used. The informations on the each inter-processor bound-aries are communicated with MPI library. Finally, the parallel performance repsented by calculating the ONERA M6 wing at transonic flow condition using CRAY T3E and C90.

An elliptic grid generation scheme using Laplace's equations guarantees the resulting grids to be crossing-free as a result of maximum principle in its analytic form. Numerical results, however, often show the grid lines overlapping each other or crossing the boundaries, especially for very sharp convex corners. The cause of this problem is investigated, and it is found that this problem can be handled by properly modifying the coefficients of transformed Laplace's equations in the computational domain.

An inverse method for the subsonic and transonic airfoil design was developed using the Euler equations. Two testcases were performed. One was a design of the supercritical airfoil, aiming to be used for the Korean mid-sized (100 passengers class) transport aircraft. The other was the design of an airfoil showing a good cruising performance (L/D ratio) in the high subsonic/transonic flow regimes. These testcases demonstrated the efficiency and the robustness of the developed method.

Reducing drag of vehicles are the main concern for the body shape designers in order to lower fuel consumption rate and to aid the driving stability. The drag of bluff bodies like transportation vehicles is mostly pressure drag due to the flow separation, which can minimized by controlling the location and size of the separation bubble. In the present study, the TURBO-3D code is incorporated with optimal algorithm based on analytical approximation method to obtain optimal afterbody shape of the MIRA Model corresponding to the lowest drag coefficient. For this purpose three mutually independent afterbody angles are chosen as design variables, while the drag coefficient is chosen as an objective function. It is demonstrated in the present study that an optimal body shape having lowest drag coefficient which is about $6\%$ lower than that of the original shape has been successfully obtained within number of iterations of the optimal design loop.

A new advanced integral reactor of 330 MWt capacity named SMART(System-integrated Modular Advanced ReacTor) is currently under development at KAERI(Korea Atomic Energy Research Institute). One of the major design features of the integral reactor is locating the steam generators(SG) inside reactor vessel and eliminating the possibility of LB LOCA(large Break Loss of Coolant Accident). Orifices are fitted at the low part of steam generator cassette to stabilize and balance coolant flow distribution in the MCP (Main Circulation Pump) pressure header. A sensitivity analysis is conducted to determine the optimum orifice size using computer code 'CFX'.

There are many moving-body problems to be solved, the solution of which necessary for proper design of flight vehicles in aerospace industry. Since a body moves relative to other bodies in the category of these problems, difficulty arises regarding both generation of computational grid around the body in motion and conservation of flow properties in the moving grid system. A few example could be store separation from the aircraft and relative vibration of multiple bodies in the high-speed flow passage. In this paper we report on the progress made in computing moving-body aerodynamics related with sabot separation characteristics. Conservative overlapping grid together with cell-merging-unmerging technique is used to solve the Euler equations for a body in high-speed motion. Carbuncle errors has to be removed before we obtain physically adequate solution. Two-dimensional application is reported here.

본 논문은 2차원 해면효과의 수치계산 결과를 정리하였다. 지면으로부터의 높이변화에 따른 점성유동장을 계산하기 위하여 지배방정식으로는 비압축성 RANS방정식을, 시간에 대하여 서는 음해법으로 프로그램을 구성하였다. 압력항은 가상압축성을 도입 4차 수치확산항을 추가하는 것에 의해 계산하였으며, 높은 레이놀즈수에서의 효과적인 계산을 위해 Baldwin-Lomax 난류 모델을 도입하였다. 해면효과가 없는 무한유중에서의 NACA-0012단면 계산결과를 실험데이터와 비교하는 것에 의해 프로그램의 타당성을 확인하였다. NACA-6409와 두께비 $4.6\%$의 날개에 대하여 해면효과를 고려한 계산을 수행하였다. 높이의 변화에 따라 계산된 무차원계수, 압력 및 속도분포는 해면효과의 특성을 잘 보여주고 있다.

This paper describes the analysis of flow and aeroacoustic field around a car-like body. The governing equations, 3-D unsteady incompressible Navier-Stokes equations, are solved with the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror, After the flow field analysis has been converged, the aerodynamic noise analysis of the side-view mirror has been performed by solving Ffowcs Williams and Hawkings equation. From the present numerical simulation, the A- and C-pillar vortex are evidently shown and the aerodynamic noise level induced by the side-view mirror is predicted to about 100dB.

An axisymmetric Navier-Stokes procedure has been developed to analyze the pulse jet flow in a ceramic filter unit for the dust dislodging process. Using Baldwin-Lomax turbulence model as a closure relationship, the SIAF(Scalar Implicit Approximate Factorization) algorithm together with the ${\delta}^k-Correction$ iterative time marching scheme is adopted to solve the unsteady compressible Navier-Stokes equations. After some validation tests, the code has been applied to solve the pulse jet flow and examine the effects of geometry and reservoir pressure condition on the pressure level inside the filter unit. To avoid dealing with the uncertainty of such factors as the cohesion of the collected dust and the adhesion of the dust to the medium and also to simplify the analysis, the filter wall is assumed to be impermeable. The results for various test cases are presented.

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.

A computer code has been developed for an analysis of conduction and convection heat transfer on a personal computer. FVM based on unstructured mesh has been employed for triangular and quadrilateral element. All variable of p, u, v, T, k and ${\varepsilon}$ has been defined on a node (not on a center of element). The code possesses the pre-and post-processor for itself to provide user-friendly interface.

In this numerical experiment, the preconditioned compressible Navier-Stokes equation is tested to analyze the laminar spray combustion. Sprayed flow field is formulated by Eulerian-Lagrangian system for the gas and liquid phases each. DSF(Deterministic Separated Flow) model was adopted for the sprays with the vortex model to describe transients of individual droplet heating. Simplified single global reaction model approximates methanol-air reaction with and without disk flame holder. The equation system is discretized by finite difference technique and time integrated by LU-SGS. Due to greatly simplified chemical reaction mechanism and the lack of experimental evidences, most of the efforts were devoted to show the applicability and robustness of preconditioned compressible flow calculation algorithm. Computation results in qualitatively reasonable combusting flow field, hence it is believed that further refinement are required to produce quantitatively accurate solutions.

The BGK equation, which is the kinetic model equation of Boltzmann equation, is solved using FDDO(finite difference with the discrete-ordinate method) to compute the rarefied flow of monatomic gas. Using reduced velocity distribution and discrete ordinate method, the scalar equation is transformed into a system of hyperbolic equations. High resolution ENO(Essentially Non-Oscillatory) scheme based on Harten-Yee's MFA(Modified Flux Approach) method with Strang-type explicit time integration is applied to solve the system equations. The calculated results are well compared with the experimental density field of NACA0012 airfoil, validating the developed computer code. Next. the computed results of circular cylinder flow for various Knudsen numbers are compared with the DSMC(Direct Simulation Monte Carlo) results by Vogenitz et al. The present scheme is found to be useful and efficient far the analysis of two-dimensional rarefied gas flows, especially in the transitional flow regime, when compared with the DSMC method.

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

This describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows of fluid behaviour with free-surface. The elliptic differential equations governing the flows have been linearized by means of finite-difference approximations, and the resulting equations have been solved via a fully-implicit iterative method. The free-surface is defined by the motion of a set of marker particles and interface behaviour was investigated by way of a 'Lagrangian' technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions or experimental results from the literature. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

A computer code has been developed for the analysis of 2D viscous flow with free surface. VOF method and higher order upwind scheme have been employed for the accurate prediction of free surface motion. Surface tension effect and axisymmetric flow can be computed by the present code.

A numerical analysis of axisymetric backward facing step main nozzle flow in air jet loom has been accomplished. To obtain basic design data for an optimum main nozzle for an air-jet loom and to predict the transonic/supersonic flow, a characteristic based upwind flux difference splitting compressible Navier-Stokes method has been used. The wall static pressure of the main nozzle and the flow velocity changes in the nozzle tube were analyzed by changing air tank pressures and acceleration tube lengths. The flow inside the nozzle experiences double choking one at the needle tip and the other at the acceleration tube exit at tank pressures over $4kg_f/cm^2$. The tank pressure $P_t$ leading to the critical condition depends on the acceleration tube length; i.e, $P_t$ is higher for longer acceleration tubes. The $P_t$ value required to bring the acceleration tube exit to the critical condition is nearly constant regardless of acceleration tube length. The round needle tip shape might lead to less total pressure loss when compared with step shape.

Numerical study of three-dimensional turbulent flow in a forward curved centrifugal fan is presented. Standard $k-{\varepsilon}$ turbulence model and non-orthogonal curvilinear coordinates are used to consider the turbulent flow field and complex geometry. Finite Volume approach is adopted for discretization scheme and structured grid system is used to help convergence. Multiblock grid system is used for flow field and divided into five domains that are inlet, outlet, impeller, tip clearance and scroll. It is assumed that the flow field is steady state and incompressible. This numerical work is performed with commercial CFD-ACE code developed by CFD Research Corporation, and the results are compared wi th the experimental data

A numerical study based on the three-dimensional Reynolds averaged Navier-Stokes equations is presented to analyze the transonic flowfield through two-stage axial compressor. Explicit four-step Runge-Kutta scheme is used for solution algorithm, and local time step and implicit residual averaging are introduced for enhancing the convergency. Artificial dissipation model is adopted to assure the stability of solution. The solver is coupled with Baldwin-Lomax model to describe turbulence. To avoid calculating the unsteady flow, a mixing process is modeled at a station between rotating and stationary blade rows. Results show a variety of important physical phenomena. Comparison of the flowfields with and without tip clearance shows that the effect is considerable in this flowfield. Comparisons with experimental data carried out to validate the calculational results show reasonable agreements. Some remedies are also suggested to improve the revealed problems.

The numerical methodology for computing the impact forces of the water entry bodies has been developed. The present method assumed the impact occurs within a very short time interval and the viscous effects do not have time enough to play a significant role in the impact forces, that is, the flow around a water-entry object was assumed as an inviscid potential flow and is solved by the source panel method. The elements fully submerged into the water are routinely treated, but the elements intersected with the effective planar free surface are redefined and reorganized to be amendable to the source panel method. To validate the present code, it has been applied to disk and ogive model and compared with experimental data. Good agreement has been obtained.