@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 {573, title = {Thermodynamics, kinetics, and mechanics of cesium sorption in cement paste: A multiscale assessment}, journal = {Physical Review Materials}, volume = {2}, year = {2018}, month = {May-30-2018}, chapter = {Article Number: 053608 }, abstract = {

Cesium-137 is a common radioactive byproduct found in nuclear spent fuel. Given its 30 year half life, its interactions with potential storage materials-such as cement paste-is of crucial importance. In this paper, simulations are used to establish the interaction of calcium silicate hydrates (C-S-H)-the mam binding phase of cement paste-with Cs at the nano- and mesoscale. Different C-S-H compositions are explored, including a range of Ca/Si ratios from 1.0 to 2.0. These calculations are based on a set of 150 atomistic models, which qualitatively and quantitatively reproduce a number of experimentally measured features of C-S-H-within limits intrinsic to the approximations imposed by classical molecular dynamics and the steps followed when building the models. A procedure where hydrated Ca2+ ions are swapped for Cs1+ ions shows that Cs adsorption in the C-S-H interlayer is preferred to Cs adsorption at the nanopore surface when Cs concentrations are lower than 0.19 Mol/kg. Interlayer sorption decreases as the Ca/Si ratio increases. The activation relaxation technique nouveau is used to access timescales out of the reach of traditional molecular dynamics (MD). It indicates that characteristic diffusion time for Cs1+ in the C-S-H interlayer is on the order of a few hours. Cs uptake in the interlayer has little impact on the elastic response of C-S-H. It leads to swelling of the C-S-H grams, but mesoscale calculations that access length scales out of the range of MD indicate that this leads to practically negligible expansive pressures for Cs concentrations relevant to nuclear waste repositories.

Figure

Figure

Figure

Figure

Figure

Figure

Figure

}, keywords = {CALCIUM-SILICATE-HYDRATE; C-S-H; MOLECULAR-DYNAMICS; CRYSTAL-STRUCTURE; MONTE-CARLO; IMMOBILIZATION; RADIONUCLIDES; SIMULATIONS; RELAXATION; COLLOIDS}, doi = {10.1103/PhysRevMaterials.2.053608}, url = {https://link.aps.org/doi/10.1103/PhysRevMaterials.2.053608}, author = {Arayro, Jack and Dufresne, Alice and Zhou, Tingtao and Katerina Ioannidou and Franz-Josef Ulm and Roland Jean-Marc Pellenq and B{\'e}land, Laurent Karim} } @article {629, title = {Time-Resolved Mechanical Spectroscopy of Soft Materials via Optimally Windowed Chirps}, journal = {Physical Review X}, volume = {8}, year = {2018}, month = {Dec-06-2018}, pages = {Article Number 041042 }, abstract = {

The ability to measure the bulk dynamic behavior of soft materials with combined time and frequency resolution is instrumental for improving our fundamental understanding of connections between the microstructural dynamics and the macroscopic mechanical response. Current state-of-the-art techniques are often limited by a compromise between resolution in the time and frequency domains, mainly due to the use of elementary input signals that have not been designed for fast time-evolving systems such as materials undergoing gelation, curing, or self-healing. In this work, we develop an optimized and robust excitation signal for time-resolved mechanical spectroscopy through the introduction of joint frequency- and amplitude-modulated exponential chirps. Inspired by the biosonar signals of bats and dolphins, we optimize the signal profile to maximize the signal-to-noise ratio while minimizing spectral leakage with a carefully designed modulation of the envelope of the chirp, obtained using a cosine-tapered window function. A combined experimental and numerical investigation reveals that there exists an optimal range of window profiles (around 10\% of the total signal length) that minimizes the error with respect to standard single-frequency sweep techniques. The minimum error is set by the noise floor of the instrument, suggesting that the accuracy of an optimally windowed-chirp (OWCh) sequence is directly comparable to that achievable with a standard frequency sweep, while the acquisition time can be reduced by up to 2 orders of magnitude, for comparable spectral content. Finally, we demonstrate the ability of this optimized signal to provide time- and frequency-resolved rheometric data by studying the fast gelation process of an acid-induced protein gel using repeated OWCh pulse sequences. The use of optimally windowed chirps enables a robust time-resolved rheological characterization of a wide range of soft materials undergoing rapid mutation and has the potential to become an invaluable rheometric tool for researchers across different disciplines.

}, doi = {10.1103/PhysRevX.8.041042}, url = {https://link.aps.org/doi/10.1103/PhysRevX.8.041042}, author = {Geri, Michela and Keshavarz, Bavand and Divoux, Thibaut and Clasen, Christian and Curtis, Daniel J. and McKinley, Gareth H.} } @article {633, title = {Two-dimensional numerical simulation of chimney fluidization in a granular medium using a combination of discrete element and lattice Boltzmann methods}, journal = {Physical Review E}, volume = {97}, year = {2018}, month = {May-10-2018}, pages = { Article Number 052902 }, abstract = {

We present here a numerical study dedicated to the fluidization of a submerged granular medium induced by a localized fluid injection. To this end, a two-dimensional (2D) model is used, coupling the lattice Boltzmann method (LBM) with the discrete element method (DEM) for a relevant description of fluid-grains interaction. An extensive investigation has been carried out to analyze the respective influences of the different parameters of our configuration, both geometrical (bed height, grain diameter, injection width) and physical (fluid viscosity, buoyancy). Compared to previous experimental works, the same qualitative features are recovered as regards the general phenomenology including transitory phase, stationary states, and hysteretic behavior. We also present quantitative findings about transient fluidization, for which several dimensionless quantities and scaling laws are proposed, and about the influence of the injection width, from localized to homogeneous fluidization. Finally, the impact of the present 2D geometry is discussed, by comparison to the real three-dimensional (3D) experiments, as well as the crucial role of the prevailing hydrodynamic regime within the expanding cavity, quantified through a cavity Reynolds number, that can presumably explain some substantial differences observed regarding upward expansion process of the fluidized zone when the fluid viscosity is changed.

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.97.052902}, url = {https://link.aps.org/doi/10.1103/PhysRevE.97.052902}, author = {Jeff Ngoma and Pierre Philippe and St{\'e}phane Bonelli and Farhang Radja{\"\i} and Jean-Yves Delenne} } @article {612, title = {Temperature distribution and heat radiation of patterned surfaces at short wavelengths}, journal = {Physical Review E}, volume = {95}, year = {2017}, month = {May-03-2017}, pages = { Article Number 052104 }, abstract = {

We analyze the equilibrium spatial distribution of surface temperatures of patterned surfaces. The surface is exposed to a constant external heat flux and has a fixed internal temperature that is coupled to the outside heat fluxes by finite heat conductivity across the surface. It is assumed that the temperatures are sufficiently high so that the thermal wavelength (a few microns at room temperature) is short compared to all geometric length scales of the surface patterns. Hence the radiosity method can be employed. A recursive multiple scattering method is developed that enables rapid convergence to equilibrium temperatures. While the temperature distributions show distinct dependence on the detailed surface shapes (cuboids and cylinder are studied), we demonstrate robust universal relations between the mean and the standard deviation of the temperature distributions and quantities that characterize overall geometric features of the surface shape.

}, issn = {2470-0045}, doi = {10.1103/PhysRevE.95.052104}, url = {http://link.aps.org/doi/10.1103/PhysRevE.95.052104}, author = {Emig, Thorsten} } @article {596, title = {Three-dimensional bonded-cell model for grain fragmentation}, journal = {Computational Particle Mechanics}, volume = {4}, year = {2017}, month = {Oct-2017}, pages = {441 - 450}, chapter = { Special Issue: SI }, abstract = {

We present a three-dimensional numerical method for the simulation of particle crushing in 3D. This model is capable of producing irregular angular fragments upon particle fragmentation while conserving the total volume. The particle is modeled as a cluster of rigid polyhedral cells generated by a Voronoi tessellation. The cells are bonded along their faces by a cohesive Tresca law with independent tensile and shear strengths and simulated by the contact dynamics method. Using this model, we analyze the mechanical response of a single particle subjected to diametral compression for varying number of cells, their degree of disorder, and intercell tensile and shear strength. In particular, we identify the functional dependence of particle strength on the intercell strengths. We find that two different regimes can be distinguished depending on whether intercell shear strength is below or above its tensile strength. In both regimes, we observe a power-law dependence of particle strength on both intercell strengths but with different exponents. The strong effect of intercell shear strength on the particle strength reflects an interlocking effect between cells. In fact, even at low tensile strength, the particle global strength can still considerably increase with intercell shear strength. We finally show that the Weibull statistics describes well the particle strength variability.

}, keywords = {Bonded-cell model; Fragmentation; Discrete element method; Contact dynamics method; Voronoi cell; Weibull statistics}, issn = {2196-4378}, doi = {10.1007/s40571-016-0129-0}, url = {http://link.springer.com/10.1007/s40571-016-0129-0}, author = {Cantor, D and Emilien Az{\'e}ma and Philippe Sornay and Farhang Radja{\"\i}} } @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 {243, title = {Topological Control on the Structural Relaxation of Atomic Networks under Stress}, journal = {Physical Review Letters}, volume = {119}, year = {2017}, month = {Jul-21-2017}, pages = {Article Number: 035502}, abstract = {

Upon loading, atomic networks can feature delayed irreversible relaxation. However, the effect of composition and structure on relaxation remains poorly understood. Herein, relying on accelerated molecular dynamics simulations and topological constraint theory, we investigate the relationship between atomic topology and stress-induced structural relaxation, by taking the example of creep deformations in calcium silicate hydrates (C-S-H), the binding phase of concrete. Under constant shear stress, C-S-H is found to feature delayed logarithmic shear deformations. We demonstrate that the propensity for relaxation is minimum for isostatic atomic networks, which are characterized by the simultaneous absence of floppy internal modes of relaxation and eigenstress. This suggests that topological nanoengineering could lead to the discovery of nonaging materials.

Figure

Figure

Figure

}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.035502}, author = {Mathieu Bauchy and Wang, Mengyi and Yu, Yingtian and Wang, Bu and Krishnan, M. Anoop and Enrico Masoero and Franz-Josef Ulm and Roland Jean-Marc Pellenq} } @article {152, title = {Transient dynamics of a 2D granular pile}, journal = {The European Physical Journal E}, volume = {38}, year = {2015}, month = {May-26-2015}, pages = {Article Number: 47}, abstract = {

We investigate by means of Contact Dynamics simulations the transient dynamics of a 2D granular pile set into motion by applying shear velocity during a short time interval to all particles. The spreading dynamics is directly controlled by the input energy whereas in recent studies of column collapse the dynamics scales with the initial potential energy of the column. As in column collapse, we observe a power-law dependence of the runout distance with respect to the input energy with nontrivial exponents. This suggests that the power-law behavior is a generic feature of granular dynamics, and the values of the exponents reflect the distribution of kinetic energy inside the material. We observe two regimes with different values of the exponents: the low-energy regime reflects the destabilization of the pile by the impact with a runout time independent of the input energy whereas the high-energy regime is governed by the input energy. We show that the evolution of the pile in the high-energy regime can be described by a characteristic decay time and the available energy after the pile is destabilized.

https://static-content.springer.com/image/art\%3A10.1140\%2Fepje\%2Fi2015-15047-x/MediaObjects/10189_2015_133_Figa_HTML.gif
}, issn = {1292-8941}, doi = {10.1140/epje/i2015-15047-x}, author = {Patrick Mutabaruka and Krishna Kumar and Kenichi Soga and Farhang Radja{\"\i} and Jean-Yves Delenne} } @proceedings {359, title = {Transient regime to fluidized chimney within a granular bed by means of a 2D DEM/LBM modeling}, journal = {IV International Conference on Particle-based Methods (PARTICLES 2015)}, volume = {PARTICLE-BASED METHODS IV-FUNDAMENTALS AND APPLICATIONS}, year = {2015}, month = {sept-30-2015}, pages = {558-566}, publisher = {International Center for Numerical Methods in Engineering (CIMNE)}, address = {SEP 28-30 2015 Barcelona, SPAIN}, abstract = {

Beyond a given threshold, an upward fluid flow at constant flowrate, injected through a small size section, is able to generate a fluidization along a vertical chimney over the entire height of a granular assembly. Fluidization is first initiated in the immediate vicinity of the injection hole and then the fluidized zone grows gradually until reaching the upper surface of the granular packing. In this work, we present numerical results on the kinetics of chimney fluidization in an immersed granular bed produced with two-dimensional simulations coupling the Discrete Element and Lattice Boltzmann Methods (DEM-LBM). A parametric study is carried out with 11 different sets of physical parameters and analyzed based on spatio-temporal diagrams. Then a dimensional analysis allows finding general scaling laws for both threshold and growth rate of the fluidized zone by use of two dimensionless numbers, namely Reynolds and Archimedes numbers, while quite simple empirical relationships can also be proposed.

}, keywords = {DEM, Fluidization, Granular Materials, LBM}, url = {https://hal.archives-ouvertes.fr/hal-01269325}, author = {Jeff Ngoma and Pierre Philippe and St{\'e}phane Bonelli and Cuellar, Pablo and Jean-Yves Delenne and Farhang Radja{\"\i}}, editor = {Onate, E and Bischoff, M and Owen, DRJ} } @article {295, title = {Tensile strength and fracture of cemented granular aggregates}, journal = {The European Physical Journal E}, volume = {35}, year = {2012}, month = {Nov-2012}, pages = { Article Number: 117}, abstract = {

Cemented granular aggregates include a broad class of geomaterials such as sedimentary rocks and some biomaterials such as the wheat endosperm. We present a 3D lattice element method for the simulation of such materials, modeled as a jammed assembly of particles bound together by a matrix partially filling the interstitial space. From extensive simulation data, we analyze the mechanical properties of aggregates subjected to tensile loading as a function of matrix volume fraction and particle-matrix adhesion. We observe a linear elastic behavior followed by a brutal failure along a fracture surface. The effective stiffness before failure increases almost linearly with the matrix volume fraction. We show that the tensile strength of the aggregates increases with both the increasing tensile strength at the particle-matrix interface and decreasing stress concentration as a function of matrix volume fraction. The proportion of broken bonds in the particle phase reveals a range of values of the particle-matrix adhesion and matrix volume fraction for which the cracks bypass the particles and hence no particle damage occurs. This limit is shown to depend on the relative toughness of the particle-matrix interface with respect to the particles.

}, issn = {1292-8941}, doi = {10.1140/epje/i2012-12117-7}, author = {Rafik Aff{\`e}s and Jean-Yves Delenne and Yann Monerie and Farhang Radja{\"\i} and Vincent Topin} } @article {93, title = {Thermodynamics of Water Confined in Porous Calcium-Silicate-Hydrates}, journal = {Langmuir}, volume = {28}, year = {2012}, month = {Aug-07-2012}, pages = {11422-11432}, type = {Article}, abstract = {
Abstract Image

Water within pores of cementitious materials plays a crucial role in the damage processes of cement pastes, particularly in the binding material comprising calcium-silicate-hydrates (C\–S\–H). Here, we employed Grand Canonical Monte Carlo simulations to investigate the properties of water confined at ambient temperature within and between C\–S\–H nanoparticles or \“grains\” as a function of the relative humidity (\%RH). We address the effect of water on the cohesion of cement pastes by computing fluid internal pressures within and between grains as a function of \%RH and intergranular separation distance, from 1 to 10 \Å. We found that, within a C\–S\–H grain and between C\–S\–H grains, pores are completely filled with water for \%RH larger than 20\%. While the cohesion of the cement paste is mainly driven by the calcium ions in the C\–S\–H, water facilitates a disjoining behavior inside a C\–S\–H grain. Between C\–S\–H grains, confined water diminishes or enhances the cohesion of the material depending on the intergranular distance. At very low \%RH, the loss of water increases the cohesion within a C\–S\–H grain and reduces the cohesion between C\–S\–H grains. These findings provide insights into the behavior of C\–S\–H in dry or high-temperature environments, with a loss of cohesion between C\–S\–H grains due to the loss of water content. Such quantification provides the necessary baseline to understand cement paste damaging upon extreme thermal, mechanical, and salt-rich environments.

}, issn = {0743-7463}, doi = {10.1021/la301738p}, author = {Patrick A Bonnaud and Qing Ji and Benoit A. Coasne and Roland Jean-Marc Pellenq and Krystyn J. Van Vliet} } @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} }