SEARCH · Engineering Papers

Results for “elastic modulus”

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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Probabilistic Evaluation of Geomechanical Risks in CO2 Storage: An Exploration of Caprock Integrity Metrics Using a Multilaminate Model

The probabilistic uncertainty assessment of geomechanical risk—specifically, caprock failure—attributable to CO2 injection, as presented in a simplified hypothetical geological model, was the focus of this study. Our approach amalgamates the implementation of a multilaminate model, the creation of a response surface model in conjunction with the Box–Behnken sampling design, the execution of associated numerical modeling experiments, and the utilization of Monte Carlo simulations. Probability distributions to encapsulate the inherent variability (elastic and mechanical properties of the caprock and reservoir) and uncertainty in prediction estimates (vertical displacement, total strain, and F value) were employed. Our findings reveal that the Young modulus of the caprock is a key factor controlling equivalent total strain but is insufficient as a stand-alone indicator of caprock integrity. It is confirmed that the caprock can accommodate significant deformation without failure, if it possesses a low Young’s modulus and high mechanical strength properties, such as the friction angle and uniaxial compressive strength. Similarly, vertical displacement was found to be an unreliable indicator for caprock integrity, as caprock failure can occur across a broad spectrum of vertical displacements, particularly when both the Young modulus and mechanical strength properties have wide ranges. This study introduces the F value as the most dependable indicator for caprock failure, although it is a theoretical attribute (the shortest distance between the Mohr circle and the nearest failure envelope used to measure the sensitivity to failure) and not physically measurable in the field. Deviatoric stress levels were found to vary based on stress regimes, with the maximum levels observed under extensive and compressive stress regimes. In conjunction with the use of the response surface method, this study demonstrates the efficacy of the multilaminate framework and the Mohr–Coulomb constitutive model in providing a simplified, yet effective, probabilistic model of the mechanical behavior of caprock failure, reducing mathematical and computational complexities.

Energy & Fuels↗

DOE OSTI2024

Surfactant templated biogenic nanoporous silica thermal insulation composite

Carbon-sequestration biogenic insulation materials offer a promising solution for sustainable and energy-efficient buildings. However, minimizing energy-intensive manufacturing processes and ensuring a uniform microstructure are essential to achieve widespread application. In this work, we introduce cetyltrimethylammonium bromide (CTAB) surfactant as both synthesis and dispersion template for nanoporous silica to create biogenic straw-based thermal insulation composites. As the amount of nanoporous silica increases to 72 wt%, the elastic insulation composite dried using the solvent exchange method exhibits a density of 0.058 g cm –3 , a thermal conductivity of 30.3 mW m –1 K –1 and a flexural modulus of 3.96 MPa, while demonstrating materials circularity and fire retardance. The use of CTAB surfactant effectively captures nanoporous silica, preventing its loss during manufacturing and enhancing the homogeneity and dispersion of the silica within the composite. In addition, the results show that the ethanol solvent exchange drying at ambient conditions provides the optimum thermal insulation performance with less energy consumption compared to freeze-drying.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Data-Driven Tailoring Optimization of Thermoset Polymers Using Ultrasonics and Machine Learning

Thermoset polymers are highly demanded for their structural robustness, thermal stability, and chemical resistance. Tailoring the properties of these polymers for high-performance applications is often preferred to designing brand-new polymers. However, the traditional destructive techniques used to characterize their properties as a function of manufacturing parameters are expensive and time-consuming. A novel non-destructive, data-driven method leveraging ultrasonics and machine learning techniques to tailor the properties of thermosets as a function of the manufacturing parameters is demonstrated. Thermoset epoxy samples with varying curing temperatures (15–40 °C) and curing agent amounts (±40%) were manufactured and tested. Their curing kinetics were monitored by determining the sound speed in the material in real time, while the longitudinal modulus of the samples was determined post-cure. Machine learning models were developed using a k-nearest neighbors algorithm. These models were implemented to predict the curing and final elastic properties using the manufacturing parameters, i.e., stoichiometry and curing temperature, and vice versa. Understanding and modeling how these parameters affect the cure kinetics and final properties will allow for efficient and reliable optimization of thermoset tailoring and manufacturing.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Selective Plasticization of Poly (ethylene oxide) (PEO) Block in Nanostructured Polystyrene– PEO– Polystyrene Triblock Copolymer Electrolytes

The plasticization of a polymer electrolyte usually promotes its ionic conductivity but decreases its storage modulus due to the increased polymer chain flexibility. Herein, we show that such a tradeoff between the ionic conductivity and the mechanical robustness of the polymer electrolyte can be alleviated by selective plasticization of the ion-conductive block, such as poly(ethylene oxide) (PEO) in a polystyrene (PS)– PEO–PS block copolymer (SEO) electrolyte using an ether type plasticizer, tetraethylene glycol dimethyl ether (TEGDME). In this work, at maximum plasticizer loading, the room temperature ionic conductivity increases by up to 3 orders, whereas the storage modulus, G' reduces to half, is still on the order of 10 2 MPa. At above the melting temperature of the PEO block, the dynamic storage modulus, G' of the plasticized membrane surpasses its dry PS-PEO-PS counterpart. Such a phenomenon results from that, a) TEGDME co-crystallizes with PEO to promote its crystallinity and hence the storage modulus, b) TEGDME swells the amorphous PEO phase to enhance the polymer chain segmental mobility and hence ionic conductivity, and c) the PS phase remains intact from TEGDME to keep the SEO elastic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Determination of Elastic Coefficients for a Composite with Crossing Steel - 20086

For a composite material with an embedding of crossing reinforcement bars, the elastic coefficients are calculated by solving a boundary-value problem (elastostatic problem). The background material is concrete whereas the bars are made of steel. The bars cross each other horizontally at a right angle. The boundary-value problem is solved numerically by using 3-D quadratic basis functions in ABAQUS{sup R}. Under external shear loading conditions, the effective elastic coefficients on the macroscale are quite larger for a cross-section in which the bars are immediately below the ones in the cross-section normal to it. The shear modulus in the off-plane direction is about one and half times larger than that in the in-plane direction. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2026

Effects of PLA/PHB Blend Ratio and Wood Flour Loading on the Melt Rheology and Thermomechanical Properties of Biobased Polymer Composites

The development of sustainable, biobased polymer composites is crucial for reducing supply chain dependence on fossil fuels. A major challenge lies in understanding and predicting the processability of these sustainable material alternatives, which directly impacts large-scale manufacturing. Here, this study systematically investigates the effect of poly(lactic acid) (PLA)/polyhydroxybutyrate (PHB) blend ratios and wood flour (WF) loadings on the thermal, rheological, and mechanical properties of the composites. Composites were prepared via melt compounding and characterized by using differential scanning calorimetry (DSC), thermogravimetric analysis (TGA), oscillatory shear melt rheology, tensile testing, and scanning electron microscopy. Results showed that increasing the PHB content lowered the glass transition temperature, while higher PHB and WF loadings decreased thermal stability. Complex viscosity decreased with increasing PHB in 0 and 10% WF, whereas it increased in 20% and 30% WF for PHB-rich blends (80/20 and 70/30 PLA/PHB) due to stronger filler interactions with flexible PHB chains. At 30% WF, PHB-containing composites also showed elastic dominance with G′ higher than G″. This behavior was confirmed by van Gurp–Palmen plots, where the phase angle decreased with increasing PHB and WF loadings. WF addition also shifted Cole–Cole plots away from semicircular terminal relaxation, while time–temperature superposition validity was maintained, confirming that polymer chain dynamics continue to dominate melt behavior. Young’s modulus increased with WF loading, while tensile strength and toughness decreased due to weaker interfacial adhesion between the matrix and the filler. Interestingly, 10% WF systems containing PHB exhibited the highest toughness among all of the filled systems, indicating synergistic reinforcement from WF fillers and PHB-induced ductility. Overall, melt rheology effectively captured the internal structure and stiffness of the composites as the formulation changed, reinforcing the concept that the PLA/PHB ratio and WF content can be tuned to adjust melt elasticity and balance mechanical property trade-offs for targeted applications.

melt rheology↗

DOE OSTI2022

Vickers hardness prediction from machine learning methods

Abstract The search for new superhard materials is of great interest for extreme industrial applications. However, the theoretical prediction of hardness is still a challenge for the scientific community, given the difficulty of modeling plastic behavior of solids. Different hardness models have been proposed over the years. Still, they are either too complicated to use, inaccurate when extrapolating to a wide variety of solids or require coding knowledge. In this investigation, we built a successful machine learning model that implements Gradient Boosting Regressor (GBR) to predict hardness and uses the mechanical properties of a solid (bulk modulus, shear modulus, Young’s modulus, and Poisson’s ratio) as input variables. The model was trained with an experimental Vickers hardness database of 143 materials, assuring various kinds of compounds. The input properties were calculated from the theoretical elastic tensor. The Materials Project’s database was explored to search for new superhard materials, and our results are in good agreement with the experimental data available. Other alternative models to compute hardness from mechanical properties are also discussed in this work. Our results are available in a free-access easy to use online application to be further used in future studies of new materials at www.hardnesscalculator.com .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Behavioral simulations and experimental evaluations of stress induced spatial nonuniformity of dynamic bulk modulus in additive friction stir deposited AA 6061

Additive friction stir deposition (AFSD) is a rapidly developing solid-state additive manufacturing (AM) technology that enables the fabrication of customized products on a relatively large scale. The unique forge-liked process offers high quality products with refined equiaxed microstructure and exceptional mechanical performance. However, in AFSD, the spatially distributed mechanical properties within a deposited layer have not been established so far, due to the lack of a suitable large-scale mapping method. In the present work, an ultrasonic dynamic-elasticity mapping method was employed to map the distribution of physical properties in a AFSD 6061 Al. Strong spatial asymmetries in elastographies were mainly due to the residual stresses. The origin of the asymmetries was investigated through thermo-mechanical simulations of the AFSD process. It was observed that majority of the asymmetries were associated with the mechanical residual stresses instead of the symmetric thermal residual stresses during linear deposition of the AFSD process. However, the thermal residual stresses can be manipulated to become asymmetric with a 90° turn in the linear deposition path, due to the difference between the levels of cooling along the edges of the inner and outer corners. Hence, the mechanical and thermo-mechanical residual stresses can be additive or subtractive depending on the direction of the turns. Here, such numerically predicted behaviors were also experimentally observed on the AFSD fabricated samples by ultrasonic elastography. Furthermore, additional preliminary studies on multi-layer AFSD fabrications show the asymmetric mechanical stresses generated during the deposition process on different layers generate constructive or destructive interference which results to spatially variable dynamic bulk modulus depending on the deposition orientations of layers.

36 MATERIALS SCIENCE↗

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 OSTI2021

X-ray diffraction and in situ pressurization of dentine apatite reveals nanocrystal modulus stiffening upon carbonate removal

Bone-like materials comprise carbonated-hydroxyapatite nanocrystals (c-Ap) embedding a fibrillar collagen matrix. The mineral particles stiffen the nanocomposite by tight attachment to the protein fibrils creating a high strength and toughness material. The nanometer dimensions of c-Ap crystals make it very challenging to measure their mechanical properties. Mineral in bony tissues such as dentine contains 2~6 wt.% carbonate with possibly different elastic properties as compared with crystalline hydroxyapatite. Here we determine strain in biogenic apatite nanocrystals by directly measuring atomic deformation in pig dentine before and after removing carbonate. Transmission electron microscopy revealed the platy 3D morphology while atom probe tomography revealed carbon inside the calcium rich domains. High-energy X-ray diffraction in combination with in situ hydrostatic pressurization quantified reversible c-Ap deformations. Additionally, crystal strains differed between annealed and ashed (decarbonated) samples, following 1 or 10 h heating at 250 °C or 550 °C respectively. Measured bulk moduli (K) and a-/c-lattice deformation ratios (η) were used to generate synthetic K syn and ηsyn identifying the most likely elastic constants C 33 and C 13 for c-Ap. These were then used to calculate the nanoparticle elastic moduli. For ashed samples, we find an average E 11 =107 GPa and E 33 =128 GPa corresponding to ~5% and ~17% stiffening of the a-/c-axes of the nanocrystals as compared with the biogenic nanocrystals in annealed samples. Ashed samples exhibit ~10% lower Poisson's ratios as compared with the 0.25~0.36 range of carbonated apatite. Carbonate in c-Ap may therefore serve for tuning local deformability within bony tissues.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗

DOE OSTI2024

Digital Multiphase Composites via Additive Manufacturing

Abstract Mechanical properties of traditional engineering materials are typically coupled to each other, presenting a challenge to practitioners with multi‐dimensional material property requirements. In this work, continuous, independent control over multiple mechanical properties is demonstrated in composite materials realized using additive manufacturing. For the first time, composites additively manufactured from rigid plastic, soft elastomer, and liquid constituents are experimentally characterized, demonstrating materials which span four orders of magnitude in modulus and two orders of magnitude in toughness. By forming analytical mappings between relative concentrations of constituents at the microscale and resulting macroscale material properties, inverse material design is enabled; the method is showcased by printing artifacts with prescribed toughness and elasticity distributions. The properties of these composites are placed in the context of biological tissues, showing they have promise as mechanically plausible tissue mimics.

Smith, Lawrence T.↗

DOE OSTI2021

A new coupled multiphase flow–finite strain deformation–fault slip framework for induced seismicity

Production of hydrocarbons and water from subsurface reservoirs are known to cause permanent deformation of the reservoir and seismicity along faults both of which are detrimental to sustainable development of natural resources. Most of the prior studies on understanding fluid flow-induced plasticity and seismicity have focused on one or the other phenomenon due to the numerical difficulty associated with simultaneous modeling of the two failure phenomena because stress and deformation evolve non-linearly in both plasticity and seismicity. However, in reservoirs undergoing long-term production, plastic failure can alter the stress paths of points on a fault such that the onset, location, and magnitude of actual seismic events can no longer be predicted by a poroelastic simulation due to inaccurate stress and deformation history. We present a computational framework for coupled multiphase flow, finite strain poroplastic deformation, and dynamic fault slip and use it to understand the impact of plastic deformation on the onset, location, and magnitude of induced fault slip events. We evaluate the impact of plasticity on reservoir pressure, deformation, induced stress, and fault slip by comparing infinitesimal strain elastic and finite strain poro-elastoplastic models. For real-world applications, we consider different scenarios where the reservoir is either mechanically weaker or stronger than the caprock. We analyze the stress evolution as a function of the change in reservoir pressure to understand the role of contrast in reservoir and caprock elastic moduli on geomechanical stability. The results show that the poroplastic reservoir exhibits larger vertical deformation and delayed slip than the poroelastic reservoir after the same amount of oil production. For the same amount of pressure drop, a reservoir with a smaller modulus than the caprock displays a larger vertical displacement and an earlier onset of both plastic failure and fault slip. For a reservoir with a larger modulus than the caprock, vertical displacement is larger on the reservoir top boundary and smaller on the ground surface, and a higher pressure drop is needed to induce plastic failure and fault slip.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Resonant ultrasound spectroscopy measurement and modeling of additively manufactured octet truss lattice cubes

In this work, the elastic characteristics of additively manufactured Ti-5553 octet truss lattice cubes were examined using resonant ultrasound spectroscopy (RUS) techniques and correlated to X-ray computed tomography measurements of strut continuity. RUS models and measurements identified an outlier lattice possessing lower effective modulus despite larger effective density. Here, experimental measurements of density and effective Young's modulus were treated as input parameters in a homogenous anisotropic continuum model. The continuum model exhibited good agreement with experimental RUS measurements, thus confirming the potential for a simplified approximation of the octet truss lattice.

42 ENGINEERING↗

DOE OSTI2022

Atoms to fibers: Identifying novel processing methods in the synthesis of pitch-based carbon fibers

Understanding and optimizing the key mechanisms used in the synthesis of pitch-based carbon fibers (CFs) are challenging, because unlike polyacrylonitrile-based CFs, the feedstock for pitch-based CFs is chemically hetero- geneous, resulting in complex fabrication leading to inconsistency in the final properties. In this work, we use molecular dynamics simulations to explore the processing and chemical phase space through a framework of CF models to identify their effects on elastic performance. The results are in excellent agreement with experiments. We find that density, followed by alignment, and functionality of the molecular constituents dictate the CF mechanical properties more strongly than their size and shape. Last, we propose a previously unexplored fabrication route for high-modulus CFs. Unlike graphitization, this results in increased sp 3 fraction, achieved via generating high-density CFs. In addition, the high sp 3 fraction leads to the fabrication of CFs with isometric compressive and tensile moduli, enabling their potential applications for compressive loading.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2022

Elevated temperature contact creep and friction of nickel-based superalloys using machine learning assisted finite element analysis

Nickel-based superalloys with superior thermochemical, mechanical, and tribological properties are highly utilized for critical components in several high temperature applications such as gas turbines and nuclear reactors. Inconel 617, in particular, is considered as one of the main candidate superalloys for tribo-components in very-high-temperature gas-cooled nuclear reactors. Recent findings indicate that this alloy grows unique surface oxide especially in a high-temperature helium environment with distinctive wear, friction, and contact properties. Here, this study investigates the high temperature contact area evolution and frictional behavior of Inconel 617 using finite element simulation and provides predictive models for the contact and friction performance at different normal loads, dwell times, and temperatures. High temperature helium-aged Inconel 617 top surface properties (up to 600 °C) are utilized along with a single asperity-based deformable elastic-plastic contact model under combined normal and tangential loading. Machine learning is used to assist the finite element results and to predict friction coefficient as well as contact area evolution. While a small difference is observed in the instantaneous friction coefficient (no dwell time) for all temperatures, friction coefficient increases considerably with dwell time. This shows that the effect of contact creep for longer dwell times significantly dominants the effect of high temperature variation in basic mechanical parameters such as modulus and yield strength. It is found that increasing temperature and dwell times lead to the friction coefficient increase, yet the dominance of dwell time effects decreases at higher temperatures and loads. While the analysis is presented for Inconel 617, the methodology is easy to be generalized and can be applied to other HT alloys.

36 MATERIALS SCIENCE↗

DOE OSTI2023

X-ray diffraction reveals two structural transitions in szomolnokite

Hydrated sulfates have been identified and studied in a wide variety of environments on Earth, Mars, and the icy satellites of the solar system. The subsurface presence of hydrous sulfur-bearing phases to any extent necessitates a better understanding of their thermodynamic and elastic properties at pressure. End-member experimental and computational data are lacking and are needed to accurately model hydrous, sulfur-bearing planetary interiors. In this work, high-pressure X-ray diffraction (XRD) and synchrotron Fourier-transform infrared (FTIR) measurements were conducted on szomolnokite (FeSO 4 ·H 2 O) up to ~83 and 24 GPa, respectively. This study finds a monoclinic-triclinic (C2/c to P1¯) structural phase transition occurring in szomolnokite between 5.0(1) and 6.6(1) GPa and a previously unknown triclinic-monoclinic (P1¯ to P2 1 ) structural transition occurring between 12.7(3) and 16.8(3) GPa. The high-pressure transition was identified by the appearance of distinct reflections in the XRD patterns that cannot be attributed to a second phase related to the dissociation of the P1¯ phase, and it is further characterized by increased H 2 O bonding within the structure. We fit third-order Birch-Murnaghan equations of state for each of the three phases identified in our data and refit published data to compare the elastic parameters of szomolnokite, kieserite (MgSO 4 ·H 2 O), and blödite (Na 2 Mg(SO 4 ) 2 ·4H 2 O). At ambient pressure, szomolnokite is less compressible than blödite and more than kieserite, but by 7 GPa both szomolnokite and kieserite have approximately the same bulk modulus, while blödite’s remains lower than both phases up to 20 GPa. Furthermore, these results indicate the stability of szomolnokite’s high-pressure monoclinic phase and the retention of water within the structure up to pressures found in planetary deep interiors.

58 GEOSCIENCES↗