Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2019

Physics-Informed CoKriging: A Gaussian-Process-Regression-Based Multifidelity Method for Data-Model Convergence

In this work, we propose a new Gaussian process regression (GPR)-based multifidelity method: physics-informed CoKriging (CoPhIK). In CoKriging-based multifidelity methods, the quantities of interest are modeled as linear combinations of multiple parameterized stationary Gaussian processes (GPs), and the hyperparameters of these GPs are estimated from data via optimization. In CoPhIK, we construct a GP representing low-fidelity data using physics-informed Kriging (PhIK), and model the discrepancy between low- and highfidelity data using a parameterized GP with hyperparameters identified via optimization. Our approach reduces the cost of optimization for inferring hyperparameters by incorporating partial physical knowledge. We prove that the physical constraints in the form of deterministic linear operators are satisfied up to an error bound. Furthermore, we combine CoPhIK with a greedy active learning algorithm for guiding the selection of additional observation locations. The eciency and accuracy of CoPhIK are demonstrated for reconstructing the partially observed modified Branin function, reconstructing the sparsely observed state of a steady state heat transport problem, and learning a conservative tracer distribution from sparse tracer concentration measurements.

Yang, Xiu↗
DOE OSTI2026

First study of the nuclear response to fast hadrons via angular correlations between pions and slow protons in electron-nucleus scattering

We report on the first measurement of angular correlations between high-energy pions and slow protons in electron-nucleus (eA) scattering, providing a new probe of how a nucleus responds to a fast-moving quark. The experiment employed the CLAS detector with a 5-GeV electron beam incident on deuterium, carbon, iron, and lead targets. For heavier nuclei, the pion-proton correlation function is more spread-out in azimuth than for lighter ones, and this effect is more pronounced in the πp channel than in earlier ππ studies. The proton-to-pion yield ratio likewise rises with nuclear mass, although the increase appears to saturate for the heaviest targets. These trends are qualitatively reproduced by state-of-the-art eA event generators, including BeAGLE, eHIJING, and GiBUU, indicating that current descriptions of target fragmentation rest on sound theoretical footing. At the same time, the precision of our data exposes model-dependent discrepancies, delineating a clear path for future improvements in the treatment of cold-nuclear matter effects in eA scattering.

Correlations↗
DOE OSTI2024

Synthesis and characterization of isotopically barcoded nickel, molybdenum, and tungsten taggants for intentional nuclear forensics

Intentional nuclear forensics is a concept wherein the deliberate addition of benign and persistent material signatures to nuclear material can be used to reduce the time between the discovery of material outside of regulatory control and determination of its original provenance. One concept within intentional nuclear forensics involves the use of perturbed stable isotopes to generate unique isotope ratio “barcodes” to encode information (e.g., production batch, location, etc.) and track material throughout the nuclear fuel cycle. Synthesis of taggant species of nickel (Ni), molybdenum (Mo), and tungsten (W) was undertaken via a double-spike mechanism, wherein two highly enriched isotopes of interest per elemental taggant were mixed to form an enriched “double-spike” which was subsequently isotopically diluted with bulk material having a natural isotopic composition. Two taggant species perturbing isotopic ratios, alpha (α) and beta (β), for each of Ni, Mo, and W were synthesized. Independent measurements of double spikes and alpha and beta taggant species agreed within uncertainty and are clearly resolvable from natural compositions. High-precision analyses were independently performed by MC-ICP-MS at two U.S. National Laboratories, with consensus values and uncertainties calculated for all samples. Observed isotopic perturbations in the final taggant species measured on the order of hundreds to thousands of permille (‰) with respect to natural for isotope ratios of interest (e.g., 60 Ni/ 58 Ni, 100 Mo/ 98 Mo, 186 W/ 183 W). Discrepancies between modeled and measured isotopic compositions were observed and are largely attributed to imprecise vendor assay values for starting materials. Using measured starting material compositions as inputs for the mixing model improved the level of agreement between predicted and measured α and β taggant isotope ratios. Overall, characterization of all taggant species demonstrates that this “barcode” concept could have viability for use in nuclear forensics. Finally, it is expected that for any two-isotope mixing array dozens of isotopic barcodes could be encoded into a material system and subsequently resolved utilizing modern mass spectrometric methods.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2019

Constraining a Historical Black Carbon Emission Inventory of the United States for 1960–2000

Abstract We present an observationally constrained United States black carbon emission inventory with explicit representation of activity and technology between 1960 and 2000. We compare measured coefficient of haze data in California and New Jersey between 1965 and 2000 with predicted concentration trends and attribute discrepancies between observations and predicted concentrations among several sources based on seasonal and weekly patterns in observations. Emission factors for sources with distinct fuel trends are then estimated by comparing fuel and concentration trends and further substantiated by in‐depth examination of emission measurements. We recommend (1) increasing emission factors for preregulation vehicles by 80–250%; (2) increasing emission factors for residential heating stoves and boilers by 70% to 200% for 1980s and before; (3) explicitly representing naturally aspired off‐road engines for 1980s and before; and (4) explicitly representing certified wood stoves after 1985. We also evaluate other possible sources for discrepancy between model and measurement, including bias in modeled meteorology, subgrid spatial heterogeneity of concentrations, and inconsistencies in reported fuel consumption. The updated U.S. emissions are higher than the a priori estimate by 80% between 1960 and 1980, totaling 690 Gg/year in 1960 and 620 Gg/year in 1970 (excluding open burning). The revised inventory shows a strongly decreasing trend that was present in the observations but missing in the a priori inventory.

Sun, T.↗
DOE OSTI2022

Transient vibration and product formation of photoexcited CS 2 measured by time-resolved x-ray scattering

Here, we have observed details of the internal motion and dissociation channels in photoexcited carbon disulfide (CS 2 ) using time-resolved x-ray scattering (TRXS). Photoexcitation of gas-phase CS 2 with a 200 nm laser pulse launches oscillatory bending and stretching motion, leading to dissociation of atomic sulfur in under a picosecond. During the first 300 fs following excitation, we observe significant changes in the vibrational frequency as well as some dissociation of the C–S bond, leading to atomic sulfur in the both 1 D and 3 P states. Beyond 1400 fs, the dissociation is consistent with primarily 3 P atomic sulfur dissociation. This channel-resolved measurement of the dissociation time is based on our analysis of the time-windowed dissociation radial velocity distribution, which is measured using the temporal Fourier transform of the TRXS data aided by a Hough transform that extracts the slopes of linear features in an image. The relative strength of the two dissociation channels reflects both their branching ratio and differences in the spread of their dissociation times. Measuring the time-resolved dissociation radial velocity distribution aids the resolution of discrepancies between models for dissociation proposed by prior photoelectron spectroscopy work.

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

Magnification, dust and time-delay constraints from the first resolved strongly lensed Type Ia supernova iPTF16geu

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply-imaged Type Ia supernova iPTF16geu, at z = 0.409, using Hubble Space Telescope (HST) observations in combination with supporting ground-based data. Multi-band photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using HST and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of $\mu = 67.8^{+2.6}_{-2.9}$, accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: -0.23 ± 0.99, -1.43 ± 0.74 and 1.36 ± 1.07 days. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: $\Delta m_1 = -3.88^{+0.07}_{-0.06}$, $\Delta m_2 = -2.99^{+0.09}_{-0.08}$, $\Delta m_3 = -2.19^{+0.14}_{-0.15}$ and $\Delta m_4 = -2.40^{+0.14}_{-0.12}$ mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro- or millilensing, in addition to the large scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cosmological implications of DESI DR2 BAO measurements in light of the latest ACT DR6 CMB data

We report cosmological results from the Dark Energy Spectroscopic Instrument (DESI) measurements of baryon acoustic oscillations (BAO) when combined with recent data from the Atacama Cosmology Telescope (ACT). By jointly analyzing ACT and Planck data and applying conservative cuts to overlapping multipole ranges, we assess how different 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 + ACT dataset combinations affect consistency with DESI. While ACT alone exhibits a tension with DESI exceeding 3⁢𝜎 within the Λ⁢ CDM model, this discrepancy is reduced when ACT is analyzed in combination with Planck. For our baseline DESI DR2 BAO + 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 P⁢R⁢4 + ACT likelihood combination, the preference for evolving dark energy over a cosmological constant is about 3⁢𝜎, increasing to over 4⁢𝜎 with the inclusion of type Ia supernova data. While the dark energy results remain quite consistent across various combinations of Planck and ACT likelihoods with those obtained by the DESI collaboration, the constraints on neutrino mass are more sensitive, ranging from ∑𝑚 𝜈 < 0.061 eV in our baseline analysis, to ∑𝑚 𝜈 < 0.077 eV (95% confidence level) in the CMB likelihood combination chosen by ACT when imposing the physical prior ∑𝑚 𝜈 > 0 eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗
DOE OSTI2025

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Characterization of the Polarization Beam Response of SPT-3G Using Point Sources

Precise measurements of cosmic microwave background (CMB) polarization require rigorous control of instrumental systematics. For the South Pole Telescope's third-generation camera (SPT-3G), which observes in three bands centered near 95, 150, and 220 GHz, accurate beam characterization is critical for interpreting the polarized mm-wave sky. We present direct measurements of SPT-3G's polarized beam response from observations of 100 bright extragalactic point sources. Previous SPT-3G power spectrum analyses introduced a phenomenological parameter, $β_{\rm pol}$, to describe the polarization preserved in beam sidelobes, and found evidence for significant depolarization from the requirement of inter-frequency polarization power spectrum consistency. Our direct measurements yield $β_{\rm pol}=0.89\pm0.10$ at 95 GHz, $1.08\pm0.10$ at 150 GHz, and $0.90\pm0.22$ at 220 GHz, indicating minimal sidelobe depolarization. We validate these results with systematic tests of posterior sampling versus bootstrap resampling, real-space versus Fourier-space analysis, temperature-to-polarization leakage handling, covariance determination, and source selection. Compared to values inferred from previous cosmological analyses, our results differ by an effective $1.3σ$. This apparent discrepancy is model dependent, because the point source analysis derives much of its $β_{\rm pol}$ constraining power from higher multipoles than the power spectrum analysis. These measurements therefore admit three explanations for the frequency-dependent residuals observed in the power spectrum analysis: a statistical fluctuation, the need for more sophisticated polarized beam models, or systematics other than beam depolarization.

de Haan, T. [KEK, Tsukuba]↗
DOE OSTI2025

SpecDis: Value Added Distance Catalog for 4 Million Stars from DESI Year-1 Data

We present the SpecDis value-added stellar distance catalog accompanying DESI Data Release 1. SpecDis trains a feed-forward neural network (NN) with Gaia parallaxes and gets the distance estimates. To build up an unbiased training sample, we do not apply selections on parallax error or signal-to-noise (S/N) of the stellar spectra, and instead, we incorporate parallax error into the loss function. Moreover, we employ principal component analysis to reduce the noise and dimensionality of stellar spectra. Validated by independent external samples of member stars with precise distances from globular clusters, dwarf galaxies, stellar streams, combined with blue horizontal branch stars, we demonstrate that our distance measurements show no significant bias up to 100 kpc, and are much more precise than Gaia parallax beyond 7 kpc. The median distance uncertainties are 23%, 19%, 11%, and 7% for S/N < 20, 20 ≤ S/N < 60, 60 ≤ S/N < 100, and S/N ≥ 100. Selecting stars with ${\mathrm{log}}\,g\lt 3.8$ and distance uncertainties smaller than 25%, we have more than 74,000 giant candidates within 50 kpc of the Galactic center and 1500 candidates beyond this distance. Additionally, we develop a Gaussian mixture model to identify unresolvable equal-mass binaries by modeling the discrepancy between the NN-predicted and the geometric absolute magnitudes from Gaia parallaxes and identify 120,000 equal-mass binary candidates. Our final catalog provides distances and distance uncertainties for >4 million stars, offering a valuable resource for Galactic astronomy.

astronomy data analysis↗
NASA NTRS1978

Scale covariant cosmology and the temperature of the earth

Geological data are used as cosmological determinants in a study of the temperature of the early earth (2.3 to 4.5 billion years ago). It is known that the energy output of the sun during that period was on the order of 30-40% lower than at present, and deduced that the mean temperature of the earth should have fallen to as low as 245 K, i.e., below the freezing point of seawater. Strong evidence exists, however, to indicate that algae (therefore liquid water) was present. To reconcile the discrepancies, a model is proposed whereby terrestrial G and M vary. It is further noted that atmosphere H2 may be a better agent than NH3 for producing a greenhouse effect.

Canuto, V.↗
NASA NTRS1987

Nitrous oxide in the tropical middle atmosphere, observed by ground-based mm-wave spectrometry

Measurements of stratospheric N2O were made from Mauna Kea in Hawaii in June 1983, and in May and June 1986, by observing thermal emission of the molecule in a rotational transition at about 1 mm wavelength. Analyses of the data yield altitude profiles in the middle and upper stratosphere. Useful measurements of N2O may be made in one to two hours. The N2O profiles agree reasonably well with model predictions and with published satellite data, though significantly more N2O is reported near the stratopause than shown by the satellite measurement, and significantly more N2O in the middle stratosphere than in one of the models. The discrepancy between these data and the satellite measurement may be due in part to variations induced by the solar cycle.

Connor, Brian J.↗
NASA NTRS1990

Comparison between reference atmosphere winds and radar winds from selected locations

Zonal and meridional 60-110-km wind profiles obtained by radar measurements at Saskatoon, Adelaide, Christchurch, Puerto Rico, and Mawson are presented graphically and compared with those from the COSPAR International Reference Atmosphere (CIRA) for 1986. Good general agreement is found below about 80 km, but above 80 km the CIRA 1986 models show discrepancies, including: (1) no spring tongue of weak westward flow at latitudes 20-70 deg; (2) too strong an eastward flow at 20-52 deg in summer; (3) too great reversal heights at 35-43 deg N in summer; and (4) too strong (by a factor of 2) summer and winter jets at 65-70 deg N.

Manson, A. H.↗
NASA NTRS1993

Measuring rapid ocean tidal earth orientation variations with very long baseline interferometry

Ocean tidal effects on universal time and polar motion (UTPM) are investigated at four nearly diurnal (K(sub 1), P(sub 1), O(sub 1), and Q(sub 1)) and four nearly semidiurnal (K(sub 2), S(sub 2), M(sub 2), and N(sub 2)) frequencies by analyzing very long baseline interferometry (VLBI) data extending from 1978 to 1992. We discuss limitations of comparisons between experiment and theory for the retograde nearly diurnal polar motion components due to their degeneracy with prograde components of the nutation model. Estimating amplitudes of contributions to the modeled VLBI observables at these eight frequencies produces a statistically highly significant improvement of 7 mm to the residuals of a fit to the observed delays. Use of such an improved UTPM model also reduces the 14-30 mm scatter of baseline lengths about a time-linear model of tectonic motion by 3-14 mm, also withhigh significance levels. A total of 28 UTPM ocean tidal amplitudes can be unambiguously estimated from the data, with resulting UTI and PM magnitudes as large as 21 micro secs and 270 microarc seconds and formal uncertainties of the order of 0.3 micro secs and 5 microarc secs for UTI and PM, respectively. Empirically determined UTPM amplitudes and phases are com1pared to values calculated theoretically by Gross from Seiler's global ocean tide model. The discrepancy between theory and experiment is larger by a factor of 3 for UTI amplitudes (9 micro secs) than for prograde PM amplitudes (42 microarc secs). The 14-year VLBI data span strongly attenuates the influence of mismodeled effects on estimated UTPM amplitudes and phases that are not coherent with the eight frequencies of interest. Magnitudes of coherent and quasi-coherent systematic errors are quantified by means of internal consistency tests. We conclude that coherent systematic effects are many times larger than the formal uncertainties and can be as large as 4 micro secs for UTI and 60 microarc secs for polar motion. On the basis of such ealistic error estimates, 22 of the 31 fitted UTPM ocean tidal amplitudes differ from zero by more than 2 sigma.

Sovers, O. J.↗