Search NASA⌕ Search

SEARCH · Search NASA

Results for “volumetric deformation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Complete Equations of State for Cyclotetramethylene Tetranitramine

Complete equations of state for the volumetric deformation of the β-polymorph of cyclotetramethylene tetranitramine (β-HMX) have been derived from its Helmholtz free energy. Dispersion-corrected density functional theory calculations were used to compute the dependence of the frequencies of the vibrational normal modes on specific volume. Here, the normal mode frequencies were used directly to generate both a tabular equation of state and an approximate model where the sum over normal mode frequencies is replaced by a set of five Debye models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cohesive instability in elastomers: insights from a crosslinked Van der Waals fluid model

Abstract The resistance to volumetric deformations displayed by polymer networks is largely due to secondary and tertiary interactions between neighboring polymer chains. These interactions are both entropic and enthalpic in nature but are fundamentally different from the entropic forces that resist shearing in these networks. In this paper, we introduce a new depiction of elastomers as a crosslinked Van der Waals fluid. Starting from first principles, we develop constitutive equations that are implemented in a continuum model as well as a discrete network model. Our models predict that the failure of polymer networks may be driven by an instability in the underlying polymer bulk ‘fluid’ or by the breaking of polymer chains, depending on the loading path taken. The results of this study indicate that material failure in elastomers exposed to a purely triaxial state, such as in a poker chip experiment, may be driven by an entirely different mode of instability than those deformed in pure shear, such as in a uniaxial tension experiment.

Lamont, Samuel C.↗
DOE OSTI2021

Stress-Dependent Chemo-Mechanical Performance of Amorphous Si Anodes for Li-Ion Batteries upon Lithiation

Alloying-type anodes are significantly governed by their chemo-mechanical performance during the electrochemical cycling. The reaction-induced huge volumetric change of these anodes may cause material degradation and failure under mechanical constraints. Here, we investigate the stress-dependent lithiation behavior of amorphous Si (a-Si) anodes using molecular dynamics simulations. It is indicated that a-Si anodes can sustain higher hydrostatic stress than biaxial/uniaxial ones without the occurrence of mechanical failure. Thermodynamic and electrochemical calculations demonstrate that although the lithiation procedure also affects the thermodynamic stability of a-Si anodes, it is mainly dominated by the external mean stresses. Compressive stress is confirmed to destabilize a-Si anodes and further trigger their capacity fading. Compared with our atomistic simulations, previous continuum models underestimate the open-cell potentials of a-Si anodes, due to their ignored large volumetric deformation at higher stresses and Li concentrations. Finally, this computational study provides the intensive atomic-level understanding of the stress-dependent lithiation behavior of a-Si anodes.

25 ENERGY STORAGE↗
DOE OSTI2021

Volumetric Density-Equalizing Reference Map with Applications

The density-equalizing map, a technique developed for cartogram creation, has been widely applied to data visualization but only for 2D applications. Here, we propose a novel method called the volumetric density-equalizing reference map (VDERM) for computing density-equalizing map for volumetric domains. Given a prescribed density distribution in a volumetric domain in $\mathbb{R}^3$, the proposed method continuously deforms the domain, with different volume elements enlarged or shrunk according to the density distribution. With the aid of the proposed method, medical and sociological data can be visualized via deformations of 3D objects. The method can also be applied to adaptive remeshing and shape modeling. Furthermore, by exploiting the time-dependent nature of the proposed method, applications to shape morphing can be easily achieved. Experimental results are presented to demonstrate the effectiveness of the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Volumetric additive manufacturing of shape memory polymers

Shape memory polymers (SMPs) capable of recovering from a deformed state through heating were 3D printed using volumetric additive manufacturing, which enabled the layerless, single-step fabrication of self-standing tripod and actuating gripper structures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High tensile alloy of copper to mitigate current collector deformation in silicon electrodes for lithium-ion batteries

The volumetric changes of silicon electrodes, along with the strong adhesive properties of certain binders, can lead to plastic deformation of the current collector and create damage in the electrode coating. Here, in this study, we report a detailed study of silicon coatings on a high-tensile alloy (HTA) foil of copper with strength over twice that of conventional copper foils. The HTA current collectors with high mechanical strength can mitigate plastic deformation upon continuous cycling. At moderate areal capacities (2.5–3 mAh cm −2 ), conventional copper foils show significant wrinkling after only a few electrochemical cycles, whereas the HTA foils remain intact. We demonstrate viability of the HTA foils in large format xx6395 pouch cells, in which the HTA current collectors remain intact even at an areal capacity of 4.5 mAh cm −2 ; in contrast, wrinkles form in conventional copper current collectors increasing the likelihood of lithium plating. Computational studies show that stresses generated during cycling of silicon electrodes are very high in the current collector and at the current collector-coating interface, explaining the wrinkling of conventional Cu foils. Our studies highlight importance of current collector to solve the electrochemical and chemo-mechanical performance challenges associated with high-loading silicon electrodes.

Chemo-mechanical degradation↗
DOE OSTI2025

Stabilization of Preternatural Barium Oxidation States as an Unexpected Byproduct of β-Decay: Discovery of a New Halide Semiconductor Alloy

137 Cs has a wide range of roles in the nuclear industry. The solid material, safely encapsulated in CsCl as 137 CsCl, is stored as fission product waste from nuclear power production and legacy waste from nuclear weapons production; it has also served as a radiation source in food and sewage irradiators as well as medical devices. However, because of the solubility of the chloride salt and the relatively high specific activity of 137 Cs, damaged or broken capsules can lead to severe radiological accidents. Safe capsule design and material recycling are complicated by the unclear structural evolution during β-decay, which remains ambiguous due to the differing oxidation states of Cs (1+) and Ba (2+). Here, in this study, we use first-principles calculations to investigate the evolving structure–property relationships of Cs 1–x Ba x Cl during β-decay. Despite the well-established 2+ formal oxidation state of alkali-earth metals, we find that Ba 1+ can be stabilized in the form of a mixed-valence alloy at low concentrations. Specifically, we identify three regimes for the β-decay of 137 Cs into CsCl: Ba-doped CsCl (Ba ≤ 14%), wherein Ba has the expected 2+ oxidation state; Cs–Ba–Cl alloys, where Ba has a mix of the usual Ba 2+ and highly unusual Ba 1+ oxidation state in the form of a quasi-disordered mixed-valence alloy (Ba = 25%); and phase separation into a CsCl + BaCl 2 + Ba (m) mechanical mixture, where Ba reverts to its expected 2+ oxidation state (Ba > 25%). Surprisingly, the Cs 0.75 Ba 0.25 Cl mixed-valence alloy is a narrow indirect band gap semiconductor (1.05 eV) despite the insulating nature of both CsCl and BaCl 2 . It also exhibits strongly excitonic polarized optical properties, has glass-like ultralow thermal conductivity (directional average of 0.21 W/mK at 300 K), and shows greater resistance to deformation under both tensile and volumetric strengths compared with the original CsCl structure (e.g., shear and Young’s modulus of 9.04 and 31.62 GPa, respectively). These findings imply that transmutation of 137 Cs leads to highly unusual chemical bonding that stabilizes Ba 1+ in local regions of the quasi-disordered Cs 0.25 Ba 0.75 Cl, resulting in anomalous physical properties. Moreover, this discovery provides valuable insight for safe nuclear waste capsule design, which can aid in preventing environmental or human exposure to radioactive materials.

Fuhr, Addis S. [Oak Ridge National Laboratory (ORN↗
DOE OSTI2024

A soft departure from jamming: the compaction of deformable granular matter under high pressures

Here, the high-pressure compaction of three-dimensional granular packings is simulated using a bonded particle model (BPM) to capture linear elastic deformation. In the model, grains are represented by a collection of point particles connected by bonds. A simple multibody interaction is introduced to control Poisson's ratio and the arrangement of particles on the surface of a grain is varied to model both high- and low-frictional grains. At low pressures, the growth in packing fraction and coordination number follows the expected behavior near jamming and exhibit friction dependence. As the pressure increases, deviations from the low-pressure power-law scaling emerge after the packing fraction grows by approximately 0.1 and results from simulations with different friction coefficients converge. These results are compared to predictions from traditional discrete element method simulations which, depending on the definition of packing fraction and coordination number, may only differ by a factor of two. As grains deform under compaction, the average volumetric strain and asphericity, a measure of the change in the shape of grains, are found to grow as power laws and depend heavily on the Poisson's ratio of the constituent solid. Larger Poisson's ratios are associated with less volumetric strain and more asphericity and the apparent power-law exponent of the asphericity may vary. The elastic properties of the packed grains are also calculated as a function of packing fraction. In particular, we find the Poisson's ratio near jamming is 1/2 but decreases to around 1/4 before rising again as systems densify.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A subdivision-stabilized B-spline mixed material point method

Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Gas Bubble Evolution in Polycrystalline UMo Fuels Under Elastic-Plastic Deformation: A Phase-Field Model With Crystal-Plasticity

In monolithic UMo fuels, the interaction between the Al cladding and large gas bubble volumetric swelling causes both elastic-plastic and creep deformation. In this work, a phase-field model of gas bubble evolution in polycrystalline UMo under elastic-plastic deformation was developed for studying the dynamic interaction between evolving gas bubble/voids and deformation. A crystal plasticity model, which assumes that the plastic strain rate is proportional to resolved shear stresses of dislocation slip systems on their slip planes, was used to describe plastic deformation in polycrystalline UMo. Xe diffusion and gas bubble evolution are driven by the minimization of chemical and deformation energies in the phase-field model, while evolving gas bubble structure was used to update the mechanical properties in the crystal plasticity model. With the developed model, we simulated the effect of gas bubble structures (different volume fractions and internal gas pressures) on stress-strain curves and the effect of local stresses on gas bubble evolution. The results show that 1) the effective Young’s modulus and yield stress decrease with the increase of gas bubble volume fraction; 2) the hardening coefficient increases with the increase of gas bubble volume fraction, especially for gas bubbles with higher internal pressure; and 3) the pressure dependence of Xe thermodynamic and kinetic properties in addition to the local stress state determine gas bubble growth or shrinkage. The simulated results can serve as a guide to improve material property models for macroscale fuel performance modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Volumetric and shear strain localization throughout triaxial compression experiments on rocks

Deformation localization is a widely observed, but rarely quantified process in the crust. Recent observations suggest that the localization of seismicity and fracture networks can help identify the approach to catastrophic failure. Here, we quantify the localization processes of the volumetric and deviatoric strain components in twelve triaxial compression experiments imaged with X-ray tomography. We capture three-dimensional images of the rock cores during triaxial compressing toward failure, and then calculate the local strain components using digital volume correlation. The divergence and curl of the incremental displacement vector field provide the volumetric and deviatoric components of the strain field. We quantify localization using the proportion of the rock occupied by high magnitudes of the volumetric and deviatoric strains, and the Gini coefficient of these high magnitude strains, which measures the deviation from a uniform process. We find that the vast majority, but not all, of the experiments experience strain localization toward failure. The rocks typically experience their maximum degree of strain localization not immediately preceding failure, but on average at 90% of the failure stress. The volumetric strain tends to localize more than the deviatoric strain. Finally, these observations support using the localization of the volumetric strain, along with the deviatoric strain, to identify the evolution of the precursory phase preceding earthquakes.

58 GEOSCIENCES↗
DOE OSTI2021

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2024

Asymptotic quasisymmetric high-beta three-dimensional magnetohydrodynamic equilibria near axisymmetry

Quasisymmetry (QS), a hidden symmetry of the magnetic field strength, is known to support nested flux surfaces and provide superior particle confinement in stellarators. In this work, we study the ideal magnetohydrodynamic (MHD) equilibrium and stability of high-beta plasma in a large-aspect-ratio stellarator. In particular, we show that the lowest-order description of a near-axisymmetric equilibrium vastly simplifies the problem of three-dimensional quasisymmetric MHD equilibria, which can be reduced to a standard elliptic Grad–Shafranov equation for the flux function. We show that any large-aspect-ratio tokamak, deformed periodically in the vertical direction, is a stellarator with approximate volumetric QS. We discuss exact analytical solutions and numerical benchmarks. Finally, we discuss the ideal ballooning and interchange stability of some of our equilibrium configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

HFIR SiC-SiC Composite Clad Tube Bowing Test: Pre-Irradiation Characterization

This report describes the pre-irradiation characterization of the SiC-SiC composite clad tube bowing experiment that will be irradiated in the High Flux Isotope Reactor (HFIR). There are concerns that SiCSiC fuel cladding in light water reactors could undergo bowing because of the non-uniform fast neutron flux profiles. This experiment is combined with the SiC-SiC channel box experiment, which has similar concerns regarding control blade movements in boiling water reactors. The HFIR experiment aims to validate the predicted deformation and stresses in the composite clad specimen after the exposure of fast neutron flux gradients. Significant radial fast neutron flux gradients that exist in the permanent reflector of HFIR were thoroughly characterized by using detailed 3D neutronic calculations. The 3D displacement damage dose rate profile and the resulting volumetric swelling in SiC were used as inputs to structural analyses to determine the predicted deformation that will affect the specimens. The tube specimens were characterized by using traditional dimensional inspection and surface profilometry to provide detailed information regarding the pre-irradiation condition as the first step to evaluate the radiation-induced deformation. Furthermore, fine engraving markers were inscribed along all outer surfaces of the specimen and mapped by using a digital microscope and a multidimensional stage. This allows marker spacings to be accurately measured and compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. The experiment successfully completed two HFIR cycle irradiations, cycles 492 and 493, which finished on July 25, 2021.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Assessing continuum plasticity postulates with grain stress and local strain measurements in triaxially compressed sand

Critical state and continuum plasticity theories have been used in research and engineering practice in soil and rock mechanics for decades. These theories rely on postulated relationships between material stresses and strains. Some classical postulates include coaxiality between stress and strain rates, stress–dilatancy relationships, and kinematic assumptions in shear bands. Although numerical and experimental data have quantified the strains and grain kinematics in such experiments, little data quantifying grain stresses are available. Here, we report the first-known grain stress and local strain measurements in triaxial compression tests on synthetic quartz sands using synchrotron X-ray tomography and 3D X-ray diffraction. We use these data to examine the micromechanics of shear banding, with a focus on coaxiality, stress-dilatancy, and kinematics within bands. Our results indicate the following: 1) elevated deviatoric stress, strain, and stress ratios in shear bands throughout experiments; 2) coaxial principal compressive stresses and strains throughout samples; 3) significant contraction along shear bands; 4) vanishing volumetric strain but nonvanishing stress fluctuations throughout samples at all stages of deformation. Our results provide some of the first-known in situ stress and strain measurements able to aid in critically evaluating postulates employed in continuum plasticity and strain localization theories for sands.

42 ENGINEERING↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

HFIR SiC Bowing Test Ready to Insert

This report describes the successful assembly of a High Flux Isotope Reactor (HFIR) irradiation experiment designed to assess radiation-induced lateral bowing of silicon carbide fiber–reinforced, silicon carbide matrix composite (SiC/SiC) components under a radial fast neutron flux gradient. Excessive bowing of a SiC/SiC channel box in a boiling water reactor could potentially interfere with control blade movements. Similar concerns exist for SiC/SiC fuel cladding in light water reactors. The experiment described herein will provide experimental validation of the structural response of a miniature SiC/SiC channel box and tube specimens with pressurized water reactor diameters during irradiation. The significant radial fast neutron flux gradients that exist in the permanent reflector of HFIR were characterized using detailed three-dimensional neutronic calculations. The three-dimensional displacement damage dose rate profile and the resulting volumetric swelling in SiC were used as inputs to structural analyses that predicted the deformation and stresses in the channel box specimen. The specimens were thoroughly characterized prior to irradiation using traditional dimensional inspection and surface profilometry so that these measurements can later be compared with similar measurements that will be made post-irradiation to determine radiation-induced deformations. Furthermore, fine engraving markers were inscribed along all outer surfaces of the specimen and mapped using a digital microscope and a three-dimensional stage. This technique allowed for accurate measurements of the marker spacings, which can be compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. The experiment was successfully assembled and is scheduled for insertion during HFIR cycle 492, which is currently scheduled to run from May 25, 2021 to June 18, 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗