Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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

Calculating Average Hot Water Mixes of Residential Plumbing Fittings Using the ANSI 301-2019 Hot Water Draw Model and National Residential Data to Estimate Hot Water Use in Showerheads and Lavatory Faucets

Available studies of hot water use percentages may not necessarily be generalized to a national level given regional differences, varying methodological approaches, and limited sample sizes. While some publicly available usage data, reports, and surveys estimate actual household usage, a lack of water use data specific to end uses make more precise savings calculations difficult. Additionally, most hot water draw models tend to focus on overall household use and may not readily estimate lavatory fixture use. However, reviewing more recent studies and standards enables the U.S. Environmental Protection Agency’s WaterSense program to update its estimates for hot water use and consumers’ corresponding energy and monetary savings to better reflect realworld conditions. This report specifically focuses on improving hot water use and savings estimates of lavatory faucets and showerheads. To estimate the hot water use of lavatory faucets and showerheads, we employed the 2015 Residential Energy Consumption Surveys (RECS) microdata alongside the ANSI 301-2019 Hot Water Draw Model. These estimates account for regional differences in hot water use, including regional cold water inlet temperatures. As a result, the refinements presented in this report are more robust, more recent, and better describe the geographic variation than previous inputs used by the WaterSense program. In addition, the approach described in this paper can be updated over time or tailored to regionally specific needs given available inputs. We conclude that hot water percentages for showers and faucets calculated using publicly available national data are close to the percentages found by regional studies and are consistent with household-level models of water use.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

DECOVALEX-2023: Task E Final Report

This is the Task E final report for DECOVALEX-2023. Task E is focused on understanding thermal, two-phase hydrological, and mechanical (TH2M) processes, especially related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and essentially non-porous in the far field (away from excavations). Investigation of the short-term (days to years) near-field (centimeters to tens of meters) behavior of salt is important for radioactive waste disposal because this early period strongly controls the amount of brine in a salt repository. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment. The main test case used in Task E is the ongoing Brine Availability Test in Salt (BATS) heater test located underground at the Waste Isolation Pilot Plant (WIPP) near Carlsbad, New Mexico, USA. The Task was divided into a series of Steps. Step 0 was an introduction to processes in salt, that included matching historical unheated brine inflow data from boreholes at WIPP and matching temperature observations during BATS heater test 1a. Step 1 included validation of models against a thermo-poroelastic analytical solution relevant to heated boreholes in salt, and two-phase flow around an excavation in salt. Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the BATS 1a heater test. There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. The key learning points from Task E have been: • Heat conduction through salt typically requires non-linear thermal conductivity (as a function of temperature), but most models do a good job matching observations, given appropriate adjustments to the applied power and some thermocouple locations. • Thermal pressurization requires coupled thermal-hydrological-mechanical (THM) responses that consider the thermal expansion of the fluid and solid phases. • Initialization of two-phase flow models around a borehole or excavation in salt are more realistically represented as “wetting up”, rather than “drying down” (i.e., the initial state after excavation is mostly dry, rather than mostly wet). • The BATS 1a heater test includes a significant release of brine after the end of heating, which requires a large increase in permeability to recreate. Task E has been a great learning experience for all the teams involved, and feedback from the modeling teams has led to changes in the design of follow-on BATS experiments, which are now ongoing underground at WIPP. There was a balance throughout the task between freedom to model phenomena how each team saw fit, and prescriptiveness in problem design to bring the modeling teams closer together to allow attribution of smaller differences between models to different modeling choices. The modeling approaches seem to go through two phases: an early phase of discovery or testing, and a later phase of refinement and improvement. In future modeling efforts, different field data could be used (e.g., BATS 2) and more time should be included in the processes for teams to make multiple model refinement or even significant changes to their conceptual model or setup, based on lessons learned from the modeling exercise.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Particle Receiver Models for Systems Analysis

Particle receivers are gaining importance in the field of Concentrating Solar Power (CSP) due to the high temperature that particles can achieve without degradation. Several researchers are studying the potential of this technology by means of system analyses, which need simple and light models of the system. This study presents two simple models for the particle receiver. The simplest model is a correlation obtained by fitting the results calculated with a more complex receiver model simulated in CFD. The other model is a 1D model, which is benchmarked against the same CFD results. Although both models achieve high coefficient of determination, R2, when compared to CFD results, the 1D model seems to provide more accurate results (especially during sunsets and sunrises). Both models are integrated into a tecno-economic model developed in previous work. The LCOE obtained with the 1D model is between 7% and 10% greater than the one obtained with the correlation.

González-Portillo, Luis F. (ORCID:0000000348643825↗
DOE OSTI2026

Constraints on White Dwarf Hydrogen Layer Masses Using Gravitational Redshifts

The hydrogen envelope is the outermost layer of a DA white dwarf; it makes up the entirety of the stellar photosphere, and yet its typical extent is difficult to model theoretically and remains poorly observationally constrained. As a result, hydrogen envelope mass is a substantial source of systematic uncertainty in the physical properties of white dwarfs, including overall masses and cooling ages. In this work, we fit a Gaussian mixture model to gravitational redshifts from high-resolution spectroscopy, paired with radius measurements from Gaia BP/RP spectra, to measure the mass–radius relation for a sample of 468 white dwarfs. Our results are in excellent agreement with the predicted mass–radius relations of state-of-the-art evolutionary models, including those from the MESA Isochrones and Stellar Tracks (MIST) library. We find that mass–radius relations such as those from MIST that assume a thick and mass-dependent hydrogen envelope are preferred by the observed probability density function over models that assume a hydrogen envelope of constant mass. Proper treatment of the evolution of white dwarf progenitors is thus important for accurately modeling the mass–radius relation. Our results indicate that gravitational redshift measurements of large samples of white dwarfs in wide binaries are promising probes of the hydrogen envelope masses of DA white dwarfs.

Astronomy and AstroPhysics↗
DOE OSTI2023

Physics-informed machine-learning model of temperature evolution under solid phase processes

We model temperature dynamics during Shear Assisted Proccess Extrusion (ShAPE), a solid phase process that plasticizes feedstock with a rotating tool and subsequently extrudes it into a consolidated tube, rod, or wire. Control of temperature is critical during ShAPE processing to avoid liquefaction, ensure smooth extrusion, and develop desired material properties in the extruded products. Accurate modeling of the complicated thermo-mechanical feedbacks between process inputs, material temperature, and heat generation presents a significant barrier to predictive modeling and process design. In particular, connecting micro-structural scale mechanisms of heat generation to macro-scale predictions of temperature can become computationally intractable. In this work we use a neural network (NN) model of heat generation to bridge this gap, by combining it with a simplified model of the temperature dynamics due to conduction and convection to capture the macro scale evolution of temperature. We inform the construction of the NN heat generation model using crystal plasticity simulations at the micro-structural scale to model the effects of process inputs on generation of heat. We achieved close fits of the temperature dynamics model to a diverse experimental data-set. Further, the relationships learned by the NN model between process inputs and heat generation showed qualitative agreement with those predicted by crystal plasticity simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Selecting samples of galaxies with fewer Fingers-of-God

The radial positions of galaxies inferred from their measured redshift appear distorted due to their peculiar velocities. We argue that the contribution from stochastic velocities — which gives rise to `Fingers-of-God' (FoG) anisotropy in the inferred maps — does not lend itself to perturbative modelling already on scales targeted by current experiments. To get around this limitation, we propose to remove FoG using data-driven indicators of their abundance that are local in nature and thus avoid selection biases. In particular, we show that the scale where the measured power spectrum quadrupole changes sign is tightly anti-correlated with both the satellite fraction and the velocity dispersion, and can thus be used to select galaxy samples with fewer FoG. In addition, we show that excluding galaxies in haloes more massive than a given mass threshold can help to discard many of the most problematic galaxies. Such selection could be achieved in practice using maps of the thermal Sunyaev-Zel'dovich distortion of the cosmic microwave background frequency spectrum. These techniques could potentially improve reconstructions of the large-scale velocity and displacement fields from the redshift-space positions of galaxies. They may also extend the reach of perturbative models for galaxy clustering, though in practice we find only marginal gains when fitting one-loop EFTofLSS models to simulations with mitigated FoG due to the relevance of other effects entering at two-loop order.

cosmological parameters from LSS↗
DOE OSTI2023

Robust cosmological inference from non-linear scales with k -th nearest neighbour statistics

ABSTRACT We present the methodology for deriving accurate and reliable cosmological constraints from non-linear scales ($\lt 50\, h^{-1}$ Mpc) with k-th nearest neighbour (kNN) statistics. We detail our methods for choosing robust minimum scale cuts and validating galaxy–halo connection models. Using cross-validation, we identify the galaxy–halo model that ensures both good fits and unbiased predictions across diverse summary statistics. We demonstrate that we can model kNNs effectively down to transverse scales of $r_{\rm p}\sim 3\, h^{-1}$ Mpc and achieve precise and unbiased constraints on the matter density and clustering amplitude, leading to a 2 per cent constraint on σ8. Our simulation-based model pipeline is resilient to varied model systematics, spanning simulation codes, halo finding, and cosmology priors. We demonstrate the effectiveness of this approach through an application to the Beyond-2p mock challenge. We propose further explorations to test more complex galaxy–halo connection models and tackle potential observational systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Prediction and Validation of Flow Properties in Porous Lattice Structures

High-porosity metal foams have been extensively studied as an attractive candidate for efficient and compact heat exchanger design. With the advancements in additive manufacturing, such foams can be manufactured with controlled topology to yield highly tailorable mechanical and transport properties. In this study, a lattice Boltzmann method (LBM)-based pore-scale model is implemented to simulate the fluid flow in additively manufactured (AM) metal foams with unit cell topologies of Cube, Face Diagonal (FD)-Cube, Tetrakaidecahedron (TKD), and Octet lattices. The pressure gradient versus average velocity profiles predicted by the LBM model were validated against in-house measurements on the AM lattice samples with the same unit cell topologies. Based on the simulation results, a novel hybrid model is proposed to accurately predict the volume averaged flow properties (permeability and inertial coefficients) of the four structures. Specifically, the linear LBM (neglecting inertial forces) is first implemented to obtain the intrinsic permeability, and then the standard LBM is applied to obtain the inertial coefficient. Convenient correlations for those flow properties as a function of porosity and fiber diameter are constructed. The effects of the AM print qualities on the flow properties are also discussed. The advantages of the hybrid model compared to the polynomial fitting approach for determining flow properties are discussed and compared quantitatively. The hybrid model and presented results are valuable for flow and thermal transport evaluation when designing new metal foams for specific applications and with different materials and topologies. Finally, the presented correlations based on pore-scale simulations can also be conveniently used in volume-averaged models to predict the macroscale flow behavior in such complex structures.

42 ENGINEERING↗
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 OSTI2026

The Influence of Shallow Subsurface Properties on Particle Motion in Acoustic-Seismic Coupling

Atmospheric acoustic waves transmit energy into the solid Earth through air-to-ground coupling. These waves are recorded by seismic sensors and provide insight into both atmospheric phenomena and subsurface properties. Interpreting these signals is often challenging because they are modulated by subsurface structure and the incidence angle of the acoustic wave. This study examines acoustic--seismic coupling generated by the 2012 Camp Minden Explosion, which was recorded by hundreds of seismoacoustic stations. We apply a novel technique to quantify the seismic particle motion, model coupled waves with a propagator matrix approach, and apply a Bayesian inversion to infer properties of the shallow subsurface. Our analysis reveals that prograde motion is widespread and focused in low shear-wave velocity regions, such as the Mississippi Embayment, and retrograde motion is more common in higher shear-velocity areas. Inversion results at some stations produce plausible subsurface models with strong waveform fits, while inversions at other sites are less successful. These results indicate prograde particle motion in air-to-ground coupled waves is more prevalent than previously recognized and may serve as a diagnostic for shallow velocity structure. Our comprehensive modeling and inversion framework provides a potential method to extract layered near-surface properties from acoustic-seismic coupling observations.

58 GEOSCIENCES↗
DOE OSTI2022

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Single-channel and single-energy partial-wave analysis with continuity improved through minimal phase constraints

Single-energy partial-wave analysis has often been applied as a way to fit data with minimal model dependence. However, remaining unconstrained, partial waves at neighboring energies will vary discontinuously because the overall amplitude phase cannot be determined through single-channel measurements. This problem can be mitigated through the use of a constraining penalty function based on an associated energy-dependent fit. However, the weight given to this constraint results in a biased fit to the data. In this paper, for the first time, we explore a constraining function which does not influence the fit to data. The constraint comes from the overall phase found in multichannel fits which, in the present study, are the Bonn-Gatchina and Jülich-Bonn multichannel analyses. The data are well reproduced and weighting of the penalty function does not influence the result. The method is applied to K⁢Λ photoproduction data and all observables can be maximally well reproduced. While the employed multichannel analyses display very different multipole amplitudes, we show that the major difference between two sets of multipoles can be related to the different overall phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Difficulties with late-time solutions for the Hubble tension

We explore the notion that cosmological models that modify the late-time expansion history cannot simultaneously fit the SH0ES Collaboration’s measurements of the Hubble constant, Dark Energy Spectroscopic Instrument baryon acoustic oscillations data, and type Ia supernova distances. Adopting a few simple phenomenological models, we quantitatively demonstrate that a satisfactory fit with a model with late-time expansion history can only be achieved if one of the following is true: (1) there is a sharp step in the absolute magnitude of type Ia supernovae at very low redshift, 𝑧 ∼0.01, or (2) the distance duality relation, 𝑑 𝐿 ⁡(𝑧) =(1+𝑧) 2 ⁢𝑑 𝐴 ⁡(𝑧), is broken. Both solutions are trivial in that they effectively decouple the calibrated type Ia supernovae measurements from other data, and this qualitatively agrees with previous work built on studying specific dark-energy models. We also identify a less effective class of late-time solutions with a transition at 𝑧 ≃0.15 that lead to a more modest improvement in fit to the data than models with a very low-𝑧 transition. Our conclusions are largely unchanged when we include surface brightness fluctuation distance measurements, with their current systematic uncertainties, to our analysis. Here, we finally illustrate our findings by studying a physical model which, when equipped with the ability to smoothly change the absolute magnitude of type Ia supernovae, partially resolves the Hubble tension.

Cosmological parameters↗
DOE OSTI2023

Quintessential cosmological tensions

Several cosmological tensions have emerged in light of recent data, most notably in the inferences of the parameters and . We explore the possibility of alleviating both these tensions simultaneously by means of the Albrecht-Skordis “quintessence” potential. The field can reduce the size of the sound horizon while concurrently suppressing the power in matter density fluctuations before it comes to dominate the energy density budget today. Interestingly, this rich set of dynamics is governed entirely by one free parameter that is of in Planck units. We find that the inferred value of can be increased, while that of can be decreased, both by compared to the case. However, ultimately the model is disfavored by Planck and BAO data alone, compared to the standard model, with a . When including large scale structure and supernova data we find . We note that historically much attention has been focused on preserving the three angular scales , , and to their values. Our work presents an example of how, while doing so indeed maintains a relatively good fit to the cosmic microwave background data for an increased number of ultrarelativistic species, it is a priori insufficient in maintaining such a fit in more general model spaces.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Analytic_Deflection_of_Asteroids_by_NEDs

This package evaluates an analytic formula for the change in velocity imparted to an asteroid by a 1 or 2 keV black body source given the radius of the asteroid, the standoff distance, the x-ray yield, the density of the asteroid, and the fit coefficients for the 1 or 2 keV black body source (provided in the module). Three variations of the formula are given; the "original" form, the "modified" form that accounts for the angle of incidence in calculating the melt depth, and a form based on an "impulse" model. Each formula is fit to simulations of spherical asteroids of uniform composition illuminated by black body x-rays.

Managan, RobertA↗
DOE OSTI2020

A Novel Relative Permeability Model for Gas and Water Flow in Hydrate-Bearing Sediments With Laboratory and Field-Scale Application

In a producing gas hydrate reservoir the effective porosity available for fluid flow constantly changes with dissociation of gas hydrate. Therefore, accurate prediction of relative permeability using legacy models (e.g. Brooks-Corey (B-C), van Genuchten, etc.) that were developed for conventional oil and gas reservoirs would require empirical parameters to be calibrated at various S h over its range of variation, but such calibrations are precluded because of lack of experimental relative permeability data. This study proposes a new relative permeability model for gas hydrate-bearing media that is a function of maximum capillary pressure, capillary entry pressure, pore size distribution index, residual saturations, hydrate saturation, and four other constants. The three novel features of the proposed model are: (i) requires fitting its six empirical parameters only once using experimental data from any single S h , and the same set of empirical parameters predict relative permeability at all S h , (ii) includes the effect of capillarity, and (iii) includes the effect of pore-size distribution. From practical standpoint, the model can be used to simulate multiphase flow in gas hydrate-bearing sediments where the proposed relative permeability can account for the evolving hydrate saturation. The proposed model is implemented in a numerical simulator and the wall time required to perform simulations using the proposed model is shown to be similar to the time it takes to run same simulations with the B-C model. The proposed model is a step forward towards achieving the goal of physically accurate modeling of multiphase flow for gas hydrate-bearing sediments that accounts for the effect of gas hydrate saturation change on relative permeability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Continuous and discrete modeling of HIV-1 decline on therapy

Mathematical models have shed light on the dynamics of HIV- 1 infection in vivo. In this paper, we generalize continuous mathematical models of drug therapy for HIV-1 by Perelson et al. (Science 271:1582–1586, 1996) and Perelson and Nelson (SIAM Rev 41:3–44, 1999) on time scales, i.e., a nonempty closed subset of real numbers in order to derive new discrete models that predict the total concentration of plasma virus as a function of time. Here, one of our main goals is to compare discrete mathematical models with the continuous model in Perelson et al. (1996) where HIV infected patients were given protease inhibitors and sampled frequently thereafter. For the comparison, we use experimental data collected in Perelson et al. (1996) and estimate the parameters such as the virion clearance rate and the rate of loss of infected cells by fitting the total concentration of plasma virus to this data set. Our results show that discrete systems describe the best fit. In the previous models of this study, the efficacy of protease inhibitor is assumed to be perfect. Motivated by Perelson and Nelson (1999), we end the paper with a mathematical model of imperfect protease inhibitor and reverse transcriptase (RT) inhibitor combination therapy of HIV-1 infection on time scales with its stability analysis.

60 APPLIED LIFE SCIENCES↗