In order to reduce the excessive numerical dissipation which is induced when a grid system is not aligned with a discontinuity, a new spatial treatment of cell-interface fluxes is introduced. The M-AUSMPW+ in this paper has the formula that has an additional procedure of re-defining transferred properties at a cell-interface, based on AUSMPW+. The newly defined transferred property could eliminate numerical dissipation effectively in non-flow aligned grid system of multi-dimensional flows.

It is well known that preconditioning methods are efficient for convergence acceleration at compressible low Mach number flows. In this study, the original Euler equations and three preconditioners nondimensionalized differently are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as flux discretization and time integration respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning one produces Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

A numerical method for simulating two-phase flows including surface force is presented. The method is based on fractional step method of finite volume formulation and the interface is tracked with PLIC VOF method. In the CSF model, as color function, f, representing the location of interface varies steeply in the interface region, we need to use smoothed function f to get accurate unit normal and the curvature. Peskin kernel is used to get smoothed function f. A spherical drop in static equilibrium and three-dimensional merging of gas bubble are tested, resulting in the validation of this method

A great deal of exhaust gas inside a combustion room goes out through exhaust pipe but residual gas, is called 'Blow by gas', enters the crankcase through a small gap between the piston and the cylinder wall. Here, if the crankcase isn't vented, this causes many bad effects such as lubricant oil contamination, corrosion by that and crankcase explosion by rising pressure. So most automobiles are constituted with a PCV(Positive Crankcase Ventilation) system to prevent previous problems. PCV valve is the most important part in this ventilation system. When companies are manufacturing new cases, engineers are designing it depending on their experiments than theoretical knowledges. Much efforts and times are needed for new development. This study will show quantitative results to increase the possibilities.

The purpose of this paper is to present essential numerical data of FCV(Flow rate Control Valve) which is used in a projectile ejection system for a underwater vehicle. A commercial CFD code is applied to analyze the 3 dimensional viscous compressible flow field within the FCV as four cases of opening rate, 25, 50, 75, and $100\%$. The flow coefficient of each cases are mainly calculated. And other characteristics such as the location of shock wave and total temperature distribution are also determined.

A three dimensional numerical analysis and experimental measurements were made for incompressible flows in a globe valve with three different types of plug. Characteristics of pressure drop across the valve were investigated for each type of plug at various opening positions. Numerical simulation results show a good agreement with experimental data.

The ballast tank of a ship is a system that realizes the required shipping condition and controls the draft of a ship. The loading/unloading of the ballast tank is frequently operated during navigation and the accurate prediction of the loading/unloading time is very important. A numerical algorithm that predicts the loading/unloading time of the ballast tank has been developed and applied to the prediction of the loading/unloading time of the ballast tank with various piping systems. This algorithm can be useful in optimizing the ballast tank system in early design stage.

Liquid injection in a Venturi Scrubber creates great effect on the dust-collection efficiency and operation cost of venturi scrubbers. We have developed a model that can numerically simulate atomization of the liquid jet in the Venturi Scrubber. This simulation consists of models on liquid column, jet surface breakup, column fracture and secondary droplet breakup. These models have been embedded in the KIVA3-V code. We have calculated such parameters as the jet penetration, jet trajectory, droplet size, velocity field and the volume flux distribution. The results are compared with the experimental data in this paper.

Numerical analysis for pressure drag on boattail afterbodies have been studied by Mach number, boattail angle and length ratio of body diameter and base diameter using CFD-FASTRAN that the commercial external flow CFD code. The numerical results have been compared with the experimental data that have been shown pressure drag reduction and supersonic turbulent flow characteristics for boattail afterbodies. And the prediction equation tot boattail base drag has been made by the numerical results about Mach number and boattail configuration parameters.

Accurate Prediction of a supersonic missile base drag continues to defy even well-rounded CFD codes. In an effort to address the accuracy and predictability of the base drags, the influence of grid system and competitive turbulence models on the base drag is analyzed. Characteristics of some turbulence models is reviewed through incompressible turbulent flow over a flat plate, and performance for the base drag prediction of several turbulence models such as Baldwin-Lomax(B-L), Spalart-Allmaras(S-A), $\kappa-\epsilon$, $\kappa-\omega$ model is assessed. When compressibility correction is injected into the S-A model, prediction accuracy of the base drag is enhanced. The NSWC wind tunnel test data are utilized for comparison of CFD and semi-empirical codes on the accuracy of base drag predictability: they are about equal, but CFD tends to perform better. It is also found that, as angle of attack of a missile with control (ins increases, even the best CFD analysis tool we have lacks the accuracy needed for the base drag prediction.

A general purpose program NUFLEX for the analysis of 3-D heat/fluid flow in complex geometry with pre/post processor have been developed, which consists of a flow solver based on FVM and a dedicated pre/post processor. The program employs a general non-orthogonal grid system and solve laminar and turbulent (lows with standard and RNG $\kappa-\epsilon$ turbulence models. NUFLEX is capable of analysing two-phase flow with topologically complex interface, turbulent diffusion combustion, solidification problems and magnetic flow. For the purpose of verification of the program and testing the applicability, several practical problems are solved and compared with the available data. Comparison of the NUFLEX results with that by the STAR-CD program has been also made for the same flow configuration and grid structure.

To analyze the performance of Wind turbine of the drag force type, 3-D RANS equations were solved by the iterative time marching method on sliding multiblock grid system. The numerical flow simulations by changing blade number and pitch angle were carried out : blade number = 15, 20 circumferentially; pitch angle = $30^{\circ},\; 50^{\circ}$ radially. The torque coefficient was also calculated.

The design procedure for centrifugal blower with high inlet resistance is not presented yet. Overall fluid dynamic performance is estimated for comparison between the case of atmospheric inlet condition and the present model. Detail information between blades is prepared by using a commercial program, SCRYU-Tetra. A centrifugal blower with large inlet pressure is adopted in an air purifier having filtering devices. As the inlet residence increases the flow rate of the system is decreased. In parallel, outlet area of the system affects the performance of the system in the sense of flow balance. Consequently, the flow balance between the inlet and outlet becomes an important parameter for the design of the scroll casing for the centrifugal blower with high inlet pressure.

The performance of the centrifugal fan in a vacuum cleaner is affected by the hydraulic loss, such as the friction loss, the recirculation loss and the impact loss etc., Those losses depend on the rotational speed of the impeller, the inlet and exit widths, the relative flow angles to the blade, the number of the blades and the geometry of the shroud and the diffuser. These parameters are complicatedly interrelated, so the experimental means in analyzing the fans are rather limited. In the present study, the flow analysis are done numerically by changing the relevant fan parameters. A commercial code, STAR-CD, is used for the calculations. It is seen from the analyses that the computational results agree well with the experimental results. The results obtained can be used for the basic design of a centrifugal fan.

In the present study, an efficient and robust implicit operator for the LU-SGS method is proposed. Numerical experiments for supersonic flow are performed to demonstrate the performance of the proposed method.

As an alternative numerical method, the lattice Boltzmann method (LBM) is used to simulate a 2-dimensional pressure driven microchannel flow which comes from frequently in MEMS problems. The flow is assumed to be isothermal ideal gas flow. The flow field is calculated with various Knudsen numbers, pressure ratios and aspect ratios of the microchannel. The LBM can show the fundamental characteristics in microchannel flow such as velocity slip and nonlinear pressure drop.

In this paper an automated hole-finding method for overset grids is introduced which uses recursive octree-cell division. A graphic program which enables the user to do the hole-cutting with ease is also introduced. Using this program it was found that a proper combination of the level of octree division and vector calculation should be used for efficient and fast hole finding.

A numerical analysis was made for the unsteady flow fields of rotor system of a Tilt-Rotor aircraft in cruise mode. The Reynolds-averaged thin-layer Navier-Stokes equations were discretized by Roe's upwind differencing scheme and integrated in time by the LU-SGS algorithm. The computational domain of the rotor system was constructed by six multi-block Chimera grids. Simulated unsteady flow fields of rotating blades were studied in several different view points.

A numerical analysis was made to investigate the intensity diminution of a simple silencer for high pressure blast flow fields. Reynolds-Averaged Wavier-Stokes equations were solved for an axisymmetric computational domain constructed by multi block Chimera grids. A blast flow field without the silencer was also calculated to validate the present numerical method. The evolution of high pressure blast flow fields was observed by depicting calculated contours of pressure and Mach number. It was found that the tested silencer could achieve 76 percent intensity diminution.

The developing Smart UAV in KARI supposes high speed flight as like a conventional plane, as well as vertical takeoff and landing as like a helicopter. Therefore, the air intake system should be designed to provide the sufficient air flow to the engine and the maximum possible total pressure recovery at the engine intake screen over a wide range of flight conditions. For this purpose, we designed the intake system using a pilot type intake model and plenum chamber In this paper, we designed the intake model and analyzed the performance of designed intake system using the general-purpose commercial CFD code, CFD-ACE+ For 3-D calculation, we generated mesh using the unstructured gird and used $\kappa-\epsilon$ turbulence model. The analysis results of the total pressure variation and the velocity distribution was illustrated in this paper. The pressure recovery and distortion coefficient at a plane coincident with the compressor inlet were calculated and streamline variation through the intake system was investigated at the worst condition as well as the standard flight condition.

CFD simulation is peformed for 2D and 3D frisbees flying at 10m/s. For convenience of simulation, rotation of 3D model is not considered. CFD results show that pitching moment makes the nose down and holes at the leading and trailing edges improve the lift characteristics of the frisbee.

CFD simulation for Smart-UAV(TR-E2Sl) is performed to analyze its aerodynamic characteristics. Base geometry and several cases, decided by control surfaces being deflected, are simulated. To obtain the better lift characteristics, the elevator should be deflected between 10o and 20o with the incidence angle of the wing 1o.

CFD code that simulates stator-rotor interactions is developed applying parallel computing method. Modified Multi-Block Grid System which enhances perpendicularity in grid and is appropriate in parallel processing is introduced and Patched Algorithm is applied in sliding interface which is caused by movement of rotor. The experimental model in the turbo-machine is composed of 11 stators and 14 rotors. Analyses on two test cases which are one stator - one rotor model and three stators - four rotors model are performed. The results of the two cases have been compared with the experimental test data.

The effect of a continuous blowing or suction on an oscillating 2-D NACA0012 airfoil was investigated numerically for the dynamic stall control. The influence of control parameter variation was also studied in the view point of aerodynamic characteristics. The result showed that the blowing control kept a higher lift drag ratio before stall angle but the dynamic stall angle was not exceed to without control result. As the slot position was closer to leading edge, the positive control effect becomes greater. The stronger jet and the smaller jet angel made more favorable roles on the control performance. In the cases of the suction, the overall control features were similar to those of the blowing, but dynamic stall angle was increased, i.e. suction was more effective to control dynamic stall. It was also founded that the suction control was showed better control effect as the slot position moves to trail edge within thirty percentage of chord length. In the simulation for the jet strength and the jet angle control, the same tendencies were observed to those of blowing cases.

The preprocessor - solver - postprocessor software for 2D/Axisymmetric CSCM Upwind Flux Difference Splitting Navier-Stokes code has been developed for undergraduate educational purpose. This computational fluid dynamics (CFD) software allows students to setup, solve, visualize and control dynamically server for their own fluid problems via Internet. The preprocessor Is capable of generating geometry and grid, initial solution data and required solver control parameters. The postprocessor shows vector plot and contour plot with different options while residual plot shows root-mean-square (RMS) error history graphically and retrieves the data from solver interactively. Special feature of the preprocessor is grid generation part which is based on MFC/Visual C++ application and FORTRAN single block grid generator process. Many users can access solver via Internet from client computers and solve desired problems using locally installed pre- and postprocessor and remote powerful solver part.

In this study, three-dimensional numerical calculations are peformed to simulate the flow and heat transfer in helically coiled tube steam generator employing a commercial CFD (Computational Fluid Dynamics) code. The problem considered herein includes the boiling phase change flow of tube side fluid and the single-phase counter-current flow of shell side hot fluid transferring heat to the tube side flow thru the tube wall. Detailed investigations are performed for both shell-side and tube-side flow fields in terms of density and volume fractions of each phase of fluids as well as for the tube wall heat transfer field in terms of heat transfer coefficients.

Thermal control of satellite propellant tank is achieved by patch heaters enabled by thermostat's behavior. It is important to attach the thermostat on the appropriate position of the propellant tank. However its position cannot be given with exact numerics because tank is spherical. Actually the position for thermostat is designated in relevant drawing approximately, therby, the engineer practices depending on his own experience and intuition. The sensitivity analysis for the position of thermostat is performed such that the influence on the thermal behavior and control of tank is examined quantatively. When assembling tank module, the reasonable performance on the thermal control is believed with possible human errors if the uncertainty in the position of thermostat is not quite large.

The monopropellant hydrazine thrusters are widely used for the satellite on-board propulsion system fulfilling various missions in space. They have outstanding features caused by the nearly unlimited restart capability and the very high credibility. The sole monopropellant thruster used at precent in nation is MRE-1 that is a standard component of NASA. It can produce 4.45 N of nominal thrust. Due to the glowing complexity with a satellite mission, the needs for thrusters of the diverse performance are being increased. The numerical simulation could give useful information to develop a new type thruster instead of the experiments performed previously. Therefore it is critical to make a reliable computer code to prepare design change of a thruster. In this paper, the performance analysis and validation of the satellite monopropellant hydrazine thruster currently used is accomplished as the preliminary study to serve valuable data for future design change.

An analysis procedure for the MCFC channel flow has been developed to predict the fuel cell performance. The channel formed by the uniformly distributed trapezoidal supports is approximated by the porous medium that yields the equivalent pressure drop. The Wavier-Stokes, energy, and species equations are solved to obtain the velocity, temperature and concentration fields for a local current density which is computed from electrochemical correlations. The fuel cell characteristics, such as the temperature, pressure, mole concentration, voltage and current density distributions, are presented and discussed.

In this study the numerical analysis on staging flow with forward ejector is conducted. The forward ejector plays a vital role in staging, which jets out from aftbody. This staging environment needs careful flow analysis for securing staging safety Present study investigates the steady inviscid staging flow phenomena with variation of separation distance. The performance index is forebody base pressure coefficients. The three dominant flow phenomena are observed according to separation distance which could be told as impinging stage, cavity vortex dominancy stage, and pure base flow characteristics stage. Impinging stage shows high thrust for forebody as one might think. However, important point is that cavity vortex dominancy stage can be more favorable for separation than impinging stage as one simply think in certain separation distance.

The HANARO, a multi-purpose research reactor of 30 MWth open-tank-in-pool type, has been under normal operation since its initial criticality in February, 1995. Many experiments should be safely performed to activate the utilization of the HANARO. A flow simulated test facility has been developed for the verification of structural integrity of those experimental facilities prior to loading In the HANARO. This test facility is composed of three major parts; a half-core structure assembly, flow circulation system and support system. The half-core structure assembly is composed of plenum, grid plate, core channel with flow tubes, chimney and dummy pool. The flow channels are to be filled with flow orifices to simulate similar flow characteristics to the HANARO. This paper describes an analysis of the flow distribution of the cote channel and compares with the test results. As results, the analysis showed similar flow characteristics compared with those in the test results.

This paper describes a supersonic separation of air-launching rocket from the mother plane. Three dimensional Euler equations were numerically solved to analyze steady/unsteady state fluid flows. To solve the Euler equations, named CFD-FASTRAN that is commercial computation code. The results of simulation clearly demonstrate effect of shock-expansion wave interaction between the rocket and the mother plane. Moreover, important influential factors at separating stage of the rocket were extracted with a comprehensive analysis. Finally, from the consideration of supersonic-separation, a guideline to safety-separation is given to the design of supersonic air-launching rocket.

A finite volume numerical model is developed for simulating non-equilibrium electroosmotic flow in micro- and nanochannels. The Guoy-Chapman model is adopted to compute the flow and electric potential. The Nernst-Planck equation is employed to trace unsteady transports of ionic species, i.e., time-dependent net charge density. A new set of boundary conditions based on surface charge density are designed rather than using the conventionally-employed zeta potential. A few issues for an efficient computation of electroosmotic flows are discussed. Representative computational examples are given to illustrate the robustness of the numerical model.

Smart UAV, which adopting tiltrotor aircraft concept, requires long endurance and high speed capability simultaneously These two contradictable flight performances are hard to meet with single wing concept and inevitably require the operation of flap system which should reveal optimal performance for each flight mode. In order to design SUAV flaperon satisfying the performance requirement, various configurations are generated and their aerodynamic performances are analyzed using numerical flow computations around flaps. Considering aerodynamic performance and manufacturing simplicity, a final flap configuration is selected.

In this study, the multidisciplinary aerodynamic-structural optimal design is carried out for the supersonic fighter wing. Through the aeroelastic analyses of the various candidate wings, the aerodynamic and structural performances are calculated such as the lift coefficient, the drag coefficient and the deformation of the wing. In general, the supersonic fighter is maneuvered under the various flight conditions and those conditions must be considered all together during the design process. The multi-point design, therefore, is deemed essential. For this purpose, supersonic dash, long cruise range and high angle of attack maneuver are selected as representative design points. Based on the calculated performances of the candidate wings, the response surfaces for the objectives and constraints are generated and the supersonic fighter wing is designed for better aerodynamic performances and less weights than the baseline. At each design point, the single-point design is performed to obtain better performances. Finally, the multi-point design is performed to improve the aerodynamic and structural performances for all design points. The optimization results of the multi-point design are compared with those of the single-point designs and analyzed in detail.

The development of a tilting train with construction of electric line on the conventional railway is required for speed-up on the conventional railway with many curving sections. For development of tilting train, the study and development of the tilting system and tilting bogie having the different mechanism with a general high speed train will play a main role for improving the technology in the field of Korean railway The study and development of the pantograph tilting mechanism in order to keep a good contact behavior between a pantograph and a contact wire by tilting a pantograph on the opposite direction of the vehicle tilting direction. In this study, we analyzed the aerodynamic characteristic of a developing pantograph on the tilting train and obtained the contact force with catenary by aerodynamic lift force by the aerodynamic analysis. We also performed the numerical analysis for design the device controlling lift force on a pantograph. From the aerodynamic simulation and parameter study for a device to control the lift force, we will suggest the various shape and the optimal shape of it corresponding to a developing tilting pantograph. The Fluent software is used for the calculation of flow profile in this study.

In gradient based optimization methods, the finite differencing which uses small perturbations in the design variables has been used to calculate the sensitivity. Recently, the automatic differentiation has been widely studied to calculate the function value and the sensitivities simultaneously. In this paper, the applicability of the automatic differentiation In the aerodynamic design optimization is studied. ADIFOR and TAPENADE are used to generate the codes which give the function value and the sensitivities for 2D compressible inviscid flows.

The conventional reliability based design optimization(RBDO) methods require high computational cost compared with the deterministic design optimization(DO) methods. To overcome the computational inefficiency of RBDO, single loop methods have been proposed. These need less function calls than that of RBDO but much more than that of DO. In this study, the approximate reliability method is proposed that the computational requirement is nearly the same as DO and the reliability accuracy is good compared with that of RBDO. Using this method, the 3-D viscous aerodynamic shape design optimization with uncertainty is performed very efficiently.

Numerical optimization method of long endurance airfoil has been performed with a RSM(Response Surface Method) for smart UAV wing design. For the base line airfoil, NACA 64621 airfoil was selected and optimized to satisfy long endurance condition for smart UAV Aerodynamic coefficients required for RSM are obtained by using 2-D Navier-Stokes solver with Spalart-Allmaras turbulence model. The optimized airfoil showed increased maximum lift and endurance factors together with reasonable thickness ratio.

A design optimization of Sphere-Cone blunt nose hypersonic flight vehicle has been studied by using upwind Navier-Stokes method and numerical optimization method. Heat transfer coefficient and drag coefficient are selected as objective function or design constraint. Control points of Bezier curve are considered as design variable.

In this study, a numerical analysis has been performed for a three-dimensional pulsatile blood flow associated with the elastic blood vessel and curved bileaflet for multiple cycles in terms of fluid-structure interaction. Here, blood has been assumed as a Newtonian, incompressible fluid. Pressure profiles have been used as boundary conditions at the ventricle and the aorta. From this analysis, the motion of the leaflet has been observed with fluttering phenomenon and rebound, and the flow fields of blood have been obtained with recirculation and regurgitation. The results can contribute to the development of design methodology for the curved bileaflet mechanical heart valve.

A kinetic theory analysis is used to study the ultra-thin gas flow field in a gas slider bering, The Boltzmann equation simplified by a collision model is solved by means of a finite difference approximation with the discrete ordinate method. Calculations are made for a flow in a micro-channel between an inclined slider and a moving disk drive platter. The results are compared well with those from the DSMC method. The present method does not suffer from statistical noise which is common in particle based methods and requires much less computational effort.

A mathematical formulation based on two-phase, two-fluid hyperbolic conservation laws is developed to investigate propagation of shock waves in one- and two-dimensions. We used a high resolution upwind scheme called the split-coefficient matrix method. Two extreme cases are computed for validation of the computer code: the states of a pure gas and a pure liquid. Computed results agreed well with the previous experimental and numerical results. It is studied how the shock wave propagation pattern is affected by the void fraction in the two-phase flow. The shock structure in a two-phase flow turned out, in fact, much deviated from the shape well known in the gas only phase.

Evaluation of turbulence models is performed for a better prediction of thermal stratification in an upper plenum of a liquid metal reactor by applying them to the experiment conducted at JNC. The turbulence models tested in the present study are the two-layer model, the $\kappa-\omega$ model, the v2-f model and the low-Reynolds number differential stress-flux model. When the algebraic flux model or differential flux model are used for treating the turbulent heat flux, there exist little differences between turbulence models in predicting the temporal variation of temperature. However, the v2-f model and the low-Reynolds number differential stress-flux model better predict the steep gradient o( temperature at the interface of thermal stratification, and only the v2-f model predicts properly the oscillation of temperature. The LES Is needed for a better prediction of the amplitude and frequency of the temperature fluctuation.