A new approach, Multi-Objective Design Exploration (MODE), is presented to address Multidisciplinary Design Optimization (MDO) problems by CFD-CSD coupling. MODE reveals the structure of the design space from the trade-off information and visualizes it as a panorama for Decision Maker. The present form of MODE consists of Kriging Model, Adaptive Range Multi Objective Genetic Algorithms, Analysis of Variance and Self-Organizing Map. The main emphasis of this approach is visual data mining. An MDO system using high fidelity simulation codes, Navier-Stokes solver and NASTRAN, has been developed and applied to a regional-jet wing design. Because the optimization system becomes very computationally expensive, only brief exploration of the design space has been performed. However, data mining result demonstrates that design knowledge can produce a good design even from the brief design exploration.

First, methods of numerical analysis of gas-particle flows is classified into micro, meso and macro scale approaches based on the concept of multi-scale mechanics. Next, the explanation moves on to discrete particle simulation where motion of individual particles is calculated numerically using the Newtonian equations of motion. The author focuses on the cases where particle-to-particle interaction has significant effects on the phenomena. Concerning the particle-to-particle interaction, two cases are considered: the one is collision-dominated flows and the other is the contact-dominated flows. To treat this interaction mathematically, techniques named DEM(Distinct Element Method) or DSMC (Direct Simulation Monte Carlo) have been developed DEM, which has been developed in the field of soil mechanics, is useful for the contact -dominated flows and DSMC method, developed in molecular gas flows, is for the collision-dominated flows. Combining DEM or DSMC with CFD (computer fluid dynamics), the discrete particle simulation becomes a more practical tool for industrial flows because not only the particle-particle interaction but particle-fluid interaction can be handled. As examples of simulations, various results are shown, such as hopper flows, particle segregation phenomena, particle mixing in a rotating drum, dense phase pneumatic conveying, spouted bed, dense phase fluidized bed, fast circulating fluidized bed and so on.

Gas Chromatography (GC) is a wisely technique used for the separation and analysis of liquid and gas sample. Separation of the sample vapors is achieved via their differential migration through a capillary column with an insert carrier gas. The identity and quantity of each vapor in the mixer can be determined from its retention time in the column and a particular property of the gas, such as thermal conductivity, which can be related to the concentration of sample vapor in the carrier gas. Therefore, the flow characteristics in the spiral gas chromatographic column are numerically investigated in this study. Especially, different pressure drop between the front and the rear of GC column with various flow rates is estimated the governing equations are derived from making using of three-dimensional Naver-Stokes equation with incompressible and laminar model due to the nature of low Reynolds number flow. Using a commercial code, FLUENT, the pressure and flow fields in GC column are calculated with various flow rates. The characteristics of thermal cycling which is one of the most important factors affecting the column efficiency and analysis time is also estimated. Furthermore, numerical analyses are also carried out by using commercial code, ANSYS, with various values of power, which is applied to the heating element located at lower GC column.

In a microfluidic chips pressure driven flow or electro-osmotic flow has been usually employed to deliver bio-samples. Flow in the chips is usually slow and the mixing performance is poor. A micro-mixer with a rapid mixing is important for practical applications. In this study a newly designed and electro-osmotic driven micro-mixer is proposed. This design is comprised of a channel and a series of metal electrodes periodically attached on the side surface. In this configuration electro-osmotic flows and the stirring effects are simulated three-dimensionally using a commercial code, CFD-ACE. Focus is given the effect on the electro-osmotic flow characteristics under the local variation of the electric field.

This article has investigated the spatial distribution of the solid particles in polymerization reactors using CFD analysis (FLUENT v. 6.2.1). The suspension of the solids in stirred reactors is affected by a number of parameters including particle diameter, vessel shape, impeller size, impeller speed, and rotating direction of stirrer. The degree of solids suspension in the vessel was quantified with a statistical average value, ${\sigma}^2$. The best stirring conditions were determined based on ${\sigma}^2$, which was found to depend on the vessel bottom shape.

인버터는 전기적으로 DC(직류)를 AC(교류)로 변환하는 역 변환 장치이며, 상용전원(AC 22V/440V)으로부터 공급된 전력을 입력받아 전압과 주파수를 가변시켜 전동기에 공급함으로써 전동기의 속도를 고효율로 이용하게 제어하는 장치이다. 본 논문에서는 550kW급 IGBT, DIODE의 발열에 의한 인버터 내부의 열 및 온도분포를 ICEPAK 상용코드를 통하여 수치적 해석을 수행하였다. 인버터의 발열은 캐패시터의 수명과 소자들의 오작동 등 많은 열적 문제를 가지고 있으며 Heat-sink, Fan, Duck등을 통하여 전도와 대류가 이루어지는 시스템이다. 인버터은 많은 파라미터들에 의해 온도가 결정되기 때문에 실험을 통한 해석은 제한적 일수 밖에 없다. 따라서 수치해석을 통하여 빠른 시간에 효율적인 열 설계를 할 수 있으며, 인버터의 크기는 최종적으로 Heat-sink의 형상에 따라 달라지므로 이를 최적화 하고 소형화 하는 작업이 필요하다. 인버터의 복잡한 내부구조상 하단부(발열원, Heat-sink, Fan등)만 수치해석을 수행하였을 때와 Full Model과는 $15^{\circ}C$ 온도 차이를 보였다. 최종적으로 인버터의 최적 열 설계를 위하여 Frame위치 변경, Heat-Sink 형상변화등 많은 수치해석을 통하여 만족할 만한 결과를 얻었다.

In HANARO, a multi-purpose research reactor of 30 MWth, the emergency water supply system consists essentially of an emergency water storage tank located in the level of about thirteen meter (13 m) above the reactor core, a three inch ('3\%') diameter water injection pipe line including injection valves from the tank to the reactor cooling inlet pipe and a test loop to do periodic system performance test. When the water level of the reactor pool comes down to the extremely low due to a loss of reactor pool water accident the emergency water stored in the tank should be fed to the core by the gravity force and at that time the design flow rate is eleven point four kilogram per second (11.4 kg/s). But it is impossible periodically to measure the injection flow rate under the emergency condition because the normal water level should be maintained during the reactor operation. This paper describes a flow network analysis to simulate the flow rate under the emergency condition. As results, it was confirmed through the analysis results that the calculated flow rate agrees with the design requirement under the emergency condition.

We have investigated the mixing behavior of a tubular heat exchanger reactor using CFD and compared its mixing performance with different type of reactors such as jet mixer and continuous stirred tank reactor (CSTR). The mixedness in each reactor was quantified introducing a statistical average value, the coefficient of variation (CoV), which is a normalized standard deviation of concentration of a component over the whole fluid domain. Through the analysis of the flow pattern and turbulent energy distribution, we suggested a simple but effective way to improve the mixing performance of the tubular heat-exchanger reactor, which include the addition of the internals and/or the increase of the recycle flow rate. It was found that the CoV value of the tubular reactor could be nearly equivalent to that of CSTR by applying those two alternatives suggested here.

This paper presents a multi-scale hybrid simulation for the design of a catalytic multi-tubular reactor with high performance. The multi-tubular reactor consists of shell and a large number of tubes in which various catalytic chemical reactions occur. To consider fluid dynamics in the shell-side and kinetics in the tube-side at the same time, commercial CFD package and process simulation tool are coupled. This hybrid approach allowed us to predict many kinds of meaningful results such as tube center temperature profile, heat transfer coefficients on the tube wall, temperature rise of cooling medium, pressure drop through shell and tube side, concentration profile of each chemical species along the tube, and so on., and to achieve the optimal reactor design.

CFD analysis has been conducted to find the two stage impeller configuration which is the most suitable for a stirred tank with an internal helical cooling coil and a cooling jacket, which is frequently used in chemical industries for highly exothermic reactions ranged from low to medium viscosity. Two typical types of impellers are considered; pitched paddle impellers and Rushton turbine impellers. Interestingly, pitched paddle impellers show a good mixing performance for multi-species, whereas Rushton turbine impellers achieve a good mixing performance for multi-phases. Besides the type of an impeller, the location of an impeller is another important factor to be considered in order to accomplish an effective mixing. The best set of types and locations of two impellers is recommended based on the coefficient of variation(CoV) value and the heat removal capability obtained from CFD results. The former is a measure to quantify the degree of mixing.

Natural convection flows in a cubical, air-filled cavity that has one pair of opposing faces isothermal at different temperatures, Th and Tc, the remaining faces having a linear variation from Tc to Th are numerically simulated by a new solution code(PowerCFD) using unstructured cell-centered method. Solutions are obtained for configurations with a Rayleigh number as high as 105 and three inclination angles ${\theta}$ of the isothermal faces from horizontal: namely ${\theta}=0$, 45 and $90^{\circ}$. Interesting features are presented in detail and comparisons are made with benchmark solutions and experimental results found in the literature. It is found that the code is capable of producing accurately the nature of the laminar convection in a cubical, air-filled cavity with differentially heated walls.

Fluid flow and heat transfer in rectangular duct system are measured and computed by commercial software of Star-CD for comparison between them. Three rectangular systems are investigated in this study. Those are a rectangular duct with 90 degree bended elbow, a rectangular duct with two branchs, and a circular cylinder in a rectangular duct. But heat transfer is studied only for last system. These investigations show us that the numerical solutions predict satisfactorily design factors (K-factor for the elbowed duct, distributions of flow rates into each branch from a duct, and Nusselt number around circular cylinder) even though there are some disagreements in velocity profiles and turbulent kinetic energy.

The present work was undertaken to numerically analyze the defrosting phenomena of windshield glass. In order to analysis the phase change from frost to water on windshield glass by discharging hot air from a defroster nozzle, the flow and the temperature field of the cabin interior, the heat transfer through the windshield glass, and the phase change of frost should be solve simultaneously. In the present work, the flow field was obtained by solving 3-D incompressible Navier-Stokes equations, and the temperature field was computed from the incompressible energy equation. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change of the driven cavity were calculated. The calculation showed a good agreement with other numerical results. Then, the present code was applied to the defrosting problem of a real automobile, and a good agreement with the experimental data was also obtained.

Three-dimensional steady incompressible laminar entry flows in a square duct of $90^{\circ}$ bend are numerically simulated by a new solution code(PowerCFD) using unstructured cell-centered method. Solutions are obtained with three unstructured grid types of hexahedron, prism and hybrid at a Reynolds number, based on the hydraulic diameter and bulk velocity, of 790. Interesting features of the flow are presented in detail. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a $90^{\circ}$ bend of a square duct with fully-developed entry flows. It is found that the code is capable of producing the nature of laminar flow in curved square duct with no grid type dependency.

A conceptual design for the movable roll vane system is done for the roll stability control of KSLV-I. The control effectiveness of the roll vanes is estimated using the numerical simulation. The hinge location is selected to minimize the torque requirement at the maximum dynamic pressure condition, and the maximum torque of 3.0 kN-m is found to be required to actuate the roll vanes for the entire range of operation. An electro-mechanical actuator system which is composed of a DC motor, the speed reducers, the battery package and the controller is designed using the given requirements, the maximum torque of 3.0 kN-m, the maximum deflection angle of 25 deg. and the maximum angular velocity of 30 deg/sec. More detailed design to make more compact and highly efficient system will be done in the future.

Air interlacing serves to protect the yarn against damage, strengthens inter-filament compactness or cohesion, and ensures fabric consistency. The air interlacing nozzle is used to introduce intermittent nips to a filament yarn so as to improve its performance in textile processing. This study investigates the effect of interlacing nozzle geometry on the interlacing process. The geometries of interlacing nozzles with multiple air inlets located across the width of a yarn channels are investigated. The basic interlacing nozzle is the yarn channel, with a perpendicular single air inlet in the middle. The yarn channel shapes are cross sections with semicircular or rectangular shapes. This paper presents three doubled sub air inlets with main air inlet and one of them is slightly inclined doubled sub air inlets with main air inlet. The compressed air coming out from the inlet hits the opposing wall of the yarn channel, divides into two branches, flows trough the top side of yarn channel, joins with the compressed air coming out from the sub air inlet and then creates two free jets at both ends of the yarn channel. The compressed air moves in the shape of two opposing directional vortices. The CFD-FASTRAN was used to perform steady simulations of impinging jet flow inside of the interlace nozzles. The vortical structure and the flow pattern such as pressure contour, particle traces, velocity vector plots inside of interlace nozzle geometry are discussed in this paper.

For the efficient reliability analysis, Bi-direction two-point approximation(BTPA) method is developed which solves shortcomings of conventional two-point approximation(TPA) methods that generate an approximate surface with low accuracy or sometimes do an unstable approximate surface. 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, the approximate reliability analysis approaches on the TPA surface are proposed. Using these FORM and SORM analysis strategies, multi-point aerodynamic-structure interacted shape design optimizations with uncertainty are performed very efficiently.

Evaluation of elliptic blending turbulence model (EBM) together with the two-layer model, shear stress transport (SST) model and elliptic relaxation model (V2-F) 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 algebraic flux model is used for treating the turbulent heat flux. There exist much differences between turbulence models in predicting the temporal variation of temperature. The V2-F model and the EBM better predict the steep gradient of temperature at the interface of thermal stratification, and the V2-F model and EBM predict properly the oscillation of temperature. The two-layer model and SST model fail to predict the temporal oscillation of temperature.

Two coupling methods for the Navier-Stokes equations and a two-equation turbulence model equations are compared. They are the strongly coupled method and the loosely coupled method. The strongly coupled method solves the Navier-Stokes equations and the two-equation turbulence model equations simultaneously, while the loosely coupled method solves the Navier-Stokes equation with the turbulence viscosity fixed and subsequently solves the turbulence model equations with all the flow quantities fixed. In this paper, performances of two coupling methods are compared for two and three-dimensional problems.

Splitting algorithms of the incompressible Navier-Stokes equations using P1P1/P2P1 finite element formulation are newly proposed. P1P1 formulation allocates velocity and pressure at the same nodes, while P2P1 formulation allocates pressure only at the vertex nodes and velocity at both the vertex and mid nodes. For comparison of the elapsed time and accuracy of the two methods, they have been applied to the well-known benchmark problems. The three cases chosen are the two-dimensional steady and unsteady flows around a fixed cylinder, decaying vortex, and impinging slot jet. It is shown that the proposed P2P1 semi-splitting method performs better than the conventional P1P1 splitting method in terms of both accuracy and computation time.

Flow control has been performed using synthetic jet on NACA23012. In order to improve aerodynamic performance, synthetic jet is located near separation paint on airfoil with leading edge droop and plain flap. The flow control using synthetic jet shows that stall characteristics and control surface performance can be improved through resizing separation vortices. Stall is delayed and stall characteristics are improved when synthetic jet is applied from separation region of leading edge droop. Control surface effectiveness is increased and lift is increased when synthetic jet applied at the flap leading edge region. The results show that aerodynamic characteristics can be improved through leading edge droop with synthetic jet at near separation and plain flap with synthetic jet at the flap leading edge. The combination of synthetic jet and simple high lift device is as good as fowler flap system.

Conventional high order interpolation schemes are limitative in several aspects mainly because they need data of neighboring cells at the reconstruction step. However, discontinuous Galerkin method and spectral volume method, two high order flux schemes which will be analyzed and compared in this paper, have an important benefit that they are not necessary to determine the flow gradients from data of neighboring cells or elements. These two schemes construct polynomial of variables within a cell so that even near wall or discontinuity, the high order does not deteriorate.

An unstructured overset mesh method has been developed for the simulation of unsteady viscous flow fields around multiple bodies in relative motion. For this purpose, a robust and fast search technique is proposed for both triangle and high-aspect quadrilateral cell elements. The interpolation boundary is defined for data communication between grid systems and an interpolation method is suggested for viscous and inviscid cell elements. This method has been applied to calculate the flow fields around 2-D airfoil including relative motion. Validation were made by comparing the predicted results with those of experiments or other researcher's numerical results. It was demonstrated that the present method is efficient and robust for the prediction of unsteady time-accurate flow fields involving multiple bodies in relative motion.

The fluid-structure interaction analysis such as a static aeroelastic analysis requires the result of each analysis as an input to other analysis. Usually the grids for the fluid analysis and the structural analysis are different, so the results should be transformed properly for each other. The Infinite Plate Spline(IPS) and the Thin Plate Spline(TPS) are used in interpolating the displacement and the pressure. In this study, such interpolation methods are compared with kriging which provides a precise response surface. The static aeroelastic analysis is performed for the supersonic flow field with shock waves and the pressure field is interpolated by the TPS and kriging. The TPS shows tendency to weaken the shock stength, whereas kriging preserves the shock strength.

Generally, Compressible Navier-Stokes codes are used to solve high mach number flows. But, Most of high mach number flows embrace low mach number flows. This phenomenon results in low convergence rate and non-physical solution in CFD analysis. So Many researchers developed preconditioning technique to solve these problems. This Study presents how to modify previous compressible N-S computer code with little changes of structure into preconditioned compressible N-S code applying Roe's Approximate Riemann Solver. And this study show developed preconditioning code is very well operated at all mach number flows.

Three-dimensional structures of unsteady detonation wave propagating through a square-shaped tube is studied using computational method and parallel processing. Inviscid fluid dynamics equations coupled with variable-${\gamma}$ formulation and simplified one-step Arrhenius chemical reaction model were analysed by a MUSCL-type TVD scheme and four stage Runge-Kutta time integration. Results in three dimension show the two unsteady detonation wave propagating mode, the Rectangular and diagonal mode of detonation wave instabilities. Two different modes of instability showed the same cell length but different cell width and the geometric similarities in smoked-foil record.

Parallel implementation and performance assessment of the grid assembly in a structured chimera grid approach is studied. The grid assembly process, involving hole cutting and searching donor, is parallelized on the PC cluster. A message passing programming model based on the MPI library is implemented using the single program multiple data(SPMD) paradigm. The coarse-grained communication is optimized with the minimized memory allocation because that the parallel grid assembly can access the decomposed geometry data in other processors by only message passing in the distributed memory system such as a PC cluster. The grid assembly workload is based on the static load balancing tied to flow solver. A goal of this work is a development of parallelized grid assembly that is suited for handling multiple moving body problems with large grid size.

A research to evaluate efficiency of design optimization was performed for aerodynamic design optimization problem in distributed computing environment. The aerodynamic analyses which take most of computational work during design optimization were divided into several jobs and allocated to associated PC clients through network. This is not a parallel process based on domain decomposition rather than a simultaneous distributed-analyses process using network-distributed computers. GBOM(gradient-based optimization method), SAO(Sequential Approximate Optimization) and RSM(Response Surface Method) were implemented to perform design optimization of transonic airfoil and to evaluate their efficiencies. One dimensional minimization followed by direction search involved in the GBOM was found an obstacle against improving efficiency of the design process in distributed computing environment. The SAO was found quite suitable for the distributed computing environment even it has a handicap of local search. The RSM is apparently the fittest for distributed computing environment, but additional trial and error works needed to enhance the reliability of the approximation model are annoying and time-consuming so that they often impair the automatic capability of design optimization and also deteriorate efficiency from the practical point of view.

A high resolution scheme for solving gas-liquid two-phase flows with cavitation is described. This scheme uses the curvilinear coordinate grid and solves the density based momentum equations for mixture of gas-liquid medium with a preconditioning method to treat both compressible and incompressible flow characteristics. The present preconditioned method is based on the Runge-Kutta explicit finite-difference scheme, and is improved by using the diagonalization, the flux difference splitting and the MUSCL-TVD schemes to save computational effort and to increase stability and resolvability, especially at gas-liquid contact surfaces. A homogeneous equilibrium cavitation model is used to treat the gas-liquid two-phase medium in cavitating flow as a locally homogeneous pseudo-single-phase medium. Therefore, it is easy to solve cavitating flow, including wave propagation, large density changes and incompressible flow characteristic at low Mach number. Some numerical results obtained by the present scheme are shown.

We have developed a numerical scheme to reproduce the unsteady flows with cavitation by the finite-difference method. The evolution of cavitation is represented by the source/sink of vapor phase in the incompressible liquid flow. The pressure-velocity coupling is based on the fractional-step method for incompressible fluid flows, in which the compressibility is taken into account through the low Mach number assumption. We applied our method for the cavitating flows in a two-dimensional cascade, which approximates the portion near the tip of inducer in liquid-fuel engine. Particular attention was focused on the influence of turbulence model in this report. Using an eddy viscosity model, although it was not an optimized one for our purpose, the agreement with the experimental observation was improved.

We developed an upwind numerical formulation based on the eigenvalues of the approximate Jacobian matrix in order to solve the hyperbolic conservation laws governing the two-fluid two-phase flow models. We obtained eight analytic eigenvalues in the two dimensions that can be used for estimate of the wave speeds essential in constructing an upwind numerical method. Two-dimensional underwater cavitation in a flow past structural shapes or by underwater explosion can be solved using this method. We present quantitative prediction of cavitation for the water tunnel wall and airfoils that has both experimental data as well as numerical results by other numerical methods and models.

For the analysis of compressible multi-phase and real gas flows, characteristic form of Roe's Riemann solver was derived using real gas equation of state. It was extended to multi component reactive system considering variable specific heat. From this study, it is known that some correction should be made for the use of existing numerical algorithm. 1) Sonic speed and characteristic variable should be corrected with real gas effect. 2) Roe's average was applicable only with the assumption of constant properties. 3) Artificial damping term and characteristic variables should be corrected but their influences may not be significant.

Wave breaking phenomenon near the fore body of a ship is numerically simulated. The ship advance with uniform velocity in calm water. For the simulation, incompressible Navier-Stokes equations and continuity equation are adopted as governing equations. The simulation is carried out in staggered variable mesh system with finite difference method. Marker and Cell(MAC) method and Marker-Density method are employed to track the free surface. Body boundary conditions are satisfied with the adoption of porosity method and no-slip condition on the hull surface. The ship model has a wedge type fore-body, and the computational domain is an appropriate region around the fore-body. The computation results are compared with some experimental results. Also the difference of the free surface tracking methods are discussed.

In order to realistically represent the complex turbulence-chemistry interaction at the partially premixed turbulent lifted flames encountered in the gas turbine combustors, the combined conserved-scalar/level-set flamelet approach has been adopted. The parallel unstructured-grid finite-volume method has been developed to maintain the geometric flexibility and computational efficiency for the solution of the physically and geometrically complex flows. Special emphasis is given to the swirl effects on the combustion characteristics of the lean-premixed gas turbine combustor. Numerical results suggest that the present approach is capable of realistically simulating the combustion characteristics for the lean-premixed gas turbine engines and the lifted turbulent jet flame with a vitiated coflow.

In designing a ventilation system of a road tunnel, a possibility of using the system as a smoke control system in case of a tunnel fire has to be considered. In the present study, a numerical simulation on ventilation system is performed considering jet fan operations and moving traffic. A fire-mode operation by reversing some fan operations in case of a tunnel fire is also simulated. The results show that ventilation operation can control the pollutants effectively, and fire-mode operation can control smoke and temperature effectively to prevent a disaster.

In this study the water including sediment is assumed that the density is different from fresh water. And the phenomenon inhaled by low pressure around the pipe is numerically simulated in two dimension. The simulation is done using finite difference method in rectangular staggered mesh system and Navier-Stokes equations and continuty equation are employed as governing equations. The method of Irregular leg lengths and stars are adopted to satisfy boundary condition of body boundary. Marker-Density method is used to calculate the density of mixed flow. Also SGS turbulence model is applied to consider vortex smaller then grid at high Renolds number. This study is to analyze inhalation phenomenon of mixed flow with sediment and to verify the numerical method for mixed flow. To verify the numerical results are compared with experimental results

Recently, the interest of the engine capacity and environment of the atmosphere is increasing, so the researches for the engine capacity have been conducted for a long time. But the internal environment of an automotive engine is very severe. A piston is exposed to combustion gas of over $2000^{\circ}C$ and strong friction is occurred by high speed motion in the cylinder. The fraction between piston and wall of the cylinder causes the increase of temperature in the engine. The temperature of the engine has an effect on the engine capacity. If the temperature is high, the capacity of the engine is low. So we have to maintain the optimum temperature. To maintain the optimum temperature, the enough flow rate of the engine oil is needed. The oil jet is used to control the flow rate of the engine oil and supply the engine oil to the piston and cylinder. The purpose of this study is to check the mass flow rate of the engine oil and the characteristics of internal flow of the oil jet. Flow pattern of the engine oil is very important because it concludes the loss in the oil jet. This study is the previous research about the oil jet and we will consider the movement of the ball check valve to get more accuracy result.

A PCV valve is a part to control the flow rate of Blowby gas in a PCV system. A PCV system re-burns Blowby gas with fuel in a combustion chamber. Some gas enters to a crankcase room through the gap between piston ring and engine cylinder wall. This gas si called 'Blowby gas'. This gas causes many problems. In environmental view, Blowby gas includes about $25\~35\%$ hydrocarbon{HC) of total generated HC in an automobile. Hydrocarbon is a very harmful pollutant element in our life. In mechanical view, Blowby gas has some reaction with lubricant oil of crankcase room. Then, this causes lubricant oil contamination, crankcase corrosion and a decrease fo engine efficiency. Consequently, Blowby gas must be eliminated from a crankcase room. In this study, we simulated internal flow characteristics in a PCV valve according to spool dynamic behavior using local remeshing method And, we programmed our sub routine to simulate a spool dynamic motion. As results, spool dynamic behavior is periodically oscillated by the relationship between fluid force and elastic force of spring. And its magnitude is linearly increased by the differential pressure between inlet and outlet. Also, as spool is largely moved, flow area is suddenly decreased at orifice. For this reason, flow velocity is rapidly decreased by viscous effect.

Following the recent trend in the automotive manufacturing technologies, the product design subject to the die casting becomes more and more complex. The requirement of the die design becomes more demanding than ever before. In some cases the product's shape can have multiple slender manifolds. In such cases, design of the inlet and outlet parts of the die is very important in the whole manufacturing process. The main issues required for the qualified products are to attain gentle and uniform flow of the molten liquid within the passages of the die. To satisfy such issues, the inlet cylinder ('bed cylinder' in this paper) must be as large as possible and simultaneously the outlet opening at the end of each passage must be as small as possible. However these in turn obviously bring additional manufacturing costs caused by re-melting of the bed cylinder and increased power due to the small outlet-openings. The purpose of this paper is to develop effective simulation methods of calculation for fluid flows in multiple columns, which mimic the actual complex design, and to get some useful information which can give some contributions to the die-casting industry. We have used a commercial code CFX in the numerical simulation. The primary parameter involved is the size of the air-vent. We will show how the very small opening of the outlet, i.e. the air-vent, can be treated with the aid of the porous model provided in the code. To check the validity of the numerical results we have also conducted a simple experiment by using water.

The unsteady supersonic flow over two- and three-Dimensional cavities has been analyzed by the integration of unsteady Reynolds-Averaged Navier-Stokes(RANS) with the k - w turbulence model. The unsteady flow is characterized by the periodicity due to the mutual relation between the shear layer and the internal flow in cavities. Numerical method is upwind TVD scheme based on the flux vector split with the Van Leer limiters, and time accuracy is used explicit 4th stage Runge-Kutta scheme. Cavity flows are Comparison of two- and three-dimensional. The cavity has a L/D ratio of 3 for two-dimensional case. and same L/D and W/D ratio is 1 for three-dimensional case. The Mach and Reynolds numbers are held constant at 1.5 and 450000 respectively. For the three-dimensional case, the flow field is observed to oscillate in the 'shear layer mode' with a feedback mechanism that follow Rossiter's formula. On the other hand, the self-sustained oscillating flow transitions to a 'wake mode' for the two-dimensional simulation, with more violent fluctuations inside the cavity.

Fully three-dimensional Large Eddy Simulation calculations of the flow past 2D cavity are conducted to study the purging of neutrally buoyant or dense miscible contaminants introduced instantaneously inside the cavity. The length to depth ratio(L/D) is 2 and Reynolds number based on the depth is 3,360. Fully developed turbulent inflow are fed at the inlet from precursor simulation of channel flow. Mean flow pattern and unsteady features are investigated based on the experimental data of Pereira and Sousa. From the study of mass exchange processes, it is found that the mechanism of removal of the contaminant is very different between the non-buoyant and buoyant cases. In the buoyant case, internal wave motion which interacts with a strong cavity vortex is dominant in the ejection mechanism of the contaminants.

SCRamjet is the key technology for hypersonic flight over mach number 6. It is characterized by very short residence time in combustor because its internal flow is supersonic. In this short time, the whole process of combustion must be done. Especially numerical study of combustor is important because air-fuel mixing rate influences the performance of combustor. Various methods of air-fuel mixing enhancement are proposed. Among these, cavity injection method is selected to study in this paper. The numerical study is conducted with the variation of the cavity length at the fixed height of unit and jet injection on the downstream of cavity.

The turbulence noise generated from blunt trailing-edge is numerically predicted by using the hydrodynamic/acoustic splitting method at the Reynolds number based on thickness of flat plate, $Re_h=1000$, and the freestream Mach number $M_o=0.2$. The turbulent flow field is simulated by incompressible large-eddy simulation and the acoustic field is predicted efficiently with the linearized perturbed compressible equations (LPCE) recently proposed by the authors. The turbulent flow characteristics are validated with the results of the previous experimental study and direct numerical simulation. The acoustic properties predicted from LPCE are compared with the solutions of analytical formulations.

An analysis is made of flow and rocket motion during a supersonic separation stage of air-launching rocket from the mother plane. Three-dimensional Euler and Navier-Stokes equations are numerically solved to analyze the steady/unsteady flow field around the rocket which is being separated from two cases of mother plane configuration: one is an idealized ogive-cylinder body and the other is a real F-4E Phantom. The simulation results clearly demonstrate the effect of shock-expansion wave interaction between the rocket and the mother plane. As a result, a design-guideline of supersonic air-launching rocket for the safe separation is proposed.

The unsteady flow analysis of staging system is conducted. This study focuses on comparing the results of two different governing equations between Euler equations and Navier-Stokes equations. The Chimera grid scheme is applied to moving simulations for unsteady flow analysis with dynamic simulation. As a result, it is certified that inviscid simulation have capabilities enough to analyze the present staging problem.

The on-board system for the air supply to the payload fairing(PLF) of a launch vehicle using both high and low pressure air was designed. The design concept was obtained from the CFD analysis of a Russian interstage air supply system, and a collector was adopted to expand the high pressure air. To verify that the on-board system would work as designed, a simplified axisymmetric computational model was made and a CFD analysis was also performed. It was found that the flow ejected from the hole of the collector expands to the Mach number of 4 and is soon retarded due to the action of viscosity. It was also found that a small gap between the low pressure duct and equipment bay wall can cause large velocity in PLF over the velocity requirement and no gap should be allowed in the design.

The vortex convection and induced flow field behind the KARI 3m x 4m LSWT gust generator was computed by using Computational Fluid Dynamics. For the accurate simulation of vortex convection, inviscid, laminar, Spalart-Allmars k-e and k-w turbulence models were tested with the NAL gust generator configuration and Spalart-Allmaras turbulence model was selected for the prediction of induced flow field behind the KARI LSWT gust generator. The wind tunnel test was also carried out at KARI LSWT and the results were compared with CFD prediction.

Viscous solutions of supersonic jet impinging on a flat wall in a confined plenum are simulated using three-dimensional Navier-Stokes solver. A confined plenum was designed for simulating the missile launch and analyzing the behavior of the exhaust plume, which were accompanied by complex flow interactions with shock and boundary layer. Concerns of this paper are to show accurate simulation of internal flow in confined plenum and to demonstrate the jet flow structure when the jet interacts with a small opening on the side. Objectives of this numerical simulation are to understand the effect of changing the plume exit area of the plenum. Pressure and temperature rise at certain position in the plenum are traced and compared with test data.

The flowfields generated by gaseous slot injection into a supersonic flow at a Mach number of 3.75 and a Reynolds number of $2.07{\times}10^7$ are simulated numerically. Fine-scale turbulence effects are represented by a two-equation(k-w SST model) closure model which includes $y^+$ effects on the turbulence model. Grid convergence index(GCI) is also considered to provide a measure of uncertainty of the grid convergence. Comparison is made with experimental data and other turbulence model in term of surface static pressure distributions, the length of the upstream separation region, and the penetration height. Results indicate that the k-w SST model correctly predicts mean surface pressure distribution and upstream separation length. However, it is also observed that the numerical simulation over predicts the pressure spike and penetration height compared with experimental data. All these results are taken within $1\%$ error band of grid convergence.

We calculate the coordinates of an axisymmetric nozzle with a central body. This nozzle ensures a transonic flow with a plane sound surface, which is orthogonal to the symmetry axis and has a wall kink at the sonic point, The Chaplygin transformation in the subsonic part of the flow leads the Dirichlet problem for a system of nonlinear equations. The definition domain of the solution in the velocity-hodograph plane is taken as a rectangle. This enables one to obtain the nozzle with a monotonic distribution of velocity along its subsonic part. In the nonlinear differential equation, the linear Chaplygin operator for plane flows is separated, which allows the iterative calculation of the solution. The supersonic part of the nozzle is calculated under the assumption that the flow at the nozzle exit is uniform and parallel to the symmetry axis; i.e., the supersonic jet outflows to the submerged space with the same pressure. The calculation is performed by the characteristic method. The exact solution of Tricomi equation for near-sonic flows with the straight sonic line is used to 'move away' the sound plane. The velocity distribution alone the supersonic part of the nozzle is also monotonic, which ensures the absence of the boundary-layer separation and, therefore, the adequacy of the ideal-gas model. calculations show that the flow in the supersonic part of the nozzle is continuous (compression shocks are absent)

The present study focuses on the flow analysis of a turbopump inducer by performing both numerical and experimental methods. The head rise, efficiency and detailed flow fields such as outlet flow angles, pressure and velocity vectors are measured and compared with the computational data. Generally a good agreement is obtained between numerical and experimental results. However, some discrepancies are observed due to complex flow structures inside the inducer. Future calculations with an advanced turbulence model and a dense computational grid needs to be performed to obtain accurate numerical solution for the detailed flow fields.

The velocity and pressure fields of a ship's propulsion mechanism of the Weis-Fogh type, in which a airfoil moves reciprocally in a channel, are studied in this paper using the advanced vortex method. The airfoil and the channel are approximated by a finite number of source and vortex panels, and the free vortices are introduced from the body surfaces. The viscous diffusion of fluid is represented using the core-spreading model to the discrete vortices. The velocity is calculated on the basis of the generalized Biot-Savart low and the pressure field is calculated from integrating the equation given by the instantaneous velocity and vorticity fields. Two-dimensional unsteady viscose flows of this propulsion mechanism are numerically clarified, and the calculated results agree well with the experimental ones.

It has been the most progressive interruption technique to use the ablation gas from the surface of PTFE nozzle driven by arc plasma during switching process in $SF_6$ gas circuit breakers. This advanced interruption technique can reduce the required mechanical energy to compress and blow the gas for extinguishing the arc plasma between the electrodes due to using the ablation effect instead. In order to consider the phenomena during calculation of switching process, it is required to confirm the principles of ablation from PTFE nozzle as well as of arc plasma during switching process. In this study, we have calculated the switching process considered the ablation of PTFE nozzle driven by arc plasma using multidisciplinary simulation technique and compared the results with the data without the ablation effect. More $50\%$ difference of pressure rise inside expansion chamber has been found from the results and it should be indispensable for this type of computational work to consider and include the ablation effect of PTFE nozzle. Further study on turbulence and radiation will be followed.