@article {620, title = {Earth concrete. Stabilization revisited}, journal = {Cement and Concrete Research}, volume = {114}, year = {2018}, month = {Dec-2018}, pages = {90 - 102}, abstract = {

Not surprisingly, with the increased awareness of environmental issues, construction with raw (crude, unbaked) earth (subsoil) is gaining renewed interest. However, it suffers from a poor image and from the difficulty to meet modern productivity standards and to pass some durability tests designed for industrial materials. The recent trend is to overcome these drawbacks by \“stabilizing\” the material most often with Portland cement (PC). Here we show that stabilization with PC is in general neither technically nor environmentally advisable. It brings only moderate mechanical improvement at a high environmental cost. Rather than massively transforming crude earth into a low quality concrete, it would be more appropriate to adapt the architectural practice and/or to look for milder ways to improve properties. In this respect, the recent successful attempts to improve the workability and the strength of raw earth by controlling the dispersion of its fine fraction seem to be particularly promising.

}, issn = {00088846}, doi = {10.1016/j.cemconres.2017.02.035}, url = {https://linkinghub.elsevier.com/retrieve/pii/S0008884616308365}, author = {Henri Van Damme and Houben, Hugo} } @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 {632, title = {On the existence and origin of sluggish diffusion in chemically disordered concentrated alloys}, journal = {Current Opinion in Solid State and Materials Science}, volume = {22}, year = {2018}, month = {Jun-2018}, pages = {65-74}, abstract = {

Concentrated single phase solid solutions, including medium- and high-entropy alloys, represent a new class of materials that have recently attracted significant interest due to exceptional functional and structural properties. Their fascinating properties are mainly attributed to the sluggish atomic-level diffusion and transport, but its controlling mechanisms are largely unknown and there is certain skepticism about its very existence. By using microsecond-scale molecular dynamics, on-the-fly and conventional kinetic Monte Carlo, we reveal the governing role of percolation effects and composition dependence of the vacancy migration energy in diffusion. Surprisingly, an increase of concentration of faster species (Fe) in face-centered cubic Ni-Fe alloy may decrease the overall atomic diffusion. Consequently, the composition dependence of tracer diffusion coefficient has a minimum near the site percolation threshold, similar to 20 at.\%Fe. We argue that this coupled percolation and composition-dependent barriers for vacancy jumps within different subsystems in medium- and high-entropy alloys leads, indeed, to the sluggish diffusion. A fast method for preselecting materials with potentially desired properties is suggested.

}, issn = {13590286}, doi = {10.1016/j.cossms.2018.05.003}, url = {https://linkinghub.elsevier.com/retrieve/pii/S135902861830038X}, author = {Osetsky, N. and B{\'e}land, Laurent Karim and Barashev, Alexander V. and Zhang, Yanwen} } @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 {321, title = {Early-Age Stress and Pressure Developments in a Wellbore Cement Liner: Application to Eccentric Geometries}, journal = {Journal of Applied Mechanics}, volume = {83}, year = {2016}, month = {Sep-2016}, pages = {Article Number: 091012}, abstract = {

This paper introduces a predictive model for the stress and pressure evolutions in a wellbore cement liner at early ages. A pressure state equation is derived that observes the coupling of the elastic changes of the solid matrix, the eigenstress developments in the solid and porespaces, and the mass consumption of water in course of the reaction. Here, the transient constitution of the solid volume necessitates advancing the mechanical state of the poroelastic cement skeleton incrementally and at constant hydration degree. Next, analytic function theory is employed to assess the localization of stresses along the steel\–cement (SC) and rock\–cement (RC) interfaces by placing the casing eccentrically with respect to the wellbore hole. Though the energy release rate due to complete debonding of either interface is only marginally influenced by the eccentricity, the risk of evolving a microcrack along the thick portion of the sheath is substantially increased. Additionally, it is observed that the risk of microannulus formation is principally affected by the pressure rebound, which is engendered by the slowing reaction rate and amplified for rock boundaries with low permeability.

}, issn = {0021-8936}, doi = {10.1115/1.4034013}, author = {Thomas Alexander Petersen and Franz-Josef Ulm} } @article {250, title = {Effect of Polydispersity of Clay Platelets on the Aggregation And Mechanical Properties of Clay at the Mesoscale}, journal = {Clays and Clay Minerals}, volume = {64}, year = {2016}, month = {Aug-2016}, pages = {425 - 437}, abstract = {

The results from mesoscale simulations of the formation and evolution of microstructure for assemblies of Na-smectite particles based on assumed size distributions of individual clay platelets are presented here. The analyses predicted particle arrangements and aggregation (i.e. platelets linked in face-face configurations) and are used to link geometric properties of the microstructure and mechanical properties of the particle assemblies. Interactions between individual ellipsoidal clay platelets are represented using the Gay-Berne potential based on atomistic simulations of the free energy between two Na-smectite clay-platelets in liquid water, following a novel coarse-graining method developed previously. The current study describes the geometric (aggregate thickness, orientation, and porosity) and elastic properties in the \‘jammed states\’ from the mesoscale simulations for selected ranges of clay particle sizes and confining pressures. The thickness of clay aggregates for monodisperse assemblies increases (with average stack thickness consisting of n = 3-8 platelets) with the diameter of the individual clay platelets and with the level of confining pressure. Aggregates break down at high confining pressures (50-300 atm) due to slippage between the platelets. Polydisperse simulations generate smaller aggregates (n = 2) and show much smaller effects of confining pressure. All assemblies show increased order with confining pressure, implying more anisotropic microstructure. The mesoscale simulations are also in good agreement with macroscopic compression behavior measured in conventional 1-D laboratory compression tests. The mesoscale assemblies exhibit cubic symmetry in elastic properties. The results for larger platelets (D = 1000 \Å) are in good agreement with Nano-indentation measurements on natural clays and shale samples.

}, issn = {0009-8604}, doi = {10.1346/CCMN.2016.0640407}, author = {Ebrahimi, Davoud and Andrew J. Whittle and Roland Jean-Marc Pellenq} } @article {217, title = {Effect of Water on Elastic and Creep Properties of Self-Standing Clay Films}, journal = {Langmuir}, volume = {32}, year = {2016}, month = {Feb-09-2016}, pages = {1370 - 1379}, abstract = {

We characterized experimentally the elastic and creep properties of thin self-standing clay films, and how their mechanical properties evolved with relative humidity and water content. The films were made of clay montmorillonite SWy-2, obtained by evaporation of a clay suspension. Three types of films were manufactured, which differed by their interlayer cation: sodium, calcium, or a mixture of sodium with calcium. The orientational order of the films was characterized by X-ray diffractometry. The films were mechanically solicited in tension, the resulting strains being measured by digital image correlation. We measured the Young\&$\#$39;s modulus and the creep over a variety of relative humidities, on a full cycle of adsorption-desorption for what concerns the Young\&$\#$39;s modulus. Increasing relative humidity made the films less stiff and made them creep more. Both the elastic and creep properties depended significantly on the interlayer cation. For the Young\&$\#$39;s modulus, this dependence must originate from a scale greater than the scale of the clay layer. Also, hysteresis disappeared when plotting the Young\&$\#$39;s modulus versus water content instead of relative humidity. Independent of interlayer cation and of relative humidity greater than 60\%, after a transient period, the creep of the films was always a logarithmic function of time. The experimental data gathered on these mesoscale systems can be of value for modelers who aim at predicting the mechanical behavior of clay-based materials (e.g., shales) at the engineering macroscopic scale from the one at the atomistic scale, for them to validate the first steps of their upscaling scheme. They provide also valuable reference data for bioinspired clay-based hybrid materials.

}, issn = {0743-7463}, doi = {10.1021/acs.langmuir.5b03431}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Bornert, Michel and Ferrage, Eric and Hubert, Fabien and Henri Van Damme} } @article {618, title = {Electromagnetic Casimir energy of a disk opposite a plane}, journal = {Physical Review A}, volume = {94}, year = {2016}, month = {Sep-15-2016}, abstract = {

Building on work by J. Meixner [Z. Naturforschung 3a, 506 (1948)], we show how to compute the exact scattering amplitude (or T-matrix) for electromagnetic scattering from a perfectly conducting disk. This calculation is a rare example of a nondiagonal T-matrix that can nonetheless be obtained in a semianalytic form. We then use this result to compute the electromagnetic Casimir interaction energy for a disk opposite a plane, for arbitrary orientation angle of the disk, for separations greater than the disk radius. We find that the proximity force approximation (PFA) significantly overestimates the Casimir energy, in the case of both the ordinary PFA, which applies when the disk is parallel to the plane, and the \“edge PFA\”, which applies when the disk is perpendicular to the plane.

}, issn = {2469-9926}, doi = {10.1103/PhysRevA.94.032509}, url = {https://link.aps.org/doi/10.1103/PhysRevA.94.032509}, author = {Emig, Thorsten and Graham, Noah} } @article {269, title = {An equation of state for granular media at the limit state of isotropic compression}, journal = {EPL (Europhysics Letters)}, volume = {114}, year = {2016}, month = {Apr-2016}, pages = {Article Number: 14004}, abstract = {

It is well believed that the volumetric entropy of Edwards captures part of the physics of granular media, but it is still unclear whether it can be applied to granular systems under mechanical stress. By working out a recent proposal by Aste, Di Matteo et al. to measure Edwards\&$\#$39; compactivity from the volume distribution of Vorono\ï or Delaunay tessellations (Phys. Rev. E, 77 (2008) 021309), and assuming that the total volume divides into elementary cells of fixed minimal volume, we derive an equation of state relating the compactivity to the packing fraction, and we show by extensive molecular-dynamics simulations that this equation and its underlying assumption describe well the volumetric aspects of both the limit state of isotropic compression and the limit state of shear (also called critical state in soil mechanics) for three-dimensional ensembles of mono-disperse spheres, for a broad range of the sliding and rolling friction coefficients. In addition, by using the limit state of isotropic compression as testing ground, we find that the compactivity, the entropy per elementary cell and the number of elementary cells per grain computed by this method are the same within statistical precision, either by using Vorono\ï, Delaunay, or centroidal Vorono\ï tessellations, allowing thus for an objective definition. This means that not only Aste\&$\#$39;s cell method is robust and suitable to measure Edwards\&$\#$39; compactivity of granular systems under mechanical stress but also the actual nature of the elementary cells might be unimportant.

}, issn = {0295-5075}, doi = {10.1209/0295-5075/114/14004}, author = {Oquendo, W. F. and Mu{\~n}oz, J. D. and Farhang Radja{\"\i}} } @article {247, title = {Evolution of organo-clay composites with respect to thermal maturity in type II organic-rich source rocks}, journal = {Geochimica et Cosmochimica Acta}, volume = {195}, year = {2016}, month = {Dec-15-2016}, pages = {68 - 83}, abstract = {

Among the proposed pathways of sequestration and progressive transformation of organic matter (OM) during burial, the sorptive protection mechanism explains the strong relationship between total organic carbon (TOC) and mineral surface area (MSA) noted in numerous black shales around the globe. The complete mechanistic framework of preservation and evolution of OM in organo-mineral complexes remains, however, an enigma. On one hand, smectite layers composing the clay fraction are known to have a major influence on available surface area. OM maturation occurs, on the other hand, concurrently with the smectite illitization that provokes the closure of the interlayer spaces. The potential of smectite layers in the sequestration and preservation of organic molecules and the fate of these molecules with respect to the smectite illitization were therefore addressed. Here, the mineralogy of three organic-rich source rocks of various maturities was characterized in regards with the geochemistry of their OM. A thorough examination of the clay minerals present in the clay matrices provided evidences of mixed layer minerals containing smectite and illite layers with an increasing illite component with respect to maturity. The comprehensive interpretation of the X-ray diffractograms and analytical electron microscopy results suggested the presence of organic molecules in the inter-particulate and possibly the interlayer spaces of the smectite-rich components in immature source rocks. This eventuality was further supported by the presence of intercalated clay-organic nanocomposites observed by transmitted electron microscopy coupled with energy dispersive spectroscopy. Textural observations also showed that the increased illite content found in the overmature sample led to the reorganization of the OM and the clay particles into nano-scale aggregates. These results clarify the geochemical mechanism beyond the reported relationship between TOC and MSA and allow generalizing it to various organic-rich source rocks. Eventually, the understanding of the close organo-clay association raises prospects on the prediction of the evolution of the mechanical properties of economically valuable source rocks with respect to maturity.

}, issn = {00167037}, doi = {10.1016/j.gca.2016.09.008}, author = {Berthonneau, Jeremie and Grauby, Olivier and Muhannad Abuhaikal and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Henri Van Damme} } @article {135, title = {Effect of Chain Length and Pore Accessibility on Alkane Adsorption in Kerogen}, journal = {Energy \& Fuels}, volume = {29}, year = {2015}, month = {Oct-28-2015}, pages = {7889 - 7896}, abstract = {
Abstract Image

Configurational-biased Grand Canonical Monte Carlo and Molecular Dynamics simulations were used to investigate the adsorption and structure of n-alkanes in realistic models of kerogen (the organic phase in gas shales). Both the effects of the n-alkane length (from methane to dodecane) and of the adsorbent porosity/density and connectivity were considered. For all n-alkanes, due to the subnano pore size of kerogen, adsorption follows a Langmuir-type adsorption isotherm as the adsorbed amount increases continuously with pressure until maximum loading is reached. While all n-alkanes get adsorbed very efficiently in low density kerogen, size exclusion in high density kerogen lead to low adsorbed amounts for molecules larger than hexane. Using Molecular Dynamics simulations reflecting the setup of actual adsorption experiments, we also probed the effect of pore accessibility on n-alkane adsorption in kerogen. All pores are connected to the surface and, hence, accessible for the low density kerogen. In contrast, most pores are isolated and inaccessible in the kerogen with the high density. We also address the effect of mesoporosity and its connection to nanoporosity on alkane adsorption in kerogen.

}, issn = {0887-0624}, doi = {10.1021/acs.energyfuels.5b02015}, author = {Kerstyn I. Falk and Roland Jean-Marc Pellenq and Franz-Josef Ulm and Benoit A. Coasne} } @article {161, title = {Effect of curvature and confinement on the Casimir-Polder interaction}, journal = {Physical Review A}, volume = {91}, year = {2015}, month = {Jan-30-2015}, abstract = {

Modifications of Casimir-Polder interactions due to confinement inside a cylindrical cavity and due to curvature in- and outside the cavity are studied. We consider a perfectly conducting cylindrical shell with a single particle (atom or macroscopic sphere) located next to its interior or exterior surface, or two atoms placed inside the shell. By employing the scattering approach, we obtain the particle-cavity interaction and the modification of the two-particle interaction due to the cavity. We consider both retardation and thermal effects. While for the atoms a dipole description is sufficient, for the macroscopic sphere we sum (numerically) over many multipole fluctuations to compute the interaction at short separations. In the latter limit we compare to the proximity approximation and a gradient expansion and find agreement. Our results indicate an confinement induced suppression of the force between atoms. General criteria for suppression and enhancement of Casimir interactions due to confinement are discussed.

}, issn = {1050-2947}, doi = {10.1103/PhysRevA.91.012516}, author = {Rodriguez-Lopez, Pablo and Emig, Thorsten and Noruzifar, Ehsan and Zandi, Roya} } @article {182, title = {Effect of Surface Texture on Freezing in Nanopores: Surface-Induced versus Homogeneous Crystallization}, journal = {Langmuir}, volume = {31}, year = {2015}, month = {Mar-10-2015}, pages = {2706 - 2713}, abstract = {
Abstract Image

Freezing of argon in ordered and disordered carbon pores of a similar diameter D \∼ 2.4 nm is investigated using extensive molecular simulations with large system sizes up to 104 atoms. While crystallization in the atomistically smooth pore consists in a surface-induced phase transition occurring at a temperature larger than the bulk, crystallization in the disordered pores, which is only partial as it is spatially restricted to the pore center, occurs through homogeneous crystallization. These results shed light on solidification in pores by showing that there is a crossover between surface-induced and homogeneous crystallization upon increasing the surface disorder of the host material. In the latter case, the Gibbs\–Thomson equation, in which crystallization is assumed to occur when the crystal size equals the pore size corrected for the thickness of the unfreezable layer at the pore surface, is in reasonable agreement with the observed freezing temperature.

}, issn = {0743-7463}, doi = {10.1021/la5041213}, author = {Benoit A. Coasne} } @article {129, title = {Effects of shape and size polydispersity on strength properties of granular materials.}, journal = {Phys Rev E Stat Nonlin Soft Matter Phys}, volume = {91}, year = {2015}, month = {Mar-18-2015}, pages = {Article Number: 032203}, abstract = {

By means of extensive contact dynamics simulations, we analyze the combined effects of polydispersity both in particle size and in particle shape, defined as the degree of shape irregularity, on the shear strength and microstructure of sheared granular materials composed of pentagonal particles. We find that the shear strength is independent of the size span, but unexpectedly, it declines with increasing shape polydispersity. At the same time, the solid fraction is an increasing function of both the size span and the shape polydispersity. Hence, the densest and loosest packings have the same shear strength. At the scale of the particles and their contacts, we analyze the connectivity of particles, force transmission, and friction mobilization as well as their anisotropies. We show that stronger forces are carried by larger particles and propped by an increasing number of small particles. The independence of shear strength with regard to size span is shown to be a consequence of contact network self-organization, with the falloff of contact anisotropy compensated by increasing force anisotropy.

}, issn = {1550-2376}, doi = {10.1103/PhysRevE.91.032203}, author = {Duc-Hanh Nguyen and Emilien Az{\'e}ma and Philippe Sornay and Farhang Radja{\"\i}} } @proceedings {343, title = {Evolution of particle size distributions in crushable granular materials}, journal = {3rd International Symposium on Geomechanics from Micro to Macro}, volume = {Geomechanics from Micro to Macro}, year = {2015}, month = {Feb-02-2015}, pages = {275 - 280}, publisher = {CRC Press}, address = {SEP 01-03-2014 Univ Cambridge, Cambridge, ENGLAND}, abstract = {

By means of the contact dynamics method together with a particle fracture model, in which the particles are cohesive aggregates of irreducible polygonal fragments, we investigate the evolution of particle size distribution in the process of uniaxial compaction of granular materials. The case of single particle breakup under compressive stress is used to test the method and the influence of discretization (number of irreducible fragments). We show that the breaking threshold of the granular assembly scales with the internal cohesion of the particles but it depends also on the initial size distribution and irregularity of polygonal particle shapes. The evolution of size distribution proceeds by consecutive periods of intense particle crushing, characterized by local shattering instability, and periods of little breaking activity. Starting with either monodisperse or power-law distribution of particle sizes, the latter evolves towards a broad distribution of the fragmented particles with a nearly power-law distribution in the range of intermediate particle sizes. Interestingly, a finite number of large particles survive despite ongoing crushing process due to the more homogeneous distribution of forces in the presence of small fragmented particles filling the pores between larger particles.

}, isbn = {978-1-138-02707-7}, doi = {10.1201/b1739510.1201/b17395-48}, author = {Duc-Hanh Nguyen and Emilien Az{\'e}ma and Farhang Radja{\"\i} and Philippe Sornay}, editor = {Kenichi Soga and Krishna Kumar and Giovanna Biscontin and Kuo, Matthew} } @article {148, title = {Experimental chemo-mechanics of early-age fracture properties of cement paste}, journal = {Cement and Concrete Research}, volume = {75}, year = {2015}, month = {Sep-2015}, pages = {42 - 52}, abstract = {

The risk of early-age fracture of cementitious materials in ever more challenging environments provides a unique opportunity to employ an experimental chemo-mechanical platform to develop functional relations between hydration degree, fracture and strength properties, assessed by isothermal calorimetry, micro-scratching, splitting and microindentation on white cement paste at various curing ages from 7\ h to 28\ days. We show that the modulus, tensile strength, fracture toughness and energy all evolve with a natural logarithmic dependence on the hydration degree. These trends are linked to the densification of the material during the hydration process, explained by compaction mechanics and free volume theory. We show that while the fracture process zone size is essentially constant during the hydration process, the ductility of the material, quantified by M/H, decreases, and is consistent with the evolution of Kc/H. Both quantities provide a convenient way to experimentally assess the fracture sensitivity of early-age cement-based materials.

}, issn = {00088846}, doi = {10.1016/j.cemconres.2015.04.004}, author = {Christian G. Hoover and Franz-Josef Ulm} } @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 {273, title = {Effect of size polydispersity versus particle shape in dense granular media}, journal = {Physical Review E}, volume = {90}, year = {2014}, month = {Jul-21-2014}, pages = {Article Number: 012202}, abstract = {

We present a detailed analysis of the morphology of granular systems composed of frictionless pentagonal particles by varying systematically both the size span and particle shape irregularity, which represent two polydispersity parameters of the system. The microstructure is characterized in terms of various statistical descriptors such as global and local packing fractions, radial distribution functions, coordination number, and fraction of floating particles. We find that the packing fraction increases with the two parameters of polydispersity, but the effect of shape polydispersity for all the investigated structural properties is significant only at low size polydispersity where the positional and/or orientational ordering of the particles prevail. We focus in more detail on the class of side/side contacts, which is the interesting feature of our system as compared to a packing of disks. We show that the proportion of such contacts has weak dependence on the polydispersity parameters. The side- side contacts do not percolate but they define clusters of increasing size as a function of size polydispersity and decreasing size as a function of shape polydispersity. The clusters have anisotropic shapes but with a decreasing aspect ratio as polydispersity increases. This feature is argued to be a consequence of strong force chains (forces above the mean), which are mainly captured by side-side contacts. Finally, the force transmission is intrinsically multiscale, with a mean force increasing linearly with particle size.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, issn = {1539-3755}, doi = {10.1103/PhysRevE.90.012202}, author = {Duc-Hanh Nguyen and {\'e}ma, Emilien and Farhang Radja{\"\i} and Philippe Sornay} } @article {101, title = {Effect of temperature on adsorption of mixtures in porous materials}, journal = {Molecular Simulation}, volume = {40}, year = {2014}, month = {Jan-02-2014}, pages = {45-51}, type = {Article}, abstract = {

The effect of temperature on the adsorption of a simple mixture (Ar/Kr) in disordered porous materials is investigated by means of molecular simulation. In the larger mesopores of porous silica glasses, capillary condensation occurs upon decreasing the temperature. At temperatures above the capillary condensation temperature, Kr is preferentially adsorbed at the pore surface and Ar adsorption occurs in regions of low Kr density. For temperatures below the capillary condensation temperature, Ar density surprisingly increases as temperature increases, the behaviour that is consistent with an over-solubility effect. In contrast, in the disordered sub-nanoporous carbon, filling of the pores occurs in a reversible and continuous way upon decreasing the temperature, owing to the small size and amorphous shape of the pores. These results show that the crossover between capillary condensation and continuous reversible filling observed for pure fluids in pores also exists for mixtures. We also show that the Kr selectivity exhibits a minimum in the disordered porous silica that is located at the capillary condensation temperature. In contrast, in the disordered porous carbon where no capillary condensation occurs, the selectivity decreases monotonically with increasing the temperature. These results shed light on low-temperature adsorption of mixtures confined in porous materials and provide a guide to design efficient phase separation processes.

}, keywords = {Adsorption, grand canonical Monte Carlo simulations, Mixtures, porous materials}, issn = {0892-7022}, doi = {10.1080/08927022.2013.852192}, author = {Kerstyn I. Falk and Benoit A. Coasne and Roland Jean-Marc Pellenq} } @article {114, title = {Elastic Properties of Swelling Clay Particles at Finite Temperature upon Hydration}, journal = {Journal of Physical Chemistry C}, volume = {118}, year = {2014}, month = {May-01-2014}, pages = {8933-8943}, type = {Article}, abstract = {
Abstract Image

The knowledge of mechanical properties of swelling clays on the scale of the clay layer is of crucial importance to build predictive hydromechanical macroscopic constitutive laws in a bottom-up approach. In this work, we computed the elastic properties of particles of a swelling clay (i.e., montmorillonite) by molecular dynamics simulations. Because of the softness of the material in the direction orthogonal to the clay layer with increasing water content, the computation of the whole stiffness tensor was not trivial: we needed to implement the elastic bath method,1 which allows us to attenuate thermal strains and compute the stiffness via strain fluctuations in isothermal\–isobaric simulations. We investigated the effect of water content, temperature, and the interlayer cation on the elastic properties of swelling clays. In particular, we showed that the out-of-plane stiffness coefficients computed at 300 K differed significantly from the same coefficients computed at 0 K. The out-of-plane coefficients were very sensitive to both temperature and water content. In contrast, the dependence of the in-plane coefficients to temperature was slight. Moreover, the decrease of the in-plane coefficients with water content could be entirely explained by the change of geometry of the system due to the swelling of the interlayer space. The interlayer cation impacted the elastic properties of montmorillonite only in the driest states.

}, issn = {1932-7447}, doi = {10.1021/jp412160e}, author = {Benoit Carrier and Matthieu Vandamme and Roland Jean-Marc Pellenq and Henri Van Damme} } @article {85, title = {Enhanced H-2 Uptake of n-Alkanes Confined in Mesoporous Materials}, journal = {Journal of Physical Chemistry C }, volume = {118}, year = {2014}, month = {Apr-30-2014}, pages = {10720-10727}, type = {Article}, abstract = {
Abstract Image

The hydrogen uptake in hybrid sorbents consisting of n-alkane solvents confined in mesoporous silica aerogel is measured at different temperatures from 273 to 313 K and pressures up to 40 bar. An apparent \“oversolubility\” effect is observed as the H2 uptake in the hybrid sorbents is much larger than that in bulk solvents. The H2 uptake in the hybrid sorbents is found to increase with increasing temperature, which suggests that the flexibility and conformation of n-alkane molecules confined in the aerogel play a crucial role; high-entropy (disordered) alkane configurations lead to the creation of numerous cavities which make it possible to solubilize a larger number of H2 molecules. This departs from adsorption-driven solubility effects for which the number of solubilized molecules decreases with increasing temperature. For a given temperature and pressure, it is found that the number of solubilized H2 molecules per unit volume increases with decreasing alkane chain length. Such an effect, which is observed for both the bulk alkanes and the alkanes confined in the silica aerogel, can be rationalized by considering the number density of CHx (x = 2 or 3) groups; for a given temperature, the latter number density decreases with decreasing alkane chain length so that the free volume available to solubilize H2 molecules increases.

}, issn = {1932-7447}, doi = {10.1021/jp411526f}, author = {St{\'e}phanie Clauzier and Linh Ngoc Ho and Pera-Titus, M. and David Farrusseng and Benoit A. Coasne} } @article {328, title = {Erratum to {\textquotedblleft}Static states of cohesive granular media{\textquotedblright}}, journal = {Journal of Mechanical Science and Technology}, volume = {28}, year = {2014}, month = {Oct-2014}, pages = {4345 - 4345}, abstract = {

vol 21, pg 1957, 2007

This paper investigates the influence of cohesion on the behavior of 2D granular media simulated by using a molecular dynamics method, involving a simple contact law with adhesion. The study considers an adhesion index which is non intrinsic but significant for interpreting the mechanical effects of applied pressure. A static state system of 4000 particles is numerically simulated by oedometrical compression under a constant force, without gravity. The results show how the geometrical texture and the network of contact forces change according to the level of adhesion. This effect is explained essentially in term of the increase in the number of the tensile contacts and by the development of the internal self-stress structure. As in the case of non cohesive granular media, a high spatial heterogeneity of the contact forces is observed.

}, issn = {1738-494X}, doi = {10.1007/s12206-014-0951-3}, author = {Preechawuttipong, I. and Peyroux, R. and Farhang Radja{\"\i} and Rangsri, W.} } @article {218, title = {Environmental scanning electron microscopy (ESEM) and nanoindentation investigation of the crack tip process zone in marble}, journal = {Acta Geotechnica}, volume = {8}, year = {2013}, month = {Jun-2013}, pages = {223 - 245}, abstract = {

This study explores the interaction between crack initiation and nanomechanical properties in the crack-tip fracture process zone of Carrara marble. Specimens with preexisting cracks were loaded in a uniaxial testing machine until the process zone appeared at the tips of the preexisting cracks. ESEM analysis reveals an increase in microcrack density in the process zone with increased loading of the specimen. Nanoindentation testing comprised of lines and grids of single nanoindentations located both near and far from the process zone shows a decrease in both indentation modulus and indentation hardness near grain boundaries in intact material, and with closeness to the process zone. Ultimately, the study confirms that the crack-tip process zone manifests itself as an area of reduced indentation hardness and indentation modulus in marble.

}, issn = {1861-1125}, doi = {10.1007/s11440-013-0213-z}, author = {Brooks, Z. and Franz-Josef Ulm and Einstein, H. H.} } @article {46, title = {ESEM Study of the Humidity-Induced Swelling of Clay Film}, journal = {Langmuir}, volume = {29}, year = {2013}, month = {Oct-02-2013}, pages = {12823{\textendash}12833}, abstract = {
Abstract Image

We measured the humidity-induced swelling of thin self-standing films of montmorillonite clay by a combination of environmental scanning electron microscopy (ESEM) and digital image correlation (DIC). The films were about 40 \μm thick. They were prepared by depositing and evaporating a suspension of clay and peeling off the highly oriented deposits. The rationale for creating such original samples was to obtain mesoscopic samples that could be used to bridge experimentally the gap between the scale of the clay layer and the engineering scale of a macroscopic clay sample. Several montmorillonite samples were used: the reference clay Swy-2, the same clay homoionized with sodium or calcium ions, and a sodium-exchanged Cloisite. The edges of the clay films were observed by ESEM at various relative humidity values between 14\% and 95\%. The ESEM images were then analyzed by DIC to measure the swelling or the shrinkage of the films. We also measured the adsorption/desorption isotherms by weighing the film samples in a humidity-controlled environment. In order to analyze our results, we compared our swelling/shrinkage and adsorption/desorption data with previously published data on the interlayer spacing obtained by X-ray diffraction and with numerical estimates of the interlayer water obtained by molecular dynamics simulation. The swelling and the hysteresis of this swelling were found to be comparable for the overall macroscopic films and for the interlayer space. The same correspondence between film and interlayer space was observed for the amount of adsorbed water. This suggests that, in the range of relative humidities values explored, the films behave like freely swelling oriented stacks of clay layers, without any significant contribution from the mesoporosity. The relevance of this result for the behavior of clayey sedimentary rocks and the differences with the behavior of nonoriented samples (powders or compacted powders) are briefly discussed.

}, doi = {10.1021/la402781p}, author = {Benoit Carrier and Wang, Linlin and Matthieu Vandamme and Roland Jean-Marc Pellenq and Bornert, Michel and Tanguy, Alexandre and Henri Van Damme} } @proceedings {350, title = {Effect of Particle Shape non-Convexity on the Rheology of Granular Media : 3D Contact Dynamics Simulations}, journal = {2nd International Conference on Particle-Based Methods - Fundamentals and Applications (Particles)}, volume = {PARTICLE-BASED METHODS II: FUNDAMENTALS AND APPLICATIONS}, year = {2012}, month = {Apr-10-2012}, pages = {427-434}, address = {OCT 26-28 2011 Barcelona, SPAIN}, abstract = {

We analyze the effect of particle shape non-convexity on the quasi-static behavior of granular materials by means of contact dynamics simulations. The particles are regular aggregates of four overlapping spheres described by a nonconvexity parameter depending on the relative positions of the particles. Several packings are first submitted to isotropic compression without friction. We find that, as in 2D, the solid fraction of isotropic packings increases with non-convexity up to a maximum value and then declines to be nearly equal to that of a packing composed of only spheres. It is also remarkable that the coordination number increases quickly and saturates so that the packings composed of grains with a high level of nonconvexity are looser but more strongly connected. Then, the quasi-static behavior, structural and force anisotropies are analyzed by subjecting each packing to a triaxial compression. We find that the shear strength increases with non-convexity. We show that this increase results from the presence of multiple contacts between trimers leading to enhanced frictional interlocking.

}, keywords = {force transmission, Granular Materials, non-convexity, particle shape, texture}, url = {https://hal.archives-ouvertes.fr/hal-00686453}, author = {Baptiste Saint-Cyr and Emilien Az{\'e}ma and Jean-Yves Delenne and Farhang Radja{\"\i} and Philippe Sornay}, editor = {Onate, E and Owen, DRJ} } @article {52, title = {Enhanced CO2 solubility in hybrid adsorbents: optimization of solid support and solvent properties for CO2 capture}, journal = {The Journal of Physical Chemistry C}, volume = {116}, year = {2012}, month = {Feb-9-2012}, pages = {3600{\textendash}3607}, abstract = {
Abstract Image

In this study, we optimize the CO2 adsorption performance of hybrid adsorbents prepared by confining physical solvents in porous solid supports. A number of prospective solid supports and physical solvents are chosen to prepare hybrid adsorbents, and are subsequently evaluated in CO2 adsorption experiments. Generally, all the hybrid adsorbents show an enhancement of CO2 solubility compared to the bulk physical solvent. However, not all the adsorbents positively display an improvement in the CO2 adsorption performance as compared with the original solids after confining the physical solvent into the solids\’ pore. The micropore blocking effect is observed in the impregnated forms of zeolite, activated carbon, silicagel, and cecagel. Furthermore, we have obtained certain requisites for a good solid support, as efficient structures should be mesoporous with large surface area. In addition, there is an optimized solvent\’s size to achieve an optimized enhanced solubility. As a result, among the candidates, N-methyl-2-pyrrolidone confined in MCM-41 and alumina are identified as the most suitable hybrid adsorbents for an effective CO2-removal application.

}, doi = {10.1021/jp2099625}, author = {Linh Ngoc Ho and Perez-Pellitero, Javier and Porcheron, Fabien and Roland Jean-Marc Pellenq} } @article {54, title = {Evidence on the Dual Nature of Aluminum in the Calcium-Silicate-Hydrates Based on Atomistic Simulations}, journal = {Journal of the American Ceramic Society}, volume = {95}, year = {2012}, month = {Jan-23-2012}, pages = {1128{\textendash}1137}, abstract = {

Hydration of tri-calcium silicate (C3S) and di-calcium silicate (C2S) precipitates calcium-silicate-hydrate (CSH) which is the bonding phase responsible for the strength of cementitious materials. Substitution of part of C3S and C2S with aluminum-containing additives alters the chemical composition of hydration products by precipitating calcium-aluminate-silicate-hydrate (CASH). Incorporation of aluminum in the molecular building blocks of CSH entails structural and chemo-mechanical consequences. These alterations can be measured through solid state nuclear magnetic resonance (NMR) experiments. By conducting a wide spectrum of atomistic simulation methods on thousands of aluminum-containing molecular CASH structures, an overall molecular approach for determination of CASH nanostructure is presented. Through detailed analysis of different order parameters, it is found that aluminum can exhibit a tetra-/penta-/octahedral behavior which is fully consistent with the recent NMR observations. This corresponds to the formation of a class of complex three-dimensional alumino-silicate skeletons with partial healing effect in the CASH nanostructure potentially increasing durability and strength of hydration products. We explored the variation of mechanical observables by increasing aluminum content in CASH structures of varying calcium to silicon ratio. Finally, deformation of CSHs and CASHs of different chemical formula in a multi-scale fashion unravels the effect of chemical composition on the strength and kinematics of deformation in this particular type of composites.

}, doi = {10.1111/j.1551-2916.2011.05058.x}, author = {Mohammad Javad Abdolhosseini Qomi and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {259, title = {Experimental determination of the fracture toughness via microscratch tests: Application to polymers, ceramics, and metals}, journal = {Journal of Materials Research}, volume = {27}, year = {2012}, month = {Jan-2012}, pages = {485 - 493}, abstract = {

This article presents a novel microscratch technique for the determination of the fracture toughness of materials from scratch data. While acoustic emission and optical imaging devices provide quantitative evidence of fracture processes during scratch tests, the technique proposed here provides a quantitative means to assess the fracture toughness from the recorded forces and depth of penetration. We apply the proposed method to a large range of materials, from soft (polymers) to hard (metal), spanning fracture toughness values over more than two orders of magnitude. The fracture toughness values so obtained are in excellent agreement with toughness values obtained for the same materials by conventional fracture tests. The fact that the proposed microscratch technique is highly reproducible, almost nondestructive, and requires only small material volumes makes this technique a powerful tool for the assessment of fracture properties for microscale materials science and engineering applications.

}, issn = {0884-2914}, doi = {10.1557/jmr.2011.402}, author = {Ange-Therese Akono and Randall, Nicholas X. and Franz-Josef Ulm} } @article {290, title = {Enhanced CO2 Solubility in Hybrid MCM-41: Molecular Simulations and Experiments}, journal = {Langmuir}, volume = {27}, year = {2011}, month = {Jul-05-2011}, pages = {8187 - 8197}, abstract = {
Abstract Image

Grand canonical Monte Carlo simulations are performed in a hybrid adsorbent model in order to interpret the CO2 solubility behavior. The hybrid adsorbent is prepared by confining a physical solvent (OMCTS) into the pores of a mimetic MCM-41 solid support. As a result, simulated adsorption isotherms of CO2 nicely match the experimental data for three distinctive systems: bulk solvent, raw MCM-41, and hybrid MCM-41. The microscopic mechanisms underlying the apparition of enhanced solubility are then clearly identified. In fact, the presence of solvent molecules favors the layering of CO2 molecules within the pores; therefore, the CO2 solubility in the hybrid adsorbent markedly increases in comparison to that found in the raw adsorbent as well as in the bulk solvent. In addition, a good understanding of confined solvents\’ properties and solid surface structures is essential to fully evaluate the efficiency of hybrid adsorbents in capturing CO2. The sorbent\–solid interactions along with the solvent molecular size\’s impact on CO2 solubility are therefore investigated in this study. We found that an ideal hybrid system should possess a weak solvent\–solid interaction but a strong solvent\–CO2 interaction. Besides, an optimal solvent size is obtained for the enhanced CO2 solubility in the hybrid system. According to the simulation results, the solvent layer builds pseudomicropores inside the mesoporous MCM-41, enabling more CO2 molecules to be absorbed under the greater influence of spatial confinement and surface interaction. In addition, the molecular sieving effect is clearly observed in the case of larger solvent molecular sizes.

}, issn = {0743-7463}, doi = {10.1021/la2012765}, author = {Linh Ngoc Ho and Perez Pellitero, Javier and Porcheron, Fabien and Roland Jean-Marc Pellenq} } @article {293, title = {Empirical force fields for complex hydrated calcio-silicate layered materials}, journal = {Phys. Chem. Chem. Phys.}, volume = {13}, year = {2010}, month = {Nov-10-2010}, pages = {1002 - 1011}, abstract = {

The use of empirical force fields is now a standard approach in predicting the properties of hydrated oxides which are omnipresent in both natural and engineering applications. Transferability of force fields to analogous hydrated oxides without rigorous investigations may result in misleading property predictions. Herein, we focus on two common empirical force fields, the simple point charge ClayFF potential and the core-shell potential to study tobermorite minerals, the most prominent family of Calcium-Silicate-Hydrates that are complex hydrated oxides. We benchmark the predictive capabilities of these force fields against first principles results. While the structural information seem to be in close agreement with DFT results, we find that for higher order properties such as elastic constants, the core-shell potential quantitatively improves upon the simple point charge model, and shows a larger degree of transferability to complex materials. In return, to remedy the deficiencies of the simple point charge potential for hydrated calcio-silicates, we suggest using both structural data and elasticity data for potential calibration, a new force field potential, CSH-FF. This re-parameterized version of ClayFF is then applied to simulating an atomistic model of cement (Pellenq et al., PNAS, 2009). We demonstrate that this force field improves the predictive capabilities of ClayFF, being considerably less computational intensive than the core-shell model.

}, issn = {1463-9076}, doi = {10.1039/C0CP00516A}, author = {Rouzbeh Shahsavari and Roland Jean-Marc Pellenq and Franz-Josef Ulm} }