SEARCH · Search NASA

Results for “modeling 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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Chilldown↗

NASA NTRSDate not supplied

Validation Experiments for Turbulent Separation over an Axisymmetric Body of Revolution

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges to validating numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires data from advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed measurements of the flow field undergoing a smooth, adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states progressing from fully attached, to incipient separation, and finally to small-scale separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configurations of the axisymmetric model host a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation, as well as a slightly larger adverse pressure gradient, inducing a small-scale region of turbulent separation. Experiments include steady pressure measurements and 2-D PIV to provide the preliminary dataset for simulation studies, which examine the effect of variable grid domains and RANS turbulence models. This is done in an effort to understand and evaluate the critical variability between solutions for the present model configurations. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and slight variability in reference parameters. Ongoing work will focus on higher fidelity experimental campaigns to obtain surface flow visualizations and Stereo-Particle Image Velocimetry (SPIV) to deliver higher spatial resolution of the three-dimensional flow field to aid turbulence modelers.

Validation↗

NASA NTRS2009

Modeling of the 22 July 2009 Storm

The magnetic storm on 22 July 2009 is the greatest storm observed since the summer of 2008 when TWINS began its stereo imaging of the magnetosphere, On 22 July 2009, the Dst dropped to nearly -80 nT at 07:00 and 10:00 UT. During the main phase and at the peak of the storm, TWINS 1 and 2 were near apogee and moving from pre-dawn to post-dawn local time. The energetic neutral atom (ENA) Imagers on the 2 spacecraft captured the storm intensification and the formation of the partial ring current. The development of this storm has been simulated using the Comprehensive Ring Current Model (CRCM) to understand and interpret the observed signatures. The CRCM reproduced the double-dip in the Dst index and the simulated ENA flux intensities agree very well with the TWINS images. However, the peak of ion flux predicted by the model is always eastward of the observed maximum by TWINS. This discrepancy posts a challenge to reexamine the physical models employed in the CRCM.

Fok, Mei-Ching↗

DOE OSTI2025

Permafrost, Peatland, and Cropland Regions Are Key to Reconciling North American Carbon Sink Estimates

Persistent discrepancies between bottom-up, terrestrial biosphere models (TBMs), and top-down, atmospheric inversions, have made it difficult to quantify the magnitude of the North American terrestrial carbon sink. Previous studies have compared aggregated continent-scale estimates of carbon fluxes from TBMs and inversions for all of North America, but this provides limited insights into finer-scale mismatches that contribute to the overall discrepancies. Here we evaluate agreement between TBM and inversion carbon flux estimates at 1° × 1° resolution to provide more direct insights into where models disagree and what underlying factors drive discrepancies. We find that the additional carbon uptake estimated by inversions, in just 16% of the area of North America, is large enough to account for the discrepancy between TBMs and inversions across the whole continent. The majority of these differences occur in permafrost, peatland, and cropland regions. In these regions, we find a higher likelihood of potential biases in the weaker sink estimates from TBMs, suggesting that the stronger sink implied by inversions is more likely to be realistic. However, the current observational coverage is insufficient for fully assessing the causes of discrepancies or the magnitude of biases in either approach. Encouragingly, improved representation of agricultural processes in a TBM led to better agreement with inversions in croplands. Efforts to accurately model cropland dynamics will help improve agreement between TBMs and inversions. Overall, this work presents a clear path for reconciling the discrepancies between inversion and TBM estimates of the North American carbon sink that have persisted for two decades.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Lagrangian displacement field estimators in cosmology

The late-time nonlinear Lagrangian displacement field is highly correlated with the initial field, so reconstructing it could enable us to extract primordial cosmological information. Our previous work carefully studied the displacement field reconstructed from the late-time density field using the iterative method proposed and found that it does not fully converge to the true, underlying displacement field (e.g., ~8 % offset at k ~ 0.2 h Mpc –1 at z = 0.6 ). We also constructed the Lagrangian perturbation theory model for the reconstructed field, but the model could not explain the discrepancy between the true and the reconstructed fields in the previous work. Here, the main sources of the discrepancy were speculated to be a numerical artifact in the displacement estimator due to the discreteness of the sample. In this paper, we develop two new estimators of the displacement fields to reduce such a numerical discreteness effect, the normalized momentum estimator and the rescaled resumed estimator. We show that the discrepancy reported is not due to the numerical artifacts. We conclude that the method cannot fully reconstruct the shape of the nonlinear displacement field at the redshift we studied, while it is still an efficient baryon acoustic oscillation reconstruction method. In parallel, by properly accounting for the UV-sensitive term in a reconstruction procedure with an effective field theory approach, we improve the theoretical model for the reconstructed displacement field, by almost 5 times, from ~15% to the level of a few percent at k ~ 0.2 h Mpc –1 at the redshift z = 0.6.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Discrete element simulation of Pebble Bed Reactors on graphics processing units

Prediction of pebble positions in a Pebble Bed Reactor (PBR) is necessary for both reactor physics and thermal hydraulics simulations as the arrangement of pebbles has a significant impact on the resulting core power, coolant flow, and fuel temperature. Knowledge of pebble movement as the fuel is cycled through the core is also critical for predicting the fuel residence time and subsequently, the fuel burnup. Simulation with the Discrete Element Method (DEM) can provide knowledge of both the fuel packing and the fuel movement during cycling. Previous works that have performed 3D full-core DEM simulation of PBRs have used simplified models that neglect reflector wall features. This work employs a graphics processing unit (GPU)-enabled DEM code, Project Chrono, to analyze the differences in pebble packing and pebble velocities between a simplified smooth PBR reflector and a more realistic reflector that includes circular wall features. Additionally, a sensitivity study is performed on the depth of the wall features to ensure that crystallization is prevented. Project Chrono is also validated for PBR cycling applications using experimental data. It is found that wall features with a depth of at least 0.5 pebble diameters significantly reduce crystallization in the near-wall region, leading to discrepancies in both packing fraction and pebble velocity in this region compared to the simplified reflector models. These discrepancies are found to lead to roughly a 5–10% difference in the prediction of the near-wall porosity and a 10% difference in the prediction of the velocity of pebbles near the wall. As a result of these discrepancies, it is suggested that future DEM simulations of PBRs include wall features to reduce modeling errors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1978

Quantitative calculations of helium ion escape fluxes from the polar ionospheres

An attempt is made to reduce the discrepancy between the earlier theoretical model and the recent experimental observations of H(+) outward fluxes for winter and summer hemispheres by improving the theoretical model. The helium photoionization cross sections used are accurate to 10%, the latest solar EUV fluxes measured by the Atmosphere Explorer satellites are incorporated, and the most recent MSIS model of the neutral atmosphere is contained in the model. A range of conditions covering solar cycle, seasonal, and geomagnetic conditions are studied. The results show a maximum H(+) escape flux of 1.4 times 10 to the 7th per sq cm per sec for solar maximum, winter, low magnetic activity conditions, which is within the scatter of the measured fluxes. The computed summer H(+) escape flux is a factor of 20 lower than the winter value. Possible reasons for the slight discrepancy between theory and experiment in summer are discussed.

Raitt, W. J.↗

DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks. Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Bridging the Gap between Cosmological Simulations with Graph Neural Networks and Domain Adaptation

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross dataset tasks (up to 28% better relative error and up to almost an order of magnitude better χ 2 ). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Domain Adaptive Graph Neural Networks for Constraining Cosmological Parameters Across Multiple Data Sets

Deep learning models have been shown to outperform methods that rely on summary statistics, like the power spectrum, in extracting information from complex cosmological data sets. However, due to differences in the subgrid physics implementation and numerical approximations across different simulation suites, models trained on data from one cosmological simulation show a drop in performance when tested on another. Similarly, models trained on any of the simulations would also likely experience a drop in performance when applied to observational data. Training on data from two different suites of the CAMELS hydrodynamic cosmological simulations, we examine the generalization capabilities of Domain Adaptive Graph Neural Networks (DA-GNNs). By utilizing GNNs, we capitalize on their capacity to capture structured scale-free cosmological information from galaxy distributions. Moreover, by including unsupervised domain adaptation via Maximum Mean Discrepancy (MMD), we enable our models to extract domain-invariant features. We demonstrate that DA-GNN achieves higher accuracy and robustness on cross-dataset tasks (up to $28\%$ better relative error and up to almost an order of magnitude better $\chi^2$). Using data visualizations, we show the effects of domain adaptation on proper latent space data alignment. This shows that DA-GNNs are a promising method for extracting domain-independent cosmological information, a vital step toward robust deep learning for real cosmic survey data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

VALIDATION, VERIFICATION, AND CALIBRATION THROUGH A CAUSAL LENS

This paper presents an alternative method based on causal inference to perform validation, verification, and calibration of simulation models. While classical validation and verification approaches focus on the identification of the associations between data elements using statistical and machine learning methods, the novel methods in this paper focus instead on the identification of causal relationships between data elements. Statistical and machine learning-based approaches are strictly data-driven, meaning that they provide quantitative comparison measures between datasets without explicitly considering the hypotheses behind them. This can lead to the erroneous conclusion that, if two data sets are close enough, then the models that generated them are similar. In addition, when experimental and simulated data differ to an extent that fails to meet the acceptance criteria, calibration techniques are used to tweak simulation model parameters to reduce the gap between the two types of data. This produces the false expectation that a simulation model will match reality. The methods presented in this paper move away from these strictly data-driven methods for validation and calibration toward more robust, model-driven methods based on causal inference. Causal inference aims to identify the possible mechanisms that might have generated data. Thus, this analysis targets the prediction of the effects when one (or more) of the identified mechanisms are altered. There are many approaches to identify, quantify and illustrate causal relationships. For the scope of this paper, directed graphs are employed as causal models. If the directed graph lacks cycles it is known as a directed acyclic graph (DAG). A node in such a graph represents an observed data element while a directed edge connecting two nodes represents a causal relationship between two variables. The developed causal methods are designed to extract causal models from simulation models and from experimental data. Causal models capture the causal relationships between data elements (e.g., simulated and experimental data). In this context, validation and verification are performed by comparing causal models. The proposed approach does not only inform system analysts on how a simulation model matches real-world data, but also identifies elements of the simulation model that should be revised when discrepancies between simulation and experimental data are observed. Through these causal methods, analysts have a means to identify the portion of the model equation(s) that are behind an edge connecting two variables. Hence, once the structural differences between causal models have been determined, model calibration can occur by changing only those model parameters that impact the identified causal relationships.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1980

Vibrational kinetics in CO electric discharge lasers - Modeling and experiments

A model of CO laser vibrational kinetics is developed, and predicted vibrational distributions are compared with measurements. The experimental distributions were obtained at various flow locations in a transverse CW discharge in supersonic (M = 3) flow. Good qualitative agreement is obtained in the comparisons, including the prediction of a total inversion at low discharge current densities. The major area of discrepancy is an observed loss in vibrational energy downstream of the discharge which is not predicted by the model. This discrepancy may be due to three-dimensional effects in the experiment which are not included in the model. Possible kinetic effects which may contribute to vibrational energy loss are also examined.

Stanton, A. C.↗

DOE OSTI2023

Secondary Pion Production using the LINAC Beam on Graphite and Tungsten Targets

The Fermilab 400-MeV Linac primary proton beam can produce precision, single-species secondary beams using a production target in the MeV Test Area experimental hall; specifically, pions from 4 up to 120 MeV kinetic energy. (These low energy pions rapidly decay into a muon beam in a secondary collection and transport beamline.) Although graphite is a standard production target material given its ultra-high temperature tolerance and low impact on the post-target primary beam, secondary pion/muon total production cross sections increase approximately as Z1/3 for positive muons and N2/3 for negative muons. Higher mass targets such as tungsten can potentially increase low-energy m+ and m- rates by factors of 3 and 8, respectively. Initial production studies have confirmed higher muon yield from heavy targets (tungsten vs carbon), but also, surprisingly, show significant differences in pion production between modern hadronic models (GENIEhad) and between GEANT and MCNP, state of the art production modeling codes. This project will further investigate the production discrepancies between these two target materials and between different physics models with high-statistics runs using the GEANT-based code, G4Beamline to better quantify the results and compare with production data which will be obtained using the MTA secondary beamline.

Ahmed, Shiza↗

NASA NTRS2003

Coulomb Stress Accumulation along the San Andreas Fault System

Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Smith, Bridget↗

NASA NTRS2020

Bennu Shape Model Validation Methods and Results

The OSIRIS-REx Independent Shape Modeling Team, staffed by the Jet Propulsion Laboratory (JPL), has used stereophotoclinometry (SPC) to produce a three-dimensional shape model of Bennu. The SPC process is informed by (but is separate from) orbit determination, so it is important to ensure that the resulting shape model is consistent with all available spacecraft tracking data. Specifically, checking the shape model’s consistency with LIDAR measurements can illuminate any discrepancies, because LIDAR measurements are highly correlated with the shape model. This study focuses on the JPL experience with these LIDAR measurements and the greater context of the shape model validation process.

Lauretta, Dante S.↗

NASA NTRS2003

Measurement of Low Amounts of Precipitable Water Vapor Using Ground-Based Millimeterwave Radiometry

Extremely dry conditions characterized by amounts of precipitable water vapor (PWV) as as 1-2 mm commonly occur in high-latitude regions during the winter months. While such atmospheres carry only a few percent of the latent heat energy compared to tropical atmospheres, the effects of low vapor amounts on the polar radiation budget - both directly through modulation of longwave radiation and indirectly through the formation of clouds - are considerable. Accurate measurements of precipitable water vapor (PWV) during such dry conditions are needed to improve polar radiation models for use in understanding and predicting change in the climatically sensitive polar regions. To this end, the strong water vapor absorption at 183.310 GHz provides a unique means of measuring low amounts of PWV. Weighting function analysis, forward model calculations based upon a 7-year radiosonde dataset, and retrieval simulations consistently predict that radiometric measurements made using several millimeter-wavelength (MMW) channels near the 183 GHz line, together with established microwave (MW) measurements at the 22.235 GHz water vapor line and -3 1 GHz atmospheric absorption window can be used to determine within 5% uncertainty the full range of PWV expected in the Arctic. This unique collective capability stands in spite of accuracy limitations stemming from uncertainties due to the sensitivity of the vertical distribution of temperature and water vapor at MMW channels. In this study the potential of MMW radiometry using the 183 GHz line for measuring low amounts of PWV is demonstrated both theoretically and experimentally. The study uses data obtained during March 1999 as part of an experiment conducted at the Department of Energy s Cloud and Radiation Testbed (CART) near Barrow, Alaska. Several radiometers from both NOAA and NASA were deployed during the experiment to provide the first combined MMW and MW ground-based data set during dry arctic conditions. Single-channel retrievals of PWV were performed using the MW and MMW data. Discrepancies in the retrieved values were found to be consistent with differences observed between measured brightness temperatures (TBs) and forward-modeled TBs based on concurrent radiosonde profiles. These discrepancies are greater than can be explained by measurement error alone and are attributed to absorption model uncertainty. We discuss here the measurements, retrieval technique, and line model discrepancies along with difficulties and potential of MMW/MW PWV measurement.

Racette, Paul E.↗

NASA NTRS1980

Estimation of the accuracy of dynamic flight-determined coefficients

This paper discusses means of assessing the accuracy of maximum likelihood parameter estimates obtained from dynamic flight data. The commonly used analytical predictors of accuracy are compared from both statistical and simplified geometric standpoints. Emphasizing practical considerations, such as modeling error, the accuracy predictions are evaluated with real and simulated data. Improved computations of the Cramer-Rao bound to correct large discrepancies caused by colored noise and modeling error are presented. This corrected Cramer-Rao bound is the best available analytical predictor of accuracy. Engineering judgement, aided by such analytical tools, is the final arbiter of accuracy estimation.

Maine, R. E.↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 2. Quantitative Moisture Dynamics Estimation Model

Abstract The long‐term containment of high‐level radioactive waste in geological disposal repositories relies on Engineered Barrier Systems (EBS), with bentonite clay emerging as a candidate material due to its unique properties. Understanding moisture dynamics within bentonite buffers is crucial for EBS performance, as it directly influences the material's swelling capacity, thermal and hydraulic conductivity, mechanical properties, and long‐term evolution under complex thermal‐hydrological‐mechanical (THM) processes. This study develops an advanced Electrical Resistivity Tomography (ERT)‐based framework to quantitatively monitor moisture dynamics under THM conditions. Our framework extends the Waxman‐Smits model to incorporate the coupled effects of temperature, water content, fluid chemistry, and mechanical changes on bentonite's electrical properties. Utilizing HotBENT‐Lab data from our companion paper, which includes electrical conductivity, CT density, and thermocouple measurements, this study offers a novel methodological framework bridging different scales of the model. Our results show that the extended model can estimate water content from ERT data, capturing spatial and temporal variations in moisture distribution within bentonite columns. However, the model tends to overestimate water content compared to CT density‐derived measurements. We address this discrepancy by incorporating a simplified swelling effect model, which improves agreement between ERT and CT density‐based water content estimates. We also discuss model limitations, including simplified treatment of swelling and micropore effects, and propose a conceptual framework for transitioning from laboratory to field applications, addressing challenges such as parameter scalability, field validation methods, and integration of diverse data sources. This ERT‐based framework can potentially advance real‐world moisture monitoring of bentonite‐based EBS in nuclear waste repositories. Plain Language Summary Safely containing high‐level radioactive waste depends on barriers made from materials like bentonite clay, which is effective because it swells and seals in the waste. To ensure these barriers work well over time, it's important to understand how moisture moves through the clay. Our study developed a new method using ERT to monitor moisture levels in bentonite under conditions that mimic those in actual storage sites, including changes in temperature, water content, and mechanical stress. This study improved an existing model to better account for how these factors affect the clay, allowing us to create more accurate moisture maps. Initially, the proposed model overestimated the amount of water in the clay, but its accuracy was improved by factoring in how the clay swells when wet. This study also identified some limitations of the model and suggested ways to adapt it for use in real‐world waste storage sites. This new approach could lead to better monitoring and safety checks for nuclear waste storage systems, helping to ensure long‐term containment. Key Points This work develops an ERT‐based framework extending the Waxman‐Smits model to monitor bentonite moisture dynamics during coupled THM processes The extended model accurately estimates water content from Electrical Resistivity Tomography data, incorporating swelling effects to improve precision This work proposes a conceptual framework for transitioning from laboratory to field applications, advancing EBS monitoring in nuclear waste repositories

Chen, Hang↗