Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariant”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

An unsupervised machine-learning checkpoint-restart algorithm using Gaussian mixtures for particle-in-cell simulations

We propose an unsupervised machine-learning checkpoint-restart (CR) algorithm for particle-in-cell (PIC) algorithms using Gaussian mixtures (GM). The algorithm compresses the particle population per spatial cell by constructing a velocity distribution function using GM. Particles are reconstructed at restart time by local resampling of the Gaussians. To guarantee fidelity of the CR process, we ensure the exact preservation of invariants such as charge, momentum, and energy for both compression and reconstruction stages, everywhere on the mesh. We also ensure the preservation of Gauss' law after particle reconstruction by exactly matching the density profile at restart time. As a result, the GM CR algorithm is shown to provide a clean, conservative restart capability while potentially affording orders of magnitude savings in input/output requirements. Here, we demonstrate the algorithm using a recently developed exactly energy- and charge-conserving PIC algorithm using both electrostatic and electromagnetic tests. The tests demonstrate not only a high-fidelity CR capability, but also its potential for enhancing the fidelity of the PIC solution for a given particle resolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Scalable semidefinite programming approach to variational embedding for quantum many-body problems

In quantum embedding theories, a quantum many-body system is divided into localized clusters of sites which are treated with an accurate ‘high-level’ theory and glued together self-consistently by a less accurate ‘low-level’ theory at the global scale. The recently introduced variational embedding approach for quantum many-body problems combines the insights of semidefinite relaxation and quantum embedding theory to provide a lower bound on the ground-state energy that improves as the cluster size is increased. The variational embedding method is formulated as a semidefinite program (SDP), which can suffer from poor computational scaling when treated with black-box solvers. Here, we exploit the interpretation of this SDP as an embedding method to develop an algorithm which alternates parallelizable local updates of the high-level quantities with updates that enforce the low-level global constraints. Moreover, we show how translation invariance in lattice systems can be exploited to reduce the complexity of projecting a key matrix to the positive semidefinite cone.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Permutation-adapted complete and independent basis for atomic cluster expansion descriptors

Atomic cluster expansion (ACE) methods provide a systematic way to describe particle local environments of arbitrary body order. For practical applications it is often required that the basis of cluster functions be symmetrized with respect to rotations and permutations. Existing methodologies yield sets of symmetrized functions that are over-complete. These methodologies thus require an additional numerical procedure, such as singular value decomposition (SVD), to eliminate redundant functions. In this work, it is shown that analytical linear relationships for subsets of cluster functions may be derived using recursion and permutation properties of generalized Wigner symbols. From these relationships, subsets (blocks) of cluster functions can be selected such that, within each block, functions are guaranteed to be linearly independent. It is conjectured that this block-wise independent set of permutation-adapted rotation and permutation invariant (PA-RPI) functions forms a complete, independent basis for ACE. Along with the first analytical proofs of block-wise linear dependence of ACE cluster functions and other theoretical arguments, numerical results are offered to demonstrate this. The utility of the method is demonstrated in the development of an ACE interatomic potential for tantalum. Using the new basis functions in combination with Bayesian compressive sensing sparse regression, some high degree descriptors are observed to persist and help achieve high-accuracy models.

Angular momentum↗
DOE OSTI2026

Quantum mechanical closure of partial differential equations with symmetries

We develop a statistical framework for the dynamical closure of spatiotemporal dynamics governed by partial differential equations. Employing the mathematical framework of quantum mechanics to embed the original classical dynamics into a quantum mechanical representation, we use the space of quantum density operators to model the unresolved degrees of freedom of the original dynamics in a statistical sense, and the framework of quantum measurement to predict their contributions to the resolved dynamics. The embedded dynamics is discretized by a positivity preserving process, leading to a compressed representation that is invariant under the dynamical symmetries of the resolved dynamics. We present a data based formulation of the closure scheme and apply it to a closure problem for the shallow water equations. The numerical results demonstrate that our closure model can accurately predict the main features of the true dynamics, including for out of sample initial conditions.

Delay embedding↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
DOE OSTI2022

Structural investigation of the temperature-stable relaxor dielectric Ba 0.8 Ca 0.2 TiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3

Aberration corrected scanning transmission electron microscopy (STEM) and electron diffraction have been used to disclose local structure and nano-chemistry in a Ca modified BaTiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3 relax or dielectric ceramic which exhibits high and near-invariant relative permittivity over a wide temperature range. High resolution, synchrotron X-ray diffraction indicated a globally cubic structure ($Pm$$\overline{3}$$m$), but direct atomic-scale imaging by STEM revealed local tetragonal distortions. Nanopolar clusters were identified from B-site atomic displacement vectors measured relative to oxygen ion positions along < 100 > and < 110 > zone axes of integrated differential phase contrast (iDPC) STEM images, highlighting cluster sizes of 2–5 nm. Chemical analysis by STEM-energy dispersive X-ray spectroscopy and full pattern refinements of X-ray powder diffraction data each implied high levels of Bi vacancies within the matrix. The possibility that these A-site vacancies modulate the nanopolar structure and promote flattening of the permittivity-temperature response in this class of dielectric is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling the processes of soil moisture in regulating microbial and carbon-nitrogen cycling

Soil carbon (C) and nitrogen (N) cycles and their complex responses to hydro-climatic forcing have gained increasing attention. While the temperature effects have been intensively studied, soil moisture response functions (SMRFs) are not well documented for various microbial and enzymatic processes due to the difficulties in directly measuring and differentiating the moisture effects on various processes. In this work we extended our C-only Microbial-ENzyme Decomposition (MEND) model to the C-N coupled MEND model with flexible element stoichiometry. Our model calibration showed good agreement between simulated and observed C:N ratios in soil organic matter and microbial biomass, as well as the ammonium and nitrate concentrations. We show that the selection of SMRFs for specific biogeochemical processes could result in significant differences in model simulated microbial and C-N processes. In particular, it is essential to account for the soil moisture effects on microbial dormancy and resuscitation, as the changes in microbial physiology under favorable or stressful conditions will exert strong controls on soil C and N dynamics. We also advocate the utilization of dynamic (time-variant) data (e.g., litter input, N deposition, soil temperature and moisture), instead of time-invariant data, to drive model simulations and analyses. Dynamic forcing data (particularly dynamic soil moisture) better represent the real-world climate and environmental conditions, which could facilitate more realistic modeling and understanding of soil C and nutrient cycling in a changing world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Multi-fidelity Fourier neural operator for fast modeling of large-scale geological carbon storage

Deep learning-based surrogate models have been widely applied in geological carbon storage (GCS) problems to accelerate the prediction of reservoir pressure and CO2 plume migration. Large amounts of data from physics-based numerical simulators are required to train a model to accurately predict the complex physical behaviors associated with this process. In practice, the available training data are always limited in large-scale 3D problems due to the high computational cost. Therefore, we propose to use a multi-fidelity Fourier neural operator (FNO) to solve large-scale GCS problems with more affordable multi-fidelity training datasets. FNO has a desirable grid-invariant property, which simplifies the transfer learning procedure between datasets with different discretization. Here, we first test the model efficacy on a GCS reservoir model being discretized into 110 k grid cells. The multi-fidelity model can predict with accuracy comparable to a high-fidelity model trained with the same amount of high-fidelity data with 81% less data generation costs. We further test the generalizability of the multi-fidelity model on a same reservoir model with a finer discretization of 1 million grid cells. This case was made more challenging by employing high-fidelity and low-fidelity datasets generated by different geostatistical models and reservoir simulators. We observe that the multi-fidelity FNO model can predict pressure fields with reasonable accuracy even when the high-fidelity data are extremely limited. The findings of this study can help for better understanding of the transferability of multi-fidelity deep learning surrogate models.

58 GEOSCIENCES↗
DOE OSTI2021

Fatigue of laser powder bed fusion processed 17-4 stainless steel using prior process exposed powder feedstock

The rapid pace of development seen in the metal additive manufacturing (AM) process of laser powder-bed fusion (LPBF) requires in-step advances in processes qualification to enable full-scale adoption. This particularly applies to quantifying how powder feedstock conditions impact end-component quality. Here this study examines how in-machine 17-4 stainless steel powder feedstocks are affected by prior LPBF processes, and how these effects impact subsequent builds. Examinations of powder morphology, chemistry, flowability, and rheology were conducted to characterize the powder conditions. The resultant effects of powder feedstock condition on produced component quasi-static tensile and high-cycle fatigue properties were analyzed. Fatigue life was analyzed using a reliability modeling approach that enabled a robust statistical comparison of life. Powder characteristics were found to evolve with powder exposure to prior LPBF processes, particularly in the extremes of powder size distribution and measures of bulk flow. No significant effects of these changes on tensile properties were observed. Reliability modeling methods, including the lognormal and Weibull distributions as well as the empirical survival function, are shown to be effective tools for modeling fatigue variability in LPBF manufactured components. Through these tools, fatigue life was found to be invariant with changes in powder condition.

42 ENGINEERING↗
DOE OSTI2024

Spherical Indentation and Implementation of S 3 /P for yield stress determination of brittle materials

Here, a mathematically transparent and robust experimental method has been developed to estimate the yield stress of brittle materials through the analysis of depth-sensing spherical indentation. Employing Hertzian contact mechanics, an elastically invariant ratio based on the simple equation S 3 / P = 6RE r 2 , (where S and P are contact stiffness and indentation load, respectively) has been derived that enables more accurate and confident determination of the transition from elastic to inelastic deformation; a transition that the yield stress dictates and represents. Using two diamond spheres with radii of 3.2 and 8.6 μm, the indentation test method and analyses are applied to two vitreous silicates: Corning's HPFS 7980® fused silica and Vitro's Starphire® soda lime silicate. The estimated yield strengths are 8.15 GPa ± 2.5% for the fused silica and 6.1 GPa ± 3.3% for the soda-lime silicate, and both were independent of indenter radius. Verification of this new experimental method is demonstrated with an as-drawn titanium by showing equivalence of measured yield stress by its spherical indentation and that from uniaxial compression testing. This method will enable easier and more confident estimation of yield stress in brittle materials - a property that historically has been elusive to measure for these materials using common laboratory mechanical test methods.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Polyconvex neural network models of thermoelasticity

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. Here, in this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

In situ measurement of three-dimensional intergranular stress localizations and grain yielding under elastoplastic axial-torsional loading

The three-dimensional grain-averaged response of solid bar samples under non-proportional (NP) elastoplastic axial-torsional loading was investigated using in situ high energy diffraction microscopy (HEDM) and companion crystal plasticity finite element (CPFE) modeling. Important stress metrics including applied shear (σ θZ ) and axial (σ ZZ ) stress tensor components, stress and stress deviator tensor invariants (I 1 , J 2 , and J 3 ), von Mises equivalent stress (σ$^{grain}_{VM}$), maximum resolved shear stress (mRSS), stress triaxiality (η), and lode angle parameter ($\barθ$) values were tracked for ~300 grains under two different loading conditions: (1) Torsion-dominated loading (low NP) and (2) Tension-torsion loading (high NP) in equiatomic NiCoCr, a representative multicomponent face-centered cubic (FCC) superalloy. Overall, significant stress localizations existed within both samples as evidenced by the radial dependence of grain-resolved σ θZ , σ$^{grain}_{VM}$, and J 2 ; by comparison, I 1 , J 3 , η, and $\barθ$ metrics did not show discernible trends within the volume. These stress localizations reveal a complex interplay between axial and shear stress components (e.g., stress coupling) resulting in grain yielding near the sample surface largely driven by shear stress, whereas internal grain yielding was largely accommodated by axial stress. Grain-resolved stress localization trends were described well by the CPFE model, although some discrepancies in magnitude occurred, particularly for volumetric stress metrics (I 1 and η) due to initial type II residual stress distributions. The superposition of initial residual stress states onto CPFE grain-resolved data significantly improved model accuracy for η. This suggests that residual stresses more strongly influence the simulation of volumetric rather than deviatoric (yield) stress metrics.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Measurement of PuO 2 film thickness by electron probe microanalysis (EPMA) calibration curve method

Plutonium (Pu) surfaces are highly reactive toward oxygen containing species and, therefore, invariably covered with oxides (e.g. PuO 2 ) during transport and handling. The actual thickness of the surface oxide may dictate if a plutonium part is suitable for a certain application. As a result, a cost-effective, quick, non-destructive, yet reliable means to measure the oxide layer thickness formed on Pu samples is desirable. In this study, the cross-sections of a series of room temperature grown oxides on Pu samples were trenched by focus ion beam (FIB) then observed by scanning electron microscopy (SEM) to measure the surface oxide thicknesses, which were then combined with the corresponding oxygen k-ratios provided by electron probe microanalysis (EPMA) to form calibration curves. Oxide thickness measurements for the calibration curves were made on samples within the typical SEM observable range for PuO 2 (35–400 nm). The portion of the calibration curve in the thinner oxide region (<35 nm) were approximated via Pouchou and Pichoir’s Φ(ρz) theory. Two specimens with micrometer-thick PuO 2 standards (one formed at room temperature and the other at higher temperature with a higher level of crystallinity) were made for the k-ratios in this study, allowing EPMA users to choose the standard that best suits their needs. If the surface corrosion is known to be PuO 2 (from the environment in which the Pu sample is stored) or if the stoichiometry of the surface oxide is confirmed by a preliminary/compliment technique, these calibration curves allow EPMA users to quickly and efficiently determine PuO 2 thicknesses from the measured oxygen k-ratios of their samples. The methodology presented in this study can also be used as a template for creating calibration curves for oxides grown on other actinides.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In situ high-energy X-ray study of deformation mechanisms in additively manufactured 316L stainless steel

A key structural material planned to be used in the core of the Transformational Challenge Reactor is the additively manufactured 316L stainless steel (AM 316L SS). While the deformation of conventionally-manufactured 316 SS have been extensively studied, the deformation mechanism of AM 316L SS is less known. Here the room-temperature tensile behavior of AM 316L SS was studied in situ using high-energy X-ray diffraction and X-ray tomography. Diffraction yielded the lattice strains and the dislocation densi-ties of the AM 316L SS, which were compared to those of a conventionally processed 316H SS. It was found that while the conventional 316H SS had three distinct deformation stages, namely stage II ather-mal hardening, stage III dynamic recovery, and stage IV, the AM 316L experienced only the stages III and IV after an initial transitional period. The plastic instability stresses (PIS), i.e. the true ultimate tensile strength, were nearly the same in both materials. The AM 316L had very low dislocation barrier strength throughout the plastic deformation, which, in combination with the invariant PIS, explained its high dis-location storage capability and the attractive combination of strength and ductility. Tomography revealed morphological changes of the built-in pores in the AM 316L during deformation and fracture, and it was concluded that large pores close to the surface might have significant roles in the necking stage.

316 stainless steel↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2025

Supported ionic liquid membranes (SILMs) with exceptional selectivity and permeability for dilute CO 2 separations

Supported ionic liquid membranes (SILMs), containing phosphonium ionic liquids with aprotic N-heterocyclic anions (AHA ILs) in an inert inorganic support, were tested under both dry and humidified (40 % RH) mixed-gas conditions down to 420 ppm CO 2 in N 2 at 35 °C. In the dry case, the best performing IL, triethyl(octyl)phosphonium 4-bromopyrazolide ([P 2228 ][4-BrPyra]) exhibited mixed-gas CO 2 permeabilities and CO 2 /N 2 permeability selectivities as high as 26,800 barrer and 7,000, respectively. In the presence of humidity, the CO 2 permeability and CO 2 /N 2 selectivity increased to 49,100 barrer and 13,200, respectively, and these are the highest reported combination in the literature. Humidity amplifies CO 2 permeabilities and CO 2 /N 2 permeability selectivities through increases in CO 2 capacity due to bicarbonate formation and through faster mobility of the mobile carrier from decreased viscosity. Furthermore, N 2 permeability stayed roughly invariant in the presence of humidity, likely from competing effects of viscosity reduction and lower N 2 solubility.

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

Effect of copper content on the tensile elongation of Al–Cu–Mn–Zr alloys: Experiments and finite element simulations

Microstructures of cast aluminum alloys used in automotive engine applications often consist of intermetallic particles that can impact the tensile elongation of these alloys. Here in this paper, we investigate the effect of intermetallic grain boundary particles on tensile elongation by fabricating a series of alloys with Cu content varying between 6.0 - 9.0 wt% in cast Al–Cu–Mn–Zr (ACMZ) type compositions. The tensile elongation of as-aged ACMZ alloys decreases monotonically with increase in Cu content. While the microstructure within the grains and yield stress of the alloy remains invariant with Cu content, the decrease in tensile elongation correlates well with increase in the size and volume fraction of grain boundary particles. Microstructural observations are combined with finite element simulations to explain the trend in tensile elongation with changing Cu content. Crack initiation is found to occur by brittle fracture of the grain boundary particles. Increase in particle size promotes crack initiation by reduction in size dependent particle fracture strength. Lower inter-particle spacing at higher particle volume fraction further facilitates crack initiation by increasing stress within the particles caused by the interaction between stress fields of neighboring particles. Increase in particle volume fraction also accelerates crack propagation through the formation of macro shear zones in the microstructure. The increase in Cu content of cast ACMZ alloys, therefore, decreases tensile elongation by promoting both crack initiation and crack propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Bayesian automated weighting of aggregated DFT, MD, and experimental data for candidate thermodynamic models of aluminum with uncertainty quantification

Atomic-scale modeling methods such as density functional theory (DFT) and molecular dynamics (MD) can predict the thermodynamic properties of materials at a lower cost than experimental measurements. However, their regular usage in thermodynamic model construction is hampered by the lack of quantitative agreement with experimental measurements and the lack of uncertainty estimates on the data. To make regular usage of this atomistic simulation data, it is important to assess whether the atomistic simulation datasets, by themselves or in combination with experimental measurements, result in the same physics-informed models best supported by experimental measurements alone. Here, models of aluminum thermodynamic properties are discussed using three data sources: atomistic calculations (DFT and MD), experiments, and a combination of atomistic calculations and experiments. The study shows that, after ensuring self-consistency in predicting key invariant points, both experimental measurements and atomistic calculations can significantly contribute to an optimal model.

36 MATERIALS SCIENCE↗