• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.216-227
• /
• 2004

A new cavitating model by using bubble size distribution based on bubbles-mass has been proposed. Both liquid and vapor phases are treated with Eulerian framework as a mixture containing minute cavitating bubbles. In addition vapor phase consists of various sizes of vapor bubbles, which are distributed to classes based on their mass. The bubble number-density for each class was solved by considering the change of the bubble-mass due to phase change as well as generation of new bubbles due to heterogeneous nucleation. In this method, the bubble-mass is treated as an independent variable, and the other dependent variables are solved in spatial coordinates and bubble-mass coordinate. Firstly, we employed this method to calculate bubble nucleation and growth in stationary super-heated liquid nitrogen, and bubble collapse in stationary sub-cooled one. In the case of bubble growth in super-heated liquid, bubble number-density of the smallest class based on its mass is increased due to the nucleation. These new bubbles grow with time, and the bubbles shift to larger class. Therefore void fraction of each class is increased due to the growth in the whole class. On the other hand, in the case of bubble collapse in sub-cooled liquid, the existing bubbles are contracted, and then they shift to smaller class. It finally becomes extinct at the smallest one. Secondly, the present method is applied to a cavitating flow around NACA00l5 foil. Liquid nitrogen and liquid oxygen are employed as working fluids. Cavitation number, $\sigma$, is fixed at 0.15, inlet velocities are changed at 5, 10, 20 and 50m/s. Inlet temperatures are 90K in case of liquid nitrogen, and 90K and 1l0K in case of liquid oxygen. 110K of oxygen is corresponding to the 90K of nitrogen because of the same relative temperature to the critical one, $T_{r}$=$T/T_c^{+}$. Cavitating flow around the NACA0015 foils was properly analyzed by using bubble size distribution. Finally, the method is applied to a cavitating flow in an inducer of the LE-7A hydrogen turbo-pump. This inducer has 3 spiral foils. However, for simplicity, 2D calculation was carried out in an unrolled channel at 0.9R cross-section. The channel moves against the fluid at a peripheral velocity corresponding to the inducer revolutions. Total inlet pressure, $Pt_{in}$, is set at l00KPa, because cavitation is not generated at a design point, $Pt_{in}$=260KPa. The bubbles occur upstream of the foils and collapse between them. Cavitating flow in the inducer was successfully predicted by using the bubble size distribution.

• 한국지구물리탐사학회:학술대회논문집
• /
• 1999.08a
• /
• pp.92-116
• /
• 1999

Tomographic approaches to image underground structure using electrical properties, can be divided into DC resistivity, electromagnetic, and radar tomography, based on the operating frequency. DC resistivity and radar tomography methods have been recently applied to site investigation for engineering purpose in Korea. This paper review these two tomography methods, through the case histories acquired in Korea. As another method of borehole radar survey, borehole radar reflection method is included, and its inherent problem and solution are discussed, how to find the azimuth angle of reflector using direction-finding-antenna. Since the velocity anisotropy of radar wave has been commonly encountered in field data, anisotropic radar tomography is discussed in this paper. In DC resistivity tomography, two subjects are focussed, electrode arrays, and borehole effect owing to the conductive fluid in borehole. Using the numerical modeling data, various kinds of electrode ways are compared, and borehole effect is illustrated. Most of the case histories presented in this paper are compared with known geology, core logging data, and/or Televiewer images.

• Journal of the Society of Naval Architects of Korea
• /
• v.55 no.3
• /
• pp.252-264
• /
• 2018

The calculation of steady-state cavitating flows around Supercavitating Underwater Bodies (SUB's), which consist of a circular disk head (cavitator), a conical fore-body, a cylindrical middle-body and either a boat-tail or a flare-tail, are carried out. To calculate the axisymmetric cavitating flow, used is a commercial computational fluid dynamics code based on the finite volume method, Fluent. From the analysis of numerical results, the cavity and drag, affected by the fore-body and tail of the SUB's, are investigated. Firstly, the effect of the fore-body shape is investigated with the same disk cavitator and a cylindrical rear-body of fixed diameter. Then with the same cavitator and a fixed fore-body, the effect of the rear-body shape is investigated. Before the cavity generated by the cavitator covers the slant of fore-bodies sufficiently, the larger the cone angle of the fore-body(i.e., the shorter the slant length), the larger the drag and the slower the development of cavity. After the cavity covers the fore-body completely so that the pressure drag component of the body is vanished, the characteristics of drag-velocity curves are identical. Also, as the tail angle is bigger, the cavity generated by the cavitator is suppressed further and the drag becomes larger. The peak of the drag appears for the flare-tail, i.e., when the tail angle is positive(+). On the contrary, the trough of the drag appears for the boat-tail, i.e., when the tail angle is negative(-). When the tail angle is 5 degrees, the peak of the drag appears at the body speed of 80m/s and the value of the drag is 43% larger than that at the design speed of 100m/s. When the tail angle is -5 degrees, the trough of the total drag appears at 75m/s and that drag is 30% smaller than that of the cavitator, which means the rest of the body has a negative drag.

• Journal of Environmental Health Sciences
• /
• v.36 no.5
• /
• pp.402-410
• /
• 2010

The objectives of this study were to evaluate the air quality surrounding an indoor swimming pool, to estimate the cancer risk based on the airborne exposure to trihalomethanes (THMs), and to examine the ventilation efficiency by Computational Fluid Dynamics (CFD). Chlorine and THMs were measured poolside, and in the staff room and reception area. The indoor swimming pool was modeled using the Airpak program, with ventilation drawings and actual survey data. Temperature, flow and mean age of the air were analyzed. Levels of chlorine poolside, and in the staff room, and reception area were $203\;{\mu}g/m^3$, $5\;{\mu}g/m^3$, and $10\;{\mu}g/m^3$, respectively. Chloroform was the dominant THM in all sampling sites and mean concentrations were $16.30\;{\mu}g/m^3$, $0.51\;{\mu}g/m^3$, and $0.06\;{\mu}g/m^3$ poolside, in the staff room and reception area, respectively. Bromodichloromethane and Dibromochloromethane levels were respectively estimated as $10.3\;{\mu}g/m^3$ and $1.7\;{\mu}g/m^3$ poolside, $1.3\;{\mu}g/m^3$ and $0.1\;{\mu}g/m^3$ in the staff room, and were not detected in the reception area. The cancer risks from inhalation exposure to THMs were estimated between $3.37{\times}10^{-7}$ and $1.84{\times}10^{-5}$. A short circulation phenomenon was observed from the supply air vents to the exhaust air vents located in the ceiling. A high temperature layer was formed within one meter of the ceiling, and a low temperature layer was formed under this layer due to the low velocity and high temperature of the supply air, and the improper locations of the supply air vents and exhaust air vents. The stagnation was evident at the above adult pool and the mean age of the air was 22 minutes. Disinfection by-products in the indoor swimming pool were present in higher concentrations than in the outdoor air. In order to increase the removal of pollutants, adjustment was required of the supply air volume and the supply/exhaust position.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.2
• /
• pp.92-108
• /
• 2017

Present work shows a numerical method to analysis of interaction analysis between solitary wave and onshore bridge. Numerical simulation is carried out by TWOPM-3D (three-dimensional one-field model for immiscible two-phase flows), which is based on Navier-Stokes solver. To do this, the solitary wave is generated numerically in numerical wave channel, and numerical results and experimental results were compared and analyzed in order to verify the applicability of force acting on an onshore bridge. From this, we discussed precisely the characteristics of horizontal and vertical forces (uplift and downward forces) changes including water level and velocity changes due to the variation of solitary wave height, water depth, onshore bridge's location and type, and number of girder. Furthermore, It is revealed that the maximum horizontal and vertical forces acting on the girder bridge show different varying properties according to the number of girder, although each maximum force acting on the girder bridge is proportional to the increasement of incident solitary wave height, and the entrained air in the fluid flow affects the vertical force highly.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.5
• /
• pp.346-353
• /
• 2016

In the design conditions of some research reactors, the siphon phenomenon can cause continuous efflux of water during pipe rupture. A siphon breaker is a safety device that can prevent water efflux effectively. However, the analysis of the siphon breaking is complicated because many variables must be included in the calculation process. For this reason, a simulation program was developed with a user-friendly GUI to analyze the siphon breaking easily. The program was developed by MFC programming using Visual Studio 2012 in Windows 8. After saving the input parameters from a user, the program proceeds with three steps of calculation using fluid mechanics formulas. Bernoulli's equation is used to calculate the velocity, quantity, water level, undershooting, pressure, loss coefficient, and factors related to the two-phase flow. The Chisholm model is used to predict the results from a real-scale experiment. The simulation results are shown in a graph, through which a user can examine the total breaking situation. It is also possible to save all of the resulting data. The program allows a user to easily confirm the status of the siphon breaking and would be helpful in the design of siphon breakers.

• Journal of Biomedical Engineering Research
• /
• v.21 no.4
• /
• pp.363-372
• /
• 2000

The present study investigated flow dynamics of a straight elastic blood vessel under sinusoidal flow conditions in order to understand influence of wall motion and impedance phase angle(time delay between pressure and flow waveforms) on wall shear stress distribution using computational fluid dynamics. For the straight elastic tube model considered in the our method of computation. The results showed that wall motion induced additional terms in the axial velocity profile and the pressure gradient. These additional terms due to wall motion reduced the amplitude of wall shear stress and also changed the mean wall shear stress. Te trend of the changes was very different depending on the impedance phase angle. As the wall shear stress increased. As the phase angle was reduced from 0$^{\circ}$to -90$^{\circ}$for ${\pm}$4% wall motion case, the mean wall shear stress decreased by 10.5% and the amplitude of wasll shear stress increased by 17.5%. Therefore, for hypertensive patients vulnerable state to atherosclerosis according to low and oscillatory shear stress theory.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.11 no.3
• /
• pp.143-150
• /
• 2011

Energy loss at manholes, often exceeding friction loss of pipes under surcharged flow, is considered as one of the major causes of inundation in urban area. Therefore, it is important to analyze the head losses at manholes, especially in case of surcharged flow. The stream characteristics were analyzed and head loss coefficients were estimated by using the computational fluid dynamics(CFD) model, FLUENT 6.3, at surcharged square manhole in this study. The CFD model was carefully assessed by comparing simulated results with the experimental ones. The study results indicate that there was good agreement between simulation model and experiment. The CFD model was proved to be capable of estimating the head loss coefficients at surcharged manholes. The head loss coefficients with variation of the ratio of manhole width(B) to inflow pipe diameter(d) and variation of the drop height at surcharged square manhole with a straight-path through were calculated using FLUENT 6.3. As the ratio of B/d increases, head loss coefficient increases. The depth and head loss coefficient at manhole were gradually increased when the drop height was more than 5cm. Therefore, the CFD model(Fluent 6.3) might be used as a tool to simulate the water depth, energy losses, and velocity distribution at surcharged square manhole.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.8
• /
• pp.449-455
• /
• 2019

The increase in impermeable pavement from recent urbanization has resulted in an increase in surface runoff. The surface runoff has also increased the burden of the existing drainage system. This drainage system has structural limitations in that the catchment area is reduced by the waste particles transported with the surface runoff. In addition, the efficiency of the drainage system is decreased. To overcome these limitations, a new type of drainage system with a drainage layer was developed and applied. In this study, various volume porosity and permeability of the lower drainage layer were simulated using ANSYS CFX, which is a three dimensional computational fluid dynamics program. The results showed that the outlet velocity of the 35% volume porosity was faster than that of the 20% and 50% cases, and there was no relationship between the volume porosity and drainage performance. The permeability of the drainage layer can be determined from the particle size of the material, and a simulation of five conditions showed that 2 mm sand grains are most suitable for workability and usability. This study suggests appropriate values of the volume porosity and particle size of the drainage layer. This consideration can be advantageous for reducing and preventing flood damage.

• Clean Technology
• /
• v.25 no.2
• /
• pp.161-167
• /
• 2019

Several types of reactors have been used during the past decade to perform fast pyrolysis of biomass. Among the developed fast pyrolysis reactors, fluidized bed reactors have been widely used in the fast pyrolysis process. In recent years, experimental studies have been conducted on the characteristics of biomass fast pyrolysis in a spouted bed reactor. The fluidization characteristics of a spouted bed reactor are influenced by particle properties, fluid jet velocity, and the structure of the core and annulus. The geometry of the spouted bed reactor is the main factor determining the structure of the core and annulus. Accordingly, to optimize the design of a spouted bed reactor, it is necessary to study the pyrolysis characteristics of biomass. However, no detailed investigations have been made of the fast pyrolysis characteristics of biomass in accordance with the geometry of the spouted bed reactor. In this study, fast pyrolysis experiments using Jatropha curcas L. seed shell cake were conducted in a conical spouted bed reactor to study the effects of reaction temperature and reactor cone angle on the product yield and pyrolysis oil quality. The highest energy yield of pyrolysis oil obtained was 63.9% with a reaction temperature of $450^{\circ}C$ and reactor cone angle of $44^{\circ}$. The results showed that the reaction temperature and reactor cone angle affected the quality of the pyrolysis oil.