@article {666, title = {Modified Szeg{\"o}{\textendash}Widom Asymptotics for Block Toeplitz Matrices with Zero Modes}, journal = {Journal of Statistical Physics}, volume = {174}, year = {2019}, month = {Jan-2019}, pages = {28 - 39}, abstract = {

The Szego-Widom theorem provides an expression for the determinant of block Toeplitz matrices in the asymptotic limit of large matrix dimension n. We show that the presence of zero modes, i.e, eigenvalues that vanish as n, ||\<1, when n, requires a modification of the Szego-Widom theorem. A new asymptotic expression for the determinant of a certain class of block Toeplitz matrices with one pair of zero modes is derived. The result is inspired by one-dimensional topological superconductors, and the relation with the latter systems is discussed.

}, issn = {0022-4715}, doi = {10.1007/s10955-018-2177-8}, url = {http://link.springer.com/10.1007/s10955-018-2177-8}, author = {Basor, E. and Dubail, Jerome and Emig, Thorsten and Santachiara, Raoul} } @article {679, title = {Multiscale poromechanics of wet cement paste}, journal = {Proceedings of the National Academy of Sciences}, volume = {2}, year = {2019}, month = {May-09-2019}, pages = {201901160}, abstract = {

Capillary effects, such as imbibition drying cycles, impact the mechanics of granular systems over time. A multiscale poromechanics framework was applied to cement paste, which is the most common building material, experiencing broad humidity variations over the lifetime of infrastructure. First, the liquid density distribution at intermediate to high relative humidity is obtained using a lattice gas density functional method together with a realistic nanogranular model of cement hydrates. The calculated adsorption/desorption isotherms and pore size distributions are discussed and compare well with nitrogen and water experiments. The standard method for pore size distribution determination from desorption data is evaluated. Second, the integration of the Korteweg liquid stress field around each cement hydrate particle provided the capillary forces at the nanoscale. The cement mesoscale structure was relaxed under the action of the capillary forces. Local irreversible deformations of the cement nanograins assembly were identified due to liquid\–solid interactions. The spatial correlations of the nonaffine displacements extend to a few tens of nanometers. Third, the Love\–Weber method provided the homogenized liquid stress at the micrometer scale. The homogenization length coincided with the spatial correlation length of nonaffine displacements. Our results on the solid response to capillary stress field suggest that the micrometer-scale texture is not affected by mild drying, while nanoscale irreversible deformations still occur. These results pave the way for understanding capillary phenomena-induced stresses in heterogeneous porous media ranging from construction materials to hydrogels and living systems.

}, issn = {0027-8424}, doi = {10.1073/pnas.1901160116}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1901160116}, author = {Zhou, Tingtao and Katerina Ioannidou and Franz-Josef Ulm and Bazant, Martin Z. and Roland Jean-Marc Pellenq} } @article {533, title = {Mechanical strength of wet particle agglomerates}, journal = {Mechanics Research Communications}, volume = {92}, year = {2018}, month = {Sep-2018}, pages = {1 - 7}, abstract = {

Using particle dynamics simulations, we investigate the strength and microstructure of agglomerates of wet frictional particles subjected to axial compression. The numerical model accounts for the cohesive and viscous effects of the binding liquid up to a debonding distance with the liquid assumed to be distributed homogeneously inside the agglomerate. We show that wet agglomerates undergo plastic deformation due to the rearrangements of primary particles during compression. The compressive strength is thus characterized by the plastic threshold before the onset of failure by the irreversible loss of wet contacts between primary particles. We find that the agglomerate plastic threshold is proportional to the characteristic cohesive stress defined from the liquid-vapor surface tension and the mean diameter of primary particles, with a prefactor that is a nearly linear function of the debonding distance and increases with size span. We analyze the agglomerate microstructure and, considering only the cohesive capillary forces at all bonds between primary particles, we propose an expression of the plastic strength as a function of the texture parameters such as the wet coordination number and packing fraction. This expression is shown to be consistent with our simulations up to a multiplicative factor reflecting the distribution of the capillary bridges. (C) 2018 Published by Elsevier Ltd.

Fig. 1. Schematic drawing of the forces acting on particle i by a contacting particle j{\textellipsis}

Fig. 2. Schematic representation of diametrical compression test (a), and snapshots of{\textellipsis}

Fig. 3. Snapshot of a granule under diametrical compression for α=1

Fig. 4. Vertical strength σzz normalized by the reference stress σc for different{\textellipsis}

Fig. 5. Normalized plastic strength of the granule for several values of the debonding{\textellipsis}

Fig. 6. Normalized plastic strength σp/σc of the granules for different values of the{\textellipsis}

Fig. 7. The initial wet coordination number Z0 of the granules for different values of{\textellipsis}

Fig. 8. The prefactor η in Eq

Fig. 9. The ZΦs in Eq

}, issn = {00936413}, doi = {10.1016/j.mechrescom.2018.07.003}, url = {https://www-sciencedirect-com.libproxy.mit.edu/science/article/pii/S0093641318301216}, author = {Vo, Thanh-Trung and Patrick Mutabaruka and Saeid Nezamabadi and Jean-Yves Delenne and Izard, Edouard and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}} } @article {630, title = {Mesoscale structure, mechanics, and transport properties of source rocks{\textquoteright} organic pore networks}, journal = {Proceedings of the National Academy of Sciences}, volume = {115}, year = {2018}, month = {Dec-04-2018}, pages = {12365-12370}, abstract = {

Organic matter is responsible for the generation of hydrocarbons during the thermal maturation of source rock formation. This geochemical process engenders a network of organic hosted pores that governs the flow of hydrocarbons from the organic matter to fractures created during the stimulation of production wells. Therefore, it can be reasonably assumed that predictions of potentially recoverable confined hydrocarbons depend on the geometry of this pore network. Here, we analyze mesoscale structures of three organic porous networks at different thermal maturities. We use electron tomography with subnanometric resolution to characterize their morphology and topology. Our 3D reconstructions confirm the formation of nanopores and reveal increasingly tortuous and connected pore networks in the process of thermal maturation. We then turn the binarized reconstructions into lattice models including information from atomistic simulations to derive mechanical and confined fluid transport properties. Specifically, we highlight the influence of adsorbed fluids on the elastic response. The resulting elastic energy concentrations are localized at the vicinity of macropores at low maturity whereas these concentrations present more homogeneous distributions at higher thermal maturities, due to pores\&$\#$39; topology. The lattice models finally allow us to capture the effect of sorption on diffusion mechanisms with a sole input of network geometry. Eventually, we corroborate the dominant impact of diffusion occurring within the connected nanopores, which constitute the limiting factor of confined hydrocarbon transport in source rocks.

}, keywords = {porous media; electron tomography; mechanics; fluid transport; mesoscale}, issn = {0027-8424}, doi = {10.1073/pnas.1808402115}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1808402115}, author = {Berthonneau, Jeremie and Ama{\"e}l Obliger and Valdenaire, Pierre-Louis and Grauby, Olivier and Ferry, Daniel and Chaudanson, Damien and Pierre E. Levitz and Kim, Jae Jin and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @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 {610, title = {A minimal power model for human running performance}, journal = {PLOS ONE}, volume = {13}, year = {2018}, month = {Nov-16-2018}, pages = {e0206645}, abstract = {

Models for human running performances of various complexities and underlying principles have been proposed, often combining data from world record performances and bio-energetic facts of human physiology. The purpose of this work is to develop a novel, minimal and universal model for human running performance that employs a relative metabolic power scale. The main component is a self-consistency relation for the time dependent maximal power output. The analytic approach presented here is the first to derive the observed logarithmic scaling between world (and other) record running speeds and times from basic principles of metabolic power supply. Our hypothesis is that various female and male record performances (world, national) and also personal best performances of individual runners for distances from 800m to the marathon are excellently described by this model. Indeed, we confirm this hypothesis with mean errors of (often much) less than 1\%. The model defines endurance in a way that demonstrates symmetry between long and short racing events that are separated by a characteristic time scale comparable to the time over which a runner can sustain maximal oxygen uptake. As an application of our model, we derive personalized characteristic race speeds for different durations and distances.

}, keywords = {ENERGY-COST; OXYGEN-UPTAKE; ATHLETIC RECORDS; ECONOMY; RUNNERS; RUN; ENDURANCE; MARATHON; EXERCISE; EXHAUSTION}, doi = {10.1371/journal.pone.0206645}, url = {http://dx.plos.org/10.1371/journal.pone.0206645}, author = {Mulligan, Matthew and Adam, Guillaume and Emig, Thorsten}, editor = {Piknova, Barbora} } @article {634, title = {Multi-periodic boundary conditions and the Contact Dynamics method}, journal = {Comptes Rendus M{\'e}canique}, volume = {346}, year = {2018}, month = {Mar-2018}, pages = {263-277}, abstract = {

For investigating the mechanical behavior of granular materials by means of the discrete element approach, it is desirable to be able to simulate representative volume elements with macroscopically homogeneous deformations. This can be achieved by means of fully periodic boundary conditions such that stresses or displacements can be applied in all space directions. We present a general framework for periodic boundary conditions in granular materials and its implementation more specifically in the Contact Dynamics method. (c) 2017 Academie des sciences. Published by Elsevier Masson SAS.

}, issn = {16310721}, doi = {10.1016/j.crme.2017.12.007}, url = {https://linkinghub.elsevier.com/retrieve/pii/S1631072117302346}, author = {Farhang Radja{\"\i}} } @article {627, title = {Multiscale modeling for bioresources and bioproducts}, journal = {Innovative Food Science \& Emerging Technologies}, volume = {46}, year = {2018}, month = {Avr-2018}, pages = {41 - 53}, abstract = {

Designing and processing complex matter and materials are key objectives of bioresource and bioproduct research. Modeling approaches targeting such systems have to account for their two main sources of complexity: their intrinsic multi-scale nature; and the variability and heterogeneity inherent to all living systems. Here we provide insight into methods developed at the Food \& Bioproduct Engineering division (CEPIA) of the French National Institute of Agricultural Research (INRA). This brief survey focuses on innovative research lines that tackle complexity by mobilizing different approaches with complementary objectives. On one hand cognitive approaches aim to uncover the basic mechanisms and laws underlying the emerging collective properties and macroscopic behavior of soft-matter and granular systems, using numerical and experimental methods borrowed from physics and mechanics. The corresponding case studies are dedicated to the structuring and phase behavior of biopolymers, powders and granular materials, and to the evolution of these structures caused by external constraints. On the other hand machine learning approaches can deal with process optimizations and outcome predictions by extracting useful information and correlations from huge datasets built from experiments at different length scales and in heterogeneous conditions. These predictive methods are illustrated in the context of cheese ripening, grape maturity prediction and bacterial production.

}, issn = {14668564}, doi = {10.1016/j.ifset.2017.09.015}, url = {https://linkinghub.elsevier.com/retrieve/pii/S1466856417302230}, author = {Barnabe, M. and Blanc, Nicolas and Chabin, T. and Jean-Yves Delenne and Duri, A. and Frank, Xavier and Hugouvieux, V. and Lutton, E. and Mabille, F. and Saeid Nezamabadi and Perrot, N. and Farhang Radja{\"\i} and Ruiz, T. and Tonda, A.} } @article {410, title = {Magnetism as indirect tool for carbon content assessment in nickel nanoparticles}, journal = {Journal of Applied Physics}, volume = {122}, year = {2017}, month = {Jul-12-2017}, pages = {213902}, abstract = {

We report a combined experimental and theoretical study to ascertain carbon solubility in nickel nanoparticles embedded into a carbon matrix via the one-pot method. This original approach is based on the experimental characterization of the magnetic properties of Ni at room temperature and Monte Carlo simulations used to calculate the magnetization as a function of C content in Ni nanoparticles. Other commonly used experimental methods fail to accurately determine the chemical analysis of these types of nanoparticles. Thus, we could assess the C content within Ni nanoparticles and it decreases from 8 to around 4\ at. \% with increasing temperature during the synthesis. This behavior could be related to the catalytic transformation of dissolved C in the Ni particles into graphite layers surrounding the particles at high temperature. The proposed approach is original and easy to implement experimentally since only magnetization measurements at room temperature are needed. Moreover, it can be extended to other types of magnetic nanoparticles dissolving carbon.

}, issn = {0021-8979}, doi = {10.1063/1.5006138}, url = {http://aip.scitation.org/doi/10.1063/1.5006138http://aip.scitation.org/doi/pdf/10.1063/1.5006138}, author = {Oumellal, Y. and Yann Magnin and {\'\i}nez de Yuso, A. and Aguiar Hualde, J. M. and Amara, H. and Paul-Boncour, V. and Matei Ghimbeu, C. and Malouche, A. and Christophe Bichara and Roland Jean-Marc Pellenq and Zlotea, C.} } @article {615, title = {Many-body heat radiation and heat transfer in the presence of a nonabsorbing background medium}, journal = {Physical Review B}, volume = {95}, year = {2017}, month = {Feb-15-2017}, abstract = {

Heat radiation and near-field radiative heat transfer can be strongly manipulated by adjusting geometrical shapes, optical properties, or the relative positions of the objects involved. Typically, these objects are considered as embedded in vacuum. By applying the methods of fluctuational electrodynamics, we derive general closed-form expressions for heat radiation and heat transfer in a system of N arbitrary objects embedded in a passive nonabsorbing background medium. Taking into account the principle of reciprocity, we explicitly prove the symmetry and positivity of transfer in any such system. Regarding applications, we find that the heat radiation of a sphere as well as the heat transfer between two parallel plates is strongly enhanced by the presence of a background medium. Regarding near- and far-field transfer through a gas like air, we show that a microscopic model (based on gas particles) and a macroscopic model (using a dielectric contrast) yield identical results. We also compare the radiative transfer through a medium like air and the energy transfer found from kinetic gas theory.

}, issn = {2469-9950}, doi = {10.1103/PhysRevB.95.085413}, url = {https://link.aps.org/doi/10.1103/PhysRevB.95.085413}, author = {M{\"u}ller, Boris and Incardone, Roberta and Antezza, Mauro and Emig, Thorsten and Kr{\"u}ger, Matthias} } @article {616, title = {Mechanisms of jamming in the Nagel-Schreckenberg model for traffic flow}, journal = {Physical Review E}, volume = {95}, year = {2017}, month = {Jan-2017}, abstract = {

We study the Nagel-Schreckenberg cellular automata model for traffic flow by both simulations and analytical techniques. To better understand the nature of the jamming transition, we analyze the fraction of stopped cars P(v=0) as a function of the mean car density. We present a simple argument that yields an estimate for the free density where jamming occurs, and show satisfying agreement with simulation results. We demonstrate that the fraction of jammed cars P(v\∈{0,1}) can be decomposed into the three factors (jamming rate, jam lifetime, and jam size) for which we derive, from random walk arguments, exponents that control their scaling close to the critical density.

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.95.012311}, url = {https://link.aps.org/doi/10.1103/PhysRevE.95.012311}, author = {Bette, Henrik M. and Habel, Lars and Emig, Thorsten and Schreckenberg, Michael} } @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 {587, title = {Methodology for Estimation of Nanoscale Hardness via Atomistic Simulations}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Dec-2017}, pages = { Article Number: UNSP 04017011 }, abstract = {

Statistical mechanics has provided powerful techniques to measure mechanical properties of materials at the nanoscale and paved the way for bottom-up computational materials design. The introduction of such techniques in civil engineering applications, namely construction and geotechnical materials, remains limited to the elastic and fracture properties. This paper presents an atomistic approach to calculate the nanoscale cohesion, friction angle, and hardness. This method is based on the application of biaxial external deformation, or stress, in the weakest crystallographic direction in the material. The onset of the failure is characterized by investigating the unloading paths from several points on the stress-strain curve. Such calculations of the failure stress along different deformation paths provide multiple failure Mohr circles in the normal-shear stress space, which is found to provide a failure envelope akin to the Mohr-Coulomb failure criterion that is widely used for the plastic analysis of granular geomaterials. The failure envelope characterizes the nanoscale cohesion and friction angle, which in conjunction with continuum mechanics can be utilized to estimate the nanoscale hardness of layered materials. Application of this method to tobermorite and Na-montmorillonite crystals yields values that are close to the experimental measurements obtained using nanoindentation and atomic force microscopy techniques. (C) 2017 American Society of Civil Engineers.

}, keywords = {Atomistic simulation; Friction; Cohesion; Hardness; Tobermorite; Montmorillonite}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000127}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000127}, author = {Qomi, M. J. Abdolhosseini and Ebrahimi, Davoud and Mathieu Bauchy and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {241, title = {Model representations of kerogen structures: An insight from density functional theory calculations and spectroscopic measurements}, journal = {Scientific Reports}, volume = {7}, year = {2017}, month = {Aug-01-2017}, pages = {Article Number: 7068}, abstract = {

Molecular structures of kerogen control hydrocarbon production in unconventional reservoirs. Significant progress has been made in developing model representations of various kerogen structures. These models have been widely used for the prediction of gas adsorption and migration in shale matrix. However, using density functional perturbation theory (DFPT) calculations and vibrational spectroscopic measurements, we here show that a large gap may still remain between the existing model representations and actual kerogen structures, therefore calling for new model development. Using DFPT, we calculated Fourier transform infrared (FTIR) spectra for six most widely used kerogen structure models. The computed spectra were then systematically compared to the FTIR absorption spectra collected for kerogen samples isolated from Mancos, Woodford and Marcellus formations representing a wide range of kerogen origin and maturation conditions. Limited agreement between the model predictions and the measurements highlights that the existing kerogen models may still miss some key features in structural representation. A combination of DFPT calculations with spectroscopic measurements may provide a useful diagnostic tool for assessing the adequacy of a proposed structural model as well as for future model development. This approach may eventually help develop comprehensive infrared (IR)-fingerprints for tracing kerogen evolution.

}, doi = {10.1038/s41598-017-07310-9}, author = {Weck, Philippe F. and Kim, Eunja and Wang, Yifeng and Kruichak, Jessica N. and Mills, Melissa M. and Matteo, Edward N. and Roland Jean-Marc Pellenq} } @article {264, title = {Modeling Granular Materials: Century-Long Research across Scales}, journal = {Journal of Engineering Mechanics}, volume = {143}, year = {2017}, month = {Jan-23-2017}, pages = {Article Number: UNSP 04017002}, abstract = {

Granular materials are the most recurrent form of solid-state matter on Earth. They challenge researchers and engineers in various fields not only because they occur with a broad variety of grain sizes, shapes and interactions in nature and industry, but also because they show a rich panoply of mechanical states. Despite this polymorphism, all these different types of soils, powders, granules, ores, pharmaceutical products, etc., are instances of the granular matter with the same least common denominator of being sandlike (psammoid in Greek), i.e.,\ solid grains interacting via frictional contacts. This review describes milestone contributions to the field of granular materials since the early elastic-plastic models developed for soils in the 1950s. The research on granular materials has grown into a vast multidisciplinary field in the 1980s with increasing focus on the microstructure and owing to new experimental tools and discrete simulation methods. It turns out that the granular texture, particle-scale kinematics, and force transmission are far more complex than presumed in early micromechanical models of granular materials. Hence, constitutive relations cannot easily be derived from the particle-scale behavior although advanced continuum models have been developed to account for anisotropy, intermediate stress, and complex loading paths. The subtle elastic properties and origins of bulk friction will be discussed, as well as the effects of particle shape and size distributions. The review covers also recent developments in macroscopic modeling such as the thermomechanical approach, anisotropic critical state theory, nonlocal modeling approach, inertial flows, and material instabilities. Finally, a brief account is given of open issues and some new frontiers and challenges in the field.

}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0001196}, author = {Farhang Radja{\"\i} and Jean-Noel~ Roux and Daouadji, Ali} } @article {606, title = {Modeling root growth in granular soils: effects of root stiffness and packing fraction}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 14013}, abstract = {

We use molecular dynamics simulations to investigate the effects of root bending stiffness and packing fraction on the path followed by a growing root in 2D packings of grains representing a soil. The root is modeled as a chain of elements that can grow in length and change their direction depending on the forces exerted by soil grains. We show that the root shape is mainly controlled by the bending stiffness of its apex. At low stiffness, the root randomly explores the pore space whereas at sufficiently high stiffness, of the order of soil hardness multiplied by mean grain size, the root follows a straight path across the soil. Between these two limits, the root shape can be characterized by the standard deviation of its re-directions at the scale of soil grains. We find that this shape parameter varies as a power-law function of the normalized bending stiffness.

}, doi = {10.1051/epjconf/201714014013}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714014013}, author = {Fakih, Mahmoud and Jean-Yves Delenne and Farhang Radja{\"\i} and Fourcaud, Thierry}, editor = {Saeid Nezamabadi and Luding, S.} } @article {265, title = {Modeling soft granular materials}, journal = {Granular Matter}, volume = {19}, year = {2017}, month = {Jun-17-2017}, pages = {Article Number: 8}, abstract = {

Soft-grain materials such as clays and other colloidal pastes share the common feature of being composed of grains that can undergo large deformations without rupture. For the simulation of such materials, we present two alternative methods: (1) an implicit formulation of the material point method (MPM), in which each grain is discretized as a collection of material points, and (2) the bonded particle model (BPM), in which each soft grain is modeled as an aggregate of rigid particles using the contact dynamics method. In the MPM, a linear elastic behavior is used for the grains. In order to allow the aggregates in the BPM to deform without breaking, we use long-range center-to-center attraction forces between the primary particles belonging to each grain together with steric repulsion at their contact points. We show that these interactions lead to a plastic behavior of the grains. Using both methods, we analyze the uniaxial compaction of 2D soft granular packings. This process is nonlinear and involves both grain rearrangements and large deformations. High packing fractions beyond the jamming state are reached as a result of grain shape change for both methods. We discuss the stress-strain and volume change behavior as well as the evolution of the connectivity of the grains. Similar textures are observed at large deformations although the BPM requires higher stress than the MPM to reach the same level of packing fraction.

}, issn = {1434-5021}, doi = {10.1007/s10035-016-0689-y}, author = {Saeid Nezamabadi and Thanh Hai Nguyen and Jean-Yves Delenne and Farhang Radja{\"\i}} } @proceedings {336, title = {Modelling Transient Dynamics of Granular Slopes: MPM and DEM}, journal = {PROCEEDINGS OF THE 1ST INTERNATIONAL CONFERENCE ON THE MATERIAL POINT METHOD (MPM 2017) Book Series: Procedia Engineering }, volume = {175}, year = {2017}, month = {Jul-17-2017}, pages = {94 - 101}, publisher = {Elsevier Ltd.}, address = {JAN 10-13, 2017, Delft, NETHERLANDS}, abstract = {

Transient granular flows, such as rock falls, debris flows, and aerial and submarine avalanches, occur very often in nature. In the geotechnical context, transient movements of large granular slopes are a substantial factor of risk due to their destructive force and the transformations the y may produce in the landscape. This paper investigates the ability of MPM, a continuum approach, to reproduce the evolution of a granular slope destabilised by an external energy source. In particular, a central issue is whether the power-law dependence of run-out distance and time observed with respect to the initial geometry or energy can be reproduced by a simple Mohr-Coulomb plastic behaviour. The effect of base friction on the run-out kinematics is studied by comparing the data obtained from the DEM and MPM simulations. The mechanism of energy dissipation is primarily through friction and the MPM is able to predict the run-out response in good agreement with the DEM simulations. At very low excitation energies, the DEM simulations show longer run-out in comparison to the MPM due to local destabilization at the flow front. At low input energies, a larger fraction of the energy is consumed in the destabilisation process, hence the amount energy available for flow is less. However, at higher input energy, where most of the energy is dissipated during the spreading phase, the run-out distance has a weak dependence on the distribution of velocity in the granular mass.

}, issn = {18777058}, doi = {10.1016/j.proeng.2017.01.032}, author = {Krishna Kumar and Kenichi Soga and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Rohe, A and Kenichi Soga and Teunissen, H} } @article {246, title = {Molecular Modeling and Adsorption Properties of Ordered Silica-Templated CMK Mesoporous Carbons}, journal = {Langmuir}, volume = {33}, year = {2017}, month = {Feb-06-2017}, pages = {2109 - 2121}, abstract = {

Realistic molecular models of silica-templated CMK-1, CMK-3, and CMK-5 carbon materials have been developed by using carbon rods and carbon pipes that were obtained by adsorbing carbon in a model MCM-41 pore. The interactions between the carbon atoms with the silica matrix were described using the PN-Traz potential, and the interaction between the carbon atoms was calculated by the reactive empirical bond order (REBO) potential. Carbon rods and pipes with different thicknesses were obtained by changing the silica\–carbon interaction strength, the temperature, and the chemical potential of carbon vapor adsorption. These equilibrium structures were further used to obtain the atomic models of CMK-1, CMK-3, and CMK-5 materials using the same symmetry as found in TEM pictures. These models are further refined and made more realistic by adding interconnections between the carbon rods and carbon pipes. We calculated the geometric pore size distribution of the different models of CMK-5 and found that the presence of interconnections results in some new features in the pore size distribution. Argon adsorption properties were investigated using GCMC simulations to characterize these materials at 77 K. We found that the presence of interconnection results greatly improves the agreement with available experimental data by shifting the capillary condensation to lower pressures. Adding interconnections also induces smoother adsorption/condensation isotherms, and desorption/evaporation curves show a sharp jump. These features reflex the complexity of the nanovoids in CMKs in terms of their pore morphology and topology.

Abstract Image

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.6b04169}, author = {Jain, Surendra Kumar and Roland Jean-Marc Pellenq and Keith E. Gubbins and Peng, Xuan} } @proceedings {337, title = {MPM with Frictional Contact for Application to Soft Particulate Materials}, journal = {PROCEEDINGS OF THE 1ST INTERNATIONAL CONFERENCE ON THE MATERIAL POINT METHOD (MPM 2017) Book Series: Procedia Engineering }, volume = {175}, year = {2017}, month = {Mar-23-2017}, pages = {141 - 147}, publisher = {Elsevier Ltd.}, address = {JAN 10-13, 2017, Delft, NETHERLANDS}, abstract = {

Soft particle materials are composed of discrete particles that can undergo large deformations without rupture. Most food products, many powders, colloidal pastes, vesicles and biological cells are soft particle systems. In order to model such materials, we present an efficient numerical approach combining an implicit formulation of the Material Point Method (MPM) and Contact Dynamics (CD) method. The MPM deals with bulk variables of an individual particle by discretizing it as a collection of material points, whereas the CD allows for the treatment of frictional contacts between particles. This model is applied for the simulation of the uniaxial compression of 2D soft-particle packings. The compaction is a nonlinear process in which new contacts are formed between particles and the contact areas increase. The change of particle shapes allows these materials to reach high packing fraction. We find that the contact specific surface, the orientation anisotropy and the aspect ratio of particles increase as a function of the packing fraction but at different rates. We also evidence the effect of friction, which favors strong stress chains and thus the elongation of particles, leading to larger values of the orientation anisotropy and the aspect ratio at a given level of packing fraction as compared to a frictionless particle packing.

}, issn = {18777058}, doi = {10.1016/j.proeng.2017.01.044}, author = {Saeid Nezamabadi and Thanh Hai Nguyen and Jean-Yves Delenne and Julien Averseng and Frank, Xavier and Farhang Radja{\"\i}}, editor = {Rohe, A and Kenichi Soga and Teunissen, H} } @article {196, title = {Mesoscale texture of cement hydrates}, journal = {Proceedings of the National Academy of Sciences}, volume = {113}, year = {2016}, month = {Feb-23-2016}, pages = {2029 - 2034}, abstract = {

Strength and other mechanical properties of cement and concrete rely upon the formation of calcium\–silicate\–hydrates (C\–S\–H) during cement hydration. Controlling structure and properties of the C\–S\–H phase is a challenge, due to the complexity of this hydration product and of the mechanisms that drive its precipitation from the ionic solution upon dissolution of cement grains in water. Departing from traditional models mostly focused on length scales above the micrometer, recent research addressed the molecular structure of C\–S\–H. However, small-angle neutron scattering, electron-microscopy imaging, and nanoindentation experiments suggest that its mesoscale organization, extending over hundreds of nanometers, may be more important. Here we unveil the C\–S\–H mesoscale texture, a crucial step to connect the fundamental scales to the macroscale of engineering properties. We use simulations that combine information of the nanoscale building units of C\–S\–H and their effective interactions, obtained from atomistic simulations and experiments, into a statistical physics framework for aggregating nanoparticles. We compute small-angle scattering intensities, pore size distributions, specific surface area, local densities, indentation modulus, and hardness of the material, providing quantitative understanding of different experimental investigations. Our results provide insight into how the heterogeneities developed during the early stages of hydration persist in the structure of C\–S\–H and impact the mechanical performance of the hardened cement paste. Unraveling such links in cement hydrates can be groundbreaking and controlling them can be the key to smarter mix designs of cementitious materials.

Fig. 1.

}, issn = {0027-8424}, doi = {10.1073/pnas.1520487113}, author = {Katerina Ioannidou and Konrad J. Krakowiak and Mathieu Bauchy and Christian G. Hoover and Enrico Masoero and Sidney Yip and Franz-Josef Ulm and Pierre E. Levitz and Roland Jean-Marc Pellenq and Emanuela Del Gado} } @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} } @article {645, title = {Magnetic nanopantograph in the SrCu 2 (BO 3 ) 2 Shastry{\textendash}Sutherland lattice}, journal = {Proceedings of the National Academy of Sciences}, volume = {112}, year = {2015}, month = {Feb-17-2015}, pages = { 1971-1976}, abstract = {

Magnetic materials having competing, i.e., frustrated, interactions can display magnetism prolific in intricate structures, discrete jumps, plateaus, and exotic spin states with increasing applied magnetic fields. When the associated elastic energy cost is not too expensive, this high potential can be enhanced by the existence of an omnipresent magnetoelastic coupling. Here we report experimental and theoretical evidence of a nonnegligible magnetoelastic coupling in one of these fascinating materials, SrCu2(BO3)(2) (SCBO). First, using pulsed-field transversal and longitudinal magnetostriction measurements we show that its physical dimensions, indeed, mimic closely its unusually rich field-induced magnetism. Second, using density functional-based calculations we find that the driving force behind the magnetoelastic coupling is the CuOCu superexchange angle that, due to the orthogonal Cu2+ dimers acting as pantographs, can shrink significantly (0.44\%) with minute (0.01\%) variations in the lattice parameters. With this original approach we also find a reduction of similar to 10\% in the intradimer exchange integral J, enough to make predictions for the highly magnetized states and the effects of applied pressure on SCBO.

}, issn = {0027-8424}, doi = {10.1073/pnas.1421414112}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1421414112}, author = {Guillaume Radtke and Andres Sa{\`u}l and Dabkowska, Hanna A. and Salamon, Myron B. and Jaime, Marcelo} } @article {117, title = {Magnetic nanopantograph in the SrCu2(BO3)2 Shastry{\textendash}Sutherland lattice}, journal = {Proceedings of the National Academy of Sciences}, volume = {112}, year = {2015}, month = {Feb-02-2015}, pages = {1971-1976}, abstract = {

Magnetic materials having competing, i.e., frustrated, interactions can display magnetism prolific in intricate structures, discrete jumps, plateaus, and exotic spin states with increasing applied magnetic fields. When the associated elastic energy cost is not too expensive, this high potential can be enhanced by the existence of an omnipresent magnetoelastic coupling. Here we report experimental and theoretical evidence of a nonnegligible magnetoelastic coupling in one of these fascinating materials, SrCu2(BO3)2 (SCBO). First, using pulsed-field transversal and longitudinal magnetostriction measurements we show that its physical dimensions, indeed, mimic closely its unusually rich field-induced magnetism. Second, using density functional-based calculations we find that the driving force behind the magnetoelastic coupling is the CuOCu^ superexchange angle that, due to the orthogonal Cu2+ dimers acting as pantographs, can shrink significantly (0.44\%) with minute (0.01\%) variations in the lattice parameters. With this original approach we also find a reduction of \∼10\% in the intradimer exchange integral J, enough to make predictions for the highly magnetized states and the effects of applied pressure on SCBO.

}, doi = {10.1073/pnas.1421414112}, author = {Guillaume Radtke and Andres Sa{\`u}l and Dabkowska, Hanna A. and Salamon, Myron B. and Jaime, Marcelo} } @article {185, title = {The Maya blue nanostructured material concept applied to colouring geopolymers}, journal = {RSC Adv.}, volume = {5}, year = {2015}, month = {Nov-2015}, pages = {98834 - 98841}, abstract = {

Maya blue is an ancient nanostructured pigment synthetized by assembling indigo, a natural dye, with palygorskite, a microfibrous clay mineral. The novelty of our approach is to mimic \“pre-Columbian nanotechnology\” and to functionalize geopolymers with a sepiolite-based hybrid organic\–inorganic nanocomposite inspired from the Maya blue. It is acid- and UV-resistant, as confirmed by the stability of Maya mural paintings over time. We synthesized analogous pigments, using methylene blue (MB) and methyl red (MR) as organic dyes and sepiolite as fibrous clay mineral. We used an aqueous and a solid-state method, both leading to encapsulation of dye monomers into the clay micropores, as confirmed by UV-vis spectroscopy. This nanostructured pigment was then included into a geopolymer matrix at room temperature. The stability of the new material to UV and acid was tested. It was confirmed that it is the prior encapsulation of the dye into sepiolite that leads to the stability of the pigment in the geopolymer matrix. This first study opens the way to numerous possibilities for functionalizing inorganic binder materials with organic elements that would be otherwise sensitive to thermal treatment in conventional ceramic processing.

}, doi = {10.1039/C5RA14076E}, author = {Ouellet-Plamondon, Claudiane and Aranda, Pilar and Favier, Aur{\'e}lie and Habert, Guillaume and Henri Van Damme and Ruiz-Hitzky, Eduardo} } @article {183, title = {Mechanism of H2O Insertion and Chemical Bond Formation in AlPO 4-54{\textperiodcentered}x H2O at High Pressure}, journal = {Journal of the American Chemical Society}, volume = {137}, year = {2015}, month = {Jan-21-2015}, pages = {584 - 587}, abstract = {

The insertion of H2O in AlPO4-54\•xH2O at high pressure was investigated by single-crystal x-ray diffraction and Monte Carlo molecular simulation. H2O molecules are concentrated, in particular near the pore walls. Upon insertion, the additional water is highly disordered. Insertion of H2O (superhydration) is found to impede pore collapse in the material, thereby strongly modifying its mechanical behavior. However, instead of stabilizing the structure with respect to amorphization, the results provide evidence for the early stages of chemical bond formation between H2O molecules and tetrahedrally coordinated aluminum, which is at the origin of the amorphization/reaction process.

}, issn = {0002-7863}, doi = {10.1021/ja511153n}, author = {Alabarse, Frederico G. and Rouquette, J{\'e}r{\^o}me and Benoit A. Coasne and Haidoux, Abel and Paulmann, Carsten and Cambon, Olivier and Haines, Julien} } @article {145, title = {Mechanism of strength reduction along the graphenization pathway}, journal = {Science Advances}, volume = {1}, year = {2015}, month = {Nov-20-2015}, pages = {e1501009 - e1501009}, abstract = {

Even though polycrystalline graphene has shown a surprisingly high tensile strength, the influence of inherent grain boundaries on such property remains unclear. We study the fracture properties of a series of polycrystalline graphene models of increasing thermodynamic stability, as obtained from a long molecular dynamics simulation at an elevated temperature. All of the models show the typical and well-documented brittle fracture behavior of polycrystalline graphene; however, a clear decrease in all fracture properties is observed with increasing annealing time. The remarkably high fracture properties obtained for the most disordered (less annealed) structures arise from the formation of many nonpropagating prefracture cracks, significantly retarding failure. The stability of these reversible cracks is due to the nonlocal character of load transfer after a bond rupture in very disordered systems. It results in an insufficient strain level on neighboring bonds to promote fracture propagation. Although polycrystallinity seems to be an unavoidable feature of chemically synthesized graphenes, these results suggest that targeting highly disordered states might be a convenient way to obtain improved mechanical properties.

}, doi = {10.1126/sciadv.1501009}, author = {Antonio Gamboa and Baptiste Farbos and Aurel, P. and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @proceedings {330, title = {The Meso-Scale Texture of Cement Hydrate Gels: Out-of-Equilibrium Evolution and Thermodynamic Driving}, journal = {10th International Conference on Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete and Concrete Structures CONCREEP 10}, volume = {CONCREEP 10: MECHANICS AND PHYSICS OF CREEP, SHRINKAGE, AND DURABILITY OF CONCRETE AND CONCRETE STRUCTURES}, year = {2015}, month = {Sep-17-2015}, pages = {34-38}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015 Vienna, AustriaReston, VA}, abstract = {

By the end of cement hydration calcium-silicate-hydrate (C-S-H) gels extends over tens and hundreds of nanometers. Their complex texture affects directly, and to a large extent, the macroscopic hygrothermal and mechanical behavior of cement. Here we review a statistical physics approach recently developed, which allows us to investigate the gel formation under the out-of-equilibrium conditions typical of cement hydration and the role of the nano-scale structure in C\—S\—H mechanics upon hardening. Our investigations have unveiled the role, in the C-S-H gels, of nano-scale structural and mechanical heterogeneities that develop due to the the far-from-equilibrium physico-chemical environment in which the material forms. A subtle interplay between the out-of-equilibrium evolution and the effective interactions emerging between the nano-scale units of the gels at different stages of the hydration process ultimately determines the mesoscale texture of cement hydrates and their material properties.

}, doi = {10.1061/978078447934610.1061/9780784479346.005}, author = {Emanuela Del Gado and Katerina Ioannidou and Enrico Masoero and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Sidney Yip}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @proceedings {358, title = {Micromechanical analysis of the surface erosion of a cohesive soil by means of a coupled LBM-DEM model}, journal = {IV International Conference on Particle-based Methods (PARTICLES 2015)}, volume = {PARTICLE-BASED METHODS IV-FUNDAMENTALS AND APPLICATIONS}, year = {2015}, month = {sept-30-2015}, pages = {519-528}, publisher = {International Center for Numerical Methods in Engineering (CIMNE)}, address = {SEP 28-30 2015 Barcelona, SPAIN}, abstract = {

The elementary mechanisms driving the ubiquitous surface erosion of cohesive geomaterials can be analysed from a micromechanical perspective combining well-known numerical techniques. Here, a coupled model combining the Discrete Element and Lattice Boltzmann methods (DEM-LBM) provides an insight into the solid-fluid interaction during the transient erosion caused by a vertical fluid jet impinging on the surface of a granular assembly. The brittle cementation providing cohesion between the solid grains is described here by means of a simple bond model with a single-parameter yield surface. The initial topology of the surface erosion tends to mimic the profile of fluid velocity directly above the soil surface. We find that both the rate of erosion and the magnitude of eroded mass depend directly on the micromechanical strength of the single solid bonds.

}, keywords = {Cohesion, DEM, LBM, Surface erosion}, url = {https://hal.archives-ouvertes.fr/hal-01269324}, author = {Cuellar, Pablo and Pierre Philippe and St{\'e}phane Bonelli and Benahmed, Nadia and Brunier-Coulin, Florian and Jeff Ngoma and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Onate, E and Bischoff, M and Owen, DRJ} } @article {155, title = {Modeling force transmission in granular materials}, journal = {Comptes Rendus Physique}, volume = {16}, year = {2015}, month = {Feb-2015}, pages = {3 - 9}, abstract = {

The probability density function of contact forces in granular materials has been extensively studied and modeled as an outstanding signature of granular microstructure. Arguing that particle environments play a fundamental role in force transmission, we analyze the effects of steric constraints with respect to force balance condition and show that each force may be considered as resulting from a balance between lower and larger forces in proportions that mainly depend on steric effects. This idea leads to a general model that predicts an analytical expression of force density with a single free parameter. This expression fits well our simulation data and generically predicts the exponential fall-off of strong forces, a small peak below the mean force and the non-zero probability of vanishingly small forces.

}, issn = {16310705}, doi = {10.1016/j.crhy.2015.01.003}, author = {Farhang Radja{\"\i}} } @article {109, title = {Magnetic interactions in 3d metal chains on Cu2X/Cu(001) (X = N, O): Comparison with corresponding unsupported chains}, journal = {Physical Review B}, volume = {90}, year = {2014}, month = {Nov-2014}, pages = {195423}, type = {Article}, abstract = {

In this work we present a systematic study of the magnetic interactions within 3d transition-metal chains adsorbed on Cu2N and Cu2O monolayers grown on Cu(001). We are interested in the particular geometric adsorption configuration which gives rise, after relaxation, to the development of diatomic TM-X (X = N, O) chains. By using density functional theory (DFT), we calculate the energy difference between the ferromagnetic and antiferromagnetic intrachain configurations for Ti, V, Cr, Mn, Fe, and Co. Both substrates give rise, with minor differences, to the same magnetic trends, the only chains which are ferromagnetic after adsorption are Cr chains. By performing similar calculations in unsupported chains and introducing a tight-binding-model Hamiltonian based on physically reasonable assumptions we reproduce the magnetic trends obtained from the DFT calculations.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.90.195423}, author = {Urdaniz, M. C. and Barral, M. A. and Llois, A. M. and Andres Sa{\`u}l} } @article {113, title = {Mesoscale properties of clay aggregates from potential of mean force representation of interactions between nanoplatelets}, journal = {Journal of Chemical Physics }, volume = {140}, year = {2014}, month = {Apr-21-2014}, pages = {Article Number: 154309}, type = {Article}, abstract = {

Face-to-face and edge-to-edge free energy interactions of Wyoming Na-montmorillonite platelets were studied by calculating potential of mean force along their center to center reaction coordinate using explicit solvent (i.e., water) molecular dynamics and free energy perturbation methods. Using a series of configurations, the Gay-Berne potential was parametrized and used to examine the mesoscale aggregation and properties of platelets that are initially random oriented under isothermal-isobaric conditions. Aggregates of clay were defined by geometrical analysis of face-to-face proximity of platelets with size distribution described by a log-normal function. The isotropy of the microstructure was assessed by computing a scalar order parameter. The number of platelets per aggregate and anisotropy of the microstructure both increases with platelet plan area. The system becomes more ordered and aggregate size increases with increasing pressure until maximum ordered state at confining pressure of 50 atm. Further increase of pressure slides platelets relative to each other leading to smaller aggregate size. The results show aggregate size of (3-8) platelets for sodium-smectite in agreement with experiments (3-10). The geometrical arrangement of aggregates affects mechanical properties of the system. The elastic properties of the meso-scale aggregate assembly are reported and compared with nanoindentation experiments. It is found that the elastic properties at this scale are close to the cubic systems. The elastic stiffness and anisotropy of the assembly increases with the size of the platelets and the level of external pressure. (C) 2014 AIP Publishing LLC.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
}, issn = {0021-9606}, doi = {10.1063/1.4870932}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @proceedings {251, title = {Mesoscale simulation of clay aggregate formation and mechanical properties}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro: Granular Matter}, year = {2014}, month = {Sep-2014}, pages = {539-544}, address = {SEP 01-03 2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

This paper proposes a novel methodology for understanding the meso-scale aggregation of clay platelets in water. We use Molecular Dynamics simulations using the CLAYFF force fields to represent the interactions between two layers of Wyoming montmorillonite (Na-smectite) in bulk water. The analyses are used to establish the potential of mean force at different spacings between the layers for edge-to-edge and face-to-face interactions. This is accomplished by finding the change in free energy as a function of the separation distance between the platelets using thermodynamic perturbation theory with a simple overlap sampling method. These nanoscale results are then used to calibrate the Gay\–Berne (GB) potential that represents each platelet as a single-site ellipsoidal body. A coarse-graining upscaling approach then uses the GB potentials and molecular dynamics to represent the meso-scale aggregation of clay platelets (at submicron length scale). Results from meso-scale simulations obtain the equilibrium/jamming configurations for mono-disperse clay platelets. The results show aggregation for a range of clay platelets dimensions and pressures with mean stack size ranging from 3 to 8 platelets. The particle assemblies become more ordered and exhibit more pronounced elastic anisotropy at higher confining pressures. The results are in good agreement with previously measured nano-indentation moduli over a wide range of clay packing densities.

}, issn = {1434-5021}, doi = {10.1007/s10035-016-0655-8}, author = {Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Andrew J. Whittle}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin} } @proceedings {99, title = {Modelling cement at fundamental scales: From atoms to engineering strength and durability}, journal = {Euro-C Conference}, volume = {COMPUTATIONAL MODELLING OF CONCRETE STRUCTURES, VOL 1}, year = {2014}, month = {2014}, pages = {139-148}, publisher = {CRC Press-Taylor \& Francis Group}, address = {MAR 24-27, 2014, St Anton am Alberg, AUSTRIA}, abstract = {

The mechanics and rheology of cement pastes are largely determined by the Calcium-Silicate-Hydrate (C-S-H) gel, which is normally the main binding phase. The C-S-H gel is an amorphous material with a nanoscale pore network. Here we consider a model structure of the gel, where a continuous molecular model of non-porous C-S-H is combined with a colloidal model that incorporates the porosity of the gel. Numerical simulations predict mechanical properties and microscopic processes associated with yielding under shear. The proposed model can now be used to extract further properties and to contribute to a broader effort aimed at describing and engineering the multi-scale structure and mechanical behaviour of the cement paste. We finally discuss the opportunities offered by our model with respect to the scientific and technological challenges of modelling and controlling ageing and creep over long time scales.

}, isbn = {978-1-138-02641-4; 978-1-315-76203-6}, doi = {10.1201/b16645-16}, author = {Enrico Masoero and Hamlin M. Jennings and Franz-Josef Ulm and Emanuela Del Gado and Hegoi Manzano and Roland Jean-Marc Pellenq and Sidney Yip}, editor = {Bicanic, N and Mang, H and Meschke, G and DeBorst, R} } @proceedings {342, title = {Modelling soft-particle materials}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro}, year = {2014}, month = {Aug-26-2014}, pages = {43-48}, publisher = {CRC Press}, address = {SEP 01-03-2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

Soft-particle materials include colloidal pastes, vesicles, many powders, microgels and suspensions. They share the common feature of being composed of particles that can undergo large deformations without rupture. For the simulation of such materials, we present a modelling approach based on an implicit formulation of the Material Point Method (MPM) interfaced with the Contact Dynamics (CD) method for the treatment of frictional contacts between particles. Each particle is discretized as a collection of material points. The information carried by the material points is projected onto a background mesh, where equations of motion are solved. The mesh solution is then used to update the material points. The implicit formulation of MPM allows for unconditional numerical stability and efficient coupling with implicit treatment of unilateral contacts and friction between the particles by the CD method. We use this model to analyse the compaction process of 2D soft-particle packings. The packing can reach high solid fractions by particle shape change and still flow plastically. The compaction is a nonlinear process in which new contacts are formed between particles and the contact areas increase. We find that the evolution of the packing fraction is a slow logarithmic function of the driving stress as a consequence of increasing contact area. We also evidence the effect of friction, which favours strong stress chains and thus the elongation of particles, leading to a larger packing fraction at a given level of compressive stress as compared to a frictionless particle packing.

}, url = {http://prodinra.inra.fr/record/370208}, author = {Saeid Nezamabadi and Farhang Radja{\"\i} and Julien Averseng}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin} } @article {650, title = {Molecular dynamics simulation of amorphous HfO 2 for resistive RAM applications}, journal = {Modelling and Simulation in Materials Science and Engineering}, volume = {22}, year = {2014}, month = {Sep-2014}, pages = { Article Number 065006 }, abstract = {

HfO2 is widely investigated as the favoured material for resistive RAM device implementation. The structural features of HfO2 play a fundamental role in the switching mechanisms governing resistive RAM operations, and a comprehensive understanding of the relation between the atomistic properties and final device behaviour is still missing. In addition, despite the fact that ultra-scaled 10 nm resistive RAM will probably be made of amorphous HfO2, a deeper investigation of the structure is necessary. In this paper, the classical molecular dynamics technique was used to investigate the disordered atomic configuration of amorphous HfO2. The influence of density on both the atomistic structure and the diffusion of O species was carefully analysed. The results achieved show that the atomistic structure of an amorphous HfO2 system is strongly affected by the density, and the amorphous system is rearranged in an atomic configuration similar to the crystalline configuration at similar densities. The diffusion of oxygen atoms increases with the decrease of the density, consistent with a less-packed atomic structure which allows for easier movement of this species.

}, issn = {0965-0393}, doi = {10.1088/0965-0393/22/6/065006}, url = {http://stacks.iop.org/0965-0393/22/i=6/a=065006?key=crossref.b7d708e3a43b9124ea539668eee311b0}, author = {Broglia, G and Guido Ori and Larcher, L and Montorsi, M} } @article {69, title = {Molecular dynamics simulations of the formation of 1D spin-valves from stretched Au-Co and Pt-Co nanowires}, journal = {J. Phys.: Condens. Matter}, volume = {26}, year = {2014}, month = {Oct-2014}, chapter = {474206}, abstract = {

We have performed molecular dynamics (MD) simulations of stretched Aux-Co1 - x and Ptx-Co1 - x nanowires to investigate the formation of bimetallic monoatomic wires between two electrodes. We have considered nanowires with two concentrations x = 0.2 and 0.8, aspect ratio of 13, a cross section of 1 nm(2) and a wide range of temperatures (from 10 to 400 K). For the MD simulations we have used a semi-empirical interatomic potential based on the second moment approximation (SMA) of the density of states to the tight-binding Hamiltonian.For Au-Co alloys, Au atoms tends to migrate towards the narrowed region to form almost pure Au wires. In the PtCo case the formed chains usually consist of Pt enriched alternating structures. The most striking result is probably the Au(0.2)-Co(0.8) alloy where pure monoatomic Au chains form between two Co electrodes constituting a potential 1D spin valve. Despite the known ease with which the 5d metals (Pt, Ir, and Au) form monoatomic chains (MACS), our results show that in the presence of Co (x = 0.2), the percentage of chain formation is higher than in the Pt and Au rich cases (x = 0.8).

}, doi = {10.1088/0953-8984/26/47/474206}, author = {Robinson Cortes-Huerto and T. Sondon and Andres Sa{\`u}l} } @article {70, title = {Multiphase equation of state for carbon addressing high pressures and temperatures}, journal = {Phys Rev. B}, volume = {89}, year = {2014}, month = {Jun-2014}, chapter = {224109}, abstract = {

We present a 5-phase equation of state for elemental carbon which addresses a wide range of density and temperature conditions: 3g/cc\<\ρ\<20g/cc,0K\<T\<\∞. The phases considered are diamond, BC8, simple cubic, simple hexagonal, and the liquid/plasma state. The solid phase free energies are constrained by density functional theory (DFT) calculations. Vibrational contributions to the free energy of each solid phase are treated within the quasiharmonic framework. The liquid free energy model is constrained by fitting to a combination of DFT molecular dynamics performed over the range 10000K\<T\<100000K, and path integral quantum Monte Carlo calculations for T\>100000K (both for \ρ between 3 and 12 g/cc, with select higher- \ρ DFT calculations as well). The liquid free energy model includes an atom-in-jellium approach to account for the effects of ionization due to temperature and pressure in the plasma state, and an ion-thermal model which includes the approach to the ideal gas limit. The precise manner in which the ideal gas limit is reached is greatly constrained by both the highest-temperature DFT data and the path integral data, forcing us to discard an ion-thermal model we had used previously in favor of a new one. Predictions are made for the principal Hugoniot and the room-temperature isotherm, and comparisons are made to recent experimental results.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevB.89.224109}, author = {L. X. Benedict and K. P. Driver and S. Hamel and B. Militzer and T. Qi and Alfredo A. Correa and Andres Sa{\`u}l and E. Schwegler} } @article {23, title = {Mesoscopic Monte Carlo simulations of microstructure and conductivity of zirconia-based ceramics}, journal = {J. Eur. Ceram. Soc.}, volume = {61}, year = {2013}, month = {Aug-2013}, pages = {5018-5025}, abstract = {

Zirconia-based refractories are ceramics that serve as furnaces for the production of high-quality speciality glasses. These materials are obtained by cooling a melt composed of different oxides and oxide precursors. The microstructure of these refractories consists of monoclinic zirconia embedded in a three-dimensional interconnected amorphous phase. Using mesoscopic Monte Carlo simulations, the present paper addresses the effect of the amount of amorphous phase and temperature on the microstructure of the material. The simulated microstructures resemble those obtained using tomography experiments. A theoretical model describing the electrical behavior of the material as a function of temperature and composition is also reported. In this model, the conductivity of the ZrO2 particles is assumed to be constant, but the overall conductivity of the sample depends on its tortuosity s. Comparison with experimental data suggests that the present model provides a realistic picture of the electrical behavior of zirconia-based ceramics and leads to quantitative predictions.

Image for unlabelled figure

}, doi = {10.1016/j.actamat.2013.04.044}, author = {P. Epicoco and Benoit A. Coasne and A. Gioia and P. Papet and I. Cabodi and M. Gaubil} } @article {283, title = {Micromechanical analysis of second order work in granular media}, journal = {Granular Matter}, volume = {15}, year = {2013}, month = {Apr-2013}, pages = {221 - 235}, abstract = {

This paper examines instabilities in granular materials from a microscopic point of view through numerical simulations conducted using a discrete element method on two three-dimensional specimens. The detection and the tracking of grain scale deformation mechanisms constitute the key point for a better understanding the failure process and puzzling out what lies behind the vanishing of the macroscopic second order work. For this purpose, the second order work from microscopic variables, involving contact force and branch vector, was introduced and tracked numerically. Then, all contacts depicting negative values of the second order work were deeply investigated, especially their spatial distribution (homogeneity, agglomeration, dispersion\…) within the specimen according to the density of the granular assembly and to the loading direction. A set of comparisons has been considered in this context in order to highlight how a specimen is populated with such contacts whether it is loaded along a direction included within the plastic tensorial zone or along a direction for which the specimen is likely to behave elastically (elastic tensorial zone). Moreover, these comparisons concerned also loading directions within the cone of instability so that links between the vanishing of both microscopic and macroscopic second order works can be established and the local mechanisms responsible for failure occurrence may be figured out.

}, issn = {1434-5021}, doi = {10.1007/s10035-013-0402-3}, author = {Hadda, Nejib and Fran{\c c}ois Nicot and Bourrier, Franck and Sibille, Luc and Farhang Radja{\"\i} and Darve, F{\'e}lix} } @article {76, title = {Molecular intermittent dynamics of interfacial water: probing adsorption and bulk confinement}, journal = {Soft Matter}, volume = {9}, year = {2013}, month = {Jul-26-2013}, pages = {8654-8663}, type = {Article}}, abstract = {

Numerous natural and manufactured systems such as colloidal suspensions, geological pore networks, catalysts, and nanofluidic devices develop a large and sometimes complex interface strongly influencing the dynamics of the fluid entrapped inside these materials. A coarse grain picture of this molecular dynamics can be considered as an intermittence of adsorption steps and bulk relocations from one point to another point of the interface. Adsorption statistics such as the adsorption time distribution and its first moment reflect the degree of interaction of the molecule with the colloidal interface. Relocation statistics strongly depend on the shape of the pore, the surface forces and the bulk confinement. In this paper, a theoretical analysis of this intermittent dynamics is presented. A direct comparison with molecular dynamics simulations is proposed in the case of liquid water confined inside a hydrophilic substrate (silica slit pore with hydroxylated surfaces) or inside a hydrophobic substrate (carbon nanotube). Analysis of this intermittent dynamics allows quantification of the level of interaction of the vicinal water with the solid interface inside an independent adsorption region in exchange with the confined bulk fluid. The possibility of experimentally probing this dynamics using NMR relaxometry is emphasized.

}, issn = {1744-683X}, doi = {10.1039/c3sm51940f}, author = {Pierre E. Levitz and Patrick A Bonnaud and Pierre-Andre Cazade and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @article {103, title = {On the molecular origin of high-pressure effects in nanoconfinement: The role of surface chemistry and roughness}, journal = {Journal of Chemical Physics}, volume = {139}, year = {2013}, month = {Oct-14-2013}, type = {Article}, abstract = {

Experiments and simulations both suggest that the pressure experienced by an adsorbed phase confined within a carbon nanoporous material can be several orders of magnitude larger than the bulk phase pressure in equilibrium with the system. To investigate this pressure enhancement, we report a molecular-simulation study of the pressure tensor of argon confined in slit-shaped nanopores with walls of various models, including carbon and silica materials. We show that the pressure is strongly enhanced by confinement, arising from the effect of strongly attractive wall forces; confinement within purely repulsive walls does not lead to such enhanced pressures. Simulations with both the Lennard-Jones and Barker-Fisher-Watts intermolecular potentials for argon-argon interactions give rise to similar results. We also show that an increase in the wall roughness significantly decreases the in-pore pressure due to its influence on the structure of the adsorbate. Finally, we demonstrate that the pressures calculated from the mechanical (direct pressure tensor calculations) and the thermodynamic (volume perturbation method) routes yield almost identical results, suggesting that both methods can be used to calculate the local pressure tensor components in the case of these planar geometries. (C) 2013 AIP Publishing LLC.

}, issn = {0021-9606}, doi = {10.1063/1.4824125}, author = {Yun Long and Jeremy C. Palmer and Benoit A. Coasne and Malgorzata Sliwinska-Bartkowiak and George Jackson and Erich A. Mueller and Keith E. Gubbins} } @article {34, title = {Molecular simulation of adsorption and transport in hierarchical porous materials}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {May-29-2013}, pages = {7864{\textendash}7875}, abstract = {

\<\>Abstract Image

Adsorption and transport in hierarchical porous solids with micro- (\∼1 nm) and mesoporosities (\>2 nm) are investigated by molecular simulation. Two models of hierarchical solids are considered: microporous materials in which mesopores are carved out (model A) and mesoporous materials in which microporous nanoparticles are inserted (model B). Adsorption isotherms for model A can be described as a linear combination of the adsorption isotherms for pure mesoporous and microporous solids. In contrast, adsorption in model B departs from adsorption in pure microporous and mesoporous solids; the inserted microporous particles act as defects, which help nucleate the liquid phase within the mesopore and shift capillary condensation toward lower pressures. As far as transport under a pressure gradient is concerned, the flux in hierarchical materials consisting of microporous solids in which mesopores are carved out obeys the Navier\–Stokes equation so that Darcy\’s law is verified within the mesopore. Moreover, the flow in such materials is larger than in a single mesopore, due to the transfer between micropores and mesopores. This nonzero velocity at the mesopore surface implies that transport in such hierarchical materials involves slippage at the mesopore surface, although the adsorbate has a strong affinity for the surface. In contrast to model A, flux in model B is smaller than in a single mesopore, as the nanoparticles act as constrictions that hinder transport. By a subtle effect arising from fast transport in the mesopores, the presence of mesopores increases the number of molecules in the microporosity in hierarchical materials and, hence, decreases the flow in the micropores (due to mass conservation). As a result, we do not observe faster diffusion in the micropores of hierarchical materials upon flow but slower diffusion, which increases the contact time between the adsorbate and the surface of the microporosity.

}, doi = {10.1021/la401228k}, author = {Benoit A. Coasne and Anne Galarneau and Gerardin, Corine and Fajula, Fran{\c c}ois and Fran{\c c}ois Villemot} } @article {45, title = {Molecular simulations of supercritical fluid permeation through disordered microporous carbons}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {Aug-13-2013}, pages = {9985{\textendash}9990}, abstract = {
Abstract Image

Fluid transport through microporous carbon-based materials is inherent in numerous applications, ranging from gas separation by carbon molecular sieves to natural gas production from coal seams and gas shales. The present study investigates the steady-state permeation of supercritical methane in response to a constant cross-membrane pressure drop. We performed dual control volume grand canonical molecular dynamics (DCV-GCMD) simulations to mimic the conditions of actual permeation experiments. To overcome arbitrary assumptions regarding the investigated porous structures, the membranes were modeled after the CS1000a and CS1000 molecular models, which are representative of real microporous carbon materials. When adsorption-induced molecular trapping (AIMT) mechanisms are negligible, we show that the permeability of the microporous material, although not significantly sensitive to the pressure gradient, monotonically decreases with temperature and reservoir pressures, consistent with diffusion theory. However, when AIMT occurs, the permeability increases with temperature in agreement with experimental data found in the literature.

}, doi = {10.1021/la402087r}, author = {Botan, Alexandru and Vermorel, Romain and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {299, title = {Mechanistic Approach to Pavement-Vehicle Interaction and Its Impact on Life-Cycle Assessment}, journal = {Transportation Research Record: Journal of the Transportation Research Board}, year = {2012}, month = {Dec-2012}, pages = {171 - 179}, abstract = {

The accuracy and the comprehensiveness of any pavement life-cycle assessment are limited by the ability of the supporting science to quantify the environmental impact. Pavement-vehicle interaction represents a significant knowledge gap that has important implications for many pavement life-cycle assessment studies. In the current study, the authors assumed that a mechanistic model that linked pavement structure and properties to fuel consumption could contribute to closing the uncertainty gap of pavement-vehicle interaction in life-cycle assessment of pavements. The simplest mechanistic pavement model, a Bernoulli-Euler beam on a viscoelastic foundation subjected to a moving load, was considered. Wave propagation properties derived from falling weight deflectometer time history data of FHWA\&$\#$39;s Long-Term Pavement Performance program were used to calibrate top-layer and substrate moduli for various asphalt and concrete systems. The model was validated against recorded deflection data. The mechanistic response was used to determine gradient force and rolling resistance to link deflection to vehicle fuel consumption. A comparison with independent field data provided realistic order-of-magnitude estimates of fuel consumption related to pavement-vehicle interaction as predicted by the model.

}, issn = {0361-1981}, doi = {10.3141/2306-20}, author = {Mehdi Akbarian and Moeini-Ardakani, Seyed and Franz-Josef Ulm and Nazzal, Munir} } @proceedings {354, title = {Microscopic Origins of Shear Strength in Packings Composed of Elongated Particles}, journal = {9th International Workshop on Bifurcation and Degradation in Geomaterials (IWBDG 2011)}, volume = {Springer Series in Geomechanics and Geoengineering - ADVANCES IN BIFURCATION AND DEGRADATION IN GEOMATERIALS}, year = {2012}, month = {Apr-11-2012}, pages = {21-27}, address = {MAY 23-26 2011 Porquerolles, FRANCE}, abstract = {

We investigate the rheology, force transmission and texture of granular materials composed of elongated particles by means of contact dynamics simulations. The particles have a rounded-cap rectangular (RCR) shape described by a single elongation parameter varying from 0 for a circular particle to 1 for an infinitely thin or long particle. We study the quasi-static behavior, structural and force anisotropies as a function of the elongation parameter for packings submitted to biaxial compression. The shear strength is found to increase linearly with this parameter whereas the solid fraction both at the initial isotropic state and in the critical state is nonmonotonous. We show that for these elongated particles a harmonic decomposition of the stress tensor provides a fairly good approximation of the internal state. Our data suggest that the increase of shear strength with reflects both enhanced friction mobilization and anisotropic particle orientation as the elongation of the particles increases.

}, keywords = {Elongated particles, Fabric properties, force transmission, Harmonic decomposition}, url = {https://hal.archives-ouvertes.fr/hal-00686728}, author = {Emilien Az{\'e}ma and Farhang Radja{\"\i}}, editor = {St{\'e}phane Bonelli and Dascalu, Cristian and Fran{\c c}ois Nicot} } @proceedings {352, title = {Modeling Porous Granular Aggregates}, journal = {9th International Workshop on Buifurcation and Degradation in Geomaterials (IWBDG 2011)}, volume = {Springer Series in Geomechanics and Geoengineering - ADVANCES IN BIFURCATION AND DEGRADATION IN GEOMATERIALS}, year = {2011}, month = {May-28-2011}, pages = {249 - 254}, publisher = {Springer Netherlands}, address = {MAY 23-26 2011 Porquerolles, FRANCE}, abstract = {

We rely on 3D simulations based on the Lattice Element Method (LEM) to analyze the failure of porous granular aggregates under tensile loading. We investigate crack growth by considering the number of broken bonds in the particle phase as a function of the matrix volume fraction and particle-matrix adhesion. Three regimes are evidenced, corresponding to no particle damage, particle abrasion and particle fragmentation, respectively. We also show that the probability density of strong stresses falls off exponentially at high particle volume fractions where a percolating network of jammed particles occurs. Decreasing the matrix volume fraction leads to increasingly broader stress distribution and hence a higher stress concentration. Our findings are in agreement with 2D results previously reported in the literature.

}, isbn = {978-94-007-1420-5}, issn = {1866-8755}, doi = {10.1007/978-94-007-1421-2_32}, author = {Rafik Aff{\`e}s and Vincent Topin and Jean-Yves Delenne and Yann Monerie}, editor = {St{\'e}phane Bonelli and Dascalu, Cristian and Fran{\c c}ois Nicot} } @article {257, title = {Multitechnique Investigation of Extruded Clay Brick Microstructure}, journal = {Journal of the American Ceramic Society}, volume = {94}, year = {2011}, month = {Sept-2011}, pages = {3012 - 3022}, abstract = {

Despite the omnipresence of clay brick as construction material since thousands of years, fundamental knowledge about the link between composition, microstructure, and mechanical performance is still scarce. In this paper, we use a variety of advanced techniques of experimental mechanics and material characterization for extruded clay brick for masonry, that range from scanning electron microscopy (SEM) coupled with energy-dispersive X-ray spectroscopy, mercury intrusion porosimetry, to instrumented nanoindentation and macroscopic strength and durability tests. We find that extruded clay brick possesses a hierarchical microstructure: depending on the firing temperature, a \“glassy\” matrix phase, which manifests itself at submicrometer scales in form of neo-crystals of mullite, spinel-type phase, and other accessory minerals, forms either a granular or a continuum matrix phase that hosts at submillimeter scale the porosity. This porous composite forms the backbone for macroscopic material performance of extruded brick, including anisotropic strength, elasticity, and water absorption behavior.

}, doi = {10.1111/jace.2011.94.issue-910.1111/j.1551-2916.2011.04484.x}, author = {Konrad J. Krakowiak and Lourenco, Paulo B. and Franz-Josef Ulm} }