• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2696-2701
• /
• 2007

For the analysis of a two-phase flow, the interaction between two phases such as the interfacial momentum or heat transfer is proportional to the interfacial area. So the interfacial area concentration (IAC) is one of the most important parameters governing the behavior of each phase. This study focuses on the development of a computational fluid dynamics (CFD) code for investigating a boiling flow with a one-group IAC transport equation. It was based on the two-fluid model and governing equations were calculated by SMAC algorithm. For checking the robustness of the developed code, the experiment of a subcooled boiling in a vertical annulus channel was analyzed to validate the capability of the IAC transport equation. As the results, the developed code was confirmed to have the capability in predicting multi-dimensional phenomena of vapor generation and propagation in a subcooled boiling.

• Geomechanics and Engineering
• /
• v.10 no.4
• /
• pp.423-436
• /
• 2016

Currently, there is a great interest in the coupling between multiphase fluid flow and geomechanical effects in hydrocarbon reservoirs and surrounding rocks. The ideal solution for this coupled problem is to introduce the geomechanical effects through the stress analysis solution and implement an algorithm, which assures that the equations governing the flow and stress analyses are obeyed in each time step. This paper deals with the implementation of a program (FORTRAN90 interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators, using a partial coupling algorithm. The explicit coupled hydro-mechanical behavior of Iranian field during depletion and $CO_2$ injection is studied using the soils consolidation procedure available in ABAQUS. Time dependent reservoir pressure fields obtained from three dimensional compositional reservoir models were transferred into finite element reservoir geomechanical models in ABAQUS as multi-phase flow in deforming reservoirs cannot be performed within ABAQUS. The FEM analysis of the reservoir showed no sign of plastic strain under production and $CO_2$ injection scenarios in any part of the reservoir and the stress paths do not show a critical behavior.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.294-298
• /
• 2003

The complex, unsteady, self-sustained pressure oscillations induced by supersonic flow past a rectangular cavity is investigated using numerical simulations. The present numerical study is performed using a parallel, multiblock solver for the two-dimensional, compressible Navier­Stokes equations. Open cavities with length-to-depth (L / D) ratio in the range 0.5 - 3.3 are considered. This paper sheds light on the cavity physics, cavity oscillatory mechanism, and the organisation of vortical structures inside the cavity. The vortex shedding phenomenon, the shear layer impingement event at the aft wall and the movement of the acoustic/compression wave within the cavity are well predicted. The vortical structures· and the source of the acoustic disturbances are found to be located near the aft wall of the cavity. With the increase in the cavity length, strong recompression of the flow near the aft wall leading to a sudden jump in the cavity form drag is observed. The estimated cavity tones are in good agreement with the available semi­empirical relation. Multiple peaks are noticed in deep and long cavities. For the present free­stream Mach number 1.71, it is observed that around L/D=2.0, the cavity oscillatory mechanism changes from the transverse to longitudinal oscillatory mode. The effects of this transition on various fluid dynamics and acoustic properties are also discussed.

• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• Journal of Aerospace System Engineering
• /
• v.12 no.2
• /
• pp.15-22
• /
• 2018

A compressible multiphase flow was investigated using a DIM consisting of seven equations, including the fifth-order MLP and a modified HLLC Riemann solver to achieve a precise interface structure of liquid and gas. The numerical methods were verified by comparing the flow structures of the high-pressure water and low-pressure air in the shock tube. A 2D air-helium shock-bubble interaction at the incident shock wave condition (Mach number 1.22) was numerically solved and verified using the experimental results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.9
• /
• pp.688-698
• /
• 2009

Using the multi-dimensional, multi-phase, nonisothermal Polymer Electrolyte Fuel Cell (PEFC) model presented in Part I, the effects of land/channel flow-field on temperature and liquid saturation distributions inside PEFCs are investigated in Part II. The focus is placed on exploring the coupled water transport and heat transfer phenomena within the nonisothermal and two-phase zone existing in the diffusion media (DM) of PEFCs. Numerical simulations are performed varying the land and channel widths and simulation results reveal that the water profile and temperature rise inside PEFCs are considerably altered by changing the land and channel widths, which indicates that oxygen supply and heat removal from the channel to the land regions and liquid water removal from the land toward the gas channels are key factors in determining the water and temperature distributions inside PEFCs. In addition, the adverse liquid saturation gradient along the thru-plane direction is predicted near the land regions by the numerical model, which is due to the vapor-phase diffusion driven by the temperature gradient in the nonisothermal two-phase DM where water evaporates at the hotter catalyst layer, diffuses as a vapor form and then condenses on the cooler land region. Therefore, the vapor phase diffusion exacerbates DM flooding near the land region, while it alleviates DM flooding near the gas channel.

• Nuclear Engineering and Technology
• /
• v.52 no.4
• /
• pp.755-763
• /
• 2020

The critical flow phenomenon has been studied because of its significant effect for design basis accidents in nuclear power plants. Transition points from thermal non-equilibrium to equilibrium are different according to the geometric effect on the critical flow. This study evaluates the uncertainty parameters of the critical flow model for analysis of DBA (Design Basis Accident) with the MARS-KS (Multi-dimensional Analysis for Reactor Safety-KINS Standard) code used as an independent regulatory assessment. The uncertainty of the critical flow model is represented by three parameters including the thermal non-equilibrium factor, discharge coefficient, and length to diameter (L/D) ratio, and their ranges are determined using large-scale Marviken test data. The uncertainty range of the thermal non-equilibrium factor is updated by the MCDA (Model Calibration through Data Assimilation) method. The updated uncertainty range is confirmed using an LBLOCA (Large Break Loss of Coolant Accident) experiment in the LOFT (Loss of Fluid Test) facility. The uncertainty ranges are also used to calculate an LBLOCA of the APR (Advanced Power Reactor) 1400 NPP (Nuclear Power Plants), focusing on the effect of the PCT (Peak Cladding Temperature). The results reveal that break flow is strongly dependent on the degree of the thermal non-equilibrium state in a ruptured pipe with a small L/D ratio. Moreover, this study provides the method to handle the thermal non-equilibrium factor, discharge coefficient, and length to diameter (L/D) ratio in the system code.

• Journal of Korean Society on Water Environment
• /
• v.33 no.4
• /
• pp.449-459
• /
• 2017

In Geum River of Korea, three multi-purpose weirs were built at the downstream of Daecheong Reservoir during the Four Major River Restoration Project (FMRRP). The weirs have altered the hydraulic characteristics of the river, and consequently transformed the large areas of flowing ecosystem to deep and wide stagnant environment. In every summer, a thermal stratification occurred near the Baekje Weir having mean depth of 4.0 m, and the surface algal blooms dominated by buoyant cyanobacteria have been frequently formed after the FMRRP. The objective of this study was to investigate the relationship between flow velocity and thermal stability of the waterbody using a three-dimensional (3D) hydrodynamic model (EFDC+) after calibration against the thermistor chain data obtained in 2014. A new Sigma-Zed vertical grid system of EFDC+ that minimize the pressure gradient errors was used to better simulate the thermodynamics of the waterbody. The model reasonably simulated the vertical profiles of the observed water temperatures. The vertical mean flow velocity and the Richardson Number (Ri) that represents the stability of waterbody were estimated for various management water levels and flow rates scenarios. The results indicated that the thermal stability of the waterbody is mostly high ($Ri{\gg}0.25$) enough to establish stratification, and largely depend on the flow velocity. The critical flow velocity that can avoid a persistent thermal stratification was found to be approximately 0.1 m/s.

• Proceedings of the Korean Vacuum Society Conference
• /
• 2014.02a
• /
• pp.221-221
• /
• 2014

The present work proposes an improved numerical simulator for design and modification of large area capacitively coupled plasma (CCP) processing chamber. CCP, as notoriously well-known, demands the tremendously huge computational cost for carrying out transient analyses in realistic multi-dimensional models, because electron dissociations take place in a much smaller time scale (${\Delta}t{\approx}10-8{\sim}10-10$) than time scale of those happened between neutrals (${\Delta}t{\approx}10-1{\sim}10-3$), due to the rf drive frequencies of external electric field. And also, for spatial discretization of electron flux (Je), exponential scheme such as Scharfetter-Gummel method needs to be used in order to alleviate the numerical stiffness and resolve exponential change of spatial distribution of electron temperature (Te) and electron number density (Ne) in the vicinity of electrodes. Due to such computational intractability, it is prohibited to simulate CCP deposition in a three-dimension within acceptable calculation runtimes (<24 h). Under the situation where process conditions require thickness non-uniformity below 5%, however, detailed flow features of reactive gases induced from three-dimensional geometric effects such as gas distribution through the perforated plates (showerhead) should be considered. Without considering plasma chemistry, we therefore simulated flow, temperature and species fields in three-dimensional geometry first, and then, based on that data, boundary conditions of two-dimensional plasma discharge model are set. In the particular case of SiH4-NH3-N2-He CCP discharge to produce deposition of SiNxHy thin film, a cylindrical showerhead electrode reactor was studied by numerical modeling of mass, momentum and energy transports for charged particles in an axi-symmetric geometry. By solving transport equations of electron and radicals simultaneously, we observed that the way how source gases are consumed in the non-isothermal flow field and such consequences on active species production were outlined as playing the leading parts in the processes. As an example of application of the model for the prediction of the deposited thickness uniformity in a 300 mm wafer plasma processing chamber, the results were compared with the experimentally measured deposition profiles along the radius of the wafer varying inter-electrode gap. The simulation results were in good agreement with experimental data.