• Nuclear Engineering and Technology
• /
• v.28 no.2
• /
• pp.103-117
• /
• 1996

A mathematical model of vapor explosion propagation is presented. The model predict two-dimensional, transient flow fields and energies of the four fluid phases of melt drop, fragmented debris, liquid coolant and vapor coolant by solving a set of governing equations with the relevant constitutive relations. These relations include melt fragmentation, coolant-phase-change, and heat and momentum exchange models. To allow thermodynamic non-equilibrium between the coolant liquid and vapor, an equation of state for oater is uniquely formulated. A multiphase code, TRACER, has been developed based on this mathematical formulation. A set of base calculations for tin/water explosions show that the model predicts the explosion propagation speed and peak pressure in a reasonable degree although the quantitative agreement relies strongly on the parameters in the constitutive relations. A set of calculations for sensitivity studies on these parameters have identified the important initial conditions and relations. These are melt fragmentation rate, momentum exchange function, heat transfer function and coolant phase change model as well as local vapor fractions and fuel fractions.

• Journal of the Korean Geotechnical Society
• /
• v.36 no.11
• /
• pp.149-156
• /
• 2020

Permeation grouting effectively enhances soil strength and decreases permeability of soil; however, the flow of grout is heavily affected by anisotropy of hydraulic conductivity in layers. Therefore, this study investigates the effect of permeability anisotropy on the effective radius of horizontal permeation grout using computational fluid dynamics (CFD). We modeled the horizontal permeation grout flow as a two-phase viscous fluid flow in porous media, and the model incorporated the chemical diffusion and the viscosity variation due to hardening. The numerical simulation reveals that the permeability anisotropy shapes the grout bulb to be elliptic and the dissolution-driven diffusion causes a gradual change in grout pore saturation at the edge of the grout bulb. For the grout pore saturations of 10%, 50% and 90%, the horizontal and vertical radii of grout bulb are estimated when the horizontal-to-vertical permeability ratio varies from 0.01 to 100, and the predictive model equations are suggested. This result contributes to more efficient design of injection strategy in formation layers with permeability anisotropy.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.48 no.9
• /
• pp.711-716
• /
• 2020

Injection visualization of heated mixed simulant droplets based on hydrocarbon fuel was performed under supercritical state environment. Mixed simulant consisted of Decane and Methylcyclohexane with different critical pressure and critical temperature. Flows injected into the supercritical state environment created droplet by Rayleigh breakup mechanism, and the Oh number and Re number were determined to confirm the breakup area. The temperature of the mixed simulant varied from Tr=0.49 to Tr=1.34. The flow rate was maintained at 0.7 to 0.8 g/s. Droplet became shorter in breakup length as heated and into a lumped form. Second droplet was formed and when Tr=1.34, the phase was not visible in the supercritical state with local unsteady flow.

• Economic and Environmental Geology
• /
• v.44 no.2
• /
• pp.161-170
• /
• 2011

The multiphase flow simulator, CHEMPS, was developed based on the fractional flow approach reported in the petroleum engineering literature considering fully three phase flow in physically and chemically heterogeneous media. It is a extension of MPS developed by Suk and Yeh (2008) to include the effect of wettability on the migration of NAPL. The fractional flow approach employs water, total liquid saturation and total pressure as the primary variables. Most existing models are limited to two-phase flow and specific boundary conditions when considering physically heterogeneous media. In addition, these models focused mainly on the water-wet media. However, in a real system, variations in wettability between water-wet and oil-wet media often occur. Furthermore, the wetting of porous media by oil can be heterogeneous, or fractional, rather than uniform due to the heterogeneous nature of the subsurface media and the factors that affect the wettability. Therefore, in this study, the chemically heterogeneous media considering fractional wettability as well as physically heterogeneous media were simulated using CHEMPS. In addition, the general boundary conditions were considered to be a combination of two types of boundaries of individual phases, flux-type and Dirichlet type boundaries.

• Journal of Energy Engineering
• /
• v.8 no.2
• /
• pp.309-316
• /
• 1999

A numerical study for the turbulent reacting flow in an orimulsion gasifier has been carried out to analyze the characteristics of chemical reaction by orimulsion droplets. In this study, our interest has been focused on the effect of oxidizer to orimulsion ratio, which is one of the key parameters of gasification operation, as well as on the distribution of chemical species. In addition, we have conducted numerical calculations to understand the effect of droplet size, spray angle and injection velocity of fuel so as to acquire the basic information on the operating range of orimulsion gasifier. The result of numerical calculations showed that the gas composition of CO and H$_2$concentrations was the highest when the oxidizer to orimulsion ratio was about 0.88 and the reactivity of orimulsion increased as the droplet size reduced with proper spray angle. Also, we have carried out the analysis on the orimulsion gasification in the 100 ton/day-scale gasifier based upon the prior study in order to obtain the basic data for the proper operating condition using orimulsion feed.

• Journal of Ocean Engineering and Technology
• /
• v.26 no.2
• /
• pp.20-25
• /
• 2012

During drilling, the precipitation velocity of cuttings within an annulus depends on the density and configuration of the cuttings, and on the density, viscosity, and rheological characteristics of the drilling fluid. In directional drilling in particular, it is difficult to adjust and control the cuttings. In contrast to vertical drilling, it is very important to evaluate the flow characteristics of a drilling flow field. However, research on the transfer features of cuttings is inadequate. In this study, in order to identify transfer features of cuttings, an experiment was performed under wide-ranging conditions by constructing a slim hole annulus ($44mm{\times}30mm$) device. In this experiment, the particle volume fraction were influenced by particle size, particle concentration within the flow, pipe rotation, flow volume, and inclination of the annulus. In addition, a mathematical formula for volumetric concentration was deduced and compared to the test results and behavior of cuttings under the other drilling condition was made to be predicted. Therefore, this study can provide meaningful data for vertical and horizontal drilling, and for directional drilling.

• Journal of the Korean Institute of Gas
• /
• v.24 no.2
• /
• pp.1-8
• /
• 2020

Centrifugal pump is a facility that transfers energy to fluid through centrifugal force, which is usually generated by rotating the impeller at high speed, and is a major process facility used in many LNG production bases such as vaporization seawater pump, industrial water and fire extinguishing pump using seawater. to be. Currently, pumps in LNG plant sites are subject to operating conditions that vary depending on the amount of supply desired by the customer for a long period of time. Pumps in particular occupy a large part of the consumption strategy at the plant site, and if the optimum operation condition is not available, it can incur enormous energy loss in long term plant operation. In order to solve this problem, it is necessary to identify the performance deterioration factor through the flow analysis and the result analysis according to the fluctuations of the pump's operating conditions and to determine the optimal operation efficiency. In order to evaluate operation efficiency through experimental techniques, considerable time and cost are incurred, such as on-site operating conditions and manufacturing of experimental equipment. If the performance of the pump is not suitable for the site, and the performance of the pump needs to be reduced, a method of changing the rotation speed or using a special liquid containing high viscosity or solids is used. Especially, in order to prevent disruptions in the operation of LNG production bases, a technology is required to satisfy the required performance conditions by processing the existing impeller of the pump within a short time. Therefore, in this study, the rotation difference of the pump was applied to the ANSYS CFX program by applying the modified 3D modeling shape. In addition, the results obtained from the flow analysis and the curve fitting toolbox of the MATLAB program were analyzed numerically to verify the outer diameter correction theory.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.4
• /
• pp.227-236
• /
• 2016

In the present experimental study, the effect of gas- and liquid-injected methods on the formation of bubble and liquid slug at the merging micro T-junction of a square microchannel with dimensions $600{\mu}m{\times}600{\mu}m$ was investigated. Nitrogen and water were used as test fluids. The superficial velocities of the liquid and gas were in the range of 0.05 - 1 m/s, and 0.1 - 1 m/s, respectively, where the Taylor flow was observed. The bubble length, liquid slug length, bubble velocity, and bubble generation frequency were measured by analyzing the images captured using a high-speed camera. Under similar inlet superficial velocity conditions, in the case of gas injection to the main channel at the merging T-junction (T_gas-liquid), the lengths of the bubble and liquid slug were longer, and the bubble generation frequency was lower than in the case of liquid injection to the main channel at the merging T-junction (T_liquid-gas). On the other hand, in both cases, the bubble velocity was almost the same. The previous correlation proposed using experimental data for T_liquid-gas had predicted the present experimental data of bubble length, bubble velocity, liquid slug length, and bubble generation frequency for T_gas-liquid to be ~24%, ~9%, ~39%, ~55%, respectively.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.35 no.5
• /
• pp.299-308
• /
• 2022

This study was conducted to predict the tendency for heat exchange and boil-off gas (BOG) in a liquefied hydrogen tank under sloshing excitation. First, athe fluid domain excited by sloshing was modeled using a multiphase-thermal flow domain in which liquid hydrogen and hydrogen gas are in the saturated state. Both the the volume of fluid (VOF) and Eulerian-based multi-phase flow methods were applied to validate the accuracy of the pressure prediction. Second, it was indirectly shown that the fluid velocity prediction could be accurate by comparing the free surface and impact pressure from the computational fluid dynamics with those from the experimental results. Thereafter, the heat ingress from the external convective heat flux was reflected on the outer surfaces of the hydrogen tank. Eulerian-based multiphase-heat flow analysis was performed for a two-dimensional Type-C cylindrical hydrogen tank under rotational sloshing motion, and an inflation technique was applied to transform the fluid domain into a computational grid model. The heat exchange and heat flux in the hydrogen liquid-gas mixture were calculated throughout the analysis,, whereas the mass transfer and vaporization models were excluded to account for the pure heat exchange between the liquid and gas in the saturated state. In addition, forced convective heat transfer by sloshing on the inner wall of the tank was not reflected so that the heat exchange in the multiphase flow of liquid and gas could only be considered. Finally, the effect of sloshing on the amount of heat exchange between liquid and gas hydrogen was discussed. Considering the heat ingress into liquid hydrogen according to the presence/absence of a sloshing excitation, the amount of heat flux and BOG were discussed for each filling ratio.

• Nuclear Engineering and Technology
• /
• v.16 no.2
• /
• pp.47-57
• /
• 1984

Mathematically-rigorous time-volume averaged conservation equations were simplified to established the differential equations of THERMIT-6S, which is a two-fluid 3-D code. The difference equations of THERMIT-6S were obtained by discretizing the proceeding set of differential equations. The spatial discretization is characterized by a first-order spatial scheme, donor cell method, and staggered mesh layout. For time discretization, a first order semi-implicit scheme treats implictly sonic terms and terms relating to local transport phenomena and explicitly convective terms. The results were linearized by the Newton-Raphson method. In order to construct the reduced pressure equation, the linearized equations were manipulated so that all variables are coupled between mesh cells through only the pressure variable. By simulating numerically the OPERA-15 experiment, it was found that THERMIT-6S is a very powerful code in predicting reactor behavior after sodium boiling including flow coastdown, reversal flow and flow oscillation.