Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plastic behavior”

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 1,207 records · Page 67

Mercury Speciation via Diffusive Gradients Thin-Films Technology

Objectives: Develop and test various diffusive gradient in thin films (DGT) samplers for mercury - SRNL is developing reactive DGTS (rDGTs). Test chemistry to differentiate methylmercury from total and/or inorganic mercury in environmental samples. Fabricate and test rDGT samplers for deployment. Deploy the rDGTs in variable settings. Mercury is a persistent-bioaccumulative-toxic environmental pollutant. Mercury exists in different species such as elemental, inorganic, and methyl. Mercury speciation determine mercury behavior and toxicity. For example, methylmercury strongly accumulates in biota such as fish resulting in potential human health impacts. DGTs are innovative samplers to measure water concentration by diffusion and capture. Types of Hg DGTs used: Total Hg, Inorganic Hg, and Methyl Hg. Deployment Sites: Upper Three Runs, Lower Three Runs, Steele Creek, East Fork Poplar Creek (SRNL Sites), Horizons Creek, NOA Creek, Bruners Creek (ORNL Sites). DGT concentration estimates closely matched alternative measurements for soluble mercury at each Oak Ridge site - DGTs did not respond to particulate mercury. DGTs provide a representative measure for biota uptake since they are left out longer which allows for an average exposure concentration. DMA-80 provided an efficient and quick analysis. Variability in site location made a difference in Hg species. DGT preparations and sampling wasn't as quick and simple as microcolumn technology. Copper reagent degraded agarose gel - more work and cleaner data are needed to assess speciation in rDGTs. Future Work: Use cross-linked polyacrylamide for the collection layer for a stronger plastic. Make a copper diffusion layer without agarose. Test copper separation in lab using realistic stream conditions such as high organic carbon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Ti–6Al–4V using crystal plasticity, calibrated with multi-scale experiments, to understand the effect of the orientation and morphology of the α and β phases on time dependent cyclic loading

Classically, crystal plasticity modeling has used a range of constitutive equations, in which the incorporation of additional physics-based relationships typically results in additional model parameters. These additional parameters need to be reliably calibrated, which often necessitates the use of a range of experimental data acquired at multiple length scales. In this work, a crystal plasticity based finite element (CPFE) model for a dual-phase Titanium alloy, Ti–6Al–4V, is developed. The α and β phases of the microstructure are explicitly modeled. The model is calibrated using a systematic optimization routine and experimental data that consist of macroscopic stress-strain curves coupled with lattice strains on different crystallographic planes for the two phases. These experimental data were obtained from in situ high energy X-ray diffraction experiments for multiple material pedigrees, with varying crystallographic orientation distribution and β volume fractions. Depending on the thermomechanical-processing route and the heat treatment used to manufacture the alloy, Ti–6Al–4V can exist in a wide number of microstructural forms, which often results in the α and β phases either having well aligned slip systems (following the Burgers orientation relationship (BOR)) or possessing no alignment of the slip systems across the interphase boundary (not following the BOR). In this study, the fully-calibrated CPFE model is used to gain a comprehensive understanding of the deformation behavior of Ti–6Al–4V, specifically, the effect of microstructures that follow the BOR (or not) on time-dependent cyclic loading (including the effects of dwell hold times).

36 MATERIALS SCIENCE↗

Comprehensive Data Analysis and Analytic Method Development for PBX9501 Material (FY2019 Annual Report of Aging and Lifetimes Program)

In the polymer bonded explosive (PBX) 9501, a binder consisting of 2.5 wt% nitroplasticizer (NP) and 2.5 wt% Estane R 5703 (Estane) is combined with 94.9 wt% of the high-explosive HMX and 0.1 wt% stabilizer. Because of its flexibility and tensile strength, this binder lowers the sensitivity and improves the manufacturability of PBX 9501. However, like many plasticizers comprised of low molecular weight components, NP has a tendency to diffuse out of the PBX 9501 matrix and can decompose at moderate temperatures into reactive byproducts, such as NO, NO 2 , H 2 O, and HNO x . Through oxidation and hydrolysis, these molecules can further degrade NP and Estane, ultimately degrading the properties of PBX 9501. While gaseous molecules can readily diffuse out of the PBX 9501 charges, the intermediates with lower volatility are more likely trapped in the condensed phase inside the charges and, in a closed system, will co-exist with Estane for an extended period of time creating an ongoing reactive environment. Understanding the rates of formation for these volatiles and intermediates during the NP degradation is therefore a critical prerequisite for understanding the long-term stability of polymeric binders used in munition systems. To work toward this ultimate goal, in the past year, we have conducted comprehensive studies on the physical properties of NP and finished three-year long aging experiment to understand the aging behavior of NP under thermal treatment. Important results are summarized in this annual report. Although NP is widely used in the DOE complex, their physical properties are rather scattered and inconsistent in the open literature. For example, there are at least two widely different values, 14.5°C and -15°C, for the melting point of NP. Although it is known that NP is a eutectic 50:50 mixture of BDNPA and BDNPF, their eutectic phase diagram is not well documented. Furthermore, the effect of temperature on the miscibility between water and eutectic BDNPA/F mixture is rarely reported. To fill these knowledge gaps, a large set of BDNPA/F mixtures with BDNPA concentration ranging from 0 to 100 wt% was analyzed using DSC techniques. In addition to determining the eutectic melt point as -25°C, a phase diagram of the BDNPA/F system was constructed from -30°C to 45°C. With this phase diagram, the phase transition temperatures and composition can be readily found for the BDNPA/F mixtures with various mass ratios.

36 MATERIALS SCIENCE↗

Report on FY 2023 Research and Development on Specially Designed Creep-fatigue Experiments on Alloy 617 in Support of Improving Creep-fatigue Evaluation Approaches

Experimental and numerical studies in support of developing the integrated Elastic–Perfectly Plastic (EPP) plus Simplified Model Test (SMT) design methodology, referred to as the EPP+SMT method, continued in FY 2023. This report focuses on the methods for extrapolating the EPP+SMT creep-fatigue life curves at long hold times and low strain ranges at elevated temperatures. In this work, the available uniaxial creep-fatigue failure data on Alloy 617 at temperatures of 950°C and 850°C were analyzed to provide a guidance on the development of the extrapolation method and the creep-fatigue failure criteria. A viscoplastic constitutive model for Alloy 617 was adopted to extrapolate the mechanical responses to low strain ranges and long hold times. A set of design curves of Alloy 617 at temperatures of 950°C, 850°C, and 800°C with tensile hold times of 1 hr, 100 hr, and 1,000 hr are developed. Furthermore, creep-fatigue testing on two notch specimen geometries, shallow-notch and sharp V-notch, on Alloy 617 was performed 950°C to understand the multiaxial stress relaxation behavior. The experimental and numerical results on the notch specimens were compared with those on the standard uniaxial smooth bar specimens. The effect of multi-axial stress state combined with elastic follow-up on the stress relaxation behavior was investigated in this report.

36 MATERIALS SCIENCE↗

A preliminary characterization of the tensile and fatigue behavior of tungsten-fiber/Waspaloy-matrix composite

A microstructural study and a preliminary characterization of the room temperature tensile and fatigue behavior of a continuous, tungsten fiber, Waspaloy-matrix composite was conducted. A heat treatment was chosen that would allow visibility of planar slip if it occurred during deformation, but would not allow growth of the reaction zone. Tensile and fatigue tests showed that the failed specimens contained transverse cracks in the fibers. The cracks that occurred in the tensile specimen were observed at the fracture surface and up to approximately 4.0 mm below the fracture surface. The crack spacing remained constant along the entire length of the cracked fibers. Conversely, the cracks that occurred in the fatigue specimen were only observed in the vicinity of the fracture surface. In instances where two fiber cracks occurred in the same plane, the matrix often necked between the two cracked fibers. Large groups of slip bands were generated in the matrix near the fiber cracks. Slip bands in the matrix of the tensile specimen were also observed in areas where there were no fiber cracks, at distances greater than 4 mm from the fracture surface. This suggests that the matrix plastically flows before fiber cracking occurs.

Corner, Ralph E.↗

Dynamic Measurements of the Structural Evolution of Material Defects at the Mesoscale

This is a technical report to be published as part of the FY2021 annual overview package. We have designed experiments and developed analytical methods using images from x-ray and electron microscopy in order to quantify the behavior of crystal dislocation defects at their native nanometer-length and micro-length scales. The results of these analyses are used to directly inform theoretical physics models and guide interpretation of dark field x-ray microscopy (DFXM) images that offer spatial resolution surpassing the resolution of beamline optics. Additionally, we developed new analytical methods to automate the alignment of optics hardware used in both synchrotron and x-ray free electron laser beamline facilities. The tools and workflows developed for this project are now being deployed as proof-of-concept in reaction history analysis (RHA) for future programmatic integration.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Partial connectomes of labeled dopaminergic circuits reveal non-synaptic communication and axonal remodeling after exposure to cocaine

Dopaminergic (DA) neurons exert profound influences on behavior including addiction. However, how DA axons communicate with target neurons and how those communications change with drug exposure remains poorly understood. We leverage cell type-specific labeling with large volume serial electron microscopy to detail DA connections in the nucleus accumbens (NAc) of the mouse (Mus musculus) before and after exposure to cocaine. We find that individual DA axons contain different varicosity types based on their vesicle contents. Spatially ordering along individual axons further suggests that varicosity types are non-randomly organized. DA axon varicosities rarely make specific synapses (<2%, 6/410), but instead are more likely to form spinule-like structures (15%, 61/410) with neighboring neurons. Days after a brief exposure to cocaine, DA axons were extensively branched relative to controls, formed blind-ended ‘bulbs’ filled with mitochondria, and were surrounded by elaborated glia. Finally, mitochondrial lengths increased by ~2.2 times relative to control only in DA axons and NAc spiny dendrites after cocaine exposure. We conclude that DA axonal transmission is unlikely to be mediated via classical synapses in the NAc and that the major locus of anatomical plasticity of DA circuits after exposure to cocaine are large-scale axonal re-arrangements with correlated changes in mitochondria.

59 BASIC BIOLOGICAL SCIENCES↗

An atomistic study connecting underlying dislocation behavior with superior mechanical properties of NiCoCr medium entropy alloy

NiCoCr-based medium-entropy alloy (MEA) with a simple face-centered cubic crystal phase exhibits excellent mechanical properties, often attributed to the synergy of multiple deformation mechanisms. However, the atomistic origin of their outstanding mechanical response, including microstructural evolution and dislocation behavior under varying strain-rates and orientation, remains unclear. In this work, we employ large-scale molecular dynamics (MD) simulations to investigate the changes in deformation mechanisms along three distinct orientations ([110], [111], [100]) under varying strain rates (1 ×10 8 /sec, 1 ×10 10 /sec, 1 ×10 12 /sec) in the NiCoCr MEA. The presence of the stair-rod and the Shockley partial dislocations under uniaxial tensile strain are found to play a key role in the formation of deformation twinning and ε-martensite, which positively correlates with strain-rate dependent dislocation analysis. These findings further establish the role of the dislocations in controlling the superior mechanical response and excellent fracture toughness of the NiCoCr MEA. Systematic transmission-electron microscopy tests performed on the [111]-oriented crystals, deformed at different strain levels, at room temperature provide clear evidence of both the extended stacking-fault and the stair rods, confirming the predicted microstructural features. Finally, this study offers key insights into the complex nucleation mechanisms of deformation twinning and ε-martensite, such as twinning – and transformation–induced plasticity (TWIP-TRIP), providing valuable guidelines for studying similar material classes.

36 MATERIALS SCIENCE↗

Modeling of melting and erosion in high velocity microparticle impacts

High velocity collisions can cause melting, a situation involving elasticity, plasticity, and phase transformations in tandem, which is very challenging for mechanical modeling. Here, this paper advances the use of smoothed particle hydrodynamics (SPH) simulations for this complex situation, investigating stainless steel microspheres impacting a tin substrate. After extensive validation against a range of mechanical and geometrical measurements of experimental impacts, the model provides a detailed, quantitative view of melt volume production over a broad range of impact velocities. It also enables direct decomposition of the impact energy partitioning, including analysis of the energy absorbed by plasticity in the impactor, the system kinetic energy, and the adiabatic heating of the substrate in addition to the melting itself. Since these quantities are not available experimentally, SPH modeling substantially enhances the quantitative understanding of such impacts. Finally, the model speaks to the phase characteristics of the ejecta formed during impact, confirming the experimental inference that the ejected material likely contains a significant fraction of solid material.

Constitutive modeling↗

Guidelines for Proof Test Analysis

These guidelines integrate state-of-the-art elastic-plastic fracture mechanics (EPFM) and proof test implementation issues into a comprehensive proof test analysis procedure in the form of a road map which identifies the types of data, fracture mechanics based parameters, and calculations needed to perform flaw screening and minimum proof load analyses of fracture critical components. Worked examples are presented to illustrate the application of the road map to proof test analysis. The state-of-the art fracture technology employed in these guidelines is based on the EPFM parameter, J, and a pictorial representation of a J fracture analysis, called the failure assessment diagram (FAD) approach. The recommended fracture technology is validated using finite element J results, and laboratory and hardware fracture test results on the nickel-based superalloy Inconel 718, the aluminum alloy 2024-T3511, and ferritic pressure vessel steels. In all cases the laboratory specimens and hardware failed by ductile mechanisms. Advanced proof test analyses involving probability analysis and multiple-cycle proof testing (MCPT) are addressed. Finally, recommendations are provided on how to account for the effects of the proof test overload on subsequent service fatigue and fracture behaviors.

Chell, G. G.↗

Relaxor Ferroelectric-Like Spatiotemporal Memory in Field-Driven Lipid Bilayers

Lipid membranes are often regarded as passive barriers, yet their nonlinear dielectric response remains poorly understood. Using all-atom molecular dynamics, we show that fully hydrated dipalmitoylphosphatidylcholine bilayers exhibit relaxor ferroelectric-like behavior under time-dependent electric fields. Unlike crystalline relaxors, which are bipolar and display little remanent polarization, lipid bilayers exhibit a unipolar polarization response: even an alternating current field produces persistent, asymmetric polarization. Furthermore, the underlying free-energy landscape contains two distinct minima, a nonpolarized state and a unipolarly polarized state, between which stochastic thermally activated transitions occur. Directionally resolved Van Hove analysis reveals pronounced anisotropy arising from out-of-plane electric dipole alignment, interleaflet coupling, and lateral polarization domains. Each field cycle nucleates polarization at distinct sites and monitors their relaxation, marking a crossover from thermal fluctuations to field-sustained polarization. Remarkably, these polarized domains persist after field removal, generating long-lived, spatially coherent dipolar patterns that encode nanoscale polarization memory. Potassium chloride amplifies these effects via dielectric screening and a modified hydration structure, enhancing electric dipole flexibility and cooperativity. Together, these results establish protein-free bilayers as nonlinear, history-dependent dielectrics capable of sustaining field-tunable electromechanical coupling, providing an emergent physical foundation for nanoscale information storage and memory phenomena reminiscent of short- and long-term plasticity in soft neuromorphic systems.

Insulators↗

NASA Contractor Report: Guidelines for Proof Test Analysis

These Guidelines integrate state-of-the-art Elastic-Plastic Fracture Mechanics (EPFM) and proof test implementation issues into a comprehensive proof test analysis procedure in the form of a Road Map which identifies the types of data, fracture mechanics based parameters, and calculations needed to perform flaw screening and minimum proof load analyses of fracture critical components. Worked examples are presented to illustrate the application of the Road Map to proof test analysis. The state-of-the-art fracture technology employed in these Guidelines is based on the EPFM parameter, J, and a pictorial representation of a J fracture analysis, called the Failure Assessment Diagram (FAD) approach. The recommended fracture technology is validated using finite element J results, and laboratory and hardware fracture test results on the nickel-based superalloy IN-718, the aluminum alloy 2024-T351 1, and ferritic pressure vessel steels. In all cases the laboratory specimens and hardware failed by ductile mechanisms. Advanced proof test analyses involving probability analysis and Multiple Cycle Proof Testing (MCPT) are addressed. Finally, recommendations are provided on to how to account for the effects of the proof test overload on subsequent service fatigue and fracture behaviors.

Chell, G. G.↗

Thin Film Solar Cells: Organic, Inorganic and Hybrid

Thin film solar cells are an important developing resource for hundreds of applications including space travel. In addition to being more cost effective than traditional single crystal silicon cells, thin film multi-crystaline cells are plastic and light weight. The plasticity of the cells allows for whole solar panels to be rolled out from reams. Organic layers are being investigated in order to increase the efficiency of the cells to create an organic / inorganic hybrid cell. The main focus of the group is a thin film inorganic cell made with the absorber CuInS2. So far the group has been successful in creating the layer from a single-source precursor. They also use a unique method of film deposition called chemical vapor deposition for this. The general makeup of the cell is a molybdenum back contact with the CuInS2 layer, then CdS, ZnO and aluminum top contacts. While working cells have been produced, the efficiency so far has been low. Along with quantum dot fabrication the side project of this that is currently being studied is adding a polymer layer to increase efficiency. The polymer that we are using is P3OT (Poly(3-octylthiopene-2,5-diyll), retroregular). Before (and if) it is added to the cell, it must be understood in itself. To do this simple diodes are being constructed to begin to look at its behavior. The P3OT is spin coated onto indium tin oxide and silver or aluminum contacts are added. This method is being studied in order to find the optimal thickness of the layer as well as other important considerations that may later affect the composition of the finished solar cell. Because the sun is the most abundant renewable, energy source that we have, it is important to learn how to harness that energy and begin to move away from our other depleted non-renewable energy sources. While traditional silicon cells currently create electricity at relatively high efficiencies, they have drawbacks such as weight and rigidness that make them unattractive especially for space applications. Thin film photovoltaics have the potential to alleviate these problems and create a cheap and efficient way to harness the power of the sun.

Dankovich, John↗

Assessing the viability of a new subsize tensile specimen geometry for evaluation of structural nuclear and additively manufactured materials

The use of subsize specimens in nuclear materials testing has been a subject of ongoing interest due to radiation safety concerns and resource conservation needs. It is also of interest in additive manufacturing (AM), also known as 3D-printing, as subsize specimens can more accurately represent the behavior of the small, intricate geometries often produced using AM. A novel, extremely small geometry, called the Subsize Teeny (SST) was recently developed and is of interest for implementation. Given its extraordinarily small size and the complexities associated with subsize specimen testing, adequate vetting of this geometry is necessary to ensure data quality. Here, a variety of unirradiated nuclear structural materials were tested in the SST geometry and compared against the well-established SSJ3 geometry. In addition, two case studies implementing the SST as a screening geometry for AM materials were also conducted. The question of SST viability was found to be highly nuanced and will often be dependent on the context or application in question. It was determined, however, that the SST is a largely invalid geometry for exceptionally coarse-grained materials or in cases where the physical defect volume equals or exceeds 0.1 % or where the specimen machining parameters result in significant surface alterations. On the other hand, it was determined that the SST may be employed with confidence if the test material is nearly or totally free of physical defects, isotropic, demonstrates homogeneous plastic deformation, and possesses a fine-grained, nearly or totally homogeneous microstructure with at least twelve slip systems.

Additive manufacturing↗

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗

Knowledge gaps for neuromorphic ionic computing

BACKGROUND Neuromorphic computing, inspired by the human brain’s ability to process information efficiently, represents a transformative approach to computation. In this Review, we explore the emerging field of neuromorphic ionic computing, which leverages ionic conduction and coupling to mimic neural processes, and identify critical knowledge gaps that must be addressed to realize its full potential. A central theme of the discussion is energy efficiency, a challenge that is both a limitation and an opportunity for this technology. Although complementary metal-oxide semiconductor (CMOS)–based neuromorphic technologies have made strides in scaling to billions of neurons and are increasingly applied in artificial intelligence and numerical computing, they remain orders of magnitude behind the human brain in terms of connectivity and energy efficiency. Neuromorphic ionic computing promises to overcome these limitations by leveraging the distinct architectural and operational principles of the brain. Our brains achieve this energy efficiency by combining several key features: using the same network elements to store and process information; using an incredibly complex and massively interconnected three-dimensional (3D) network of locally active elements that enables sparsity, robustness in the presence of noise, adaptation, and life-long learning; computing at comparatively low voltage and frequency; and last, taking advantage of a plethora of ions and small molecules as information carriers. Here, we propose that ionic computing systems can take advantage of similar features to achieve substantial gains in energy efficiency. ADVANCES Since the first reports of neuromorphic ionic behavior in nanofluidic channels, we have witnessed an explosion of reports that used ionic devices to produce synaptomimetic behaviors. However, achieving the goals of ionic computing requires not only implementation of much more sophisticated device functionality but also overcoming fundamental barriers in materials science, device architecture, and system integration. Current ionic devices, even those incorporating state-of-the-art materials, still suffer from limited functionality and stability, which restrict their performance and increase energy demands. Developing new materials with enhanced ionic properties is essential to overcome these limitations. Similarly, the design of neuromorphic devices must evolve to leverage the particular advantages of ionic processes. Existing architectures often follow a single-information-carrier logic of conventional electronics or are constructed of mesoscale fluidics, failing to capitalize on the energy-efficient mechanisms inherent to ionic systems or implement the multiple-information-carrier paradigm. Current neuromorphic chips focus on large-scale networks of analog memory elements based on mechanisms such as charge trap (flash), filamentary, phase change, or spin, which are built on top of a network of artificial CMOS neurons. Although such prototype networks have achieved impressive performance, it is difficult to envision how they can implement the key features such as massive connectivity, sophisticated plasticity, adaptability, sparsity, and “multichromatic” computing. Although small-scale devices have demonstrated promising results, integrating them, maintaining energy efficiency, and implementing temperature control as systems grow in complexity and size to computationally relevant scale remain major hurdles. Furthermore, interfacing neuromorphic ionic devices with existing computing technologies presents technical and conceptual challenges that will require innovative approaches that combine insights from neuroscience, materials science, and engineering. OUTLOOK Despite these challenges, the potential impact of neuromorphic ionic computing is profound with potential applications ranging from artificial intelligence to robotics and beyond. We also argue that neuromorphic ionic computing systems should not, at least in the beginning, compete with CMOS technologies but rather should focus on applications that require extreme energy efficiency with chemical and/or biological compatibility, such as biomedical applications (for example, brain-computer interfaces), environmental monitoring, and agricultural and food applications. Ultimately, this Review highlights the crucial role of interdisciplinary collaboration in advancing the field. Neuromorphic ionic computing is not merely a technological innovation; it represents a substantial step toward sustainable computation, aligning with the growing demand for energy-conscious solutions in a world that is increasingly reliant on data and computation.

Neuromorphic↗

Computational thermomechanics for crystalline rock. Part II: Chemo-damage-plasticity and healing in strongly anisotropic polycrystals

We present a thermal–mechanical–chemical-phase field model that captures the multi-physical coupling effects of precipitation creeping, crystal plasticity, anisotropic fracture, and crack healing in polycrystalline rock at various temperature and strain-rate regimes. This model is solved via a fast Fourier transfer solver with an operator-split algorithm to update displacement, temperature and phase field, and chemical concentration incrementally. In nuclear waste disposal in salt formation, brine inside the crystal salt may migrate along the grain boundary and cracks due to the gradient of interfacial energy and pressure. This migration has a significant implication on the permeability evolution, creep deformation, and crack healing within rock salt but is difficult to incorporate implicitly via effective medium theories compared with computational homogenization. As such, we introduce a thermodynamic framework and a corresponding computational implementation that explicitly captures the brine diffusion along the grain boundary and crack at the grain scale. Meanwhile, the anisotropic fracture and healing are captured via a high-order phase field that represents the regularized crack region in which a newly derived non-monotonic driving force is used to capture the fracture and healing due to the solution–precipitation. Numerical examples are presented to demonstrate the capacity of the thermodynamic framework to capture the multiphysics material behaviors of rock salt.

42 ENGINEERING↗

Experimental Shock Transformation of Gypsum to Anhydrite: A New Low Pressure Regime Shock Indicator

The shock behavior of gypsum is important in understanding the Cretaceous/Paleogene event and other terrestrial impacts that contain evaporite sediments in their targets (e.g., Mars Exploration Rover Spirit detected sulfate at Gusev crater, [1]). Most interest focuses on issues of devolatilization to quantify the production of SO2 to better understand its role in generating a temporary atmosphere and its effects on climate and biota [2,3]. Kondo and Ahrens [4] measured induced radiation emitted from single crystal gypsum shocked to 30 and 40 GPa. They observed greybody emission spectra corresponding to temperatures in the range of 3,000 to 4,000 K that are a factor of 2 to 10 times greater than calculated pressure-density energy equation of state temperatures (Hugoniot) and are high enough to melt gypsum. Chen et al. [5] reported results of shock experiments on anhydrite, gypsum, and mixtures of these phases with silica. Their observations indicated little or no devolatilization of anhydrite shocked to 42 GPa and that the fraction of sulfur, by mass, that degassed is approx.10(exp -2) of theoretical prediction. In another report of shock experiments on calcite, anhydrite, and gypsum, Badjukov et al. [6] observed only intensive plastic deformation in anhydrite shock loaded at 63 GPa, and gypsum converted to anhydrite when shock loaded at 56 GPa but have not experimentally shocked gypsum in a step-wise manner to constrain possible incipient transformation effects. Schmitt and Hornemann [7] shock loaded anhydrite and quartz to a peak pressure of 60 GPa and report the platy anhydrite grains were completely pseudomorphed by small crystallized anhydrite grains. However, no evidence of interaction between the two phases could be observed and they suggested that recrystallization of anhydrite grains is the result of a solid-state transformation. They concluded that significant decomposition of anhydrite requires shock pressures higher than 60 GPa. Gupta et al. [8] reanalyzed the calcite and anhydrite shock wave experiments of Yang [9] using improved equations of state of porous materials and vaporized products. They determined the pressures for incipient and complete vaporization to be 32.5 and 122 GPa for anhydrite GPa which is a factor of 2 to 3 lower than reported earlier by Yang [9]. These studies are not in agreement regarding the onset of sulfate decomposition and documentation of shock effects in gypsum is incomplete.

Bell, Mary S.↗