Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Damage models”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A directional microcrack damage theory for brittle solids based on continuous hyperplasticity

We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.

Materials Science↗
DOE OSTI2024

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2020

Targeted inhibition of gut bacterial β-glucuronidase activity enhances anticancer drug efficacy

Irinotecan treats a range of solid tumors, but its effectiveness is severely limited by gastrointestinal (GI) tract toxicity caused by gut bacterial β-glucuronidase (GUS) enzymes. Targeted bacterial GUS inhibitors have been shown to partially alleviate irinotecan-induced GI tract damage and resultant diarrhea in mice. Here, we unravel the mechanistic basis for GI protection by gut microbial GUS inhibitors using in vivo models. We use in vitro, in fimo, and in vivo models to determine whether GUS inhibition alters the anticancer efficacy of irinotecan. We demonstrate that a single dose of irinotecan increases GI bacterial GUS activity in 1 d and reduces intestinal epithelial cell proliferation in 5 d, both blocked by a single dose of a GUS inhibitor. In a tumor xenograft model, GUS inhibition prevents intestinal toxicity and maintains the antitumor efficacy of irinotecan. Remarkably, GUS inhibitor also effectively blocks the striking irinotecan-induced bloom of Enterobacteriaceae in immune-deficient mice. In a genetically engineered mouse model of cancer, GUS inhibition alleviates gut damage, improves survival, and does not alter gut microbial composition; however, by allowing dose intensification, it dramatically improves irinotecan’s effectiveness, reducing tumors to a fraction of that achieved by irinotecan alone, while simultaneously promoting epithelial regeneration. These results indicate that targeted gut microbial enzyme inhibitors can improve cancer chemotherapeutic outcomes by protecting the gut epithelium from microbial dysbiosis and proliferative crypt damage.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Predicting Cement-Sheath Integrity with Consideration of Initial State of Stress and Thermoporoelastic Effects

In conventional wellbore-integrity analysis, the cement sheath’s initial state of stress and transient thermoporoelastic effects are often neglected. However, the initial state of stress is prerequisite information for accurately predicting the safe operating conditions that prevent a cemented well from being damaged. In addition, transient thermoporoelastic effects can have a profound effect on when damage will occur. Here in this paper, we propose a model that includes these effects to predict the safe operating pressures and temperatures that will prevent cement-sheath failure. For the initial state of stress, we proposed an empirical model using measurements. Subsequent stress changes are evaluated by a fully coupled transient thermoporoelastic model to analyze the mechanical behavior of the cement sheath. We predict the safe operating envelope (SOE) for shear, tensile, and debonding cement-sheath failures caused by pressure and temperature perturbations after the cement sets. Our model predicts that pore pressure is a key factor for cement failure, especially for rapid temperature changes. If the formation is low permeability, the transient pore pressures are amplified, causing the risk of damage to increase. Compared with conventional thermoelastic models, the thermoporoelastic model predicts a smaller SOE when heating the internal casing fluid and a larger envelope when cooling the internal casing fluid. Finally, the heating rate was considered with respect to field applications. The heating rate was also considered, and slower heating/cooling rates can prevent damage to the cement sheath. Finally, the thermoporoelastic model was applied to explain several laboratory and field experiments and achieved good matches.

02 PETROLEUM↗
DOE OSTI2021

NFD Methodology Report

The NetFlow Dynamics (NFD) model was developed for estimating the availability of a commodity supplied by a national- or regional-scale infrastructure following unexpected disruption of one or more of its components. The large scope of the disruptions of interest produce changes in availability lasting days to weeks. Consequently, the model does not resolve daily variations in system state and does not include the long-term processes that cause infrastructures to evolve as assets are added and removed according to owners ’planning decisions. NFD simulates fluid flow, including petroleum and other incompressible fluids, as well as natural gas and other compressible fluids, through pipeline networks characterized by limits on transmission capacity and storage. It was designed to enable efficient exploration of possible transmission system responses to large-scale disruptions lasting for days or longer. The model formulation reflects constraints on transmission and storage capacity imposed by the physical system assets. Those capacity limits are input parameters and are not derived from more basic system properties such as pipeline diameters and compressor power. A system’s response to a large disruption is controlled by operational decisions as well as damage to physical assets. The NFD model formulation allows users to efficiently consider alternative scenarios about the way remaining capacity might be used so that the analysis result appropriately reflects uncertainties about operator response.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Postearthquake Damage Mapping via Remote Sensing: Lessons From the 2023 Türkiye Disaster

This review addresses the urgent need for scalable, accurate, and reproducible remote sensing solutions following the February 2023 Türkiye earthquakes. It synthesizes the contributions of five peer-reviewed studies published in the IEEE JSTARS Special Issue on postearthquake damage and risk assessment. These studies cover areas such as damage classification with deep learning, fusion of multisource remote sensing data, creation of benchmark datasets, detailed damage mapping, and analysis of geophysical signals using outgoing longwave radiation. The article summarizes the methodological approaches and the practical relevance of the reviewed studies for detecting, evaluating, and quantifying damage, and outlines key challenges, including model generalization, class ambiguity, and data integration. It also discusses emerging trends, including explainable artificial intelligence, multimodal data fusion, and open-data platforms. This synthesis provides a foundation for building robust, interpretable, and real-time disaster response systems and aims to guide future research in earthquake-related Earth observation and rapid damage assessment.

Taskin, Gulsen [Istanbul Technical University] (OR↗
DOE OSTI2023

Virtual sensing of wind turbine hub loads and drivetrain fatigue damage, Virtuelle Sensoren für die Messung von Hauptwellenlasten und Ermüdungsschäden im Antriebstrang von Windenergieanlagen

Abstract: This paper presents a Digital Twin for virtual sensing of wind turbine aerodynamic hub loads, as well as monitoring the accumulated fatigue damage and remaining useful life in drivetrain bearings based on measurements of the Supervisory Control and Data Acquisition (SCADA) and the drivetrain condition monitoring system (CMS). The aerodynamic load estimation is realized with data-driven regression models, while the estimation of local bearing loads and damage is conducted with physics-based, analytical models. Field measurements of the DOE 1.5 research turbine are used for model training and validation. The results show low errors of 6.4% and 1.1% in the predicted damage at the main and the generator side high-speed bearing respectively. Zusammenfassung: In diesem Aufsatz wird ein digitaler Zwilling für Windenergieanlagen vorgestellt, welcher die virtuelle Erfassung der Hauptwellenlasten und die Zustandsüberwachung von Ermüdungschäden und der verbleibende Nutzungsdauer der Antriebsstranglager ermöglicht. Der digital Zwilling nutzt Messdaten des Supervisory Control and Data Acquisition (SCADA) Systems und des Zustandsüberwachungssystems des Antriebsstranges (CMS). Die Berechnung der Hauptwellenlasten ist mit datenbasierten Regressionsmodellen umgesetzt, während die Berechnung der Lagerkräfte und der Ermüdungsschaden mit physikbasierten Modelle durchgeführt wird. Für die Modellentwicklung und -validierung werden Feldmessdaten der DOE 1.5MW Turbine eingesetzt. Die Abweichungen in den Ermüdungsschäden am Hauptwellenlager und am Generatorwellenlager betragen lediglich 6,4% beziehungsweise 1,1%.

17 WIND ENERGY↗
DOE OSTI2024

Statistical analysis of displacement damage in small devices from neutron and ion irradiation

Modern semiconductor devices, such as gate-all-around nanosheet field-effect transistors (GAA NS FETs), are smaller than displacement damage cascades from fission neutrons. In this regime, device failure may occur through low-probability single events, rather than by parametric degradation previously seen in larger devices. Here, we present a statistical model that predicts the probability of a damage event in a small device and the probability distribution of the magnitude, i.e., number of displacements within the device, from each event. The model is developed first for neutron irradiation and then for energetic ion irradiation. The model is consistent with results from recent experiments in which lithium-ion irradiation produced stepwise increases in subthreshold current in GAA NS FETs.

Wampler, W. R.↗
DOE OSTI2025

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗
DOE OSTI2024

Bayes_Opt-SWMM: A Gaussian process-based Bayesian optimization tool for real-time flood modeling with SWMM

Real-time flood model plays a pivotal role in averting urban flood damage, particularly when there is minimal lead time for preparatory measures. However, urban flood modeling in real-time often contends with inherent uncertainties arising from input data uncertainty and parameter ambiguities. Here this study introduces a real-time calibration (RTC) tool called Bayes_Opt-SWMM, specifically tailored for real-time urban flood modeling and uncertainty optimization. This tool leverages the Gaussian process-based Bayesian optimization algorithm and interfaces seamlessly with the Stormwater Management Model (SWMM). It integrates real-time model forcing data and flood monitoring collected through sensors and gauges which are strategically placed within critical locations of urban drainage systems. Our approach hinges on the Surrogate Model based Uncertainty Optimization (SMUO) concept, providing an avenue for enhancing real-time flood modeling. Bayes_Opt-SWMM runs the optimization process using a surrogate model called Gaussian Process emulator with two inference methods: (1) the Gaussian Process (GP) model and (2) Markov Chain Monte Carlo (MCMC) algorithm in GP model (GP_MCMC). Furthermore, three acquisition functions, namely Expected Improvement (EI), Maximum Probability of Improvement (MPI), and Lower Confidence Bound (LCB), facilitate optimal parameter fitting within the surrogate models. The efficiency of GP-based surrogate models in learning SWMM model parameters, leads to an improved uncertainty quantification and accelerated real-time flood modeling in urban areas. Overall, Bayes_Opt-SWMM emerges as a cost-effective and valuable tool for real-time flood modeling and monitoring, with significant potential for managing intelligent storm water systems in urban environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

High Temperature Erosion Modeling in Particle Based CSP Systems

1. Introduction. Wear and erosion damage of materials from solid particle and surface interactions is a major issue in various industries. Although more common in mining and oil and gas production, erosion is becoming a critical issue in renewable technologies as well such as particle based concentrated solar thermal power (CSP) systems. In particle based CSP systems, solid particles are used to absorb solar energy and as thermal storage. However, these particles may cause significant amount of wear to system components while through the system at high temperatures. This damage can be costly, and therefore, requires a greater understanding of solid particle erosion in CSP systems. Throughout the years, models and tools have been developed to predict and control erosion in industries such as oil and gas production. However, these erosion models and erosion prediction tools have been mainly developed based on erosion data for much higher velocities and lower temperatures, i.e. the operating conditions for which the existing models have been built are not pertinent to those expected in CSP system. It is known that erosion depends on many factors such as material properties, erodent particle properties, and particle impact speed and angle. However, the effect of temperature on erosion is not vastly investigated. In this work, a temperature-based correlation is introduced that will be used to modify the existing erosion models to predict erosion rates at velocities and temperatures relevant to Gen3 CSP systems. The existing models are first validated against erosion experiments run at low temperatures and low velocities. Subsequently, a temperature correction term is developed that can extend the existing models to high temperatures, based on the available experimental data showing the effect of temperature on impact erosion. 2. Erosion Models. Continuing sequence of impacts from solid particles on surfaces would result in loss of material due to mechanical interaction between solid surface and particles. Erosion equations have been developed to predict erosion under different conditions and erosion mechanism, including cutting and deformation erosion. These models are mainly either mechanistic, empirical, and semi-mechanistic models, when the latter combine the theory of the erosion mechanism and particle motion with the available experimental data. One of the first empirical correlation was introduced by American Petroleum Institute (API) Recommended Practice (RP) 14E [1]. This correlation was very conservative, and several improved empirical and semi-mechanistic equations were developed in the years following. More recently, the models introduced by Erosion/Corrosion Research Center (E/CRC) at the University of Tulsa and Arabnejad. et al. [2] are frequently used in the literature and industry, as they account for many parameters affecting erosion including particle impact speed and angle, material density and hardness, and particle size and shape. However, both these models were developed based on data at relatively high velocities and also do not take into account any mechanical changes in material associated with thermal cycling at high temperatures. 3. Erosion Prediction at High Temperatures. In this work, the accuracy of Arabnejad et al. model and E/CRC model is investigated to predict impact erosion at conditions relevant to Gen3 CSP systems. Two significant deviations expected in Gen3 systems compared to operating conditions used to build these models are particle impact velocities and system temperature. As a first step, the performance of these models was validated at low velocities and temperatures. Impact erosion experiments were conducted on SS316 coupons using HSP 40/70 ceramic particles. Assuming particle rate of 1 kg/s/m, 0.0254 m of particle curtain thickness, and particle velocity of 1.5 m/s in the system, an overall erosion of 0.59 mm/year was calculated. The experimental results were subsequently compared to those from computational simulations and erosion of 0.47 mm/year was obtained using the Arabnejad et al. model. The results indicated that the models work well at low velocities and low temperature conditions. To account for temperature effect, a mathematical correlation was developed using data published by DUCOM [3] for Inconel 600 eroded by alumina particles at high velocities. The correlation was applied to both Arabnejad et. al model and the E/CRC model. The prediction results from these modified models were within 20% of this experimental data. Predictions of of erosion by the modified E/CRC model at three different temperatures are made. Similar to the calculation at low temperature, a typical CSP system with 1 kg/s/m of particle flow rate per unit length of the particle curtain and a curtain thickness of 0.0254 m is assumed. The annual thickness loss calculations were run for different particle impact velocities, assuming a uniform particle impact area equal to the cross-sectional area of the curtain (particle-particle interactions and dispersion of particles are not considered in obtaining the results). It is observed that erosion increases exponentially as temperature increases. Furthermore, it is also noted that, erosion changes non-linearly with impact velocity. Currently, experiments are also being conducted to measure erosion of SS316 at 800 ºC with HSP 40/70 particles at low impact velocities. We expect to use the results from high temperature testing to further improve the temperature correlation function. Similar models are also being developed for abrasion erosion resulting from particle sliding along the surfaces as well as attrition from particle to particle and particle to surface interactions. References. [1] Institute, A. P. (1991). API Recommended Practice for Design and Installation of Offshore Production Platform Piping System, API RP 14E. [2] Arabnejad, H., Mansouri, A., Shirazi, S. A., and McLaury, B. S. (2015a). Development of mechanistic erosion equation for solid particles. Wear, 332–333, 1044–1050. http://doi.org/10.1016/j.wear.2015.01.031. [3] https://ducom.com/high-temperature-erosion-evaluating-sample-wear/

14 SOLAR ENERGY↗
DOE OSTI2025

Effect of Loading Frequency on Fatigue Behavior of Additively Manufactured Short Fiber Thermoplastics via Infrared Thermography

Fracture, Fatigue, Failure, Damage Evolution and Thermomechanics & Infrared Imaging, Volume 4 of the Proceedings of the 2025 SEM Annual Conference & Exposition on Experimental and Applied Mechanics, the fourth volume of five from the Conference, brings together contributions to this important area of research and engineering. The collection presents early findings and case studies on a wide range of areas, including: • 3D Deformation and Failure • In-situ Techniques and Microscale Effects on Mechanical Behavior • Vibration and Cyclic Fatigue • Integration of Models and Experiments • Interface-Mediated Damage and Failure • Fracture and Fatigue in Brittle Materials • Material Characterizations Using Thermography • Fatigue, Damage & Fracture Evaluation Using Infrared Thermography

Pathak, Pharindra [Auburn University]↗
DOE OSTI2021

Accelerated Creep Test Qualification of Creep-Resistance Using the Wilshire–Cano–Stewart Constitutive Model and Stepped Isostress Method

In this paper, a qualification of accelerated creep-resistance of Inconel 718 is assessed using the novel Wilshire–Cano–Stewart (WCS) model and the stepped isostress method (SSM) and predictions are made to conventional creep data. Conventional creep testing is a long-term continuous process; in fact, the ASME B&PV III requires that 10,000+ h of experiments must be conducted to each heat for materials employed in boilers and/or pressure vessel components. This process is costly and not feasible for rapid development of new materials. As an alternative, accelerated creep testing techniques have been developed to reduce the time needed to characterize the creep resistance of materials. Most techniques are based upon the time-temperature-stress superposition principle that predicts minimum-creep-strain-rate (MCSR) and stress-rupture behaviors but lack the ability to predict creep deformation and consider deformation mechanisms that occur for experiments of longer duration. The SSM has been developed, which enables the prediction of creep deformation response as well as reduce the time needed for qualification of materials. The SSM approach has been successful for polymer, polymeric composites, and recently has been introduced for metals. In this study, the WCS constitutive model, calibrated to SSM test data, qualifies the creep resistance of Inconel 718 at 750 °C and predictions are compared to conventional creep testing data. The WCS model has proven to make long-term predictions for stress-rupture, MCSR, creep deformation, and damage in metallic materials. The SSM varies stress levels after time interval adding damage to the material, which can be tracked by the WCS model. The SSM data is calibrated into the model and the WCS model generates realistic predictions of stress-rupture, MSCR, damage, and creep deformation. The calibrated material constants are used to generate predictions of stress-rupture and are postaudit validated using the National Institute of Material Science database. Similarly, the MCSR predictions are compared from previous studies. Finally, the creep deformation predictions are compared with real data and is determined that the results are well in between the expected boundaries. Material characterization and mechanical properties can be determined at a faster rate and with a more cost-effective method. This is beneficial for multiple applications such as in additive manufacturing, composites, spacecraft, and industrial gas turbines.

42 ENGINEERING↗
DOE OSTI2020

Microstructural characterization of cold-worked 316 stainless steel flux thimble tubes irradiated up to 100 dpa in a commercial Pressurized Water Reactor

Two flux thimble tubes (FTT) made of 15% cold-worked 316 stainless steel (SS) were harvested from Ringhals Pressurized Water Reactor (PWR) Unit 2, with peak damages of 76 and 100 displacements per atom (dpa) after 29 and 34 years’ service, respectively. Specimens sectioned from parent tubes were comprehensively characterized with nominal damage levels of ~0, ~41, ~74, 76, and 100 dpa at a nominal temperature range of 285-323 °C. Both FTTs contained helium and hydrogen gases as transmutation products. The helium follows a production rate of ~9.8 appm/dpa, while environmental factors complicate hydrogen production obscuring an exact H/dpa ratio. Irradiation-induced dislocation loops, nano-cavities, solute clusters, and microsegregation were all observed. The dislocation loops and nano-cavities indicated saturation at 41 dpa. The solute clusters continued to evolve with Ni-Si clusters formed at 41 dpa, and Ni-Si-Mn-P clusters formed at 74 and 100 dpa, but neither clusters exhibited distinct diffraction patterns at any damage levels. Solute clusters were observed to frequently be co-located with dislocation loops, but fully decorated loops were rarely detected. Significant radiation-induced segregation (RIS) was observed around grain boundaries at all damage levels. The modified inverse Kirkendall (MIK) model captured the RIS behavior of major elements. Large cavities within or around an Mn-S rich region were observed for the first time. Through all the damage levels, void swelling is always below 0.05%, making significant dimensional change unlikely in core internals when used at similar conditions. Meanwhile, the role of overwhelming nanocavities, presumably helium bubbles, should be considered in other potential degradation mechanisms, including irradiation-assisted stress corrosion cracking, embrittlement, and loss of fracture toughness, which remain the concerns for extended operation of nuclear power plants.

Post Irradiation Examination, Isotope Dilution Mas↗
DOE OSTI2022

Multiscale mechanical design of the lightweight, stiff, and damage-tolerant cuttlebone: A computational study

Cuttlebone, the endoskeleton of cuttlefish, offers an intriguing biological structural model for designing low-density cellular ceramics with high stiffness and damage tolerance. Cuttlebone is highly porous (porosity ~93%) and lightweight (density less than 20% of seawater), constructed mainly by brittle aragonite (95 wt%), but capable of sustaining hydrostatic water pressures over 20 atmospheres and exhibits energy absorption capability under compression comparable to many metallic foams (~4.4 kJ/kg). Here, in this work, we computationally investigate how such remarkable mechanical efficiency is enabled by the multiscale structure of cuttlebone. Using the common cuttlefish, Sepia Officinalis, as a model system, we first conducted high-resolution synchrotron micro-computed tomography (µ-CT) and quantified the cuttlebone's multiscale geometry, including the 3D asymmetric shape of individual walls, the wall assembly patterns, and the long-range structural gradient of walls across the entire cuttlebone (ca. 38 chambers). The acquired 3D structural information enables systematic finite-element simulations, which further reveal the multiscale mechanical design of cuttlebone: at the wall level, wall asymmetry provides optimized energy absorption while maintaining high structural stiffness; at the chamber level, variation of walls (number, pattern, and waviness amplitude) contributes to progressive damage; at the entire skeletal level, the gradient of chamber heights tailors the local mechanical anisotropy of the cuttlebone for reduced stress concentration. Our results provide integrated insights into understanding the cuttlebone's multiscale mechanical design and provide useful knowledge for the designs of lightweight cellular ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radiation damage and dose limits in serial synchrotron crystallography at cryo- and room temperatures

Radiation damage limits the accuracy of macromolecular structures in X-ray crystallography. Cryogenic (cryo-) cooling reduces the global radiation damage rate and, therefore, became the method of choice over the past decades. The recent advent of serial crystallography, which spreads the absorbed energy over many crystals, thereby reducing damage, has rendered room temperature (RT) data collection more practical and also extendable to microcrystals, both enabling and requiring the study of specific and global radiation damage at RT. Here, we performed sequential serial raster-scanning crystallography using a microfocused synchrotron beam that allowed for the collection of two series of 40 and 90 full datasets at 2- and 1.9-Å resolution at a dose rate of 40.3 MGy/s on hen egg white lysozyme (HEWL) crystals at RT and cryotemperature, respectively. The diffraction intensity halved its initial value at average doses ( D 1/2 ) of 0.57 and 15.3 MGy at RT and 100 K, respectively. Specific radiation damage at RT was observed at disulfide bonds but not at acidic residues, increasing and then apparently reversing, a peculiar behavior that can be modeled by accounting for differential diffraction intensity decay due to the nonuniform illumination by the X-ray beam. Specific damage to disulfide bonds is evident early on at RT and proceeds at a fivefold higher rate than global damage. The decay modeling suggests it is advisable not to exceed a dose of 0.38 MGy per dataset in static and time-resolved synchrotron crystallography experiments at RT. This rough yardstick might change for proteins other than HEWL and at resolutions other than 2 Å.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗