• Journal of Soil and Groundwater Environment
• /
• v.20 no.2
• /
• pp.41-46
• /
• 2015

The drag of the excess charge in an electrical double layer at the solid fluid interface due to water flow induces the streaming current, i.e., the streaming potential (SP). Here we introduce a sandbox experiment to study this hydroelectric coupling in case of a tracer test. An acrylic tank was filled up with homogeneous sand as a sand aquifer, and the upstream and downstream reservoirs were connected to the sand aquifer to control the hydraulic gradient. Under a steady-state water flow condition, a tracer test was performed in the sandbox with the help of peristaltic pump, and tracer samples were collected from the same interval of five screened wells in the sandbox. During the tracer test, SP signals resulting from the distribution of 20 nonpolarizable electrodes were measured at the top of the tank by a multichannel meter. The results showed that there were changes in the observed SP after injection of tracer, which indicated that the SP was likely to be related to the solute transport.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.22 no.11
• /
• pp.727-734
• /
• 2010

The present study numerically investigates the motion of a solid body suspended in the square enclosure with natural convection. A two-dimensional circular cylinder levitated thermally has been simulated by using thermal lattice Boltzmann method(TLBM) with the direct-forcing immersed boundary method. To deal with the ascending, falling or levitation of a circular cylinder in natural convection, the immersed boundary method is expanded and coupled with the TLBM. The circular cylinder is located at the bottom of a square enclosure with no restriction on the motion and freely migrates due to the Boussinesq approximation which is employed for the coupling between the flow and temperature fields. For different density ratio between the cylinder and the fluid, the motion characteristics of the circular cylinder for various Grashof numbers have been carried out. The Prandtl number is fixed as 0.7.

• Structural Engineering and Mechanics
• /
• v.59 no.1
• /
• pp.37-57
• /
• 2016

Reservoir geomechanics can play an important role in hydrocarbon recovery mechanism. In $CO_2$-EOR process, reservoir geomechanics analysis is concerned with the simultaneous study of fluid flow and the mechanical response of the reservoir under $CO_2$ injection. Accurate prediction of geomechanical effects during $CO_2$ injection will assist in modeling the Carbon dioxide recovery process and making a better design of process and production equipment. This paper deals with the implementation of a program (FORTRAN 90 interface code), which was developed to couple conventional reservoir (ECLIPSE) and geomechanical (ABAQUS) simulators, using a partial coupling algorithm. A geomechanics reservoir partially coupled approach is presented that allows to iteratively take the impact of geomechanics into account in the fluid flow calculations and therefore performs a better prediction of the process. The proposed approach is illustrated on a realistic field case. The reservoir geomechanics coupled models show that in the case of lower maximum bottom hole injection pressure, the cumulative oil production is more than other scenarios. Moreover at the high injection pressures, the production rates will not change with the injection bottom hole pressure variations. Also the FEM analysis of the reservoir showed that at $CO_2$ injection pressure of 11000 Psi the plastic strain has been occurred in the some parts of the reservoir and the related stress path show a critical behavior.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Journal of the Korean Chemical Society
• /
• v.42 no.6
• /
• pp.618-628
• /
• 1998

Pure CaO, Mn-doped CaO, Mn/CaO, and K/CaO catalysts were prepared and tested as catalysts for the oxidative coupling of methane in the temperature range of 600 to 800$^{\circ}C$ to investigate the effects of Mn- and K-addition on the catalytic activity of calcium oxide. To characterize the catalysts, X-ray powder diffraction(XRD), XPS, SEM, DSC, and TG analyses were performed. The catalytic reaction was carried out in a single-pass flow reactor using on-line gas chromatography system. Normalized reaction conditions were generally $p(CH_4)/p(O_2)=250$ Torr/50 Torr, total feed flow rate=30 mL/min, and 1 atm of total pressure with He being used as diluent gas. Among the catalysts tested, 6.3 mol% Mn-doped CaO catalyst showed the best $C_2$ yield of 8.0% with a selectivity of 43.2% at 775$^{\circ}C$. The $C_2$ selectivity increased on lightly doped CaO catalysts, while decreased on heavily doped CaO([Mn] > 6.3 mol%) catalysts. 6 wt.% Mn/CaO and 6 wt.% K/CaO catalysts showed the $C_2$ selectivities of 13.2% and 30.9%, respectively, for the reaction. Electrical conductivities of CaO and Mn-doped CaO were measured in the temperature range of 500 to 1000$^{\circ}C$ at Po2's of $10^{-3}\; to\;10^{-1}\;atm.$ The electrical conductivity was decreased with Mn-doping and increased with increasing $P0_2$in the range of $10^{-3}\;to\;10^{-1}\;atm,$ indicating the specimens to be p-type semiconductors. It was suggested that the interstitial oxygen ions formed near the surface can activate methane and the formation of interstitial oxygen ions was discussed on the basis of solid-state chemistry.