Search NASA⌕ Search

SEARCH · Search NASA

Results for “extrapolation method”

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 289 records · Page 16

Monotonic Gaussian Process for Physics-Constrained Machine Learning With Materials Science Applications

Physics-constrained machine learning is emerging as an important topic in the field of machine learning for physics. One of the most significant advantages of incorporating physics constraints into machine learning methods is that the resulting model requires significantly less data to train. By incorporating physical rules into the machine learning formulation itself, the predictions are expected to be physically plausible. Gaussian process (GP) is perhaps one of the most common methods in machine learning for small datasets. In this paper, we investigate the possibility of constraining a GP formulation with monotonicity on three different material datasets, where one experimental and two computational datasets are used. The monotonic GP is compared against the regular GP, where a significant reduction in the posterior variance is observed. The monotonic GP is strictly monotonic in the interpolation regime, but in the extrapolation regime, the monotonic effect starts fading away as one goes beyond the training dataset. Imposing monotonicity on the GP comes at a small accuracy cost, compared to the regular GP. The monotonic GP is perhaps most useful in applications where data are scarce and noisy, and monotonicity is supported by strong physical evidence.

36 MATERIALS SCIENCE↗

Synergy of semiempirical models and machine learning in computational chemistry

Catalyzed by enormous success in the industrial sector, many research programs have been exploring data-driven, machine learning approaches. Performance can be poor when the model is extrapolated to new regions of chemical space, e.g., new bonding types, new many-body interactions. Another important limitation is the spatial locality assumption in model architecture, and this limitation cannot be overcome with larger or more diverse datasets. The outlined challenges are primarily associated with the lack of electronic structure information in surrogate models such as interatomic potentials. Given the fast development of machine learning and computational chemistry methods, we expect some limitations of surrogate models to be addressed in the near future; nevertheless spatial locality assumption will likely remain a limiting factor for their transferability. Here, we suggest focusing on an equally important effort—design of physics-informed models that leverage the domain knowledge and employ machine learning only as a corrective tool. In the context of material science, we will focus on semi-empirical quantum mechanics, using machine learning to predict corrections to the reduced-order Hamiltonian model parameters. The resulting models are broadly applicable, retain the speed of semiempirical chemistry, and frequently achieve accuracy on par with much more expensive ab initio calculations. These early results indicate that future work, in which machine learning and quantum chemistry methods are developed jointly, may provide the best of all worlds for chemistry applications that demand both high accuracy and high numerical efficiency.

36 MATERIALS SCIENCE↗

Validation of Fault Displacements from Dynamic Rupture Simulations against the Observations from the 1992 Landers Earthquake

ABSTRACT Coseismic fault displacements in large earthquakes have caused significant damage to structures and lifelines on and near fault lines. Coseismic displacements represent a real threat, especially to distributed infrastructure systems. For infrastructure systems that can not avoid active faults, engineering displacement demands are defined using probabilistic fault-displacement hazard analyses (PFDHA). However, PFDHA models are sparse and poorly constrained partly due to the scarcity of detailed fault-displacement observations. Advancements in dynamic rupture simulation methods make them an attractive approach to address this important issue. Because fault displacements can be simulated for various geologic conditions as constrained by current knowledge about earthquake processes, they can be used to supplement the observation datasets. In addition to providing on-fault displacements, when used with appropriate constitutive models for the bulk medium, they can capture off-fault distributed inelastic deformations as well. For viable extrapolation, simulations must first be validated against data. In this article, we summarize the calibration and validation of the dynamic rupture model against the observations of the well-documented 1992 Landers earthquake. We defined a preferred model that reproduces several first-order fault-displacement metrics such as the on-fault partition of the total displacement, the mean fault-zone width, and the location of the peak displacement. Simulated ground motions consistent with the observations ensure that all physics important to modeling have been properly parameterized. For the extrapolation, we generated a suite of dynamic rupture models to quantify expected fault-displacement metrics, their intercorrelations, and magnitude dependencies, which are in part supported by the Landers and other recent earthquakes. Our validation and extrapolation exercise paves the way for using dynamic rupture modeling to quantitatively address fault-displacement hazard on a broader scale. The results are promising and are expected to be useful to inform PFDHA model development.

Geochemistry & Geophysics↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗

Inferring Heat Flow in Laser Absorption Regions using Diagnostic Magnetic Fields

This research was motivated by the desire to obtain a measurement of heat-flow in overdense regions of a laser ablated plasma for the first time using a high energy proton beam. For sufficiently high energies, the proton beam can pass through overdense regions and only be deflected by electric and magnetic fields. As the transport of magnetic fields in a plasma are intrinsically linked to the transport of thermal energy, it was posited that such a measurement could be used to infer a heat-flow. While the project did not prove the method to be infeasible, difficulties with understanding the proton radiographs without an imposed magnetic field made the final aim of this project unreached. Instead, however, a closer understanding of proton radiographs of laser-solid interactions without an imposed field has been obtained. 5 shot days on The OMEGA Laser Facility have been awarded (3 through ICF and 2 as ride-along) to further investigate some of these features. Extrapolating current results to hohlraum conditions has the magnetic energy being a significant (10%) fraction of the laser drive energy, which has the potential to explain the NIF hohlraum drive deficit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Inferring Heat Flow in Laser Absorption Regions using Diagnostic Magnetic Fields: Abstract

This research was motivated by the desire to obtain a measurement of heat-flow in overdense regions of a laser ablated plasma for the first time using a high energy proton beam. For sufficiently high energies, the proton beam can pass through overdense regions and only be deflected by electric and magnetic fields. As the transport of magnetic fields in a plasma are intrinsically linked to the transport of thermal energy, it was posited that such a measurement could be used to infer a heat-flow. While the project did not prove the method to be infeasible, difficulties with understanding the proton radiographs without an imposed magnetic field made the final aim of this project unreached. Instead, however, a closer understanding of proton radiographs of laser-solid interactions without an imposed field has been obtained. 5 shot days on The OMEGA Laser Facility have been awarded (3 through ICF and 2 as ride-along) to further investigate some of these features. Extrapolating current results to hohlraum conditions has the magnetic energy being a significant (10%) fraction of the laser drive energy, which has the potential to explain the NIF hohlraum drive deficit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collaborative Research: Properties and Dynamics of the Shallow Crust (Final Report)

Ground motions recorded at one location are often extrapolated to nearby regions within a given radius. Here, we explore the appropriateness of spatial extrapolation using data from seven small aperture seismic network deployments in southern California. Six of these deployments are linear arrays of 4-13 stations, and one is a 2D array of 13 stations at Pinyon Flats Observatory. The spatial footprint array diameters are 3 km or less, and each array was operational for a year or more. From our base catalog (M2.5+ earthquakes; 4038 events; September 2010 - June 2023), automated methods remove temporally overprinted waveforms from nearby events (< 5 km) in quick succession (< 5 min) and data with nonviable waveforms. These 200 samples per second data are filtered at 0.5-25 Hz and must have signal-tonoise ratios (SNRs) of 2.5+. Peak ground acceleration (PGA) and peak ground velocity (PGV) are derived individually from the maximum absolute values of each of the 3-component waveforms (vertical, northsouth, and east-west). Five of the seven arrays traverse the San Jacinto fault, and two do not. Ground motion observations are compared with theoretical estimates from Abrahamson et al., 2014. On average, arrays deployed within and across fault zones consistently record ground motions above theoretical expectations, whereas off-fault arrays record ground motions at or slightly below theoretical expectations. We attribute these differences to site conditions because these trends prevail for the full data suites. For each network and each individual channel, the coefficient of variation indicates that the standard deviations are ~30±6% of the mean. Exploring relative ground motion contributions from all three channels (ternary plots), as expected, most data show that vertical ground motions are attenuated compared to horizontal ones. However, this is not always the case for RA array data, where vertical motions can be ~2-3 times larger than horizontal motions for select events near Cahuilla, CA. These anomalously high vertical motions are focal mechanism-related. These results suggest that ternary plots created using only a small amount of data can be used as a data quality metric and a tool to find anomalous features in three-component data.

58 GEOSCIENCES↗

Fast & accurate emulation of two-body scattering observables without wave functions

We combine Newton's variational method with ideas from eigenvector continuation to construct a fast & accurate emulator for two-body scattering observables. The emulator will facilitate the application of rigorous statistical methods for interactions that depend smoothly on a set of free parameters. Our approach begins with a trial K or T matrix constructed from a small number of exact solutions to the Lippmann–Schwinger equation. Subsequent emulation only requires operations on small matrices. We provide several applications to short-range potentials with and without the Coulomb interaction and partial-wave coupling. It is shown that the emulator can accurately extrapolate far from the support of the training data. When used to emulate the neutron-proton cross section with a modern chiral interaction as a function of 26 free parameters, it reproduces the exact calculation with negligible error and provides an over 300x improvement in CPU time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Identification of Life Models for Li-Ion Batteries Using Penalized Regression and Bilevel Optimization

Reduced-order physics-based life models are extremely useful for rapidly predicting battery state-of-health and for simulating battery lifetime in arbitrary aging conditions. However, identification of well-parameterized models is difficult. This is because, for maximum usefulness in predicting lifetime under a variety of conditions, aging test data exhibits many degradation mechanisms, which all need to be accurately modeled. However, because aging tests are time-consuming and expensive, especially for large-format batteries, a minimum of tests are conducted while probing many stress factors. Building a well-parameterized model is then very challenging: an under-parameterized model will neglect critical degradation modes, and an over-parameterized model will extrapolate poorly to new testing conditions. To complicate this matter, the functional form of the model for any individual degradation rate can be very difficult to identify. In this work, the statistical tools of penalized regression and bilevel optimization are used to help identify both the functional forms of and optimize the parameters of reduced-order life models, accelerating identification of robust models. Model robustness is demonstrated through traditional statistics methods of cross-validation and sensitivity analysis, uncertainty quantification through bootstrap resampling and Monte-Carlo simulation, and simulation of real-world use cases.

47 OTHER INSTRUMENTATION↗

Hydroxyl radical yields in the heavy ion radiolysis of water

Scavenging experiments coupled with Monte Carlo track structure and kinetics simulations have been used to determining and understand the yields of hydroxyl radicals in the radiolysis of water with protons and carbon ions. Carbon dioxide production from aerated formic acid solutions with concentrations ranging from 10-3 to 1 mol./dm3 was used as a probe of the time dependence of the hydroxyl radical yield from about 5 ns to 1 ?s. Numerical extrapolation of the results at slightly different ion energies in combination with data from gamma radiolysis experiments allows a systematic analysis of both track average and track segment yields, with considerable track chemistry found to occur on the nanosecond to microsecond time scale. Track chemistry calculations utilizing Monte Carlo track structure simulations and nonhomogeneous independent reaction times diffusion-reaction modeling of product yields reproduce experimental observations. Subsequently, the track simulation methods were used to predict hydroxyl radical kinetics in deaerated and aerated water at neutral pH. Taken with literature data for helium ion beam studies, the results offer insight into the potential for indirect radiation induced damage to biological systems for light ions with very different track structures.

38 - RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCL↗

Model fusion with physics-guided machine learning: Projection-based reduced-order modeling

The unprecedented amount of data generated from experiments, field observations, and large-scale numerical simulations at a wide range of spatiotemporal scales has enabled the rapid advancement of data-driven and especially deep learning models in the field of fluid mechanics. Although these methods are proven successful for many applications, there is a grand challenge of improving their generalizability. This is particularly essential when data-driven models are employed within outer-loop applications like optimization. In this work, we put forth a physics-guided machine learning (PGML) framework that leverages the interpretable physics-based model with a deep learning model. Leveraging a concatenated neural network design from multi-modal data sources, the PGML framework is capable of enhancing the generalizability of data-driven models and effectively protects against or inform about the inaccurate predictions resulting from extrapolation. We apply the PGML framework as a novel model fusion approach combining the physics-based Galerkin projection model and long- to short-term memory (LSTM) network for parametric model order reduction of fluid flows. We demonstrate the improved generalizability of the PGML framework against a purely data-driven approach through the injection of physics features into intermediate LSTM layers. Our quantitative analysis shows that the overall model uncertainty can be reduced through the PGML approach, especially for test data coming from a distribution different than the training data. Moreover, we demonstrate that our approach can be used as an inverse diagnostic tool providing a confidence score associated with models and observations. The proposed framework also allows for multi-fidelity computing by making use of low-fidelity models in the online deployment of quantified data-driven models.

42 ENGINEERING↗

Finite-volume pionless effective field theory for few-nucleon systems with differentiable programming

Finite-volume pionless effective field theory provides an efficient framework for the extrapolation of nuclear spectra and matrix elements calculated at finite volume in lattice QCD to infinite volume, and to nuclei with larger atomic number. In this work, it is demonstrated how this framework may be implemented via a set of correlated Gaussian wave functions optimized using differentiable programming and via solution of a generalized eigenvalue problem. This approach is shown to be significantly more efficient than a stochastic implementation of the variational method based on the same form of correlated Gaussian wave functions, yielding comparably accurate representations of the ground-state wave functions with an order of magnitude fewer terms. The efficiency of representation allows such calculations to be extended to larger systems than in previous work. Further, the method is demonstrated through calculations of the binding energies of nuclei with atomic number A ϵ {2,3,4} in finite volume, matched to lattice QCD calculations at quark masses corresponding to m π = 806 MeV, and infinite-volume effective field theory calculations of A ϵ {2,3,4,5,6} systems based on this matching.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

On the applicability of the Redlich-Kister framework for viscosity estimation of molten halide salt mixtures

For molten halide salt mixtures already being utilized or under consideration for carbon-free energy production systems, it is crucial that their viscosity is well understood so that system thermal hydraulics can be reliably assessed. Because of the difficulty in accurately measuring molten halide viscosity and the sheer size of the matrix of possible higher order salt mixtures that may be of interest to the energy industry, there are several gaps in the quantified understanding of molten halide viscosity across this matrix. As such, both first-principles and semi-empirical modeling techniques may be crucial for rapidly assessing this broad, complex compositional domain. Herein, the Redlich-Kister framework is applied to assess the feasibility of broadly interpolating and estimating the viscosity of several pseudobinary and pseudoternary molten halide salt systems that may be of key interest to the energy industry. The framework is based on the assumption that an ideal component and a nonideal component collectively describe the viscosity as a function of composition and temperature for a given molten halide system. Three different ideal models were considered for the ideal component, including Grunburg-Nissan, Katti-Chaudhri, and Gambill methods. Regarding the pseudobinary interpolations, the Redlich-Kister models with either the Grunburg-Nissan or Katti-Chaudhri models as the ideal component resulted in either highly (average error less than 5%) or reasonably (average error less than 15%) accurate interpolations of pseudobinary halide viscosity; BeF 2 - or UF 4 -bearing salts tended to result in reasonably accurate interpolations, whereas other pseudobinary mixtures tended to show high accuracy. Regarding the pseudoternary extrapolations, the Redlich-Kister framework shows reasonable success at estimating the extent to which a pseudoternary system may indicate deviations from ideal Grunburg-Nissan mixing, where discrepancies with comparative experimental data generally stay within 30%. Finally, the primary reasons identified for such discrepancies are (1) inaccuracy in the underlying experimental data, (2) different complexation behavior in the higher order systems compared to the pseudobinary subsystems, and (3) extrapolation into temperatures too far out of the domain, which is valid for the underlying experimental data feeding the Redlich-Kister model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear masses learned from a probabilistic neural network

Machine learning methods and uncertainty quantification have been gaining interest throughout the last several years in low-energy nuclear physics. In particular, Gaussian processes and Bayesian neural networks have increasingly been applied to improve mass model predictions while providing well-quantified uncertainties. In this work, we use the probabilistic Mixture Density Network (MDN) to directly predict the mass excess of the 2016 Atomic Mass Evaluation within the range of measured data, and we extrapolate the inferred models beyond available experimental data. The MDN provides not only mean values but also full posterior distributions both within the training set and extrapolated testing set. We show that the addition of physical information to the feature space increases the accuracy of the match to the training data as well as provides for more physically meaningful extrapolations beyond the the limits of experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗

Fast and converged classical simulations of evidence for the utility of quantum computing before fault tolerance

A recent quantum simulation of observables of the kicked Ising model on 127 qubits implemented circuits that exceed the capabilities of exact classical simulation. We show that several approximate classical methods, based on sparse Pauli dynamics and tensor network algorithms, can simulate these observables orders of magnitude faster than the quantum experiment and can also be systematically converged beyond the experimental accuracy. Our most accurate technique combines a mixed Schrödinger and Heisenberg tensor network representation with the Bethe free entropy relation of belief propagation to compute expectation values with an effective wave function–operator sandwich bond dimension >16,000,000, achieving an absolute accuracy, without extrapolation, in the observables of <0.01, which is converged for many practical purposes. We thereby identify inaccuracies in the experimental extrapolations and suggest how future experiments can be implemented to increase the classical hardness.

Science & Technology - Other Topics↗

Towards determination of the strong coupling $α_s(m_Z)$ from four-flavor lattice QCD using the continuous $β$-function method

The precise value of the strong coupling $α_s(m_{Z})$ at the $Z$-boson mass $m_{Z}$ is essential for high-energy phenomenology and precision tests of quantum chromodynamics (QCD). We present the status of a program targeting a $\sim 0.3\%$ determination of $α_s(m_{Z})$ using the renormalization group $β$-function in the infinite volume gradient flow scheme based on lattice QCD simulations of degenerate four-flavor highly improved staggered quark (HISQ) ensembles. In particular, we analyze both tree-level cutoff effects and finite-mass effects. We also outline the next steps of the analysis, including the infinite-volume and continuum extrapolations required for a precise determination of $α_s(m_Z)$.

Mandlecha, Yash (ORCID:000000020587962X)↗