• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.43-59
• /
• 2018

The presence of the pore fluid strongly influences the reponse of the soil subjected to external loading and in many cases increases the risk of final failure. In this paper, we propose the use of a discrete beam lattice model with the aim to investigate the coupling effects of the solid and fluid phase on the response and failure mechanisms in the saturated soil. The discrete cohesive link lattice model used in this paper, is based on inelastic Timoshenko beam finite elements with enhanced kinematics in axial and transverse direction. The coupling equations for the soil-pore fluid interaction are derived from Terzaghi's principle of effective stresses, Biot's porous media theory and Darcy's law for fluid flow through porous media. The application of the model in soil mechanics is illustrated through several numerical simulations.

• Coupled systems mechanics
• /
• v.7 no.6
• /
• pp.649-668
• /
• 2018

In this paper, we present a numerical model for fluid-structure interaction between structure built of porous media and acoustic fluid, which provides both pore pressure inside porous media and hydrodynamic pressures and hydrodynamic forces exerted on the upstream face of the structure in an unified manner and simplifies fluid-structure interaction problems. The first original feature of the proposed model concerns the structure built of saturated porous medium whose response is obtained with coupled discrete beam lattice model, which is based on Voronoi cell representation with cohesive links as linear elastic Timoshenko beam finite elements. The motion of the pore fluid is governed by Darcy's law, and the coupling between the solid phase and the pore fluid is introduced in the model through Biot's porous media theory. The pore pressure field is discretized with CST (Constant Strain Triangle) finite elements, which coincide with Delaunay triangles. By exploiting Hammer quadrature rule for numerical integration on CST elements, and duality property between Voronoi diagram and Delaunay triangulation, the numerical implementation of the coupling results with an additional pore pressure degree of freedom placed at each node of a Timoshenko beam finite element. The second original point of the model concerns the motion of the outside fluid which is modeled with mixed displacement/pressure based formulation. The chosen finite element representations of the structure response and the outside fluid motion ensures for the structure and fluid finite elements to be connected directly at the common nodes at the fluid-structure interface, because they share both the displacement and the pressure degrees of freedom. Numerical simulations presented in this paper show an excellent agreement between the numerically obtained results and the analytical solutions.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.06a
• /
• pp.309-313
• /
• 2000

A structural-acoustic coupling problem involving fluid in a cavity divided with flexible walls and porous materials is investigated in this paper. In many practical problems, the use of finite elements to discretize the fluid region leads to large stiffness and mass matrices. But, since the acoustic boundary element discretization requires to put elements only on the surface of structure, the size of matrices is reduced considerably. Here, we developed a numerical analysis program for the structural-acoustic coupling problems of the multi-region cavity, using boundary elements for the fluid regions and finite elements for the structure. By considering sound transmission through layered systems placed in a cavity, the accuracy of the coupled acoustical-structural finite element model has been verified by comparing its transmission loss predictions with analytical sloutions. Example problems are included to investigate the characteristics of the multi-region structural-acoustic coupling system with porous material.

• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.83-101
• /
• 2008

A derivation is given of two-scale models that are able to describe deformation and flow in a fluid-saturated and progressively fracturing porous medium. From the micromechanics of the flow in the cavity, identities are derived that couple the local momentum and the mass balances to the governing equations for a fluid-saturated porous medium, which are assumed to hold on the macroscopic scale. By exploiting the partition-of-unity property of the finite element shape functions, the position and direction of the fractures are independent from the underlying discretization. The finite element equations are derived for this two-scale approach and integrated over time. The resulting discrete equations are nonlinear due to the cohesive crack model and the nonlinearity of the coupling terms. A consistent linearization is given for use within a Newton-Raphson iterative procedure. Finally, examples are given to show the versatility and the efficiency of the approach.

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

• Journal of the Korean GEO-environmental Society
• /
• v.4 no.1
• /
• pp.5-10
• /
• 2003

The solids and the fluids in porous media have a relative velocity to each other. Due to physically and chemically different material properties and their relative velocity, the behavior of saturated porous media is extremely complicated. Thus, in order to describe and clarify the deformation behavior of saturated porous media, constitutive models for deformation of porous media coupling several effects need to be developed in frame of Arbitrary Lagrangian Eulerian(ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian elements, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media which are considered for the behavior of the solids and the fluids. For this reason, mass balance equations for saturated porous media are derived here in ALE description frames. ALE formulations of mass conservation for the solid phase and the fluid phase are expressed. Then, linear momentum balance equation for porous media as multiphase media is expressed.

• Journal of the Korean Geotechnical Society
• /
• v.24 no.2
• /
• pp.67-75
• /
• 2008

It is interesting to note that distinct element method has been used extensively to model the response of micro and discontinuous behavior in geomechanics. Impressive advances related to response of distinct particles have been conducted and there were difficulties in considering fluid effect simultaneously. Current distinct element methods are progressively developed to solve particle-fluid coupling focused on fluid flow through soil, rock or porous medium. In this research, numerical simulations of fluid injection into particulate materials were conducted to observe cavity initiation and propagation using distinct element method. After generation of initial particles and wall elements, confining stress was applied by servo-control method. The fluid scheme solves the continuity and Navior-Stokes equations numerically, then derives pressure and velocity vectors for fixed grid by considering the existence of particles within the fluid cell. Fluid was injected as 7-step into the assembly in the x-direction from the inlet located at the center of the left boundary under confining stress condition, $0.1MP{\alpha}\;and\;0.5MP{\alpha}$, respectively. For each simulation, movement of particles, flow rate, fluid velocity, pressure history, wall stress including cavity initiation and propagation by interaction of flulid-paricles were analyzed.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.125-145
• /
• 2020

In this paper, we present a 3D thermo-hydro-mechanical coupled discrete beam lattice model of structure built of the nonisothermal saturated poro-plastic medium subjected to mechanical loads and nonstationary heat transfer conditions. The proposed model is based on Voronoi cell representation of the domain with cohesive links represented as inelastic Timoshenko beam finite elements enhanced with additional kinematics in terms of embedded strong discontinuities in axial and both transverse directions. The enhanced Timoshenko beam finite element is capable of modeling crack formation in mode I, mode II and mode III. Mode I relates to crack opening, mode II relates to in-plane crack sliding, and mode III relates to the out-of-plane shear sliding. The pore fluid flow and heat flow in the proposed model are governed by Darcy's law and Fourier's law for heat conduction, respectively. The pore pressure field and temperature field are approximated with linear tetrahedral finite elements. By exploiting nodal point quadrature rule for numerical integration on tetrahedral finite elements and duality property between Voronoi diagram and Delaunay tetrahedralization, the numerical implementation of the coupling results with additional pore pressure and temperature degrees of freedom placed at each node of a Timoshenko beam finite element. The results of several numerical simulations are presented and discussed.

• Tunnel and Underground Space
• /
• v.34 no.2
• /
• pp.127-142
• /
• 2024

The robustness of a numerical method means that its computational performance is maintained under various modeling conditions. New numerical methods or codes need to be assessed for robustness through benchmark testing. The TOUGH-FLAC modeling approach has been applied to various fields such as subsurface carbon dioxide storage, geological disposal of spent nuclear fuel, and geothermal development both domestically and internationally, and the modeling validity has been examined by comparing the results with experimental measurements and other numerical codes. In the present study, a benchmark test of the TOUGH-FLAC approach was performed based on a coupled thermal-hydro-mechanical behavior problem with an analytical solution. The analytical solution is related to the temperature, pore water pressure, and mechanical behavior of a fully saturated porous medium that is subjected to a point heat source. The robustness of the TOUGH-FLAC approach was evaluated by comparing the analytical solution with the results of numerical simulation. Additionally, the effects of thermal-hydro-mechanical coupling terms, fluid phase change, and timestep on the computation of coupled behavior were investigated.

• Journal of the Korean Geotechnical Society
• /
• v.29 no.8
• /
• pp.65-73
• /
• 2013

Ground-coupled heat pump system has attracted attention as a promising renewable energy technology due to its improving energy efficiency and eco-friendly mechanism for space cooling and heating. Pipes buried in the ground play a role of direct thermal interaction between circulating fluid inside the pipe and surrounding soils in the geothermal exchange system. However, both complexities of turbulent flow coupling thermal-hydraulic phenomena and very long aspect ratio of the pipe make it difficult to model the heat exchange system directly. Energy balance for fluid flow inside the pipe was derived to model thermal-hydraulic phenomena, and one-dimensional pipe element was proposed through Galerkin formation and time integration of the equation. Developed element is combined to pre-developed FEM code for THM phenomena in porous media. Numerical results of Thermal Response Test showed that line-source model overestimates equivalent thermal conductivity of surrounding soils due to thermal interaction between adjacent pipes and finite length of the pipe. Thus, inverse analysis for the TRT simulation was conducted to present optimal transformation matrix with utmost convergence.