SEARCH · Engineering Papers

Results for “tensor sampling”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

NASA NTRSDate not supplied

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

NASA NTRS2014

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗

DOE OSTI2020

Complete Strain Mapping of Nanosheets of Tantalum Disulfide

Quasi-two-dimensional (quasi-2D) materials hold promise for future electronics because of their unique band structuresthat result in electronic and mechanical properties sensitive to crystal strains in all three dimensions. Quantifying crystal strain is a prerequisite to correlating it with the performance of the device, and calls for high resolution but spatially resolved rapid characterization methods. Here we show that using fly-scan nano X-ray diffraction we can accomplish a tensile strain sensitivity below 0.001% with a spatial resolution of better than 80 nm over a spatial extent of 100 µm on quasi 2D flakes of 1T-TaS2. Coherent diffraction patterns were collected from a ~100 nm thick sheet of 1T-TaS2 by scanning 12keV focused X-ray beam across and rotating the sample. We demonstrate that the strain distribution around micron and sub-micron sized ‘bubbles’ that are present in the sample may be reconstructed from these images. The experiments use state of the art synchrotron instrumentation, and will allow rapid and non-intrusive strain mapping of thin film samples and electronic devices based on quasi 2D materials.

nano X-ray diffraction↗

DOE OSTI2021

Multi-modal imaging of a single mouse brain over five orders of magnitude of resolution

Mammalian neurons operate at length scales spanning six orders of magnitude; they project millimeters to centimeters across brain regions, are composed of micrometer-scale-diameter myelinated axons, and ultimately form nanometer scale synapses. Capturing these anatomical features across that breadth of scale has required imaging samples with multiple independent imaging modalities. Translating between the different modalities, however, requires imaging the same brain with each. Here, we imaged the same postmortem mouse brain over five orders of spatial resolution using MRI, whole brain micrometer-scale synchrotron x-ray tomography ($\mu$CT), and large volume automated serial electron microscopy. Using this pipeline, we can track individual myelinated axons previously relegated to axon bundles in diffusion tensor MRI or arbitrarily trace neurons and their processes brain-wide and identify individual synapses on them. This pipeline provides both an unprecedented look across a single brain's multi-scaled organization as well as a vehicle for studying the brain's multi-scale pathologies.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Overview of Hydraulic Fracturing Test Site 2 in the Permian Delaware Basin (HFTS-2)

Here, the Hydraulic Fracturing Test Site 2 (HFTS-2) is a large collaborative field-based R&D program in the Permian Delaware Basin, funded by the US Department of Energy through the National Energy Technology Laboratory (NETL) and the E&P industry, with support from academia. The projects' main objective is to improve the understating of the hydraulic fracturing process through utilization of advanced diagnostics and collection of through-fracture cores to provide undisputable evidence and attributes of the created hydraulic fractures. At the HFTS-2, in excess of $30 million was used to perform hydraulic fracturing research focusing on the Wolfcamp formation at a field site hosted and operated by Occidental. In addition to the research data collected by the project, Occidental provided a significant amount of background data for about a dozen existing wells in the test area as well as access to previously collected core. Additional technical and laboratory support was provided by the program members. Building on learnings and unanswered questions from HFTS-1 in the Permian Midland basin, the HFTS-2 used eight new producing wells and two existing (parent) wells to perform hydraulic fracturing research. Multiple science wells were drilled to sample and characterize the subsurface, including the collection of 540 feet of core in a vertical pilot hole and 948 feet of high-angle through-fracture core. The project installed permanent fiber optic cables in 3 wells to monitor near wellbore signals during fracturing and to collect cross-well strain measurements. Additional advanced diagnostics included a significant formation evaluation program on the vertical whole core, multi array moment tensor inversion capable microseismic survey, multi-well time-lapse geochemistry analysis, analysis of proppant distribution in producing child and slant core well, and others. We will provide an overview of the HFTS-2 project, including list of the consortium members, details of the test site, experiments performed, and technologies tested.

58 GEOSCIENCES↗

NASA NTRSDate not supplied

Enhanced Quantification of 3D Woven Composites Via Fourier Analysis and Structure Tensors Applied to CT Scans

NASA's development of advanced 3D woven carbon composites, such as HEEET and 3MDCP, represents a significant leap in materials designed to endure extreme atmospheric entry conditions. However, the internal complexity of these materials introduces a variety of structures and novel defect types. Traditional characterization methods rely on labor-intensive, manual examination of X-ray CT scans, demanding considerable time and expertise. This process, while thorough, is not sustainable for large-scale analysis and defect tracking. Addressing this critical challenge, this presentation introduces an innovative suite of automated tools and analytical techniques, pioneered by the Mars Sample Return (MSR) Earth Entry System (EES) Thermal Protection System (TPS) project. These advancements offer a more efficient and accurate approach for quantifying CT data and characterizing the defects of 3D woven composites.

Magnus Haw↗

NASA NTRSDate not supplied

Quantitative Characterization of 3D Woven Composites Via Fourier Analysis and Structure Tensors Applied to CT Scans

3D Woven↗

DOE OSTI2024

Exploiting universal nonlocal dispersion in optically active materials for spectro-polarimetric computational imaging

Recent years have seen significant advancements in exploring novel light-matter interactions such as hyperbolic dispersion within natural crystals. However, current studies have predominantly concentrated on local optical response of materials characterized by a dielectric tensor without spatial dispersion. Here, we investigate the nonlocal response in optically-active crystals with screw symmetries, revealing their lossless, super-dispersive properties compared to traditional optical response functions. We leverage this universal nonlocal dispersion, i.e. the dispersion of optical rotatory power, to explore a novel spectral de-multiplexing scheme compared to conventional gratings, prisms and metasurfaces. We design and demonstrate an ‘Nonlocal-Cam’ - a camera that exploits nonlocal dispersion through sampling of polarized spectral states and the application of computational spectral reconstruction algorithms. The Nonlocal-Cam captures information in both laboratory and outdoor field experiments which is unavailable to traditional intensity cameras - the spectral texture of polarization. Merging the fields of nonlocal electrodynamics and computational imaging, our work paves the way for exploiting nonlocal optics of optically active materials in a variety of applications, from biological microscopy to physics-driven machine vision and remote sensing.

Wang, Xueji [Purdue Univ., West Lafayette, IN (Uni↗

DOE OSTI2020

Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows

Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Investigating Quantum Materials with Half-Polarized Diffraction and magnetic PDF analysis at the HB-2A Neutron Powder Diffractometer

Local magnetic ordering and anisotropy is often central to the emergent behavior and subsequent functional properties in quantum materials and beyond. Neutron powder diffraction provides a straightforward yet extremely powerful technique for quantitative measurements of microscopic magnetic properties. The HB-2A powder diffractometer located at the High Flux Isotope Reactor in ORNL is traditionally utilized for long-range magnetic structure determination. Recently these capabilities have been extended to include methods aimed at accessing local magnetism: Half- polarized neutron powder diffraction (pNPD) and magnetic pair distribution function (mPDF) analysis. These two distinct techniques are possible on HB-2A due to the versatility of the instrument’s reciprocal space coverage, resolution and novel ultra-low temperature multi-sample changers that operate down to dilution refrigerator temperatures. This provides unique capabilities not found on any powder diffraction instrument and is particularly well suited to investigations of magnetic quantum materials. The development and implementation of these techniques will be discussed with a series of science case examples ranging from geometric frustrated magnets to magnetic metal-organic frameworks. Data reduction and analysis tools will be presented that enable the extraction of the local site susceptibility tensor and local spin-spin correlations in real space. Finally, potential combinations of these techniques in the form of half-polarized magnetic pair distribution function (pmPDF) analysis will be considered. Looking forward, HB-2A is undergoing a detector upgrade that will be in the user program by 2026. This will offer an order of magnitude increase in count rates to further aid the development of these often low signal measurements and provide new scientific capabilities.

Neutron Scattering↗

DOE OSTI2025

Machine learning interatomic potential for predicting the thermal properties of uranium nitride

We present a combined computational and experimental investigation of the thermal properties of uranium nitride (UN), focusing on the development of a machine learning interatomic potential (MLIP) using the moment tensor potential framework. The MLIP was trained on density functional theory (DFT) data and validated against various quantities including energies, forces, elastic constants, phonon dispersion, and defect formation energies, achieving excellent agreement with DFT calculations, prior experimental results, and our thermal conductivity measurement. The potential was then employed in molecular dynamics simulations to predict key thermal properties such as melting point, thermal expansion, specific heat, and lattice thermal conductivity. To further assess model accuracy, we fabricated a UN sample and performed new thermal conductivity measurements representative of single-crystal properties, which showed strong agreement with the MLIP predictions. This work confirms the reliability and predictive capability of the developed potential for determining the thermal properties of UN.

36 - MATERIALS SCIENCE↗

NASA NTRS2019

An Experimental Investigation of a Wing-Fuselage Junction Model in the NASA Langley 14- by 22-Foot Subsonic Tunnel

Current turbulence models, such as those employed in Reynolds-averaged Navier-Stokes CFD, are unable to reliably predict the onset and extent of the three-dimensional separated flow that typically occurs in wing-fuselage junctions. To critically assess, as well as to improve upon, existing turbulence models, experimental validation-quality flow-field data in the junction region is needed. In this report, we present an overview of experimental measurements on a wing-fuselage junction model that addresses this need. The experimental measurements were performed in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The model was a full-span wing-fuselage body that was configured with truncated DLR-F6 wings, both with and without leading-edge extensions at the wing root. The model was tested at a fixed chord Reynolds number of 2.4 million, and angles-of-attack ranging from -10 degrees to +10 degrees were considered. Flow-field measurements were performed with a pair of miniature laser Doppler velocimetry (LDV) probes that were housed inside the model and attached to three-axis traverse systems. One LDV probe was used to measure the separated flow field in the trailing-edge junction region. The other LDV probe was alternately used to measure the flow field in the leading-edge region of the wing and to measure the incoming fuselage boundary layer well upstream of the leading edge. Both LDV probes provided measurements from which all three mean velocity components, all six independent components of the Reynolds-stress tensor, and all ten independent components of the velocity triple products were calculated. In addition to the flow-field measurements, static and dynamic pressures were measured at selected locations on the wings and fuselage of the model, infrared imaging was used to characterize boundary-layer transition, oil-flow visualization was used to visualize the separated flow in the leading- and trailing-edge regions of the wing, and unsteady shear stress was measured at limited locations using capacitive shear-stress sensors. Sample results from the measurement techniques employed during the test are presented and discussed.

Kegerise, Michael A.↗

DOE OSTI2021

QM7-X, a comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for ≈4.2 million equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures—comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers)—as well as 100 non-equilibrium structural variations thereof to reach a total of ≈4.2 million molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Complex deuteron NMR signals

To determine the spin polarization of deuterons, nuclear magnetic resonance (NMR) is used. This is necessary for polarized targets, such as for the upcoming A zz and b 1 experiment at Jefferson Lab. NMR measures the impedance of a solenoid around a deuterated sample. Although the impedance is a complex value, conventionally only the real part of the impedance has been used for this purpose. However, often the tune is not precisely real, meaning the signal has at least some small imaginary portion. This conventionally has been dealt with by an offset parameter, such as Dulya’s false asymmetry method. For vector polarization, this suffices, as the tuning is factored into the overall error of the results, and for a small phase angle doesn’t make much of a difference. However, for tensor polarization, the exact lineshape of the signal is quite significant, and treating the impedance as complex during analysis removes the need for an offset parameter. As a result, it also provides more accurate results, as the conventional false asymmetry method over- or underestimates polarization, depending on the sign of the phase angle.

McClellan, Michael [University of New Hampshire, D↗

DOE OSTI2021

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

DOE OSTI2023

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Physics Division Strategic Plan Fiscal Years 2020-2024

The vision of the Physics Division (PHY) at Argonne National Laboratory is to continue enhancing its role as a world-leading institution in basic nuclear physics research and its applications. Key to this vision is for PHY to continue to safely and effectively operate and evolve the capabilities of the Argonne Tandem Linac Accelerator System (ATLAS) facility to best serve its users. ATLAS is the Department of Energy (DOE) accelerator facility for low-energy nuclear physics research. The research carried out by PHY covers many themes in contemporary science but can be distilled into five areas of focus. These themes are interconnected, with continuously evolving synergies between the various groups and facilities in PHY and the broader Laboratory. This five-year strategic plan serves to illustrate our current capabilities and new directions related to these five areas of focus (ATLAS also develops an independent strategic plan): Accelerator research and design. The goal of this theme is to design, fabricate, test, and implement novel accelerator systems, with a focus on high-intensity ion and electron systems. These R&D activities have led to enhancements in the ATLAS accelerator system. Among others, the division’s accelerator systems are in use or planned for use at Fermilab, the Advanced Photon Source, and in a broad variety of applications at the Facility for Rare Isotope Beams (FRIB). Atom trapping and fundamental symmetries. The goal of this theme is to explore and exploit the uses of advanced laser cooling and trapping techniques to manipulate atoms. There are three main areas of focus. First is in the application of the atom trap trace analysis (ATTA) technique for age determination of groundwater and ice by using radio-krypton dating. The next two trapping-based programs involve tests of fundamental symmetries in nature: the cooling and trapping of radium-225 with the aim of determining limits on an observation of its electric dipole moment and precision measurements of the beta decay properties of helium-6 to set limits on the tensor coupling constant. This is complemented by measurements of similar properties in lithium-8 and boron-8. Nuclear astrophysics. The goal of this theme is to enhance our understanding of how elements are created in the universe via explosive nucleosynthesis and how stars evolve. To meet the challenges of this theme, many of the capabilities of ATLAS are being enhanced, including the development of new beams through a new in-flight separator (RAISOR) and the anticipated neutron-generator upgrade of the Californium Rare Isotope Breeder Upgrade facility (nuCARIBU). These are coupled to state-of-the-art instruments such as the Canadian Penning Trap, Helical Orbit Spectrometer, Gammasphere, GRETINA, Multi-Sampling Ionization Chamber, the Fragment Mass Analyzer, and Argonne Gas-filled Fragment Analyzer, and a new low-background experimental area for decay studies. Nuclear structure. The goal of this theme is to understand the structure of nuclei, both stable and radioactive, in terms of single-particle properties, their shapes, and the dynamics governing reactions between them. These include questions such as what are the limits of nuclear stability and what are the properties of super heavy nuclei. As with the nuclear astrophysics theme, the capabilities of ATLAS, guided by the ATLAS user community, are continuously being enhanced to this end. The Division has strategic initiatives to play a leading role in the development of instrumentation and research programs at the FRIB. Quantum chromodynamics (QCD) and hadron physics. The goal of this theme is to lead major research programs focused on revealing the quark and gluon structure of protons, neutrons, nuclei, and short-lived mesons and baryons. These involve major programs at Jefferson Lab, Fermilab, and smaller facilities. They are complemented by theoretical endeavors centered on the question of how hadrons and their properties emerge from QCD. Out of these activities arise strategic initiatives to play a major role in the forthcoming Electron-Ion Collider (EIC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗