A study has been made of the condition to maintaining Taylor-Proudman column flows in a compressible rotating fluid, which is driven by small mechanical and/or thermal perturbations imposing on the container wall in the basic state of isothermal rigid body rotation. The Rossby and system Ekman numbers are assumed to be very small. The Taylor-Proudman column flow can be produced when energy parameter, e, becomes constant on the whole flow region. Energy balance concept, related to energy parameter, and its physical interpretation are given with comprehensive discussions.

Bifurcation of unstable symmetric flow patterns to stable asymmetric ones in laminar sudden-expansion flow has been numerically investigated. Computations were carried out for an expansion ratio of 3 and over a range of the flow Reynolds numbers by using numerical methods of second-order time accuracy and a fractional-step method that guarantees divergence-free flowfields at all times. The critical Reynolds number above which bifurcation of pitchfork type to asymmetric flow pattern takes place is lower in a flow with a higher expansion ratio, in agreement with the previously reported results. The bifurcation diagrams show that the bifurcation takes place at a Reynolds number, $Re_c = 86.3$, higher than the value that has been reported. The lower critical Reynolds number may be due to deficiencies in their computations which employed SIMPLE-type relaxation methods rather than the initial-value approach of the present study. Characteristics of the flow development during the transition to asymmetric stable flow have been investigated by using spectral analysis of the velocity signals obtained by the simulations.

The thermocapillary flow in a differentially heated cylindrical liquid bridge under steady rotation of the hot disk is considered in the limit of zero capillary number. Steady flow states and their three-dimensional stability are calculated numerically. A linear stability analysis reveals that the most dangerous perturbations are oscillatory with azimuthal wavenumber m=1 or m=2 depending on the parameters.

This experimental study concerns the characteristics of a transitional flow in a concentric annulus with a diameter ratio of 0.52, whose outer cylinder is stationary and inner one rotating. The pressure losses and skin-friction coefficients have been measured for the fully developed flow of a 0.2 % aqueous solution of sodium carbomethyl cellulose (CMC) at a inner cylinder rotational speed of $0{\sim}600$ rpm. The transitional flow has been examined by the measurement of pressure losses, to reveal the relation of the Reynolds numbers with the skin-friction coefficients, in the laminar and transitional flow regimes. The occurrence of transition has been checked by the gradient change of pressure losses and skin-friction coefficient with respect to the Reynolds numbers. The increasing rate of skin-friction coefficient due to the rotation is uniform for laminar flow regime, whereas it is suddenly reduced for transitional flow regime and, then, it is gradually declined for turbulent flow regime. Consequently, the critical(axial-flow) Reynolds number decrease as the rotational speed increases. Thus, the rotation of inner cylinder promotes the early occurrence of transition due to the onset of taylor vortices.

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 for the fundamental understanding of vapor explosion phenomena. Therefore, KAERI launched a real material experiment called TROI using 20 kg of UO2 and ZrO2 to investigate the vapor explosion phenomena. In addition, 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 number of instruments were used to measure the physical change occurring during the vapor explosion. In this experiment, the vapor explosion generated by molten fuel water interaction is visualized using high speed camera and the pressure behavior accompanying the explosion is investigated.

Korea Atomic Energy Research Institute (KAERI) launched an intermediate scale steam explosion experiment named Test for Real cOrium Interaction with water (TROI) using reactor material to investigate whether the corium would lead to energetic steam explosion when interacted with cold water at low pressure. The melt-water interaction is confined in a pressure vessel with the multi-dimensional fuel and water pool geometry. The cold crucible technology, where the mixture of powder in a water-cooled cage is heated by high frequency induction, is employed. In this paper, results of the first series of tests ($TROI-1{\sim}5$) were discussed. The ZrO2 jets with 5kg mass and 5cm diameter were poured into the 67cm deep water pool at $30{\sim}95^{\circ}C$. Either spontaneous steam explosions or quenching was observed. The morphology of debris and pressure wave profiles clearly indicates the each case.

An experimental study was carried out to identify the various regimes of natural convective boiling and to determine the Critical Heat Flux(CHF) on a 70mm square surface which is inclined at $180^{\circ}$(upward), $90^{\circ}, \;45^{\circ}$. The heater block made of copper with cartridge heaters is submerged in a water tank with windows for visualization. As the heat flux increases from $100kW/m^2$ to $1.1MW/m^2$, the heat transfer regime migrates from the nucleate boiling to film boiling and results in a rapid heat up of the heater block. An explosive vapor generation on the heated surface, whose size and frequency are characterized by the heat flux, is visualized by using a digital camcorder with $512{\times}512$ pixel size at 30fps.

It is well known that an air bubble trapped in water emits light at its collapse robustly with a proper forcing amplitude of ultrasound. Instability mechanism which causes deviation from sphericity of bubble wall was investigated theoretically. The rapid change of the bubble wall velocity which is both dependent on the forcing amplitude, was found to be a major factor of instability of the interface. The Rayleigh-Taylor instability which occurs when rapid acceleration is directed from the lighter towards the heavier fluid is found to be not related to the instability of the sonoluminescing gas bubble. A good agreement between the calculation results and experimental data is found.

The present article deals with the numerical calculations for the inter-spray impingement of two diesel sprays under the conditions of high injection pressure. The Wave model involving the cavitation effect inside the nozzle was used for describing the atomization process. In particular, a hybrid model for drop collision was newly suggested in this study and compared with the O'Rourk's model, which has been widely used for diesel sprays. The impingement angles of 60 and 90 degrees were considered for simulation of non-evaporative diesel sprays. The calculated results for tip penetration were compared with experimental data and the Sauter Mean Diameter(SMD) characteristic was analyzed. It was concluded that the hybrid model slightly shows better agreement with experimental data than the O'Rourke's model. However, the more elaborate study should be needed for better understanding of spray-to-spray impingement phenomena.

Interfacial pressure jump terms based on the physics of phasic interface and bubble dynamics are introduced into the momentum equations of the two-fluid model for bubbly flow. The pressure discontinuity across the phasic interface due to the surface tension force is expressed as the function of fluid bulk moduli and bubble radius. The consequence is that we obtain from the system of equations the real eigenvalues representing the void-fraction propagation speed and the pressure wave speed in terms of the bubble diameter. Inversely, we obtain an analytic closure relation for the radius of bubbles in the bubbly flow by using the kinematic wave speed given empirically in the literature. It is remarkable to see that the present mechanistic model using this practical bubble radius can indeed represent both the mathematical well-posedness and the physical wave speeds in the bubbly flow.

In this study, the velocity characteristics of liquid elements formed by two impinging jets is analysed using double pulse image capturing technique. For the droplets formed by low speed impinging jets, the droplet velocities are higher with smaller azimuthal and impingement angle. The maximum droplet velocities are about 25 % lower than jet velocity. With an increase of azimuthal angle, the shedding angles increases but remains lower than azimuthal angle. The velocities of ligaments formed by high speed impinging jets gradually decreases with an increase of azimuthal angle. The maximum ligament velocities are about 40 % lower than jet velocity. Higher impingement angles produce lower ligament velocities. The shedding angles of ligament almost increases with the same value of azimuthal angle, which implies that the moving direction of ligaments is radial from the origin as the impingement point.

The numerical simulation of the particle dispersion in the vortical flows provides insight into the mechanism of particle-fluid interaction. The simulation results show that the mixing layers are characterized by the large-scale vortical structures undergoing pairing process. The particle dispersion is strongly influenced by the large-scale structures and the particle sizes. The analysis shows that the mixing layers grows like a step-function.

The size distributions of electrospray droplets from the Taylor cone in cone-jet mode are directly measured by using a freezing method and a transmission electron microscope (TEM) image processing technique. These results are compared with the data obtained by an aerodynamic size spectrometer (TSI Aerosizer DSP). The use of glycerol seeded with NaI and a freezing method make it possible to sample droplets with their original sizes preserved. Since pictures of droplets are taken with TEM with very low vapor pressure of the solution, evaporation is suppressed by freezing. For liquid flow rates below 1 nl/sec, the measured droplet diameters by the TEM image processing technique and the aerosizer are in the range of 0.25 to $0.32{\mu}m$ and 0.30 to $0.40{\mu}m$, respectively. Comparing the TEM data with the aerosizer measurements, it has been revealed that the TEM image processing technique can afford more accurate values of droplet size distributions in the submicron range of 0.1 to $0.4{\mu}m$.

The effect of ambient gas (steam) condensation on swirl spray characteristics were studied experimentally for low subcooling condition of the liquid. The configuration of the liquid(water) sheet and the breakup modes were examined. Also variation of the discharge coefficient, breakup length, local and the cross-sectional area-averaged SMD of droplets with the liquid flow(injection) rate were obtained. The perforation breakup mode appears dominant with condensation while the aerodynamic wave breakup mode is dominant without condensation(in the air environment). The discharge coefficient, breakup length and the mean drop sizes decrease in a same manner with increasing of the liquid flow rate for both cases(with and without condensation). The condensation effects are insignificant with the discharge coefficient. However, the local and cross-sectional area-averaged SMD are larger and the breakup length becomes shorter in the steam environment. The spray angle predicted from the volumetric flux distribution along the radial direction of the sprays in the steam environment becomes larger with condensation.

We present a new interface tracking method for computing two-phase flow. This method is based on VOF method coupled with Level set method. The method is verified to calculate an interfacial curvature accurately as well as to achieve volume conservation during the whole computation period. We apply the present method to calculate a falling drop. The calculated shape and terminal velocity of the falling drop showed good agreement with the data reported in the literature. Also, the present method was proven to be applicable to drop-wall collision phenomenon.

This experimental study was conducted to figure out the characteristics of convective heat transfer non boiling vertical downward flow with polymer additives. This experiment was studied in diameter, 800mm heating length and $1{\times}10^5 W/m^2$ heat flux. The polymer concentration ranged 0ppm to 500ppm with corresponding from superficial liquid velocity 1.25m/s to 2.5m/s in non bo vertical up and downward flow. Experimental results show that the characteristics of convective transfer was a strong function of polymer concentration and it has decreased with increasing polymer concentration in non boiling up and vertical downward flow.

Experimental and analytical researches have been conducted on the twin-fluid atomizers for better droplet breakup during the past decades. But, the studies on the disintegration mechanism still present a great challenge to understand the drop behavior and breakup structure. In an effort to describe the aerodynamic behavior of the sprays issuing from the internal mixing counter-swirling nozzle, the spatial distribution of axial (U) radial (V) and tangential (W) components of droplet velocities are investigated across the radial distance at several axial locations of Z=30, 50, 80, 120 and 170mm, respectively. Experiments were conducted for the liquid flow rates which was kept constant at 7.95 g/s and the air injection pressures were varied from 20 kPa to 140 kPa. Counter-swirling internal mixing nozzles manufactured at angles of $15^{\circ},\;30^{\circ},\;45^{\circ}$ and $60^{\circ}$ the central axis with axi-symmetric tangential-drilled holes was considered. The distributions of velocities and turbulence intensities are comparatively analyzed. PDPA is installed to specify spray flows, which have been conducted along the axial downstream distance from the nozzle exit. Ten thousand of sampling data was collected at each point with time limits of 30 second. 3-D automatic traversing system is used to control the exact measurement. It is observed that the sprays with all swirl angle have the maximum SMD for on air injection pressure of 20 kPa and 140 kPa with centerline, respectively. The nozzle with swirl angle of $60^{\circ}$ has vest performance.

The atomization characteristics of an annular liquid (water) sheet of small radius with a core gas (air) flow were studied. Different sizes of annular gaps (0.2, 0.4 and 0.8 mm) were tested to find the effect of liquid sheet thickness on SMD. The inner diameter of the gas port for the core gas flow was 4 mm. Cross-section averaged SMD was measured for various liquid and gas velocities. Regions of the SMD decrease with the increase of the liquid velocity always existed regardless of the liquid sheet thickness. This attributes to the transition of the flow patterns of spray and also to the aerodynamic interaction between the atomizing gas and the ripples on the liquid sheet surface.

In this study, we analysed tree surface flow by using the experimental and numerical method with a different surfactant concentration. We compared numerical solution with experimental results for one-dimensional model. The result shows that in general the tree surface velocity can well be reproduced by the one-dimensional model for various surfactant concentration.

Natural convection of a magnetic fluid is different from that of Newtonian fluids because magnetic body force exists in an addition to gravity and buoyancy. In this paper, natural convection of a magnetic fluids(W-40) in a cubic cavity is examined by numerical and experimental method. One side wall was kept at a constant temperature($25^{\circ}C$), and the opposite side wall was also held at a constant but lower temperature($20^{\circ}C$). Under above conditions, various magnitudes of the magnetic fields were applied up. GSMAC scheme is used for a numerical method, and the thermo-sensitive liquid crystal film(R20C5A) is utilized in order to visualize wall-temperature distributions as an experimental method. This study has resulted in the following fact that the natural convection of a magnetic fluids is controlled by the direction and intensity of the magnetic fields.

A computational investigation of the effect of the electromagnetic force(or Lorentz force) on the flow behavior around a circular cylinder, a typical model of bluff bodies, is conducted. Two-dimensional unsteady flow computation for $Re=10^2$ is carried out using a numerical method of finite difference approximation in a curvilinear body-fitted coordinate system by solving the momentum equations including the Lorentz force as a body force. The effect of the spatial variations of the Lorentz forcing region and forcing direction along the cylinder circumference is investigated. The numerical results show that the Lorentz force can effectively suppress the flow separation and oscillation of the lift force of the circular cylinder cross-flow, leading to the reduction of the drag.

The effects of pleating parameters on the performance of cylindrical cartridge air filters was experimentally studied. The tested filters are 150 mm in outer diameter and 700 mm in length. As the pleat count increases, the pressure drop across a cartridge filter is decreased for a constant pleat depth and flow rate. This is due to the increased filtration area which brings out less pressure loss of the filter. However, in the case of the filters having the pleat depth of 40 mm, the inner diameter of the cartridges is so small that the pressure loss is greatly increased.

A new type baghouse with tangential inlet and inner tube was developed and it's performances were evaluated. Experiments with variable shapes of inner tube were performed to suggest an optimum shape of inner tube which might decrease the dust loading onto bag filter surface and the amount of reentrained particles by internal gas flow. The dust loading was lower when inner tubes parallel to outer casing exist. When an inner tube covering around the lower portion of bag filters was used, the dust particles detached from the filter surface by cleaning process were reentraind by internal recirculating gas flow.

In the present study, flow visualizations, hot-wire velocity measurements and computational fluid dynamics were performed in order to determine complicated air flow characteristics in a neonate incubator. In this study, following conclusions can be made: (1) The flow visualization technique developed in the present study revealed an enough qualitative information for the flow field in the neonate incubator. Flow structures in a neonate incubator with a realistic three-dimensional shape was successfully visualized the present study. (2) Results from the flow visualization were relatively in good agreements with those obtained from the computational fluid dynamics. (3) Velocities very near the neonate measured by the hot-wire anemometer were relevant to those obtained from the computational fluid dynamics. (4) Temperatures were higher at the neck region and the medial aspect of both thighs, but lower in both extremities. (5) Small vortices between the neonate and the mattress might interfere with convective and evaporative heat transfers on the neonate's surface. In the fluid dynamic aspect, it is important to eliminate the formation of these small vortices for the design of incubator chamber.

This paper reports on the development of high aspect ratio structure and 3-D integrated process for MEMS-based micro gas turbines. To manufacture high aspect ratio structures, Deep Reactive Ion Etching (DRIE) process have been developed and optimized. Specially, in this study, structures with aspect ratios greater than 10 were fabricated. Also, wafer direct bonding and Infra-Red (IR) camera bonding inspection systems have been developed. Moreover, using glass/silicon wafer direct bonding, we optimized the 3-D integrated process.

In a clothes dryer, various thermo-fluid flow phenomena occur such as the heat and mass transfer in the process of removing moisture from clothes, the flow field generated by the fan, and the various flow characteristics from the complex flow paths. The study and understanding of such phenomena is an important factor in increasing the performance of dryers. In this study, as part of a dryer research, the flow field inside a vented dryer was measured using PIV, which the result will be used as the basic material in analyzing the various flow phenomena.

A shrouded rotating disk airflow has a simple figure on geometric basis, but has various and complicated forms of flow. This flow type can be applied to many turbo devices such as information storage device(optical disk). Circumferential velocity frequency in the middle plane between disk and shroud wall is measured using laser Doppler velocimeter. Solid body region of flow was founded when low Reynolds number relatively. Through the informations of the experimental results. we could examine the number and distribution of the vortices. When Reynolds number $3.80{\times}10^5$ there is a dominant frequency of which vortices number is 5.

Rotating disks are used in various machines such as floppy disks, hard disk, turbines and circular sawblades. The problems of vibrations of rotating disks are important in improving these machines. Many investigators have dealt with these problem. Specially, vibrations of a rotating flexible disk taking into account the effect of air is difficult problem in simulation. The governing equation of a rotating flexible disk coupled to the surrounding fluid is investigated by a simple mathematical model. And several important parameters concerned with the stability of a rotating flexible disk are defined. Coupling strength between air and rotating flexible disk is proportional to square of disk radius directly and square root of the all of bending rigidity, disk density and thickness inversely. Lift-to-damping coefficient has relation to the onset of disk flutter.

The applications of magnetic fluid can be normally made by 1) using changes of a property of matter caused by applied magnetic field; 2) preserving magnetic fluid at a certain position or in a magnetic fluid keeping the body in a floating condition; 3) controlling the flow of magnetic fluid by means of magnetic field. However, these are usually made by using their methods together. In this study, the natural convection flow of a magnetic fluid in annular pipes is experimentally analyzed. High temperature is kept constantly inside of a circular pipe of experimental model, on the other hand, low temperature is kept constantly outside of it. In experiments, several cases are carried out in order to clarify the fluence of direction and intensity of magnetic fields on the natural convection of magnetic fluid. Therefore magnetic fields are applied in various intensity and up and down directions by permanent magnets.

In this paper, we have described the microchannel fluid behavior in a slot between rotating curvilinear surfaces of revolution using micropolar fluid theory. ]n order to solve this problem, we have used boundary layer equations and applied non-zero values of the microrotation vector on the wall. The results are compared with the corresponding flow problems for Newtonian fluid. Results show that both the velocity distribution and the microrotation component distribution for a micropolar fluid are lower than that of a Newtonian fluid.

The characteristics of heat transfer on the fire clay with microwave heating are numerically investigated using finite element method. The modelled regular hexahedron chamber($50cm{\times}50cm{\times}50cm$) filled with air consists of vertical heat source and sink walls, a fire clay model, and adiabatic plates at the top and bottom walls. With different geometrical aspect ratios of the fire clay model, the heat energy distribution is throughly investigated. The optimal shape of the fire clay for given chamber geometry and microwave power is analyzed.

Electro-Rheological(ER) fluid are suspensions which show an abrupt increase in rheological properties under electric fields. ER effects arise from electrostatic forces between the starch particles dispersed in the electrically insulating silicone oil, induced when an electric field is applied. Yield stress of the fluids were measured on the couette cell type rheometer as a function of electric fields. This paper presents performance analyses of four types of the two parallel-plate. Which have different electrode length and width but same electrode area. On the basis of the pressure drop and flow rate analysis. Four types of the two parallel-plate are designed and manufactured. Using ER fluid, it is possible to directly interface between electric signals and fluid power without moving parts.

The characteristics of the fluid flow and mass transfer in a vertical atmospheric pressure chemical vapor deposition (APCVD) are numerically studied. In order to get the optimal process parameters for the uniformity of deposition on a substrate, Navier-Stokes and energy equations have been solved for the pressure, mass-flow rate and temperature distribution in a CVD reactor. Results show that the thermal boundary condition at the reactor wall has an important effect in the formation of buoyancy-driven secondary cell when radiation effect is considered. Results also show that reduction of the buoyancy effect on the heated reactor improves the uniformity of deposition.

The objective of the present study is to evaluate the effects of the stenting shapes on flow velocity and wall shear stress in angulated coronary stenosis by computer simulation. Coronary angiogram and Doppler ultrasound measurement in the patients with angulated coronary stenosis were obtained. Inlet wave velocity distribution obtained from in vivo intracoronary Doppler data was used for the numerical simulation. Spatial pattern of blood flow velocity and recirculation area were drawn through out the selected segment of coronary models. Wall shear stresses in the intracoronary stent models were calculated from three-dimensional computer simulation. A negative shear stress region, which is consistent with re-circulation area on flow pattern, was noted on the inner wall of post-stenotic area of pre-stenting model. The negative shear stress was disappeared after stenting. Shear stress in the post-stenting model was markedly reduced up to about two orders of magnitude compared to that of the pre-stenting model.

This study was obtained the pressure distribution of a falling body that is deadrise angle $0^{\circ}$ and deadrise angle $5^{\circ}$ upon a water surface by the experiment with the impact machine. The theoretical equation was obtained the air region and the interface and the water region which devide 3 parties between the body and the water surface for an investigation of the complete phenomena. Pressure distributions and histories compare favorably with available experimental data. The numerical results are similar to the experimental results for the impact force type with $Fo(1+cos{\pi}t/tc)$.

An adiabatic counter-current vertical two-phase flow of air and water in narrow rectangular channels with offset strip fm was investigated experimentally. Tests were systematically performed with downward liquid superficial velocities and upward gas velocities covering 0 to 0.06 m/s and 0 to 2.5 m/s ranges, respectively. Two-phase flow regimes were classified by examining the video images of flow patterns in transparent test sections of 760 mm long and 100 mm wide channel with gaps of 3.0 and 5.0 mm. The channel average void fraction was measured by the quick-closing valve method. Unlike the flow regimes in the channels without fin, where bubbly, slug, chum, and annular flow were identified, only bubbly and chum flow regimes were found for the channels with offset strip fin. However the existence of fin in the channels showed negligible effects on the void fraction. Instead counter-current flow limitations were found to happen at lower air superficial velocity once offset strip fin was introduced in narrow rectangular channels.

The CHFG (Critical Heat Flux in Gap) test results have been evaluated to quantify the critical power in hemispherical narrow gaps and a new correlation has been developed. The CHFG test results have shown that increases in the gap thickness and pressure lead to an increase in critical power. The pressure effect on the critical power was found to be much milder than predictions by CHF correlations of other researches. From the CHFG test results, a new correlation on critical power in the hemispherical gap has been developed using the non-dimensional parameters as follows: $$\frac{qCHF}{{\rho}g^hfg}{\cdot}4\sqrt{\frac{{\rho}_g^2}{g{\sigma}{\Delta}{\rho}}=\frac{0.1042}{1+0.1375({\rho}g/{\rho}l)^{0.21}(D/s)}$$ The developed correlation has been expanded to apply the spherical geometry using the Siemens/KWU's correlation.

Vortices terminating at free surface have been investigated extensively. Most of investigations, however, are focused on surface parallel vortices and little has been known about surface normal vortex or columnar vortex. Visualized experimental results utilizing LIF technique are discussed for the purpose of characterization of columnar vortex interacting with a clean and a contaminated free surfaces and a solid body interface in the present investigation. The results reveal that surface tension changes due to surface contamination although bulk viscosity remains constant and eventually the behavior of a columnar vortex interacting with a contaminated free surface and a solid body interface are totally different from the clean free surface case.

The flow around free end of a finite circular cylinder(FC) embedded in an atmospheric boundary layer has been investigated experimentally. The experiments were carried out in a closed-return type subsonic wind tunnel with varying aspect ratio of the finite cylinder mounted vertically on a flat plate. The wake structures behind a 2-D cylinder and a finite cylinder located in a uniform flow were also measured for comparison. Reynolds number based on the cylinder diameter was about Re=20,000. A hot-wire anemometer was employed to measure the wake velocity and the mean pressure distributions on the cylinder surface were also measured. The flow past the FC free end shows a complicated three-dimensional wake structure and flow phenomenon is quite different from that of 2-D cylinder. The three-dimensional flow structure was attributed to the downwashing counter rotating vortices separated from the FC free end. As the FC aspect ratio decreases, the vortex shedding frequency is decreased and the vortex formation length is increased compared to that of 2-D cylinder. Due to the descending counter-rotating twin-vortex, in the region near the FC free end, regular vortex shedding from the cylinder is suppressed and the vortex formation region is hardly established. In the wake center region, the mean velocity for the FC located in atmospheric boundary layer has large velocity deficit, compared to that of uniform flow.

Control of drag force on a circular cylinder using a detached splitter plate is numerically studied for laminar flow. A splitter plate with the same length as the cylinder diameter(d) is placed horizontally in the wake region. Its position is described by the gap ratio(G/d), where G represents the gap between the cylinder base point and the leading edge of the plate. The drag varies with the gap ratio; it has the minimum value at a certain gap ratio for each Reynolds number. The drag sharply increases past the optimum gap ratio; this seems to be related to the sudden change in the bubble size in the wake region. This trend is consistent with the experimental observation currently available in case of turbulent flow. It is also found that the net drag coefficient significantly depends on the variation of base suction coefficient.

High-resolution simulations using vortex methods have been performed for simulating unsteady viscous flows around an impulsively started square cylinder. In order to investigate the phenomenon from laminar to transition flow, simulations are performed for Reynolds numbers 25, 50, 150 and 250. At extremely low Reynolds number, flow around a square cylinder is known to separate at the trailing edges rather than the leading edges. With an increase of Reynolds number, the flow separation at the leading edges will be developed. The main flow characteristics of developing recirculation region and separations from leading and trailing edges are studied with the unsteady behavior of the wake after the cylinder starts impulsively. A notable change in the flow evolution is found at Re=150, that is, it is shown that the flow separations begin at both leading and trailing edges of the square cylinder. On the other hand, when Re=250, the strong secondary vorticity from the rear surfaces of the square cylinder increases the drag coefficient as the primary vortex layer is pushed outwards. The comparisons between results of the present study and experimental data show a good consistency.

This study presents a feedback control methodology for suppression of the vortex shedding from a circular cylinder in a uniform flow. A rotational oscillation is applied as a controlled forcing and the lift coefficient ($C_L$) is used as a feedback signal. A feedback control concept is made based on the phase relation between the rotation velocity and $C_L$ at 'lock-on', The phase between the forcing and the vortex formation is changed $180^{\circ}$ from the phase of enhancing the lock-on state. This concept is examined by solving the Van del Pol equation. The results are satisfactory.

Three dimensional flow characteristics in a liquid fuel ramjet combustor are investigated using PIV method. The combustors have two rectangular inlets that form 90 degree each other. Three guide vane is installed in each rectangular inlet to improve the flow stability. We made three cases of test combustors in which those inlet angles are 30 degree, 45 degree and 60 degree. Each combustor easily changes the size of combustor's recirculation zone with the replacement of combustors dome. The experiments are performed in the water tunnel test with the same Reynolds number in the case of Mach 0.3 at inlet. PIV software is developed to measure the flow field in the combustor and the accuracy of developed PIV program is verified with rotating disk experiment and standard data. The experimental results show that the two main streams from rectangular inlet collide near the plane of symmetry and generate two large longitudinal vortex, A large and complex three-dimensional recirculating flow is measured in the recirculation zone.

The KSR-III Main Propulsion System configuration of the liquid oxygen (LOX) feed line is analyzed. This feed line includes a tighter radius and cavitation venturi for flow mass flow-rate passive control. There were concerns that these configurations might generate a great flow distortion at the engine interface. Also both the pressure drop at the feed line and any presence of separation area are a great concern according to the propellant flow. To resolve these issues, a computational fluid dynamic analysis was conducted to determine the flow field in the LOX feed lines.

The CSCM Upwind method and Material Transport Analysis (MTA) have been used to predict the thermal response and ablation rate for non-charring material to be used as thermal protection material (TPM) in KSR-III test rocket nozzle. The thermal boundary conditions such as cold wall heat-transfer rate and recovery enthalpy for MTA code are obtained from the upwind Navier-Stokes solution procedure. The heat transfer rate and temperature variations at rocket nozzle wall were studied with shape change of the nozzle surface as time goes by. The surface recession was severely occurred at nozzle throat and this affected nozzle performance such as thrust coefficient substantially.

This paper presents the water-hammer effect due to the rapid opening and closing of isolation valve and thruster valve in the spacecraft propulsion system. The single propellant feed system was modeled to investigate the maximum peak pressure due to the water-hammer effect. The test parameters are tank supply pressure, shape and throat length of orifice and line length. Kerosene was used as the inert simulant propellant liquid instead of hydrazine. As downstream line length after isolation valve increased from 1.5 to 2.5m, the maximum line-filling water-hammer peak pressure decreased, but the average time interval between peak pressures increased. The maximum line-filling water-hammer peak pressure with orifice was lower than without orifice, and the maximum line-filling water-hammer peak pressure with orifice at the back of isolation valve was lower than with orifice in front of isolation valve. Without orifice, the maximum water-hammer peak pressure due to the rapid opening and closing of the thruster valve was about 126% of tank supply pressure. With orifice, it decreased. As orifice throat length increased, it decreased. The maximum water-hammer peak pressure due to the rapid closing of the thruster valve with converging-diverging orifice was lower than normal orifice. It was found that the orifice as a means of pressure drop was very effective to reduce the water hammer peak pressure at the thruster valve. The results of this study can be used for the design of spacecraft liquid propulsion feed system.

Satellite propulsion system is employed for orbit transfer, orbit correction, and attitude control. The monopropellant feeding system in the low-earth-orbit satellite blowdowns fuel to the thrust chamber. The thrust produced by the thruster depends on fuel amount flowed into the combustion chamber. If the thruster valve be given on-off signal from on-board commander in the satellite, valve will be opened or closed. When the thrusters fire fuel flows through opened thruster valve. Instantaneous stoppage of flow in according to valve actuation produces transient pressure due to pressure wave. This paper describes transient pressure predictions of the KOMPSAT2 propulsion system resulting from latching valve and thrust control valve operations. The time-dependent set of the fluid mass and momentum equations are calculated by Method of Characteristics (MOC).

The current study describes experimental and computational work on the passive control of the steady and unsteady condensation shock waves, which are generated in a transonic nozzle. The bleed slots are installed on the contoured wall of the transonic nozzle in order to control the magnitude of the condensation shock wave and its oscillations. For computations, a droplet growth equation is incorporated into the two-dimensional Navier-Stokes equation systems. Computations are carried out using a third-order MUSCL type TVD finite-difference scheme with a second-order tractional time step. Baldwin-Lomax turbulence model is employed to close the governing equations. An experiment using an indraft transonic wind tunnel is made to validate the computational results. The current computations represented well the experimental flows. From both the experimental and computational results it is found that the magnitude of the condensation shock wave in the bleed slotted nozzle is significantly reduced, compared with no passive control of solid wall. The oscillations of the condensation shock wave are successfully suppressed by a bleed slot system.

This paper presents the computational method for analyzing the compressible flow fields in a high voltage gas circuit breaker. There are many difficult problems in analyzing the gas flow in GCB due to complex geometry, moving boundary, shock wave and so on. In particular, the distortion problem of the grid due to the movement of moving parts can be worked out by the fixed grid technique. Numerical simulations are based on a fully implicit finite volume method of the compressible Reynolds-averaged Navier-Stokes equations to obtain the pressure, density, and velocity through the entire interruption process. The presented method is applied to the real circuit breaker model and the pressure in front of the piston is good agreement with the experimental one.

The thin-layer Navier-Stokes equations are solved for the hypersonic flow over blunt cone configurations with applications to laminar as well as turbulent flows. The equations are expressed in the forms of flux-vector splitting and explicit algorithm. The upwind schemes of Steger-Warming and van Leer are investigated in their ability to accurately predict the heating loads along the surface of the body. A comparison with the second order extensions of these schemes is made and a hybrid scheme incorporating a combination of central differencing and flux-vector-splitting is presented. This scheme is also investigated in its ability to accurately predict heat transfer distributions.

This study depicts an experimental work of the impulsive wave discharged from the open end of several kinds of right-angle bend pipes, which are attached to the open end of a simple shock tube. The weak normal shock wave with its magnitude of Mach number from 1.03 to 1.20 is employed to obtain the impulsive wave propagating outside the open end of the bend pipes. The experimental data of the magnitude of the impulsive wave and its propagation directivity are analyzed to characterize the impulsive waves discharged from the right-angle bend pipes and compared with those from a straight pipe. The results obtained show that a right-angle miter bend considerably reduces the magnitude of the impulsive wave and its directivity toward to the pipe axis, compared with the straight pipe. It is believed that the right angle miter bend pipe can playa role of passive control against the impulsive wave.

Supersonic coaxial, axisymmetric, jets issuing from various kinds of dual coaxial nozzles were experimentally investigated. Four different kinds of coaxial, dual nozzles were employed to characterize the major features of the supersonic, coaxial, dual jets. Two convergent-divergent supersonic nozzles with an impinging angle in the jet axis of the annular jets were designed to have the Mach number 2.0 and used to compare the coaxial jet flows with those discharging from two sonic nozzles. The primary pressure ratio was changed in the range from 4.0 to 10.0 and the assistant jet ratio from 1.0 to 4.0. The results obtained show that the assistant jets from the annular nozzle affect the coaxial jet flows and an increase of both the primary jet pressure ratio and assistant jet pressure ratio produces longer supersonic length of the dual, coaxial jet.

In order to investigate the operating characteristics of a supersonic steam ejector, the axisymmetric, compressible, Reynolds-averaged, Navier-Stokes computations are performed using a finite volume method. The secondary and back pressures of the ejector system with a second throat are changed to investigate their effects on the suction mass flow. Three operation modes of the steam ejector system, the critical mode, subcritical mode and back flow mode, are discussed to predict the critical suction mass flow. The present computations are validated with some experimental results. The secondary and back pressures of the supersonic steam ejector significantly affect the critical suction mass flow. The present computations predict the experimented critical mass flow with fairly good accuracy. A good correlation is obtained for the critical suction mass flow. The present results show that provided the primary nozzle configuration and secondary pressure are known, we can predict the critical mass flow with good accuracy.

There are many tunnels located apart short distance from each other in mountainous country like Korea. Serial tunnels are connected by snow shelter in many cases. This study presents some countermeasures against micro-pressure wave at the tunnel exit using snow shelter. Through 1/60 scale model laboratory test, we find that snow shelter with 3.6m slit and slit cover show the effect of reducing the micro pressure wave to about 50.6%.

The purpose of present study is to investigate for reducing pressure fluctuations in the case of installing the air-shafts on the side wall of the tunnel with small cross-sectional area on conventional line. Experiments were performed with a 1/61-scale moving model rig for the tunnel of 0.764 km length in the condition of tunnel cross-section area of $28 m^2$. According to the results, the maximum pressure fluctuation is reduced by 45 % for 19 air-shafts. This results have the speed-up effects of about 33.4 km/h for the train running in tunnel.

The present study addresses the open end correction associated with the reflection and discharge phenomena of a weak shock wave from an open end of a duct. The open end correction of the weak shock wave is investigated experimentally and by numerical computation. An experiment is made using a simple shock tube with an open end, and computation is performed to simulate the experimental flow field using the unsteady, axisymmetric, compressible, flow governing equations. The results obtained show that an open end correction should be involved for shock wave discharge and reflection problems generated from the exit of the duct with an open end baffle plate. With a baffle plate less than three times the duct diameter, it is found that the open end correction is a function of both the diameter of the baffle plate and normal shock wave magnitude. However, for a baffle plate larger than three times the duct diameter, it is independent of the baffle plate diameter. The present computations predict the results of shock tube experiment with good accuracy. A new empirical equation for prediction of the open end correction is found for the weak shock reflection and discharge phenomena occurring at the open end of the duct with and without a baffle plate.

This study addresses a computational work of the impulsive wave which is discharged from the open end of a pipe. An initial compression wave inside the pipe is assumed to propagate toward atmosphere. The over pressure and wave-length of the initial compression wave are changed to investigate the characteristic values of the impulsive wave. The second order total variation diminishing (TVD) scheme is employed to solve the axisymmetric, compressible, unsteady Euler equations. The relationship between the initial compression wave form and impulsive wave is characterized in terms of the peak pressure of the impulsive wave and its directivity. The results obtained show that for the initial compression wave of a large wave-length the peak pressure of the impulsive wave does not depend on the over pressure of the initial compression wave, but for the initial compression wave of a very short wave-length, like a shock wave, the peak pressure of the impulsive wave is increased with an increase in the over pressure of the initial compression wave. The directivity of the impulsive wave to the pipe axis becomes significant with a decrease in the wave-length of the initial compression wave.

The present study describes a practical estimation procedure about the pantograph under several severe aerodynamic load conditions. As the operating speed of the Korean Train Express(KTX) reaches 350km/h, structural safety at various conditions should be examined at the design stage. In the present study, a compact and reliable procedure is developed to get aerodynamic loads on each part of the pantograph regarding the typhoon condition, the train/tunnel interaction, the train/train interaction and the side wind condition. In the estimation procedure, 3-dimensional steady and unsteady CFD simulation around the high speed train facilitates assigning the external local flow condition around the pantograph. The procedure is verified using the results of the low speed wind tunnel test at JARI and applied to 7 flow conditions and 4 operation configurations.

The paper presents an improved way of aerodynamic quality in Korean Very High Speed Railway, The pantograph model being under development dissatisfies the required grade of aerodynamic lift force. So the present work proposes modified configurations of panheads to maintain consistent aerodynamic characteristics. Analysis has been performed using commercial CFD program. Simulation based analysis has been conducted with two different models. One is to attach the thin plate on the crossbar and the other is the use of trapezoidal cross section in contact strip. Various length of thin plate is simulated for flow velocity and acceptable value of plate length is selected which satisfy the necessary average lift force. Aerodynamic variation on the panheads strip is studied.

When a train runs into a tunnel at high-speed, aerodynamic drag suddenly increases and the booming noise is generated at the exit of tunnel. The noise shape is very important to reduce the aerodynamic drag in tunnel as well as on open ground, and the micro-pressure wave that is a source of booming noise is dependent on nose shape, especially on area distribution. In this study, the nose shape has been optimized employing the response surface methodology and the axi-symmetric compressible Navier-Stokes equations. The optimal designs have been executed imposing various conditions of the aerodynamic drag and the micro-pressure wave on object functions. The results show that the multi-objective design was successful to decrease micro-pressure wave and aerodynamic drag of trains.

In the present study, flow characteristics of turbulent oscillatory flow in a square-sectional $180^{\circ}$ curved duct are investigated experimentally. In order to measure wall shear stress and pressure distributions, experimental studies for air flow are conducted in a square-sectional $180^{\circ}$ curved duct by using the LDV system with the data acquisition and the processing system. The wall shear stress measuring point bend angle of the $150^{\circ}$ and pressure distribution of the inlet (${\phi}=0^{\circ}$) to the outlet (${\phi}=180^{\circ}$) at $10^{\circ}$ intervals of the duct. The results obtained from the experimentation are summarized as follows: A wall shear stress value in an inner wall is larger than that in an outer wall, except for the phase angle (${\omega}t/{\pi}/6$) of 3, because of the intensity of secondary flow. The pressure distributions are the largest in accelerating and decelerating regions at the bend angle(${\phi}$) of $90^{\circ}$ and pressure difference of inner and outer walls is the largest before and after the ${\phi}=90^{\circ}$.

A finite element code based on P2P1 tetra element has been developed for the large eddy simulation (LES) of turbulent flows around a complex geometry. Fractional 4-step algorithm is employed to obtain time accurate solution since it is less expensive than the integrated formulation, in which the velocity and pressure fields are solved at the same time. Crank-Nicolson method is used for second order temporal discretization and Galerkin method is adopted for spatial discretization. For very high Reynolds number flows, which would require a formidable number of nodes to resolve the flow field, SUPG (Streamline Upwind Petrov-Galerkin) method is applied to the quadratic interpolation function for velocity variables, Noting that the calculation of intrinsic time scale is very complicated when using SUPG for quadratic tetra element of velocity variables, the present study uses a unique intrinsic time scale proposed by Codina et al. since it makes the present three-dimensional unstructured code much simpler in terms of implementing SUPG. In order to see the effect of numerical diffusion caused by using an upwind scheme (SUPG), those obtained from P2P1 Galerkin method and P2P1 Petrov-Galerkin approach are compared for the flow around a sphere at some Reynolds number. Smagorinsky model is adopted as subgrid scale models in the context of P2P1 finite element method. As a benchmark problem for code validation, turbulent flows around a sphere and a MIRA model have been studied at various Reynolds numbers.

A numerical study is performed to investigate the interaction between subsonic axial turbine blade boundary layer and periodically oncoming rotor induced wakes. An implicit scheme for solving the compressible Navier-Stokes equation is developed, which adopts a 4th-order compact difference for spatial discretiztion, a 2nd order Crank-Nicolson scheme for temporal discretization and the dynamic eddy viscosity model as the subgrid scale model. The efficiency and the accuracy of the proposed method are verified by applying to some benchmark problems such as laminar cylinder flow, laminar airfoil cascade flow and a transitional flat plate boundary layer flow. Computational results show good agreements with previous experimental and numerical results. Finally, flow through a stator cascade is simulated at $Re = 7.5{\times}10^5$ without free-stream turbulence intensity. The velocity fields and skin friction coefficients in the transitional region show similar trends with previous boundary layer natural transition.

Pressure drop at a helical square duct orifice is numerically evaluated. The orifice is installed at the entrance of a once-through steam generator tube to suppress flow instabilities. The calculated results are compared with the available experimental correlations, and showed good agreement. Effects of curvature ratio and characteristics of the secondary flow with Reynolds number are reported. Through the numerical simulations, pressure drop mechanisms were well understood inside the compact and complicated orifice geometry.

An account of second-order fractional-step methods and boundary conditions for the incompressible Navier-Stokes equations is presented. The present work has aimed at (i) identification and analysis of all possible splitting methods of second-order splitting accuracy; and (ii) determination of consistent boundary conditions that yield second-order accurate solutions. It has been found that only three types (D, P and M) of splitting methods called the canonical methods are non-degenerate so that all other second-order splitting schemes are either degenerate or equivalent to them. Investigation of the properties of the canonical methods indicates that a method of type D is recommended for computations in which the zero divergence is preferred, while a method of type P is better suited to the cases when highly-accurate pressure is more desirable. The consistent boundary conditions on the tentative velocity and pressure have been determined by a procedure that consists of approximation of the split equations and the boundary limit of the result. The pressure boundary condition is independent of the type of fractional-step methods. The consistent boundary conditions on the tentative velocity were determined in terms of the natural boundary condition and derivatives of quantities available at the current timestep (to be evaluated by extrapolation). Second-order fractional-step methods that admit the zero pressure-gradient boundary condition have been derived. The boundary condition on the new tentative velocity becomes greatly simplified due to improved accuracy built in the transformation.

An analytical and numerical examination of second-order fractional-step methods and boundary condition for the incompressible Navier-Stokes equations is presented. In this study, the compatibility condition for pressure Poisson equation and its boundary conditions, stability, and numerical accuracy of canonical fractional-step methods has been investigated. It has been found that satisfaction of compatibility condition depends on tentative velocity and pressure boundary condition, and that the compatible boundary conditions for type D method and approximately compatible boundary conditions for type P method are proper for divergence-free velocity for type D and approximately divergence-free for type P method. Instability of canonical fractional-step methods is induced by approximation of implicit viscous term with explicit terms, and the stability criteria have been founded with simple model problems and numerical experiments of cavity flow and Taylor vortex flow. The numerical accuracy of canonical fractional-step methods with its consistent boundary conditions shows second-order accuracy except $D_{MM}$ condition, which make approximately first-order accuracy due to weak coupling of boundary conditions.

This paper presents the numerical prediction of sound generated by viscous flow past a circular cylinder. The two dimensional flow field is predicted using FEM based Reynolds-averaged Navier-Stokes solver, and the calculated unsteady fluid field values are utilized by an acoustic code that implements Ffowcs Willianms-Hawkings(FW-H) equation. The integration surface used in acoustic analysis is extended from the cylinder surface to permeable surfaces. The 2D based CFD calculations overpredict the acoustic amplitude, however, if adequate correlation length is used, the predicted acoustic amplitude agrees well with experiment. The predictions using extended integral surface in FW-H equation show results that contain the characteristics of quadrupole - volume integration - noise term, and do not vary seriously with the integral surface location.

A two dimensional hierarchical elements are investigated for a use on the incompressible flow computation. The construction of hierarchical elements are explained through the tensor product of 1-D hierarchical functions, and a systematic treatment of essential boundary values has been developed for the degrees of freedom corresponding to higher order terms. The numerical study for the poisson problem showed that the present scheme can increase the convergence and accuracy of finite element solutions, and can be more efficient than the standard first order with many elements. Also, for Stokes and cavity flow cases, solutions from hierarchical elements showed better resolutions and future promises for higher order solutions.

This paper introduces a pressure correction method for microflow computation. Conventional CFD methods with no slip boundary condition fail to predict the rarefaction effect of the wall when simulating gas microflows in the slip-flow regime. Pressure correction method with an appropriate slip boundary condition is an efficient tool in analyzing microscale flows. The present unstructured SIMPLE algorithm adopts both the classical Maxwell boundary condition and Langmuir boundary condition proposed by Myong. The simulation results of microchannel flows show that the proposed method has an effective predictive capability for microscale flows.

In order to investigate the mixing of dissimilar grades during the arbitrary grade transition in bloom caster, a computational model has been developed. The model is fully transient and consists of two sub models, which account for mixing in the bloom tundish, mixing in the strand. The developed model was verified using concentration histories measured on 1 : 1 scale bloom tundish water model. The result of numerical model showed good agreement with the experimental results of water model. By using this numerical model, the mixing of dissimilar grades in bloom caster has been simulated. As that result, the characteristics of the steel mixing in the bloom tundish and strand was showed and the amount of the intermixed grade bloom was predicted.

The flow between two concentric cylinders, with the inner one rotating and with an imposed pressure-driven axial flow, is studied using numerical simulation. This study considers the identical flow geometry as in the experiments of Wereley and Lueptow[Phys. Fluid, 11 (12), 1999]. They carried out experiments using PIV to measure the velocity fields in a meridional plane of the annulus in detail. When an axial flow is imposed, the critical Taylor number is increased. The axial flow stabilizes the flow field and decreases the torque required to rotate the inner cylinder. The velocity vector fields obtained also show the same flow features found in the experiments of Wereley and Lueptow.

An analytical study is made of transient adjustment process of an initially stationary, stably-stratified fluid in a square. The boundary walls are highly-conducting. The overall Rayleigh number $R_a$ is large. Considerations are given to both opposing (${\delta}w/{\delta}T>0$) and cooperating (${\delta}w/{\delta}T<0$) configurations. The flow character in opposing configuration can be classified into (a) a forced-convection dominaut mode (${\delta}w/{\delta}T>1/\sqrt{2}$), (b) a buoyancy-convection dominant mode ($0<{\delta}w/{\delta}T<1/\sqrt{2}$), and (c) a static mode (${\delta}w/{\delta}T{\cong}1/\sqrt{2}$). Global evolutionary processes are depicted. and physical rationalizations are provided.

Comprehensive numerical computations are made of a homogenous spin-up in a cylindrical cavity with a time-dependent rotation rate. Numerical solutions are acquired to the governing axisymmetric cylindrical Navier-Stokes equation. A rotation rate formula is ${\Omega}_f={\Omega}_i+{\Delta}{\Omega}(1-{\exp}(-t/t_c))$. If $t_c$ is large, it implies that a rotation change rate is small. The Ekman number, E, is set to $10^{-4}$ and the aspect ratio, R/H, fixed to I. For a linear spin-up(${\epsilon}<<$), the major contributor to spin-up in the interior is not viscous-diffusion term but inviscid term, especially Coriolis term, though $t_c$ is very large. The viscous-diffusion term only works near sidewall. But for spin-up from rest, when $t_c$ is very large, viscous-diffusion term affects interior area as well as sidewall, initially. So azimuthal velocity of interior for large $t_c$ appears faster than that of interior for relatively small $t_c$. However, the viscous-diffusion term of interior decreases as time increases. Instead, inviscid term appears in the interior.

This paper presents the methodology for thermal hydraulic analysis of Pressurized Water Reactor (PWR) steam generators. Topics include porous media approach, governing equations, physical models and correlations for solid-to-fluid interaction and heat transfer and numerical solution scheme. Some details about the ATHOS3 code currently used widely for thermal hydraulic analysis of PWR steam generators in the industry are presented. The ATHOS3 code is applied to the thermal hydraulic analysis of steam generator in the Korea YGN 3&4 nuclear power plant and the computed results are presented.

It is very important that piping system shall meet the optimum design requirement as predicted in designing system. If the piping system proved not to meet the requirement in commission it shall be redesigned and corrected till the required capacity is satisfied. which costs much expense. The objective of this study is to understand steady 3-dimensional flow phenomena in a process piping system numerically. 3-dimensional numerical simulations with standard $k-{\epsilon}$ model were carried out by using ALGOR code for three cases of Reynolds number. 2500, 3500 and 4500, based on inlet pipe diameter and three cases of inflow air temperature, $20^{\circ}C,\;50^{\circ}C$ and $100^{\circ}C$.

This paper addresses a numerical method for predicting transient temperature distributions in the wall of a curved pipe subjected to internally thermal stratification flow. A simple and convenient numerical method of treating the unsteady conjugate heat transfer in the non-orthogonal coordinate systems is presented. The proposed method is implemented in a finite volume thermal-hydraulic computer code based on a cell-centered, non-staggered grid arrangement, the SIMPLEC algorithm, a higher-order bounded convection scheme, and the modified version of momentum interpolation method. Calculations are performed for the transient evolution of thermal stratification in two curved pipes, where the one has thick wall and the other has so thin wall that its presence can be negligible in the heat transfer analysis. The predicted results show that the thermally stratified flow and transient conjugate heat transfer in a curved pipe with a finite wall thickness can be satisfactorily analyzed by the present numerical method, and that the neglect of wall thickness in the prediction of pipe wall temperature distributions can provide unacceptably distorted results.

A CFD study was conducted to evaluate the nuclear fuel assembly coolant mixing that is promoted by the flow-mixing vanes on the grid spacer. Four mixing vanes (split vane, swirl vane, twisted vane, hybrid vane) were chosen in this study. A single subchannel of one grid span is modeled using the flow symmetry. The three mixing vanes other than swirl vane generate a large crossflow between the subchannels and a skewed elliptic swirling flow in the subchannel near the grid spacer. The swirl vane induces a circular swirling flow in the subchannel and a negligible crossflow. The split vane and the twisted vane were predicted to result in relatively larger pressure drop across the grid spacer. Since the average turbulent kinetic energy in the subchannel rapidly decreases to a fully developed level downstream of the spacer, turbulent mixing caused by the mixing vanes appears to be not as effective as swirling flow mixing in the subchannel. In summary, the CFD analysis represented the overall characteristics of coolant mixing well in a nuclear fuel assembly with the flow mixing vanes on the grid spacer. The CFD study is therefore quite useful for the development of an advanced flow-mixing vane.

This paper introduces some useful Java technologies for utilizing the Internet in numerical analysis, and suggests one architecture performing efficient numerical analysis in the Intranet by using them. The present work has verified it's possibility by implementing some parts of this architecture with two easy examples. One is based on Servlet-Applet communication, JDBC and Swing. The other is adding multi-threads, file transfer and Java Remote Method Invocation to the former. Through this work it has been intended to make the base for the later advanced and practical research that will include efficiency estimates of this architecture and deal with advanced load balancing.

When supersonic jet impinges on wall from the nozzle, complex flow pattern appears such as Mach disc, expansion fan, and jet boundary. The numerical computation of this supersonic jet is important on flame deflecctor design for launch space especially. In this paper, we analyzed supersonic jet structure impinging on deflector wall using three dimensional steady and unsteady compressible equation and showed temperature and pressure distribution on the wall surface. As a result, some dominant factors of jet flows are discussed for conceptual design of flame deflector.

The objective of the paper was to measure the pressure drop and to investigate the recirculation region of the conical orifices used in Kwang-yang Iron & Steel Company. The flow field with water used as a working fluid was the turbulent flow for Reynolds number of $2{\times}10^4$. The effective parameters for the pressure drop and the recirculation region were the conical orifice's inclined angle (${\theta}$) against the wall, the interval(L) between orifices, the relative angle of rotation(${\alpha}$) of the orifices, the shape of the orifice's hole(circle, rectangle, triangle) having the same area. It was found that the shape of the orifice's hole affected the pressure drop and the flow field a lot, But the other parameters did not make much differences to the pressure drop. The PISO algorithm with FLUENT code was employed.

This paper presents how redevelopment of the boundary layer in a backward-facing step flow is affected by boundary conditions imposed on velocity at the inlet, top and exit of the flow. A two-dimensional, laminar, incompressible flow over a backward-facing step with an open top boundary has been computed by using numerical methods of second-order time and spatial accuracy and a fractional-step method that guarantees a divergence-free velocity field at all time. The inlet velocity profile above the step is of Blasius type. Along the top boundary, shear-tree and Dirichlet conditions on the streamwise velocity were considered and at the exit fully-developed and convective boundary conditions were examined. (The vertical velocity at all boundaries were assumed to be zero explicitly or implicitly.) From the computed flow fields, the reattachment on the bottom side of shear layer separated from the tip of the step and succeeding redevelopment of the boundary layer were investigated.

In this study, spin-up flows in a rectangular container are analysed by using three-dimensional computation. In the numerical computation, we use the parallel computer system of PC-cluster type. We compared our results with those obtained by two-dimensional computation. Effect of velocity and vorticity on the flow is studied. The result shows that two-dimensional solution is in good agreement with the 3-D result. Attention is given to the region where the 3-D flow is significant.

Vortex methods were originally conceived as a tool to model the evolution of unsteady, incompressible, high Reynolds number flows of engineering interest. Recently various methods have been proposed for simulating the diffusion in vortex methods for two-dimensional incompressible flows. We test the diffusion schemes of vortex methods. In this paper we directly compare the particle strength exchange scheme with the vorticity redistribution scheme in tenus of their accuracy and computational efficiency. Comparisons between both viscous models described are presented for short-time runs of impulsively started flows past a circular cylinder for Reynolds number of 60. The particle strength exchange scheme has been shown more accurate and efficient than the vorticity redistribution scheme.

Because of the ozone layer depletion and global wanning, new alternative refrigerants are being developed. Among them, HFC refrigerants are thought promising, but some European countries are arguing that these refrigerants are also harmful to the global wanning. Therefore, natural refrigerants should be considered as an eventual alternative in refrigerators and heat pumps. In the present study, the supercritical gas cooling process are computationally analysed by employing various turbulence models of carbon dioxide in a trans critical refrigeration cycle. The gas cooling process near the critical point experiences a drastic change in thermodynamic and transport properties, thus the heat transfer characteristics would be different from those of two or single phases. Based on the computational results, the correlations to estimate the near-critical heat transfer characteristics will are obtained.

Variations of temperature field in a Hele-Shaw convection cell (HSC) were measured using a holographic interferometry with varying Rayleigh number. Experimental results show a steady flow pattern at low Rayleigh numbers and a time-dependent periodic flow at high Rayleigh numbers. Especially, the period of oscillation at $Ra = 6.35{\times}10^6$ was 62 seconds. Two different measurement methods of holographic interferometry, double-exposure method and real-time method, were employed to measure the temperature field variations of HSC convective flow. In the double-exposure method, unwanted waves can be eliminated and reconstruction images are clear, but transient flow structure cannot be observed clearly. On the other hand, transient flow can be observed and reconstructed well using the real-time method. However, the fringe patterns reconstructed by the real-time method contain more noise, compared with the double-exposure method. The two holographic interferometer techniques employed complementary in this study were proved to be useful for analyzing the temperature field variations of unsteady thermal fluid flows.

A high resolution digital cinematic Particle Image Velocimetry(PIV) has been developed. The system consists of a high speed CCD camera, a continuous Ar-ion laser and a computer with camera controller. To improve the spatial resolution, we adopt a Recursive Technique for velocity interrogation. At first, we obtain a velocity vector for a larger interrogation window size based on the conventional two-frame cross-correlation PIV analysis using the FFT algorithm. Based on the knowing velocity information, more spatially resolved velocity vectors are obtained in the next iteration step with smaller interrogation windows. The correct velocity vector at the first step is found to be critical, so we apply a Multiple Correlation Validation(MCV) technique in order to decrease the spurious vectors. The MCV technique turns out to improve SNR(Signal to Noise Ratio) of the correlation table. The developed cinematic PIV method has been applied to the measurement of the unsteady flow characteristics of a Rushton turbine mixer. A total of 3,245 instantaneous velocity vectors were successfully obtained with 4 ms time resolution. The acquired spatial resolution corresponds the performance of the conventional high resolution digital PIV system using a $1K{\times}1K$ CCD camera.

A non-intrusive Planar Laser-Induced Fluorescence(PLIF) technique was applied to study the turbulent mixing process in a Rushton turbine reactor. Instantaneous and ensemble averaged concentration fields are obtained by measuring the fluorescence intensity of Rhodamine B tracer excited by a thin Nd:Yag laser sheet illuminating the whole center plane of the stirred tank. The gray level images captured by a 14-bit cooled CCD camera can be transformed to the local concentration values using a calibration matrix. The dye injection point was selected at the tank wall with three quarter height (3/4H) from the tank bottom to observe the mixing characteristics in upper bulk flow region. There exist distinct two time scales: the rapid decay of mean concentration in each region after the dye infusion reflects the large scale mixing while the followed slow decay reveals the small scale mixing. The temporal change of concentration probability functions conjectures the two sequential processes in the batch type mixing. An inactive column of water existed above the impeller disk, in which the fluid rotates with the shaft but is isolated from the mean bulk flow.

A GA(Genetic Algorithm) based 3D-PTV technique has been developed. The measurement system consists of three CCD cameras, Ar-ion laser, an image grabber and a host computer. The fundamental of the developed technique was based on that one-to-one correspondence is found between two tracer particles selected at two different image frames taking advantage of combinatorial optimization of the genetic algorithm. The fitness function controlling reproductive success in the genetic algorithm was expressed by a kind of continuum theory on the sparsely distributed particles in space. In order to verify the capability of the constructed measurement system, a performance test was made using the LES data set of an impinging jet. The developed 3D-PTV system was applied to the measurement of flow characteristics of the wake of a circular cylinder.

A new stereoscopic PIV is developed using two CCD cameras, stereoscopic photogrammetry, and a 3D-PTV principle. The wake of a circular cylinder is measured by the developed stereoscopic PIV technique. The B mode vortical structure of the wake over the Reynolds number 300 is clearly seen by the developed technique. The arrangement of the two cameras is based on angular position. The calibration of cameras and the pair-matching of the three-dimensional velocity vectors are based on 3D-PTV technique.

The objective of the present study is to investigate the characteristics for flow and wall shear stress in the aneurysm which is a local dilatation of the blood vessel. The numerical simulation using the commercial software for the laminar and steady flow were carried out over the diameter ratios(ratio of maximum diameter of aneurysm to the diameter of blood vessel) ranging from 1.5 to 2.5 and Reynolds number ranging from 900 to 1800. It was shown that a recirculating vortex occupied the entire bulge with its core located closer to the distal end of the bulge and the strength of vortex increased with increase of the Reynolds number and diameter ratio. Especially, for the Reynolds number of 1800 and diameter ratio of 2.5, the very weak secondary recirculating flow was produced at the left upper of the aneurysm. The position of a maximum wall shear stress was the distal end of the aneurysm(z=18mm) regardless of the Reynolds number and diameter ratios. But the maximum values of the wall shear stress increased in proportion to the increase of Reynolds number and diameter ratio.

Researchers have invectigated nasal flow both numerically and experimentally for centuries. Experimental studies most have suffered from various limitations necessary to allow the measurements to be obtained with available equipment. Nasal airflow can be subdivided into two interrelated categories; nasal airflow resistance and heat and mass transfer between the air stream and the walls of the nasal cavity. In this study, thanks to a new method for model casting by a combination of Rapid prototyping and Solidification of clear silicone, a transparent rectangular box containing the complex nasal cavity is made for PIV experiments. The CBC PIV algorithm is used for analysis.

An experimental study is performed to investigate the characteristics of vortex shedding behind a circular cylinder with serrated fins using hot-wire anemometer. Strouhal numbers which are calculated using outer diameter of a circular cylinder with serrated fins are higher than that of a circular cylinder. Fin thickness and pitch are closely related with vortex shedding frequency and play increasing or decreasing vortex shedding after transient Reynolds numbers. Strouhal numbers using effective diameters which are proposed in this paper agree with that of a circular cylinder. After transient Reynolds number, a trend of Strouhal number can be estimated by checking the ratio of effective diameter to inner diameter.

When a liquid is supplied through a nozzle onto a relatively nonwetting inclined solid surface, a narrow rivulet forms. This work provides novel physical insights into the following phenomena in the rivulet flow that have not been well understood to date. Firstly, the fundamental mechanism behind the transition of a linear rivulet to a droplet flow is investigated. The experiments show that the droplet flow emerges due to the necking of a liquid thread near the nozzle. Based on the observation, it is argued that when the retraction velocity of a liquid thread exceeds its axial velocity, the bifurcation of the liquid thread occurs, and this argument is experimentally verified. Secondly, a discussion on the curved motion of a meandering rivulet is given. This study proposes the contact angle hysteresis as a primary origin of the centripetal force that enables the rivulet's curved motion A simple scaling analysis based on this assumption predicts a radius of curvature which agrees with the experimental observation.

A numerical method using FLUENT code was employed to investigate fluid drag and lift forces on a cylinder in a group of circular cylinders, subjected to a uniform cross flow. The cylinders can be arranged in tandem or in a staggered arrangements relative to the free stream flow. A vortex street behind the cylinder pairs or jets between the cylinders forms according to the arrangements. Vibration on a cylinder can occurs due to vortex shedding, fluid-elastic stiffness and wake galloping. The flow is first investigated and then the forces acting on the cylinder are calculated. The lift and drag forces on an elastically mounted cylinder in the wake of an upstream fixed cylinder arise from the mean flow plus velocity and pressure gradients in the wake. The analytical results of two staggered cylinder were compared with the existing experimental ones for validation of the present method. The analytical results of the forces were in good agreement with the experimental ones. The present method can be used for the analysis of the fluid induced vibration where the group of circular cylinders are subjected to a cross flow.

본 연구를 통해 초음속 전투기 날개의 공력-구조를 동시에 고려한 다학제간 설계를 수행하였다. 공력해석을 위해 사용된 3 차원 Euler Code는 수렴 속도를 개선하기 위해 Multigrid를 적용하였으며, 3차의 transfinite interpolation을 사용하여 O-H type의 공력해석 격자계를 생성하였다. 구조 분야는 절점당 54개의 자유도를 가지는 9 절점 쉘 혼합 유한요소(9-node shell mixed finite element)를 사용하여 해석을 수행하였다. 설계변수는 공력쪽으로 날개의 평면형상에 관련된 변수 3개, 구조쪽은 날개 윗면과 아래면의 표피두께에 관련된 4개의 설계변수 사용하였으며, D-optimality 조건을 만족시키는 실험점들에 대해 공력해석과 구조해석이 연동된 정적 공탄성 해석을 수행한 후, 반응면 기법을 이용하여 목적함수와 제약조건에 대한 반응면을 구성하였다. 단일점 설계를 수행한 후 이를 바탕으로 3개의 설계점을 동시에 고려한 다점 설계를 수행하였으며, 공력만을 고려한 설계 결과와 공력-구조를 동시에 고려한 다학제간 설계결과의 비교를 통해 다학제간 설계의 타당성과 우수성을 입증하였다.

Thermal mass flow meter(TMF) and thermal mass flow controller(MFC) were used to measure and to control the mass flow rate of gases. TMF and MFC were designed for specified working pressure and gas. For the case of different working pressure and gases, the flow rate measurement accuracy decreased dramatically. In this study, a TMF and MFC was tested with three different gases and pressure range from 0.2 MPa up to 1.0 MPa. Effect of specific heat causes to increase flow measurement error as much as ratio of specific heat compared with reference gas. Changing of pressure causes to increase flow rate measurement error about -0.2% as the working pressure decreased 0.1 MPa. Response time of MFC was below 3.12 s for the case of increasing of flow rate. But the response time was increased up to 6.92 s for the case of decreasing of flow rate. When the solenoid valve was fully closed, a initial delay time of output of MFC was increased up to 1.36 s.

Measurement uncertainty should be evaluated according to ISO/IEC 17025. In Flow measurement area, uncertainty evaluation scheme was applied to the reference flow meter, sonic nozzle. Uncertainty was calculated by evaluating various uncertainty factors affected in flow measurement. The expanded uncertainty of the sonic nozzle was 0.21 % (confidence level of 95 %). This evaluation example will be useful in flow measurement uncertainty determination of other flow meters.

Based on the flow characteristics around a piercing cylinder, a free surface-flow velocitmetry which can be used in extremely harsh environment such as molten steel flow was developed. The velocimetry is consisted of finite length cylinder, load detecting elastic plate, electric signal transducer and data acquisition H/W and S/W. Using such a velocimetry, two velocity measurement schemes were established which one is flow resistance detecting scheme and the other is Karman Vortex frequency detecting scheme. For calibration of each scheme, realistic flow water model was used and in followings, detailed calibration processes were explained.

For the evaluation of operability of MOV(Motor Operated Valve), the precision prediction of thrust/torque acting on the valve is important. In this paper, the analytical prediction method of thrust/torque was proposed. The design basis stem thrust calculation typically considers the followings: Packing thrust, Stem rejection load, design basis differential pressure load. In general, test results show that temperature, pressure, fluid type, and differential pressure, independently and combination, all have an effect on the friction factor. The prediction results of thrust/torque are well agrement with dynamic test results.

A performance test facility for turbopump inducer cavitation was developed and the inducer cavitation performance tests were performed. Major components of the performance test facility are driving unit, test section, piping, water tank, and data acquisition and control system. The maximum of testing capability of this facility are as follows: flow rate - 30kg/s; pressure - 13 bar; rotational speed 10,000rpm. This cavitation test facility is characterized by the booster pump installed at the outlet of the pump that extends the flow rate range, and by the pressure control system that makes the line pressure down to vapor pressure. The vacuum pump is used for removing the dissolved air in the water as well as the line pressure. Performance tests were carried out and preliminary data of test model inducer were obtained. The cavitation performance test and cavitation bubble flow visualization were also made. This facility is originally designed for turbopump inducer performance test and cavitation test. However it can be applied to the pump impeller performance test in the future with little modification.

The inducers in liquid-rocket engines are to increase the inlet pressure of the pump to avoid any malfunction due to cavitation. Inducers are typically designed to be operated with some amount of cavitation for the compactness of the turbopump system. Also, inducers are designed to produce low headrise to prevent the decrease of the overall pump efficiency due to the low efficiency of inducers. In the present paper, a computational study on the hydrodynamic behavior of the inducer for the rocket-engine turbopump are presented including the effect of the mass flow rate under the constant rotational speed. As the mass flow rate is decreased, the inducer showed better performance with strong back flows which may have deleterious effects upon the anti-cavitation ability. But the adopted inducer showed very low headrise with high volume flow rates, which may be caused by the small passage area near the trailing edge. The modified version of the present inducer is proposed and numerically evaluated, which in turn showed better results.

The internal flow in the rocket pump inducer of LE-7 engine for H-II rocket was predicted at design and off-design flow rates using CFD code, CFX- Tascflow. In this numerical study, the performance curve of inducer coressponding to flow rates variation and the internal flow in the front of blade leading edge show good agreement between the calculations and the measurements. Backflow is appeared at suction side of leadinge edge tip, and this region is extended to upstream as flowrate decrease. Because of backflow, pressure loss coressponding to meridinal coordinate occupy 50% from inlet domain to leading edge. By this phenomena, pressure loss in front of blade leading edge take a great effect to inducer performance.

Performance and cavitation characteristics of a turbo pump inducer which is designed for oxygen pump is experimentally investigated. It is found that the static performance of inducer is dependent on the location of inlet pressure measurement and that enough distance from inducer is ensured for accurate evaluation of performance. With the increase of flow rate. NPSH tends to decrease as opposed to pump characteristics, which seems to be due to the cavitation region.

The direct simulation Monte Carlo(DSMC) method is applied to investigate steady and unsteady flow fields of a single-stage disk-type drag pump. Two different kinds of pumps are considered: the first one is a rotor-rotor combination, and the second one is a rotor-stator combination. The pumping channels are cut on a rotor and stator. The rotor and stator have 10 Archimedes' spiral blades, respectively. In the present DSMC method, the variable hard sphere model is used as a molecular model, and the no time counter method is employed as a collision sampling technique. For simulation of diatomic gas flows, the Borgnakke-Larsen phenomenological model is adopted to redistribute the translational and internal energies. The DSMC results are in good agreement with the experimental data.

This paper describes aerodynamic preliminary design performance prediction and flow analysis for centrifugal compressor of the marine middle engine turbocharger. The performance characteristics of turbocharger compressor are investigated at various operating conditions using mass flow rate and revolution speed, and computational flow analysis for impeller and diffuser at design point are performed. Preliminary design results correspond to actual compressor geometric values comparatively by applying modified slip factor. Performance prediction and flow analysis results show good agreement with experiments. Therefore, this will provide the performance prediction in preliminary design, and help to increase the design capability for optimized impeller.

The performance test of a centrifugal compressor for APU(Auxiliary Power Unit) which is developed by the collaborative research of KARI and Samsung TechWin has been conducted. The investigated compressor consists of a curved inlet, a centrifugal impeller, a channel diffuser and a plenum chamber. The experiments were carried out in an open-loop centrifugal compressor test rig driven by a turbine. For three different diffusers, overall performance data were obtained at 80%, 90% and 97% of design speed. For the initially designed wedge-type diffuser, test results showed that the compressor was operated at a higher mass flow rate than the design requirement. By reducing the diffuser throat area, the compressor operating range was shifted to lower mass flow rate range. The test result of redesigned wedge-type diffuser showed high pressure loss. To reduce the diffuser loss, diffuser inlet radius was increased and airfoil-type of diffuser was adopted. This airfoil-type diffuser showed reasonal results in terms of design requirement.

This paper is concerned with the effects of periodically approaching upstream wakes on cascade-flow and loss. The reduced frequency of the periodic wakes was varied in the narrow range from 0.5 to 0.7. According to a wake-passing through the cascade, two velocity deficits appear near the boundary layer contour in the downstream from about 60% chord-length. The first velocity deficit is caused by a periodic wake and the second one appears after some delayed time. The second velocity deficit may be interpreted as the results of reattachment of flow-separation. The higher reduced frequency decreases the duration time of separation appearance and the lesser loss of aerofoil is resulted.

For the saving of energy consumption, it is necessary that the high efficient energy transfer machine. We found the optimum operating condition of screw water chiller during this performance test. The specifications of compressor for this performance test are as follows. Built-in volume ratio; Vi=2.4, using refrigerant; HCFC-22, nominal motor output; 60kW(80Hp). We found the test results on water cooled chiller standard test condition that the evaporating capacity is increased 8% maximum and 5.6% under same operating condition. and COP is increased 11 % max. and 8% under same operating condition with change of the evaporating inlet/outlet condition. For the results of these test, we can estimate the optimum operating condition of water cooled type screw chiller.

This study has been conducted to design the high efficiency centrifugal compressor for a HFC-134a. The 2-stage centrifugal compressor consists of inlet guide vanes, two impellers with splitters, a deswirler, a vaneless diffuser and a volute casing. We have designed the high efficiency centrifugal compressor by applying the repeated design procedure including a meanline design, a 3D geometry generation of 1st and 2nd impellers etc. Also, a fluid dynamic calculation of impellers and deswirler have been conducting using a commercial code STAR-CD.

The present work has carried out experimental study on a cross-flow fan system with a simplified vortex wall scroll casing. A cross-flow fan test rig was constructed to obtain pressure rise and volume flow rate for various fan operating conditions. The performance estimation is using a wind tunnel with a motor driven damper for flow rate control and flows are quantitatively visualized by light scattering system with a pulsed laser. Min focus on the visualization is finding a eccentric vortex inside a fan which is a major factor reducing fan efficiency. Comprehensive engineering data are prepared for industrial applications and show a good agreement with a prior work by experimental measurements.

In this study, three-dimensional viscous flow analysis and optimization are presented for the design of a mixed-flow fan. Steady, imcompressible, three-dimensional Reynolds averaged Navier-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. Upwind difference scheme is used for the discretization of the convective term and SIMPLEC algorithm is used as a velocity-pressure correction procedure. The computational results are compared with the results obtained by TASCflow. For the numerical optimization of the design, objective function is defined as a ratio of generation of the turbulent energy to pressure head. Sweep angles are used as design variables.

The three-dimensional spatial structures of impeller flow created by a six bladed Rushton turbine have identified based on the volumetric velocity information from multi-plane stereoscopic PIV measurements. A total of 10 planes with 2 mm space with a 50 mm by 64 mm size of the field of view were targeted. To reduce the depth of focus, we adopted an angle offset configuration which satisfied the Scheimpflug condition. The distortion compensation procedure was utilized during the in situ calibration. Phase-locked instantaneous data were ensemble averaged and interpolated in order to obtain mean 3-D, volumetric velocity fields on a 60 degree sector of a cylindrical ring volume enclosing the turbine blade. Using the equi-vorticity surface rendering, the spatial structure of the trailing vortices was clearly demonstrated. Detail flow characteristics of the radial jet reported in previous studies of mixer flows were easily identified.

The performance and the internal flow of the impeller of the centrifugal chiller compressor with refrigerant R12 as working fluid were studied numerically, using CFD code, CFX-Tascflow, which is commercially available. In this numerical study, the thermodynamic and transport properties of the refrigerant gas were generated by the property program of NIST and linked with main program to extend the capability of the code to refrigerant gases. Numerical study was applied to several mass flow rates near the design mass flow rate at constant rotating speed. Overall performance and flow characteristics of the impeller at impeller exit were investigated. The results were physically reasonable and showed good agreement with experimental measurement at the design flow rate.

The characteristics of flow in dividing regions are precise, therefore their classification is very important not only in industry but also in hydrodynamics. By now, many studies of flow in dividing regions have been performed, but flow characteristics that use visualization in dividing regions have not been studied. The present study of the PIV and the CFD exhibit average velocity distributions, kinetic energy distributions and total pressure distributions etc of the total flow field due to the development of the accurate visualization optical laser and of optical equipment. Also, PIV is accurate with the flows characteristics of the dividing region as continuous analysis is done using input equipment. The study analyzes average velocity vector field, average kinetic energy, x-axis stress distributions, average and total pressure distributions of dividing regions with flow for visualization of the PIV and the CFD measurement in a dividing rectangular duct.

The paper presented some results of a experimental study of natural convection in partitioned 2D square enclosure. The square enclosure consist of two adiabatic vertical walls and the upper cold and the lower hot walls. A partition is positioned perpendicularly at the center of left vertical insulated wall The PIV mesaurements were performed with the variations of the partition length and inclination of enclosure. The working fluid is water with a Prandtl number of 6.996 at $20^{\circ}C$ temperature. A captured images were calculated by using a Cross-Correlation(Multi-frame/Single-exposure) method.

The Holographic Particle Velocimetry system can be a promising optical tool for the measurements of three dimensional particle velocities. In this research, the optical system for double pulse holographic recording and reconstruction of particle images was developed. Validation experiments for the developed system were conducted measuring the velocities of glass beads on a rotating disk. Uncertainty analysis was performed to identify the sources of all relevant errors and to evaluate their magnitude. The measurement results of distance between glass beads, size, and velocities of them using holographic method compared reasonalbly well with the known values within acceptable range of errors.

Two-dimensional, angle-resolved LDV(Laser Doppler Velocimetry) measurements of the turbulent rotating flow field in a confined cylinder have been performed. The configurations of interest are flows between a rotating upper disk with a rod attached by a disk or impeller(${\theta}= 90^{\circ},\;45^{\circ}$) and a stationary lower disk in a confined cylinder. The mean flow velocity as well as the turbulent intensity of the flow field have been measured. The results show that the flow is strongly dependent on the position of the impellers or the disk, negligibly affected by the Reynolds number in turbulent flow. It is observed that the mixing effect of the axial flow impeller(${\theta}= 45^{\circ}$) is better than that of the radial flow impeller(${\theta}= 90^{\circ}$) or a disk.

The purpose of this study is to compare the time mean velocity distribution, the time mean kinetic energy, and the time mean turbulence intensity between vertical and horizontal flow fields in a coaxial circular pipe by PIV measurement. Experiments are performed at a Reynolds number 2,000, measuring regions divided as the section regions A, B, C, D in flow fields. The angle of the high-frequency ultrasonic is selected in the direction of $45^{\circ}$ to the flow axes and it is reflected several times. In results, it is clarified that the effect of gravity is given in the vertical flow field compared with the horizontal flow field and the ultrasonic affects the turbulence enhancement. And kinetic energy and turbulence intensity with ultrasonic are shown slightly bigger than those in flow field without it.

Heat generated from the electronic parts in PDP is undesirable physical properties. To attain optimal arrangement of the electronic parts in PDP, thermal flows in PDP should be analyzed. PIV measurement has been made to quantify the characteristics of the inner flows and outer flows of an actual PDP. The quantity of heat flux from PDP has been estimated using the PIV results. Measurement system consists of Ar-ion laser, CCD camera and an image grabber installed on a host computer.

The purpose of the present study is to measure the viscosity of liquid in the capillary tube viscometer using the unsteady flow concept. The capillary tube viscometer is consisted of a small cylindrical reservoir, capillary tubes, and the mass flow rate measuring system interfaced with computer. Two capillary tubes with 1.152 and 3.002 mm I.D. are used to determine the diameter effects on the viscosity measurements. The instantaneous shear rate and gravitational driving force in the capillary tube are determined by measuring the mass flow rate through the capillary tube instantaneously. The measured viscosities of water and aqueous Separan solution are in good agreement with the reported experimental data.

A study on the ultrasonic effect for turbulence enhancement is carried out in the horizontal flow field of a coaxial circular pipe. A large transparent acryl tank is made to perform several experiments for the above research. The front flow field from jet exit is divided as 4 measuring regions to observe characteristics of the above flow field according to those with and without ultrasonic. An ultrasonic transducer with 2MHz high frequency is used to give them the ultrasonic forcing. Characteristics such as the velocity distribution, the kinetic energy and the turbulence intensity are visualized, observed, examined and considered at Re No. 2000. In results, it is clarified that the ultrasonic increases the turbulence enhancement. And the optimum and harmonious intensity suited to the power of flow is needed to maximize the turbulence enhancement.

The objective of the present study is to visualize the pulsatile flow fields by using three-dimensional computer simulation and the PIV system. A closed flow loop system was built for the steady and unsteady experiments. The Harvard pulsatile pump was used to generate the pulsatile pressure and velocity waveforms. Conifer powder as the tracing particles was added to water to visualize the flow field. Two consecutive particle images were captured by a CCO camera for the image processing at several cross section. The range validation and the area interpolation methods were used to obtain the final velocity vectors with high accuracy. The finite volume predictions were used to analyze three-dimensional flow patterns in the bifurcation model. The results of the PIV experiment and the computer simulation are in good agreement and the results show the recirculation zones and formation of the paired secondary flow distal to the apex of the bifurcated model. The results also show that the branch flow is pushed strongly to the inner wall due to the inertial force effect and helical motions are generated as the flow proceeds toward the outer wall.

This paper presents unsteady computational investigations and wind tunnel tests on the flow field around a square cylinder with a gap between the body and the ground plane. Two-dimensional unsteady, incompressible Navier-Stokes codes are developed for the computation of the viscous turbulent flows. By computing the flow around a square cylinder without ground effect, three two-equation turbulence models are evaluated and the developed code is validated. The results show a good agreement with experimental values and other computational results. Critical gap height at which the formation of Karman vortex streets is interrupted, is demonstrated and another transition regime is pointed out

Some multidimentional generalizations of the Fokker-Planck Equation used by Friedrich and Peinke for description of a turbulent cascade was solved by A.A.Donkov, A.D.Donkov, and G.I.Grancharova. The solutions are two types, isotropic and anisotropic diffusion case. We introduce their methods to solve the Equation and solutions. Furthermore we get the more generalized exact solution as combination of two cases and plot to compare those to experimental results for the isotropic case.

An experimental study was undertaken to investigate the effect of coriolis force for the turbulent flow at low Reynolds numbers in a rotating straight square duct. The study was carried out using a hot-wire anemometer. The flow Reynolds number based on the hydraulic diameter ranged from 4,000 to 18,000 and Rotation number ranged from 0 to 0.196. At Re=9000, developing turbulent flow was calculated for mean velocity and Reynolds stress. Pressure coefficient and energy dissipation spectrum were also calculated.

The flow structures of turbulent shear layer behind oil fences with different tip configurations were investigated experimentally using flow visualization and PIV velocity field measurement. An oil fence was installed in a circulating water channel and the flow structure around the fence tip was mainly analyzed in this experiment. The four tip configurations tested in this experiment are knife edge; semi-circle edge, circular edge and rectangular edge. The 300 instantaneous velocity fields were measured using the single-frame PIV system and they were ensemble averaged to give the mean velocity field and spatial distribution of turbulent statistics. Free stream velocity was fixed at 10ms/sec and the corresponding Reynolds number based on the fence height was Re=4000. As a result, for the oil fence with rectangular edge, the streamwise velocity component was decreased. On the other hand it was increased for the oil fence with circular edge. For all four fences tested in this study, general flow pattern of the lower shear layer is analogous but the upper layer shows difference depending on the tip configurations. The oil fence with circular edge has more diffusive upper shear layer than that of the others. The shear layer of the oil fence with rectangular edge has relatively thin thickness. The oil fence with circular edge was found to be proper shape for tandem fence.

In this study, we used In-line orifice mixer for efficient chemicals mixing in water treatment. The method of using In-line orifice mixer has been already proved the improvement of water treatment efficiency. Numerical study was performed using FLUENT, a commercial code, to standard design and production of effective In-line orifice mixer. As variable for exactly standardizing, a proper ratio between an outer diameter of cone and a diameter of pipe, a distance between cone and orifice, a determination of orifice diameter for an optimal mixing, a distance between injection nozzle's position and cone, Numerical study has been performed for optimal standard and analyzed flow field on a basis of turbulent intensity in an orifice downstream.

Choosing the most suitable mathematical model and relating this to turbulent tangential tensions model are very important in the investigations of turbulent two-phase flow. This paper considers two-fluid scheme. According to it, two phases have their own densities, velocities, and temperatures at any spatial point and at any moment. The equations of motion and heat transfer for each phase are linked with the forces of interaction between two phases. These forces are considered as predominant for the flow. As a closure in the system of motion equations, one modification of $K - {\epsilon}$ turbulent model is worked out. The modification uses two equations for turbulent kinetic energy of the phases and one - for the turbulent energy loss of main phase. This model can be set as a $K_g - K_p -{\epsilon}$ model. The modified model has been tested for both a two-phase non-isothermal flat jet and axially symmetrical jet. The numerical results are compared with the reference data revealing a good agreement between them.

In this study, a program for the off-design performance prediction of multi-stage axial-compressors is developed based on stage-stacking method. To account for the increased losses at off-design conditions, generalized performance curve is applied. The purpose of this study is to investigate the influence of the choice of generalized performance curve and stator exit angle. For this purpose, we tested various generalized performance curves and stator exit angles. In conclusion, Muir's pressure coefficient curve gives a good prediction results regardless of the efficiency curve for a low-stage compressors. On the other hand, for high-stage compressors, The combination of Muir's pressure coefficient curve and Stone's efficiency curve gives a optimistic results. Stator exit angle has a small effect on overall performance curve.

In this study, a program to design a multi-stage axial compressor is developed wi th mean-line analysis and vortex methods. In a preliminary design stage, a method. to design in a short time is needed and mean-line analysis is usually used for this purpose. Arbitrary pressure ratio and reaction can be assigned to generate overall geometry and several vortex methods are adopted to consider the radial distribution of velocity and reaction. The variation of performance, when we use free vortex, forced vortex, and exponential method, is compared and discussed.

The objective of present study is to find the interaction between volute and impeller of the centrifugal turbomachinery with rectangular cross-sectional volute. Flow measurement were taken in shrouded impeller with 12 backward type blades by using a five-hole pressure probe. The measurements are carried out in 7 flow rate, respectively. Primary function of a centrifugal turbo machinery volute is to serve the flow from the impeller and diffuser to pipe system. For the off-design conditions, Influence of pressure distortion was shown by these measurements. As a result, It has caused the decrease of total efficiency of centrifugal turbomachinery. We have also taken data to design volute by these measurements.

The effects of casing shapes on the interaction of the impeller and volute in a small-size turbo-compressor are investigated. Numerical analysis is conducted for the compressor with circular and single volute casings from inlet to discharge nozzle. In order to predict the flow pattern inside the entire impeller, vaneless diffuser, and casing, calculations with a multiple frame of reference method between the rotating and stationery parts of the domain are carried out. For incompressible turbulent flow fields, the continuity and three-dimensional time-averaged Navier-Stokes equations are employed. To predict the performance of two types of casings, the static pressure and loss coefficients are obtained with various flow rates. Also, static pressure distributions around casings are studied for different casing shapes, which are very important to predict the distribution of radial load.

The flow characteristics of the vaned diffuser were complicated with geometric shapes. We have studied the effects of various vaned diffuser configurations, such as divergence angles and rectangular and conical cross-section shapes. Numerical analyses are carried out for the diffuser and casing. The pressure recovery coefficient was calculated to estimate the performance of the diffuser, and then compared with the measure data. Results show that the shapes and the divergence angles of the diffuser strongly influence on the performance of the small-size turbo-compressor.

Mean velocity and Reynolds stress components of the developing turbulent flows in a rotating 90 degree bend with square cross-section were measured by a hot-wire anemometer. Effects of the centrifugal and Coriolis forces generated by the curvature and rotation of bend on the mean motion and turbulence structures are investigated experimentally. Results show that the Coriolis force associated with the rotation of the bend may act both through the mean motion and turbulent structures, thereby changing the pressure fields, mean and turbulent velocities distributions.

A numerical study on a quantitative analogy of the fully developed flow between in a straight square duct rotating about an axis perpendicular to that of the duct and a stationary curved duct of square cross-section is carried out. In order to clarify the similarity of two turbulent flows, the dimensionless parameters $K_{TR} = Re^{1/4}/\sqrt{Ro}$ and the Rossby number, Ro, in a rotating straight duct flow were used as a set corresponding to $K_{TC} = Re^{1/4}/\sqrt{{\lambda}}$ and curvature ratio, ${\lambda}$, in a stationary curved duct flow so that they have the same dynamical meaning as $K_{LR} = Re/\sqrt{Ro}$ and $K_{LC} = Re/\sqrt{{\lambda}}$ of the fully developed laminar flows. For the large values of Ro or A, it is shown that the flow field satisfies the asymptotic invariance property: there are strong quantitative similarities between the two flows such as flow patterns and friction factors for the same values of $K_L$ and $K_T$.