This article presents the unique characteristics of cavitation, those are very high pressure and very high temperature even in a very short time. Such the high-pressure causes the destruction of material, which sometimes brings a severe problem to fluid machinery. However, if we look the cavitation from the different direction, such the high pressure can be used to various applications. The author presents two examples of the utilization of cavitation recently done at the Toyo University. Those are 'dispersion of spilled oil' and 'killing planktons'.

The issues associated with the application of CFD for ship design are addressed. It is quite certain that the CFD tools are very useful in evaluating hull forms a prior to traditional towing tank tests. However, the time-consuming pre-processing is an obstacle in the daily application of CFD tools to improve hull forms. The accuracy of computational modeling without sacrificing the usability of CFD system is also to be assessed. The wave generation is still predicted by using potential panel methods, while velocity profiles entering into propeller plane is solved using turbulent flow solvers. The choice of turbulence model is a key to predict nominal wake distribution within acceptable accuracy. The experimental data for CFD validation are invaluable to improve physical and numerical modeling. Other applications of CFD for ship design than hull form improvement are also given. It is certain that CFD can be a cost-effective tool for the design of new and better ships.

Human circulatory system between heart and tissue is not directly connected in normal condition but mandatory to go through the capillary system in order to fulfill its physiologic aim to deliver oxygen and nutrients, etc. to the tissue and retrieve used blood together with waste products from the tissue properly. When abnormal connection between arterial and venous system (AV fistula), these two circulatory systems respond differently to the hemodynamic impact of this abnormal connection between high pressure (artery) and low pressure (vein) system. Depending upon the location and/or degree (e.g. size and flow) of fistulous condition, each circulatory system exerts different compensatory hemodynamic response to this newly developed abnormal inter-relationship between two systems in order to minimize its hemodynamic impact to own system of different hemodynamic characteristics. Pump action of the heart can assist the failing arterial system directly to maintain arterial circulation against newly established low peripheral resistance by the AV fistula during the compensation period, while it affects venous system in negative way with increased venous loading. However, the negative impact of increased heart action to the venous system is partly compensated by the lymphatic system which is the third circulatory system to assist venous system independently with different hemodynamics. The lymphatic system with own unique Iymphodynamics based on peristaltic circulation from low resistance to high resistance condition, also increases its circulation to assist the compensation of overloaded venous system. Once these compensation mechanisms should fail to fight to newly established hemodynamic condition due to this abnormal AV connection, each system start to show different physiologic ${\underline{de}compensation}$ including heart and lymphatic system. The vicious cycle of decompensation between arterial and vein, two circulatory system affecting each other by mutually negative way steadily progresses to show series of hemodynamic change throughout entire circulation system altogether including heart. Clinical outcome of AV fistula from the compensated status to decompensated status is closely affected by various biological and mechanical factors to make the hemodynmic status more complicated. Proper understanding of these crucial biomechanical factors iii particular on hemodyanmic point of view is mandatory for the advanced assessment of biomechanical impact of AV fistula, since this new advanced concept of AY fistula based on blomechanical information will be able to improve clinical control of the complicated AV fistula, either congenital or acquired.

In this paper, we present two successful results from active controls of flows over a circular cylinder and a sphere for drag reduction. The Reynolds number range considered for the flow over a circular cylinder is 40-3900 based on the free-stream velocity and cylinder diameter, whereas for the flow over a sphere it is $10^{5}$ based on the free-stream velocity and sphere diameter. The successful active control methods are a distributed (spatially periodic) forcing and a high-frequency (time periodic) forcing. With these control methods, the mean drag and lift fluctuations decrease and vortical structures are significantly modified. For example, the time-periodic forcing at a high frequency (larger than 20 times the vortex shedding frequency) produces $50{\%}$ drag reduction for the flow over a sphere at $Re=10^{5}$. The distributed forcing applied to the flow over a circular cylinder results in a significant drag reduction at all the Reynolds numbers investigated.

Non-uniform and unsteady inflow into a Horizontal Axis Wind Turbine (HAWT) brings about an asymmetric flow field on the rotor plane and an unsteady aerodynamic load on the blades. In the present paper effects of yawed inflow and wind shear are analyzed by an inviscid aerodynamic model based on the asymptotic acceleration potential method. In the analysis the rotor blades are represented by spanwise and chordwise pressure distribution composed of analytical first-order asymptotic solutions for the Laplace equation. As the actual wind field experienced by a HAWT is turbulent, the effects of the turbulence are also examined using the Veers' model.

An evaluation of one algebraic and two one-equation eddy viscosity-transport turbulence closure models as implemented to the CFDS(Characteristic Flux Difference Splitting) scheme is presented for the efficient computation of the turbulent flow. Comparisons of Baldwin-Lomax model as algebraic turbulence model and Baldwin-Barth and Spalart-Allmaras model as one-equation turbulence model are presented for three test cases for 3-dimensional flow. The numerical result of the CFDS schem is examined through comparison with the experimental data.

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. In this study, 3-dimensional unsteady incompressible Wavier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade' scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

An efficient implicit multigrid method is presented for the Navier-Stokes and $k-{\omega}$ turbulence equations. Freezing and limiting strategies are applied to improve the robustness and convergence of the multigrid method. The eddy viscosity and strongly nonlinear production terms of turbulence are frozen in the coarser grids by passing down the values without update of them. The turbulence equations together with the Navier-Stokes equations, however, are consecutively solved on the coarser grids in a loosely coupled fashion. A simple limit for k is also introduced to circumvent slow-down of convergence. Numerical results for the unseparated and separated transonic airfoil flows show that all computations converge well without any robustness problem and the computing time is reduced to a factor of about 3 by the present multigrid method.

The convergence acceleration methods for the compressible Wavier-Stokes equations are studied ,which are multigrid method and implicit preconditioned multistage time stepping method. In this paper, the performance of implicit preconditioning methods are studied for the full-coarsening multigrid methods on the high Reynolds number compressible flow computations. The effect of numerical flux on the convergence are investigated for the inviscid and viscous calculations.

The nuclear transmutation technology to Incinerate the long lived radioactive nuclides and produce energy during the incineration process is believed to be one or the best solutions. HYPER(${\underline{HY}}brid {\underline{P}}ower {\underline{E}}xtraction {\underline{R}}$eactor)is the accelerator driven transmutation system which is being developed by KAERI(Korea Atomic Energy Research Institute). Lead-bismuth(Pb-Bi) is adopted as a coolant and spallation target material. In this paper, we performed the thermal-hydraulic analysis of HYPER target using the commercial code FLUENT, and also calculated thermal and mechanical stress of the beam window using the commercial code ANSYS. It is found that there is an optimum value for the window diameter and the maximum allowable beam current can be increased to 17.3 mA for the inner diameter of windows, 40 cm. Finally, the other shapes such as uniform or scanned beam were considered. The results of FLUENT calculations show that the uniform type is preferable to the other shapes of the beam in terms of the window and target cooling and the maximum window temperature is lower than that of the parabolic beam by $58 ^{\circ}C$ for the beam current, 13 mA.

The HANARO, multi-purpose research reactor, 30 MWth open-tank-in-pool type, is under 24 MWth of power operation since it reached to the initial critical in February, 1995. Many useful experiments should be safely performed to activate the utilization of the HANARO, but there is a radioactive risk of using the HANARO. To reduce the risk, a test facility, which is not reacted by nuclear fuel, is being developed to simulate similar flow characteristics with the HANARO. This paper describes the computational flow analysis to determine each shape of simulating fuels for simulating the flow similarities of 36 elements hexagonal fuels assembly and 18 elements circulating fuels assembly loaded in HANARO. The shares of orifices were determined by the trial and error method and the structural integrities of them were verified by the finite element method assuming that the flow rate and pressure differences of reactor core are constant. The analysis results will be verified with the results of the flow test to be performed after the installation of this test facility.

The numerical simulation of a 271-rod fuel assembly of nuclear Liquid-Metal Fast Reactor (LMFR) with an infernal blockage has been carried out. Internal blockage within a subassembly is addressed in the safety assessment because it potentially has very serious consequences for the reactor as a whole. Three dimensional calculations were performed using the SABRE4 computer code for the range of blockage positions and sizes to investigate the seriousness and detectability of the internal blockage. The magnitude and location of the peak temperatures together with the temperature distribution at the subassembly exit were calculated in order to look at the potential for damage within the subassembly, and the possibility of blockage detection. The analysis result shows that the 6-subchannel blockage causes large temperature rise within a assembly with practically no change in mixed mean temperature at the assembly exit.

When a cold HPSI (High Pressure Safety Injection) fluid associated with an design basis accident, such as LOCA (Loss of Coolant Accident), enters the cold legs of a stagnated primary coolant loop, thermal stratification phenomena will arise due to incomplete mixing. If the stratified flow enters a reactor pressure vessel downcomer, severe thermal stresses are created in a radiation embrittled vessel wall by local overcooling. Previous thermal-mixing analyses have assumed that the thermal stratification phenomena generated in stagnated loop of a partially stagnated coolant loop are neutralized in the vessel downcomer by strong flow from unstagnated loop. On the basis of these reasons, this paper presents the thermal-mixing analysis results in order to identify the fact that the cold plume generated in the vessel downcomer due to the thermal stratification phenomena of the stagnated loop is affected by the strong flow of the unstagnated loop.

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

In the present study, solution algorithms for the computation of unsteady flows on an unstructured mesh are presented. Dual time stepping is incorporated to achieve the 2-nd order temporal accuracy while reducing the linearization and the factorization errors associated with a linear solver. Hence, any time step can be used by only considering physical phenomena. Gauss-Seidel scheme is used to solve linear system of equations. Rigid motion and spring analogy method fur moving mesh are all considered and compared. Special treatments of spring analogy for high aspect ratio cells are presented. Finally, numerical results for oscillating wing are compared with experimental data.

The present study introduces some useful methods of implementing the Internet numerical analysis program with existing numerical codes for utilizing the Internet environment. The Internet gives developers good environment for development and release. Several methods were suggested, and some of them were implemented with an existing numerical code named SOLA-VOF, a computational fluid dynamics program to solve two-dimensional transient flow problems with free surface. This was reconstructed with Java technologies and compared with the original one. Java technologies have been applied to development of Internet applications for a long time. The objective of this work is to contrive methods of implementing Internet numerical analysis program with existing numerical codes and confirm the possibility of them. Methods using the applet-servlet communication were suggested and implemented. In addition, the Java web services with XML was introduced, which makes possible the cooperation of components. Although the concept has been suggested and developed for business applications, it can also be used for engineering softwares. Therefore, this study will be a preparation for numerical analysis to participate in engineering web services.

In the present study a pre/post-processor program has been developed to be used with a general CFD code. This program is capable of performing the basic functions of the pre/post-processing, which include mesh generation and post processing plots. Also through perspective projection, this program can be used to check the quality of generated mesh by moving around inside the mesh. The smoke visualization can be also performed with the present program to visualize the smoke behavior in the case of fire simulation. The examples of the program execution are given in paper.

Fluid flow in a rectangular duct system are measured by W laser doppler velocity meter, and also computed by commercial software of STAR-CD for comparison between then First, for a rectangular duct with 90 degree metered elbow, the fluid flow with Reynolds numbs's of 1,508 is predicted by assumption of both laminar and turbulent models. But, even though the Reynolds number is less than 2,300-3,000, the computation by turbulent model is close to the experimental data. Moeover, the computation by turbulent model for Reynolds number of 11,751 also predicts the experimental data satisfactorily. Second, for a rectangular duct with two branch ducts, the ratios between flow rates in the two branches are invariant to Reynolds number according to both of numerical and experimental results.

In this paper, flow on the rear region of a butterfly valve was analysed by using numerical and experimental methods. The butterfly-valve disk angle is changed as 0-60 degree and the uniform flow velocity was fixed In this experiment. It was shown that the numerical results are similar to the experimental results. General discussions are given to the flow-pattern change upon the disk angle of the valve.

Three-dimensional pipe flows with elbows and tees for few different pipe fittings are calculated to estimate the effect of swirling flow on measuring accuracy of orifice flow meter. It is evaluated how the pressure difference across the orifice is dependent on the length of upstream straight pipe in a branch and how swirl intensity, swirl angel and axial velocity distribution affect the measuring error of orifice flowmeter. From the results, it is found that, regardless of flow rate specified in this calculation, the effect of the straight pipe length can be neglected for the lengths larger than thirty diameters although there still remain significant swirl at the orifice

Numerical simulations are conducted for laminar flow past a sphere rotating In the transverse direction, in order to investigate the effect of the rotation on the characteristics of flow over a sphere. The Reynolds numbers considered are Re=100, 250 and 300 based on the free-stream velocity and the sphere diameter, and the rotational speeds are in the range of $0{\leq}{\omega}{\leq}1$, where ${\omega}^{\ast}$ is the maximum velocity on the sphere surface normalized by the free-stream velocity. At ${\omega}^{\ast}=0$ (without rotation), the flow past the sphere experiences steady axisymmeoy, steady planar-symmetry and unsteady planar-symmetry, respectively, at Re=100, 250 and 300. However, with rotation, the flow becomes planar-symmetric for all the cases investigated and the symmetry plane is orthogonal to the axis of the rotation. The flow is also steady or unsteady depending on both the Reynolds number and the rotational speed, and the vortical structures behind the sphere are significantly modified by the rotation. For example, at Re=300, hairpin vortices completely disappear in the wake at ${\omega}^{\ast}=0.4\;and\;0.6$, and at ${\omega}^{\ast}=1$ vortical structures of a high frequency are newly generated due to the shear layer instability. It is also shown that with increasing rotational speed, the time-averaged drag and lift coefficients increase monotonically.

Computational study has been undertaken to investigate the aerodynamic influence of side jet on a supersonic missile and to find a similarity condition between the flight condition and the wind tunnel testing. Tasks were performed to validate the existing Raytheon test body with side jet, to simulate the flow inside the supersonic wind tunnel, and to compare the flow fields between the missile in free flight and that in the wind tunnel. Then sub-scale model of body-tail configuration was analyzed to estimate the influence of the side jet on the missile components. It Is found that the influence of side Jet is not as significant on the tail region as on the body surface and a simple algebraic formula for aerodynamic coefficients accounting for the side jet as a point force may be cautiously utilized in setting up control logic.

It is important to develop new effective technologies for increasing the interruption capacity and reducing the size of a GCB (Gas Circuit Breaker). It is not easy to test the real GCB model in practice as in theory. Therefore, a simulation tool based on a CFD (Computational Fluid Dynamics) algorithm has been developed to facilitate an optimization of the interrupter. But the choice of grid is not at all trivial in the complicated geometries like a GCB. In this paper, we have applied a CFD-CAD integration using Cartesian cut-cell method, which is one of the grid generation techniques for dealing with complex and multi-component geometries.

Three-dimensional flow analyses for two different ratios of radius to height of sedimentation bed are implemented to evaluate the effect of blockage ratio of center feed wall and angle of distributor on sedimentation efficiency, and to find the optimal value of those parameters. Sedimentation efficiencies for three different shapes are compared with and without rotation speed. And then, five different combinations of blockage ratio of center feed wall and angle of distributor are compared It reveals that the effect of blockage ratio of center fled wall and angle of distributor is considerable to sedimentation efficiency while rotation effect can be neglected and $0.55 and 33^{\circ}$for blockage ratio of center food wall and angle of distributor, respectively, ive the best sedimentation efficiency.

The present study describes a numerical analysis for simulation of the sloshing of flows with free-surface which contained in a rectangular tank The SOLA-VOF (Volume of fluid) method uses a fixed mesh for calculating the motion of flow and the free-surface. This Eulerian approach enables the VOF method to use only a small amount of computer memory for simulating sloshing problems with complicated free-surface contours. The VOF function, representing the volume fraction of a cell occupied by the fluid, is calculated for each cells, which gives the locating of the free-surface filling any some fraction of cells with fluid. Using SOLA-VOF method, the study describes visualization about simulation of the sloshing of flows and damping effect by baffle. Translation and pitching motion of the forms have been investigated The time-dependent changes of free-surface height are used for visualization subject to several conditions such as fluid height horizontal acceleration, sinusoidal motion, and viscosity. The free-surface heights were used for comparing wall-force, which is caused by sloshing of flows. Baffle was Installed to reduce the force on the wall by sloshing of flows. Damping effects was extensively expressed under the conditions such as baffle shape and position.

This paper presents an effective numerical method for analyzing three-dimensional unsteady conjugate heat transfer problems of a curved pipe subjected to infernally thermal stratification. In the present numerical analyses, the thermally stratified flows in the pipe are simulated using the standard $k-{\varepsilon}$turbulent model and the unsteady conjugate heat transfer is treated numerically with a simple and convenient numerical technique. The unsteady conjugate heat transfer analysis method is implemented in a finite volume thermal-hydraulic computer code based on a non-staggered grid arrangement, SIMPLEC algorithm and higher-order bounded convection scheme. Numerical calculations have been performed far the two cases of thermally stratified pipe flows where the surging directions are opposite each other i.e. In-surge and out-surge. The results show that the present numerical analysis method is effective to solve the unsteady flow and conjugate heat transfer in a curved pipe subjected to infernally thermal stratification.

The present study deals with CFD analysis of 'The vortex generators on plastic plate heat exchanger'. When a vortex generator is placed on the heat transfer surface, the flow gets more complex because it entails complicated three-dimensional flows such as separation, reattachment, and recirculation. CFX-5.4, a commercial code utilizing unstructured mesh, has been used as a computational method for solving RANS(Reynolds-Averaged Wavier-Stokes) equations, and the applied turbulence model is $k-{\varepsilon}$ model. In addition, those computational analyses were implemented under various conditions , with or without the vortex generator between two plates, the number, form and the size of vortex generator, and different attack of angle. From the calculated temperature, velocity and pressure distribution, vorticity, wall heat flux and so on under those conditions, this study shows the effect of vortex on heat transfer.

The purpose of this study is to investigate the flow analysis for the optimum design of flat burner using the surface combustion technique. The numerical analysis far flows were compared with the experimental results in order to get a basic database for establishing the integrated numerical method of combustion flows. As results, we knew the current numerical calculation method had confidence with error range under $6{\%}$ in compare with the experiments.

The present study deals with CFD analysis of a plastic heat exchanger with corrugated wall. This exchanger has sinusoidal corrugations, and the flow through the exchanger is three dimensional. In addition, CFX-5.4, a commercial code utilizing unstructured mesh, was used as a computational method for solving RANS(Reynolds-Averaged Navier-Stokes) equations, and the applied turbulence model is $k-{\varepsilon}$ model. The factors to affect the efficiency of a plastic heat exchanger are heat conductivity, flow characteristics and so on. For those two factors, heat conductivity is fixed by the wall material. Therefore, the How along the corrugation affects the efficiency more, provided the same material. In conclusion, the heat transfer enhancement of a plastic heat exchanger with corrugated wall can be recognized from the flow characteristics such as velocity streamline, local heat transfer coefficient, velocity contour, and pressure contour. To confirm the results, both of the measured and the computational data for pressure loss were compared with each other, and they were identical.

In general, Liquid Injection Thrust Vector Control(LITVC) is accomplished by injecting a liquid into the supersonic exhaust flow through holes in the wall of the propulsion nozzle. This injection flow field is highly complicated and detailed flow physics associated with the secondary flow injection should be known far the practical design and use of the LITVC system. The present study aims at understanding the LTTVC flow field and obtaining fundamental design parameters for LITVC. The experimentations were performed in a supersonic blow-down wind tunnel. Compressed, dry air was used for both the main exhaust and injection flows but the pressures of these two flows were controlled independently. The location of the injection holes was changed and the pressures of the two streams were also changed between 2.0 and 15.0 bar. The effectiveness of LITVC was discussed in details using the results of the pressure measurements and flow visualizations

Very frequently the compressible flow in an extinction nozzle of gas circuit breaker is simulated under no arc assumption, which can be reasonable for both high and low current breakings. In the present study, computations are performed to investigate the major features of the compressible flows inside the arc extinction nozzle of gas circuit breaker. A fully implicit finite volume scheme is applied to solve the two-dimensional, steady, compressible, Wavier-Stokes equations. The computed results are validated with the previous experimental data available. Several types of turbulence models are explored to reasonably predict the complicated flows inside the arc extinction nozzle. The obtained results show that the shock wave boundary layer interaction inside the nozzle significantly influences the whole performance of the gas breaker.

Supersonic jet flow has been applied to many various industrial applications of manufacturing fields. Such a supersonic jet is generally classified by three flow patterns, depending on the flow state at nozzle exit, that is, under-, correctly- and over-expanded flows. Of these three flows, the correctly-expanded supersonic jet is most frequently used since it provides a maximum performance of a flow device. However detailed information on what conditions are the Jet correctly expanded at the exit of nozzle is not well known. In the current study, computations are applied to the axisymmetric, compressible, Navier-Stokes equations. The design Mach number used are 2.0,1.2 and 2.6. The computational results obtained are compared with the previous experimental ones. A theoretical analysis is conducted to predict the major features of the correctly-expanded jet. The results show that the jet core length is increased as Mach number is increased.

Compressible gas flow through a convergent-divergent nozzle is choked at the nozzle throat under a certain critical pressure ratio, and then being no longer dependent on the pressure change in the downstream flow field. In practical, the flow field at the divergent part of the critical nozzle can affect the effective critical pressure ratio. In order to investigate details of flow field through a critical nozzle, the present study solves the axisymmetric, compressible, Wavier-Stokes equations. The diameter of the nozzle throat is D=8.26mm and the half angle of the diffuser is changed between $2^{\circ}\;and\;10^{\circ}$ Computational results are compared with the previous experimental ones. The results obtained show that the divergence angle is significantly influences the critical pressure ratio and the present computations predict the experimented discharge coefficient and critical pressure ratio with a good accuracy. It is also found that a nozzle with the half angle of $4^{\circ}$ nearly predicts the theoretical critical pressure ratio.

High-speed ejection of burnt gases from the resonator cavity is essential for performance optimization of the chemical laser system. Additionally, to maintain the population of lasing species at a level for maximum optical power, the pressure within the cavity must be of order of 10 torr. In the present study, a small-scale ejector was designed and built for parametric study of its performance. High-pressure air was used as a motive gas. Measurements include schlieren visualization and pressure distribution trace near the ejector nozzle and along the diffuser downstream of the ejector. preliminary tests showed performance of the ejector is a function of parameters including mass flow rate and stagnation pressure of the motive gas, ejector nozzle area ratio, throat area of the diffuser downstream of the ejector.

An experimental model and a conceptual model are investigated in this paper with both shock tube experiment and numerical technique. The shock-vortex interaction generated by this model is visualized with various methods: holographic interferometry, shodowgraphy, and numerical computation. In terms of shock dynamics, there are two meaningful physics in the present problem. They are reflective wave from the slip layer at the vortex edge and transmitted shock penetrating the vortex core. The discussion in this study is mainly focused on the two kinds of waves contributing to the quadrupolar pressure distribution around the vortex center during the interaction.

Recently the micro flows have been received much attention in the applications with regard to Micro Electro Mechanical Systems(MEMS). Such flows are governed by relatively large viscous forces, compared with inetia force, and frequently specified by slip wall boundary conditions. In the present study, computations are applied to investigate the compressible viscous flows through an extremely small channel, and theoretical analyses are conducted using the Fanno flow theory. The axisymmetic, compressible, Wavier-Stokes equations are numerically solved by a fully implicit finite implicit method. The predicted results are validated with previous experimental data available. The results obtained show that for small Reynolds numbers, the viscous frictional force significantly influences the compressible micro channel flows.

In general, a single gas flow through a converging nozzle is choked when the pressure communications between the downstream and upstream flowfields are broken by the sonic condition of Mach number, M=1. A similar phenomenon may occur In two streams of different stagnation properties flowing side by side in a converging nozzle. In this case, the limiting condition of M=1 for flow choking is no longer applied to such a compound compressible flow. The compound choking phenomenon can be explained by means of a compound sound wave at the nozzle exit. In order to detail the flow characteristics involved in such a compound choking of the two streams, the two-dimensional, compressible, Wavier-Stokes equations have been solved using a fully implicit finite volume method and compared with the results of the one-dimensional theoretical analysis. The computational and theoretical results show that the compound sound wave can reasonably explain the compound choking phenomenon of the two streams in the convergent flow channel.

The twin impulse wave leads to very complicated flow fields, such as Mach stem, spherical waves, and vortex ring. The twin impulse wave discharged from the exits of the two tubes placed in parallel is investigated to understand detailed flow physics associated with the twin impulse wave, compared with those in a single impulse wave. In the current study, the merging phenomena and propagation characteristics of the impulse waves are investigated using a shock tube experiment and by numerical computations. The Harten-Yee's total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The Mach number $M_{s}$, of incident shock wave is changed below 1.5 and the distance between two-parallel tubes, L/d, is changed from 1.2 to 4.0. In the shock tube experiment, the twin impulse waves are visualized by a Schlieren optical system for the purpose of validation of computational work. The results obtained show that on the symmetric axis between two parallel tubes, the peak pressure produced by the twin-impulse waves and its location strongly depend upon the distance between two parallel tubes, L/d and the incident shock Mach number, $M_{s}$. The predicted Schlieren images represent the measured twin-impulse wave with a good accuracy.

In general, TVC(Thermal Vapor Compressor) is used to boost/compress a low pressure vapor to a higher pressure for further utilization. The one-dimensional method is simple and reasonably accurate, but cannot realize the detail as like the back flow and recirculation in the mixing chamber, viscous shear effect, and etc. In this study, the axisymmetric How simulations have been performed to reveal the detailed flow characteristics for the various ejector shapes. The Navier-Stokes and energy equations are solved together with the continuity equation In the compressible flow fields. The standard $k-{\epsilon}$ model is selected for the turbulence modeling. The commercial computational fluid dynamic code FLUENT software is used for the simulation. The results contain the entrainment ratio under the various motive, suction and discharge pressure conditions. The numerical results are compared with the experimental data, and the comparison shows the good agreement. The three different flow regimes (double chocking, single chocking and back flow) have been clearly distinguished according to each boundary pressure values. Also the effects of the various shape variables (nozzle position, nozzle outlet diameter, mixing tube diameter, mixing tube converging angle, and etc.) are quantitatively discussed.

New S-shaped aeroelastic mesoflaps are utilized to control normal shock/boundary-layer interactions. New generation of the mesoflaps is designed f3r a better rigidness and a good flow uniformity across the ulteractions. ,Major advantages of the mesoflap system can be a better total pressure recovery downstream of the interactions due to the lambda shock structure over the flap system, and a rehabilitation of the thickened boundary layer due to bleeding through a cavity underneath the flap system. Skin friction has been measured downstream of the interactions, using the laser interferometer skin friction (LISF) meter, which optically detects the rate of thinning of an oil film applied to the test surface. Various flap-thicknesses of the S-shaped mesoflap arrays are tested, and the results are compared to the solid-wall reference case. Overall, not much difference in the level of skin friction is noticed for the S-shaped flap arrays of various thicknesses, and its level is lower than the skin friction downstream of the solid-wall interaction

In the present study, a passive control method, using the porous wall and cavity system, is applied to the shock wave/boundary layer interactions in transonic moist air flow. The two-dimensional, unsteady, compressible Navier-Stokes equations, which are fully coupled with a droplet growth equation, are solved by the third-order MUSCL type TVD finite difference scheme. Baldwind-Lomax turbulence model is employed to close the governing equations. In order to investigate the effectiveness of the present control method, the total pressure losses of the flow and the time-dependent behaviour of shock motions are analyzed in detail. The computed results show that the present passive control method considerably reduces the total pressure losses due to the shock/boundary layer interaction in transonic moist air flow and suppresses the unsteady shock wave motions over the airfoil, as well. It is also found that the location of the porous ventilation significantly influences the control effectiveness.

This paper describes the dynamics of the weak shock wave propagating inside some kinds of branched pipe bends. Computations are carried out by solving the two-dimensional, compressible, unsteady Euler Equations. The second-order TVD(Total Variation Diminishing) scheme is employed to discretize the governing equations. For computations, two types of branched pipe($90^{\circ}$ branch,$45^{\circ}$ branch) with a diameter of D are used. The incident normal shock wave is assumed at D upstream of the pipe bend entrance, and its Mach number is changed between 1.1 and 2.4. The flow fields are numerically visualized by using the pressure contours and computed schlieren images. The comparison with the experimental data performed for the purpose of validation of computational work. Reflection and diffraction of the propagating shock wave are clarified. The present computations predicted the experimented flow field with a good accuracy.

When a plane shockwave reflects ken a concave wall, it is focused at a certain location, resulting in extremely high local pressure and temperature. This focusing is due to a nonlinear phenomenon of shock wave. The focusing phenomenon has been extensively applied to many diverse folds of engineering and medical treatment as well. In the current study, the focusing of shock wave over a reflector is numerically investigated using a CFD method. The Harten-Yee total variation diminishing (TVD) scheme is used to solve the unsteady, two-dimensional, compressible, Euler equations. The incident shock wave Mach number $M_{s}\;of\;1.1{\~}l.3$ is applied to the parabolic reflectors with several different depths. Detailed focusing characteristics of the shock wave are investigated in terms of peak pressure, gasdynamic and geometrical foci. The results obtained are compared with the previous experimental results. The results obtained show that the peak pressure of shock wave focusing and its location strongly depend on the magnitude of the incident shock wave and depth of parabolic reflector. It is also found that depending up on the depth of parabolic reflector, the weak shock wave focusing process can classified into three distinct patterns : the reflected shock waves do not intersect each other before and after focusing, the reflected shock waves do not intersect each other before focusing, but intersect after focusing, and the reflected shock waves intersect each other before and after focusing. The predicted Schlieren images represent the measured shock wave focusing with a good accuracy.

The purpose of this study is to develop a tool for the analysis of the cavitating flow around trans-cavitating marine propellers. In this study, a linearized super-cavitation theory was applied in order to analyze the performance of the 2-dimensional foils. The numerical results correlated very well with experimental data. The trans-cavitating propellers, manufactured and tasted in KRISO, is selected to validate the lifting surface procedure. For a TCP with a Johnson's five term section, the comparison between the numerical prediction and experiments is fairly good and promising. The new lifting surface procedure, developed and validated with 2-D foils and a TCP, is generally considered applicable to the practical design of the trans-cavitating propeller with Johnson's five term section

The study on the design of pre-swirl stator vanes is performed. The pre-swirl stator vanes is an energy-saving device to improve propulsive performance by providing pre-swirl to the propeller inflow. The theoretical background and the design conditions for pre-swirl stator vanes are presented. The flow characteristics around the pre-swirl stator vanes attached ship hull are analyzed through the experimental method. The technique to determine the optimum location, angle and the number of stator vane is investigated and applied it to 310,000 TBW VLCC The flow velocities are measured using 5-hole Pilot tubes at the condition with and without a propeller.

CFD calculations are carried out to investigate the turbulent flow characteristics inside the duct of marine waterjet propulsors. The Reynolds-averaged Wavier-Stokes equations are solved using a finite-volume method. Standard $k-{\varepsilon}$ model and realizable $k-{\varepsilon}$ model are evaluated with an existing experimental data. Multi-block grid topology is adopted to describe the details of complex duct geometry. The present numerical methods are applied to the preliminary duct design of new waterjet propulsor system. Four different influx conditions are simulated to find out pressure and velocity distribution inside the intake duct. Attention is also paid upon the possible flow separation inside the waterjet duct. It is found that CFD tools can be used for the initial evaluation of inflow condition into the impeller of waterjet propulsor system.

The objective of present paper is to apply a stereoscopic PIV(Particle Image Velocimetry) techiique for measuring the 3 dimensional flow structure of turbulent wake behind a marine propeller with 5 blades. It is essential to measure 3-components velocity fields for the investigation of complicated near-wake behind the propeller. The out-of-plane velocity component was measured using the particle images captured by two CCD cameras in the angular displacement configuration.400 instantaneous velocity fields were measured for each of few different blade phases of $0^{\circ},\;18^{\circ},\;36^{\circ}\;and\;54^{\circ}$. They were ensemble averaged to investigate the spatial evolution of the propeller wake in the region ranged from the trailing edge to the region of one propeller diameter(D) downstream. The phase-averaged velocity fields show the viscous wake formed by the boundary layers developed along the blade surfaces. Tip vortices were formed periodically and the slipstream contraction occurs in the near-wake region. The out-of-plane velocity component has large values at the tip and trailing votices. With going downstream, the axial turbulence intensity and the strength of tip vortices were decreased due to the visous dissipation, turbulence diffusion and blade-to-blade interaction. The blade wake traveling at higher speed with respect to the tip vortex overtakes and interacts with tip vortices formed from the previous blade. Tip vortices are separated from the wake and show oscillating trajectory

In this paper, the submarine model, called DARPA SUBOFF model, has been numerically analyzed to investigate the aerodynamic forces variation in terms of angle of attacks and yaw angles. The SUBOFF model is consisted of the three parts : axisymmetric body, fairwater, and four symmetric stern appendages. Three dimensional unsteady incompressible Wavier-Stokes equation was used on curvilinear multi-block grid system. To validate the present code, the SUBOFF tare hull and an ellipsoid at angle of attacks of $10^{\circ}\;and\;30^{\circ}$ were simulated and a good agreement with experiments was obtained. After the code validation, the flows over SUBOFF model were simulated at three different angle of attacks and yaw angles. The variation of aerodynamic forces in terms of angle of attack and yaw angle were calculated. Also, to understand the flow features around a submarine with variation of yaw and attack angle, the pressure contours and streamlines were plotted.

The relativistic theory has not been properly taken up by the marine hydrodynamicists. To take on a relativistic view, we confine ourselves to a simple vector case of a wave train in spacetime, to be shown to represent a sound wave or a surface wave, and bring in an observer who is travelling on another platform. We are interested in relative position of each event on these two worldlines. It, then, will be shown that the velocity, the acceleration, the encounter frequency, the group velocity, and the time and the space distance between the wave and the observer on the worldlines should all be derivable in principle. This is interpreted to mean that we really have the relativistic events taking place with different values of time dilation in the sense of 'spacetime', and that the well-known ${\lceil}special Theory of Relativity{\rfloor}$ applies just as well in hydrodynamic waves.

In this paper, roll damping coefficients for a non-conventional cross section, which is herein named as 'step' model, are investigated numerically and experimentally. Experiments are extensively carried out to estimate the roll damping coefficients. Numerical estimations are also made with the help of numerical codes. For convenience, the roll damping is divided into wave-making component and viscous component. The wave-making component is determined using a potential code and the viscous component using a viscous flow code, in which the fluid domain is taken as unbounded. In order to validate the present approach, a typical cross section with bilge is considered and our results are compared with published data. The comparison shows a good agreement qualitatively. For the step model, numerical results are compared well with experimental data besides some quantitative discrepancies at a certain range of frequency. It is thought that the discrepancy might be caused by the ignorance of the free surface in viscous computations. It is found in the case of the step model that not only the viscous component but also the wave component increases considerably compared to the section with bilge.

This paper describes the hydrodynamic characteristics of a test-bed AUV SNUUV-I constructed at Seoul National University. The main purpose of the AUV is to carry out fundamental control and hydrodynamic experiments. Its configuration is basically a long cylinder of 1.35m in length and 0.25m in diameter with delta-type wings near its rear end. On the edge of each wing, a thruster of 1/4HP is mounted, which is used for both drive and turn the vehicle for horizontal movement as the output control power is varied. A pair of control surfaces installed near its font part generates pitch moments for vertical movement. The 6 DOF mathematical model of SNUUV-I contains hydrodynamic forces and moments expressed in terms of a set of hydrodynamic coefficients. These coefficients can be classified into linear damping coefficients, linear inertial coefficients and nonlinear damping coefficients. It is important to estimate the exact value of these coefficients to control the vehicle precisely. Among these, the linear coefficients are known to affect the motion of the vehicle dominantly. The linear damping coefficients are estimated by using Extended Kalman Filter. The responses of the vehicle to input signals are used to estimate the hydrodynamic coefficients, which can be inferred from output signals measured from an IMU (inertial motion unit) sensor, while the linear inertial coefficients are calculated by a potential code. By using these coefficients estimated as described above, a simulation program is constructed using Matlab.

Numerous studies have been performed to develop an analytical model that can predict the reflection of regular or irregular waves from a perforated-wall caisson breakwater. Though such irregular wave models as Suh et at. (2001) become available, regular wave models are still in extensive use because of their simplicity. In the present study, using the regular wave model of Fuggazza and Natale(1992), the reflection of irregular waves from a perforated-wall caisson breakwater was calculated in several different methods. First, the regular wave model was re-validated by the hydraulic model tests. Though the model somewhat over-predicted the reflection coefficients at larger values and under-predicted them at smaller values, overall agreement was pretty good between calculation and measurement. Then, the regular wave model was applied to calculate the irregular wave reflection in the experiments of Suh et at.(2001) and Bennett et al. (1992). In applying the regular wave model to irregular wave reflection, several different methods were used. The results showed that it is the most reasonable to use the regular wave model repeatedly for each frequency component of the irregular wave specuum with the root-mean-squared wave height for all the frequencies .

Since Berkhoff proposed the mild-slope equation in 1972, it has widely been used for calculation of shallow water wave transformation. Recently, it was extended to give an extended mild-slope equation, which includes the bottom slope squared term and bottom curvature term so as to be capable of modeling wave transformation on rapidly varying topography. These equations were derived by integrating the Laplace equation vertically. In the present study, we develop a finite element model to solve the Laplace equation directly while keeping the same computational efficiency as the mild-slope equation. This model assumes the vertical variation of wave potential as a cosine hyperbolic function as done in the derivation of the mild-slope equation, and the Galerkin method is used to discretize . The computational domain was discretized with proper finite elements, while the radiation condition at infinity was treated by introducing the concept of an infinite element. The upper boundary condition can be either free surface or a solid structure. The applicability of the developed model was verified through example analyses of two-dimensional wave reflection and transmission. .

In this paper, flow patterns around two vertical cylinders were investigated numerically and experimentally Flow patterns between cylinders are very complex and interative. Changing gaps between cylinders the flow patterns are measured and numerically studied at a fixed comming velocity. The experiment has beeen conducted in circulating water channel with PIV system and manometer, and numerical analysis has been made by F.D.M. and multi block method. The results can be applied in the understanding and design of multiple pile array structures.

The assessment of PIV to measure the mean velocity and turbulence was carried out over a train of fixed two-dimensional dunes. The agreement between the PIV and LDV is good enough even in regions of flow reversals and high shear. Though limited in the wall normal direction field-of-view, PIV provides instantaneous flow fields, which reveal the complex nature of flow over dunes, as well as more sophisticated analyses such as two-point space correlation and quadrant analysis with a reasonable accuracy.

A technique of mixing enhancement in an underexpanded sonic round jet is studied with fully expanded jet Mach number 1.5. Tonal sound, jet screech can be produced at some underexpansion pressure ratio in a sonic jet. Since the jet screech excites the initial Jet shear layer to change the flow, a reflector which focuses the jet screech near the nozzle lip is designed. The reflector has an elliptic shape of which two foci are located near the nozzle lip and the jet screech source region. Jet screech tone near the nozzle lip increases with the elliptic reflector and spreading of the jet largely increases. It is concluded that mixing enhancement of the jet with the elliptic reflector is attributed to large scale structures which are initially excited by the increased jet screech.

Supersonic jet impingement on a flat plate has been investigated to show the flow physics for different jet heights and to demonstrate the adequacy of the characteristics-based flux-difference Wavier-Stokes code Current study also compares the steady-state solutions obtained with variable CFL number for different grid spacing with the time-accurate unsteady solutions using the inner iterations, displaying a good agreement between the two sets of numerical solutions. The unsteady nature of wall fluctuations due to bouncing of the plate shock is also uncovered for high pressure ratios. The methodology is then applied to a complex vertical launcher system where the jet plume hits the bottom wail, deflects into the plenum and eventually exits through the vertical uptake. Flow structures within vertical launcher system are captured and solutions are partially verified against the flight test data. Present jet impingement study thus shows the usefulness of CFD in designing a complex structure and predicting flow behavior within such a system.

For the purpose of Thermal Protection Material design problem, a numerical analysis of axisymmetric high temperature supersonic impinging jet flows of exhaust gas from combustor on curved surfaces has been accomplished. A modifed CSCM Upwind Navier-Stokes method which is able to cure the carbuncle Phenomena has been developed to study strong shock wave structure and thermodynamic wall properties such as pressure and heat transfer rate on various curved surfaces. The results show that the maximum heat transfer rate which is the most important parameter affecting thermo-chemical surface ablation on the plate did not occur at the center of jet impingement, but rather on a circle slightly away from the center of impingement and the shear stress distribution alone the wall is similar to the wall heat transfer late distribution.

The high-pressurized spray nozzle is used f3r special washing and cutting with strong impact force. The performance of this nozzle, which focused on spray penetration and radial dispersion, was mainly investigated to maximize the momentum and minimize the flow loss. Hence, our experimental research was conducted by changing the aspect ratio ranging from 0 to 3 with nozzle outlet of 1.1. The spray trajectory far high-pressurized water was experimentally investigated using PDPA diagnostics, which was available at spray downstream region. As the spray at upstream near the nozzle exit did not show the improved disintegration. The results showed empirical correlation with regard to non-dimensional axial velocity distribution, spray penetration, and radial spreading rate with photographic visualization.

In the present study, we perform LES of turbulent flow and temperature fields in a circular impinging jet at Re=5000 for two cases of H/D=2 and 6 (H denotes the distance between the jet exit and flat plate, and D does the diameter of the jet exit). In the case of H/D=2, the regular vortical structures observed in free jet do not exist because of the smaller distance than the potential core. The Nusselt number on the wall is largest at $r/D{\cong}10.67$ where vortex rings Impinge. At $r/D{\cong}1.5{\~}2.0$, the vortex rings induce the secondary vortices, resulting in a secondary peak in the Nusselt number there. In the case of H/D=6, the vortex rings change into three-dimensional vortical structures and the small-scale vortices impinge on the flat plate. The increase of turbulent intensity due to small-scale vortices results in the largest Nusselt number at the stagnation point.

Large eddy simulation of a circular jet at the Reynolds number of 10000 is performed to investigate turbulence suppression effect with single frequency excitation at the non-dimensional frequency of 0.017. Instantaneous flow fields show that, with excitation, naturally occurring energetic vortices are suppressed through earlier saturation and breakdown of the shear layer vortices into fine grained turbulence. Due to the excitation, the Reynolds stresses are larger for the excited case near the jet and turbulence suppression begins afterward. The Reynolds normal stresses show largest suppression in the shear layer near the jet and in the centerline further downstream, while the Reynolds shear stress shows largest suppression in the shear layer at all the downstream locations.

Vortex lock-on or resonance in the flow behind a circular cylinder is visualized by a time-resolved PIV when a single frequency oscillation is superimposed on the mean incident velocity. Measurements are made of the $K{\'{a}}rm{\'{a}}n$ vortices in the wake-transition regime at the Reynolds number 360. Basically, natural shedding state is observed to compare with lock-on state. Wake motion by the change of the shedding frequency of lock-on state is investigated. When lock-on occurs, the vortex shedding frequency is found to be half the oscillation frequency as expected from previous experiments. The physical flow phenomena of natural shedding and lock-on states are analyzed with physical parameters of recirculation and vortex formation region. Consequently, it is found that the change of wake bubble plays an important role in the flow at the lock-on state. Vortex formation region is also actively changed like recirculation region as the lock-on occurs. Therefore, it is deduced that the recirculation region is closely related with the vortex formation region.

Turbulent properties of the wake of a circular cylinder were measured The diameter of the cylinder is l0mm and the Reynolds number is 420. A new 3D-PTY system was constructed and a genetic algorithm (GA) was introduced in order to increase the number of instantaneous three-dimensional velocity vectors. In the GA two fitness functions were introduced in order to enhance the correspondences of the particles. The measurement system consists of three CCD cameras, Ar-ion laser, an image grabber and a host computer. More than 3000 instantaneous three-dimensional velocity vectors were obtained by the system. The database of the turbulent properties of the circular cylinder was constructed by the constructed 3D-PTV system.

The determination of thrust is very important in hypersonic air-breathing propulsion design and evaluation. Because of the short flow-residence time in the combustor, the evaluation of engine performance is strongly influenced upon the engine thrust. Conventional methods to determine the thrust is using thrust stand or force measurement system. However, these methods cannot be applied to the case where thrust stands are impractical, such as free jet testing of engines, and model combustor. With this reason, the thrust determination method from measured pilot pressure is considered and evaluated.

Liquid rocket engines using liquefied natural gas (LNG) or methane as a fuel is known to have several good characteristics, such as high specific impulse compared to other hydrocarbon fuels, environment-friendly exhaust gas, low production cost, and re-usability with low soot generation in the cooling channel. In this study, experimental combustion chambers capable of using LNC and $CH_{4}$ are being researched through experimental firing tests, and within easy range of eyes' inspection, there are the periodical existence of soot or discoloration in the chamber wall surface. This result means that mixture ratio of oxidizer and fuel fluctuates periodically between outer-row injectors in the mixing head in the circumferential direction. Therefore, based on this phenomenon, the variation of mixture ratio near the chamber wall caused by the spill pattern of a shear coaxial injector was analyzed quantitatively and the thermal heat flux Into the cooling channel is modified. Then, the calculated and modified results are compared with the measured ones.

Application of combustion stabilization devices such as baffle and acoustic cavity to liquid propellant rocket engine is investigated to suppress high-frequency combustion instability, i.e., acoustic instability. First, these damping devices are designed based on linear damping theory. As a principal design parameter, damping factor is considered and calculated numerically in the chambers with various specifications of these devices. Next, the unbaffled chambers with/without acoustic cavities are tested experimentally for several operating conditions. The unbaffled chamber shows the specific stability characteristics depending on the operating condition and has small dynamic stability margin. The most hazardous frequency is clearly identified through Fast Fourier Transform. As a result, the acoustic cavity with the present design has little stabilization effect in this specific chamber. Finally, stability rating tests are conducted with the baffled chamber, where evident combustion stabilization is observed, which indicates sufficient damping effect. Thrust loss caused by baffle installation is about $2{\%}$.

The KSR-III developed by KARI is the first rocket vehicle which is adopting the liquid propellant rocket engine system in Korea. Not only the engine itself, but also the propellant feeding system is one of the most important component in liquid rocket vehicle. In this paper, the authors are intended to introduce the multi-purpose test facility(PTA-II Test Facility) which is constructed for the variety of tests on KSR-III feeding system(single component tests, verification tests, cold flow tests and combustion tests). With the results of these tests, we can identify the characteristics of rocket feeding system and decide the optimum setting values of feeding system for the successful flight.

Combustion instability, which is one of the most undesirable phenomena in the development of liquid Propellant rocket engine, can cause serious damage to the rocket itself, and must be evaded by all means. Unfortunately, KSR-III rocket went through the combustion instability during engine start at the propulsion test article No.2. To resolve the problem, time sequence has been changed, and the baffle system has been applied. In consequence of the change, stable combustion was achieved.

This paper is investigated construction of the Scramjet test facility and test method of Scramjet engine combustor model. Scramiet engine combustor model test was performed at Lab C-16BK CIAM (Central Institute of Aviation Motors) at Tyraevo in Moscow. The velocity of flow in the combustion chamber equal to Mach number 2.49 with single hole fuel spray nozzle injector and test duration equal to 7 seconds. Therefore In this paper is showed high altitude test method of Scramjet combustor model and the proper structure of combustor with single hole fuel spray nozzle.

The purpose of this study is to investigate the combustion characteristics of the flat burner using metal fiber. By changing the excess air ratio which affects strongly to the combustion characteristics, the flame stability range in the metal fiber were found and the range of the blue flame mode and radiant mode were distinguished and $NOx\;and\;O_{2}$ distributions in combustor were measured at two conditions which represents the blue flame mode and radiant mode.

Recently, gas turbines for power generation adopt multistage DLN(Dry Low NOx) type combustion, where diffusion combustion is applied at low load and, with increase in load, the combustion mode is changed to lean premixed combustion to reduce NOx emissive concentration. However, during the mode changeover from diffusion to premixed flame, unfavorable phenomena, such as flashback, high amplitude combustion oscillations, or thermal damage of combustor parts could frequently occur. In the present study, to apply for the analysis of such unfavorable phenomena, three-dimensional CFD investigations are carried out to compare the detailed flow characteristics and temperature distribution inside the gas turbine combustor before and after combustion mode changeover. The fuel considered here is pure methane gas. A standard $k-{\varepsilon}$ turbulence model with wall function and a P-N type radiation heat transfer model, have been utilized. To analyze the complex geometric effects of combustor parts on combustion characteristics, fuel nozzles, a swirl vane f3r fuel-air mixing, and cooling air holes on the combustor liner wall, are included in this simulation.

The mixing of fuel and oxidizer and each mass distribution of unlike split triplet(F-O-O-F) injector for liquid rocket which are known to affect the combustion efficiency significantly, has been investigated using PLLIF technique. Experiment is conducted to investigate the effect of mixture ratio(MiR), momentum ratio(MoR) and impinging angle on mixing efficiency. The mixing efficiency, which is introduced by Rupe, gives the global chracteristic of spray mixing. Experiment far comparison with triplet injector(F-O-F), which has same momentum ratio and impinging angle with split triplet injector is conducted.

'Backhole' is an extra empty volume where is located behind the tangential entries at the rear par of the vortex chamber in the swirl coaxial injector. With the backhole, there are three major hydraulic characteristics. First, mass flow rate is increased about $15{\%}$ compared with the case without the backhole. Second, with the backhole, the center region of the injected flow has more large volume than that of without the backhole. The last, some range of the cone angle can be controlled by the backhole Experiments are conducted by using a PDPA apparatus, a mechanical patternator, stroboscopic photography and etc. With the backhole, based on cold-flow tests, the model swirl injector has some Improvement in its performance.

This study describes the external spray characteristics of deflector nozzle such as the breakup procedures of liquid sheet, spray angle, discharge coefficient and bubble behaviors of spray and SMD at deflector nozzle. In order to visualize the spray behaviors shadow graphy technique were used. According to the increase of injection pressure, development of the spray passes through the dribbling, distoted jet, closed bubble due to the contraction by surface tension forces, the bubble opens into hollow tulip shape, and the curved surface straightened to form a conical sheet like as the simplex swirl atomizer. Spray cone angle was nearly 90 deg. Variations of SMD were examined in order to describe the dependency of SMD on the injection pressure and orifice diameter. The shape of deflector and oriffice diameter had an effect on the discharge coefficient.

In the analysis of the mass distribution and SMD (Sauter Mean Diameter), planar laser imaging technique Is a convenient and useful one when compared to the mechanical patternator or PDPA (Phase Doppler Particle Analyzer). But for the accurate usage of this laser technique, it is required to compensate the signal attenuation and to find the attenuation coefficients. In the present research, we considered effects of the spray distributions on the attenuation coefficient and improved the accuracy of planar laser imaging technique by the reduction of errors issuing from the signal attenuation.

Characteristics of the twin fluid spray with ultrasonic forcing are examined in order to obtain the high efficiency of cold fog spray of the automatic pest control machine which has been widely used for the equipment cultivation recently. An electrostrictive vibrator of PZT BLT and a magnetostrictive vibrator of ${\pi}type$ with 28 kHz are applied as the ultrasonic transducer. All experiments are made and observed in 4 methods of spray ; a conventional spray method without ultrasonic forcing, an indirect vibration method with ultrasonic forcing, an improving duality method by ultrasonic forced within liquid, and a combined use method with both of the indirect vibration method and the improving quality method. In results, It was clarified that the ultrasonic effects the atomization of spray droplets and its efficiency is about $10{\%}$ and especially much more in the case of the combined use method.

Liquid atomization by means of a spinning cup is widely used as a device for combustion, in cooling and spray drying. In this study, the blowing type rotary atomizer was experimental carried out the investigations on the characteristics of the blowing type rotary atomizer which is an air flow energy of blower instead of an electric motor most commonly used to a driven energy. The analysis on the rotary cup speed, air velocity with the blower conditions was performed and also the drop size was measured using LDPA. It was tried to analyzed on air-nozzle size and liquid flowrate as the result. It was found that the increase of the relative velocity between liquid and air improve significantly atomization liquid, and decrease of the liquid flowrate improved the maximum drop size though the mean drop size is really the same.

A two-dimensional numerical method for inviscid two-fluid flows with evolution of density interface is developed, and an initially stationary two-dimensional fluid sheet surrounded by another fluid is studied. The Interface between two fluids is modeled as a vertex sheet, and the flow field u÷th the evolution of interface is solved by using vortex-in-cell/front-tracking method. The edge of the sheet Is pulled back into the sheet due to surface tension and a blob is formed at the edge. This blob and fluid sheet are connected by a thin neck. In the inviscid limit, such process of the blob and neck formation is examined in detail and their kinematic characteristics are summarized with dimensionless parameters. The edge recedes at $V=1.06({\sigma}/{\rho}h)^{0.5}$ and the capillary wave Propagating into the fluid sheet must be considered for bettor understanding of the edge receding.

In this paper a systematic design approach to determine the optimum size (height) of circulating fluidized bed heat exchanger for exhaust gas heat recovery is prososed. Unlike the convensional heat exchangers where the length of the heat exchanger section is not very much emphasized, the vertical length of heat exchanger tube in the case of fluidized bed heat exchangers is important because this length determines the time interval during which particles reside and transfer heat in the heat exchanger section. For particles initial conditions are nearly stationary, accelerating particles motion should be considered rather than simply assuming fully developed condition. A way to estimate optimum tube length at different fluid velocity and particle sizes is suggested based on the required conditioning time for heat transfer from the flue gas to solid particles.

A grease filter is used to remove grease generated from a cooking appliance in a kitchen. This numerical study has been conducted to investigate the particle collection characteristics of a grease filter having nominal flowrate of $100m^{3}/h$. The flow field and particle trajectories in the grease filter with a flow chamber were simulated by using the commercial code of STAR-CD. The pressure drop of a grease filter rapidly increases with increasing the air flowrate. The numerical values of the pressure drop are slightly lower than the experimental values when the air flowrates are 50, 75, and $100m^{3}/h$. The particle collection efficiency of a grease filter increases with increasing the particle diameter and the air flowrate, which means that the inertial impaction is a dominant particle removal mechanism in a grease filter The cut-off diameter of the tested grease filter representing $50-{\%}$ collection efficiency is about $11.6{\mu}m$ for water droplets at $100m^{3}/h$.

The effect of vibration on the viscosity of a shear-thinning fluid was investigated with a newly designed pressure-scanning capillary viscometer. The viscometer was designed to measure non-Newtonian viscosity continuously over a range of shear rates at a time. Low frequency vibration was applied perpendicularly to the direction of the flow. The effect of the transversal vibration was investigated for both Newtonian fluids and non-Newtonian fluids. The experimental results showed that the vibration had no effect on the viscosity of the Newtonian fluids. However, the vibration caused a significant reduction of the shear-thinning fluid viscosity. The viscosity reduction was strongly dependent on both vibration frequency and shear rate. In addition, the viscosity reduction was affected by the amplitude of vibration, and, the bigger amplitude applied, the more viscosity reduction occurred.

Hydraulic valve control the pressure and the How of fluid by the hydraulic oil transfered from pump but the ER fluid consists of solid particles of micrometer in size and insulating oil so in the general hydraulic valve. We invented ER-Valve using ER fluid as working fluid. The ER fluid, working fluid of ER-Valve is a functional fluid to represent the feature of fluid according to strength of electric field. In this research we made our own 4 types of plate type ER-Valve which has same surface but different width and length and then we conducted performance test. We measured flow rate and pressure drop of fluid which is flowing in the ER-Valve according to the electric field strength to conduct this test. We modeling ER-Valve relating to ER-Valve system and yield shear stress according to the strength of electric field. We used the pressure drop according to the strength of electric field by differential pressure gauge in the our own made ER-Valve. This test reviewed experimental the special changes of ER-Fluid in the steady flow condition.

We present numerical and experimental result of the rotating flows inside a rectangular container under a background rotation. In numerical computation, a parallel-computation technique with MPI is implemented. Flow visualization and PIV measurement are also performed to obtain velocity fields at the free surface. Through a series of numerical and experimental works, we aim to clarify, if any, the fundamental reasons of discrepancy between the two-dimensional computation and the experimental measurement, which was detected in the previous study for the same flow model. Specifically, we check if the various assumptions prerequisite for the validity of the classical Ekman pumping law are satisfied for periodic flows under a background rotation.

A theoretical study is made of the steady flow of a compressible fluid in a rapidly rotating finite cylinder. Flow is generated by imposing mechanical and/or thermal disturbances at the rotating endwall disks. Both the Ekman and Rossby numbers are small. A detailed consideration is given to the energy budget for a control volume in the Ekman boundary layer. A combination of physical variables, which is termed the energy contents, consisting of temperature and modified angular momentum, emerges to be relevant. The distinguishing features of a compressible fluid, in contrast to those of an incompressible fluid, are noted. For the Taylor-Proudman column to be sustained, in the interior, it is shown that the net energy transport between the solid disk wall and the interior fluid should vanish. Physical rationalizations are facilitated by resorting to the concept of the afore-stated energy content.

This paper describes a detailed experimental Investigation of heat transfer In a reciprocating rectangular channel fitted with rib structures with particular reference to the design of a piston for marine propulsive diesel engine. The parametric test matrix involves Reynolds number, reciprocating frequency, and reciprocating radius, respectively, in the ranges, $1,000\;{\~}\;6,000,\;1.7\;{\~}\;2.5\;Hz,\;and\;7\;{\~}\;15cm$ with four different rib arrangements. The rib arrangements have considerable influences on the heat transfer in the reciprocating channel due to the modified vortex flow structure. The experimental data confirm that the increases in the heat transfer can be seen in order of Case (a), Case (d), Case (c), and Case (b)

The present experimental and numerical investigations are performed for the characteristics of transitional flow in a concentric annulus with a diameter ratio of 0.52, whose outer cylinder is stationary and inner one rotating. The flow field of an annulus has been numerically solved using a finite volume method. The pressure losses and Skin-friction coefficients have been measured for the fully developed flow of water and $0.2{\%}$ aqueous solution of sodium carboximethy1 cellulose (CMC), respectively at inner cylinder rotational speed of $0{\~}600rpm$. The transitional flow has been examined by the measurement of pressure losses to reveal the relation of the Reynolds and Rossby numbers with the skin-friction coefficients. The occurrence of transition has been checked by the gradient changes of pressure losses and skin-friction coefficients with respect to the Reynolds numbers. Consequently the critical(axial-flow) Reynolds number decreases as the rotational speed increases. Thus, the rotation of inner cylinder promotes the early occurrence of transition due to the onset of Taylor vortices.

A matched asymptotic analysis is conducted for a compressible rotating flow in a cylindrical container when a mechanical and/or a thermal disturbance is imposed on the wall. The system Ekman number is assumed to be very small. The conditions for the Taylor-Proudman column in the interior, which were also given in the companion paper Park & Hyun, 2002) by means of the energy balancing analysis, have been re-derived. The concept of the variable, the energy content $e[{\equiv}T+2 {\alpha}^2 {\gamma}{\nu}]$, is reformulated, and its effectiveness in characterizing the energy transport mechanism is delineated. It is seen that, under the condition of the Taylor-Proudman column, numerous admissible theoretical solutions for interior flow exist with an associated wail boundary condition. Some canonical examples are illustrated with comprehensive physical descriptions. The differential heating problem on the top and bottom endwall disks is revisited by using the concept of the energy content. The results are shown to be in line with the previous findings.

The Flow patterns around two cylinders in various arrangements were studied by a discrete vortex method. The flow for the surface of each cylinder was represented by arranging bound vortices at adequate intervals. The viscous diffusion of fluid was represented by the random walk method. The vortex distributions, streaklines, timelines and velocity vectors around two cylinders were calculated for centre-to-centre pitch ratios of $P/D=1.5 and 2.5$, attack angles of ${\alpha}=0^{\circ},\;30^{\circ},\;60^{\circ}\;and\;90^{\circ}$, and Reynolds number of Re=1200. The results of simulation correspond to the photographs by flow visualization and the flow intereference between two cylinders in various arrangements was clearly visualized by a numerical simulation.

The time-development of the wake vortices of the unsteady viscous flow past a semicircular cylinder is simulated using the vortex particle methods for direct numerical simulations(DNS). The early wake behaviour of the flow behind an impulsively started a semicircualr cylinder is evaluated for a range of Reynolds numbers between 60 and 200 with opposite body configurations respectively. The diffusion scheme based on the particle strength exchange(PSE) is used to account far the viscous effect accurately. And the vorticity generation algorithm to enforce the no-slip boundary conditions is employed. In order to redistribute particles efficiently on the distorted Lagrangian grid the particle distribution technique is adaptively revised, while maintaining the uniform resolution. The results of the simulations are compared to other experimental results.

Planar laser induced fluorescence(PLIF) has been widely used to obtain two dimensional fuel distribution. Preliminary investigation was performed to measure quantitative air excess ratio distribution in an engine fueled with LPG. It is known that fluorescence signal from acetone as a fluorescent tracer is less sensitive to oxygen quenching than other dopants. Acetone was excited by KrF excimer laser (248nm) and its fluorescence image was acquired by ICCD camera with a cut-of filter to suppress Mie scattering from the laser light. For the purpose of quantifying PLIF signal, an image processing method including the correction of laser sheet beam profile was suggested. Raw images were divided by each intensity of laser energy and profile of laser sheet beam. Inhomogeneous fluorescence images scaled with the reference data, which was taken by a calibration process, were converted to air excess ratio distribution. This investigation showed instantaneous quantitative measurement of planar air excess ratio distribution for gaseous fuel.

The effect of distance between ninety degree elbow close to upstream face of orifice plate and orifice plate on discharge coefficient was investigated. The distributions of discharge coefficient and differential pressure caused by elbow and short upstream straight length were examined and modified discharge coefficient was introduced. The results presented in this study could be useful when orifice plate is installed under the condition of simple flow disturbance element and short upstream straight length.

An electromagnetic flowmeter(EMF) was developed and its characteristics were compared with a commercial EMF. The measuring tube of the EMF was basically designed as 100 m diameter. A signal processing circuit was also developed for generating the magnetic field and converting the flow signal to flowrate and flow quantity. To obtain a more stable and reliable flow signal, the double magnetizing frequency was adopted for magnetizing the coil of the EMF For the characterization of EU, the uncertainty of calibrator was estimated within ${\pm}0.5{\%}$, the method of estimation was in accordance with the ISO 17025 recommendation. It was found that the flow signals between the electrodes were about ${\pm}60\;-\;{\pm}300{\mu}V$, which were sufficient for the discrimination of the flowmeter and protecting the noise. The test result against the calibrator, showed the good linearity in the range of $3 m^{3}/h\;and\;70m^{3}/h$. A commercialized design of the EMF will be technically more competitive in domestic and foreign market.

The design of industrial arc plasma systems is still largely based on trial and error although the situation is rapidly improving because of the available computational power at a cost which is still fast coming down. The desire to predict the behavior of arc plasma system, thus reducing the development cost, has been the motivation of arc research. To interrupt fault current, the most enormous duty of a circuit breaker, is achieved by separating two contacts in a interruption medium, $SF_{6}$ gas or air etc., and arc plasma is inevitably established between the contacts. The arc must be controlled and interrupted at an appropriate current zero. In order to analyze arc behavior in $SF_{6}$ gas circuit breakers, a numerical calculation method combined with flow field and electromagnetic field has been developed. The method has been applied to model arc generated in the Aachen nozzle and compared the results with the experimental results. Next, we have simulated the unsteady flow characteristics to be induced by arcing of AC cycle, and conformed that the method can predict arc behavior in account of thermal transport to $SF_{6}$ gas around the arc, such as increase of arc voltage near current zero and dependency of arc radius on arc current to maintain constant arc current density.

Under the assumption of potential flow, free-surface flows around a 2-dimensional hydrofoil are calculated by high-order spoctral/boundary-integral method. This method is one of the most efficient numerical methods by which the nonlinear interactions between hydrofoil and free-surface can be simulated in time-domain. Comparisons of the calculated free-surface profiles with other experimental results show relatively good agreements. As another example, free-surface flow generated by the heaving and translating hydrofoil is calculated and discussed.

The drift force acting on a floating OWC chamber in waves is studied taking account of fluctuating air pressure in the air chamber. A velocity potential in the water due to the free surface oscillating pressure patch is added to the conventional radiation-diffraction potential problem. The potential problem inside the chamber is formulated by making use of the Green integral equation associated with the Rankine Green function while the outer problem with the Kelvin Green function. The drift forces as well as the chamber motions are calculated taking account of the air pressure in the chamber.

This paper presents the results of a study investigating methods of interpretation of wave directionality based on wavelet transforms. Two-dimensional discrete wavelet was used for the analysis. The proposed scheme utilizes a single frame of ocean waves to detect their directionality. This fact is striking considering the fact that traditional methods require long time histories of ocean wave elevation measured at various locations. The developed schemes were applied to the data generated from numerical simulations and video images to test the efficiency of the proposed scheme in detecting the directionality of ocean waves.

The free surface flow problem has been one of the most interesting and challenging topic in the area of the naval ship hydrodynamics and ocean engineering field. The problem has been treated mainly in the scope of the potential theory and its governing equation is well known Laplace equation. But in general, the exact solution to the problem is very difficult to obtain because of the nonlinearlity of the free surface boundary condition. Thus the linearized free surface problem has been treated often in the past. But as the computational power increases, there is a growing trend to solve the fully nonlinear free surface problem numerically. In the present study, a time-dependent finite element method is developed to solve the problem. The initial-boundary problem is formulated and replaced by an equivalent variational formulation. Specifically, the computations are made for a highly nonlinear flow phenomena behind a transom stern ship and a vertical strut piercing the free surface.

The researches of a two-phase atomizers have been carried out in the field of automotive and aerospace industries in order to improve the atomization performance of the liquid droplets ejecting from these nozzles. The smaller droplets have the advantages of the reduction of environmental pollution matter and effective use of energy through the improvement of heat and mass transfer efficiency. Thus, to propose the basic information of two-phase flow, an internal mixing atomizer was designed, its shape factor was 0.6 and the liquid feeding hole was positioned at the center of the mixing tube which was used to mix the air and liquid. The experimental work was performed in the field after the nozzle exit orifice. The measurement of the liquid droplets was made by PDPA system. This system can measure the velocity and size of the droplets simultaneously. The number of the droplets used in this calculation was set to 10,000. The flow patterns were regulated by ALR (Air to Liquid mass Ratio). ALR was varied from 0.1024 to 0.3238 depending on the mass flow rate of the air. The analysis of sampling data was mainly focused on the spray characteristics such as flow characteristics distributions, half-width of spray, RMS, and turbulent kinetic energy with ALR.

In many industrial applications, multiphase flow analysis is the norm rather than an exception as compared to more-conventional single-phase investigation. This paper describes the implementation of the multiphase flow simulation capability in the general purpose CFD software AVL FIRE/SWIFT. The governing equations are discretized based on a finite volume method (FVM) suitable fur very complex geometry, The pressure field is obtained using the SIMPLE algorithm. Depending on the characteristics of the multiphase flow to be examined, the user can choose either the two-fluid model or an explicit interface-tracking model based on the Volume-of-Fluid approach. For truly 'multi'-phase flow problems, it is also possible to apply a hybrid model where certain phases are explicitly tracked while the other phases are handled by the two fluid model. In order to demonstrate the capability of the method, applications to the Taylor bubble flow simulations are presented.

Submerged gas-injected system can be applied to various industrial field such as metallurgical and chemical processes, So this study aims at presenting the relevant relationship between gas phase and liquid phase in a gas-injected bath. In a cylinderical bath, local gas volume fraction and bubble frequency were measured by electroconductivity probe and oscilloscope. The temperature of each phase was measured using thermocouple and data acquisition system. In vertical gas injection system, gas-liquid two phase plume was formed, being symmetry to the axial direction of injection nozzle and in a shape of con. Lacal gas-liquid flow becomes irregular around the injection nozzle due to kinetic energy of gas and the flow variables show radical change at the vicinity of gas(air) injection nozzle As most of the kinetic energy of gas was transferred to liquid in this region, liquid started to circulate. In this reason, this region was defined as 'developing flow region' The Bubble was taking a form of churn flow at the vicinity of nozzle. Sometimes smaller bubbles formed by the collapse of bubbles were observed. The gas injected into liquid bath lost its kinetic energy and then was governed by the effect of buoyancy. In this region the bubbles which lost their kinetic energy move upward with relatively uniform velocity and separate. Near the gas nozzle, gas concentration was the highest. But it started to decrease as the axial distance increased, showing a Gaussian distribution.

Vapor explosion is one of the most important problems encountered in severe accident management of nuclear power plants. In spite of many efforts, a lot of questions still remain. So, KAERI launched a real experimental program called TROI using $UO_{2}$ and $ZrO_{2}$ to investigate the vapor explosion. Besides TROI tests, a small-scale experiment with molten-tin/water system was performed to quantify the characteristics of vapor explosion and to understand the phenomenology of vapor explosion. A vapor explosion was observed while the amount of air bubble and water temperature were systematically varied The mass and temperature of tin are $50\;g\;and\;150^{\circ}C$, respectively. Water temperature is set to $24^{\circ}C\;and\;50^{\circ}C$. The void fraction of air bubble ranges from $0\;to\;10\;{\%}$. The strength of vapor explosion was measured using dynamic pressure sensors attached in reactor tube wall. as a function of void fraction. In addition, a high speed video filming up to 1,000 flame/sec was taken in order to visually investigate the behavior of the vapor explosion .

The flow characteristics of ice slurry which was made from $6.5{\%}$ ethylene glycol-water solution flowing in the special pipings including the enlargement, the contraction and the orifice were experimentally investigated. The flow patterns and the pressure drops were measured in acrylic pipes when the fraction of ice were varied from $0\;to\;30{\%}$. The pressure drop behavior of the contraction and the orifice appears to be similar to that of the elbow pipe, since these piping may provide similar frictional resistance to the elbow. In the mean while, the pressure drop increased unexpectedly high with the Ice fraction in the enlargement pipe. It seems that the onset of sharp increase in the pressure drop depends on the flowing time as well as the ice fraction and the flow rate.

A number of the experiments on the phenomenon in which the thermal energy was transferred from a hot sphere to the surrounding water through the film boiling process had been conducted. As the sphere only carried the thermal energy associated with its initially high temperature but did not contain any other thermal source, the film boiling was only driven by the decreasing temperature of the sphere and, thus, was time dependent. The results from the experiments showed that the temperature of the sphere was slowly decreased in the beginning. This corresponded to the period in which the sphere was penetrating the water surface. Later, when the sphere was fully submerged and the transition film boiling was observed over the whole surface, the temperature of the sphere was decreased relatively much faster. In the last stage, the temperature of the sphere was again slowly decreased. This was considered caused by the relatively low temperature of the sphere, which reduced and later ceased the film boiling process. In addition, the estimation of the departure rate of the steam bubbles from the film layer was also correlated for the experiments.

A series of steam explosion experiments using real core materials of $ZrO_2$ and corium(a mixture of $ZrO_{2}\;and\;UO_{2}$) has been performed to evaluate the risk of steam explosion load in nuclear power plants. Surprisingly, spontaneous steam explosions are observed far both materials, which have been thought to be inexplosive so far. The dynamic pressure and morphology of the debris clearly indicate the evidence of an explosion. The experimental results also indicate that $ZrO_2$ is more explosive than corium.

The CHFG(Critical Heat Flux in Cap) test results have been evaluated to quantify the critical powers in hemispherical narrow gaps and Park's correlation has been verified. For verification of the Park's correlation, the critical powers in hemispherical narrow gaps have been measured using water at a pressure of 0.1-0.3 WPa for the hemispherical gap thickness of 1.0mm, and the heated top diameter of 238mm The measured critical powers were compared with Park's and Monde's correlation. As a results, the Park's correlation was most accurate of other correlations, however, estimated the critical powers somewhat larger than measured ones. So, it was suggested that the diameter factor in Park's correlation should be revised.

Numerical study on the advection in a microchannel with multi blocks has been performed. the microchannel is simulated three-dimensionally using a commercial fluid analysis code, FLUENT 6.0. Geometric factors of blocks are used in analysing the mixing effect. the numerical results show that the height of blocks in the channel is a key factor. Further study is required to investigate mixing effect of the microchamel with various shapes.

A three-dimensional inlet flow structure inside a microfluidic element has been investigated using a micro-PIV(particle image velocimetry) measurement as well as a numerical analysis. The present study employs a state-of-art micro-PIV system which consists of epi-fluorescence microscope, 620nm diameter fluorescent seed particles and an 8-bit megapixel CCD camera. For the numerical analysis, a commercial software CFD-ACE+(V6.6) was employed for comparison with experimental data. Fixed pressure boundary condition and a 39900 structured grid system was used for numerical analysis. Velocity vector fields with a resolution of $6.7{\times}6.7{\mu}m$ has been obtained, and the attention has been paid on the effect of varying measurement conditions of particle diameter and particle concentration on the resulting PIV results. In this study, the microfluidic elements were fabricated on plastic chips by means of MEMS processes and a subsequent melding process.

A custom micro-PIV optics assembly has been used to measure the flow field inside a T-shaped microchannel. The micro-PIV system consists of microscope objectives of various magnifications, a dichroic cube, and an 8-bit CCD camera. Fluorescent particles of diameters 620nm have been used with a Nd:YAG laser and color filters. A programmable syringe pump with Teflon tubings were used to inject particle-seeded distilled water into the channel at flow rates of $420,\;40,\;60{\mu}L/hr$. The microchannels are fabricated with PDMS with a silicon mold, then $O_2-ion$ bonded onto a slide glass. Results show differences in flow characteristics and resolution according to fluid injection rates, and magnifications, respectively. The results show PIV results with vector-to-vector distances of $2{\mu}m$ with 32 pixel-square interrogation windows at $50{\%}$ overlap.

The present parer outlines the system identification and vibration control performance of air traffic control tower of Yangyang international airport with tuned mass damper(TMD). From the free vibration test, natural frequency, damping ratio and mode shape of tower are obtained and these values are compared with the values from numerical analysis. In the vibration control test to evaluate the vibration control performance, equivalent damping ratio increased by tuned mass damper are obtained in case the TMD is operated as passive mode. Damping ratio of tower evaluated from free vibration test is about $1.0{\%}$. It is very low value than damping ratio recommended in general code. Damping ratio of passive mode is about $5{\%}$. These equivalent damping ratio increased by TMD is enough to enhance the serviceability of tower structure under wind load.

The excessive wind-induced motion of tall buildings most frequently result from vortex shedding induced across-wind oscillations. This form of excitation is most pronounced far relatively flexible, lightweight and lightly damped structure, e.g. tall building. This paper discusses aerodynamic means for mitigating the across-wind vortex shedding induced in such situations. Emphasis is on the change of the building cross section to design the building with openings from side to side which provide pressure equalization and tend to reduced the effectiveness of across-wind forces by reducing their magnitudes and disrupting their spatial correlation. Wind tunnel test have been carried out on the Kumoh National University of Technology using rigid models with twenty-four kinds of opening shapes. Form these results, the effective opening shape, size and location for building to reducing wind-induced vortex shedding and responses are pointed out.

In this study, a new passive aerodynamic control method is proposed. Control plate which is oscillated by TMD-like mechanism makes flutter stabilizing airflow. Effectiveness of proposed model is verified by experimental and analytical study. In addition, various parameters of the proposed system are investigated. Applicability to long span bridge is also examined. According to the research results, proposed model is very effective in suppressing flutter, and it also shows remarkable robustness.

In this study, a turbulence simulation is carried out in a suction type wind tunnel using grids, where turbulent flows with various turbulence intensity are successfully produced by the change of grid size, arrangement of grids and settling position, respectively. Response tests of rectangular cylinder models with aspect ratio of 2 and 4 are carried out in smooth flow and generated turbulent flows. Additionally, two types of fairing are considered such as right triangle and regular triangle. The effects of wind velocity fluctuations and fairing are discussed on vortex-induced oscillation.

Wake downstream of an object in the stratified flow has been of long-standing interest in fluid dynamics because of its similarity to geophysical flow over topographical terrains and more recently, concerns about the wake left behind a body moving through the ocean thermocline. Decades of studies of geophysical flow have unveiled that the flow downstream of obstacles in stratified flow consists of attached wake and strong internal waves, or separated, fluctuating wake and persistent late wakes, all of which depend on the flow conditions. Among unique and interesting characteristics of the stratified flow past obstacles is the generation of coherent vortex structure in the late wake far downstream of the object. Without the density stratification, the flow field downstream becomes undisturbed after relatively fast diminishing of the near wake. However, no matter how small the stratification is, the flow field downstream self-develops coherent vortex structures even after diminishing of the near wake. This paper present a computational approach to simulate the generation mechanism of the coherent vortex and analysis of the vortical structure.

The purpose of this 3-D numerical simulation is to calculate and examine a 500 kW Horizontal Axis Wind Turbine (HAWT) power performance and 3-D rotor flow characteristics, which are compared to calculation data from Delft University. The experimental approach, which has been the main method of investigation, appears to be reaching its limits, the cost increasing relate with the size of wind turbines. Hence, the use of Computational Fluid Dynamics (CFD) techniques and Navier-Stokes solvers is considered a very serious contender. We has used the CFD software package CFX-TASCflow as a modeling tool to predict the power performance and 3-D flow characteristics of a wind turbine on the basis of its geometry and operating data. The wind turbine with 40m diameters rotor, it was scaled to compare with the calculation data from delft university. The HAWT, which has eight-rpm variations are investigated respectively. The pitch angle is $+0.5^{\circ}$and wind speed is fixed at 5m/s. The tip speed ratio (TSR) of the HAWT ranging from 2.89 to 9.63.

In the present paper, the practicability of Commentray on Wind Loads of Korean Standard Design Loads fur Buildings, which has been enacted in Minstry of Construction and Transportation in June 2000, is validated by using the meteorological data measured at Homi-Cape, Pohang. Assuming that the vertical wind profile follows the power-law in a quasi-steady state, wind profile exponents are calculated by seasons by using wind data as collected through four monitoring towers. According to the Commentray on Wind Loads, Pohang is classified with its exposure category being B and the wind profile exponent being 0.22, while it is identified that the average wind profile exponent as calculated in this study is 0.26. Also, in this paper, a cross-correlation method is suggested in order to identify any meteorological correlation between measurement sites quantitatively.

Mixing zone characteristics of the wastewater effluents discharged from an ocean outfall have been investigated using an outfall mixing zone analysis method. This analysis was applied to Sokcho ocean outfall case and was used to simulate the concentration distribution of wastewater effluents. Hourly concentration distributions of outfall effluent discharges were averaged for the period of 15 days. The results show that near-field dilution was relatively high with the dilution of 130 for the winter season.

Many active and passive flow control methods have been studied since decades, but there are only few works about flow control methods using ion wind. This paper presents an experimental study on the wake control behind a circular cylinder using ion wind, a bulk motion of neutral molecules driven by locally ionized air of corona discharge. Experiments are done f3r different electrohydrodynamic numbers - the ratio of an electrical body farce to a fluid Inertial force - from 0 to 2 and for the Reynolds number ranging from $4{\times}10^3\;to\;8{\times}10^3$. Pressure distributions over a cylinder surface are measured and flow visualizations are carried out by smoke wire method. Flow visualizations confirm that ion wind affects significantly the wake structure behind a circular cylinder and pressure drag could be dramatically reduced by the superimposing ion wind.

In order to examine the possibility of using a cavity as a passive device for reduction of skin friction and heat transfer, an intensive parametric study over a broad range of the cavity depth and length at different Reynolds numbers is performed for both laminar and turbulent boundary layers in the present study. Direct and large eddy simulation techniques are used for turbulent boundary layers at low and moderate Reynolds numbers, respectively. for both laminar and turbulent boundary layers over a cavity, a flow oscillation occurs due to the shear layer instability when the cavity depth and length are sufficiently large and it plays an important role in the determination of drag and heat-transfer increase or decrease. For a cavity sufficiently small to suppress the flow oscillation, both the total drag and heat transfer are reduced. Therefore, the applicability of a cavity as a passive device for reduction of drag and heat transfer is fully confirmed in the present study. Scaling based on the wall shear rate of the incoming boundary layer is also proposed and it is found to be valid in steady flow over a cavity.

The offects of a small control rod (d=3mm) located near a main circular cylinder on the drag reduction and wake structure modification were investigated. The location of the small control rod mounted on a rod-like spring is self-adjusting according to the wake structure far optimal control of the flow around the main cylinder. The experiments were carried out at the Reynolds numbers based on the cylinder diameter (D=50mm) in the range $Re_{D}=1{\times}10^4{\~}6{\times}10^4$. Mean velocity and turbulent statistics were measured with varying the angle along the cylinder circumference ${\Theta}=15^{\circ},\;30^{\circ},\;45^{\circ}$ and the distance between the main and control rods L =0.7, 1. Compared with the bare cylinder, the main circular cylinder with the fixed and self-adjusting rods reduced drag coefficient by $10{\%}$ at the angle of ${\Theta}=45^{\circ}$. For the main cylinder with self-adjusting rot as the Reynolds numbers increase, the streamwise mean velocity is increased, however, the turbulence intensity is decreased. In addition, the control rods tested in this study are effective at higher Reynolds number than at lower Reynolds number.

A new stereoscopic PIV has been constructed using two CCD cameras, stereoscopic photogrammetry, and 3D-PTV principle. The capability of the developed stereoscopic PIV was verified through a test on the Standard Images which are provided on the website of VSJ. The arrangement of the two cameras was based on angular position For the acquisition of 3D velocity vectors, 2D velocity vectors were obtained using the gray-level cross-correlation method from the two camera images and they were matched stereoscopically. The wake of a circular cylinder was measured and turbulent properties were presented.

Maximum drag reduction (MDR) in turbulent channel flow by polymer additives is studied by direct numerical simulation. An Oldroyd-B model is adopted to express the polymer stress because it is believed that MDR is closely related to the elasticity of the polymeric liquids. The Reynolds number based on the bulk velocity and the channel height is 40000. MDR in the present study is $44{\%}$ and this is in a good agreement with the Virk's asymptote. Turbulence statistics are also in good agreements with the experimental observation. In the 'large drag reduction', the decrease of turbulent kinetic energy is compensated by the increase of energy transfer from the polymer to the flow. Therefore, MDR is a dynamic equilibrium state of the energy transfer between the polymer and the flow.

Riblets with longitudinal grooves along the streamwise direction have been used as an effective flow control technique for drag reduction. A flexible micro-riblet with v-grooves of peak-to-peak spacing of $300{\mu}m$ was made using a MEMS fabrication process of PDMS replica. The flexible micro-riblet was attached on the whole surface of a NACA0012 airfoil with which grooves are aligned with the streamwise direction. The riblet surface reduces drag coefficient about $7.9{\%}\;at\;U_o=3.3m/s$, however, it increases drag about $8{\%}\;at\;U_o=7.0m/s$, compared with the smooth airfoil without riblets. The near wake has been investigated experimentally far the cases of drag reduction ($U_o\;=\;3.3 m/s$) and drag increase ($U_o\;=\;7 m/s$). Five hundred instantaneous velocity fields were measured for each experimental condition using the cross-correlation PIV velocity field measurement technique. The instantaneous velocity fields were ensemble averaged to get spatial distribution of turbulent statistics such as turbulent kinetic energy. The experimental results were compared with those of a smooth airfoil under the same flow condition. The micro-riblet surface influences the near wake flow structure largely, especially in the region near the body surface

Injection of microbubbles and/or polymer solution has been known to be a promising method for the reduction of frictional drag of water-borne vehicles. Naval Architects have been interested in friction drag reduction technology, since the friction drag of a commercial ship can be over $70{\%}$ of total resistance. The reduction of friction drag is also important for autonomous underwater vehicles and naval submarines to improve their durability and survivability In this study two sets of experiments were carried out for the friction drag reduction of 2-D channel wall and flat plate in the circulating water channels in Chungnam National University. Preliminary results from the experiments are presented and discussed.

The real chip and similarity model were used to investigate the thermal behavior and velocity distribution of air from the heat source with the location and the amount of heat experimentally and numerically, and compared. The heat generated in the block is not cooled by convection and show the high temperature by the stagnation of heat flow. After maintaining the high temperature of block by the natural convection, the sudden drop of temperature with the air flow was shown in the channel but the decreasing rate was small with the time. The inward block was effected by infinitesimal air flow generated between block and channel and outward block was effected by the entry condition.

To assess the electrostatic interaction of surfaces at the triple contact line, the electrostatic field is analyzed by using the finite element method. The Helmholtz free energy is used as a functional which should be minimized under an equilibrium condition. The numerical results are compared with the nonlinear analytical solution for a two-dimensional charged interface and linear solution for a wedge shaped geometry, which shows fairly good agreement. The method is applied to the analysis of electrostatic influence on the contact angle on a charged substrate. The excess free energy found to increase drastically as the contact angle approaches to zero. This excess free energy Plays an opposite role to the Primary electrocapillary effect, as the contact angle gets smaller. This enables an alternative explanation for the contact-angle saturation phenomenon occurring in electrical control of surface tension and contact angle.

This article deals with the introduction of longitudinal instability of liquid rocket (pogo) and the analytical results on the frequency responses of KSR-III propulsion feeding system. Both the stiffness of bellows and the cavitation volume of venturi affect the frequency response of the feeding system. Especially, bellows has a great roll to reduce the natural frequency of the feeding system. Also, oxidizer and fuel feeding systems of the KSR-III have natural frequencies of ${\~}280Hz\;and\;{\~}90Hz$, respectively.

The present study has been carried out for understanding of flows over an evaporator in cooling system of water. Main emphasis is to decide the flow patterns in order to enhance the flow rate with low noise level. Two cases aye examined for comparison of flows; one is blowing system and the other is suction style with respect to Inn system. Present methodology for this work is PIV (Particle Image Velocimetry) techniques and Rot-wire anemometer for velocity measurements and wind tunnel for performance of the present fan. Consequently, it is found that flows passing evaporator and other components for cooling are more effective than the suction flows. Flow details with performance of fan system are also presented for proper explanation of the conclusion.

Electrorheological fluids(ERFs) show a rapid and reversible increase in viscosity by applied electric field. It is called the electrorheological effect (ER effect). The reason for ER effect is the induction of an electric dipole in each particle, leading to the formation of clusters in the direction of the field, which resist fluid flow. Generally, the behavior of ER fluids has been modeled on those of Bingham fluids. But there are some differences between Bingham fluids and ER fluids. The visualization of ER fliuds are presented and ER effects by the forming, growing and breaking of clusters are discussed. In the low shear rate area, the pressure drop is measured by a pressure sensor and the formation of ER particles is visualized by video camera. The reason for the nonlinear behavior of ER fluids at low shear rate is explained through results of visualization.

An algorithm of 3-D particle image velocimetry(3D-PIV) was developed for the measurement of 3-D velocity Held of complex flows. The measurement system consists of two or three CCD camera and one RGB image grabber. Flows size is $1500{\times}100{\times}180(mm)$, particle is Nylon12(1mm) and illuminator is Hollogen type lamp(100w). The stereo photogrammetry is adopted for the three dimensional geometrical mesurement of tracer particle. For the stereo-pair matching, the camera parameters should be decide in advance by a camera calibration. Camera parameter calculation equation is collinearity equation. In order to calculate the particle 3-D position based on the stereo photograrnrnetry, the eleven parameters of each camera should be obtained by the calibration of the camera. Epipolar line is used for stereo pair matching. The 3-D position of particle is calculated from the three camera parameters, centers of projection of the three cameras, and photographic coordinates of a particle, which is based on the collinear condition. To find velocity vector used 3-D position data of the first frame and the second frame. To extract error vector applied continuity equation. This study developed of various 3D-PIV animation technique.

A new stereoscopic PTV was developed using two CCD cameras, stereoscopic photogrammetry based on a 3D-PTV principle. The arrangement of the two cameras was based on angular position. The calibration of cameras and the pair-matching of the three-dimensional velocity vectors were based on Genetic Algorithm based 3D-PTV technique. The constructed Stereoscopic W technique was tested on the standard images of the impinged jet proposed by VSJ. The results on the turbulent properties of the jet obtained by the constructed system showed a good agreement with the original LES data.

The flow structure and heat transfer characteristics of a turbulent buoyant jet were investigated experimentally. The instantaneous temperature and velocity fields in the near field were measured using a two-frame PIV and PLIF techniques. A thin light sheet illuminated a two-dimensional cross section of the buoyant jet in which Rhodamine B was added as a fluorescent dye. The intensity variations of LIF signal from Rhodamine B molecules scattered by the laser light were captured by a CCD camera after passing an optical filter. By ensemble averaging the instantaneous temperature and velocity fields, the mean temperature and velocity fields as well as the spatial distributions of turbulent statistics were obtained. The results show the flow structure and convective heat transfer of the developing shear layer in the near field.

Stereoscopic particle image velocimetry is a measurement technique to acquire of three dimensional velocity field by two cameras. With a laser sheet illumination, the third velocity component can be deduced by out-of-plane velocity components using a stereoscopic matching method. Industrial fluid flows are almost three dimensional turbulent flows, so it is necessary to use the stereoscopic PIV measurement method. However the existing stereoscopic PIV system seems hard to use since it is very expensive and complex. In this study we have developed a Stereoscopic Miniature PIV(MPIV) system based on the concept of the Miniature PIV system which we have already developed. In this paper, we address the design and some first experimental results of the stereoscopic PIV system. The Stereoscopic MPIV system features relatively modest performances, but is considerably smaller, cheaper and easy to handle. The proposed Stereoscopic MPIV system uses two one-chip-only CMOS cameras with digital output. Only two other chips are needed, one for a buffer memory and one for an interfacing logic that controls the system. Images are transferred to a personal computer (PC) via its standard parallel port. No extra hardware is required (in particular, no frame grabber board is needed).

In the present study, kinetic energy of laminar steady flow in the exit region connected to the square-sectional $180^{\circ}$curved duct was investigated experimentally. The experimental study for air flows was conducted to measure kinetic energy distributions by using the Particle Image Velocimetry(PIV) system with the data acquisition and processing system of Cactus 2000 software. The results obtained from experimental studies are summarized as follows : (1) The critical Reynolds number for a change from laminar steady flow to transitional steadt flow was about 1910, in the 50 region of dimensionless axial position (x/Dh) whirh was considered as a fully developed flow region. (2) Maximum kinetic energy of laminar steady flow was gradually increased as the Reynolds number increased.

The objective of this study is to investigate characteristics of flow by the Rod shape and the choice of the turbulent intensity enhancement section. The Rod was setup vertically to the way of a nozzle exit flow and nozzle diameter is 17mm. Rod height is 5mm and its shapes are square, triangle, and circle. Characteristics of fluid such as velocity vector distribution, kinetic energy, turbulent intensity, and etc. were visualized, observed, and considered at 3 kinds of Re No. such as 2000, 3000, and 4000. The characteristics of flow field were investigated in each case of the distance rate from the nozzle exit to impinging plate(H/B=8, 10). The temperature of water is $20^{\circ}E$ and the measurement region divided by 3 sections(I, II, III). The nozzle diameter is 17mm. As the experimental result by PIV measurement, scale of the vector profile showed a tendency to an unbalance parabola distribution as increasing of the Re No. When the impinging plates such as square, triangle, and circle shape are installed respectively in front of the flow accelerated, rod shape of the highest velocity vector is circle shape and rod shape of the highest turbulent Intensity is square shape.

In coflow jets with relatively large size nozzle and low fuel jet velocity, the buoyancy effect arises from the density difference between fuel and air streams. The present study investigated the behavior of such a buoyant cold Jet both numerically and experimentally, especially when the fuel stream has higher density than air. It has been demonstrated that the cold jet has a circular cone shape since upwardly injected fuel jet decelerates and forms a stagnation region, when the fuel jet was composed of propane highly diluted with nitrogen. When the fuel was moderately diluted, numerical results showed the Kelvin-Helmholtz type instability along the mixing layer of the jet. The stagnation height increases nonlihearly with fuel jet velocity with the power of approximately 1.64.

A heating flow discharged from a 4-way ceiling type indoor unit has been investigated using a PIV(particle image velocimetry) system For the PIV measurements, an experimental model of 1/10 scale with a transparent room was devised by satisfying the Archimedes number, which is generally used in case that the forced convection has the similar magnitude as the natural convection. To optimize the heating flow, several vane angles and vane control algorithms of cross and right angle controls were considered. Regarding the vane angle, the experimental results show that it should be less than $30^{\circ}$ to avoid re-suction flows which decrease the performance of the air-conditioner. At the vane angle of $30^{\circ}$, applying open/close control gives nae to more uniform distribution of the heating flow than without control. Especially, the cross-control seems to be satisfactory for the thermal comfort.

Flow through turbomachinery has a very complex structure and is intrinsically unsteady. Especially, recent design trend to turbomachinery with short axial spacing makes the flow extremely complex due to the interaction between stator and rotor. Therefore, it is very necessary to clearly understand the complex flow structure to obtain the high efficiency turbomachinery. So, in this paper, the effects of axial spacing on the unsteady secondary flow performance in the one stage turbine are investigated by three-dimensional unsteady flow analysis. The three-dimensional solver is parallelized using domain decomposition and Message Passing Interface(MPI) standard to overcome the limitation of memory and the CPU time in three-dimensional unsteady calculation. A sliding mesh interface approach has been implemented to exchange flow information between blade rows.

The turbine performance test of an axial-type turbine is carried out with various axial gap distances between the stator and rotor. The turbine is operated at the low pressure and speed, and the degree of reaction is 0.373 at the mean radius. The axial-type turbine consists of ons-stage and 3-dimensional blades. The chord length of rotor is 28.2mm and mean diameter of turbine is 257.56mm. The power of turbo-blower for input power is 30kW and mass flow rate is $340m^3/min\;at\;290mmAq$ static-pressure. The RPM and output power are controlled by a dynamometer connected directly to the turbine shaft. The axial gap distances are changed from a quarter to two times of stator axial chord length, and performance curves are obtained with 7 different axial gaps. The efficiency is dropped about $5{\%}$ of its highest value due to the variation of axial gap on the same non-dimensional mass flow rate and RPM, and experimental results show that the optimum axial gap is 1.0-1.5Cx.

Turbines have been known to be particularly susceptible to flow-induced self-excited vibration. In such vibrations, direct damping and cross stiffness effects of aerodynamic forces determine rotordynamic stability. In axial turbines with eccentric shrouded rotors, the non-uniform sealing gap causes azimuthal non-uniformities in the seal gland pressure and the turbine torque which destabilize the rotor system. Previously, research efforts focused solely on either the seal flow or the unshrouded turbine passge flow. Recently, a model for flow in a turbine with a statically offset shrouded rotor has been developed and some stiffness predictions have been obtained. The model couples the seal flow to the passage flow and uses a small perturbation approach to determine nonaxiymmetric flow conditions. The model uses basic conservation laws. Input parameters include aerodynamic parameters (e.g. flow coefficient, reaction, and work coefficient); geometric parameters (e.g. sealing gap, depth of seal gland, seal pitch, annulus height); and a prescribed rotor offset. Thus, aerodynamic stiffness predictions have been obtained. However, aerodynamic damping (i.e. unsteady aerodynamic) effects caused by a whirling turbine has not yet been examined. Therefore, this paper presents a new unsteady model to predict the unsteady flow field due to a whirling shrouded rotor in turbines. From unsteady perturbations in velocity and pressure at various whirling frequencies, not only stiffness but also damping effects of aerodynamic forces can be obtained. Furthermore, relative contributions of seal gland pressure asymmetry and turbine torque asymmetry are presented.

The objective of this study is to clarify the detailed performances of the impulse type radial turbine and to present an optimum configuration of the turbine. The impulse type radial turbine has been manufactured and investigated experimentally under steady and sinusoidally oscillating flow conditions by model testing. Then, the starting characteristics under sinusoidally flow conditions have been evaluated by a numerical simulation using a quasi-steady analysis. As a result, the running and starting characteristics of the impulse type radial turbine for wave energy conversion have been clarified. Furthermore, the recommended configuration is presented, especially for setting angles of inner and outer guide vanes.

A numerical study on the use of the momentum interpolation mettled for flows with a large body force is presented. The inherent problems of the momentum interpolation method are discussed first. Numerical experiments are performed for a typical flow involving a large body force. The tact that the momentum interpolation method may result in physically unrealistic solutions is demonstrated. Numerical experiments changing the numerical grid have shown that a simple way of removing the physically unrealistic solution is a proper grid refinement where there is a large pressure gradient. An effective way of specifying the pressure and pressure correction at the boundary by a local mass conservation near the boundary is proposed, and it is shown that this method can effectively remove the inherent problem of the specification of pressure and pressure correction at the boundary when one uses the momentum interpolation method.

In this study, simulation of weak shock waves are peformed by a two-dimensional thermal fluid or compressible fluid model of the lattice Boltzmann method. The shock wave represents an abrupt change in fluids properties, in which finite variations in pressure, internal energies, and density occur over the shock thickness. The characteristics of the proposed model with a simple distribution function is verified by calculation of the sound speeds, and the shock tube problem. The reflection of a weak shock wave by wedge propagating in a channel is performed. The results agree well with those by finite difference method or by experiment. In the simulation of unsteady shock wave diffraction around a sharp corner, we show a flow field of vortical structure near the comer.

For the reduction of fuel consumption of high speed, the aerodynamic drag must be reduced. In early vehicle design process, it is very important to have information about aerodynamic characteristics of design models. In this phase CFD methods are usually used to predict the aerodynamic forces. But commercial programs using turbulence models cannot give a good agreement with experimental result and have also problems with convergence. PowerFLOW employs a new technology called DIGITAL PHYSICS, which provides a different approach to simulating fluids. DIGITAL PHYSICS uses a lattice-based approach (extended from lattice-gas and lattice-Boltzmann methods) where time, space and velocity are discrete. This discrete system represents the Wavier-Stokes continuum behavior without the numerical instability Issues of traditional CFD solvers, such as convergence. In this paper, aerodynamic performance of vehicles are simulated using PowerFLOW by Exa and results are compared with experimental wind tunnel data.

In the flow analysis around a bluffbody such as road vehicles, drag reduction has been of the primary concern mainly due to the effect on fuel economy. To reduce the drag, which is mostly due to the pressure difference caused by the flow separation, the location of the separation and eddy sizes are controlled. However, less attention has been given to the effect of the lift. The effect of lift may cause the driving stability problem of the vehicle at high speed white heavy downward effect of lift together with the vehicle weight may require more power to drive the vehicle forward. It is considered worthwhile to pursue the optimal design of the low drag tail shape of the MIRA model while taking the lift effect into account, even though it is considered as a reference. To this end, a commercial multi-objective optimization code, FRONTIER, Is used together with the CFD code, STAR-CD. It is hoped that the results will provide more insight into the flow field around the bluffbody as transportation means.

In the present study, flow simulations of a solution pump fer an absorption chiller are performed. The results are compared to the experimental data. Since the cavitation is more likely to occur in a solution pump due to Its operation under vacuum condition, and the cavitation was not considered in the present computations, the computed and experimental results show large discrepancies. For more accurate performance prediction of a solution pump, a cavitation model is required in the flow simulation. Flows through an inducer are also studied to see the effect of design parameters on performance characteristics. It is shown from the results that, if not properly designed, recirculation legion may exist near the hub region of the Inducer, and the suction surface may experience higher pressure than the pressure surface of the inducer, which may deteriorate the performance.

A flow characteristics of automotive oil pump of gerotor type has been investigated numerically. For the simulation of gerotor that have different rotating velocity at inner rotor and outer rotor, node expansion and contraction method was adopted. ASI (Arbitrary Sliding Interfaces) method was also applied at interface between rotating gerotor and stationary volute. The present results showed good agreement with the experimental data.

In this paper, the pumping performance of the disk-type drag pump which works in the outlet pressure range from 4 to 0.001 Torr is studied experimentally. The pumping characteristics of various drag pumps are performed. The inlet pressures are measured for various outlet pressures of the test pump. The flow-meter method is adopted to calculate the pumping speed. Compression ratios and pumping speeds for the nitrogen gas are measured. The present experimental data show the leak-limited value of the compression ratio in the molecular transition region. The rotational speed of the pump is 24,000rpm. The inlet pressures are measured for various outlet pressures of the test pump. The ultimate Pressures for zero throughput are measured for three-stage, two-stage and single-stage disk-type, respectively.

The periphery pump(or regenerative pump) has been generally applied in the automotive fuel pump due to their low specific speed(high heads and small flow rate) with stable performance curves. In this study, the performance prediction of side channel type periphery pumps has been developed. The prediction of the circulatory flow rate is based on the consideration of the centrifugal force field in the side-channel and in the impeller vane grooves. For the determination of performance curve(head-flow rate), momentum exchange theory is used. The effects of various geometric parameters and loss coefficients used in the performance prediction method on the head and efficiency are discussed and the results were compared with experimental data.

Pump casing has blockage effects on Impeller flow in a centrifugal pump such that the flow field around volute tongue has quite large change when the impeller rotates. A double suction pump is widely used in industrial world because it has lower NPSH required than a single suction pump. Thus, in this study, the interaction between impeller and volute casing has been investigated by using CFD for a double-suction centrifugal pump. Quasi-steady method and full pump model has been employed for the numerical calculation.

A commercial CFD code is applied to analyze the 3-D viscous flow field within vertical multi-stage centrifugal pump including impeller of centrifugal pump with 6 blades and guide vain with 11 blades. The numerical analysis of vertical multi-stage centrifugal pump is performed by changing flow rate from $8\;to\;26\;m^{3}/h$ at the constant 3500rpm. The characteristics such as total pressure coefficient, total head, water horse power, power efficiency are represented according to flow rate changing. In the future, we will need to perform flow calculation of vertical multi-stage centrifugal pump by considering meridional shape of impeller.

Flows in a regenerative pump were calculated for several flow-rates, using the CFX-Tascflow. The calculated results show the vortex structure in the impeller and side channel. The predicted performance shows considerable discrepancy form the measured values for low flow rates. Main source of the difference is the leakage flow of pump strongly affecting the performance of pump. A simple correlation was proposed using calculated leakage flows through the parametric calculations of the simplified passage.

In general, the function of intake structure, whether it be a open channel, a fully wetted tunnel, a sump or a tank, is to supply an evenly distributed flow to a pump station. An even distribution of flow, characterized by strong local flow, can result in formation of surface or submerged vortices, and with certain low values of submergence, may introduce air into pump, causing a reduction of capacity and efficiency, an increase in vibration and additional noise. Uneven flow distribution can also increase or decrease the power consumption with a change in total developed head. To avoid these sump problems pump station designers are considered intake structure dimensions, such as approaching upstream, baffle size, sump width, width of pump cell and so on. From this background, flow characteristics of intake within sump are Investigated numerically to obtain the optimal sump design data. The sump model is designed in accordance with HI code.

Numerical simulation using well-known commercial software Fine/Turbo is applied to the analysis of the aerodynamic performance for the supersonic turbine system with partial admission nozzle. Calculation was performed for coupled system of nozzle and blades using mixing plane method. In addition, calculations were also performed for the blades alone to investigate the effect of the performance variation with blade profile. These computational results are compared with the experiments. The agreement between the prediction and the experiment was found to be satisfactory..

Numerical analysis and its inverse design process of 2nd stage of JT8D aircraft engine is described. One of the most important factors that affect the performante of turbomachine is secondary flow in the blade passage, so that the performance of turbomachine can be improved by controlling secondary flow. In this paper, as a method to control secondary flow, commercial inverse design program, TurboDesign is used. Meridional derivative of angular momentum is selected as a parameter to control blade leading in this program, To validate inverse designed model, computational analysis is applied which includes rotor-stator-interaction. In this paper, CFB results of both original and inverse designed model are compared to examine how much the performance improves without reduction of work output.

As machines become smaller and faster multileaf foil bearings are used to overcome the problems with heat, friction and wear Systems with foil bearings do not need a separate system for lubrication. These bearings are self acting and are therefore green systems. Until now, there have been many studies on the structural and dynamical performances. Therefore the object of the present study is to predict the flow and structural characteristics by using the Fluid/structure interaction method. The increase in RPM led to the increase in pressure, temperature difference, maximum velocity, Mach number, shear stress and torque. In the case of 90,000 RPM effects such as choking led to a non-lineararity in the system. Also the effect of eccentricity ratio was observed and showed that eccentricity increased the maximum pressure and the density difference while decreasing the shear stress and torque.

Rotating stall in vaneless diffusers of centrifugal compressor occurs in the diffuser wall due to flow separation at large inlet flow angle. For this reason, the critical inlet flow angles are suggested by several researchers. Beyond this critical angle, flow separates in the diffuser, and develops into rotating stall. This paper studied this critical flow angle. Rotating stall is measured through eight fast-response pressure transducers which are equally spaced around the circumference at the inlet and exit of a vaneless diffuser. Experiments are done from 20000rpm to 60000rpm for the diffuser stall. Two-cell structure which rotates at $6{\~}l0{\%}$ of impeller speed is fully developed at $20000{\~}40000rpm$, and three-cell structure which rotates at $7{\~}9{\%}$ of impeller speed is fully developed at $50000{\~}60000rpm$. This paper shows that the critical inlet flow angle is not constant but related with tip speed of impeller. As tip speed increases, so does the critical inlet flow angle.

Turbo-fan for ceiling cassette type air conditioner doesn't operate in general volute. It is operated by porous material, heat exchanger. Heat exchanger increases resistance of air conditioning system and disturbs exit-flow of impeller. Therefore it has some influences on impeller capacity. In this study, we want to how that influence of exchanger on impeller capacity for ceiling cassette type air conditioner. To research, we made circular case that didn't have asymmetric part unlike rectangular case. With and without heat exchanger we measured total pressure and static pressure of impeller and three-dimensional rear flow field From the result, a turbo fan , installed in the 35mm back of fan and operated in heat exchanger, experienced $2{\%}{\~}5{\%}$% total pressure loss over all flow rate. With heat exchanger impeller efficiency decrease as flow rate decrease when flow rate coefficient was below 0.18. Especially when flow rate coefficient was below 0.12, there was $20{\%}{\~}30{\%}$ decrease of impeller efficiency.

It has been recognized that the flow in the blade passage of an axial turbomachinery rotor is very complex and is influenced by various flow phenomena, of which the tip leakage flow passing through the gap between rotor blade tip and casing plays a significant role. The losses produced due to the existence of the clearance have been known to be a large contributor of the rotor overall losses. Despite several experimental studies on non-rotating blade in the cascade configuration, and on actual rotating blades, the detailed nature of the complex flow phenomena associated with tip leakage, however, remains largely unresolved. Thus, a single-stage compressor test rig was built and measurements were taken at upstream and downstream of the rotor of this compressor at the aerodynamics laboratory of University of New South Wales. A five-hole probe and a hot-wire probe were used to measure mean and fluctuating flow parameters. The results show that tip leakage losses rise rapidly beyond tip gap of 0.01 Furthermore, the present project also identifies the regions in the wake behind the rotor of the axial compressor where such losses are concentrated. These results should be useful in the better design of rotors for improved performance of axial compressor.

The object of present study is to find the flow characteristics of centrifugal turbomachinery according to the volute shape change. The experiments were carried out for two model; volute-A designed by free vortex conditions and volute-B designed for use in high mass flow rate conditions using the volute-A test results, Flow measurements were taken in shrouded impeller with 12 backward type blades by using a five-hole pilot-tube and carried out in 4 flow rate, $Q/Q_d\;=0,43,\;1.0,\;1.27,\;1,47$, respectively, For volute-B, we found that pressure distribution was more uniform at high flow rate and from $Q/Q_d\;:\;0,43\;to\;Q/Q_d\;:\;1,20$, losses decreased and efficiency increased compare with volute-A.

In the present study a mesh generation program is developed for the flow analysis of industrial mixers. With this program one can select various design parameters such as impeller types, vessel type, number of baffles, number of blades, blade pitched angle, stages of impeller, rotational speed etc. Post processing capabilities are also imbedded in the program. Along with brief explanation of the program, examples of flow simulation for various type of mixers by using the program are also presented to show effectiveness of the program. It is expected that this program can be used to understand the effect of design parameters on the performance of particular type of the mixer, and hence to achieve the optimal design.

The study of the flow characteristics in two types of turbo-fans for a vacuum cleaner was performed in a previous study. In present study an analysis of a new modified model to reduce fan noise was performed by using CFD. The characteristics of three models calculated for various rotating speeds and flow rates are obtained and compared with available measured data. The results show that the modified model gives stable flow characteristics in operating range than the original model, while both models show similar performance characteristics at the range of high flow rate. Since in the modified model it takes much longer for an impeller blade to pass a diffuser blade than in the original model, and thus the peak pressure at BPF can be relieved, it is anticipated that the modified model gives much lower noise level with similar performance than the original one, which remains to be verified by unsteady computation and measurements.

In this study, three-dimensional imcompressible viscous flow analysis and optimization using response surface method are presented for the design of a jet fan. Steady, imcompressible, three-dimensional Reynolds averaged Wavier-Stokes equations are used as governing equations, and standard $k-{\varepsilon}$ turbulence model is chosen as a turbulence model. Governimg equations are discretized using finite volume method. Sweep angles are used as design variables for the shape optimization of the impeller in response surface method. The experimental points which are needed to construct response surface are obtained from the D-optimal design and finally the shape of impeller Is achieved from using a numerical optimization for the response surface which is obtained from CFD.

A tubular centrifugal fin is designed by using various methods of analysis and design. A preliminary design method based on empirical optimum curves for centrifugal fin is used to determine the geometric parameters for tubular centrifugal fan. And, Quasi-3D streamline curvature duct-flow analysis is used to provide the primary position of streamlines and spanwise distribution of flow angle f3r generation of blade geometry based on S1 surface. Three-dimensional CFD solution then is obtained to optimize the blade design. Constriction of flow path in the region of impeller, backward swept blade, and central cone, which are introduced to improve the design, successfully remove or suppress the vortices downstream of the impeller.

Experimental study was conducted to reveal the flow mechanism under rotating stall in an axial flow fan. For this study unsteady pressure was measured using high frequency pressure transducers mounted on the casing wall of rotor passage and total pressure fields were measured at the rotor upstream and downstream. The measured pressure signal was analyzed by Wavelet Transform and Double Phase Locked Averaging Technique. From the result of unsteady pressure field of the casing wall, one period of rotating stall was divided into three zones and the flow characteristics on each zone were described in detail. The pressure field was also analyzed in terms of the pressure distribution along pressure side and suction side of blade tip. From the result of total pressure fields at inlet and outlet of the rotor, the useful information on the characteristics of the stall cell in radial direction was provided.

A thermoacoustic refrigerator, using a resonant standing acoustic wave, has been built. It employs a loudspeaker to sustain a standing wave in a resonance tube. The stack of plates is 0.075mm-thick plastic sheet(Mylar), spirally wound to make a 40mm-diameter assembly 80mm long. and two heat exchanger are installed. The refrigerant is Helium gas. The description of thermoacoustic refrigerator are presented. The directly measured performance of the resonator. Tc/Th ( temperature ratio ) and $COP/COP_{c}$, (normalized by carnot's COP ) are plotted versus the thermal load Q applied to the cold end, for various speaker power.

Study of electric vehicle is popular with automobile company. However, battery cooling problem has delayed development of electric vehicle. Lifetime of electric vehicle's battery depends on the cooling effect for the battery tray. One model was simulated by 3-D, steady state, incompressible, k-e turbulent model simulation. It is found that flow inlet, outlet and inlet position are very important design parameters.

The present paper investigated the effect of ultrasonic vibrations on the melting process of a phase-change material (PCM). Furthermore, the present study considered constant heat-flux boundary conditions unlike many of the previous researches, which had adopted constant wall-temperature conditions. Therefore in the study, modified dimensionless numbers such as Stefan and Rayleigh were adopted to represent heat transfer results. The experimental results revealed that ultrasonic vibrations accompanied the effects like agitation, acoustic streaming, cavitation, and oscillating fluid motion, accelerating the melting process as much as 2.5 times, compared with the result of natural melting (i. e., the case without ultrasonic vibration). Such effects are believed to be a prime mechanism in the overall melting process when ultrasonic vibrations were applied. Subsequently, energy could be saved by applying the ultrasonic vibrations to the natural melting In addition, various time-wise dimensionless numbers provided a conclusive evidence of the important role of the ultrasonic vibrations on the melting phenomena of the PCM.

The present paper presents the experimental results for pressure drop in inclined tube bundles located in a rectangular duct. Measurements are made for pressure drop in triangular and rotated triangular tube arrays having P/d ratio of 1.6 and inclination angles of 30,45,60 and 90 degrees. The Reynolds number based on the free stream velocity and tube diameter ranges from $8{\times}10^2\;to\;6.3{\times}10^{4}$. The experimental results show that the magnitude of dimensionless pressure drop decreases significantly when the inclined angle is less than 45 degree. The measured data are compared with two existing correlations available in the literatures. The ESDU correlation agrees well with the present data far the triangular arrays. But some discrepancies are observed for the rotated triangular arrays when the inclined angles are 30 and 45 degrees. The Idel'chik correlation generally agrees well with the measured data for the rotated triangular arrays except for the inclined angle of 30 degree. The Idel'chik correlation needs modification for the triangular arrays. The modified Idel'chik correlation agrees well with the measured data within $10{\%}$. It is found that the present measured data can be applied to the evaluation and modification of previous correlations.

Hard disk drived (HDD) in computer are used extensively as data storage capacity. The trend in the computer industry to produce smaller disk drives rotating at higher speeds requires an improved understanding of fluid motion in the space between disks. The distribution of pressure disturbance on disks has relation to flow structure. To investigate the flow structure, time-resolved hot-wire measurements of the circumferential velocity component were obtained for the flow between the center pair of four disks of common radius $R_2$ coretating at angular velocity ${\Omega}$ in a fixed cylindrical enclosure. Hot-wire supporter acts as an obstruction in this case. The effects of rotating speed and size of hot-wire supporter diameter between disks on the flow driven by disks were investigated. Velocity spectra at the fixed space were measured to obtain the structure of inner and outer region in flow field.

Vehicle's airtight integrity is a very important factor that greatly affects passenger's habitability. However, when a door is closed, the pressure in the passenger compartment increases due to the vehicle's airtight integrity. That pressurizes the eardrums of the passenger, and makes passenger unpleasant. Thus, in this study, the configurations of air ventilation hasve been investigated to reduce pressure in the passenger compartment. Truck cab is utilized to measure the pressure in the passenger compartment. Various kinds of air ventilations are considered to find out optimized pressure in truck cab when a door is closed.

Pressure compensating temperature control valve(TCV) is one of the important control devices, which is used to maintain the constant temperature of working fluid in power and chemical plants. The ratio of cylinder hole diameters of inlet and outlet is the main design parameters of TCV. So this needs to be investigated to improve the function of control of temperature and void fraction. In this study, numerical analysis is carried out with various ratios of cylinder hole diameters of the inlet and outlet in the TCV. Especial1y, the distribution of the static pressure Is investigated to calculate the new coefficient($C_{\upsilon}$) and resistance coefficient(K). The governing equations are derived from making using of three-dimensional Naver-Stokes equations with standard $k-{\varepsilon}$ turbulence model and SIMPLE algorithm. Using a commercial code, PHOENICS, pressure and flow fields in TCV are calculated with different inlet and outlet diameters of the cylinder hole for cold and hot water passages.

In order to develop the water hammering cleaner system for removing scale and slime in inner metal or non-metal piping wall, the flow characteristics are investigated by numerical and experimental methods. The air bubbles in the piping systems as a shock wave are formed and transferred with the water flow in the piping. The governing equations are derived from making using of three-dimensional Wavier-Stokes equations with the standard $k-{\varepsilon}$ turbulence model and SIMPLE algorithm. Pressure distributions in the pipeline are calculated for different air supply pressures. Also, we prepared some experimental results of the pressure differences for various air supply times.

Energy dissipations in a general PHE flow are the compounded effects of the piled corrugate geometries and its wall pressure and temperature distributions. In addition, although the exchangers are substantial pieces of engineering equipment, they are composed of a very large number of nominally identical and small geometrical elements. In the present numerical study, the three-dimensionally complicated energy dissipation fields and those wall-shape-induced flow destabilization are investigated in the cross-corrugated passages, which result in high energy transports with comparatively low pressure drop. We revealed the critical conditions as $Re=157.3 for the wall-shape-induced flow destabilization in a general PHE element by initial value method, or shooting method, and compare its value to that of analytical solution of plane Poiseille flow, two-dimensional grooved flow and so on. We also observed the detailed variation of flow field and energy transportation with changes in time and flow variables such as Reynolds number. Lastly, we considered the flow natural frequency, or Strouhal number, with variation of hydrodynamic conditions for the best use of active control, such as forced mass flow rate pulsative flow, to enhance energy transportation.

The hydrodynamic instability of the three-dimensional boundary-layer over a rotating disk has been numerically investigated for three cases flows using linear stability theory (i.e. Rossby number, Ro = -1, 0, and 1). Detailed numerical values of the disturbance wave number, wave frequency, azimuth angle, radius (Reynolds number, Re) and other characteristics have been calculated for $K{\acute{a}}rm{\acute{a}}n$, Ekman and $B{\"{o}}ewadt$ boundary-layer flows. Neutral curves for these flows are presented. Presented are the neutral stability results concerning the two instability modes (Type I and Type II) by using a two-point boundary value problem code COLUEW that was based upon the adaptive orthogonal collocation method using B-spline. The prediction from the present results on both instability modes among the three cases agrees with the previously known numerical and experimental data well.

Multimode boundary-layer transition on a NACA0012 airfoil is experimentally investigated under periodically passing wakes and the moderate level of free-stream turbulence. The periodic wakes are generated by rotating circular cylinders clockwise or counterclockwise around the airfoil. The free-stream turbulence is produced by a grid upstream of the rotating cylinder, and its intensity(Tu) at the leading edge of the airfoil is $0.5\;or\;3.5\;{\%}$. The Reynolds number ($Re_c$) based on chord length (C) of the alrfoil is $2.0{\times}10^5$, and Strouhal number ($St_c$) of the passing wake is about 0.7. Time- and phase-averaged streamwise mean velocities and turbulence fluctuations are measured with a single hot-wire probe, and especially, the corresponding wall skin friction is evaluated using a computational Preston tube method. The wake-passing orientation changes pressure distribution on the airfoil in a different manner irrespective of the free-stream turbulence. Regardless of free-stream turbulence level, turbulent patches for the receding wakes propagate more rapidly than those for the approaching wake because adverse pressure gradient becomes larger. The patch under the high free-stream turbulence ($Tu=3.5{\%}$) grows more greatly in laminar-like regions compared with that under the low background turbulence ($Tu=0.5{\%}$) in laminar regions. The former, however, does not greatly change the original turbulence level in the very near-wall region while the latter does it. At further downstream, the former interacts vigorously with high environmental turbulence inside the pre-existing transitional boundary layer and gradually lose his identification, whereas the latter keep growing in the laminar boundary layer. The calmed region is more clearly observed under the lower free-stream turbulence level and for the receding wakes. The calmed region delays the breakdown further downstream and stabilizes more the boundary layer.

The most unstable situation of laminar plane Poiseuille flow for transition to turbulence is investigated by using a pseudo-spectral method. A number of various disturbance modes are tested and it is found that the flow is the most unstable when it is disturbed by an oblique wave with an angle of $29.7^{\circ}$.

Three-dimensional blood flow in the sac of the KTAH(Korean total artificial heart) is simulated using fluid-structure interaction model. The aim of this study is to delineate the three-dimensional unsteady-blood flow in the sac of KTAH. Incompressible viscous flow is assumed for blood using the assumption of Newtonian fluid. The numerical method employed in this study is the finite element software called ADINA. Fluid-structure interaction model between blood and sac is utilized to represent the deformation of the sac by the rigid moving actuator. Three-dimensional geometry of cactus type KTAH is chosen for numerical model with prescribed pressure boundary condition on the sac surface. Blood flow is generated by the motion of moving actuator and strongly interacts with the solid material surrounding blood. High shear stress is observed mainly near the inlet and outlet of the sac.

The dispersion of an aerosol bolus in acinus is analyzed numerically. Model geometry is a straight duct surrounded by an axisymmetric semicircular annulus which is expanding or contracting with breathing. Unsteady Wavier-Stokes equation is solved by CFX-F3D, an FVM commercial code and the trajectory of massless particle Is computed by Lagrangian method. For steady flow with no wall motion, mean velocity of aerosol bolus in alveolated duct is a little smaller than that in straight duct and dispersion in alveolated duct is comparable with the dispersion in straight tube. For expanding duct mean velocity of aerosol bolus approaches half of that in straight tube and effective diffusivity is smaller than that of straight tube. For contracting duct mean velocity of aerosol bolus becomes slightly larger than that in straight tube and effective diffusivity is comparable with the case of straight tube.

Shunt valves implanted in the subcutaneous tissue of brain to treat patient with hydrocephalus were numerically simulated to investigate influence of pressure pulsation on their flow control characteristics. Shunt valves are subjected to pressure variation since ventricles enclosing the brain are under pressure pulsation rather than uniform pressure due to blood pressure variation. We modeled flow orifice through shunt valve and imposed pulsating pressure and valve diaphragm movement to compute flow through the valve. The results of our study indicated that flow rate increased by $40{\%}$ by introducing pressure pulsation and diaphragm movement on the shunt valve. Our results demonstrate the pressure-flow control characteristics of shunt valves unplanted above human brain may be quite different from the characteristics obtained by syringe pump test with uniform pressure and no diaphragm movement.

The two major problems related to the blood flow in a floating type polymer valve are thrombus formation and hemolysis. It is well known that the shear stress in the fluid and flow separation around the valve are blamed for such disastrous phenomena. In this viewpoint, through study of the flow field around the valve is imperative to improve design of the valve. The aim of this study is to investigate the fluid flow around a floating type polymer valve. The numerical method employed in this study is the finite element software called ADINA. Incompressible viscous flow is assumed for blood using the assumption of Newtonian fluid. In this study, two prominent features of the axisymmetric flow around the floating type polymer valve are observed: jet-like flows observed near the gap between the conduit and the valve, and recirculating flow downstream of the valve. We also provided a detailed description of shear stress field according to the variation of flow conditions. The shear stress in fluid has its maximum value near the gap between the valve and the conduit.

Hemodynamics of the TPLS(Twin Pulse Life Support System) is numerically investigated to delineate the possibility of hemolysis in blood. Computational method employing finite element algorithm is utilized to solve the blood flow of the sac squeezed by moving actuator. We assume that the blood flow interacts with the sac material which is activated by the rigid body motion of the actuator. Valve dynamics at the ends of the sac is simplified as on/off type motion. We compute the transient viscous flow in the two-dimensional geometry of the blood sac. Incompressible laminar flow is simulated on the assumption of Newtonian fluid. Blood velocity has a step gradient near the throat of the sac formed by the moving actuator. According to the decrease of the gap size of blood passage, the magnitude of shear stress in the blood is dramatically increased. Numerical solutions show that the maximum value of shear stress in the blood flow in TPLS is relatively smaller than that of the roller type ECLS.

Catheters are used to measure translesional pressure gradients in the stenosed coronary arteries. Catheter insertions during coronary angioplasty cause flow obstructions. A narrowed flow cross section with catheter present effectively introduced a tighter stenosis than the enlarged residual stenoses after balloon angiplasty. In general, the form of stenoses are no uniform. Sometimes, these are occurred the irregularly eccentricity If the analyses are conducted for uniform stenosed artery, the results will be underestimated. Thus, in this study, three dimensional computer simulations are conducted to investigate the flow blockage effects due to the catheter insertion during the coronary angioplasty. The results for the eccentric stenosed artery are compared with those of the concentric stenosed artery.

Cancer angiogenesis is simulated using a two dimensional model. Governing equation of angiogenesis is a TAE (Tumor angiogenesis factor) conservation equation in time and space. A stochastic process model is utilized to simulate vessel formation, proliferation, and migration to a cancer pellet. Numerical results are presented especially in case of growing cancer.

A numerical study of a hemodynamical model for the tumor angiogenesis is carried out. The tumor angiogenesis process is comprised of a sequence of events; secretion of tumor angiogenesis factor(TAF) from the solid tumor, degradation of the basement membrane of nearby blood vessels, migration and proliferation of the endothelial cells. The model takes into account the effect of TAF concentration and endothelial cell density, and their conservation equations are represented as a set of one-dimensional initial boundary value problems. These equations are discretized by using a finite difference method in which the second order schemes both in time and in space are used. The effects of the parameters contained in the model are Investigated extensively through the numerical simulation of the discretized model. The result for the typical case compares very well with the known result.

Aging effect on the cardiovascular circulation is simulated by lumped parameter model. Aging phenomena can be hemodynamically explained as (1) the increase of flow resistance induced by remodeling of artery vessels and increased viscosity of blood and (2) the reduction of the vessel capacitance caused by arteriosclerosis. Appropriate physiological parameters are evaluated from the clinical data of adults and old men. Simulation results well explain the hypertension with aging of cardiovascular system.

The present study is to evaluate the performances of flow velocity and wall shear stress in the stenosed coronary artery using human in vivo hemodynamic Parameters and computer simulation. Initial and follow-up coronary angiographics in the patients with angulated coronary stenosis are performed. Follow-up coronary angiogram demonstrated significant difference in the percent of diameter in the stenosed coronary between two groups ($Group\;1:\;40.3{\%},\;Group\;2:\;25.5{\%}$). Flow-velocity wave obtained from in vivo intracoronary Doppler ultrasound data is used for the boundary condition for the computer simulation. Spatial and temporal variations of flow velocity vector and recirculation area are drawn throughout the selected segment of coronary models. The WSS of pre- and post-intracoronary stenting are calculated from three-dimensional computer simulation. Then negative shear stresses area on 3D simulation we noted on the inner wall of the post-stenotic area before stenting. The negative WSS is disappeared after stenting. High spatial and temporal WSS before stenting fell into within physiologic WSS after stenting. This finding was prominent in Model 2. The present study suggest that hemodynamic forces exerted by pulsatile coronary circulation termed WSS might affect on the evolution of atherosclerosis within the angulated vascular curvature. The local recirculation area which has low or negative WSS, might lead to progression of atherosclerosis.

The aim of this work is to analyze changes in cerebral hemodynamics and intracranial pressure mediated by cerebral blood flow challenges in patients with acute heart arrest. Lumped parameter model with feedback mechanism is utilized to simulate the hemodynamics of brain blood flow in case 40 min T-PLS operation is applied to patients of cardiac arrest. Numerical solutions show that cerebral blood flow and perfusion pressure in patients of cardiac arrest are sharply recovered in the initial state of T-PLS operation.

Computational modeling is presented to investigate the hemodynamics of coronary circulation for the cardiac arrest model during the operation of an ECLS called T-PLS. A lumped parameter model is utilized to delineate the hemodynamics of coronary circulation. For verification of the present method, we compute the coronary circulation for normal case and compare those results with the existing data. Numerical results on the cardiac arrest model show that T-PLS can increase LAD flow significantly.

A separate type thermosyphon can be utilized as a cooling device of electronic equipments (such as CPU of a personal computer or notebook). This study was carried out to investigate the cooling effect of separate type thermosyphon and to find the adequate parameters affecting the separate type thermosyphon. The heat transfer characteristics of separate type thermosyphon were obtained from experimental results. A $50{\times}50{\times}2 mm$ heat source was copied after CPU for the experiments. The results indicate that the device is capable of dissipating 60W of thermal energy and keeping the heat plate surface temperature under 50'E and the device can transfer heat from the evaporator to the condenser through natural circulation (without any external driving forces). Some transport phenomena of the working fluid and the heat transfer characteristics of the loop were observed in the experiments and are discussed in detail below.

This paper reports an experimental study around a module about forced air flow by blower($35{\times}35{\times}6mm^3$) in portable PC(10mm high, 200mm wide, and 235mm long). The channel inlet flow velocity has been varied between 0.26, 0.52 and 0.78m/s. The power input to the module is 4Wthis report, particular attention is directed to the fluid flow and adiabatic wall temperature($T_(ad)$) around a module which is under fluid mechanical and thermal influences of the module. The fluid flow around a module was visualized using PIV system. Liquid crystal thernography is used to determine the adiabatic wall temperature around a heated module on an acrylic board. Plots of $T_(ad)$ (or F) show marked effects of dispersion of thermal wake near the module.

The inside configuration of intake nozzle of vacuum cleaner greatly affects the dust-collection efficiency and acoustic-noise effect generated from flow separation Interaction between high-speed flow and internal structure. In order to improve the performance of the vacuum cleaner, flow fields inside the intake nozzles were investigated using flow visualization and FIV (Particle Image Vetocimetry) technique. The measurement to aerodynamic power, suction efficiency and noise level were also carried out. Valuable information was obtained from the experiments, revealing how to modify the intake nozzle. In this paper, the results of visualization, velocity distribution of flow fields, aerodynamic power, suction efficiency and noise level are discussed.

To develop and improve a GIS(Gas Insulated Switchgear), the prediction of the pressure in puffer cylinder and the flow between the nozzle and the moving electrode within GIS is very important.The leading companies in GIS business issue the results of the study of flow within GIS including arc plasma. In this study, the characteristics of the flow of the GIS developed by HHI(Hyundai Heavy Industries Co. Ltd.) was investigated. To simulate the compressible flow of GIS, the CFX, a commercial CFD code, was used. With moving grid method, the movement of piston and electrode was simulated. The moving grid method was superior to the method of varying the property of cells to move an obsticle, in stability and convergencce of solution. The calculated maximum pressure within the puffer cylinder was matched with experimental data within $5{\%}$ error. The oscilation of pressure in GIS after the movement of electrode was well predicted.

The supersonic, dual, coaxial jet impinging upon a vertical flat plate has recently been applied to a variety of industrial manufacturing processes, since it has several advantages over a conventional supersonic impinging jet. In the present study, experimentation is carried out to investigate the effects of the impinging angle of the annular flow and the design Mach number on the flow field formed over the vertical flat plate. A convergent-divergent nozzle is used to obtain the inner jet flow, its design Mach number being changed between $1.0\;and\;2.0$. The outer annular nozzle has a constant area of the Mach number of 1.0, and its impinging angle of $0^{\circ}\;and\;20^{\circ}$. The primary jet pressure ratio is changed in the range from 6.0 to 10.0 and for the annular flow, the assistant jet pressure ratio is changed from 1.0 to 4.0. The distance between the dual, coaxial nozzle and flat plate is also changed. Detailed pressure measurements are conducted along the axis of the jet and on the flat plate as well. The impinging coaxial Jet flows are visualized using the Schlieren and Shadow optical methods. The results show that the flow field on the plate is not strongly dependent only on the primary and assistant pressure ratios but also the impinging angle of the annular nozzle.

Experimental investigation of the flow field of supersonic cavity is described. In this research, supersonic cavity is used in chemical laser system. For efficient laser, downstream flow after cavity need to be uniform and clear for pressure recovery system. In previous research, it's known that there's oscillation In cavity and is due to Mach number and L/D ratio. A strong recompression occurs at the after wall and the flow is visibly unsteady. Cavity flow in this research is of the open type, that is, length-to-depth ratio $L/D<10\;at\;M\;=\;3$. Experiment is done with pressure measurement by piezo-type sensor and visualization by Schlirern method. The time-dependent experimental result is compared with computation.

An attempt to reduce supersonic Jet noise is carried out by using two steady microjets in a round jet. The jet is issued from a round sonic nozzle with an exit diameter of 10mm. Two micro-nozzles with an inside diameter of 1mm each are installed on the exit plane with an off-axis angle of $45^{\circ}$. Far-field noise was measured at a location 40 diameters off the jet axis. The angles between a microphone and the jet axis are $45^{\circ}\;and\;90^{\circ}$. For an injection rate less than $1{\%}$ of the main jet, screech tones were completely suppressed by the microjets. The reduction in the ovelall sound pressure levels were $2.4\;and\;2.7\;dB\;for\;90^{\circ}\;and\;45^{\circ}$ directions, respectively. The enhancement of mixing/spreading of the jet by the microjet was negligible. The reduction of noise is probably due to distorted shock cell structures and/or broken large scale vortical structures by the microjets.

The SFU(Solid Fuel Ram-Jet) propulsion is attractive for projectiles because of the combination of high propulsive performance and low system complexity more than conventional projectiles. The Objective of this research was to characterize the inlet aerodynamic characteristics (centerbody & pilot type) in SFRJ. Diffuser static pressure & combustion chamber pressure was tested and the AoA was changed $0^{\circ}\;and\;4^{\circ}$ at Mach number of 3.0 for performance estimate. The performance study of inlet was carried out with the Schlieren system and Supersonic cold-flow system. A Computational fluid dynamic solution is applied internal flow of inlet and the solutions are compared with experimental results.

The physical phenomena of rectangular open cavity are numerically investigated in this paper Two-dimensional cavity problems with laminar boundary layers in upstream are simulated by using the compressible Wavier-Stokes equations. The high-order and high-resolution numerical schemes are used for the evaluation of spatial derivatives and the time integration. Cross-correlation is used to analyze the characteristics of wave propagation along time and spatial. Sudden phase shifting of 90 degrees is appeared near downstream edge, and this is coincident with the phase lag suggested in original Rossiter's equation. The results give a further understanding of the physical phenomenon of noise generation, and the resonance of flow and acoustic in cavity. Moreover, modified Rossiter's equation, which is more accurate and can be applied in various conditions, is suggested. The distance from the point of vortex generation to the point of vortex collapsing acts as effective distance of cavity resonance, and the phase difference between the point of vortex collapsing and the point of acoustic source acts as phase lag. The mechanism of acoustic generation is fully understood in this paper. The mechanism of acoustic generation is fully understood in this paper.

Automobile aeroacoustics Is a developing area of technology where experimental and theoretical tools are being continuously refined to understand, analyze and modify the noise-generating mechanisms in the vehicle flow. Main sources of ground vehicle exterior noise are the tires (tire/road interaction) and the unsteady flow field around the vehicle. In this study, the sound source localization of a rolling tire was applied to the measurement of radiated sound by using an acoustic mirror system. A possible flow pattern that develops is suggested based on detailed wind tunnel investigations with a rotating wheel in contact with a moving belt.

A BEM is highly efficient method in the sense of economic computation. However, boundary integration is not easy for the complex and moving surface e.g. in a rotating blade. Thus, Kirchhoff surface is designed in an effort to overcome the difficulty resulting from complex boundary conditions. A Kirchhoff surface is a fictitious surface which envelopes acoustic sources of main concern. Acoustic sources may be distributed on each Kirchhoff surface element depending on its acoustic characteristics. In this study, an axial fan is assumed to have loading noise as a dominant source. Dipole sources can be computed based on the FW-H equation. Acoustic field is then computed by changing Kirchhoff surface on which near-field is implemented, to analyze the effect of Kirchhoff surface on it.

The performance and noise characteristics of the design parameters of a cross-flow fan are investigated by computational methods. The incompressible Wavier-Stokes equations in moving coordinates are time-accurately solved for obtaining the pressure fluctuations due to the aerodynamic interactions between the impeller blades and the stabilizer, and sound pressure is then computed by the Ffowcs Williams-Hawkings equation. Design parameters of the cross-flow fan include blade setting angle, exit-diffusion angle, and stabilizer installation angle. Also, an optimization of the aforementioned design parameters has been peformed using the Taguchi method.

Highly sweep leading edge extensions(LEX) applied to delta wings have greatly improved the subsonic maneuverability of contemporary fighters. In this study, systematic approach by PIV experimental method within a circulating water channel was adopted to study the fundamental characteristics of induced vortex generation, development and its breakdown appearing on a delta wing model with or without LEX in terms of four angles of attack($15^{\circ},\;20^{\circ},\;25^{\circ},\;30^{\circ}$) and six measuring sections of chord length($30{\%},\;40{\%},\;50{\%},\;60{\%},\;70{\%},\;80{\%}$). Sideslip effect in case of the LEX was also studied for two sideslip(yaw) angles($5^{\circ},\;10^{\circ}$) at one angle of attack(20). Distribution of time-averaged velocity vectors and vorticity over the delta wing model were compared along the chord length direction. Quantitative comparison of the maximum vorticity featuring the induced pressure distribution were also conducted to clarify the significance of the LEX existence. Animation presentation in velocity distribution was also implemented to reveal the effect of LEX with wing vortex interaction.

In order to analyse the mechanism of a flow tone around a cavity, the correlations between the flow in the cavity and the boundary layer flow in front of the cavity are studied experimentally in this paper. The instability In the boundary layer forms the vortex at the front edge of the cavity and the flow tone is occurred by the vortex breakdown at the rear edge of the cavity Therefore, the boundary layer measurement is important in the cavity flow control. We measure the velocity of the boundary layer at the entrance of the cavity using hot-wire anemometry and the flow tone around the cavity by microphone. The boundary layer characteristic is changed by the various angle of the flap on the front edge of the cavity, while it is less influenced by the ratio of length and depth of the cavity.

We examine the problem in which porous/viscoelastic compliant thin plates are subject to pressure fluctuations under transitional or turbulent boundary layer. Measurements are presented of the frequency spectra of the near-field pressure and radiated sound by compliant surface. A porous plate consisting of 5mm thick, open-cell foam with fabric covering and a viscoelastic painted plate of 1mm thick over an acoustic board of 4m thick were placed over a rigid surface in an anechoic wind tunnel. Streamwise velocity and wall pressure measurements were shown to highly attenuate the convective wall pressure energy when the convective wavenumber ($k_{ch}$) was 3.0 or more. The sound source localization on the compliant walls is applied to the measurement of radiated sound by using an acoustic mirror system.

The unsteady turbulent flow characteristics of NACA4406 airfoil at close proximity to the pound are numerically investigated, especially focused on the noise generation mechanism near the blunt trailing edge. The unsteady two-dimensional compressible Wavier-Stokes equations with a Spalart-Allmaras turbulence closure model are solved by the 6th-order compact scheme and the 4th-order Runge-Kutta scheme. The computation shows a noise generation by a feedback mechanism at the blunt tailing edge, where the acoustic-fluidic coupling occurs between the wall-reflected sound waves and the periodically disturbed turbulent shear layer.

In this work, numerical experiments ave conducted to find out the optimal shape of flapping-airfoil using thickness variation airfoils. In the previous study of flapping-airfoil, we had found that the thrust efficiency of thicker airfoil is better than thinner one, but the latter has higher thrust coefficient. Therefore, we have combined thin(NACA0009) and thick(NACA0015)airfoil to overcome these demerits of each airfoil. Using this combined airfoil, we can achieve acceptable aerodynamic performances from thrust efficiency and coefficient points of view. In order to computational study, we have used parallel-implemented incompressible Wavier-Stokes solver. Computational results show how to design leading and trailing edge shapes.

Cavity is inevitably included in automobile vehicle configuration. The complex unsteady flow and sound waves generated by the cavity are very important issues because of the involved fluid dynamics and the practical importance in the field of aerodynamics. The LES method used is a conventional one with Smagorinsky eddy-viscosity model and the computational grid is small enough to be handled by workstation-level computers. LES can successfully simulate of cavity noise analysis.

This paper dispicts the vortical flow characteristics over a delta wing using a computational analysis for the purpose of investigating and visualizing the effect of the angle of attack and fee stream velocity on the low-speed delta wing aerodynamics. Computations are applied to the full, 3-dimensional, compressible, Navier-Stokes Equations. In computations, the free stream velocity is changed between 20m/s and 60m/s and the angle of attack of the delta wing is changed between $16^{\circ}\;and\;28^{\circ}$. For the correct prediction of the major features associated with the delta wing vortex flows, various turbulence models are tested. The standard $k-{\varepsilon}$ turbulence model predict well the vertical flows over the delta wing. Computational results are compared with the previous experimental ones. It is found that the present CFD results predict the vortical flow characteristics over the delta wing, and with an increase in the free steam velocity, the leading edge vortex moves outboard and its streangth is increased.

Korea Aerospace Research Institute has been developing a 50-meter class airship to demonstrate the technology to be used in the development of a stratospheric airship, and importance of accurate prediction of the dynamic behavior of an airship before flight test is widely conceived. The added mass has large impact on the dynamic characteristics of an airship unlike for an airplane and the added mass of the airship with empennage is predicted in this paper. At first, the usability of the strip theory is examined which integrates the analytic two dimensional results in the cross section along the longitudinal axis. A panel method with the surface distribution of sources is developed and its validity is also examined. Finally, the panel method with both source and doublet distributions is implemented, and it is validated and applied for the calculation of the added mass of a 50-meter class airship. Using the methods developed, the influence of empennage and control surface deflection on the added mass property of the airship is studied.

The dispersion of Lagrangian fluid particles in a turbulent channel flow is studied by a direct numerical simulation. Four points Hermite interpolation in the homogeneous direction and Chebyshev polynomials in the inhomogeneous direction is adopted by assesing the acceleration of fluid particles. In order to characterize the inhomogeneous Lagrangian statistics, accurate single particle Lagrangian statistics are obtained along the wall normal direction. Integral time scales of Lagrangian velocity can be normalized by Eulerian mean shear stresses.

The suitability of high-order accurate, central and upwind-biased compact difference schemes is evaluated for the large-eddy simulations of flows in complex geometry. Two flow geometries are considered: channel and circular cylinder. The effects of numerical dissipation and aliasing error on the evaluation of subgrid scale stress are investigated by extending the analysis by Ghosal (1) to centered and upwind compact schemes. It is shown that the failure of upwind schemes mainly comes from the aliasing error.

Using a mathematical theory, we show that the optimality condition of a turbulent diffuser with maximum pressure recovery at the exit is zero shear stress along the wall. The optimal diffuser shape is designed through iterative procedures by using the $k-{\varepsilon}-{\nu}^{2}-f$ turbulence model for flow simulation. The Reynolds number based on the bulk mean velocity and the channel height at the diffuser entrance is 18,000. We also perform large eddy simulation to validate the shape design results and investigate the flow characteristics near the zero-skin friction wall. Results from large eddy simulation show that the skin friction is slightly higher than zero but is still very small as compared to that of the flat plate boundary layer flow Although the time-averaged wall shear stress is slightly above zero along the diffuser wall, instantaneous flow reversals occur intermittently. The streamwise mein velocity shows an asymptotic behavior of the half-power-law near the wall where the skin friction is close to zero.

The elliptic relaxation model(ERM) with the inhomogeneous correction intermediate between near wall with and far from the wall. The source of the ERM usually was appled quasi-homogeneous pressure-strain correlation in homogeneous situations. This formulation was easily applied to the linear model or non-linear pressure-strain model. It is observed that the boundary conditions of the relaxation operator dominate the homogeneous pressure-strain model in the near wall region. While looking at high-Reynolds number flows, it was found necessary to modify the effect of the relaxation operator throughout the log region by accounting for gradients of the flatness variable and turbulent length scales. These effects are kinematic blocking of the wall normal velocity fluctuation and pressure reflections from the surface. This model is wall distances and unit vectors which make the model applicable to flows boundary by a complex geometry. Inhomogeneous correction model is computed inertial and non-inertial channel flow These are compared DNS(Kim et at., Kristofffrsen & Andersson) for channel flow. The present model could be predicted well for rotating flows.

It is necessary for the numerical simulation of 3-dimensional incompressible isotropic decaying turbulence to construct 3-dimensional initial velocity field which resembles the fully developed turbulence. Although the previous velocity field generation method proposed by Rogallo(1981) satisfies continuity equation and 3-dimensional energy spectrum, it has limitation, as indicated in his paper, that it does not produce the higher velocity moments(e. g. velocity derivative skewness) characteristic of real turbulence. In this study, a new velocity field generation method which is able to control velocity derivative skewness of initial velocity field is proposed. Brief descriptions of the new method and a few parameters which is used to control velocity derivative skewness are given. A large eddy simulation(LES) of isotropic decaying turbulence using dynamic subgrid-scale model is carried out to evaluate the performance of the initial velocity field generated by the new method. It was shown that the resolved turbulent kinetic energy decay curve and the resolved enstrophy decay curve from the initial field of new method were more realistic than those from the initial field of Rogallo's method. It was found that the dynamic model coefficient from the former was initially half the stationary value and experienced relatively short transition period, though that from the latter was initially zero and experienced relatively longer transition period.

Experiments were conducted to show the characteristics of the flow on unventilated parallel plane jets. Measurements of mean velocity components and turbulent intensities were carried out with a particle image velocimetry. The measurements range of these experiments was Reynolds number of 5300 based on the nozzle width and the cases of nozzle-to-nozzle distance were 4, 6, 8, and 10 times the width of the nozzle. Results show that a recirculation zone with a sub-atmospheric static pressure was bounded by the inner shear layers of the individual jets and the nozzles plate. The positions, where maximum value of mean turbulent intensities and mean turbulent kinetic energy show, were at the same position with the merging point. The spread of jets in the merging region increases more rapidly than that of jets in the converging and the combined region. As nozzle-to-nozzle distances were increased, it was shown that merging and combined lengths were shorter.

A direct numerical simulation is performed for turbulent heat transfer in a concentric annulus at $Re_{Dh}=8900\;and\;Pr=0.71$ for two radius ratios ($R_{1}/R_{2}=0.1\;and\;0.5$) and $q^{\ast}=1.0$. Main emphasis is placed on the transverse curvature effect on near-wall turbulent thermal structures. Near-wall turbulent structures close to the inner and outer walls are scrutinized by computing the lower-order statistics. The fluctuating temperature variance and turbulent heat flux budgets are illustrated to confirm the results of the lower-order statistics. The present numerical results show that the turbulent structures near the outer wall are more activated than those near the inner wall, which may be attributed to the different vortex regeneration processes between the inner and outer walls.

An analytic approach is attempted to predict the amplification of turbulence in compressible flows experiencing one-dimensional and axisymmetric bulk dilatation. The variations of vortex radius and vorticity are calculated, and then the amplification of turbulence is obtained from them by tracking three representative vortices. For a one-dimensionally compressed flow, the present analysis slightly underestimates the amplification of velocity fluctuations and turbulent kinetic energy, relative to that of rapid distortion theory in the solenoidal limit. For an axisymmetrically distorted flow, the amplification of velocity fluctuations and turbulent kinetic energy depend not only on the density ratio but also on the ratio of streamwise mean velocities, which represents streamwise vortex contraction/stretching. In all flows considered, the amplification of turbulence is dictated by the mean density ratio. In the axisymmetric flow, streamwise vortex stretching/contraction, however, alters the amplification slightly.

An experiment has been performed using a facility, which simulates the safety depressurization system (SDS) and in-containment refueling water storage tank (IRWST) of APR1400, an advanced PWR being developed in Korea, to investigate the dynamic load resulting from the blowdown of steam from a steam generator through a sparser. The influence of the key parameters, such as air mass, steam pressure, submergence, valve opening time, and pool temperature, on frequency and peak toads was investigated. The blowdown phenomenon was analyzed to find out the real cause of the initiation of bubble oscillation and discrepancy in frequencies between the experiment and calculation by conventional equation for bubble oscillation. The cause of significant damping was discussed and is presumed to be the highly tortuous flow path around bubble. The Rayleigh-Plesset equation, which is modified by introducing method of image, reasonably reproduces the bubble oscillation in a confined tank. Right after the completion of air discharge the steam discharge immediately follows and it condenses abruptly to provide low-pressure pocket. It may contribute to the negative maximum being greater than positive maximum. The subsequently discharging steam does not play as at the driving force anymore.

Some open jet wind tunnels have been operating under limitations due to large pressure fluctuations at some wind tunnel speeds. The Hyundai Aero-acoustic full scale Wind Tunnel (HAWT), which was completed in 1999, shows that most of the specifications were fulfilled but wind tunnel pulsations at some wind speeds were observed. Hyundai Motor Company started the wind tunnel modification in order to solve this problem in 2001. After the modification work the amplitude of pressure fluctuation was reduced and below required level over full wind speed range. Aero-acoustic performance, e.g. background noise, as well as aerodynamic performance were improved after this work.

Unlike ordinary train, the HST(High Speed Train) is operated at a very high speed, which may cause pressure transient problems when the HST is passing through a station. In the present study, the wind pressure caused by the passing HST was measured in the Cheonan HST station and compared with the numerical simulations. For the measurement, the HST was passing through the station at speeds of 240 km/h north bound and 150 km/h south bound. MEMS based differential pressure transducers are used to measure pressure variation at various locations in the station. It is shown from the results that measured data are in good agreement with CFD simulation with moving mesh technique for the train movement. With the present validation of CFD simulation, the CFD simulation may effectively aid the design of future HST station.