• Journal of Korea Soil Environment Society
• /
• v.1 no.1
• /
• pp.89-102
• /
• 1996

An integrated model is presented to describe underground flow and mass transport, using a multicomponent multiphase approach. The comprehensive governing equation is derived considering mass and force balances of chemical species over four phases(water, oil, air, and soil) in a schematic elementary volume. Compact and systemati notations of relevant variables and equations are introduced to facilitate the inclusion of complex migration and transformation processes, and variable spatial dimensions. The resulting nonlinear system is solved by a multidimensional finite element code. The developed code with dynamic array allocation, is sufficiently flexible to work across a wide spectrum of computers, including an IBM ES 9000/900 vector facility, SP2 cluster machine, Unix workstations and PCs, for one-, two and three-dimensional problems. To reduce the computation time and storage requirements, the system equations are decoupled and solved using a banded global matrix solver, with the vector and parallel processing on the IBM 9000. To avoide the numerical oscillations of the nonlinear problems in the case of convective dominant transport, the techniques of upstream weighting, mass lumping, and elementary-wise parameter evaluation are applied. The instability and convergence criteria of the nonlinear problems are studied for the one-dimensional analogue of FEM and FDM. Modeling capacity is presented in the simulation of three dimensional composite multiphase TCE migration. Comprehesive simulation feature of the code is presented in a companion paper of this issue for the specific groundwater or flow and contamination problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.1
• /
• pp.95-101
• /
• 2014

In this study, an ISM (immersed solid method) was used for investigating the mass flow rate and efficiency of an involute-gear pump featuring very high rotational speed. For considering circulation flow at the gear pump and housing, fluid flow was assumed as turbulent, and the rotational speed of the gear pump increased under the condition of constant pressure at both the inlet and outlet. The efficiency and mass flow rate of the gear pump were studied by varying its rotational speed and the clearance between the gear tip and the housing. In the simulation results, as the rotational speed were increased, the average mass flow rate and efficiency increased. Furthermore, as the clearance between the gear tip and the housing was increased, the average mass flow rate and efficiency decreased. The efficiency was 85.11, 90.94, and 93.62 at rotational speeds of 6,000 rpm, 8,000 rpm, and 10,000 rpm, respectively, under the condition that there was no clearance. In addition, the efficiency was 93.62, 93.29, and 92.74 at clearances of 0 m, 0.00001 m, and 0.00003 m 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.

• 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 Society of Naval Architects of Korea
• /
• v.35 no.3
• /
• pp.1-13
• /
• 1998

In the present paper, the hydrodynamic characteristics of three dimensional hydrofoils moving with a constant speed below the free surface using a higher-order boundary element method based on 9-node Lagrangian curvilinear elements are investigated. A bi-quadratic spline scheme is employed to improve the numerical results on the free surface. To validate the present scheme, the calculated results are compared with the analytic solutions for a submerged sphere and a spheroid showing a good agreement. For the validation of the hydrofoil study, the computed lift and drag of a hydrofoil having $NACA64_{1}A412$ section with aspect ratio(A.R.) of 4 are compared with the experimental data by Wadlin et al.[28]. The comparison covers a number of variations of angle of attack and submergence depth. Then, using an A.R. hydrofoil with NACA0012 section, the free surface on the lift and drag are investigated and these are compared with the previous results. The wave elevations and patterns created by the aforementioned submerged bodies are also investigated with Froude numbers and submergences.

• Journal of Bio-Environment Control
• /
• v.15 no.2
• /
• pp.138-148
• /
• 2006

This study was carried out to evaluate the structural safety fur both the attached wind-protecting wall in greenhouse and the detached one installed outside. Regarding the attached wind-protecting wall in greenhouse, the analysis was conducted by doing a fluid-structure coupled field analysis using both CFX-5.7 and ANSYS 8.1 and also under the design condition of an instantaneous maximum wind velocity of $30.9m{\cdot}s^{-1}$. Three kinds of the width ranged from 30 to 90cm were considered in this study. With regard to the detached wind-protecting wall, the structural saffty was analyzed under the pressure difference of 1,117 Pa which corresponded to a wind velocity of $50m{\cdot}s^{-1}$ and the analytical results were also compared with theoretical ones. The result showed that there was little difference in the distribution of velocity overall and total pressure on the lateral side according to the width of the attached wind-protecting wall, but greenhouse with wind-protecting widths of 30 to 60cm has been reinforced to the extent of about 11% when compared with the case of being without the wall. The result also showed that the detached wind-protecting wall with a main-column interval of 3m was not stable so that it was necessary for the detached wind-protecting wall to be adequately reinforced to secure structural stability. Finally, there was great difference between analytical results and theoretical studies. The difference meant that there was some possibility of including errors when a theoretical study was done in three dimensional structure.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

• Solar Energy
• /
• v.12 no.3
• /
• pp.70-83
• /
• 1992

The interaction between the surface radiation and the mixed convection transport from an isolated thermal source, with a uniform surface heat flux input and located in a rectangular enclosure, is stuied numerically. The enclosure simulates a practical system such an air cooled electric device, where an air-stream flows through the openings on the two vertical walls. The heat source represents an electric component located in such an enclosure. The size of this cavity is $0.1[m]{\times}0.1[m]$. The inlet velocity is assumed as 0.07[m/s] and the inlet temperature is maintained as $27^{\circ}C$. The inflow is kept at a fixed position. Laminar, two dimensional flow is assumed, and the problem lies in the mixed convection regime, governed by buoyancy force and surface readiation. The significant variables include the location of the out-flow opening, of the heat source and the wall emissivity. The basic nature of the resulting interaction betwwn the externally induced air stream and the buoyancy-driven flow generated by the source is investigated. As a result, the best location of the heat source to make the active heat transfer is 0.075[m] from the left wall on the floor. The trends observed are also discussed in terms of heat removal from practical systems such as electric circuitry.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.9
• /
• pp.126-132
• /
• 2016

This study numerically investigates the characteristics of chemical reactions and thermal deformation in a steam reformer. These phenomena are significantly affected by the high-temperature burner gas and the process gas conditions. Because the high temperature of the burner gas ranges from 800 to 1000 K, the reformer tubes undergo substantial thermal deformation, eventually resulting in structural failure. Thus, it is necessary to understand the characteristics of the reaction and thermal deformation under the operating conditions to evaluate the reformer tubes for sustainable, stable operation. Extensive numerical simulations were carried out using commercial CFD code (ANSYS FLUENT/MECHANICA Ver. 13.0) while considering three-dimensional turbulent flows and combined heat transfer including conduction, convection, and radiation. Structural analysis considering conjugated heat transfer between solid tubes and fluid flows was conducted using the Fluid-Solid Interaction (FSI) method. The results show that when the injection temperature of the process gas and burner gas decreased, the hydrogen production rate decreased significantly, and thermal deformation decreased by at least 15 to 20%.