Effective Potentials and Elastic Properties in the Lattice-Element Method: Isotropy and Transverse Isotropy

International audienceLattice 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

Mesoscale Poroelasticity of Heterogeneous Media

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 μ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

Effect of Confinement on Capillary Phase Transition in Granular Aggregates

Using a 3D mean-field lattice-gas model, we analyze the effect of confinement on the nature of capillary phase transition in granular aggregates with varying disorder and their inverse porous structures obtained by interchanging particles and pores. Surprisingly, the confinement effects are found to be much less pronounced in granular aggregates as opposed to porous structures. We show that this discrepancy can be understood in terms of the surface-surface correlation length with a connected path through the fluid domain, suggesting that this length captures the true degree of confinement. We also find that the liquid-gas phase transition in these porous materials is of second order nature near capillary critical temperature, which is shown to represent a true critical temperature, i.e., independent of the degree of disorder and the nature of the solid matrix, discrete or continuous. The critical exponents estimated here from finite-size scaling analysis suggest that this transition belongs to the 3D random field Ising model universality class as hypothesized by F. Brochard and P.G. de Gennes, with the underlying random fields induced by local disorder in fluid-solid interactions

The Effect of Confinement on Capillary Phase Transition In Granular Aggregates

Utilizing a 3D mean-field lattice-gas model, we analyze the effect of confinement on the nature of capillary phase transition in granular aggregates with varying disorder and their inverse porous structures obtained by interchanging particles and pores. Surprisingly, the confinement effects are found to be much less pronounced in granular aggregates as opposed to porous structures. We show that this discrepancy can be understood in terms of the surface-surface correlation length with a connected path through the fluid domain, suggesting that this length captures the true degree of confinement. We also find that the liquid-gas phase transition in these porous materials is of second order nature near capillary critical temperature, which is shown to represent a true critical temperature, i.e. independent of the degree of disorder and the nature of solid matrix, discrete or continuous. The critical exponents estimated here from finite-size scaling analysis suggest that this transition belongs to the 3D random field Ising model universality class as hypothesized by P.G. de Gennes, with the underlying random fields induced by local disorder in fluid-solid interactions

Microporoelastic modeling of organic-rich shales

Thesis: S.M., Massachusetts Institute of Technology, Department of Civil and Environmental Engineering, 2015.Cataloged from PDF version of thesis.Includes bibliographical references (pages 175-183).Due to their abundance, organic-rich shales are playing a critical role in re-defining the world's energy landscape leading to shifts in global geopolitics. However, technical challenges and environmental concerns continue to contribute to the slow growth of organic-rich shale exploration and exploitation worldwide. The engineering and scientific challenges arise from the extremely heterogeneous and anisotropic nature of these naturally occurring geo-composites at multiple length scales. Specifically, the anisotropic poroelastic behavior of organic-rich shales becomes of critical importance for petroleum engineers. Thus, the focus of this thesis is to capture mechanisms of first-order contribution to the effective anisotropic poroelasticity of organic-rich shales which can pave the way for more efficient and effective exploration and exploitation. We introduce an original approach for micromechanical modeling of organic-rich shales which accounts for the effect of organic maturity on the overall anisotropic poroelasticity through morphology considerations. This morphology contribution is captured by means of an effective media theory that bridges the gap between immature and mature systems through the choice of the system's micro-texture; namely a matrix-inclusion morphology (Mori-Tanaka) for immature systems and a polycrystal/ granular morphology for mature systems. Also, we show that interfaces play a role on the effective elasticity of mature organic-rich shales. The models are calibrated by means of ultrasonic pulse velocity measurements of elastic properties and validated by means of lab measured nanoindentation data. Sensitivity analyses using Spearman's Partial Rank Correlation Coefficient show the importance of porosity and Total Organic Carbon (TOC) as key input parameters for accurate model predictions. These models' developments provide a mean to define a "unique" set of clay elasticity. They also highlight the importance of the depositional environment, burial and diagenetic processes on overall mechanical and poromechanical behavior of organic-rich shales.by Siavash Khosh Sokhan Monfared.S.M

Mechanical basis and topological routes to cell elimination

Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remains largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusion with cell–cell and cell–substrate interactions in a flat monolayer. Independent tuning of cell–cell versus cell–substrate adhesion forces reveals that extrusion events can be distinctly linked to defects in nematic and hexatic orders associated with cellular arrangements. Specifically, we show that by increasing relative cell–cell adhesion forces the cell monolayer can switch between the collective tendency towards fivefold, hexatic, disclinations relative to half-integer, nematic, defects for extruding a cell. We unify our findings by accessing three-dimensional mechanical stress fields to show that an extrusion event acts as a mechanism to relieve localized stress concentration

Mechanics of Live Cell Elimination

Cell layers eliminate unwanted cells through the extrusion process, which underlines healthy versus flawed tissue behaviors. Although several biochemical pathways have been identified, the underlying mechanical basis including the forces involved in cellular extrusion remain largely unexplored. Utilizing a phase-field model of a three-dimensional cell layer, we study the interplay of cell extrusion with cell-cell and cell-substrate interactions, in a monolayer. Independent tuning of cell-cell versus cell-substrate adhesion forces in the model reveals that a higher cell-substrate adhesion leads to a lower number of total extrusion events. We find extrusion events to be linked to both half-integer topological defects in the orientation field of the cells and to five-fold disclinations in cellular arrangements. We also show that increasing the relative cell-cell adhesion forces translates into a higher likelihood for an extrusion event to be associated with a five-fold disclination and a weaker correlation with +1/2 topological defects. We unify our findings by accessing mechanical stress fields: an extrusion event acts as a mechanism to relieve localized stress concentration

Disorder-induced stiffness degradation of highly disordered porous materials

International audienceThe 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

A methodology to calibrate and to validate effective solid potentials of heterogeneous porous media from computed tomography scans and laboratory-measured nanoindentation data

International audienceBuilt 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