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

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.

• Korean Chemical Engineering Research
• /
• v.53 no.6
• /
• pp.814-817
• /
• 2015

A new ionic mass transfer correlation is derived for the fluid-saturated, horizontal porous layer. Darcy-Forchheimer model is used to explain characteristics of fluid motion. Based on the microscales of turbulence a backbone mass transfer relation is derived as a function of the Darcy-Rayleigh number, $Ra_D$ and the porous medium Schmidt number, $Sc_p$. For the Darcy's limit of $Sc_p{\gg}Ra_D$, the Sherwood number, Sh is a function of $Ra_D$ only. However, for the region of high $Ra_D$, Sh can be related with $Ra_DSc_p$. Based on the present backbone equation and the electrochemical mass transfer experiments which are electro plating or electroless plating, the new ionic mass transfer correlation is suggested in the porous media.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2005.03a
• /
• pp.47-51
• /
• 2005

If a structure is founded on the ground saturated with pore water, then the ground should be modeled as a saturated two-phase porous medium for accurate earthquake response analysis. In this study, an axisymmetric transmitting boundary hyperelement is developed for modeling of far field of the ground using u-U formulation for water-saturated transversely isotropic layered stratum. The developed hyperelement is verified by comparing the dynamic stiffness of rigid circular foundation on water-saturated isotropic layered stratum with the case of using equivalent single-phase medium model.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2006.03a
• /
• pp.345-350
• /
• 2006

If a structure is founded on the ground saturated with pore water, then the ground should be modeled as a saturated two-phase porous medium for accurate earthquake response analysis. In this study, a 3-dimensional transmitting boundary is developed for modeling of far field using u-U formulation for water-saturated transversely isotropic layered stratum. The developed transmitting boundary is verified by comparing the dynamic stiffness of rigid square foundation on water-saturated isotropic layered stratum with the case of using equivalent single-phase medium model.

• Korean Chemical Engineering Research
• /
• v.57 no.5
• /
• pp.723-729
• /
• 2019

A theoretical analysis was conducted of convective instability driven by buoyancy forces under transient temperature fields in an annular porous medium bounded by coaxial vertical cylinders. Darcy's law and Boussinesq approximation are used to explain the characteristics of fluid motion and linear stability theory is employed to predict the onset of buoyancy-driven motion. The linear stability equations are derived in a global domain, and then cast into in a self-similar domain. Using a spectral expansion method, the stability equations are reformed as a system of ordinary differential equations and solved analytically and numerically. The critical Darcy-Rayleigh number is founded as a function of the radius ratio. Also, the onset time and corresponding wavelength are obtained for the various cases. The critical time becomes smaller with increasing the Darcy-Rayleigh number and follows the asymptotic relation derived in the infinite horizontal porous layer.

• Geotechnical Engineering
• /
• v.13 no.1
• /
• pp.5-18
• /
• 1997

In order to make reliable ground shock predictions in saturated geological media, it is necessary to use multi -phase material models and numerical codes. This paper presents the results of theoretical study of the fundamental behavior of multi-phase porous media subjected to high dynanlic loadings, and deals with the development of numerical code MPDAP with JWL(Jones-Wilkins-Lee) model, which is capable of considering the kinds and characters of explosives. To check the global equilhorium equations of the numerical code, we carried out some verifications. In the cases of the elastic spherical wave propagation in a single phase medium, one-dimensional linear ronsolidation, and one timensional wave propagation in saturated linear elastic soils and rocks, the results calculated by MPDAP show close agreement with closed-form solutions or numerical solutions generated with two phase code.

• Geomechanics and Engineering
• /
• v.15 no.2
• /
• pp.793-803
• /
• 2018

Based on the theory of porous media, an interaction system of a floating pile and a saturated soil in cylindrical coordinates subjected to vertical harmonic load is presented in this paper. The surrounding soil is separated into two distinct layers. The upper soil layer above the level of pile base is described as a saturated viscoelastic medium and the lower soil layer is idealized as equivalent spring-dashpot elements with complex stiffness. Considering the cylindrically symmetry and the pile-soil compatibility condition of the interaction system, a frequency-domain analytical solution for dynamic impedance of the floating pile embedded in saturated viscoelastic soil is also derived, and reduced to verify it with existing solutions. An extensive parametric analysis has been conducted to reveal the effects of the impedance of the lower soil base, the interaction coefficient and the damping coefficient of the saturated viscoelastic soil layer on the vertical vibration of the pile-soil interaction system. It is shown that the vertical dynamic impedance of the floating pile significantly depends on the real stiffness of the impedance of the lower soil base, but is less sensitive to its dynamic damping variation; the behavior of the pile in poro-visco-elastic soils is totally different with that in single-phase elastic soils due to the existence of pore liquid; the effect of the interaction coefficient of solid and liquid on the pile-soil system is limited.

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

• Journal of the Earthquake Engineering Society of Korea
• /
• v.9 no.4 s.44
• /
• pp.67-74
• /
• 2005

This paper presents a new infinite element for modeling far-field region in dynamic analysis of a fluid-saturated two-phase medium. The infinite element method combined to the infinite element method has been effectively applied to several engineering problems where the full space or half-space medium should be modeled. However, the currently available infinite element for dynamic analysis of two-phase porous medium has a limitation that Pl and P2 waves can only be Included in shape function expressing behavior ol the body. In this paper, the infinite element method is extended to simulate arbitrary number of multi-component waves. For this purpose, the far-field of the porous medium is assumed to be a layered half-space, while the near-field Includes structures as well as irregular soil medium. The accuracy and effectiveness of the proposed element have demonstrated using 1-D and 2-D wave propagation problems.