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. The HANARO is composed af inlet plenum, grid plate, core channel with flow tubes and chimney. The reactor core channel is located at about twelve meters (12 m) depth of the reactor pool and cooled by the upward flow that the coolant enters the lower inlet of the plenum, rises up through the grid plate and the core channel and comes out from the outlet of chimney. A fission moly guide tube is extended from the reactor core to the top of the reactor chimney for easily loading a fission moly target under the reactor normal operation. But active coolant through the core can be quickly raised up to the top of the chimney through the guide tube by jet flow. This paper describes an analytical analysis that is the study of the flow behavior through the guide tube under reactor normal operation and unloading the target. As results, it was conformed through the analysis results that the flow rate, reduced to about fourteen kilogram per second (14 kg/s) from the original flow rate of sixteen point three kilogram per second (16.3 kg/s) did not show the guide tube jet.

In this study, the shallow water equations on the sphere is simulated by the proposed method which has high spatial resolution and is based on the CIP(Cubic Interpolated Pseudoparticle) method. The governing equations are approximated on the longitude-latitudinal coordinate system. To avoid the problems resulting from the convergence of the meridians toward high-latitude and singularities on the poles, the semi-Lagrangian and finite volume method are employed. in addition, the CIP method is employed to solve the advection equations and is extended to apply on the equations in the polar coordinate system. The mathematical formulation and numerical results are also described. To verify of the efficiency, accuracy and capability of proposed algorithm, the standard test cases proposed bv Williamson et. al are simulated and the results are compared with other results. As a result, it is found that the present scheme gives a good properties in preserving shapes of solution and settles the pole problems in solving the shallow water equations on the sphere.

Numerical investigations have been performed to examine the effects of the computational grids on the prediction of the flow characteristics inside the turbine cascades. Three kinds of grid system based on H-type grid are applied to the high-turning transonic turbine rotor blades and comparisons with the experimental data and the numerical results of each grid structure have been done. In addition, the grid sensitivity on the estimation of the blade performances has been investigated.

In this paper, we present the theoretical and numerical results of scavenge characteristics in a small subchamber of an HCCI(Homogeneous Charge Compression Ignition) engine. Two theoretical models are proposed in prediction of the scavenge time and the efficiency; one is the non-mixing model in which the input gas(CH4) and the existing gas(air) do not mix at all, and the other is the fully-mixed model in which the two gases are assumed to mix completely before ejection. Focus is also given to the effect on the scavenge performance of the size of the chamber outlet.

To evaluate LBM we performed the simulation of the unsteady two dimensional flow over a square cylinder in a channel in moderate Reynolds number range, $100\~500$ by using LBM and Fractional-Step method. Frist of all we compared LBM solution of Poiseuille flow applied Farout and periodic boundary conditions with the analytical solution to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured 500x100 grids. Prescribed maximum velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback boundary condition was applied to the channel and the cylinder waifs. The flow patterns and vortex shedding strouhal numbers were compared with previous research results. The flow patterns by LBM were in agreement with the flow pattern by fractional step method. Furthermore the strouhal number computed by LBM simulation result was more accurate than that of fractional step method through the comparison of the previous research results.

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

This work discusses the numerical analysis, the design and experimental test of the SMA actuator along with its capabilities and limitations. Convective heating and cooling using water actuate the SMA(Shape memory alloy) element of the actuator. The fuel such as propane, having a high energy density, is used as the energy source for the SMA actuator in order to increase power and energy density of the system, and thus in order to obviate the need for electrical power supplies such as batteries. The system is composed of a pump, valves, bellows, a heater(burner), control unit and a displacement amplification device. The experimental test of the SMA actuator system results in 150 MPa stress(force : 1560 N) with $3\%$ strain and 0.5 Hz. actuation frequency. The actuation frequency is compared with the prediction obtained from numerical analysis. For the designed SMA actuator system, the results of numerical analysis were utilized in determining design parameters and operating conditions.

For the improvement of aerodynamic performance of the turbine blade in a turbopump for the liquid rocket engine, the optimization of turbine profile shape has been studied. The turbine in a turbopump in this study is a partial admission of impulse type, which has twelve nozzles and supersonic inflow. Due to the separated nozzles and supersonic expansion, the flow field becomes complicate and shows oblique shocks and flow separation. To increase the blade power, redesign ol the blade shape using CFD and optimization methods was attempted. The turbine cascade shape was represented by four design parameters. For optimization, a genetic algorithm based upon non-gradient search hue been selected as an optimizer. As a result, the final blade has about 4 percent more blade power than the initial shape.

For study on the unsteady wall interference effect, flows around a forced oscillating airfoil in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The Spalart-Allmaras one-equation model is employed for the turbulence effect. The computed results of the oscillating airfoil having a thin wake showed that the lift curve slope is increased and the magnitude of hysteresis loop is reduced by the interference effects. Since the vortex around the airfoil is generated and convected downstream faster than the free-air condition, the phase of lift, drag and pitching moment coefficients was shifted. The pressure on the test section wall shows harmonic terms having the oscillating frequency contained in the wail effect.

A fluid transient analysis for the propellant flow in a monopropellant propulsion system is conducted by using the method of characteristics(MOC). It reviews algebraic simultaneous equations method and Cramer's rule method utilized to drive the compatible and characteristic equations to understand MOC extensively. The identification of fluid transient phenomena of propulsion system of Koreasat 1 is carried out through parametric studies. The valve response time is one of the dominant parameters governing the fluid transient phenomena. The results show that the shorter closing time induces the greater pressure response amplitude. And it shows that the installation of in-line orifice is effectively to limit the fluid transients in rapid valve response time and at high pressure. But it seems that the effect of orifice weakens at slow valve response time and at low pressures.