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. 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.
Elevated temperature and pressure in the earth's subsurface alters the permeability of salt formations, due to changing properties of the salt-brine interface. Molecular dynamics (MD) simulations are used to investigate the mechanisms of temperature and pressure dependence of liquid-solid interfacial tensions of NaCl, KCl, and NaCl-KCl brines in contact with (100) salt surfaces. Salt-brine dihedral angles vary between 55 and 76° across the temperature (300-450 K) and pressure range (0-150 MPa) evaluated. Temperature-dependent brine composition results in elevated dihedral angles of 65-80°, which falls above the reported salt percolation threshold of 60°. Mixed NaCl-KCl brine compositions increased this effect. Elevated temperatures excluded dissolved Na+ ions from the interface, causing the strong temperature dependence of the liquid-solid interfacial tension and the resulting dihedral angle. Therefore, at higher temperature, pressure, and brine concentrations Na-Cl systems may underpredict the dihedral angle. Higher dihedral angles in more realistic mixed brine systems maintain low permeability of salt formations due to changes in the structure and energetics of the salt-brine interface.
Formation of zeolite supported Ag0 clusters depends on a combination of thermodynamically stable atomic configurations, charge balance considerations, and mobility of species on the surface and within pores. Periodic density functional theory (DFT) calculations were performed to evaluate how the location of Al in the mordenite (MOR) framework and humidity control Ag0 nanocluster formation. Four Al framework sites were studied (T1-T4) and the Al positions in the framework were identified by the shifts in the differential Al⋯Al pair distribution function (PDF). Furthermore, structural information about the Ag0 nanoclusters, such as dangling bonds, can be identified by Ag⋯Ag PDF data. For Ag0 formation in vacuum MOR structures with a Si:Al ratio of 5:1 with Al in the T1 position resulted in the most framework flexibility and the lowest Ag0 nanocluster charge, indicating the best result for formation of charge neutral nanoclusters. When water is present, Al in the T3 and T4 positions results in the formation of the smallest average Ag0 nanoclusters plus greater expansion of the O-T-O bond angle than in vacuum, indicating easier diffusion of the Ag0 nanoclusters to the surface. The presence of Al in 4-membered rings and in pairs indicates favorable MOR structures for formation of single Ag atoms, despite the existence of synthesis challenges. Therefore, Al in the T2 position is the least favorable for Ag0 nanocluster formation in both vacuum and in the presence of water. Al in the T1, T3, and T4 positions provides beneficial effects through framework flexibility and changes in nanocluster size or charge that can be leveraged for design of zeolites for formation of metallic nanoclusters.
Brittle materials, such as cement, compose major portions of built infrastructure and are vulnerable to degradation and fracture from chemo-mechanical effects. Currently, methods of modeling infrastructure do not account for the presence of a reactive environment, such as water, on the acceleration of failure. Here, we have developed methodologies and models of concrete and cement fracture that account for varying material properties, such as strength, shrinkage, and fracture toughness due to degradation or hydration. The models have been incorporated into peridynamics, non-local continuum mechanics methodology, that can model intersecting and branching brittle fracture that occurs in multicomponent brittle materials, such as concrete. Through development of new peridynamic capabilities, decalcification of cement and differential shrinkage in clay-cement composites have been evaluated, along with exemplar problems in nuclear waste cannisters and wellbores. We have developed methods to simulate multiphase phenomena in cement and cement-composite materials for energy and infrastructure applications.
Modeling the degradation of cement-based infrastructure due to aqueous environmental conditions continues to be a challenge. In order to develop a capability to predict concrete infrastructure failure due to chemical degradation, we created a chemomechanical model of the effects of long-term water exposure on cement paste. The model couples the mechanical static equilibrium balance with reactive–diffusive transport and incorporates fracture and failure via peridynamics (a meshless simulation method). The model includes fundamental aspects of degradation of ordinary Portland cement (OPC) paste, including the observed softening, reduced toughness, and shrinkage of the cement paste, and increased reactivity and transport with water induced degradation. This version of the model focuses on the first stage of cement paste decalcification, the dissolution of portlandite. Given unknowns in the cement paste degradation process and the cost of uncertainty quantification (UQ), we adopt a minimally complex model in two dimensions (2D) in order to perform sensitivity analysis and UQ. We calibrate the model to existing experimental data using simulations of common tests such as flexure, compression and diffusion. Then we calculate the global sensitivity and uncertainty of predicted failure times based on variation of eleven unique and fundamental material properties. We observed particularly strong sensitivities to the diffusion coefficient, the reaction rate, and the shrinkage with degradation. Also, the predicted time of first fracture is highly correlated with the time to total failure in compression, which implies fracture can indicate impending degradation induced failure; however, the distributions of the two events overlap so the lead time may be minimal. Extension of the model to include the multiple reactions that describe complete degradation, viscous relaxation, post-peak load mechanisms, and to three dimensions to explore the interactions of complex fracture patterns evoked by more realistic geometry is straightforward and ongoing.
