@article {673, title = {Bitumen@SiO2 core-shell particles green synthesis towards flowable powdered bitumen and their binder applications}, journal = {Colloids and Surfaces A: Physicochemical and Engineering Aspects}, volume = {570}, year = {2019}, month = {Mars-20-2019}, pages = {531 - 543}, abstract = {

Mineralization of occulated bitumen emulsions allows producing bitumen powders, also referred to as Bitumen@SiO2 hybrid particles, that can resist storage in hard conditions, with a load up to 8 kPa at 50 \°C during 20 h, without alteration of their morphology or release of their bitumen core. Thanks to the occulation, nal Bitumen@SiO2 powders are not highly pulverulent and allow transportation and handling of bitumen at ambient temperature without the need of water as it is the case for bitumen emulsions. This \“water free\” storage par- ticularity exhibits the double advantage of reducing the cost of transport and of avoiding inhibition of bitumen adhesion at mineral aggregates through the presence of water, therefore rendering the bitumen its high quality of binder. The silica shell is responsible for the high storage stability and, when mixed with aggregates for end use of asphalt, as for example in pavement applications, serves as mineral charge to the bitumen, modifying the bitumen ow properties. Indeed, when heated at 100 \°C and crushed, the capsules break releasing thus the bitumen. This suspension made of bitumen and silica shell fragments exhibits an interesting shear-thinning behavior that can be attributed to the presence of the remnant fragments. Indeed, when shell fragments have been transferred into a silicon oil, initially exhibiting a Newtonian behavior, the mixture of fragments and oil also exhibited a shear thinning rheogram demonstrating the e ect of the silica shell fragments. The bene t of such a behavior modi cation is the possibility to reduce or avoid the usual addition of llers. Passive adhesion has been tested by generating asphalt concretes with Bitumen@SiO2 powders and showed high stability.

}, issn = {09277757}, doi = {10.1016/j.colsurfa.2019.03.048}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0927775719302420}, author = {Anaclet, Pauline and Depardieu, Martin and Laurichesse, Eric and Julien, Viviane and Hung, Yvong and Schmitt, Veronique and Backov, Renal} } @article {665, title = {Cutting floating single-walled carbon nanotubes with a {\textquoteleft}CO2 blade{\textquoteright}}, journal = {Carbon}, volume = {143}, year = {2019}, month = {Mar-2019}, pages = {481 - 486}, abstract = {

Uniform ultra-short (\< 100 nm) single-walled carbon nanotubes (SWCNTs) with properties of rich chemistry, high dispersity and easy manipulation, are of great importance for constructing novel nanostructures and highly integrated miniature devices. Herein, we report a recent breakthrough in cutting long SWCNTs into uniform ultra-short segments in gas-phase using CO2 as a \&$\#$39;blade.\&$\#$39; SWCNTs synthesized by an aerosol reactor are directly introduced into a second reactor for gaseous cutting. The quality and the number concentration of shortened SWCNTs are enhanced after the gas-phase cutting process. Moreover, the growth and cutting of SWCNTs are accomplished in a continuous gas-phase process, thus allows direct dry deposition of ultra-short SWCNTs as individual macromolecules or thin films onto various substrates for multiple applications. (C) 2018 Elsevier Ltd. All rights reserved.

}, issn = {00086223}, doi = {10.1016/j.carbon.2018.11.035}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008622318310558}, author = {Tian, Ying and Wei, Nan and Laiho, Patrik and Ahmad, Saeed and Yann Magnin and Liao, Yongping and Christophe Bichara and Jiang, Hua and Kauppinen, Esko I.} } @article {669, title = {Fiber plucking by molecular motors yields large emergent contractility in stiff biopolymer networks}, journal = {Soft Matter}, volume = {15}, year = {2019}, month = {Feb-21-2020}, pages = {1481 - 1487}, abstract = {

The mechanical properties of the cell depend crucially on the tension of its cytoskeleton, a biopolymer network that is put under stress by active motor proteins. While the fibrous nature of the network is known to strongly affect the transmission of these forces to the cellular scale, our understanding of this process remains incomplete. Here we investigate the transmission of forces through the network at the individual filament level, and show that active forces can be geometrically amplified as a transverse motor-generated force \"plucks\&$\#$39;\&$\#$39; the fiber and induces a nonlinear tension. In stiff and densely connected networks, this tension results in large network-wide tensile stresses that far exceed the expectation drawn from a linear elastic theory. This amplification mechanism competes with a recently characterized network-level amplification due to fiber buckling, suggesting that that fiber networks provide several distinct pathways for living systems to amplify their molecular forces.

}, issn = {1744-683X}, doi = {10.1039/C8SM00979A}, url = {http://xlink.rsc.org/?DOI=C8SM00979Ahttp://pubs.rsc.org/en/content/articlepdf/2019/SM/C8SM00979A}, author = {Ronceray, Pierre and Broedersz, Chase P. and Lenz, Martin} } @article {668, title = {Spontaneous rotation can stabilise ordered chiral active fluids}, journal = {Nature Communications}, volume = {10}, year = {2019}, month = {Feb-22-2019}, pages = { Article Number: 920 }, abstract = {

Active hydrodynamic theories are a powerful tool to study the emergent ordered phases of internally driven particles such as bird flocks, bacterial suspension and their artificial analogues. While theories of orientationally ordered phases are by now well established, the effect of chirality on these phases is much less studied. In this paper, we present a complete dynamical theory of orientationally ordered chiral particles in two-dimensional incompressible systems. We show that phase-coherent states of rotating chiral particles are remarkably stable in both momentum-conserved and non-conserved systems in contrast to their non-rotating counterparts. Furthermore, defect separation-which drives chaotic flows in non-rotating active fluids-is suppressed by intrinsic rotation of chiral active particles. We thus establish chirality as a source of dramatic stabilisation in active systems, which could be key in interpreting the collective behaviors of some biological tissues, cytoskeletal systems and collections of bacteria.

}, doi = {10.1038/s41467-019-08914-7}, url = {http://www.nature.com/articles/s41467-019-08914-7}, author = {Maitra, Ananyo and Lenz, Martin} } @article {655, title = {Stress-dependent amplification of active forces in nonlinear elastic media}, journal = {Soft Matter}, volume = {15}, year = {2019}, month = {Jan-14-2019}, pages = {331 - 338}, abstract = {

The production of mechanical stresses in living organisms largely relies on localized, force-generating active units embedded in filamentous matrices. Numerical simulations of discrete fiber networks with fixed boundaries have shown that buckling in the matrix dramatically amplifies the resulting active stresses. Here we extend this result to a continuum elastic medium prone to buckling subjected to an arbitrary external stress, and derive analytical expressions for the active, nonlinear constitutive relations characterizing the full active medium. Inserting these relations into popular active gel descriptions of living tissues and the cytoskeleton will enable investigations into nonlinear regimes previously inaccessible due to the phenomenological nature of these theories.

}, issn = {1744-683X}, doi = {10.1039/C8SM00949J}, url = {http://xlink.rsc.org/?DOI=C8SM00949J}, author = {Ronceray, Pierre and Broedersz, Chase P. and Lenz, Martin} } @article {577, title = {Entropy-driven stability of chiral single-walled carbon nanotubes}, journal = {Science}, volume = {362}, year = {2018}, month = {Oct-12-2018}, pages = {212 - 215}, abstract = {

Single-walled carbon nanotubes are hollow cylinders that can grow centimeters long via carbon incorporation at the interface with a catalyst. They display semiconducting or metallic characteristics, depending on their helicity, which is determined during their growth. To support the quest for a selective synthesis, we develop a thermodynamic model that relates the tube-catalyst interfacial energies, temperature, and the resulting tube chirality. We show that nanotubes can grow chiral because of the configurational entropy of their nanometer-sized edge, thus explaining experimentally observed temperature evolutions of chiral distributions. Taking the chemical nature of the catalyst into account through interfacial energies, we derive structural maps and phase diagrams that will guide a rational choice of a catalyst and growth parameters toward a better selectivity.

}, issn = {0036-8075}, doi = {10.1126/science.aat6228}, url = {http://www.sciencemag.org/lookup/doi/10.1126/science.aat6228}, author = {Yann Magnin and Amara, Hakim and Ducastelle, Francois and Loiseau, Annick and Christophe Bichara} } @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 {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 {631, title = {A nonequilibrium force can stabilize 2D active nematics}, journal = {Proceedings of the National Academy of Sciences}, volume = {115}, year = {2018}, month = {Jul-03-2018}, pages = {6934-6939}, abstract = {

Suspensions of actively driven anisotropic objects exhibit distinctively nonequilibrium behaviors, and current theories predict that they are incapable of sustaining orientational order at high activity. By contrast, here we show that nematic suspensions on a substrate can display order at arbitrarily high activity due to a previously unreported, potentially stabilizing active force. This force moreover emerges inevitably in theories of active orientable fluids under geometric confinement. The resulting nonequilibrium ordered phase displays robust giant number fluctuations that cannot be suppressed even by an incompressible solvent. Our results apply to virtually all experimental assays used to investigate the active nematic ordering of self-propelled colloids, bacterial suspensions, and the cytoskeleton and have testable implications in interpreting their nonequilibrium behaviors.

}, issn = {0027-8424}, doi = {10.1073/pnas.1720607115}, url = {http://www.pnas.org/lookup/doi/10.1073/pnas.1720607115}, author = {Maitra, Ananyo and Srivastava, Pragya and Marchetti, M. Cristina and Lintuvuori, Juho S. and Ramaswamy, Sriram and Lenz, Martin} } @article {570, title = {Poroelasticity of Methane-Loaded Mature and Immature Kerogen from Molecular Simulations}, journal = {Langmuir}, volume = {34}, year = {2018}, month = {Nov-13-2018}, pages = {13766 - 13780}, abstract = {

While hydrocarbon expulsion from kerogen is certainly the key step in shale oil/gas recovery, the poromechanical couplings governing this desorption process, taking place under a significant pressure gradient, are still poorly understood. Especially, most molecular simulation investigations of hydrocarbon adsorption and transport in kerogen have so far been performed under the rigid matrix approximation, implying that the pore space is independent of pressure, temperature, and fluid loading, or in other words, neglecting poromechanics. Here, using two hydrogenated porous carbon models as proxies for immature and overmature kerogen, that is, highly aliphatic hydrogen-rich vs highly aromatic hydrogen-poor models, we perform an extensive molecular-dynamics-based investigation of the evolution of the poroelastic properties of those matrices with respect to temperature, external pressure, and methane loading as a prototype alkane molecule. The rigid matrix approximation is shown to hold reasonably well for overmature kerogen even though accounting for flexibility has allowed us to observe the well-known small volume contraction at low fluid loading and temperature. Our results demonstrate that immature kerogen is highly deformable. Within the ranges of conditions considered in this work, its density can double and its accessible porosity (to a methane molecule) can increase from 0 to \∼30\%. We also show that these deformations are significantly nonaffine (i.e., nonhomogeneous), especially upon fluid adsorption or desorption.

}, keywords = {ADSORPTION-INDUCED DEFORMATION; ORGANIC-MATTER; COMPETITIVE ADSORPTION; NANOPOROUS MEDIA; POROUS CARBONS; GLASSY-CARBON; COAL; TRANSPORT; SORPTION; MODEL}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.8b02534}, url = {http://pubs.acs.org/doi/10.1021/acs.langmuir.8b02534}, author = {Obliger, Ama{\"e}l and Valdenaire, Pierre-Louis and Capit, Nicolas and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Jean-Marc Leyssale} } @article {535, title = {Thermalizing and Damping in Structural Dynamics}, journal = {Journal of Applied Mechanics}, volume = {85}, year = {2018}, month = {Aug-2018}, abstract = {

Structural damping, that is the presence of a velocity dependent dissipative term in the equation of motion, is rationalized as a thermalization process between a structure (here a beam) and an outside bath (understood in a broad sense as a system property). This is achieved via the introduction of the kinetic temperature of structures and formalized by means of an extended Lagrangian formulation of a structure in contact with an outside bath at a given temperature. Using the Nose-Hoover thermostat, the heat exchange rate between structure and bath is identified as a mass damping coefficient, which evolves in time in function of the kinetic energy/temperature history exhibited by the structure. By way of application to a simple beam structure subjected to eigen-vibrations and dynamic buckling, commonality and differences of the Nose-Hoover beam theory with constant mass damping models are shown, which permit a handshake between classical damping models and statistical mechanics-based thermalization models. The solid foundation of these thermalization models in statistical physics provides new insights into stability and instability for engineering structures. Specifically, since two systems are considered in (thermodynamic) equilibrium when they have the same temperature, we show in the case of dynamic buckling that a persistent steady-state difference in kinetic temperature between structure and bath is but indicative of the instability of the system. This shows that the kinetic temperature can serve as a structural order parameter to identify and comprehend failure of structures, possibly well beyond the elastic stability considered here.

PDF FULL TEXT

jam_085_08_081001_f001.png

jam_085_08_081001_f002.png

jam_085_08_081001_f003.png

}, keywords = {structural dynamics; damping; Nose-Hoover Bath; kinetic temperature; dynamic buckling}, issn = {0021-8936}, doi = {10.1115/1.4040080}, url = {http://appliedmechanics.asmedigitalcollection.asme.org.libproxy.mit.edu/article.aspx?articleid=2680040}, author = {Arghavan Louhghalam and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {318, title = {Carbon management of infrastructure performance: Integrated big data analytics and pavement-vehicle-interactions}, journal = {Journal of Cleaner Production}, volume = {142}, year = {2017}, month = {Jan-20-2017}, pages = {956 - 964 Part: 2}, abstract = {

As a crucial part of the transportation system, roadway network provides mobility to the society and is vital for the economy. At the same time it contributes significantly to the environmental footprint during its construction, operation and maintenance. Hence, the sustainable development of our Nation\&$\#$39;s roadway system requires quantitative means to link infrastructure performance to lifecycle energy use and greenhouse gas emissions. Recent developments in mechanistic models of roughness- and deflection-induced pavement-vehicle interaction aim at providing such engineering estimates. Herein, it is demonstrated that these models when implemented at a network scale are a powerful basis for big data analytics of excess-energy consumption and carbon dioxide emissions by integrating spatially and temporally varying road conditions, pavement properties, traffic loads and climatic conditions. A novel ranking algorithm is proposed, that allows upscaling of the local carbon dioxide emissions due to pavement vehicle interaction to the size of state-wide or national sustainability goals. Implemented for 5157 lane-miles of the interstate highway system in the State of Virginia, sections contributing significantly to carbon dioxide emissions are identified. It is shown that the proposed ranking algorithm based on the inferred emission that exhibits a power-law distribution, provides the shortest path for greenhouse gas emissions savings per maintenance at network scale. That is, maintaining a few lane miles allows for a significant synergetic improvement of both infrastructure performance and environmental impact of the interstate network and helps transportation agencies in making economic and environmentally sustainable decisions.

Unlabelled figure

Fig.1. Flow of network analysis; the inputs and the output of deflection- and{\textellipsis}

Fig.2. Spatial distribution of roughness- and deflection-induced excess CO2 emissions{\textellipsis}

Fig.3. Total annual excess CO2 emissions (Tons) for a 7-year period (2007{\textendash}2013) in{\textellipsis}

Fig.4. Excess CO2 emissions in function of its rank

Fig.5. Comparison of different road selection strategies for maintenance in terms of{\textellipsis}

}, issn = {09596526}, doi = {10.1016/j.jclepro.2016.06.198}, author = {Arghavan Louhghalam and Mehdi Akbarian and Franz-Josef Ulm} } @article {600, title = {Cohesive strength of iron ore granules}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 08020}, abstract = {

We present an experimental and numerical investigation of the mechanical strength of crude iron ore (Hematite) granules in which capillary bonds between primary particles are the source of internal cohesion. The strength is measured by subjecting the granules to vertical compression between two plates. We show that the behavior of the granules is ductile with a well-defined plastic threshold which increases with the amount of water. It is found that the compressive strength scales with capillary cohesion with a pre-factor that is nearly independent of size polydispersity for the investigated range of parameters but increases with friction coefficient between primary particles. This weak dependence may be attributed to the class of fine particles which, due to their large number, behaves as a cohesive matrix that controls the strength of the granule.

}, doi = {10.1051/epjconf/201714008020}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714008020}, author = {Contreras, Rafael Jaimes and Berger, Nicolas and Izard, Edouard and Douce, Jean-Fran{\c c}ois and Koltsov, Alexey and Jean-Yves Delenne and Emilien Az{\'e}ma and Saeid Nezamabadi and van Loo, Fr{\'e}d{\'e}ric and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}}, editor = {Luding, S.} } @article {602, title = {Compaction of granular materials composed of deformable particles}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 05013}, abstract = {

In soft particle materials such as metallic powders the particles can undergo large deformations without rupture. The large elastic or plastic deformations of the particles are expected to strongly affect the mechanical properties of these materials compared to hard particle materials more often considered in research on granular materials. Herein, two numerical approaches are proposed for the simulation of soft granular systems: (i) an implicit formulation of the Material Point Method (MPM) combined with the Contact Dynamics (CD) method to deal with contact interactions, and (i) Bonded Particle Model (BPM), in which each deformable particle is modeled as an aggregate of rigid primary particles using the CD method. These two approaches allow us to simulate the compaction of an assembly of elastic or plastic particles. By analyzing the uniaxial compaction of 2D soft particle packings, we investigate the effects of particle shape change on the stress-strain relationship and volume change behavior as well as the evolution of the microstructure.

}, doi = {10.1051/epjconf/201714005013}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714005013}, author = {Thanh Hai Nguyen and Saeid Nezamabadi and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Luding, S.} } @article {239, title = {Disorder-induced stiffness degradation of highly disordered porous materials}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {106}, year = {2017}, month = {Sep-2017}, chapter = {207-228}, abstract = {

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

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

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

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

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

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

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

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

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

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

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

Fig. 11. Two-phase periodic porous solid

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

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

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

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

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

Fig. D.17. High porosity limit geometry

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

}, doi = {10.1016/j.jmps.2017.05.008}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {601, title = {Effect of particle size distribution on 3D packings of spherical particles}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 02030}, abstract = {

We use molecular dynamics simulations of frictionless spherical particles to investigate a class of polydisperse granular materials in which the particle size distribution is uniform in particle volumes. The particles are assembled in a box by uniaxial compaction under the action of a constant stress. Due to the absence of friction and the nature of size distribution, the generated packings have the highest packing fraction at a given size span, defined as the ratio of the largest size to the smallest size. We find that, up to = 5, the packing fraction is a nearly linear function of . While the coordination number is nearly constant due to the isostatic nature of the packings, we show that the connectivity of the particles evolves with . In particular, the proportion of particles with 4 contacts represents the largest proportion of particles mostly of small size. We argue that this particular class of particles occurs as a result of the high stability of local configurations in which a small particle is stuck by four larger particles.

}, doi = {10.1051/epjconf/201714002030}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714002030}, author = {Taiebat, Mahdi and Patrick Mutabaruka and Roland Jean-Marc Pellenq and Farhang Radja{\"\i}}, editor = {Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {591, title = {Effective Potentials and Elastic Properties in the Lattice-Element Method: Isotropy and Transverse Isotropy}, journal = {Journal of Nanomechanics and Micromechanics}, volume = {7}, year = {2017}, month = {Sep-2017}, pages = { Article Number: UNSP 04017007 }, abstract = {

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

}, issn = {2153-5434}, doi = {10.1061/(ASCE)NM.2153-5477.0000125}, url = {http://ascelibrary.org/doi/10.1061/\%28ASCE\%29NM.2153-5477.0000125}, author = {Hadrien Laubie and Monfared, Siavash and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {588, title = {From cellulose to kerogen: molecular simulation of a geological process}, journal = {Chemical Science}, volume = {8}, year = {2017}, month = {Dec-01-2017}, pages = {8325-8335}, abstract = {

The process by which organic matter decomposes deep underground to form petroleum and its underlying kerogen matrix has so far remained a no man\&$\#$39;s land to theoreticians, largely because of the geological (Myears) timescale associated with the process. Using reactive molecular dynamics and an accelerated simulation framework, the replica exchange molecular dynamics method, we simulate the full transformation of cellulose into kerogen and its associated fluid phase under prevailing geological conditions. We observe in sequence the fragmentation of the cellulose crystal and production of water, the development of an unsaturated aliphatic macromolecular phase and its aromatization. The composition of the solid residue along the maturation pathway strictly follows what is observed for natural type III kerogen and for artificially matured samples under confined conditions. After expulsion of the fluid phase, the obtained microporous kerogen possesses the structure, texture, density, porosity and stiffness observed for mature type III kerogen and a microporous carbon obtained by saccharose pyrolysis at low temperature. As expected for this variety of precursor, the main resulting hydrocarbon is methane. The present work thus demonstrates that molecular simulations can now be used to assess, almost quantitatively, such complex chemical processes as petrogenesis in fossil reservoirs and, more generally, the possible conversion of any natural product into bio-sourced materials and/or fuel.

}, issn = {2041-6520}, doi = {10.1039/C7SC03466K}, url = {http://xlink.rsc.org/?DOI=C7SC03466K}, author = {Atmani, Lea and Christophe Bichara and Roland Jean-Marc Pellenq and Henri Van Damme and Adri CT Van Duin and Raza, Zamaan and Truflandier, Lionel A. and Ama{\"e}l Obliger and Kralert, Paul G. and Franz-Josef Ulm and Jean-Marc Leyssale} } @article {595, title = {Heat-induced aging of agar solutions: Impact on the structural and mechanical properties of agar gels}, journal = {Food Hydrocolloids}, volume = {64}, year = {2017}, month = {Mar-2017}, pages = {59 - 69}, abstract = {

Thermoreversible gels are commonly prepared by cooling down to ambient temperature, aqueous polymer solutions first brought to a boil for a short time. Here, we study the effect of the latter duration on the subsequent gelation of an agar solution. Increasing the incubation time at high temperature (T= 80 degrees C) from hours to a few days leads to the hydrolysis and the oxidation of the polysaccharides, and results in the decrease of the viscosity and the pH of the solution. Samples withdrew at different incubation times are cooled down to form gels whose structure and mechanical properties are systematically determined. Cryoelectron microscopy and X-ray diffraction experiments reveal that agar gels formed from solutions of increasing incubation times, display a coarser microstructure composed of micron-sized foils which result from the condensation of the polysaccharides, and contrast with the fibrous-like microstructure of gels prepared from a fresh agar solution. Along with structural changes, a prolonged incubation time of the agar solution at T = 80 degrees C leads to weaker gels than those made from fresher solutions, and extensive macro-indentation experiments coupled to direct visualization show that the gel rupture scenario turns from brittle to ductile-like as the incubation time increases. Our study suggests that the incubation time of agar solutions at high temperature could be used as an external control parameter to tune the mechanical properties of agar-based gels. (C) 2016 Elsevier Ltd. All rights reserved.

}, keywords = {Agar; Aging; Gelation; Debonding; Indentation; Fracture}, issn = {0268005X}, doi = {10.1016/j.foodhyd.2016.10.020}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0268005X16305793}, author = {Mao, Bosi and Bentaleb, Ahmed and Louerat, Fr{\'e}d{\'e}ric and Divoux, Thibaut and Snabre, Patrick} } @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 {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 {623, title = {Nano-granular texture of cement hydrates}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 15027}, abstract = {

Mechanical behavior of concrete crucially depends on cement hydrates, the \“glue\” of cement. The design of high performance and more environmentally friendly cements demands a deeper understanding of the formation of the multiscale structure of cement hydrates, when they precipitate and densify. We investigate the precipitation and setting of nano-grains of cement hydrates using a combination of Monte Carlo and Molecular Dynamics numerical simulations and study their texture from nano up to the micron scale. We characterize the texture of cement hydrates using the local volume fraction distribution, the pore size distribution, the scattering intensity and the chord length distribution and we compare them with experiments. Our nano-granular model provides cement structure with realistic texture and mechanics and can be further used to investigate degradation mechanisms.

}, doi = {10.1051/epjconf/201714015027}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714015027}, author = {Katerina Ioannidou and Franz-Josef Ulm and Pierre E. Levitz and Emanuela Del Gado and Roland Jean-Marc Pellenq}, editor = {Farhang Radja{\"\i} and Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {604, title = {Numerical insight into the micromechanics of jet erosion of a cohesive granular material}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 15017}, abstract = {

Here we investigate the physical mechanisms behind the surface erosion of a cohesive granular soil induced by an impinging jet by means of numerical simulations coupling fluid and grains at the microscale. The 2D numerical model combines the Discrete Element and Lattice Boltzmann methods (DEM-LBM) and accounts for the granular cohesion with a contact model featuring a paraboloidal yield surface. Here we review first the hydrodynamical conditions imposed by the fluid jet on a solid granular packing, turning then the attention to the impact of cohesion on the erosion kinetics. Finally, the use of an additional subcritical debonding damage model based on the work of Silvani and co-workers provides a novel insight into the internal solicitation of the cohesive granular sample by the impinging jet.

}, doi = {10.1051/epjconf/201714015017}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714015017}, author = {Cuellar, Pablo and Benseghier, Zeyd and Luu, Li-Hua and St{\'e}phane Bonelli and Jean-Yves Delenne and Farhang Radja{\"\i} and Pierre Philippe}, editor = {Saeid Nezamabadi and Luding, S.} } @article {608, title = {Numerical modeling of the tensile strength of a biological granular aggregate: Effect of the particle size distribution}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 08013}, abstract = {

Wheat grains can be considered as a natural cemented granular material. They are milled under high forces to produce food products such as flour. The major part of the grain is the so-called starchy endosperm. It contains stiff starch granules, which show a multi-modal size distribution, and a softer protein matrix that surrounds the granules. Experimental milling studies and numerical simulations are going hand in hand to better understand the fragmentation behavior of this biological material and to improve milling performance. We present a numerical study of the effect of granule size distribution on the strength of such a cemented granular material. Samples of bi-modal starch granule size distribution were created and submitted to uniaxial tension, using a peridynamics method. We show that, when compared to the effects of starch-protein interface adhesion and voids, the granule size distribution has a limited effect on the samples\’ yield stress.

}, doi = {10.1051/epjconf/201714008013}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714008013}, author = {Heinze, Karsta and Frank, Xavier and Val{\'e}rie Lullien-Pellerina and George, Matthieu and Farhang Radja{\"\i} and Jean-Yves Delenne}, editor = {Saeid Nezamabadi and Luding, S.} } @article {603, title = {Numerical study of the failure of materials embedding soft to hard particles}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 02029}, abstract = {

In this study, we use a bond-based peridynamic approach to investigate the mechanical strength and cracking of composite materials with spherical inclusions. The total volume fraction of particles and the particle-matrix toughness ratio were varied to cover a range of soft to hard inclusions. The mean particle damage was characterized together with crack patterns at a sub-particle scale. Three types of crack patterns are identified depending on the toughness ratio.

}, doi = {10.1051/epjconf/201714002029}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714002029}, author = {Frank, Xavier and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Saeid Nezamabadi and Luding, S.} } @article {597, title = {Peridynamics simulation of the comminution of particles containing microcraks}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 07018}, abstract = {

In this study, we rely on a \’bond-based\’ peridynamic approach to investigate the strength and failure of 2D particles containing a collection of 1D microcracks. The mechanical tests were performed on disks under diametral compression. In an extensive parametric study, the distribution of microcracks was varied for different particle sizes. The evolution of yield stress with diameter and the probability of failure in terms of Weibull distributions are investigated in detail. Finally, by means of a floodfill algorithm, we analyze the variation of the mean fragment size as a function of the density of defects.

}, doi = {10.1051/epjconf/201714007018}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714007018}, author = {Blanc, Nicolas and Frank, Xavier and Mayer-Laigle, Claire and Farhang Radja{\"\i} and Jean-Yves Delenne}, editor = {Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {622, title = {The Potential of Mean Force concept for bridging (length and time) scales in the modeling of complex porous materials}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 01009}, abstract = {

We introduce the concept of Potential of Mean Force, PMF, as a way to implement upscaling modeling from the nano-scale to micron-scale. A PMF is a free energy function representing in an effective way the interactions between objects (cement hydrates, clay platelets, etc.) at thermodynamics conditions. The PMF is therefore the key piece of information allowing to coarse-grained Physical-Chemistry information in a meso-scale model formulation. The use of PMF offers a huge computational advantage as it allows a straight up-scaling to the meso-scale while keeping essential interactions information that are the hallmark of Physical-Chemistry processes. Such a coarse-grained modeling integrates atomistic response into inter-particle potentials that fully propagate molecular scale information all the way to the meso-scale.

}, doi = {10.1051/epjconf/201714001009}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714001009}, author = {Katerina Ioannidou and Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq}, editor = {Farhang Radja{\"\i} and Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {242, title = {A potential-of-mean-force approach for fracture mechanics of heterogeneous materials using the lattice element method}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {105}, year = {2017}, month = {Aug-2017}, pages = {116 - 130}, abstract = {

Fracture of heterogeneous materials has emerged as a critical issue in many engineering applications, ranging from subsurface energy to biomedical applications, and requires a rational framework that allows linking local fracture processes with global fracture descriptors such as the energy release rate, fracture energy and fracture toughness. This is achieved here by means of a local and a global potential-of-mean-force (PMF) inspired Lattice Element Method (LEM) approach. In the local approach, fracture-strength criteria derived from the effective interaction potentials between mass points are shown to exhibit a scaling commensurable with the energy dissipation of fracture processes. In the global PMF-approach, fracture is considered as a sequence of equilibrium states associated with minimum potential energy states analogous to Griffith\’s approach. It is found that this global approach has much in common with a Grand Canonical Monte Carlo (GCMC) approach, in which mass points are randomly removed following a maximum dissipation criterion until the energy release rate reaches the fracture energy. The duality of the two approaches is illustrated through the application of the PMF-inspired LEM for fracture propagation in a homogeneous linear elastic solid using different means of evaluating the energy release rate. Finally, by application of the method to a textbook example of fracture propagation in a heterogeneous material, it is shown that the proposed PMF-inspired LEM approach captures some well-known toughening mechanisms related to fracture energy contrast, elasticity contrast and crack deflection in the considered two-phase layered composite material.

}, issn = {00225096}, doi = {10.1016/j.jmps.2017.05.006}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {590, title = {Production of H 2 by water radiolysis in cement paste under electron irradiation: A joint experimental and theoretical study}, journal = {Cement and Concrete Research}, volume = {100}, year = {2017}, month = {Oct-2017}, pages = {110 - 118}, abstract = {

Long-term confinement of nuclear waste is one of the main challenges faced by the nuclear industry. Fission products such as Sr-90 and Cs-137, both beta(-) emitters known to induce serious health hazards, represent the largest fraction of nuclear waste. Cement is a good candidate to store them, provided it can resist the effects of irradiation over time. Here, we have investigated the effects of beta(-) decay on cement by performing electron irradiation experiments on different samples. We show that H-2 production in cement, the main effect of water radiolysis, depends strongly on composition and relative humidity. First-principles calculations indicate that the water-rich interlayer regions with Ca2+ ions act as electron traps that promote the formation of H-2. They also show that holes localize in water-rich regions in low Ca content samples and are then able to participate in H-2 production. This work provides new understanding of radiolysis effects in cements.

}, issn = {00088846}, doi = {10.1016/j.cemconres.2017.05.022}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008884617302065}, author = {Le Caer, Sophie and Dezerald, Lucile and Boukari, Khaoula and Laine, Maxime and Taupin, s{\'e}bastien and Kavanagh, Ryan M. and Johnston, Conrad S.N. and Foy, Eddy and Charpentier, Thibault and Konrad J. Krakowiak and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Tribello, Gareth A. and Kohanoff, Jorge J. and Andres Sa{\`u}l} } @article {607, title = {Scaling behavior of immersed granular flows}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 09044}, abstract = {

The shear behavior of granular materials immersed in a viscous fluid depends on fluid properties (viscosity, density), particle properties (size, density) and boundary conditions (shear rate, confining pressure). Using computational fluid dynamics simulations coupled with molecular dynamics for granular flow, and exploring a broad range of the values of parameters, we show that the parameter space can be reduced to a single parameter that controls the packing fraction and effective friction coefficient. This control parameter is a modified inertial number that incorporates viscous effects.

}, doi = {10.1051/epjconf/201714009044}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714009044}, author = {L. Amarsid and Jean-Yves Delenne and Patrick Mutabaruka and Yann Monerie and Perales, F. and Farhang Radja{\"\i}}, editor = {Saeid Nezamabadi and Luding, S.} } @article {605, title = {Small solar system bodies as granular systems}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 14011}, abstract = {

Asteroids and other Small Solar System Bodies (SSSBs) are currently of great scientific and even industrial interest. Asteroids exist as the permanent record of the formation of the Solar System and therefore hold many clues to its understanding as a whole, as well as insights into the formation of planetary bodies. Additionally, SSSBs are being investigated in the context of impact risks for the Earth, space situational awareness and their possible industrial exploitation (asteroid mining). In all these aspects, the knowledge of the geophysical characteristics of SSSB surface and internal structure are of great importance. Given their size, constitution, and the evidence that many SSSBs are not simple monoliths, these bodies should be studied and modelled as self-gravitating granular systems in general, or as granular systems in micro-gravity environments in particular contexts. As such, the study of the geophysical characteristics of SSSBs is a multi-disciplinary effort that lies at the crossroads between Granular Mechanics, Celestial Mechanics, Soil Mechanics, Aerospace Engineering and Computer Sciences.

}, doi = {10.1051/epjconf/201714014011}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714014011}, author = {Hestroffer, Daniel and Adriano Campo Bagat{\'\i}n and Losert, Wolfgang and Opsomer, Eric and S{\'a}nchez, Paul and Scheeres, Daniel J. and Staron, Lydie and Taberlet, Nicolas and Yano, Hajime and Eggl, Siegfried and Lecomte, Charles-Edouard and Murdoch, Naomi and Farhang Radja{\"\i} and Richardson, Derek C. and Salazar, Marcos and Schwartz, Stephen R. and Tanga, Paolo}, editor = {Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @book {636, title = {Springer Series in Geomechanics and GeoengineeringAdvances in Laboratory Testing and Modelling of Soils and Shales (ATMSS)Measurement of Mechanical Properties of Thin Clay Films and Comparison with Molecular Simulations}, series = {Springer Series in Geomechanics and Geoengineering}, year = {2017}, pages = {78-84}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

Here, we focus on the hydromechanical behavior of self-standing clay films with a thickness of a few dozen microns. We measure their elastic and creep properties and how those properties depend on the interlayer cation and on the relative humidity (or water content). Those experimental results are compared with the elastic and creep behavior of nanometric clay particles, which we characterize by molecular simulations. Significant qualitative differences between the behavior of the clay films and that of the clay particles are observed, which suggests that the hydromechanical behavior of the clay films is significantly impacted by their mesostructure (i.e., by how the clay particles or tactoids are arranged in space). Upscaling the hydromechanical behavior of the clay films from that of the clay particles may be challenging.

}, isbn = {978-3-319-52772-7}, issn = {1866-8755}, doi = {10.1007/978-3-319-52773-410.1007/978-3-319-52773-4_8}, url = {http://link.springer.com/10.1007/978-3-319-52773-4}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Henri Van Damme}, editor = {Ferrari, Alessio and Laloui, Lyesse} } @article {599, title = {Strength of wet agglomerates of spherical particles: effects of friction and size distribution}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 08021}, abstract = {

We investigate the mechanical behavior of wet granular agglomerates composed of spherical particles by means of molecular dynamics simulations. The capillary cohesion force is modeled as an attraction force at the contact between two particles and expressed as an explicit function of the gap and volume of the liquid bridge. We are interested in the effect of the friction coefficient between primary particles. The agglomerates are subjected to diametrical compression tests. We find that the deformation is ductile involving particle rearrangements. However, a well-defined stress peak is observed and the peak stress is used as a measure of the compressive strength of the agglomerate. The strength increases with friction coefficient but levels off at friction coefficients above 0.4. Furthermore, the compressive strength is an increasing function of particle size span.

}, doi = {10.1051/epjconf/201714008021}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714008021}, author = {Vo, Thanh-Trung and Patrick Mutabaruka and Jean-Yves Delenne and Saeid Nezamabadi and Farhang Radja{\"\i}}, editor = {Luding, S.} } @article {240, title = {Stress Transmission and Failure in Disordered Porous Media}, journal = {PHYSICAL REVIEW LETTERS}, volume = {119}, year = {2017}, month = {Aug-14-2017}, pages = {Article Number: 075501}, abstract = {

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

Figure

Figure

Figure

Figure

}, doi = {10.1103/PhysRevLett.119.075501}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and et al} } @article {593, title = {Stress Transmission and Failure in Disordered Porous Media}, journal = {Physical Review Letters}, volume = {119}, year = {2017}, month = {Aug-14-2017}, pages = { Article Number: 075501 }, abstract = {

By means of extensive lattice-element simulations, we investigate stress transmission and its relation with failure properties in increasingly disordered porous systems. We observe a non-Gaussian broadening of stress probability density functions under tensile loading with increasing porosity and disorder, revealing a gradual transition from a state governed by single-pore stress concentration to a state controlled by multipore interactions and metric disorder. This effect is captured by the excess kurtosis of stress distributions and shown to be nicely correlated with the second moment of local porosity fluctuations, which appears thus as a (dis)order parameter for the system. By generating statistical ensembles of porous textures with varying porosity and disorder, we derive a general expression for the fracture stress as a decreasing function of porosity and disorder. Focusing on critical sites where the local stress is above the global fracture threshold, we also analyze the transition to failure in terms of a coarse-graining length. These findings provide a general framework which can also be more generally applied to multiphase and structural heterogeneous materials.

}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.075501}, url = {https://link.aps.org/doi/10.1103/PhysRevLett.119.075501}, author = {Hadrien Laubie and Farhang Radja{\"\i} and Roland Jean-Marc Pellenq and Franz-Josef Ulm} } @article {609, title = {A time-dependent atomistic reconstruction of severe irradiation damage and associated property changes in nuclear graphite}, journal = {Carbon}, volume = {120}, year = {2017}, month = {Aug-2017}, pages = {111 - 120}, abstract = {

Wheat grains can be considered as a natural cemented granular material. They are milled under high forces to produce food products such as flour. The major part of the grain is the so-called starchy endosperm. It contains stiff starch granules, which show a multi-modal size distribution, and a softer protein matrix that surrounds the granules. Experimental milling studies and numerical simulations are going hand in hand to better understand the fragmentation behavior of this biological material and to improve milling performance. We present a numerical study of the effect of granule size distribution on the strength of such a cemented granular material. Samples of bi-modal starch granule size distribution were created and submitted to uniaxial tension, using a peridynamics method. We show that, when compared to the effects of starch-protein interface adhesion and voids, the granule size distribution has a limited effect on the samples\’ yield stress.

}, keywords = {TRANSMISSION ELECTRON-MICROSCOPY; MOLECULAR-DYNAMICS SIMULATIONS; HIGHLY ANISOTROPIC PYROCARBONS; THERMAL-CONDUCTIVITY; NEUTRON-IRRADIATION; CARBON; GRAPHENE; DEFECTS; ENERGY; NANOSTRUCTURE}, issn = {00086223}, doi = {10.1016/j.carbon.2017.05.009}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008622317304542}, author = {Baptiste Farbos and Freeman, Helen and Hardcastle, Trevor and Da Costa, Jean-Pierre and Brydson, Rik and Scott, Andrew J. and Weisbecker, Patrick and Germain Christian and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {598, title = {Wall roughness and nonlinear velocity profiles in granular shear flows}, journal = {EPJ Web of Conferences}, volume = {140}, year = {2017}, month = {Jun-30-2017}, pages = {Article Number 03090}, abstract = {

Inhomogeneous velocity profiles in granular flows are well known from both experiments and simulations, and considered as a hallmark of nonlocal behavior. By means of extensive contact dynamics simulations, we show that the sigmoidal velocity profiles in 2D flows of rigid disks are controlled by the roughness of driving boundary walls. We find that the velocity profile becomes linear for a critical value of wall roughness up to an exponential decay close to the walls with a characteristic length that does not depend on the flow thickness and rate. We describe the velocity profiles by introducing a state parameter that carries wall perturbation. By assuming that the local shear rate is a linear function of the state parameter, we obtain an analytical expression that fits velocity profiles. In this model, the nonlinear velocity profiles are explained in terms of the effects of wall roughness as boundary condition for the state parameter.

}, doi = {10.1051/epjconf/201714003090}, url = {http://www.epj-conferences.org/10.1051/epjconf/201714003090}, author = {Schuhmacher, Paul and Farhang Radja{\"\i} and St{\'e}phane Roux}, editor = {Saeid Nezamabadi and Luding, S. and Jean-Yves Delenne} } @article {639, title = {Activated desorption at heterogeneous interfaces and long-time kinetics of hydrocarbon recovery from nanoporous media}, journal = {Nature Communications}, volume = {7}, year = {2016}, month = {Jun-2016}, pages = { Article Number 11890 }, abstract = {

Hydrocarbon recovery from unconventional reservoirs (shale gas) is debated due to its environmental impact and uncertainties on its predictability. But a lack of scientific knowledge impedes the proposal of reliable alternatives. The requirement of hydrofracking, fast recovery decay and ultra-low permeability-inherent to their nanoporosity-are specificities of these reservoirs, which challenge existing frameworks. Here we use molecular simulation and statistical models to show that recovery is hampered by interfacial effects at the wet kerogen surface. Recovery is shown to be thermally activated with an energy barrier modelled from the interface wetting properties. We build a statistical model of the recovery kinetics with a two-regime decline that is consistent with published data: a short time decay, consistent with Darcy description, followed by a fast algebraic decay resulting from increasingly unreachable energy barriers. Replacing water by CO2 or propane eliminates the barriers, therefore raising hopes for clean/efficient recovery.

}, doi = {10.1038/ncomms11890}, url = {http://www.nature.com/articles/ncomms11890}, author = {Lee, Thomas and Lyd{\'e}ric Bocquet and Benoit A. Coasne} } @article {267, title = {Bottom-up model for understanding the effects of wheat endosperm microstructure on its mechanical strength}, journal = {Journal of Food Engineering}, volume = {190}, year = {2016}, month = {Dec-2016}, pages = {40 - 47}, abstract = {

Wheat flours are essential ingredients of daily food products like bread, cookies or pastries. Their quality depends on the milling process and mechanical strength of wheat grains. Although it is well known that the strength and rupture of grains are strongly controlled by the endosperm microstructure, the respective roles of the starch and polymer volume fractions and their adhesion are not yet fully understood. This typical biological microstructure can be modeled as a cemented granular material, where the two size populations of starch granules (large:A-type, small:B-type) are the particles, and the protein matrix, which partially fills the space between granules, plays the role of a cement. This structural model of wheat endosperm is used, together with mechanical characteristics of starch and proteins obtained by means of Atomic Force Microscopy (AFM) measurements, to simulate the mechanical behavior and breakage of wheat endosperm in milling process. We find that the porosity outweighs the effect of other parameters for the elastic modulus, which declines as a nearly linear function of porosity. We also show that the tensile strength is an increasing function of the amount and connectivity of starch granules with increasing concentration of stresses along chains of granules. This effect is more significant at low porosity where stress distribution is mainly controlled by the contact network between starch granules. This effect explains why the protein content is not fully correlated to vitreousness, and samples of similar protein content can be different in vitreosity. Finally, we find that the starch-granule adhesion strongly affects the tensile strength whereas the effect of starch volume fraction appears mainly at high interface adhesion, which is the case of hard type wheat grains.

Fig.1. (a) Scanning Electron Micrograph of the fracture surface of a soft wheat{\textellipsis}

Fig.2. Illustration of the lattice approach for one step: (1) Initial lattice state{\textellipsis}

Fig.3. (a) Bidisperse granular sample generated with DEM for the highest starch volume{\textellipsis}

Fig.4. Example of samples submitted to tensile tests: (a) Maps of vertical stress{\textellipsis}

Fig.5. Stress-strain plots showing the effect of different values of the starch{\textellipsis}

Fig.6. Effective Young{\textquoteright}s modulus as a function of (a) protein volume fraction ρp and{\textellipsis}

}, issn = {02608774}, doi = {10.1016/j.jfoodeng.2016.06.009}, author = {Chichti, Emna and Val{\'e}rie Lullien-Pellerina and George, Matthieu and Farhang Radja{\"\i} and Rafik Aff{\`e}s and Jean-Yves Delenne} } @proceedings {338, title = {Contribution of mechanical factors to the variability of root architecture: Quantifying the past history of interaction forces between growing roots and soil grains}, journal = {2016 IEEE INTERNATIONAL CONFERENCE ON FUNCTIONAL-STRUCTURAL PLANT GROWTH MODELING, SIMULATION, VISUALIZATION AND APPLICATIONS (FSPMA)}, year = {2016}, month = {Nov-11-2016}, pages = {52-60}, publisher = {I E E E}, address = {NOV 07-11, 2016, Qingdao, PEOPLES R CHINA }, abstract = {

The relation between a growing root and the soil movement has often been under-estimated. The present work aims to determine how grains in granular soils are reorganized by the action of growing roots, and in turn how the resulting forces acting on root tips modify their development. For this purpose, we have developed a 2D Discrete Element Model (DEM) able to compute a numerical growth of a single root inside a granular medium, taking into account the grain-grain and the root-grain contact forces during the growth. First in silico simulations were carried out in order to : 1-quantify the influence of the granular structure (grain diameter distribution and gaps) and root mechanical properties (root bending stiffness) on the evolution of reaction forces applied to a single root during its growth; 2-highlight \“group effects\”, e.g. how the reorganization of grains and their interaction forces due to a given growing root can affect the mechanical signal perceived by its near neighbours; 3-investigate how the presence of initial channels within the granular medium can effect the growth trajectory and minimize the resistance to penetration. All simulations were carried out assuming that root growth direction was only driven by external forces. Simlation results allowed the extraction of general physical laws that will be used further to provide mechanoperceptive indicators and analyze experimental data provided by phenotyping platforms. The final objective will be to quantify the response of plants to mechanical stresses in terms of root elongation rate, root straightness and ramification.

}, doi = {10.1109/FSPMA.2016.7818288}, author = {Fakih, Mahmoud and Jean-Yves Delenne and Farhang Radja{\"\i} and Fourcaud, Thierry} } @article {621, title = {The crucial effect of early-stage gelation on the mechanical properties of cement hydrates}, journal = {Nature Communications}, volume = {7}, year = {2016}, month = {Jul-15-2016}, pages = {Article number 12106}, abstract = {

Gelation and densification of calcium\–silicate\–hydrate take place during cement hydration. Both processes are crucial for the development of cement strength, and for the long-term evolution of concrete structures. However, the physicochemical environment evolves during cement formation, making it difficult to disentangle what factors are crucial for the mechanical properties. Here we use Monte Carlo and Molecular Dynamics simulations to study a coarse-grained model of cement formation, and investigate the equilibrium and arrested states. We can correlate the various structures with the time evolution of the interactions between the nano-hydrates during the preparation of cement. The novel emerging picture is that the changes of the physicochemical environment, which dictate the evolution of the effective interactions, specifically favour the early gel formation and its continuous densification. Our observations help us understand how cement attains its unique strength and may help in the rational design of the properties of cement and related materials.

}, doi = {10.1038/ncomms12106}, url = {http://www.nature.com/articles/ncomms12106}, author = {Katerina Ioannidou and Matej, Kandu{\v c} and Li, Lunna and Frenkel, Daan and Dobnikar, Jure and Emanuela Del Gado} } @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 {320, title = {Risk of Pavement Fracture due to Eigenstresses at Early Ages and Beyond}, journal = {Journal of Engineering Mechanics}, volume = {142}, year = {2016}, month = {Dec-2016}, pages = {Article Number: 04016105 }, abstract = {

Tensile cracks significantly affect the durability of concrete pavements, leading to an increase in the costs of maintenance and rehabilitation. A model is developed that relates thermal, chemical, and hygral evolutions at small scales due to different distress mechanisms to the risk of fracture at the structural scale. The method is based on application of linear elastic fracture mechanics (LEFM) to eigenstresses that develop in infinite- and finite-length beams on an elastic foundation that represents the subgrade. Axial and bending contributions to the energy release rate are determined for a worst-case scenario of an entirely cracked pavement section in functions of material properties, structural dimensions, and eigenstress forces and moments. By way of example, the model is used to study the risk of fracture of concrete pavements due to two different mechanisms: (1)\ autogeneous shrinkage at early ages of placing the concrete and (2)\ thermal cycles at the short term and long term after a temperature change. In addition, scaling relationships are developed that provide insight into the improvement of different structural and material properties for minimizing the risk of fracture.

}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0001167}, author = {Arghavan Louhghalam and Franz-Josef Ulm} } @article {136, title = {Capturing material toughness by molecular simulation: accounting for large yielding effects and limits}, journal = {International Journal of Fracture}, volume = {194}, year = {2015}, month = {Aug-2015}, pages = {149 - 167}, abstract = {

The inherent computational cost of molecular simulations limits their use to the study of nanometric systems with potentially strong size effects. In the case of fracture mechanics, size effects due to yielding at the crack tip can affect strongly the mechanical response of small systems. In this paper we consider two examples: a silica crystal for which yielding is limited to a few atoms at the crack tip, and a nanoporous polymer for which the process zone is about one order of magnitude larger. We perform molecular simulations of fracture of those materials and investigate in particular the system and crack size effects. The simulated systems are periodic with an initial crack. Quasi-static loading is achieved by increasing the system size in the direction orthogonal to the crack while maintaining a constant temperature. As expected, the behaviors of the two materials are significantly different. We show that the behavior of the silica crystal is reasonably well described by the classical framework of linear elastic fracture mechanics (LEFM). Therefore, one can easily upscale engineering fracture properties from molecular simulation results. In contrast, LEFM fails capturing the behavior of the polymer and we propose an alternative analysis based on cohesive crack zone models. We show that with a linear decreasing cohesive law, this alternative approach captures well the behavior of the polymer. Using this cohesive law, one can anticipate the mechanical behavior at larger scale and assess engineering fracture properties. Thus, despite the large yielding of the polymer at the scale of the molecular simulation, the cohesive zone analysis offers a proper upscaling methodology.

}, issn = {0376-9429}, doi = {10.1007/s10704-015-0045-y}, author = {Brochard, Laurent and Gy{\"o}rgy Hantal and Hadrien Laubie and Franz-Josef Ulm and Roland Jean-Marc Pellenq and Benoit A. Coasne} } @article {137, title = {Fracture toughness of calcium{\textendash}silicate{\textendash}hydrate from molecular dynamics simulations}, journal = {Journal of Non-Crystalline Solids}, volume = {419}, year = {2015}, month = {Jul-01-2015}, pages = {58 - 64}, abstract = {

Concrete is the most widely manufactured material in the world. Its binding phase, calcium\–silicate\–hydrate (C\–S\–H), is responsible for its mechanical properties and has an atomic structure fairly similar to that of usual calcium silicate glasses, which makes it appealing to study this material with tools and theories traditionally used for non-crystalline solids. Here, following this idea, we use molecular dynamics simulations to evaluate the fracture toughness of C\–S\–H, inaccessible experimentally. This allows us to discuss the brittleness of the material at the atomic scale. We show that, at this scale, C\–S\–H breaks in a ductile way, which prevents one from using methods based on linear elastic fracture mechanics. Knowledge of the fracture properties of C\–S\–H at the atomic scale opens the way for an upscaling approach to the design of tougher cement paste, which would allow for the design of slender environment-friendly infrastructures, requiring less material.

}, issn = {00223093}, doi = {10.1016/j.jnoncrysol.2015.03.031}, author = {Mathieu Bauchy and Hadrien Laubie and Mohammad Javad Abdolhosseini Qomi and Christian G. Hoover and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @proceedings {335, title = {Hydration Kinetics and Gel Morphology of C-S-H}, journal = {10th International Conference on Mechanics and Physics of Creep, Shrinkage, and Durability of Concrete and Concrete StructuresCONCREEP 10}, volume = {CONCREEP 10: MECHANICS AND PHYSICS OF CREEP, SHRINKAGE, AND DURABILITY OF CONCRETE AND CONCRETE STRUCTURES}, year = {2015}, month = {Sep-17-2015}, pages = {565-573}, publisher = {American Society of Civil Engineers}, address = {September 21{\textendash}23, 2015, Vienna, AustriaReston, VA}, abstract = {

Calcium-silicate hydrate (C-S-H) is the main binder in cement and concrete. It starts forming from the early stages of cement hydration and it progressively densifies as cement sets. C-S-H nanoscale building blocks form a cohesive gel, whose structure and mechanics are still poorly understood, in spite of its practical importance. Here we review a statistical physics approach recently developed, which allows us to investigate the C-S-H gel formation under the out-of-equilibrium conditions typical of cement hydration. Our approach is based on colloidal particles, precipitating in the pore solution and interacting with effective forces associated to the ionic environment. We present the evolution of the space filling of C-S-H with different particle interactions and compare them with experimental data at different lime concentrations. Moreover, we discuss the structural features of C-S-H in the mesoscale in terms of the scattering intensity. The comparison of our early stage C-S-H structures with small angle neutron scattering (SANS) experiments shows that long range spatial correlations and structural heterogeneties that develop in that early stages of hydration persist also in the hardened paste.

}, doi = {10.1061/978078447934610.1061/9780784479346.069}, author = {Katerina Ioannidou and Enrico Masoero and Pierre E. Levitz and Roland Jean-Marc Pellenq and Emanuela Del Gado}, editor = {Hellmich, Christian and Pichler, Bernhard and Kollegger, Johann} } @article {146, title = {Investigating carbon materials nanostructure using image orientation statistics}, journal = {Carbon}, volume = {84}, year = {2015}, month = {Apr-2015}, pages = {160 - 173}, abstract = {

A new characterization method of the lattice fringe images of turbostratic carbons is proposed. This method is based on the computation of their orientation field without explicit detection of fringes. It allows meaningful insights into the material nanostructure and nanotexture at several scales, either qualitatively or quantitatively. The calculation of pairwise spatial statistics of the orientation field at short distance provides measurements of the coherence lengths along any direction, in particular along and orthogonally to the layers. These statistics also allow representing orientation coherence patterns typical of the observed nanostructure. At larger distances, the mean disorientation of the fringes is computed and information about the homogeneity of the sample is obtained. An experimental validation is carried out on various artificial images and an application to the characterization of four bulk turbostratic carbons is provided.

}, issn = {00086223}, doi = {10.1016/j.carbon.2014.11.048}, author = {J. Philippe Da-Costa and Weisbecker, P. and Baptiste Farbos and Jean-Marc Leyssale and G{\'e}rard L. Vignoles and Germain, C.} } @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} } @article {143, title = {Nanoscale elasticity of highly anisotropic pyrocarbons}, journal = {Carbon}, volume = {94}, year = {2015}, month = {Nov-2015}, pages = {285 - 294}, abstract = {

We report on the elastic properties of high-textured laminar pyrocarbons (PyCs) as obtained at the nanoscale using atomistic simulations on realistic models of the rough laminar (RL) PyC, the regenerative laminar (ReL) PyC, and the ReL PyC heat treated at temperatures up to 1700 \°C. The purely longitudinal ( and ) and transverse ( and ) elastic properties of such materials have values of around 50\–75\% of those found for graphite. Conversely, cross longitudinal-transverse properties ( and ) are much larger in PyCs than in graphite (up to around one order of magnitude for ). Our results also show that stiffness decreases with the hydrogen content, that longitudinal properties increase with the extent of graphene domains and decrease with the misorientation of the layers, and that transverse and cross properties significantly increase with the amount of interlayer cross-links present in the materials. Comparison to known experimental data indicates that texture effect at superior scales play an important role in the effective macroscopic elasticity constants. Finally, from a materials perspective and in agreement with the structural evolution, the elastic properties of the ReL PyC evolve with heat treatment towards those of the RL PyC.

}, issn = {00086223}, doi = {10.1016/j.carbon.2015.06.060}, author = {Baptiste Farbos and J. Philippe Da-Costa and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {178, title = {Polarized Water Wires under Confinement in Chiral Channels}, journal = {The Journal of Physical Chemistry B}, volume = {119}, year = {2015}, month = {Jul-16-2015}, pages = {8707 - 8717}, abstract = {
Abstract Image

The alignment of water molecules along chiral pores may activate proton/ion conduction along dipolar hydrophilic pathways. Here we show that a simple synthetic \“T-channel\” forms a directional pore with its carbonyl moieties solvated by chiral helical water wires. Atom-scale simulations and experimental crystallographic assays reveal a dynamical structure of water and electrolyte solutions (alkali metal chlorides) confined in these organic T-channels. Oscillations in the dipole orientation, which correspond to alternative ordering (dipole up\–dipole down) of the water molecules with a period of about 4.2 \Å (imposed by the distance between two successive carbonyl groups) are observed. When ions are added to the system, despite the strong Coulombic water/ion interaction, confined water remains significantly ordered in the T-channel and still exhibits surface-induced polarization. Cation permeation can be achieved through alternated hydration\–dehydration occurring along strongly oriented water wires. The T-channel, which exhibits chirality with strong water orientation, provides an opportunity to unravel novel water-channel systems that share many interesting properties of biomolecular systems.

}, issn = {1520-6106}, doi = {10.1021/acs.jpcb.5b03322}, author = {Barboiu, Mihail and Pierre-Andre Cazade and Le Duc, Yann and Yves-Marie Legrand and Arie Van Der Lee and Benoit A. Coasne} } @article {181, title = {Predicting Adsorption on Bare and Modified Silica Surfaces}, journal = {The Journal of Physical Chemistry C}, volume = {119}, year = {2015}, month = {Mar-19-2015}, pages = {6009 - 6017}, abstract = {
Abstract Image

We show that Derjaguin\’s theory of adsorption can be used to predict adsorption on bare and modified surfaces using parameters available to simple experiments. Using experiment and molecular simulation of adsorption of various gases on hydroxylated, methylated, and trifluoromethylated silica, this simple parametrization of Derjaguin\’s model allows predicting adsorption on any functionalized surface using a minimum set of parameters such as the heat of vaporization of the adsorbate and the Henry constant of the adsorption isotherm. This general yet simple scheme constitutes a powerful tool as it avoids having to carry out tedious and complex adsorption measurements.

}, issn = {1932-7447}, doi = {10.1021/jp511726a}, author = {L{\'e}pinay, Matthieu and Broussous, Lucile and Licitra, Christophe and Bertin, Fran{\c c}ois and Rouessac, Vincent and Ayral, Andre and Benoit A. Coasne} } @article {144, title = {On the prediction of graphene{\textquoteright}s elastic properties with reactive empirical bond order potentials}, journal = {Carbon}, volume = {89}, year = {2015}, month = {Aug-2015}, pages = {176 - 187}, abstract = {

The elastic properties of graphene as described by the reactive empirical bond order potential are studied through uniaxial tensile tests calculations at both zero temperature, with a conjugate gradient approach, and room temperature, with molecular dynamics simulations. A perfect linear elastic behavior is observed at 0\ K up to \≈0.1\% strain. The Young\’s modulus and Poisson\’s ratio obtained with this potential are of \≈730\ GPa and 0.39, respectively, with little chirality effects. These values differ significantly from former estimations, much closer to experimental values. We show that these former values have certainly been obtained by neglecting the effect of atomic relaxation, leading to a severe inaccuracy. At larger strains, an extended apparent linear domain is observed in the stress\–strain curves, which is relevant to Young\’s modulus calculations at finite temperature. Our molecular dynamics simulations at 300\ K have allowed obtaining the following, chirality dependent, apparent Young\’s moduli, 860 and 761\ GPa, and Poisson\’s ratios, 0.12 and 0.23, for armchair and zigzag loadings, respectively.

}, issn = {00086223}, doi = {10.1016/j.carbon.2015.03.035}, author = {Antonio Gamboa and G{\'e}rard L. Vignoles and Jean-Marc Leyssale} } @article {174, title = {Probing the microporosity of low-k organosilica films: MP and t-plot methods applied to ellipsometric porosimetry data}, journal = {Microporous and Mesoporous Materials}, volume = {217}, year = {2015}, month = {Nov-15-2015}, pages = {119 - 124}, abstract = {

Ellipsometric porosimetry (EP) experiments are performed to obtain the adsorption isotherms of water, methanol, and toluene on pristine and damaged SiOCH porous materials. The use of gaseous adsorbates with different polarities, sizes, and surface tensions enables us to probe their affinity with such organosilica surfaces. Using reference t-curves obtained from Statistical Mechanics molecular simulations, we discuss the ability of the MicroPore analysis (MP) method to accurately estimate micropore sizes from EP measurements by comparing them with the mean pore sizes obtained using positron annihilation lifetime spectroscopy (PALS) and grazing incidence small angle X-ray scattering (GISAXS). We also report accessible microporous volumes estimated from the t-plot method used with the reference t-curves obtained by means of molecular simulation. We show that EP characterization of microporous films combined with the MP and t-plot methods can be improved by taking into account the effect of the chemical nature of the pore surface on the variation of the adsorbed thickness (t-curve).

Image for unlabelled figure

}, issn = {13871811}, doi = {10.1016/j.micromeso.2015.05.050}, author = {L{\'e}pinay, Matthieu and Broussous, Lucile and Licitra, Christophe and Bertin, Fran{\c c}ois and Rouessac, Vincent and Ayral, Andre and Benoit A. Coasne} } @article {147, title = {Roughness-Induced Vehicle Energy Dissipation: Statistical Analysis and Scaling}, journal = {Journal of Engineering Mechanics}, volume = {141}, year = {2015}, month = {Nov-2015}, pages = {Number: 04015046 }, abstract = {

The energy dissipated in a vehicle suspension system due to road roughness affects rolling resistance and the resulting fuel consumption and greenhouse gas emission. The key parameters driving this dissipation mechanism are identified via dimensional analysis. A mechanistic model is proposed that relates vehicle dynamic properties and road roughness statistics to vehicle dissipated energy and thus fuel consumption. A scaling relationship between the dissipated energy and the most commonly used road roughness index, the International Roughness Index (IRI), is also established. It is shown that the dissipated energy scales with IRI squared and scaling of dissipation with vehicle speed V depends on road waviness number w in the form of Vw\−2. The effect of marginal probability distribution of the road roughness profile on dissipated energy is examined. It is shown that although the marginal distribution of the road profile does not affect the identified scaling relationships, the multiplicative factor in these relationships does change from one distribution to another. As an example of practical application, the model is calibrated with the empirical HDM-4 model for different vehicle classes.

}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0000944}, author = {Arghavan Louhghalam and Tootkaboni, Mazdak and Franz-Josef Ulm} } @book {647, title = {Springer Series in Materials ScienceMolecular Dynamics Simulations of Disordered MaterialsFirst-Principles Modeling of Binary Chalcogenides: Recent Accomplishments and New Achievements}, volume = {215}, year = {2015}, pages = {313 - 344}, publisher = {Springer International Publishing}, organization = {Springer International Publishing}, address = {Cham}, abstract = {

This contribution is focussed on a set of first-principles molecular dynamics results obtained over the past fifteen years for disordered chalcogenides. In the first part, we sketch and review the historical premises underlying research efforts devoted to the understanding of structural properties in liquid and glassy GexSe1-x systems. We stress the importance of selecting well performing exchange-correlation functionals (within density functional theory) to achieve a correct description of short and intermediate range order. In the second part, we provide a specific, comparative example of structural analysis for chalcogenide GeX4 systems differing by the chemical identity of the X atom. We are able to demonstrate that the correct account of differences between the coordination environments of the two corresponding glasses requires system sizes substantially larger than similar to 100 atoms.

}, isbn = {978-3-319-15674-3}, issn = {0933-033X}, doi = {10.1007/978-3-319-15675-010.1007/978-3-319-15675-0_12}, url = {http://link.springer.com/10.1007/978-3-319-15675-0}, author = {Bouzid, A. and Le Roux, S. and Guido Ori and Tugene, Christine and Boero, M. and Carlo Massobrio}, editor = {Carlo Massobrio and Du, Jincheng and Bernasconi, Marco and Salmon, Philip S.} } @article {107, title = {An artificial primitive mimic of the Gramicidin-A channel}, journal = {Nature Communications}, volume = {5}, year = {2014}, month = {Jun-26-2014}, type = {Article}, abstract = {

Gramicidin A (gA) is the simplest known natural channel, and important progress in improving conduction activity has previously been obtained with modified natural gAs. However, simple artificial systems mimicking the gA functions are unknown. Here we show that gA can be mimicked using a simple synthetic triazole or {\textquoteleft}T-channel\&$\#$39; forming compound (TCT), having similar constitutional functions as the natural gAs. As in gA channels, the carbonyl moieties of the TCT, which point toward the T-channel core and surround the transport direction, are solvated by water. The net-dipolar alignment of water molecules along the chiral pore surfaces influences the conduction of protons/ions, envisioned to diffuse along dipolar hydrophilic pathways. Theoretical simulations and experimental assays reveal that the conduction through the T-channel, similar to that in gA, presents proton/water conduction, cation/anion selectivity and large open channel-conductance states. T-channels-associating supramolecular chirality with dipolar water alignment-represent an artificial primitive mimic of gA.

Molecular structures and crystal packing.

Molecular simulations of water molecules and ions confined within T-channel.

Bilayer membrane transport experiments.

Single channel conductance of T-channels under {\textquoteleft}early stage{\textquoteright} conductance records.

}, issn = {2041-1723}, doi = {10.1038/ncomms5142}, author = {Barboiu, Mihail and Le Duc, Yann and Gilles, Arnaud and Pierre-Andre Cazade and Michau, Mathieu and Yves-Marie Legrand and Arie Van Der Lee and Benoit A. Coasne and Parvizi, Paria and Post, Joshua and Fyles, Thomas} } @article {84, title = {Atomic-scale modelling of elastic and failure properties of clays}, journal = {Molecular Physics}, volume = {112}, year = {2014}, month = {May-19-2014}, pages = {1294-1305}, type = {Article}, abstract = {

The elastic and failure properties of a typical clay, illite, are investigated using molecular simulation. We employ a reactive (ReaxFF) and a non-reactive (ClayFF) force field to assess the elastic properties of the clay. As far as failure is concerned, ReaxFF was used throughout the study; however, some calculations were also performed with ClayFF. A crack parallel to the clay layers is found to have low fracture resistance when submitted to a tensile loading perpendicular to the crack. The mechanism of both yield and fracture failures is decohesion in the interlayer space. In contrast, under shear loading, the nanoscale failure mechanism is a stick-slip between clay layers. No fracture propagation is observed as the clay layers slide on top of each other. The low fracture resistance in mode I and the stick-slip failure in mode II are both the consequence of the lack of chemical bonds between clay layers where the cohesion is provided by non-covalent interactions. This work, which provides a description of the failure of clays at the microscopic scale, is the first step towards describing the failure of clays at a larger scale where the polycrystalline distribution of clay grains must be taken into account.

}, keywords = {clay, elastic properties, fracture, reactive molecular simulation}, issn = {0026-8976}, doi = {10.1080/00268976.2014.897393}, author = {Gy{\"o}rgy Hantal and Brochard, Laurent and Hadrien Laubie and Ebrahimi, Davoud and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @article {89, title = {Effect of Gas Adsorption on Acoustic Wave Propagation in MFI Zeolite Membrane Materials: Experiment and Molecular Simulation}, journal = {Langmuir}, volume = {30}, year = {2014}, month = {Aug-04-2014}, pages = {10336-10343}, type = {Article}, abstract = {
Abstract Image

The present study reports on the development of a characterization method of porous membrane materials which consists of considering their acoustic properties upon gas adsorption. Using acoustic microscopy experiments and atomistic molecular simulations for helium adsorbed in a silicalite-1 zeolite membrane layer, we showed that acoustic wave propagation could be used, in principle, for controlling the membranes operando. Molecular simulations, which were found to fit experimental data, showed that the compressional modulus of the composite system consisting of silicalite-1 with adsorbed He increases linearly with the He adsorbed amount while its shear modulus remains constant in a large range of applied pressures. These results suggest that the longitudinal and Rayleigh wave velocities (VL and VR) depend on the He adsorbed amount whereas the transverse wave velocity VT remains constant.

}, issn = {0743-7463}, doi = {10.1021/la502182k}, author = {Manga, Etoungh D. and Blasco, Hugues and J. Philippe Da-Costa and Drobek, Martin and Ayral, Andre and Le Clezio, Emmanuel and Despaux, Gilles and Benoit A. Coasne and Julbe, Anne} } @article {313, title = {Fl{\"u}gge{\textquoteright}s Conjecture: Dissipation- versus Deflection-Induced Pavement{\textendash}Vehicle Interactions}, journal = {Journal of Engineering Mechanics}, volume = {140}, year = {2014}, month = {Aug-2014}, pages = {Article Number: 04014053 }, abstract = {

The dissipation occurring below a moving tire in steady-state conditions in contact with a viscoelastic pavement is expressed using two different reference frames: a fixed observer attached to the pavement and a moving observer attached to the pavement\–tire contact surface. The first approach is commonly referred to as dissipation-induced pavement\–vehicle interaction (PVI), the second as deflection-induced PVI. Based on the principle of frame independence, it is shown that both approaches are strictly equal, from a thermodynamic point of view, and thus predict the same amount of dissipated energy. This equivalence is illustrated through application to two pavement systems: a viscoelastic beam and a viscoelastic plate both resting on an elastic foundation. The amount of dissipated energy in the pavement structure needs to be supplied by the vehicle to maintain constant speed, thus contributing to the rolling resistance, associated excess fuel consumption, and greenhouse gas emissions. The model here proposed can be used to quantify the dissipated energy and contribute to the development of engineering methods for the sustainable design of pavements.

}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0000754}, author = {Arghavan Louhghalam and Mehdi Akbarian and Franz-Josef Ulm} } @article {312, title = {Irwin׳s conjecture: Crack shape adaptability in transversely isotropic solids}, journal = {Journal of the Mechanics and Physics of Solids}, volume = {68}, year = {2014}, month = {Aug-2014}, pages = {1 - 13}, abstract = {

The planar crack propagation problem of a flat elliptical crack embedded in a brittle elastic anisotropic solid is investigated. We introduce the concept of crack shape adaptability: the ability of three-dimensional planar cracks to shape with the mechanical properties of a cracked body. A criterion based on the principle of maximum dissipation is suggested in order to determine the most stable elliptical shape. This criterion is applied to the specific case of vertical cracks in transversely isotropic solids. It is shown that contrary to the isotropic case, the circular shape (i.e. penny-shaped cracks) is not the most stable one. Upon propagation, the crack first grows non-self-similarly before it reaches a stable shape. This stable shape can be approximated by an ellipse of an aspect ratio that varies with the degree of elastic anisotropy. By way of example, we apply the so-derived crack shape adaptability criterion to shale materials. For this class of materials it is shown that once the stable shape is reached, the crack propagates at a higher rate in the horizontal direction than in the vertical direction. We also comment on the possible implications of these findings for hydraulic fracturing operations.

Fig.1. Top view of an elliptical crack (en, et and eb are the unit vectors normal to{\textellipsis}

Fig.2. Definition of the θ angle for a penny-shaped crack

Fig.3. Definition of the α-angle for an elliptical vertical crack (θ=π/2)

Fig.4. Zoom at the crack tip (vertical crack)

Fig.5. Comparison of the H(ϕ)/H2 obtained numerically {\textendash} solid line {\textendash} and its{\textellipsis}

Fig.6. Comparison of the H(ϕ)/H2 obtained numerically {\textendash} solid line {\textendash} and its{\textellipsis}

}, issn = {00225096}, doi = {10.1016/j.jmps.2014.03.004}, author = {Hadrien Laubie and Franz-Josef Ulm} } @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 {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} } @proceedings {348, title = {Pavement infrastructures footprint: The impact of pavement properties on vehicle fuel consumption}, journal = {Euro-C Conference}, volume = {COMPUTATIONAL MODELLING OF CONCRETE STRUCTURES}, year = {2014}, month = {Mar-2014}, pages = {1051-1058}, address = {MAR 24-27 2014 St Anton am Alberg, AUSTRIA}, abstract = {

A novel mechanistic model based on an infinite beam on elastic foundation is developed to
quantify the impact of pavement structural and material properties on pavement deflection and consequently on
vehicle fuel consumption. The model can also account for the effect of temperature and vehicle speed on fuel
consumption. A simplified expression for evaluating the energy dissipation for practical purposes is proposed
and used to investigate the impact of various pavement design systems on fuel consumption. GPS (General Pave-
ment Studies) sections from the FHWA\’s Long Term Pavement Performance program (FHWA 2011) are used
for this study. These sections consist of asphalt concrete (AC), portland cement concrete (PCC) and composite
pavements. The model quantifies the impact of temperature and vehicle speed on the fuel consumption and
confirms that those impacts are negligible for PCC and significant for AC pavements due to their viscoelasticity.

}, url = {https://cshub.mit.edu/sites/default/files/documents/louhghalam-Euro-C-2014.pdf}, author = {Arghavan Louhghalam and Mehdi Akbarian and Franz-Josef Ulm}, editor = {Bicanic, N and Mang, H and Meschke, G} } @article {327, title = {Plane-Strain Crack Problem in Transversely Isotropic Solids for Hydraulic Fracturing Applications}, journal = {Journal of Engineering Mechanics}, volume = {140}, year = {2014}, month = {Dec-2014}, pages = {Article Number: 04014092 }, abstract = {

This paper aims at understanding and predicting how pressurized cracks propagate in anisotropic brittle solids, a situation frequently encountered in hydraulic fracturing. Special attention is paid to transverse isotropy, often used to model shale. Although the theory of linear elastic fracture mechanics of anisotropic solids is well established at present, this paper shows that the application of Muskhelishvili\’s formalism to Lekhnitskii\’s anisotropic complex potentials provides a powerful tool to extend the validity of the classical tools of isotropic fluid-driven crack models to the anisotropic case, provided that the appropriate elastic constants are used. These elastic constants are identified and derived in closed form for transversely isotropic solids. The constants are shown to be directly related to quantities easily measured in a laboratory at macroscopic scale through indentation tests and acoustic measurements. Moreover, several crack-kinking criteria are compared. Contrary to the isotropic case, the crack-kinking criteria are not consistent among themselves, even in the case of a pure pressure loading. The orientation at which it is easier to propagate an already existing crack is sought. A critical crack length, below which this crack orientation is the one of minimal stiffness felt by the crack, is identified.

}, issn = {0733-9399}, doi = {10.1061/(ASCE)EM.1943-7889.0000807}, author = {Hadrien Laubie and Franz-Josef Ulm} } @article {87, title = {Pressure effects in confined nanophases}, journal = {Molecular Simulation}, volume = {40}, year = {2014}, month = {Feb-05-2014}, pages = {721-730}, type = {Article}, abstract = {

In this article, we review how pressure effects in pores affect both the physics of the confined fluid and the properties of the host porous material. Molecular simulations in which high-pressure effects were observed are first discussed; we will see how the strong dependence on bulk phase pressure of the freezing temperature of a fluid confined in nanopores can be explained by important variations of the pressure within the pore. We then discuss recent works in which direct calculations of the pressure tensor of fluids confined in pores provide evidence for large pressure enhancements. Finally, practical applications of these pressure effects in which gas adsorption in microporous solids (pore size \<2nm) was found to enhance their mechanical properties by increasing the elastic modulus by a factor 4 are discussed.

}, keywords = {confinement, high-pressure effects, mechanical properties, porous materials}, issn = {0892-7022}, doi = {10.1080/08927022.2013.829227}, author = {Benoit A. Coasne and Yun Long and Keith E. Gubbins} } @article {149, title = {Roughness-induced Pavement-Vehicle Interactions Key Parameters and Impact on Vehicle Fuel Consumption}, journal = {Transportation Research Record: Journal of the Transportation Research Board}, volume = {2525}, year = {2014}, month = {Dec-30-2014}, pages = {62 - 70}, abstract = {

Pavement roughness affects rolling resistance and thus vehicle fuel consumption. When a vehicle travels at constant speed on an uneven road surface, the mechanical work dissipated in the vehicle\’s suspension system is compensated by vehicle engine power and results in excess fuel consumption. This dissipation depends on both road roughness and vehicle dynamic characteristics. This paper proposes, calibrates, and implements a mechanistic model for roughness-induced dissipation. The distinguishing feature of the model is its combination of a thermodynamic quantity (energy dissipation) with results from random vibration theory to identify the governing parameters that drive the excess fuel consumption caused by pavement roughness, namely, the international roughness index (IRI) and the waviness number, w (a power spectral density parameter). It is shown through sensitivity analysis that the sensitivity of model output, that is, excess fuel consumption, to the waviness number is significant and comparable to that of IRI. Thus, introducing the waviness number as a second roughness index, in addition to IRI, allows a more accurate quantification of the impact of surface characteristics on vehicle fuel consumption and the corresponding greenhouse gas emissions. This aspect is illustrated by application of the roughness\–fuel consumption model to two road profiles extracted from FHWA\’s Long-Term Pavement Performance database.

}, issn = {0361-1981}, doi = {10.3141/2525-07}, author = {Arghavan Louhghalam and Mehdi Akbarian and Franz-Josef Ulm} } @article {304, title = {Scaling Relationships of Dissipation-Induced Pavement-Vehicle Interactions}, journal = {Transportation Research Record: Journal of the Transportation Research Board}, volume = {2457}, year = {2014}, month = {Jan-2014}, pages = {95 - 104}, abstract = {

Rolling resistance is one of the key factors that affect the fuel efficiency of the national pavement system. In addition to pavement texture and pavement roughness, the dissipation of mechanical work provided by the vehicle because of viscous deformation within the pavement structure has been recognized as a relevant factor contributing to the environmental footprint of pavement systems. This dissipation depends on material and structural parameters that can be optimized to increase the fuel efficiency of pavements. Identifying the key material and structural parameters that drive this dissipation is the focus of this paper. This identification is achieved by a combination of dimensional analysis and model-based simulations of the dissipation of a viscoelastic beam on an elastic foundation. For linear viscoelastic systems, the dissipation is found to scale with the square of the vehicle weight and with the inverse of the viscous relaxation time, in addition to distinct power relations of top-layer stiffness, thickness, and subgrade modulus. These scaling relations can be used by pavement engineers to reduce such pavement-inherent dissipation mechanisms and increase the fuel efficiency of a pavement design. An example shows the application of these scaling relations with data extracted from FHWA\&$\#$39;s Long-Term Pavement Performance database for seven road classes. The scaling relations provide a means for evaluating the performance of the various road classes in terms of the fuel efficiency related to dissipation.

}, issn = {0361-1981}, doi = {10.3141/2457-10}, author = {Arghavan Louhghalam and Mehdi Akbarian and Franz-Josef Ulm} } @article {79, title = {Structure-property relationships of water adsorption in metal-organic frameworks}, journal = {New Journal of Chemistry}, volume = {38}, year = {2014}, month = {Apr-16-2014}, pages = {3102-3111}, type = {Article}, abstract = {

A set of 15 metal\–organic frameworks (MIL-53, MIL-68, MIL-125, UiO-66, ZIF) exhibiting different pore size, morphology, and surface chemistry is used to unravel the numerous behaviors of water adsorption at room temperature in this class of materials. Outstanding \“S\”-shaped (type V) adsorption isotherms are observed for MIL-68 type solids. We show that the underlying mechanism of water adsorption can be rationalized using a simple set of three parameters: the Henry constant (i.e. the slope of the adsorption pressure in the low pressure range), the pressure at which pore filling occurs, and the maximum water adsorption capacity. While the Henry constant and pore filling pressure mostly depend on the affinity of water for the surface chemistry and on pore size, respectively, these two parameters are correlated as they both reflect different aspects of the hydrophobicity\–hydrophilicity of the material. For a given type of porous structure, the functionalization of the material by hydrophilic moieties such as hydrogen bonding groups (amine or aldehyde) systematically leads to an increase in the Henry constant concomitantly with a decrease in the pore filling pressure. As for the adsorption mechanism, we show that, for a given temperature, there is a critical diameter (Dc \∼ 20 \Å for water at room temperature) above which pore filling occurs through irreversible capillary condensation accompanied by capillary hysteresis loops. Below this critical diameter, pore filling is continuous and reversible unless the material exhibits some adsorption-induced flexibility.

Graphical abstract: Structure{\textendash}property relationships of water adsorption in metal{\textendash}organic frameworks
}, issn = {1144-0546}, doi = {10.1039/c4nj00076e}, author = {Canivet, Jerome and Bonnefoy, Jonathan and Daniel, Cecile and Legrand, Alexandre and Benoit A. Coasne and David Farrusseng} } @article {77, title = {First-principles molecular dynamics study of glassy GeS2: Atomic structure and bonding properties}, journal = {Physical Review B}, volume = {88}, year = {2013}, month = {Nov-04-2013}, type = {Article}, abstract = {

The structure of glassy GeS2 is studied in the framework of density functional theory, by using a fully self-consistent first-principles molecular dynamics (FPMD) scheme. A comparative analysis is performed with previous molecular dynamics data obtained within the Harris functional (HFMD) total energy approach. The calculated total neutron structure factor exhibits an unprecedented agreement with the experimental counterpart. In particular, the height of the first sharp diffraction peak (FSDP) improves considerably upon the HFMD results. Both the Ge and the S subnetworks are affected by a consistent number of miscoordinations, coexisting with the main tetrahedral structural motif. Glassy GeS2 features a short-range order quite similar to the one found in glassy GeSe2, a notable exception being the larger number of edge-sharing connections. An electronic structure localization analysis, based on the Wannier functions formalism, provides evidence of a more enhanced ionic character in glassy GeS2 when compared to glassy GeS2.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.88.174201}, author = {Celino, M. and Le Roux, S. and Guido Ori and Benoit A. Coasne and Bouzid, A. and Boero, M. and Carlo Massobrio} } @article {284, title = {An interdisciplinary approach towards improved understanding of soil deformation during compaction}, journal = {Soil and Tillage Research}, volume = {128}, year = {2013}, month = {Apr-2013}, pages = {61 - 80}, abstract = {

Soil compaction not only reduces available pore volume in which fluids are stored, but it alters the arrangement of soil constituents and pore geometry, thereby adversely impacting fluid transport and a range of soil ecological functions. Quantitative understanding of stress transmission and deformation processes in arable soils remains limited. Yet such knowledge is essential for better predictions of effects of soil management practices such as agricultural field traffic on soil functioning. Concepts and theory used in agricultural soil mechanics (soil compaction and soil tillage) are often adopted from conventional soil mechanics (e.g. foundation engineering). However, in contrast with standard geotechnical applications, undesired stresses applied by agricultural tyres/tracks are highly dynamic and last for very short times. Moreover, arable soils are typically unsaturated and contain important secondary structures (e.g. aggregates), factors important for affecting their soil mechanical behaviour. Mechanical processes in porous media are not only of concern in soil mechanics, but also in other fields including geophysics and granular material science. Despite similarity of basic mechanical processes, theoretical frameworks often differ and reflect disciplinary focus. We review concepts from different but complementary fields concerned with porous media mechanics and highlight opportunities for synergistic advances in understanding deformation and compaction of arable soils. We highlight the important role of technological advances in non-destructive measurement methods at pore (X-ray tomography) and soil profile (seismic) scales that not only offer new insights into soil architecture and enable visualization of soil deformation, but are becoming instrumental in the development and validation of new soil compaction models. The integration of concepts underlying dynamic processes that modify soil pore spaces and bulk properties will improve the understanding of how soil management affect vital soil mechanical, hydraulic and ecological functions supporting plant growth.

}, issn = {01671987}, doi = {10.1016/j.still.2012.10.004}, author = {Keller, T. and Lamand, M. and Peth, S. and Berli, M. and Jean-Yves Delenne and Baumgarten, W. and Rabbel, W. and Farhang Radja{\"\i} and Rajchenbach, J. and Selvadurai, A.P.S. and Or, D.} } @article {25, title = {Large pressure enhancements in confined phases: effect of surface chemistry and roughness}, journal = {J. Chem. Phys.}, volume = {139}, year = {2013}, month = {Oct-10-2013}, chapter = {144701}, 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.

Image of FIG. 1.
Image of FIG. 2.
Image of FIG. 3.
Image of FIG. 4.
}, 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 {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 {279, title = {Nano-mechanical properties of starch and gluten biopolymers from atomic force microscopy}, journal = {European Polymer Journal}, volume = {49}, year = {2013}, month = {Dec-2013}, pages = {3788 - 3795}, abstract = {

An original method based on atomic force microscopy (AFM) in contact mode was developed to abrade progressively the surface of tablets made of starch or gluten polymers isolated from wheat. The volume of the material removed by the tip was estimated from the analysis of successive topographic images of the surface, and the shear force was measured by keeping a constant normal force. Our data together with a simple tribological model provide clear evidence for a higher hardness and shear strength of starch compared to gluten. Gluten appears to have mechanical properties close to soft materials, such as talc, whereas starch displays higher hardness close to calcite. Our results are in a better agreement with structural properties of gluten (complex protein network) and starch (granular and semi-cristalline structure) than earlier studies by micro-indentation. This work shows that the AFM scratching method is relevant for the characterization of any polymer surface, in particular in application to materials made of different polymers at the nano-scale.

Graphical abstract

Unlabelled figure
}, issn = {00143057}, doi = {10.1016/j.eurpolymj.2013.08.024}, author = {Chichti, Emna and George, Matthieu and Jean-Yves Delenne and Farhang Radja{\"\i} and Val{\'e}rie Lullien-Pellerina} } @article {285, title = {Quasistatic rheology and microstructural description of sheared granular materials composed of platy particles}, journal = {Physical Review E}, volume = {87}, year = {2013}, month = {Mar-28-2013}, pages = {Article Number: 032206}, abstract = {

This is the first paper of a series devoted to the micro-mechanical modeling of clayey soils, by means of discrete element simulations. We specifically focus here on the effect of the platy shape of particles by reducing the interactions between particles to mechanical contact forces (i.e., neither electrostatic repulsion nor van der Waals forces are taken into account). The particles are three-dimensional square plates, approximated as spheropolyhedra. Several samples composed of particles of different levels of platyness (related to the ratio of length to thickness) were numerically prepared and sheared up to large deformations. We analyzed the shear strength, packing fraction, orientation of the particles, connectivity, fabric of the interactions network, and interaction forces as functions of the platyness. We find that both the mechanical behavior and microstructure are strongly dependent on the degree of platyness. The principal underlying phenomenon is the alignment of particle faces along a particular direction. This ordering phenomenon, which emerges even for shapes that deviate only slightly from that of a sphere, enhances the ability of the packing to develop an anisotropic structure leading to large shear strength, especially as a consequence of the fabric and mobilization of friction forces. Moreover, the connectivity of the packings and their packing fraction also evolve with the platyness. In particular, the packing fraction evolves in a nonmonotonic fashion, as observed in other granular materials composed of elongated or angular particles.

}, issn = {1539-3755}, doi = {10.1103/PhysRevE.87.032206}, author = {Boton, Mauricio and Emilien Az{\'e}ma and Estrada, Nicolas and Farhang Radja{\"\i} and Lizcano, Arcesio} } @proceedings {363, title = {Shear strength and microstructure of 3D assemblies of platy particles}, journal = {7th International Conference on Micromechanics of Granular Media (Powders and Grains)}, volume = {Book Series: AIP Conference Proceedings POWDERS AND GRAINS 2013}, year = {2013}, month = {Jul-09-2013}, pages = {519-522}, address = {JUL 08-12 2013 Sydney, AUSTRALIA}, abstract = {

As a first step towards particle-scale modeling of clayey soils, we investigate the mechanical behavior and microstructure of assemblies of three-dimensional rectangular platy particles by means of the discrete element method. Several samples composed of particles of different levels of platyness (ratio of width to thickness) were numerically prepared and sheared to large deformations. We analyze the shear strength, packing fraction, connectivity, contact and force anisotropies, and mobilization of friction forces as functions of platyness.We find that both the mechanical behavior and microstructure are strongly dependent on the degree of platyness. This happens, in particular, because of the alignment of particle faces along a particular direction. Additionally, as observed for other granular materials with complex shapes, the packing fraction passes by a peak value before decreasing for larger values of platyness.

}, keywords = {Clays, DEM, Microstructure, Shear strength}, url = {https://hal.archives-ouvertes.fr/hal-00842787}, author = {Boton, Mauricio and Emilien Az{\'e}ma and Estrada, Nicolas and Farhang Radja{\"\i} and Lizcano, Arcesio}, editor = {Yu, A and Dong, K and Yang, R} } @article {60, title = {Simple views on surface stress and surface energy concepts}, journal = {Advances in Natural Sciences: Nanoscience and Nanotechnology}, volume = {5}, year = {2013}, month = {Nov-27-2013}, pages = {013002}, abstract = {

Some aspects of the thermodynamics and mechanics of solid surfaces, in particular with respect to surface stress and surface energy, are reviewed. The purpose is to enlighten the deep differences between these two physical quantities. We consider successively the case of atomic flat surfaces and the case of vicinal surfaces characterized by surface stress discontinuities. Finally, experimental examples, concerning Si surfaces, are described.

}, doi = {10.1088/2043-6262/5/1/013002}, author = {M{\"u}ller, Pierre and Andres Sa{\`u}l and Leroy, Fr{\'e}d{\'e}ric} } @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} } @article {287, title = {Reply to {\textquotedblleft}Comment on {\textquoteleft}Determination of the bulk melting temperature of nickel using Monte Carlo simulations: Inaccuracy of extrapolation from cluster melting temperatures{\textquoteright} {\textquotedblright}}, journal = {Physical Review B}, volume = {84}, year = {2011}, month = {Sep-26-2011}, pages = { Article Number: 096103}, abstract = {

The thermodynamic integration (TI) method for determining the melting temperature has to be performed with great care to ensure accuracy. Since the slopes of the Gibbs free energy curves as a function of temperature (or pressure) are usually rather close, small shifts in the free energy curves can lead to a large change in the calculated melting temperature. In this Reply to a Comment by Harvey and Gheribi [ Phys. Rev. B 84 096102 (2011)] on our paper [ Phys. Rev. B 81 064112 (2010)], we show that the error in the calculated melting temperature of nickel induced by some technical imperfections in our application of the TI method, as noticed and described in the Comment, is small compared to the error bar given in our paper. We also clarify supposed doubts about the accuracy of the TI method based on Johnson\&$\#$39;s free energy expression and parametrization [ Johnson et al. Mol. Phys. 78 591 (1993)] of the Lennard-Jones fluid.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.84.096103}, author = {Los, J. H. and Roland Jean-Marc Pellenq} } @article {288, title = {Tight binding within the fourth moment approximation: Efficient implementation and application to liquid Ni droplet diffusion on graphene}, journal = {Physical Review B}, volume = {84}, year = {2011}, month = {Aug-31-2011}, pages = { Article Number: 085455}, abstract = {

Application of the fourth moment approximation (FMA) to the local density of states within a tight binding description to build a reactive, interatomic interaction potential for use in large scale molecular simulations, is a logical and significant step forward to improve the second moment approximation, standing at the basis of several, widely used (semi-)empirical interatomic interaction models. In this paper we present a sufficiently detailed description of the FMA and its technical implications, containing the essential elements for an efficient implementation in a simulation code. Using a recent, existing FMA-based model for C-Ni systems, we investigated the size dependence of the diffusion of a liquid Ni cluster on a graphene sheet and find a power law dependence of the diffusion constant on the cluster size (number of cluster atoms) with an exponent very close to \−2/3, equal to a previously found exponent for the relatively fast diffusion of solid clusters on a substrate with incommensurate lattice matching. The cluster diffusion exponent gives rise to a specific contribution to the cluster growth law, which is due to cluster coalescence. This is confirmed by a simulation for Ni cluster growth on graphene, which shows that cluster coalescence dominates the initial stage of growth, overruling Oswald ripening.

}, issn = {1098-0121}, doi = {10.1103/PhysRevB.84.085455}, author = {Los, J. H. and Christophe Bichara and Roland Jean-Marc Pellenq} }