Search NASA⌕ Search

SEARCH · Search NASA

Results for “material parameters”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neutron scattering off Weyl semimetals

We present how to detect type-I Weyl nodes in a material by inelastic neutron scattering. Such an experiment first of all allows one to determine the dispersion of the Weyl fermions. We extend the reasoning to produce a quantitative test of the Weyl equation, taking into account realistic anisotropic properties. These anisotropies are mostly contained in the form of the emergent magnetic moment of the excitations, which determines how they couple to the neutrons. Although there are many material parameters, we find several quantitative predictions that are universal and demonstrate that the excitations are described by solutions to the Weyl equation. The anisotropic coupling between electrons and neutrons implies that even fully unpolarized neutrons can reveal the spin-momentum locking of the Weyl fermions because the neutrons will couple to some components of the Weyl fermion pseudospin more strongly. On the other hand, in an experiment with polarized neutrons, the scattered neutron beam remains fully polarized in a direction that varies as a function of momentum transfer (within the range of validity of the Weyl equation). This allows measurement of the chirality of Weyl fermions for inversion-symmetric nodes. Furthermore, in this work we estimate that the scattering rate may be large enough for such experiments to be practical; in particular, the magnetic moment may be larger than the ordinary Bohr magneton, compensating for a small density of states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Smart material based multilayered microbeam structures for spatial self-deployment and reconfiguration: A residual stress approach

Alleviation of the potentially damaging effects induced by residual stresses was comprehensively investigated in previous research. Here, this paper, however, presents a spatially self-deployable and reconfigurable multilayered microbeam which takes advantage of residual stresses and shape memory effects. Reconfigurable mechanism of a typical four-layered microbeam composed of Pt\Ni 50 Ti 50 \Ni 50 Ti 50 \Pt is introduced, followed by analytical modeling of the maximum distance of the self-deployed gap as functions of variable structural and material parameters, including compressive residual stress in Pt layers and tensile residual stress in Ni 50 Ti 50 layers. Analytical solutions given by the static model agree well with the results obtained via finite element models (FEMs). Fabrication, characterization, and in-situ experiments were carried out to validate the feasibility of deployment of the as-released four-layered microbeam. The maximum distance of the gap was measured to be 41.39 μm at 20 °C, which could be increased to 51.73 μm thanks to controllable reconfiguration driven by shape memory effects. Theoretical analysis of such self-deployment and reconfiguration suggested a tensile residual stress increase by 52 MPa in Ni 50 Ti 50 layers. The multilayered microbeam structure with capabilities of self-deployment and reconfiguration offers great potential for various emerging applications, such as micro robotics, medical drug delivery devices, and intelligent chip scale spacecraft.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum efficiency of InGaN–GaN multi-quantum well solar cells: Experimental characterization and modeling

InGaN-based multi-quantum well (MQW) solar cells are promising devices for photovoltaics (e.g., for tandem solar cells and concentrator systems), space applications, and wireless power transfer. In order to improve the efficiency of these devices, the factors limiting their efficiency and stability must be investigated in detail. Due to the complexity of a MQW structure, compared with a simple pn junction, modeling the spectral response of these solar cells is not straightforward, and ad hoc methodologies must be implemented. In this paper, we propose a model, based on material parameters and closed-formula equations, that describes the shape of the quantum efficiency of InGaN/GaN MQW solar cells, by taking into account the layer thickness, the temperature dependence of the absorption coefficient, and quantum confinement effects. We demonstrate (i) that the proposed model can effectively reproduce the spectral response of the cells; in addition, (ii) we prove that the bulk p-GaN layer absorbs radiation, but the carriers photogenerated in this region do not significantly contribute to device current. Finally, we show that (iii) by increasing the temperature, there is a redshift of the absorption edge due to bandgap narrowing, which can be described by Varshni law and is taken into account by the model, and a lowering in the extraction efficiency due to the increase in recombination (mostly Shockley–Read–Hall) inside the quantum wells, which is also visible by decreasing light intensity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A molecular dynamics survey of grain boundary energy in uranium dioxide and cerium dioxide

Uranium dioxide (UO 2 ) is the primary fuel material that is used in current nuclear reactors. As one of the most fundamental material parameters, grain boundary (GB) energy strongly influences many fuel properties, and the influences depend on the characters and properties of individual GBs. Using molecular dynamics simulations, a high throughput survey of GB energy in UO 2 was carried out for the purpose of elucidating the roles of GB geometry such as misorientation and inclination, as well as the bonding nature of UO 2 , in affecting GB energy. GB energies in CeO 2 were calculated as well for comparison with UO2toinvestigate the generality of GB energy anisotropy in fluorite phase oxides. The results show significant GB energy anisotropy in both UO 2 and CeO 2 that is associated with the cubic symmetry of the fluorite structure. More interestingly, the GB anisotropy is found to be dependent not only on the crystal structure but also the ionic bonding. As such, the GB energy anisotropy in fluorite oxides has significant differences compared with that in fcc metals. We report that the data obtained and the increased knowledge on GB anisotropy will facilitate GB engineering for nuclear fuels with improved properties

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Durable and Low-Cost Fractal Structured Multifunctional Coatings for Next-Generation CSP

Fractal coatings on various low-cost and high-temperature metallic substrates such as SS316, SS347, In800H, In740H, In625, Haynes 230, etc. were fabricated employing our patented novel electrodeposition process at various deposition potentials and deposition times, as well as chemical etching. The best combination of the materials and processing conditions were co-optimized for the reliable and repeatable fabrication of the multifunctional coatings. Surface morphologies were examined using SEM and AFM, and fractal characterization was done using a 3D optical surface profiler. Fractal dimension was used to identify successful fabrication and material parameters. The mechanical durability of the coatings was evaluated using several standardized tests, such as sand abrasion (ASTM D968) and adhesion (ASTM D3359), to fully characterize and demonstrate the mechanical durability of the fabricated surfaces. HTF immersion tests on the fractal multifunctional coatings were conducted based on a modification of the ASTM D870, which is established for static water immersion testing. To calculate the corrosion rates and descaled weight loss of the samples subjected to thermal exposure for different time durations (100 h, 300 h, 500 h, and optionally 750 h) to estimate the corrosion rate. While protecting the substrates from corrosion is an important goal, an allied goals is also to examine the post-corrosion temperature stability and properties of the molten salts themselves due to the dissolution of the constituents of the alloys and the corrosion products from the coatings that cause “fouling” of the molten salts. Fouling characteristics of the different molten salt HTFs exposed to coated and uncoated substrates were investigated. Finally, considering a prototype molten salt/sCO2 heat exchanger, cost/performance model is developed to demonstrate that the coatings on low-cost alloys such as stainless steel provide lower levelized cost of heat exchanger compared to high nickel content Ha230.

14 SOLAR ENERGY↗
DOE OSTI2023

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. In this report we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO 2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revisiting the Role of Mobile Ion Size for the Activation Barrier of Charge Transport in Single Ion Conducting Polymers

The development of polymer electrolytes for energy storage and conversion technologies requires a fundamental understanding of the material parameters controlling the energy barriers for ion transport. In glassy polymers, these activation barriers are usually extracted by using Arrhenius procedures. However, our recent studies of single ion conducting polymers reveal that this traditional Arrhenius description provides anomalously small prefactors, an issue that is widely disregarded in the literature. Here, to address it, an alternative approach was introduced for extracting the effective activation barrier E by imposing a physically valid limit (∼10 –13 s) for the ion hopping attempt time. Two consequences of this recent approach are that in polymer electrolytes (i) E is significantly (∼30–40%) lower than the barrier estimated using traditional Arrhenius fits and (ii) E displays significant temperature dependence even in the glassy state. Under these circumstances, we are revisiting the role of ion size, glass transition temperature T g , and dielectric and elastic constants of the polymer matrix on the effective activation barriers for ion transport as extracted using the recently proposed approach and highlight the differences from the picture that previously emerged based on simplistic apparent Arrhenius analysis. To this end, we investigate cation transport in three families of single ion conducting trifluoromethane sulfonimide-based polyanions with varied T g . Our results indicate that E decreases with mobile cation size and is highly sensitive to changes in the dielectric permittivity of the matrix, even for large cations. These insights call for revisions of many earlier results based on apparent Arrhenius fits as the proposed approach can provide more accurate guidance for the design of polymer electrolytes with enhanced ionic conductivity.

energy barrier↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2021

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗
DOE OSTI2024

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗
DOE OSTI2020

An analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure

The objective of this paper is to develop an analytical expression for temperature in a thermodynamically consistent model with a Mie-Gruneisen equation for pressure. This expression for the temperature is obtained for both constant and variable specific heat models and can be used in common models for the material yield strength that include thermal softening due to increase in temperature. Additionally, the material parameters in the model can easily be calibrated to match known experimental data and the model has been used to determine the value of the shock pressure at the onset of melting during shock loading.

42 ENGINEERING↗
DOE OSTI2020

Microstructural characterization of dendritic evolution using two-point statistics

Characterizing large four-dimensional materials datasets is difficult due to the presence of complex mi- crostructures and time-varying length scales. We showcase the use of two-point statistics as an efficient and un-biased way of extracting materials parameters from an Al-Cu alloy during solidification. The evo- lution of dendrite primary arm thickness, average secondary arm spacing, and average tip-to-tip spacing were tracked using two-point Pearson auto-correlations of scaled mean curvatures. Insights into compet- itive side-branching are also reported. We show both visually and quantitatively that most length scales change rapidly during early stages of dendritic growth, but slow as diffusion fields of dendrites overlap.

Aluminum alloys↗
DOE OSTI2020

Modeling of thermal pressurization in tight claystone using sequential THM coupling: Benchmarking and validation against in-situ heating experiments in COx claystone

We apply thermoporoelasticity and a sequentially coupling technique for modeling thermally-driven coupled Thermo-Hydro-Mechanical (THM) processes in tight claystone. A THM benchmark case with a corresponding analytic solution for thermoporoelasticity under a constant heat loading verifies the model. Thereafter, two in situ heating experiments are simulated for model validation: a smaller-scale heating experiment (TED experiment) and a larger-scale experiment (ALC experiment) in Callovo-Oxfordian (COx) claystone at the Meuse/Haute-Marne underground research laboratory in France. The model exhibits good performance to match the observed temperature and pore pressure evolution for the smaller-scale TED experiment. For the larger-scale ALC experiment, general trends of thermal-pressurization are captured in the modeling, but pressure is underestimated at some monitoring points during cool-down. This indicates that the THM response in the field may be affected by the variability of rock's properties or irreversible or time-dependent mechanical processes that are not included in the current thermoporoelastic model. The main contributions of this work are as follows: (1) we verify and validate the numerical simulator, TOUGH-FLAC, to be a valuable coupled THM modeling tool; (2) prove that the laboratory determined material parameters can be used as reference values for upscaling experiments. However, to better identify and quantify THM processes with modeling of in situ tests, more emphasize should be dedicated to obtaining high-quality mechanical deformation data.

58 GEOSCIENCES↗