SEARCH · Search NASA

Results for “deformation function”

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 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Effect of proteoglycans at interfaces as related to location, architecture, and mechanical cues

Covalently bound functional GAGs orchestrate tissue mechanics through time-dependent characteristics. The role of specific glycosaminoglycans (GAGs) at the ligament-cementum and cementum-dentin interfaces within a human periodontal complex were examined. Matrix swelling and resistance to compression under health and modeled diseased states was investigated. The presence of keratin sulfate (KS) and chondroitin sulfate (CS) GAGs at the ligament-cementum and cementum-dentin interfaces in human molars (N = 5) was illustrated by using enzymes, atomic force microscopy (AFM), and AFM-based nanoindentation. Furthermore, the change in physical characteristics of modeled diseased states through sequential digestion of keratin sulfate (KS) and chondroitin sulfate (CS) GAGs was investigated. One-way ANOVA tests with P < 0.05 were performed to determine significant differences between groups. Additionally, the presence of mineral within the seemingly hygroscopic interfaces was investigated using transmission electron microscopy. Immunohistochemistry (N = 3) indicated presence of biglycan and fibromodulin small leucine rich proteoglycans at the interfaces. Digestion of matrices with enzymes confirmed the presence of KS and CS GAGs at the interfaces by illustrating a change in ti ssue architecture and mechanics. A significant increase in height (nm), decrease in elastic modulus (GPa), and tissue deformation rate (nm/s) of the PDL-C attachment site (215 ± 63-424 ± 94 nm; 1.5 ± 0.7-0.4 ± 0.2 GPa; 21 ± 7-48 ± 22 nm/s), and cementum-dentin interface (122 ± 69-360 ± 159 nm; 2.9 ± 1.3-0.7 ± 0.3 GPa; 18 ± 4-30 ± 6 nm/s) was observed. The sequential removal of GAGs indicated loss in intricate structural hierarchy of hygroscopic interfaces. From a mechanics perspective, GAGs provide tissue recovery/resilience. Our results provide insights into the role of GAGs toward conserved tooth movement in the socket in response to mechanical loads, and modulation of potentially deleterious strain at tissue interfaces.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Observation of nanoparticle coalescence during core-shell metallic nanowire growth in colloids via nanoscale imaging

The surface morphology and shape of crystalline nanowires significantly influence their functional properties, including phonon transport, electrocatalytic performance, to name but a few. However, the kinetic pathways driving these morphological changes remain underexplored due to challenges in real-space and real-time imaging at single-particle and atomic resolutions. This study investigates the dynamics of shell (Au, Pd, Pt, Fe, Cu, Ni) deposition on AuAg alloy seed nanowires during core-shell formation. By using chiral/non-chiral seed nanowires, advanced imaging techniques, including liquid-phase transmission electron microscopy (LPTEM), cryogenic TEM, and three-dimensional electron tomography, a three-step deposition process is revealed: heterogeneous nucleation, nanoparticle attachment, and coalescence. It is found that colloidal Ostwald ripening, metal reactivity, and deposition amount modulate nanoparticle size and surface roughness, shaping final morphologies. Noble metal nanoparticles (Au, Ag, Pd, Pt) coalesce with seed nanowire along the 〈111〉 direction, distinct from that of other metals. These findings are consistent across different metals, including Ru, Cu, Fe, and Ni, highlighting the hypothesis of these processes in nanowire formation. These findings enhance traditional crystallographic theories and provide a framework for designing nanowire morphology. Additionally, our imaging techniques may be applied to investigate phenomena like electrodeposition, dendrite growth in batteries, and membrane deformation.

Yang, Dahai↗

DOE OSTI2024

Methodological guidelines on concrete degradation based on predictive models and the release of MOSAIC for industry use

In light-water reactors (LWRs), the concrete biological shield (CBS) designates the concrete structure directly facing the reactor pressure vessel (RPV). The primary function of the CBS is to protect equipment and personnel from the neutron and gamma radiation exiting the RPV. Most CBSs in operation in the United States also provide structural support to the RPV. Based on a long-term structural analysis, the structural integrity of the CBS to transfer the in-service passive load to the foundation, dynamic loading during an earthquake, and thermal loading during a loss-of-coolant accident (LOCA) must be maintained. Based on test reactor data for accelerated conditions, the mechanical properties of irradiated concrete are affected when exposed to neutron fluence higher than ≈ 10 19 n.cm -2 (E > 0.1 MeV). At 80 years of operation, the fluence estimates at the surface of the CBS range between 1×10 19 and 7×10 19 n.cm -2 (E > 0.1 MeV): the fluence value depends mostly on the design. Although fluence attenuates through the CBS by virtue of the shielding properties of concrete, about 5% to 10% of the CBS wall depth is affected by high fluence. The main degradation mechanism is the so-called radiation-induced volumetric expansion (RIVE) caused by the amorphization of aggregate-forming minerals. RIVE causes large deformation of up to 18% in pure quartz, and it also reduces the structural properties of concrete. The structural effects of long-term exposure to fast neutron irradiation in LWRs was studied in the report entitled “Assessment of the Effect of the Irradiation-Induced Degradation on the In-Service Structural Performance of the Concrete Biological Shields” (M3LW-23OR0403013). This report illustrates that (1) irradiation-induced damage expands further than the region subject to a cumulative fluence of 10 19 n.cm -2 (E > 0.1 MeV) to reach the steel reinforcement location, and (2) the irradiation-induced damage region shows no substantial residual bearing capacity and should not be accounted for in the integrity assessment of the CBS.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2019

Static subspace approximation for the evaluation of G 0 W 0 quasiparticle energies within a sum-over-bands approach

Many-body perturbation theory within the GW approach has been established as a quantitatively accurate approach for predicting the quasiparticle and excited-state properties of a wide variety of materials. However, the successful application of the method is often complicated by the computational complexity associated with the evaluation and inversion of the frequency-dependent dielectric matrix ε(ω). Here, we describe an approach to speed up the evaluation of the frequency-dependent part of ε(ω) in the traditional sum-over-states GW framework based on the low-rank approximation of the static dielectric matrix, a technique often used in GW implementations that are based on a starting mean field within density-functional perturbation theory. We show that the overall accuracy of the approach, independently from other calculation parameters, is solely determined by the threshold on the eigenvalues of the static dielectric matrix, ε(ω = 0) , and that it can yield orders-of-magnitude speed-ups in full-frequency G W calculations. In this work, we validate our implementation with several benchmark calculations ranging from bulk materials to systems with reduced dimensionality, and show that this technique allows one not only to study larger systems, but also to carefully consider the convergence of computationally demanding systems, such as ZnO, without relying on plasmon-pole models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Evidence for Rigid Triaxial Deformation in Ge 76 from a Model-Independent Analysis

An extensive, model-independent analysis of the nature of triaxial deformation in 76 Ge, a candidate for neutrinoless double-beta (0νββ) decay, was carried out following multistep Coulomb excitation. Shape parameters deduced on the basis of a rotational-invariant sum-rule analysis provided considerable insight into the underlying collectivity of the ground-state and γ bands. Both sequences were determined to be characterized by the same β and γ deformation parameter values. In addition, compelling evidence for low-spin, rigid triaxial deformation in 76 Ge was obtained for the first time from the analysis of the statistical fluctuations of the quadrupole asymmetry deduced from the measured E 2 matrix elements. Furthermore, these newly determined shape parameters are important input and constraints for calculations aimed at providing, with suitable accuracy, the nuclear matrix elements relevant to 0νββ.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2017

Estimation of Core Recovery Time with TRACE

Following the severe accident in Fukushima Daiichi Nuclear Power Plant in Japan, safety related concerns about the zirconium based claddings used in current light water reactors have increased. High corrosion resistance and low neutron absorption cross section make zirconium a suitable cladding material for normal nuclear plant operational conditions. However, for nuclear accident scenarios when safety functions of the plant are lost, Zr claddings undergo exothermic reaction with steam and exacerbate the accident conditions. The exothermic reaction increases the temperature of the cladding and accelerates core melting in severe accident progression. And the generation of explosive hydrogen gas accumulates and challenges the structural integrity. The accident tolerant fuel (ATF) concept is the fuel and cladding system that gives a longer time margin, compared to standard Zr claddings, to cope with the accident when core cooling systems are lost in a nuclear power plant. The ATF cladding should have slower oxidation kinetics and less hydrogen generation than Zr to alleviate consequences of cladding interaction with steam at high temperatures. The maximum time available to mitigate progression of severe accident is an important measure to evaluate benefits of ATF claddings. For this work, the maximum time available to mitigate progression of severe accident is selected as the time at which cladding melting reaction starts. Melting reaction is selected as an upper limit because of the degradation of coolable geometry. Maintaining the coolable geometry is a regulation criteria for both loss of coolant (LOCA) analyses and outside the LOCA analyses. Once the melting of the cladding starts, it exacerbates the deformation of the fuel assembly and hottest rods within the core cannot be cooled by water reflood and results in fuel melting progression. The aim of this work is to evaluate time margin gain of ATF cladding in a severe accident scenario of LOCA without safety injections and station blackout considering oxidation kinetics only. Also it is assumed that safety injections are recovered before melting of the claddings start and the sensitivity of safety injection actuation time on core recovery is analyzed. The core recovery time is the time when makeup water from safety injections quench the fuel rods and halt the progression of the severe accident. This time is estimated by initiating safety injections at different times and observing the hot rod behavior which represents 1/5. of the core in this work. It is considered that core is not recoverable when 20% of the core is damaged.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Accelerated corrosion testing of cold spray coatings on 304L in chloride environments

Cold spray is an advanced metal manufacturing technique applied across many fields for a wide range of functions. Low heat input and compressive stresses induced into the substrate by the cold spray process makes it a promising choice for protective corrosion resistant coatings. One potential application for cold spray is as a protective coating against corrosion for spent nuclear fuel (SNF) interim dry storage canisters. As these canisters are currently stored at interim storage locations longer than originally intended, chloride induced stress corrosion cracking has been identified as a high priority knowledge gap, specifically with respect to prolonging or extending canister lifetimes ( Teague et al., 2019 ). The high deployability of cold spray, for which nozzles have been developed for application in constrained spaces, in conjunction with beneficial properties inherent to cold spray makes this a good candidate for a corrosion protection coating on SNF canisters. This work explores a pathway to rapidly down-select cold spray coatings for canisters by focusing on the corrosion properties. Specifically, this study examines the corrosion protection abilities of nickel and nickel-based alloy cold spray coatings on 304 L stainless steel in chloride rich environments through electrochemical scans and ferric chloride pitting tests (ASTM G48 Method A). It was shown that the porosity of the coating, the processing gas, material selection, and deformation in the substrate all impact the corrosion behavior of cold spray coatings and are areas where optimization could reduce potential materials degradation, enabling enhanced coatings development.

Karasz, Erin K.↗

DOE OSTI2019

Role of Well Operations and Multiphase Geomechanics in Controlling Fault Stability During CO 2 Storage and Enhanced Oil Recovery

CO 2 injection in active oil fields is a technology proposed for industrializing carbon storage operations. However, monitoring CO 2 migration, oil recovery, and mechanical deformation within caprock and on faults is a challenging problem in large-scale storage-enhanced oil recovery operations due to intersecting faults, multiphase flow effects, and wells with complex production-injection schedule. We develop and demonstrate a methodology based on our coupled modeling framework to monitor the movement of CO 2 , hydrocarbons and water, and the associated evolution in mechanical stability of faults during CO 2 storage-enhanced oil recovery in the Farnsworth Unit oil field in Texas, United States. The methodology honors geological, geophysical, and production-injection data acquired in the field over six decades. Differential depletion from hydraulically isolated fault compartments followed by water and CO 2 injection-induced overpressure causes volumetric contraction and expansion of the reservoir, and changes in the total and effective stresses in the overburden-reservoir-underburden complex. CO 2 migrates upward to accumulate near top of the geologic structure, and water migrates downward to pressurize the faults. Three-dimensional changes in the pressure and stress fields in the system lead to changes in the shear and effective normal tractions on three major faults compartmentalizing the field. Evolution in fault tractions is used to compute the evolution in the Coulomb failure function of the faults to quantify induced-slip tendency under production and injection. We explain the spatial heterogeneity and time variability of Coulomb failure function in terms of well location heterogeneity and well rate variability.

58 GEOSCIENCES↗

DOE OSTI2016

Radiation forces and torque on a rigid elliptical cylinder in acoustical plane progressive and (quasi)standing waves with arbitrary incidence

This paper presents two key contributions; the first concerns the development of analytical expressions for the axial and transverse acoustic radiation forces exerted on a 2D rigid elliptical cylinder placed in the field of plane progressive, quasi-standing, or standing waves with arbitrary incidence. The second emphasis is on the acoustic radiation torque per length. The rigid elliptical cylinder case is important to be considered as a first-order approximation of the behavior of a cylindrical fluid column trapped in air because of the significant acoustic impedance mismatch at the particle boundary. Based on the rigorous partial-wave series expansion method in cylindrical coordinates, non-dimensional acoustic radiation force and torque functions are derived and defined in terms of the scattering coefficients of the elliptic cylinder. A coupled system of linear equations is obtained after applying the Neumann boundary condition for an immovable surface in a non-viscous fluid and solved numerically by matrix inversion after performing a single numerical integration procedure. Computational results for the non-dimensional force components and torque, showing the transition from the progressive to the (equi-amplitude) standing wave behavior, are performed with particular emphasis on the aspect ratio a/b, where a and b are the semi-axes of the ellipse, the dimensionless size parameter, as well as the angle of incidence ranging from end-on to broadside incidence. The results show that the elliptical geometry has a direct influence on the radiation force and torque, so that the standard theory for circular cylinders (at normal incidence) leads to significant miscalculations when the cylinder cross section becomes non-circular. Moreover, the elliptical cylinder experiences, in addition to the acoustic radiation force, a radiation torque that vanishes for the circular cylinder case. The application of the formalism presented here may be extended to other 2D surfaces of arbitrary shape, such as Chebyshev cylindrical particles with a small deformation, stadiums (with oval shape), or other non-circular geometries.

42 ENGINEERING↗

DOE OSTI2012

Second order gauge invariant measure of a tidally deformed black hole

In this paper, a Lagrangian perturbation theory for the second order treatment of small disturbances of the event horizon in Schwarzchild black holes is introduced. The issue of gauge invariance in the context of general relativistic theory is also discussed. The developments of this paper is a logical continuation of the calculations presented in [1], in which the first order coordinate dependance of the intrinsic and exterinsic geometry of the horizon is examined and the first order gauge invariance of the intrinsic geometry of the horizon is shown. In context of second order perturbation theory, It is shown that the rate of the expansion of the congruence of the horizon generators is invariant under a second order reparametrization; so it can be considered as a measure of tidal perturbation. A generally non-vanishing expression for this observable, which accomodates tidal perturbations and implies nonlinear response of the horizon, is also presented.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Quasiparticle random-phase approximation calculations for M 1 transitions with the noniterative finite-amplitude method and application to neutron radiative capture cross sections

Here, we derive the equations of quasiparticle random-phase approximation (QRPA) based on the finite amplitude method with the Hartree-Fock+Bardeen–Cooper–Schrieffer (HF+BCS) single-particle states, and calculate the magnetic dipole (M⁢1) transition for deformed gadolinium isotopes. Our QRPA calculation shows both large spin-flip transitions in the 5 to 10 MeV excitation energy and the low energy orbital transition that would correspond to the M⁢1 scissors mode observed in nuclear experiments. Then, we calculate neutron capture reactions based on the statistical Hauser-Feshbach theory with the photoabsorption cross sections of even-even nuclei given by QRPA. We find that the capture cross section is enhanced due to the contribution from the low energy M⁢1 transition although the calculated capture cross section still underestimates the experimental data. This issue in the calculated capture cross section could be improved by uncertainties of the low energy E⁢1 transition neglected in our QRPA calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Failure Analysis for Molten Salt Thermal Energy Storage Tanks for In-Service CSP Plants

Thermal Energy Storage (TES) is a fundamental component in concentrating solar power (CSP) plants to increase the plant's dispatchability, capacity factor, while reducing the levelized cost of electricity. In central receivers CSP plants, nitrate molten salts have been used for several years for operation temperatures of up to 565 degrees C. Despite many efforts to advance nitrate salt to higher operation temperatures (even considering a replacement with molten chloride salts) to achieve higher energy conversion efficiencies, the 565 degrees C temperature is currently considered the state-of-the art. Although molten salt tanks have been broadly deployed in commercial CSP plants worldwide, several failures have been reported in these tanks after a few months or years of operation, causing significant economic loss and mistrust in CSP technologies. Most of these failures are associated with the infancy of the technology and multiple issues related to tank design, fabrication, commissioning, and aggressive operation. A technical standard dedicated to the design and fabrication of molten nitrate TES tanks does not exist today. Current in-service molten salt tanks have been generally designed based on the American Petroleum Institute's (API) 650 and ASME Section II standards. The API 650 code provides guidelines for dimensions and fabrication for oil storage tanks up to 260 degrees C. The ASME standard provides allowable stress values for various materials at a range of temperatures and conditions. Both standards seem to be limited for molten salt TES tanks where high temperatures, thermal cycling, and transient conditions are expected. In 2020, NREL released the Concentrating Solar Power Best Practices Study (NREL/TP-5500-75763) that summarized multiple issues in CSP plants, along with potential alternatives and recommendations to address those issues based on information collected from participants representing about 80% of operating CSP plants in the world. One of the recommendations from this study was the development of accurate and validated models to evaluate the plant's transient operation, capable of capturing the effect of short-term clouds and operator response, while being flexible in being adapted to various spatial and temporal resource data. The "Failure Analysis for Molten Salt Thermal Energy Tanks for In-Service CSP Plants" project was inspired on this recommendation and was focused on (1) the development and validation of a physics-based model for a representative, commercial-scale molten salt tank, (2) performing simulations to evaluate the behavior of the tank as a function of typical plant operation conditions, (3) understanding tank failures mechanisms, (4) determining the residual stress and distortion in the tank floor after welding fabrication and evaluating their impact in the stresses developed in the tank during operation, (5) assessing the impact of key operation parameters on the temperature and stress distribution, (6) conduct a preliminary evaluation of design features to reduce stress and improve tank's reliability, and (7) estimate the tank's service life based on the stress developed under diverse operation scenarios. From the analysis conducted in the project and presented in this report, it was found that maximum stresses surpassing the yield strength point of the stainless steel (SS) 347H are developed on the tank floor near the perimeter. These large stresses are strongly influenced by the initial residual stresses and distortion of the tank floor after welding fabrication. During operation, large stresses are developed in the tank floor at high operation temperatures with large salt inventory levels during transient operation. High stresses are also related to elevated temperature gradients in the tank floor that could be attributed to insufficient mixing within the salt inflow and the salt inventory. Based on the analysis, creep is the predominant failure mechanism. However, the large stress levels could favor the plastic deformation into buckles, and crack formation due to stress relaxation cracking during cycle operation. A lifetime below 3 years was estimated for the typical plant operation conditions and a specific initial residual stress and deformation distribution of the tank floor. The estimated life agrees with the service time to failure reported in several commercial molten salt tanks. Desing and operation guidelines can be extracted from the analysis presented in this report, which could be adopted by tank manufacturers and CSP operators to advance toward an ultimate solution for tank failures by reducing residual and operational stresses to achieve a tank service life of more than 30 years. Addressing failures in molten salt TES tanks is fundamental for the CSP industry's survivability, but it is also important for other industrial and power generation applications using this technology, including nuclear and concentrating solar thermal.

14 SOLAR ENERGY↗

DOE OSTI2019

GW170817 implications on the frequency and damping time of $f$-mode oscillations of neutron stars

Within a minimum model for neutron stars consisting of nucleons, electrons, and muons at $β$ equilibrium using about a dozen equations of state (EOS) from microscopic nuclear many-body theories and 40000 EOSs randomly generated using an explicitly isospin-dependent parametric EOS model for high-density neutron-rich nucleonic matter within its currently known uncertainty range, we study correlations among the $f$-mode frequency, its damping time and the tidal deformability as well as the compactness of neutron stars. Except for quark stars, both the $f$-mode frequency and damping time of canonical neutron stars are found to scale with the tidal deformability independent of the EOSs used. Applying the constraint on the tidal deformability of canonical neutron stars $Λ_{1.4} = 190^{+390}_{-120}$ extracted by the LIGO+VIRGO Collaborations from their improved analyses of the GW170817 event, the $f$-mode frequency and its damping time of canonical neutron stars are limited to 1.67–2.18 kHz and 0.155–0.255 s, respectively, providing a useful guidance for the ongoing search for gravitational waves from the $f$-mode oscillations of isolated neutron stars. Finally, assuming either or both the $f$-mode frequency and its damping time will be measured precisely in future observations with advanced gravitational wave detectors, we discuss how information about the mass and/or radius as well as the still rather elusive nuclear symmetry energies at suprasaturation densities may be extracted.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2019

Shape-preserving machining produces gradient nanolaminate medium entropy alloys with high strain hardening capability

A high density of grain boundaries can potentially increase structural materials' strength, but at the expense of losing the materials' strain hardening ability at high flow stress levels. However, endowing materials with grain size gradients and a high density of internal interfaces can simultaneously increase the strength and strain hardening ability. This applies particularly for through-thickness gradients of nanoscale interface structures. In this work, we apply a machining method that produces metals with nanoscale interface gradients. Conventional bulk plastic deformation such as rolling, a process applied annually to about 2 billion tons of material, aims to reduce the metal thickness. We have modified this process by introducing severe strain path changes, realized by leading the sheet through a U-turn while preserving its shape, an approach known as ‘hard turning’. Furthermore, we applied this process at both room temperature and 77 K to a NiCrCo medium entropy alloy. Micropillar compression was conducted to evaluate the mechanical response. After hard turning at room temperature, the surface microstructure obtained a ~50% increase in yield stress (0.9 GPa) over the original state with homogeneous grain size (0.4 GPa), but the initial strain hardening rate did not show significant improvement. However, after hard turning at 77 k, the gradient nanolaminate structure tripled in yield stress and more than doubled its initial strain hardening rate. The improvements were achieved by introducing a specific microstructure that consists of gradient nanolaminates in the form of nanospaced twins and martensite in the face center cubic (fcc) phase. This microstructure was formed only at cryogenic temperature. It was found after turning at room temperature that only nanospaced twins were present in the fcc phase inside nanolaminates that had formed at the surface. The origin of the enhanced strain hardening mechanism was studied. Joint density functional theory (DFT) and axial next nearest neighbor Ising (ANNNI) models were used to explain the temperature-dependent phase formation of the NiCrCo nanolaminate at the surface of the hard-turned material.

36 MATERIALS SCIENCE↗

DOE OSTI2015

WE-AB-BRA-02: Development of Biomechanical Models to Describe Dose-Volume Response to Liver Stereotactic Body Radiation Therapy (SBRT) Patients

Purpose: Previous studies have shown that radiotherapy treatment for liver metastases causes marked liver hypertrophy in areas receiving low dose and atrophy/fibrosis in areas receiving high dose. The purpose of this work is to develop and evaluate a biomechanical model-based dose-response model to describe these liver responses to SBRT. Methods: In this retrospective study, a biomechanical model-based deformable registration algorithm, Morfeus, was expanded to include dose-based boundary conditions. Liver and tumor volumes were contoured on the planning images and CT/MR images three months post-RT and converted to finite element models. A thermal expansion-based relationship correlating the delivered dose and volume response was generated from 22 patients previously treated. This coefficient, combined with the planned dose, was applied as an additional boundary condition to describe the volumetric response of the liver of an additional cohort of metastatic liver patients treated with SBRT. The accuracy of the model was evaluated based on overall volumetric liver comparisons and the target registration error (TRE) using the average deviations in positions of identified vascular bifurcations on each set of registered images, with a target accuracy of the 2.5mm isotropic dose grid (vector dimension 4.3mm). Results: The thermal expansion coefficient models the volumetric change of the liver to within 3%. The accuracy of Morfeus with dose-expansion boundary conditions a TRE of 5.7±2.8mm compared to 11.2±3.7mm using rigid registration and 8.9±0.28mm using Morfeus with only spatial boundary conditions. Conclusion: A biomechanical model has been developed to describe the volumetric and spatial response of the liver to SBRT. This work will enable the improvement of correlating functional imaging with delivered dose, the mapping of the delivered dose from one treatment onto the planning images for a subsequent treatment, and will further provide information to assist with the biological characterization of patients’ response to radiation.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Calibration of elastoplastic constitutive model parameters from full-field data with automatic differentiation-based sensitivities

Here, we present a framework for calibration of parameters in elastoplastic constitutive models that is based on the use of automatic differentiation (AD). The model calibration problem is posed as a partial differential equation-constrained optimization problem where a finite element (FE) model of the coupled equilibrium equation and constitutive model evolution equations serves as the constraint. The objective function quantifies the mismatch between the displacement predicted by the FE model and full-field digital image correlation data, and the optimization problem is solved using gradient-based optimization algorithms. Forward and adjoint sensitivities are used to compute the gradient at considerably less cost than its calculation from finite difference approximations. Through the use of AD, we need only to write the constraints in terms of AD objects, where all of the derivatives required for the forward and inverse problems are obtained by appropriately seeding and evaluating these quantities. We present three numerical examples that verify the correctness of the gradient, demonstrate the AD approach's parallel computation capabilities via application to a large-scale FE model, and highlight the formulation's ease of extensibility to other classes of constitutive models.

42 ENGINEERING↗

DOE OSTI2022

A Novel Anatomic Reconstruction for Posterior Tibialis Tendon in Treatment of Flexible Adult‐Acquired Flatfoot Deformity

Objective To present a novel approach for the anatomic reconstruction of the posterior tibialis tendon (PTT) in restoring plantar insertions and evaluate its efficiency in treating flexible adult‐acquired flatfoot deformity (AAFD) caused by PTT dysfunction. Methods For AAFD treatment, a novel PTT reconstruction method was presented. The current study involved 16 patients, including three men, and 13 women, from August 2017 to July 2019. The mean age was 43.2 ± 15.1 years (21–64 years). The innovative PTT repair method was used on all patients. The treatment involved performing a traditional Flexor Digitorum Longus (FDL) transfer in the navicular tuberosity and suturing the plantar insertions to FDL as tension was applied to tighten the plantar structures of the foot. The results were retrospectively analyzed. The clinical outcome was assessed using the pain visual analogue scale (VAS), the satisfaction VAS, and the American Orthopedic Foot and Ankle Society ankle–hindfoot scale (AOFAS‐AH). Isokinetic testing was performed using a dynamometer at 60°/s and 120°/s for inversion/eversion and plantarflexion/dorsiflexion, respectively, to determine the mean peak torque. Radiographic measurements were employed to assess the outcomes. Results Bone surgeries combined with the modified anatomic PTT reconstruction were performed on patients with medializing calcaneal osteotomy in 12 (75%) patients and subtalar joint fusion in four (25%) patients. The branch linking to the plantar insertions was detected in every case, with an average width of 3.5 ± 0.8 mm (3.1–4.3 mm). All patients were followed up for the mean of 16.8 ± 1.8 months (range, 15–20 months). The average postoperative functional scores, including pain VAS, satisfaction VAS, total AOFAS‐AH, and all AOFAS‐AH sub‐scales, steadily improved during the follow‐up. In the last follow‐up, isokinetic testing revealed no loss of plantarflexion strength ( p = 0.350 and 0.098) and significant improvement in the inversion strength ( p = 0.007 and 0.008) in the operated ankles at 60°/s and 120°/s. Radiographic outcomes, particularly the talar head uncovering, improved significantly after more than a year ( p < 0.001 for all). Conclusions The novel technique for PTT reconstruction in restoring the plantar insertions serves as an effective procedure in treating AAFD caused by PTT dysfunction in terms of delivering a consistent improvement in ankle inversion strength, medial longitudinal arch restoring, and satisfactory clinical outcomes.

Hu, Yong↗