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 847 records · Page 47

Machine Learning for LBNF Beam Diagnostics

This paper focuses on developing a machine learning model for predicting initial beam parameters for the Long Baseline Neutrino Facility (LBNF) beamline using downstream muon monitor data. Parameters such as proton beam position on target, sigma on target, focusing horn current, and focusing horn tilt are parameters we anticipate to be predictable based on the muon monitors. Uncertainty in initial beam condition measurements are a major contributor to uncertainty in downstream flux, and over operation time beam misalignment can occur [1]. A machine learning model has promise to detect anomalies along the beamline based on discrepancies between predicted configurations and measured configurations, and thus can expedite error detection and handling. A PyTorch neural network is defined, trained, and tested. The developed model currently does not provide reliable predictions, with the lowest loss being 0.09.. Further steps to improve the model’s accuracy are discussed, as well as future plans to detect anomalous beam conditions using a digital twin.

O'Brien, Bridget [Fermilab]↗

Windblown sand on Venus - Preliminary results of laboratory simulations

Small particles and winds of sufficient strength to move them have been detected from Venera and Pioneer-Venus data and suggest the existence of aeolian processes on Venus. The Venus wind tunnel (VWT) was fabricated in order to investigate the behavior of windblown particles in a simulated Venusian environment. Preliminary results show that sand-size material is readily entrained at the wind speeds detected on Venus and that saltating grains achieve velocities closely matching those of the wind. Measurements of saltation threshold and particle flux for various particle sizes have been compared with theoretical models which were developed by extrapolation of findings from Martian and terrestrial simulations. Results are in general agreement with theory, although certain discrepancies are apparent which may be attributed to experimental and/or theoretical-modeling procedures. Present findings enable a better understanding of Venusian surface processes and suggest that aeolian processes are important in the geological evolution of Venus.

Greeley, R.↗

Toward an improved nonlocal thermodynamic equilibrium model for more predictive simulations of ignition scale hohlraums

Recently, nonlocal thermodynamic equilibrium (NLTE) modeling has been identified as the primary reason for discrepant predictions of the peak neutron production time in indirectly driven inertial confinement fusion (ICF) platforms. It has also been observed that predictions of collisional excitation rates differ by as much as 50% from measurements. Theoretical uncertainties in dielectronic recombination rates have also been posited as possibly contributing to errors in NLTE predictions. This work examines the impact of multipliers on collisional excitation and dielectronic recombination rates on simulations of a directly driven gold sphere and an indirect drive ICF implosion. It is found that multipliers on the collisional excitation rates have a strong impact on radiant intensity and electron temperature and a weaker impact on ionization state, whereas multipliers on dielectronic recombinations rates strongly impact ionization state with a smaller impact on radiant intensity and electron temperature. A self-consistent NLTE model which places multipliers on differing transitions, as motivated by experimental measurements and more detailed atomic physics predictions, improves agreement but does not completely eliminate discrepancies with measurements of the radiant intensity within the 2–4 keV spectral range.

Farmer, W. A. [Lawrence Livermore National Laborat↗

Measurement of the production cross section of prompt $\Xi ^0_{\textrm{c}}$ baryons in p–Pb collisions at $\sqrt{s_{{\textrm{NN}}}}=5.02$ TeV

The transverse momentum (p T ) differential production cross section of the promptly produced charm-strange baryon $\Xi ^0_{\textrm{c}}$ (and its charge conjugate $\Xi ^0_{\textrm{c}}$) is measured at midrapidity via its hadronic decay into π + $\Xi$ – in p–Pb collisions at a centre-of-mass energy per nucleon–nucleon collision $\sqrt{s_{{\textrm{NN}}}}=5.02$ TeV with the ALICE detector at the LHC. The $\Xi ^0_{\textrm{c}}$ nuclear modification factor (R pPb ), calculated from the cross sections in pp and p–Pb collisions, is presented and compared with the R pPb of Λ$^{+}_{c}$ baryons. The ratios between the p T -differential production cross section of $\Xi ^0_{\textrm{c}}$ baryons and those of D 0 mesons and Λ$^{+}_{c}$ baryons are also reported and compared with results at forward and backward rapidity from the LHCb Collaboration. The measurements of the production cross section of prompt $\Xi ^0_{\textrm{c}}$ baryons are compared with a model based on perturbative QCD calculations of charm-quark production cross sections, which includes only cold nuclear matter effects in p–Pb collisions, and underestimates the measurement by a factor of about 50. This discrepancy is reduced when the data is compared with a model that includes string formation beyond leading-colour approximation or in which hadronisation is implemented via quark coalescence. The p T -integrated cross section of prompt $\Xi ^0_{\textrm{c}}$-baryon production at midrapidity extrapolated down to p T = 0 is also reported. These measurements offer insights and constraints for theoretical calculations of the hadronisation process. Additionally, they provide inputs for the calculation of the charm production cross section in p–Pb collisions at midrapidity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The use of ground magnetograms to time the onset of magnetospheric substorms

A critical review of magnetospheric substorm timing techniques is presented. Attention is given to substorm timing with auroral zone magnetograms, substorm onset timing with midlatitude magnetograms, multiple onset substorms, and a phenomenological model for multiple onset substorms. Emphasis is placed on some pitfalls that have led some investigators to err in substorm onset determinations. In particular, multiple onset substorms appear to account for many of the reported discrepancies in substorm timing. It is suggested that present idealized models of substorm morphology must be modified to include the multiple onset substorms.

Mcpherron, R. L.↗

Benchmarking boron carbide equation of state using computation and experiment

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 – 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are ≲ 5 % near the compression maximum at 1– 2 × 10 6 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 – 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Microscopic origin of temperature-dependent magnetism in spin-orbit-coupled transition metal compounds

A few 4 d and 5 d transition metal compounds with various electron fillings were recently found to exhibit magnetic susceptibilities χ and magnetic moments that deviate from the well-established Kotani model. This model has been considered for decades to be the canonical expression for descriing the temperature dependence of magnetism in systems with nonnegligible spin-orbit coupling effects. In this paper, we uncover the origin of such discrepancies and determine the applicability and limitations of the Kotani model by calculating the temperature dependence of the magnetic moments of a series of 4 d (Ru-based) and 5 d (W-based) systems at different electron fillings. For this purpose, we perform exact diagonalization of -derived relativistic multiorbital Hubbard models on finite clusters and compute their magnetic susceptibilities. Comparison with experimentally measured magnetic properties indicates that contributions such as a temperature-independent χ 0 background, crystal field effects, Coulomb and Hund's couplings, and intersite interactions—not included in the Kotani model—are especially crucial for correctly describing the temperature dependence of χ and magnetic moments at various electron fillings in these systems. Based on our results, we propose a generalized approach beyond the Kotani model to accurately describe their magnetism. Published by the American Physical Society 2025

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Classical bump Cepheids - Reconciliation of theory with observations

New models of classical bump Cepheids lying along the Cepheid ridge line in the H-R diagram are computed to obtain a more accurate set of predicted surface velocity curves. Acceptable Cepheid models are found to be 0.5 mag brighter at a fixed stellar mass than standard evolutionary models. All the previously obtained manifestations of the Cepheid mass discrepancy, except for the anomalous double-mode Cepheid masses, are resolved, at least in a statistical sense, by accepting the implied luminosity increase. The need for the luminosity increase is independently inferred from standard evolutionary models of intermediate-mass and high-mass giant and supergiant stars from their location in the observational H-R diagram.

Carson, T. Richard↗

Wind gust models derived from field data

Wind data measured during a field experiment were used to verify the analytical model of wind gusts. Good coincidence was observed; the only discrepancy occurred for the azimuth error in the front and back winds, where the simulated errors were smaller than the measured ones. This happened because of the assumption of the spatial coherence of the wind gust model, which generated a symmetric antenna load and, in consequence, a low azimuth servo error. This result indicates a need for upgrading the wind gust model to a spatially incoherent one that will reflect the real gusts in a more accurate manner. In order to design a controller with wind disturbance rejection properties, the wind disturbance should be known at the input to the antenna rate loop model. The second task, therefore, consists of developing a digital filter that simulates the wind gusts at the antenna rate input. This filter matches the spectrum of the measured servo errors. In this scenario, the wind gusts are generated by introducing white noise to the filter input.

Gawronski, W.↗

Fluidization of Group A Glass Particles: Experiments and Preliminary Validation

This work is aimed at providing reliable high-quality data from fluidization experiments. Geldart Group A glass beads are used as bed material. The fluidizing medium is air at atmospheric pressure which enters the system at 40% relative humidity. The air inlet velocity is varied over a broad range from 2 U mf to 26 U mf , where the minimum fluidization velocity, U mf , is 0.608 cm/s. There is not a significant change in the mean values of differential pressure; however, there are noticeable trends in their standard deviation values. Preliminary validation studies are performed using the open-source software Multiphase Flow with Interphase eXchanges (MFiX) developed by the National Energy Technology Laboratory (NETL). The results from Eulerian-Eulerian (MFiX-TFM, two-fluid model) and Eulerian-Lagrangian (MFiX-PIC, particle-in-cell) formulations are compared. The work presented in this report highlights the applicability of MFiX-PIC for large-scale applications, and reveals that PIC simulation offers a suitable alternative to TFM simulation when effects due to multiple components or polydispersity are significant. PIC modeling provides a readily extendible framework for such systems. However, PIC models do rely on an empirical closure for inter-particle stress. Some solution discrepancies were evident in MFiX-PIC predictions which require further analysis. As such, the study underlines the need for a systematic approach to quantify the sources of numerical uncertainty within the PIC model. However, even with known empiricism, MFiX-PIC has demonstrated considerable potential in analyzing gas-solid systems. The results obtained are encouraging and warrant further investigation to improve the predictive capability of MFiX-PIC.

36 MATERIALS SCIENCE↗

Development of a Discrepancy Checker for the Digital Twin in a Supervisory Control System for a Thermal Energy Delivery System

Defined as a virtual representation of a physical object, process, or service, and used to support real-world decision-making, a digital twin (DT) can be utilized to combine classical and novel frameworks in sensors, state predictions, and multi-input/multi-output systems, and to enable optimal autonomous operations. However, a DT’s usefulness largely depends on its ability to adequately mirror the state of its physical counterpart, and this adequacy should be reflected by the level of uncertainty in the underlying simulation models when estimating and predicting quantities of interest (QOIs). Moreover, simulation models in a DT may involve multiple fidelities of representations—ranging from physics-based models to data-driven ones—but classical uncertainty quantification (UQ) methods struggle to handle numerous uncertainty sources, nor are they designed for real-time applications. This work presents a UQ-based discrepancy checking and diagnosis tool for a DT-based supervisory control system applied to a thermal energy delivery system (TEDS) at Idaho National Laboratory. The discrepancy checker was developed using metadata from an automated DT development process, and these metadata included different combinations of physical model forms and model parameters, training data and hyperparameters for surrogate models, and design parameters for supervisory control systems. Next, correlations between the uncertainty results and the metadata were established and then applied to the DT operations. The discrepancy checker evaluates the discrepancies between model predictions from virtual and sensor measurements and backtraces them to the corresponding major sources of uncertainty. The discrepancy checker showed reasonable performance in detecting discrepancies and diagnosing sources of uncertainty in testing scenarios.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Salt and ion transport in a series of crosslinked AMPS/PEGDA hydrogel membranes

Salt permeability and ionic conductivity are critical properties for membranes in water purification and energy applications. Both properties depend on individual ion sorption and diffusion coefficients, which are significantly influenced by polymer chemical and physical parameters such as fixed charge concentration and membrane water content. However, systematic studies connecting polymer structure to ion transport properties are still lacking. In this study, a series of uncharged and charged membranes were synthesized using poly(ethylene glycol) diacrylate (PEGDA) as a cross-linker and 2-acrylamido-2-methyl-1-propanesulfonic acid (AMPS) as a charged monomer. Membrane fixed charge concentration and water uptake were systematically varied by adjusting AMPS content in the pre-polymerization mixture. Salt sorption and permeability coefficients and ionic conductivity of these membranes were measured as a function of NaCl solution concentration (0.01–1 M). Combining the solution-diffusion model and Nernst-Planck equation, individual ion diffusion coefficients were calculated. Experimental Na + diffusion coefficients for all materials were well described by the Mackie and Meares tortuosity model, highlighting the strong influence of water content on ion diffusivity in both uncharged and charged polymers. Model predictions for Cl – diffusion coefficients agree reasonably well with experimental values, with some deviations occurring in more highly charged membranes. Finally, this discrepancy might result from interactions not captured by the Mackie and Meares model (e.g., fixed charge-ion interactions).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Regional and Global Climate Response to Anthropogenic SO2 Emissions from China in Three Climate Models

We use the HadGEM3-GA4, CESM1, and GISS ModelE2 climate models to investigate the global and regional aerosol burden, radiative flux, and surface temperature responses to removing anthropogenic sulfur dioxide (SO2) emissions from China. We find that the models differ by up to a factor of 6 in the simulated change in aerosol optical depth (AOD) and shortwave radiative flux over China that results from reduced sulfate aerosol, leading to a large range of magnitudes in the regional and global temperature responses. Two of the three models simulate a near-ubiquitous hemispheric warming due to the regional SO2 removal, with similarities in the local and remote pattern of response, but overall with a substantially different magnitude. The third model simulates almost no significant temperature response. We attribute the discrepancies in the response to a combination of substantial differences in the chemical conversion of SO2 to sulfate, translation of sulfate mass into AOD, cloud radiative interactions, and differences in the radiative forcing efficiency of sulfate aerosol in the models. The model with the strongest response (HadGEM3-GA4) compares best with observations of AOD regionally, however the other two models compare similarly (albeit poorly) and still disagree substantially in their simulated climate response, indicating that total AOD observations are far from sufficient to determine which model response is more plausible. Our results highlight that there remains a large uncertainty in the representation of both aerosol chemistry as well as direct and indirect aerosol radiative effects in current climate models, and reinforces that caution must be applied when interpreting the results of modelling studies of aerosol influences on climate. Model studies that implicate aerosols in climate responses should ideally explore a range of radiative forcing strengths representative of this uncertainty, in addition to thoroughly evaluating the models used against observations.

sulfates↗

Validation of seismic velocity models in southern California with full-waveform simulations

SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.

58 GEOSCIENCES↗

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Graphical interface between the CIRSSE testbed and CimStation software with MCS/CTOS

This research is concerned with developing a graphical simulation of the testbed at the Center for Intelligent Robotic Systems for Space Exploration (CIRSSE) and the interface which allows for communication between the two. Such an interface is useful in telerobotic operations, and as a functional interaction tool for testbed users. Creating a simulated model of a real world system, generates inevitable calibration discrepancies between them. This thesis gives a brief overview of the work done to date in the area of workcell representation and communication, describes the development of the CIRSSE interface, and gives a direction for future work in the area of system calibration. The CimStation software used for development of this interface, is a highly versatile robotic workcell simulation package which has been programmed for this application with a scale graphical model of the testbed, and supporting interface menu code. A need for this tool has been identified for the reasons of path previewing, as a window on teleoperation and for calibration of simulated vs. real world models. The interface allows information (i.e., joint angles) generated by CimStation to be sent as motion goal positions to the testbed robots. An option of the interface has been established such that joint angle information generated by supporting testbed algorithms (i.e., TG, collision avoidance) can be piped through CimStation as a visual preview of the path.

Hron, Anna B.↗

Spark plasma sintered, MoNbTi-based multi-principal element alloys with Cr, V, and Zr

MoNbTi, MoNbTiZr, CrMoNbTiZr, and MoNbTiVZr multi-principal element alloys (MPEAs) were fabricated via spark plasma sintering (SPS) and investigated for use in high-strength applications. The fabrication method by SPS and powder metallurgy differs from those presented in the prior literature, where most MPEAs are fabricated using arc melting (AM) methods. Here cryogenic milling was used to maximize potential defect sinks (grain boundaries) for radiation resistance and to increase ultimate tensile strength through the Hall-Petch effect. SPS was chosen for consolidation in order to maintain a fine-grained structure during densification. Each alloy was characterized using x-ray diffraction and scanning electron microscopy for phase identification and compositional homogeneity. The base ternary alloy MoNbTi presented a predominantly single BCC system with minor cubic phases. The introduction of additional alloying elements—Zr, V, and Cr—heightened the phase complexity and increased the fractions of a secondary BCC phase and an HCP phase from Zr. The addition of Cr induced a larger fraction of the Laves phase to form. The addition of V caused the precipitation of small Mo inclusions. Thermodynamic analysis was performed to understand the separation of phases in each alloy. Discrepancies among the phase predictions generated by thermodynamic models, phases previously presented in the literature, and the characterization results suggest that MPEA fabrication methods, especially solid-state methods, require significant investigation to ensure that alloys can remain stable throughout their anticipated service lifetimes.

36 MATERIALS SCIENCE↗

Phonon relaxation time reconstruction from transient thermal grating experiments and comparison with density functional theory predictions

Phonon relaxation time and free path distributions are reconstructed from experimental measurements on a two-dimensional transient thermal grating and compared with density functional theory (DFT) results for silicon. The reconstruction is performed using the inverse problem formulation of Forghani et al. [Phys. Rev. B 94, 155439 (2016)]. The discrepancies observed between reconstructed and DFT results are analyzed in terms of the ability of each set of data to reproduce the experimental temperature relaxation profiles; the reconstructed data are found to predict temperature profiles in closer agreement with the experimental data than the DFT results, possibly due to discrepancies between the actual material and the idealized model studied in the DFT calculations. The reconstructed phonon properties accurately predict temperature relaxation profiles at grating length scales smaller than those spanned by the experimental data. Finally, this is a very important feature since in a variety of experimental setups, including the one providing the data in the present study, measurements are not available at all scales spanned by the material free paths.

42 ENGINEERING↗