Search NASA⌕ Search

SEARCH · Search NASA

Results for “parametric modeling”

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

Chamber Flux and Porewater Concentration of CH4, CO2 and N2O, 2018, Columbia River bank at the Hanford site, WA, USA

This is the corresponding observations data for the paper: Jorge A Villa; Garrett J Smith; Yang Ju; Lupita Renteria; Jordan C Angle; Evan Arntzen; Samuel F Harding; Huiying Ren; Xingyuan Chen; Audrey H Sawyer; Emily B Graham; James C Stegen; Kelly C Wrighton; Gil Bohrer (2020) Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river. Science of the Total Environment.Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow, intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during the low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

LaSDI: Parametric Latent Space Dynamics Identification

Enabling fast and accurate physical simulations with data has become an important area of computational physics to aid in inverse problems, design-optimization, uncertainty quantification, and other various decision-making applications. This paper presents a data-driven framework for parametric latent space dynamics identification procedure that enables fast and accurate simulations. The parametric model is achieved by building a set of local latent space model and designing an interaction among them. An individual local latent space dynamics model achieves accurate solution in a trust region. By letting the set of trust region to cover the whole parameter space, our model shows an increase in accuracy with an increase in training data. Herein we introduce two different types of interaction mechanisms, i.e., point-wise and region-based approach. Both linear and nonlinear data compression techniques are used. We illustrate the framework of Latent Space Dynamics Identification (LaSDI) enable a fast and accurate solution process on various partial differential equations, i.e., Burgers’ equations, radial advection problem, and nonlinear heat conduction problem, achieving 0 (100)x speed-up and 0 (1)% relative error with respect to the corresponding full order models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗
DOE OSTI2023

Investigation of methane and ethane pyrolysis with highly active and durable iron-alumina catalyst to produce hydrogen and valuable nano carbons: Continuous fluidized bed tests and reaction rate analysis

The goal of this research work was to evaluate a novel iron-alumina catalyst for methane and ethane pyrolysis with high conversion rates to H 2 with a long life time to address two major issues reported with previous catalysts. Fluidized bed tests demonstrated continuous hydrogen production while collecting carbon containing fines. During the 230 h fluid bed test the methane conversion rate to hydrogen at 700 °C remained above 80% for 160 h while creating valuable nano carbons. Similarly high catalyst lifetimes with high methane conversions have not been reported before. Ethane was more active than methane for pyrolysis with 100% conversion to H 2 . Transmission electron microscopy and Raman analysis of the carbon showed single and multiwall nanotubes and conglomerate nano-fibers. Ethane did not impact the allotrope of carbon formed. TGA was used to create a data array that was used for kinetic analysis and reaction model parametric regressions. Ethane had a higher rate of pyrolysis that was 2.3–2.8 times that of methane. Ethane pyrolysis was found to be more sensitive to temperature with a higher activation energy of 62 kJ/mol as compared to methane at 43 kJ/mol. The frequency factor, A i , was significantly higher for ethane at 0.73 [[kmol/m 3 -sec]/[N/m2] n ] compared to 0.03 for methane on the same basis. Both gases had a functional partial pressure order dependence of ~0.6 based on the regression analysis.

08 HYDROGEN↗
DOE OSTI2023

Predicting the creep-rupture lifetime of a cast austenitic stainless steel using Larson-Miller and Wilshire parametric approaches

An experimental dataset of just over 100 creep tests of a cast austenitic stainless steel, CF8C-Plus, was analyzed by two temperature-compensated parametric models (Larson-Miller, Wilshire et al.) to predict long-term lifetimes as functions of temperature and stress. The dataset and associated regression analyses showed greater scatter than typically found in recent similar studies of wrought Ni-based alloys by the same two models and was attributed to the microstructural inhomogeneity of the cast stainless steel. Qualitatively, the Larson-Miller formalism showed greater lifetime prediction accuracy than the Wilshire approach, with the latter model's predictive ability being particularly degraded by the presence of two very significant outlier results. This observation suggests that the Larson-Miller approach is more robust when treating rupture-time datasets that show particularly wide experimental scatter. Despite the differences in the overall predictive ability, both models yielded similar predictions of the applied stress at which CF8C-Plus would have a creep-limited lifetime of 100,000 h when loaded below the yield point.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geospatial analysis of freight accessibility and job attraction: The role of interstate ramps, airports, ports, and rail

The number of jobs within an industry is significantly influenced by geographical location, with transportation infrastructure playing a key role. While previous research has largely focused on how access to jobs affects employment, less attention has been given to how transportation infrastructure impacts business operations and job attraction. Here, this study addresses this gap by examining how the ease of transporting products to key transportation facilities affects job numbers in freight-intensive industries. Using job data from the Longitudinal Employment Household Dynamics dataset at the Census Tract level, we applied a non-parametric model to assess the impact of proximity to interstate ramps, rail intermodals, ports, and airports. Our analysis revealed that closer transportation infrastructure generally has a greater impact on employment. Specifically, interstate ramps are crucial for attracting jobs, particularly in rural areas, while airport proximity is essential for industries dealing with high-value, time-sensitive goods, as seen notably in Massachusetts. The importance of transportation facilities varies considerably across states and industries. The findings and method in this study can be used by transportation agencies for freight planning.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Effects of applied stress and grain size on creep-rupture lifetime prediction for Haynes 282 alloy

A multiple-heat set of tensile and creep-rupture data generated from a single laboratory was used to establish and evaluate predictive equations for the average creep-limited lifetime of an aged wrought Haynes 282 alloy using two different time- and temperature-compensated parametric models (those of Larson and Miller and Wilshire et al.). The range of temperatures and stresses used to produce the creep-rupture data used in the analyses revealed that two distinct regimes, above and below the proportional limit, were necessary for more accurate predictions of creep lifetimes over the range of conditions studied. Furthermore, a 2x difference in grain size amongst the three heats revealed a breakdown in the relationship between ultimate tensile strength and creep-rupture time, thus impacting the predictive ability of the Wilshire model. Finally, lifetime predictions based on the formulations from the respective models were compared to other available Haynes 282 datasets and found to reasonably replicate the experimental trends established in other studies.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A Conditional Autoencoder for Galaxy Photometric Parameter Estimation

Astronomical photometric surveys routinely image billions of galaxies, and traditionally infer the parameters of a parametric model for each galaxy. This approach has served us well, but the computational expense of deriving a full posterior probability distribution function is a challenge for increasingly ambitious surveys. In this paper, we use deep learning methods to characterize galaxy images, training a conditional autoencoder on mock data. The autoencoder can reconstruct and denoise galaxy images via a latent space engineered to include semantically meaningful parameters, such as brightness, location, size, and shape. Our model recovers galaxy fluxes and shapes on mock data with a lower variance than the Hyper Suprime-Cam photometry pipeline, and returns reasonable answers even for inputs outside the range of its training data. When applied to data in the training range, the regression errors on all extracted parameters are nearly unbiased with a variance near the Cramr-Rao bound.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Accelerating full-waveform inversion using source stacking: synthetic experiments at the global scale in a realistic 3-D earth model

SUMMARY The spectral element method is currently the method of choice for computing accurate synthetic seismic wavefields in realistic 3-D earth models at the global scale. However, it requires significantly more computational time, compared to normal mode-based approximate methods. Source stacking, whereby multiple earthquake sources are aligned on their origin time and simultaneously triggered, can reduce the computational costs by several orders of magnitude. We present the results of synthetic tests performed on a realistic radially anisotropic 3-D model, slightly modified from model SEMUCB-WM1 with three component synthetic waveform ‘data’ for a duration of 10 000 s, and filtered at periods longer than 60 s, for a set of 273 events and 515 stations. We consider two definitions of the misfit function, one based on the stacked records at individual stations and another based on station-pair cross-correlations of the stacked records. The inverse step is performed using a Gauss–Newton approach where the gradient and Hessian are computed using normal mode perturbation theory. We investigate the retrieval of radially anisotropic long wavelength structure in the upper mantle in the depth range 100–800 km, after fixing the crust and uppermost mantle structure constrained by fundamental mode Love and Rayleigh wave dispersion data. The results show good performance using both definitions of the misfit function, even in the presence of realistic noise, with degraded amplitudes of lateral variations in the anisotropic parameter ξ. Interestingly, we show that we can retrieve the long wavelength structure in the upper mantle, when considering one or the other of three portions of the cross-correlation time series, corresponding to where we expect the energy from surface wave overtone, fundamental mode or a mixture of the two to be dominant, respectively. We also considered the issue of missing data, by randomly removing a successively larger proportion of the available synthetic data. We replace the missing data by synthetics computed in the current 3-D model using normal mode perturbation theory. The inversion results degrade with the proportion of missing data, especially for ξ, and we find that a data availability of 45 per cent or more leads to acceptable results. We also present a strategy for grouping events and stations to minimize the number of missing data in each group. This leads to an increased number of computations but can be significantly more efficient than conventional single-event-at-a-time inversion. We apply the grouping strategy to a real picking scenario, and show promising resolution capability despite the use of fewer waveforms and uneven ray path distribution. Source stacking approach can be used to rapidly obtain a starting 3-D model for more conventional full-waveform inversion at higher resolution, and to investigate assumptions made in the inversion, such as trade-offs between isotropic, anisotropic or anelastic structure, different model parametrizations or how crustal structure is accounted for.

Geochemistry & Geophysics↗
DOE OSTI2022

Dark Energy Survey Year 3 results: marginalization over redshift distribution uncertainties using ranking of discrete realizations

ABSTRACT Cosmological information from weak lensing surveys is maximized by sorting source galaxies into tomographic redshift subsamples. Any uncertainties on these redshift distributions must be correctly propagated into the cosmological results. We present hyperrank, a new method for marginalizing over redshift distribution uncertainties, using discrete samples from the space of all possible redshift distributions, improving over simple parametrized models. In hyperrank, the set of proposed redshift distributions is ranked according to a small (between one and four) number of summary values, which are then sampled, along with other nuisance parameters and cosmological parameters in the Monte Carlo chain used for inference. This approach can be regarded as a general method for marginalizing over discrete realizations of data vector variation with nuisance parameters, which can consequently be sampled separately from the main parameters of interest, allowing for increased computational efficiency. We focus on the case of weak lensing cosmic shear analyses and demonstrate our method using simulations made for the Dark Energy Survey (DES). We show that the method can correctly and efficiently marginalize over a wide range of models for the redshift distribution uncertainty. Finally, we compare hyperrank to the common mean-shifting method of marginalizing over redshift uncertainty, validating that this simpler model is sufficient for use in the DES Year 3 cosmology results presented in companion papers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Systematic exploration of heavy element nucleosynthesis in protomagnetar outflows

ABSTRACT We study the nucleosynthesis products in neutrino-driven winds from rapidly rotating, highly magnetized and misaligned protomagnetars using the nuclear reaction network SkyNet. We adopt a semi-analytic parametrized model for the protomagnetar and systematically study the capabilities of its neutrino-driven wind for synthesizing nuclei and eventually producing ultra-high energy cosmic rays (UHECRs). We find that for neutron-rich outflows (Ye < 0.5), synthesis of heavy elements ($\overline{A}\sim 20-65$) is possible during the first $\sim 10\, {\rm s}$ of the outflow, but these nuclei are subjected to composition-altering photodisintegration during the epoch of particle acceleration at the dissipation radii. However, after the first $\sim 10\, {\rm s}$ of the outflow, nucleosynthesis reaches lighter elements ($\overline{A}\sim 10-50$) that are not subjected to subsequent photodisintegration. For proton-rich (Ye ≥ 0.5) outflows, synthesis is more limited ($\overline{A}\sim 4-15$). These suggest that while protomagnetars typically do not synthesize nuclei heavier than second r-process peak elements, they are intriguing sources of intermediate/heavy mass UHECRs. For all configurations, the most rapidly rotating protomagnetars are more conducive for nucleosynthesis with a weaker dependence on the magnetic field strength.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Quasielastic electromagnetic scattering cross sections and world data comparisons in the GENIE Monte Carlo event generator

The usage of Monte Carlo neutrino event generators (MC$\textit{ν}$EGs) is a norm within the high-energy $\textit{ν}$ scattering community. The relevance of quasielastic (QE) energy regimes to $\textit{ν}$ oscillation experiments implies that accurate calculations of $\textit{ν}$A cross sections in this regime will be a key contributor to reducing the systematic uncertainties affecting the extraction of oscillation parameters. In spite of this, many MC$\textit{ν}$EGs utilize highly phenomenological, parametrized models of QE scattering cross sections. Moreover, a culture of validation of MC$\textit{ν}$EGs against prolific electron ($\textit{e}$) scattering data has been historically lacking. In this work, we implement new eA cross sections obtained from nuclear ab initio approaches in GENIE, the primary MC$\textit{ν}$EG utilized by the FNAL community. In particular, we utilize results from quantum MC methods which solve the many-body nuclear problem in the short-time approximation (STA), allowing consistent retention of two-nucleon dynamics which are crucial to explain available nuclear electromagnetic (electroweak) data over a wide range of energy and momentum transfers. This new implementation in GENIE is fully tested against the world QE electromagnetic data, finding agreement with available data below ~2 GeV of beam energy with the aid of a scaling function formalism. The STA is currently limited to study $\textit{A}$ ≤ 12 nuclei, however, its semi-inclusive multibody identity components are exportable to other many-body computational techniques such as auxiliary field diffusion MC which can reach $\textit{A}$ ≤ 40 systems while continuing to realize the factorization contained within the STA’s multinucleon dynamics. Together, these developments promise to make future experiments such as DUNE more accurate in their assessment of MC$\textit{ν}$EG systematics, $\textit{ν}$ properties, and potentially empower the discovery of physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Small-scale clumping at recombination and the Hubble tension

Despite the success of the standard ΛCDM model of cosmology, recent data improvements have made tensions emerge between low- and high-redshift observables, most importantly in determinations of the Hubble constant H 0 and the (rescaled) clustering amplitude S 8 . The high-redshift data, from the cosmic microwave background (CMB), crucially relies on recombination physics for its interpretation. In this work, we study how small-scale baryon inhomogeneities (i.e., clumping) can affect recombination and consider whether they can relieve both the H 0 and S 8 tensions. Such small-scale clumping, which may be caused by primordial magnetic fields or baryon isocurvature below kpc scales, enhances the recombination rate even when averaged over larger scales, shifting recombination to earlier times. We introduce a flexible clumping model, parametrized via three spatial zones with free densities and volume fractions, and use it to study the impact of clumping on CMB observables. We find that increasing H 0 decreases both Ω m and S 8 , which alleviates the S 8 tension. On the other hand, the shift in Ω m is disfavored by the low-z baryon-acoustic-oscillations measurements. We find that the clumping parameters that can change the CMB sound horizon enough to explain the H 0 tension also alter the damping tail, so they are disfavored by current Planck 2018 data. We test how the CMB damping-tail information rules out changes to recombination by first removing ℓ > 1000 multipoles in Planck data, where we find that clumping could resolve the H 0 tension. Furthermore, we make predictions for future CMB experiments, as their improved damping-tail precision can better constrain departures from standard recombination. Both the Simons Observatory and CMB-S4 will provide decisive evidence for or against clumping as a resolution to the H 0 tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Neutrino-nucleus CC 0 π cross-section tuning in GENIE v3

This article summarizes the state of the art of v μ and $\bar{v}_μ$ CC0π cross-section measurements on carbon and argon and discusses the relevant nuclear models, parametrizations and uncertainties in GENIE v3. The CC0π event topology is common in experiments at a few-GeV energy range. Although its main contribution comes from quasi-elastic interactions, this topology is still not well understood. The GENIE global analysis framework is exploited to analyze CC0π datasets from MiniBooNE, T2K and MINERvA. A partial tune for each experiment is performed, providing a common base for the discussion of tensions between datasets. The results offer an improved description of nuclear CC0π datasets as well as data-driven uncertainties for each experiment. This work is a step towards a GENIE global tune that improves our understanding of neutrino interactions on nuclei. It follows from earlier GENIE work on the analysis of neutrino scattering datasets on hydrogen and deuterium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Long-exposure NuSTAR constraints on decaying dark matter in the Galactic halo

Here, we present two complementary NuSTAR x-ray searches for keV-scale dark matter decaying to monoenergetic photons in the Milky Way halo. In the first, we utilize the known intensity pattern of unfocused stray light across the detector planes—the dominant source of photons from diffuse sources—to separate astrophysical emission from internal instrument backgrounds using ~7- Ms/detector deep blank-sky exposures. In the second, we present an updated parametric model of the full NuSTAR instrument background, allowing us to leverage the statistical power of an independent ~20-Ms/detector stacked exposures spread across the sky. Finding no evidence of anomalous x-ray lines using either method, we set limits on the active-sterile mixing angle sin 2 ⁡(2⁢θ) for sterile-neutrino masses 6–40 keV. The first key result is that we strongly disfavor a ~7-keV sterile neutrino decaying into a 3.5-keV photon. The second is that we derive leading limits on sterile neutrinos with masses ~15–18 keV and ~25–40 keV, reaching or extending below the big bang nucleosynthesis limit. In combination with previous results, the parameter space for the neutrino minimal standard model is now nearly closed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Correlated anharmonicity and dynamic disorder control carrier transport in halide perovskites

Halide pervoskites are an important class of semiconducting materials that hold great promise for optoelectronic applications. In this work we investigate the relationship between vibrational anharmonicity and dynamic disorder in this class of solids. Via a multiscale model parametrized from first-principles calculations, we demonstrate that the non-Gaussian lattice motion in halide perovskites is microscopically connected to the dynamic disorder of overlap fluctuations among electronic states. This connection allows us to rationalize the emergent differences in temperature-dependent mobilities of prototypical MAPbI 3 and MAPbBr 3 compounds across structural phase transitions, in agreement with experimental findings. In conclusion, our analysis suggests that the details of vibrational anharmonicity and dynamic disorder can complement known predictors of electronic conductivity and can provide structure-property guidelines for the tuning of carrier transport characteristics in anharmonic semiconductors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Multiwavelength anomalous X-ray diffraction for combined imaging of atomic displacement and strain

The X-ray Bragg coherent diffractive imaging (CDI) technique assumes that the structure factor holds constant over the measured crystal. This approximation breaks down for materials exhibiting variations in the unit-cell configuration, such as piezo- and ferroelectrics. In that case, the strain field cannot be reliably determined from the reconstruction because the lattice deformation and the structure factor contribute concomitantly. Proposed here is a solution to this problem achieved by combining Bragg CDI and the multiwavelength anomalous diffraction approach that measures a Friedel pair of reflections at two different photon energies near an absorption edge. Comparing the obtained reconstructions with a parametric model that includes calculating the scattering amplitude as a function of wavelength and the unit-cell configuration, the contributions of the lattice deformation and the structure factor are separated. Simulations of the ferroelectric material BaTiO 3 demonstrate the possibility of simultaneous probing of the strain and displacement of the Ti atoms. The proposed method opens up an opportunity to apply coherent X-ray diffraction for nanoscale-resolved 3D mapping of polarization domains in micro- and nanocrystals.

lattice deformation↗
DOE OSTI2020

Parametric Design Study of a Power Electronics Package for Improving Solder Joint Reliability

The lifetime of a power electronics package is, to a large extent, determined by the reliability of its bonded interfaces under major loading conditions in an operational environment. Based on the application-level requirements of the package, bonded interface material is selected based on the results obtained from accelerated tests such as thermal cycling and power cycling. In addition to evaluating the reliability through accelerated tests, it is important to consider the impact of other component layers on the thermomechanical performance of the interface material, both from a material and geometric perspective. The co-efficient of thermal expansion (CTE) mismatch introduced by the use of different materials within a package and its structural design plays a critical role in determining the interface material reliability. In this paper, we present the results of a parametric modeling study of a power electronics package under thermal cycling in which the materials and geometric design of the different component layers were varied with respect to a baseline design to understand their impact on the reliability of the interface material. We chose the volume-averaged strain energy density per cycle computed at the corner region of the interface material as the metric for the reliability comparisons. Our results indicate that in addition to the CTE mismatch, the stiffness of individual component layers has a major impact on reliability. Among the different baseplates that we studied, aluminum silicon-carbide baseplates offered superior reliability over their copper and aluminum counterparts. We also found that the magnitude of the impact of stiffness variation—amongst the adjacent layers—on the reliability of the solder joint is dependent on the inherent CTE mismatch between the layers.

47 OTHER INSTRUMENTATION↗