• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.499-504
• /
• 2001

A validation of a 3D CFD model for predicting local subcooling of moderator in the vicinity of calandria tubes in a CANDU reactor is performed. The small scale moderator experiments performed at Sheridan Park Experimental Laboratory(SPEL) in Ontario, Canada[1] is used for the validation. Also a comparison is made between previous CFD analyses based on 2DMOTH and PHOENICS, and the current model analysis for the same SPEL experiment. For the current model, a set of grid structures for the same geometry as the experimental test section is generated and the momentum, heat and continuity equations are solved by CFX-4.3, a CFD code developed by AEA technology. The matrix of calandria tubes is simplified by the porous media approach. The standard $k-\varepsilon$ turbulence model associated with logarithmic wall treatment and SIMPLEC algorithm on the body fitted grid are used and buoyancy effects are accounted for by the Boussinesq approximation. For the test conditions simulated in this study, the flow pattern identified is a buoyancy-dominated flow, which is generated by the interaction between the dominant buoyancy force by heating and inertial momentum forces by the inlet jets. As a result, the current CFD moderator analysis model predicts the moderator temperature reasonably, and the maximum error against the experimental data is kept at less than $2.0^{\circ}C$ over the whole domain. The simulated velocity field matches with the visualization of SPEL experiments quite well.

• 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.

• Proceedings of the Korean Radioactive Waste Society Conference
• /
• 2005.06a
• /
• pp.231-238
• /
• 2005

On the basis of flow resistance theory the conceptual model and related mathematical descriptions is proposed for resistance modeling of groundwater flow in CPM(continuum Porous medium), DFN(discrete fracture network) and fractured-porous medium. The proposed model is developed on the basis of finite volume method assuming steady-state, constant density groundwater flow. The basic approach of the method is to evaluate inter-block flow resistance values for a staggered grid arrangement, i.e. fluxes are stored at cell walls and scalars at cell centers. The balance of forces, i.e. the Darcy law, is utilized for each control volume centered around the point where the velocity component is stored. The transmissivity (or permeability) at the interface is assumed to be the harmonic average of neighboring blocks. Flow resistance theory was utilized to relate the fluxes between the grid blocks with residual pressures. The flow within porous medium is described by three dimensional equations and that within an individual fracture is described by a two dimensional equivalent of the flow equations for a porous medium. Newly proposed models would contribute to develop flow simulation techniques with various matrix characteristics.

• Geomechanics and Engineering
• /
• v.23 no.4
• /
• pp.327-338
• /
• 2020

The role of precise prediction of subsurface fluids and discrimination among them cannot be ignored in reservoir characterization and petroleum prospecting. A suitable rock physics model should be build for the extraction of valuable information form seismic data. The main intent of current work is to present a rock physics model to analyze the characteristics of seismic wave propagating through a cracked porous rock saturated by a three phase fluid. Furthermore, the influence on wave characteristics due to variation in saturation of water, oil and gas were also analyzed for oil and water as wet cases. With this approach the objective to explore wave attenuation and dispersion due to wave induce fluid flow (WIFF) at seismic and sub-seismic frequencies can be precisely achieved. We accomplished our proposed approach by using BISQ equations and by applying appropriate boundary conditions to incorporate heterogeneity due to saturation of three immiscible fluids forming a layered system. To authenticate the proposed methodology, we compared our results with White's mesoscopic theory and with the results obtained by using Biot's poroelastic relations. The outcomes reveals that, at low frequencies seismic wave characteristics are in good agreement with White's mesoscopic theory, however a slight increase in attenuation at seismic frequencies is because of the squirt flow. Moreover, our work crop up as a practical tool for the development of rock physical theories with the intention to identify and estimate properties of different fluids from seismic data.

• Journal of computational fluids engineering
• /
• v.17 no.4
• /
• pp.41-48
• /
• 2012

A numerical analysis of thermal stratification in the upper plenum of the MONJU fast breeder reactor was performed. Calculations were performed for a 1/6 simplified model of the MONJU reactor using the commercial code, CFX-13. To better resolve the geometrically complex upper core structure of the MONJU reactor, the porous media approach was adopted for the simulation. First, a steady state solution was obtained and the transient solutions were then obtained for the turbine trip test conducted in December 1995. The time dependent inlet conditions for the mass flow rate and temperature were provided by JAEA. Good agreement with the experimental data was observed for steady state solution. The numerical solution of the transient analysis shows the formation of thermal stratification within the upper plenum of the reactor vessel during the turbine trip test. The temporal variations of temperature were predicted accurately by the present method in the initial rapid coastdown period (~300 seconds). However, transient numerical solutions show a faster thermal mixing than that observed in the experiment after the initial coastdown period. A nearly homogenization of the temperature field in the upper plenum is predicted after about 900 seconds, which is a much shorter-term thermal stratification than the experimental data indicates. This discrepancy is due to the shortcoming of the turbulence models available in the CFX-13 code for a natural convection flow with thermal stratification.

• Structural Engineering and Mechanics
• /
• v.66 no.4
• /
• pp.505-513
• /
• 2018

Steam flooding is widely used in heavy oil reservoir with coupling effects among the formation temperature change, fluid flow and solid deformation. The effective stress, porosity and permeability in this process can be affected by the multi-physical coupling of thermal, hydraulic and mechanical processes (THM), resulting in a complex interaction of geomechanical effects and multiphase flow in the porous media. Quantification of the state of deformation and stress in the reservoir is therefore essential for the correct prediction of reservoir efficiency and productivity. This paper presents a coupled fluid flow, thermal and geomechanical model employing a program (MATLAB interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators for coupled THM processes in multiphase reservoir modeling. In each simulation cycle, time dependent reservoir pressure and temperature 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 and new porosity and permeability are obtained using volumetric strains for the next analysis step. Finally, the proposed approach is illustrated on a complex coupled problem related to steam flooding in an oil reservoir. The reservoir coupled study showed that permeability and porosity increase during the injection scenario and increasing rate around injection wells exceed those of other similar comparable cases. Also, during injection, the uplift occurred very fast just above the injection wells resulting in plastic deformation.

• Tunnel and Underground Space
• /
• v.32 no.2
• /
• pp.144-159
• /
• 2022

The present study developed a sequential approach-based numerical simulator for modeling coupled thermal-hydrological-mechanical (THM) processes in the ground and investigated the computational performance of the coupling analysis algorithm. The present sequential approach linked the two different solvers: an open-source numerical code, OpenGeoSys for solving the thermal and hydrological processes in porous media and a commercial code, FLAC3D for solving the geomechanical response of the ground. A benchmark test of the developed simulator was carried out using a THM problem where an analytical solution is given. The benchmark problem involves the coupled behavior (variations in temperature, pore pressure, stress, and deformation with time) of a fully saturated porous medium which is subject to a point heat source. The results of the analytical solution and numerical simulation were compared and the validity of the numerical simulator was investigated.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2001.10a
• /
• pp.181-185
• /
• 2001

Three-dimensional finite-element-based numerical model of turbulent flow, heat transfer, macroscopic solidification and inclusion trajectory in a continuos steel slab caster was developed Turbulence was incorporated using the Improved Low-Re turbulence model with positive preserving approach. The mushy region was modeled as the porous media with average effective viscosity. A series of simulations was carried out to investigate the effects of the casting speed, the slab size, the delivered superheat the immersion depth of the SEN on the transport phenomena. In the absence of any known experimental data related to velocity profiles, the numerical predictions of the solidified profile on a caster was compared with breakouts data and a good agreement was found.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.10a
• /
• pp.106-111
• /
• 2001

This paper presents the numerical methodology of ATHOS3 code 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. The ATHOS3 code is applied to the thermal hydraulic analysis of steam generator in the Korea Kori Unit-1 nuclear power plant and the computed results are presented.

• 한국신재생에너지학회:학술대회논문집
• /
• 2006.11a
• /
• pp.340-343
• /
• 2006

In this paper, heat and mass transfer characteristics through experimental and numerical study are extensively investigated in steam reform ins reactor under given operating conditions. In order to get simulated data at outlet of the reformer, heterogeneous reactor model is incorporated. As the reaction also takes place in porous media, two medium approach is used to take into account thermally non-equilibrium phenomena between catalyst and bulk gas. In steam reforming reaction, heat transfer issue is so significant that geometrical configuration study is also conducted.