@article {658, title = {A methodology to calibrate and to validate effective solid potentials of heterogeneous porous media from computed tomography scans and laboratory-measured nanoindentation data}, journal = {Acta Geotechnica}, volume = {13}, year = {2018}, month = {Jan-12-2018}, pages = {1369 - 1394}, abstract = {

Built on the framework of effective interaction potentials using lattice element method, a methodology to calibrate and to validate the elasticity of solid constituents in heterogeneous porous media from experimentally measured nanoindentation moduli and imported scans from advanced imaging techniques is presented. Applied to computed tomography (CT) scans of two organic-rich shales, spatial variations of effective interaction potentials prove instrumental in capturing the effective elastic behavior of highly heterogeneous materials via the first two cumulants of experimentally measured distributions of nanoindentation moduli. After calibration and validation steps while implicitly accounting for mesoscale texture effects via CT scans, Biot poroelastic coefficients are simulated. Analysis of stress percolation suggests contrasting pathways for load transmission, a reflection of microtextural differences in the studied cases. This methodology to calibrate elastic energy content of real materials from advanced imaging techniques and experimental measurements paves the way to study other phenomena such as wave propagation and fracture while providing a platform to fine-tune effective behavior of materials given advancements in additive manufacturing and machine learning algorithms .

}, issn = {1861-1125}, doi = {10.1007/s11440-018-0687-9}, url = {http://link.springer.com/10.1007/s11440-018-0687-9}, author = {Monfared, Siavash and Hadrien Laubie and Farhang Radja{\"\i} and Hubler, Mija H. and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {239, title = {Disorder-induced stiffness degradation of highly disordered porous materials}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {106}, year = {2017}, month = {Sep-2017}, chapter = {207-228}, abstract = {

The effective mechanical behavior of multiphase solid materials is generally modeled by means of homogenization techniques that account for phase volume fractions and elastic moduli without considering the spatial distribution of the different phases. By means of extensive numerical simulations of randomly generated porous materials using the lattice element method, the role of local textural properties on the effective elastic properties of disordered porous materials is investigated and compared with different continuum micromechanics-based models. It is found that the pronounced disorder-induced stiffness degradation originates from stress concentrations around pore clusters in highly disordered porous materials. We identify a single disorder parameter, \φsa, which combines a measure of the spatial disorder of pores (the clustering index, sa) with the pore volume fraction (the porosity, \φ) to scale the disorder-induced stiffness degradation. Thus, we conclude that the classical continuum micromechanics models with one spherical pore phase, due to their underlying homogeneity assumption fall short of addressing the clustering effect, unless additional texture information is introduced, e.g. in form of the shift of the percolation threshold with disorder, or other functional relations between volume fractions and spatial disorder; as illustrated herein for a differential scheme model representative of a two-phase (solid\–pore) composite model material.

Fig. 1. Two-dimensional porous media: (a) ordered system, (b) type 1 disorder (pores{\textellipsis}

Fig. 2. Two-point probability functions and associated microstructures of four systems{\textellipsis}

Fig. 3. Probability Density Function (PDF) of the local porosity (φa) and associated{\textellipsis}

Fig. 4. (a) Degrees of freedom of a link element between points i and j, (b) D3Q18 unit{\textellipsis}

Fig. 5. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity{\textellipsis}

Fig. 6. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity

Fig. 7. Ordered system under uniaxial strain, (a) schematic stress map, the darker the{\textellipsis}

Fig. 8. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of the porosity

Fig. 9. Probability Density Function (PDF) of the normalized stress (σxx/〈σxx〉) in the{\textellipsis}

Fig. 10. Dimensionless disordered Young{\textquoteright}s modulus: F(φsa)=Eeff(φ,sa)/Eordered(φ) as{\textellipsis}

Fig. 11. Two-phase periodic porous solid

Fig. 12. Dimensionless disordered Young{\textquoteright}s modulus: FMT(φsa)=Eeff(φ,sa)/EMT(φ) as{\textellipsis}

Fig. 13. Integration path for the lower bound differential scheme

Fig. 14. Dimensionless effective Young{\textquoteright}s modulus: Eeff(φ)/Es as function of porosity the

Fig. 15. Dimensionless disordered Young{\textquoteright}s modulus: FMT(φsa)=Eeff(φ,sa)/EMT(φ) as{\textellipsis}

Fig. B.16. (a) Porosity as a function of the number of pores

Fig. D.17. High porosity limit geometry

Fig. E.18. Calculation of the local porosity in a square observation window

}, doi = {10.1016/j.jmps.2017.05.008}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {591, title = {Effective Potentials and Elastic Properties in the Lattice-Element Method: Isotropy and Transverse Isotropy}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Sep-2017}, pages = { Article Number: UNSP 04017007 }, abstract = {

Lattice approaches have emerged as a powerful tool to capture the effective mechanical behavior of heterogeneous materials using harmonic interactions inspired from beam-type stretch and rotational interactions between a discrete number of mass points. In this paper, the lattice element method (LEM) is reformulated within the conceptual framework of empirical force fields employed at the lattice scale. Within this framework, because classical harmonic formulations are but a Taylor expansion of nonharmonic potential expressions, they can be used to model both the linear and the nonlinear response of discretized material systems. Specifically, closed-form calibration procedures for such interaction potentials are derived for both the isotropic and the transverse isotropic elastic cases on cubic lattices, in the form of linear relations between effective elasticity properties and energy parameters that define the interactions. The relevance of the approach is shown by an application to the classical Griffith crack problem. In particular, it is shown that continuum-scale quantities of linear-elastic fracture mechanics, such as stress intensity factors (SIFs), are well captured by the method, which by its very discrete nature removes geometric discontinuities that provoke stress singularities in the continuum case. With its strengths and limitations thus defined, the proposed LEM is well suited for the study of multiphase materials whose microtextural information is obtained by, e.g., X-ray micro-computed tomography. (c) 2017 American Society of Civil Engineers.

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000125}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000125}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {586, title = {Mesoscale Poroelasticity of Heterogeneous Media}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Oct-19-2017}, pages = {Article Number: UNSP 04017016 }, abstract = {

The poromechanics of heterogeneous media is reformulated in a discrete framework using the lattice element method (LEM) that accounts for the presence of interfaces as well as local microtextural and elastic variations. The exchange of mechanical information between pore and solid(s) is captured by means of force field potentials for these domains, which eliminate the requirement of scale separability of continuum-based poromechanics approaches. In congruence with mu VT and NPT ensembles of statistical mechanics, discrete expressions for Biot poroelastic coefficients are derived. Considering harmonic-type interaction potentials for each link, analytical expressions for both isotropic and transversely isotropic effective elasticity are presented. The theory is validated against continuum-based expressions of Biot poroelastic coefficients for porous media with isotropic and transversely isotropic elastic solid behavior. (C) 2017 American Society of Civil Engineers.

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000136}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000136}, author = {Monfared, Siavash and Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {325, title = {A molecular informed poroelastic model for organic-rich, naturally occurring porous geocomposites}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {88}, year = {2016}, month = {Mar-2016}, pages = {186 - 203}, abstract = {

Molecular simulation results on organic maturity (mature and immature kerogen as the two asymptotic cases) are introduced into a continuum micromechanics based model for organic-rich shales. Through a fundamental functional relationship that constrains microporous kerogen density and elasticity variable spaces and within the framework of effective media theory; the model bridges the gap between asymptotic cases of organic maturity with texture as the overriding theme, specifically a matrix/inclusion (Mori\–Tanaka) texture for immature systems and a granular (self-consistent) texture for mature ones. The utility of the molecular results merged into a continuum framework is demonstrated by estimating kerogen\&$\#$39;s microporosity () from nanoindentation measurements. The effect of burial and diagenetic processes on the effective poroelasticity of these porous, naturally occurring geocomposites are captured by introduction of imperfect interfaces. Finally, the performance of the model is fully characterized by ranking the normalized contribution of uncertainty of input to the overall behavior and parameters of interest to geophysicists and geomechanicians such as degree of anisotropy and in situ stresses.

Fig.1. Clay type distribution in the studied samples

Fig.2. Schematics for the multi-scale maturity dependent model for organic-rich shales

Fig.3. A Scanning Electron Microscope (SEM) image on Haynesville sample

Fig.4. Calibration quality check by comparing predicted macroscopic values (drained{\textellipsis}

Fig.5. Validation by means of nanoindentation for Woodford on the left and Haynesville{\textellipsis}

Fig.6. Bulk modulus of kerogen against density for two of the samples studied by{\textellipsis}

}, issn = {00225096}, doi = {10.1016/j.jmps.2015.12.006}, author = {Monfared, Siavash and Franz-Josef Ulm} }