Eduardo T Lima Junior¹, Wilson S Venturini², Ahmed Benallal³

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.3&4, pp. 169-184, 2011, DOI:10.3970/cmes.2011.078.169

Abstract Some considerations on the numerical analysis of brittle rocks are presented in this paper. The rock is taken as a poro-elastic domain, in full-saturated condition, based on the Biot's Theory. The solid matrix of this porous medium is considered to be susceptible to isotropic damage occurrence. An implicit boundary element method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and two-dimensional elastostatics problems. The integration over boundary elements is evaluated by using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the… More >

P. Maghoul¹, B. Gatmiri^1,2, D. Duhamel¹

CMES-Computer Modeling in Engineering & Sciences, Vol.78, No.1, pp. 51-76, 2011, DOI:10.3970/cmes.2011.078.051

Abstract This paper aims at obtaining boundary integral formulations as well as three dimensional(3D) fundamental solutions for unsaturated soils under dynamic loadings for the first time. The boundary integral equations are derived via the use of the weighted residuals method in a way that permits an easy discretization and implementation in a Boundary Element code. Also, the associated 3D fundamental solutions for such deformable porous medium are derived in Laplace transform domain using the method of Hérmander. The derived results are verified analytically by comparison with the previously introduced corresponding fundamental solutions in elastodynamic limiting case. These solutions can be used,… More >

A.B. Tatomir^1,2, A.Szymkiewicz³, H. Class¹, R. Helmig¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 81-112, 2011, DOI:10.3970/cmes.2011.077.081

Abstract We present a two phase flow conceptual model, the corresponding simulator (2pMINC) and a workflow for large-scale fractured reservoirs, based on a continuum fracture approach which uses the multiple interacting continua (MINC) method complemented with an improved upscaling technique. The complex transient behavior of the flow processes in fractured porous media is captured by subgridding the coarse blocks in nested volume elements which have effective properties calculated from the detailed representation of the fracture system. In this way, we keep a physically based approach, preserve the accuracy of the model, avoid the common use of empirically derived transfer functions and… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.3, pp. 227-248, 2010, DOI:10.3970/cmes.2010.066.227

Abstract In this work, a meshless method based on the local Petrov-Galerkin approach is proposed for the solution of pore-dynamic problems considering elastic and elastoplastic materials. Formulations adopting the Heaviside step function as the test functions in the local weak form are considered. The moving least-square method is used for the approximation of physical quantities in the local integral equations. After spatial discretization is carried out, a nonlinear system of time-domain ordinary differential equations is obtained. This system is solved by Newmark/Newton-Raphson techniques. The present work is based on the u-p formulation and the incognita fields of the coupled analysis in… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.61, No.2, pp. 177-200, 2010, DOI:10.3970/cmes.2010.061.177

Abstract In this work, meshless methods based on the local Petrov-Galerkin approach are employed for the time-domain dynamic analysis of porous media. For the spatial discretization of the pore-dynamic model, MLPG formulations adopting Gaussian weight functions as test functions are considered, as well as the moving least square method is used to approximate the incognita fields. For time discretization, the generalized Newmark method is adopted. The present work is based on the u-p formulation and the incognita fields of the coupled analysis in focus are the solid skeleton displacements and the interstitial fluid pore pressures. Independent spatial discretization is considered for… More >

Ehsan Jabbari¹, Behrouz Gatmiri^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.1, pp. 31-44, 2007, DOI:10.3970/cmes.2007.018.031

Abstract In this paper after a discussion about the evolution of the unsaturated soils' governing differential equations and a brief history of the Green's functions for porous media, the governing equations, i.e., the mathematical model in the presence of heat effects are presented and simplified so as the derivation of the associated Green's functions be in the realm of possibility. The thermal two- and three-dimensional, full- and half-space Green's functions for unsaturated porous media, although in a relatively simplified form, are being introduced for the first time, following the previous works of the authors. The derived Green's functions have been demonstrated… More >

M.R. Arab^1,2, E. Semma³, B. Pateyron¹, M. El Ganaoui¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.2, pp. 161-176, 2009, DOI:10.3970/fdmp.2009.005.161

Abstract In this work, flows in porous media are simulated by using a Lattice Boltzmann Method (LBM). A model D2Q9 with a single collision operator is proposed. This method is applied on 2D digital images obtained by a Scanning Electron Microscope technique (SEM), and followed by a special treatment in order to obtain an image of synthesis that is finally read by the numerical code. The first results tested on two-dimensional configurations show the reliability of this strategy in simulating with a good accuracy phenomena of heat and mass transport. The numerical study is extended to the prediction of physical parameters… More >

Aouachria Z¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.2, pp. 137-148, 2009, DOI:10.3970/fdmp.2009.005.137

Abstract This work treats heat and mass transfer by natural convection along a vertical wall in porous media imbibed by fluid, using an integral method. The problem governing parameters are the buoyancy ratio,N, and the Lewis number, Le The results for the local Nusselt and Sherwood numbers are presented for a large range of these parameters. The concentration and thermal boundary layer thickness are also determined. We observe that our results are in good agreement with those obtained by Bejan and Khair (1985). More >

L. Bennamoun¹, A. Belhamri²

FDMP-Fluid Dynamics & Materials Processing, Vol.4, No.4, pp. 221-230, 2008, DOI:10.3970/fdmp.2008.004.221

Abstract This work focuses on tyipical heat and mass transfer phenomena during the processing of products in the context of the packed-bed drying method (products arranged in thick layers into dryers working in forced convection mode). The dryers are modeled as porous media at the macroscopic level. The simulations are carried out using the mass, momentum (written in the framework of the Darcy's law approximation) and energy equations applied for the different components. A diffusion model based on Fick's law is also used to take into account the drying kinetics. This approach allows monitoring of the variations of humidity and temperature… More >

D. Srinivasacharya¹, G. Swamy Reddy¹

FDMP-Fluid Dynamics & Materials Processing, Vol.9, No.3, pp. 291-305, 2013, DOI:10.3970/fdmp.2013.009.291

Abstract Free convection and related heat and mass transfer along a vertical plate embedded in a power-law fluid saturated Darcy porous medium with thermal and solutal stratification effects is studied. The governing partial differential equations are transformed into ordinary differential equations using similarity transformations and then solved numerically by means of a shooting method. The variations of non-dimensional velocity, temperature and concentration are presented graphically for various values of the power-law index, and of the thermal and solutal stratification parameters. In addition, the heat and mass transfer rates are tabulated for different values of the governing nondimensional numbers. More >