Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomistic fracture”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Artificial neural network potentials for mechanics and fracture dynamics of two-dimensional crystals **

Understanding the mechanics and failure of materials at the nanoscale is critical for their engineering and applications. The accurate atomistic modeling of brittle failure with crack propagation in covalent crystals requires a quantum mechanics-based description of individual bond-breaking events. Artificial neural network potentials (NNPs) have emerged to overcome the traditional, physics-based modeling tradeoff between accuracy and accessible time and length scales. Previous studies have shown successful applications of NNPs for describing the structure and dynamics of molecular systems and amorphous or liquid phases of materials. However, their application to deformation and failure processes in materials is still uncommon. In this study, we discuss the apparent limitations of NNPs for the description of deformation and fracture under loadings and propose a way to generate and select training data for their employment in simulations of deformation and fracture simulations of crystals. We applied the proposed approach to 2D crystalline graphene, utilizing the density-functional tight-binding method for more efficient and extensive data generation in place of density functional theory. Then, we explored how the data selection affects the accuracy of the developed artificial NNPs. It revealed that NNP’s reliability should not only be measured based on the total energy and atomic force comparisons for reference structures but also utilize comparisons for physical properties, e.g. stress–strain curves and geometric deformation. In sharp contrast to popular reactive bond order potentials, our optimized NNP predicts straight crack propagation in graphene along both armchair and zigzag (ZZ) lattice directions, as well as higher fracture toughness of ZZ edge direction. Our study provides significant insight into crack propagation mechanisms on atomic scales and highlights strategies for NNP developments of broader materials.

2D materials↗
DOE OSTI2020

Atomistic Study of the Plastic Deformation of Transition Metals and High Entropy Alloys (Final Report)

The overarching goal of the research performed in the studies supported by the DOE Basic Energy Sciences has been investigation of mechanical properties of materials employed in technologies related to energy production and conversion using atomic level computer modeling. The materials studied are body-centered-cubic (BCC) transition metals, including α-iron, their alloys and most recently also compositionally complex, so called high entropy alloys (HEAs), which are believed to become future materials in high temperature applications. The goal of the research has been to achieve atomic level understanding of the deformation and fracture of these materials starting from the electronic structure, coarse grain the material energetics into potentials describing interaction between atoms and proceed via atomistic modeling and mesoscopic models to investigation of the glide of dislocations and thus advance the theory of plastic deformation of studied materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2022

Atomistic Materials Modeling of High-Pressure Hydrogen Interactions in Ethylene Propylene Diene Monomer (EPDM) Rubber

Elastomeric rubbers serve a vital role as sealing materials in the hydrogen storage and transport infrastructure. With applications including O-rings and hose-liners, these components are exposed to pressurized hydrogen at a range of temperatures, cycling rates, and pressure extremes. Cyclic (de)pressurization is known to degrade these materials through the process of cavitation. This readily visible failure mode occurs as a fracture or rupture of the material and is due to the oversaturated gas localizing to form gas bubbles. Computational modeling in the Hydrogen Materials Compatibility Program (H-Mat), co-led by Sandia National Laboratories and Pacific Northwest National Laboratory, employs multi-scale simulation efforts to build a predictive understanding of hydrogen-induced damage in materials. Modeling efforts within the project aim to provide insight into how to formulate materials that are less sensitive to high-pressure hydrogen-induced failure. In this document, we summarize results from atomistic molecular dynamics simulations, which make predictive assessments of the effects of compositional variations in the commonly used elastomer, ethylene propylene diene monomer (EPDM).

08 HYDROGEN↗
DOE OSTI2022

Calcium silicate hydrate colloid at different humidities: Microstructure, deformation mechanism, and mechanical properties

Calcium silicate hydrate (C-S-H) is the main hydration product of Portland cement, while fundamental questions persist regarding its colloidal structure and water content. In this work, we developed a full atomistic C-S-H colloid model and investigated the physical and mechanical properties at different gel water concentrations. Water confined in nanometer gel pores degrades the mechanical properties of C-S-H because of its hydrolysis effect. A “ductile to brittle” transition was observed in the C-S-H colloid model with increasing gel water contents, evolving from transgranular failure to intergranular fracture. The gel water confined in nanometer-sized pores displays glassy characteristics, and it acts as a lubricant and facilitates structure distortions during the deformation. We elucidated the topology effect of the hydrogen bond network on the evolution of the viscoelasticity of C-S-H at different humidities. In conclusion, this work sheds light on the role of water in C-S-H colloids from nanostructure perspectives.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Atomistic simulation of brittle-to-ductile transition in silicon carbide embedded with nano-sized helium bubbles

In this study, the tensile response of cubic silicon carbide (SiC) bulk containing cavities (voids and He bubbles) has been investigated using molecular dynamic simulations. The formation of cavities in SiC leads to a significant degradation in the mechanical properties of SiC with more influence on material fracture than initial elastic deformation. The brittle-to-ductile transition occurs in cavity-embedded SiC as the pressure in He bubbles increases. This is associated with the deformation mechanism that bond breaking at a low He bubble pressure transfers to extensive dislocation activities at a higher He bubble pressure. The cavities can effectively concentrate stress around them in the direction perpendicular to the tension, which leads to preferred cracking in the region with a higher tensile stress. The failure mechanism as revealed by this study improves understanding of property degradation in SiC that may be useful for applications of SiC in advanced nuclear energy systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Water Weakening of Calcium Oxide

Chemomechanical processes such as water weakening can control the permeability and deformation of rocks and manmade materials. Here, atomistic modeling and nanomechanical experiments were used to identify molecular origins of chemomechanical effects in calcium oxide (CaO) and its effect on observed elastic, plastic, and brittle deformation. Classical molecular dynamics simulations with the bond order-based reactive force-field ReaxFF were used to assess brittle fracture. In the presence of water, CaO fractured earlier and more often during quasi-static loading, with a calculated reduction in fracture toughness of ~80% associated with changes in the stress distribution around the crack tip. Experimentally, elastic and plastic deformation of CaO surfaces exposed to water was assessed experimentally using in situ liquid nanoindentation. Nanoindentation showed that following exposure to water, the contact hardness decreased by 1–2 orders of magnitude and decreased the modulus by 2–3 orders of magnitude due to surface hydroxylation. Furthermore, the strong chemomechanical effects on the mechanical processes in CaO suggests that minerals with similar structures may exhibit comparable effects, influencing the stability of cements and geomaterials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Evaluating the Role of Metastable Surfaces in Mechanochemical Reduction of Molybdenum Oxide

Mechanochemistry and mechanocatalysis are gaining increasing attention as environmentally friendly chemical processes because of their solvent-free nature and scalability. Significant effort has been devoted for studying continuum-scale phenomena in mechanochemistry, such as temperature and pressure gradients, but the atomic-scale mechanisms remain relatively unexplored. In this work, we focus on the mechanochemical reduction of MoO3 as a case study. We use experimental techniques to determine the mechanochemical reduction conditions and density functional theory (DFT) simulations to establish an atomistic framework for identifying the metastable surfaces that are most likely to enable this process. Our results show that metastable surfaces can significantly lower or remove thermodynamic barriers for surface reduction and that kinetic energy from milling can facilitate the formation of metastable surfaces that have high surface fracture energies and are not thermally accessible. These findings indicate that metastable surfaces are an important aspect of mechanochemistry along with hot spots and other continuum-scale phenomena.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Polymer-cement composites with adhesion and re-adhesion (healing) to casing capability for geothermal wellbore applications

Deterioration of cement/casing adhesion in wellbore scenarios can result in unwanted and potentially harmful leakage with the potential of serious repair costs. In this work, we explore the use of self-healing polymers added to conventional wellbore cements as a way to bring about self-healing and readhering (to steel casing) properties to the composite material. The polymers are pH resistant and seem to improve the cement integrity after exposure to typical chemical and thermal stresses encountered under geothermal wellbore conditions. We find that addition of about 10-15 wt% of polymer to the cement visually increases its resistance to fracturing from exposure to geothermal conditions, while the adhesive strength of cement/stainless steel increases with curing time for a period of about 10 days. Self-healing capability was demonstrated by permeability analysis showing that polymer-cement composites reduce flow by 50-70% at cement bulk and at the cement/steel interface. Use of atomistic simulations imply that these polymers have good wetting properties on the steel surfaces. Analysis of the interactions between steel/polymer and cement/polymer show that they are complementary, resulting in a wider range of bonding patterns. Cracks are likely to expose under-coordinated sites that result in more bonding interactions, which agrees well with the permeability measurements showing high degree of healed cracks and healed (cement-steel) interfacial gaps together with an overall increased in structural integrity of these advanced polymer-cement composite materials.

Rod, Kenton A.↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗
DOE OSTI2025

Ab Initio Molecular Dynamics Insights into Stress Corrosion Cracking and Dissolution of Metal Oxides

Oxide phases such as α-Fe2O3 (hematite) and α-Al2O3 (corundum) are highly insoluble in water; however, subcritical crack growth has been observed in humidity nonetheless. Chemically induced bond breaking at the crack tip appears unlikely due to sterically hindered molecular transport. The molecular mechanics of a crack in corundum with a reactive force field reveal minimal lattice trapping, leading to bond breaking before sufficient space opens for water transport. To address this, we model a pre-built blunt crack with space for H2O molecule adsorption at the tip and show that it reduces fracture toughness by lowering the critical J-integral. Then, we explore stress-enhanced dissolution to understand the mechanism of crack tip blunting in the oxide/water system. Density functional theory combined with metadynamics was employed to describe atomic dissolution from flat hematite and corundum surfaces in pure water. Strain accelerates dissolution, stabilizing intermediate states with broken bonds before full atom detachment, while the free energy profile of unstrained surfaces is almost monotonic. The atomistic calculations provided input for a kinetic model, predicting the shape evolution of a blunt crack tip, which displays three distinct regimes: (i) dissolution primarily away from the tip, (ii) enhanced blunting near but not at the apex, and (iii) sharpening near the apex. The transition between regimes occurs at a low strain, highlighting the critical role of water in the subcritical crack growth of oxide scales, with dissolution as the fundamental microscopic mechanism behind this process.

Chemistry↗
DOE OSTI2021

Atomistic Mechanics of Torn Back Folded Edges of Triangular Voids in Monolayer WS 2

Triangular nanovoids in 2D materials transition metal dichalcogenides have vertex points that cause stress concentration and lead to sharp crack propagation and failure. Here, the atomistic mechanics of back folding around triangular nanovoids in monolayer WS 2 sheets is examined. Combining atomic-resolution images from annular dark-field scanning transmission electron microscopy with reactive molecular modelling, it is revealed that the folding edge formation has statistical preferences under geometric conditions based on the orientation mismatch. It is further investigated how loading directions and strong interlayer friction, interplay with WS 2 lattice's crack preference, govern the deformation and fracture pattern around folding edges. Finally, these results provide fundamental insights into the combination of fracture and folding in flexible monolayer crystals and the resultant Moiré lattices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elastic Properties and Deformation Mechanisms in the van der Waals Single-Crystalline Indium Selenide

The bulk van der Waals (vdW) single-crystalline indium selenide (β-InSe) displays exceptional plasticity in a layered crystalline form at both micro- and macroscale. However, the nanoscale origin of plasticity remains unclear. Herein, an atomic-level study on the deformation mechanisms of InSe by using first-principles calculations is reported. Remarkable anisotropic elasticity is observed in the vdW InSe layered crystal, and the stiffness is dramatically softened because of the vdW gap in the layered structure. The simulations capture the distinct fracture modes in the uniaxial tensile deformation, depending on the loading directions—brittle fracture in the [100] and [110] directions while ductile failure in the [001] direction. The InSe layered crystal structure exhibits superplastic deformability under uniaxial compression. Different transition pathways, including interlayer tangling, amorphization, and cross-linking, are tracked along respective deformation directions. The unprecedented plasticity of InSe layered crystals can be attributed to the phase transition coupled with interlayer gliding and cross-layer dislocation slipping. Here, this study deepens our understanding of the deformation mechanisms of layered materials at the atomistic level and provides insights into tailoring material properties for low-dimensional material design based on its deformation mechanisms.

2D materials↗
DOE OSTI2025

Atomistic Modeling of the Interphase Structure between Alumina and Al12Fe2Cr Quasicrystal Approximant

Molecular dynamics simulations are performed to study the stability of an fcc metallic interlayer phase between an iron aluminide quasicrystal approximant (QCA) phase (84.2 Al, 5.7 Fe, 5.4 Ni, and 4.7 at.% Cr) and aluminum oxides (α-Al 2 O 3 and γ-Al 2 O 3 ). The presence of the interlayer phase was experimentally observed in some regions of the inner diameter surface of a TPBAR cladding tube. The simulations employed existing many-body (MEAM and EAM) interatomic potentials to describe atomic interactions in the interlayer and QCA phases, and Buckingham+Coulomb pair interactions for the alumina and between the alumina and the QCA (or fcc). Interface systems of QCA/alumina, fec/alumina, and QCA/fcc are modeled and the fracture energy of the interfaces is calculated. Based on the fracture energies, for the γ-Al 2 O 3 , both MEAM and EAM suggest that the fcc interlayer phase would not form. For the α-Al 2 O 3 , the MEAM predicts that the fcc phase would form, however, the EAM predicts that it forms only if the α-Al 2 O 3 is Al-terminated at the interface. For the O-terminated α-Al 2 O 3 , the EAM predicts the fcc phase would not form, consistent with the results for γ-Al 2 O 3 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Mapping the structural–mechanical landscape of amorphous carbon with ReaxFF molecular dynamics

We use ReaxFF molecular dynamics (MD) to investigate the relationship between structural and mechanical properties in bulk and nanostructured amorphous carbon (a-C). The liquid-quench MD method is used to generate isotropic bulk samples with mass densities ranging from 0.96 to 3.29 g/cm3. Structural analysis identifies two types of structures with distinct short- and medium-range order: lower-density sp2-dominated a-C, which is characterized by a bimodal ring-size distribution, and higher-density sp3-dominated tetrahedral amorphous carbon (ta-C), exhibiting a unimodal ring-size distribution. Stress–strain MD simulations and analysis reveal how an atomistic structure impacts elastic properties and post-yield atomic rearrangements. All stretched structures demonstrate elastic isotropy and plasticity driven by a ring-size expansion mechanism reflected in changes in ring statistics. The plastic region is substantially larger in ta-C than in a-C due to the post-yield shift from sp3 to sp2 C dominant bonding. In both a-C and ta-C, ultimate failure occurs when a reactive crack, traversed by long sp chains, forms and propagates predominantly perpendicular to the direction of the applied strain. Oxygen infiltration into the fractured region significantly reduces stress resistance, primarily through the early rupture of long sp chains. MD simulations and analysis are extended to a-C slabs, a-C nanotubes, and partially a-C nanotubes. The latter nanostructure highlights the differences between the elastically isotropic a-C walls, which develop circumferential cracking, and the crystalline walls, which tear along crystallographic directions. These results provide a strong foundation for further computational characterization of a-C materials.

Dernov, A. (ORCID:0009000004220973)↗
DOE OSTI2023

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗
DOE OSTI2020

Structural Implications of Interfacial Hydrogen Bonding in Hydrated Wyoming-Montmorillonite Clay

Montmorillonite (MMT) clay - a layered porous nanomaterial used as seals in engineered waste containment barriers for spent nuclear fuel - adopts discrete hydration/swelling states depending upon surrounding water and ion activities and confining pressure. The structure of nanoconfined water and charge-balancing counterions in the clay mineral interlayers dictate the swelling and mechanical behavior of MMT, so a molecular model for this clay with high structural fidelity is required to accurately predict the reliability of long-term nuclear waste storage. Here, we present a molecular model for MMT that is based on high resolution transmission electron microscopy of Wyoming-MMT single crystals. Imaging data unambiguously show a cis-vacant arrangement of structural hydroxyl groups in the octahedral sheet, whereas existing molecular models assume a centrosymmetric trans-vacant configuration for MMT. Furthermore, using atomistic simulations, we find that the cis-vacant arrangement of structural hydroxyl groups significantly affects the structure of adsorbed water yielding a larger population of hydrogen bonds with bridging oxygens on the tetrahedral sheet and weak hydrogen bonding between the hydroxyl groups in the octahedral sheet and water in the clay mineral interlayers. As a result, water adsorbed in the interlayer is more "ice-like", with stronger ordering and lower density, although the diffusivity of the interlayer species is not significantly diminished. Our improved structural model for MMT provides insight into the energetics of water adsorption, which ultimately dictates its pore- to macro-scale swelling, transport, and fracture properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗