Search NASA⌕ Search

SEARCH · Search NASA

Results for “convergence diagnostics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Revisiting Source Convergence Diagnostics in the KENO Monte Carlo Neutron Transport Codes [Abstract]

Monte Carlo criticality transport codes, which rely on the power iteration procedure, are a fundamental tool for nuclear criticality safety practitioners in assessing the neutron multiplication factor (k eff ) for problems involving fissile material. In these calculations, ensuring the convergence of both the fission source distributions and the k eff estimate for accurate results is crucial. However, a converged k eff estimate does not necessarily mean the fission source distribution is also converged because the fission source and flux distribution may continue to evolve even after k eff convergence. Therefore, most Monte Carlo transport criticality codes now offer various diagnostic tests to assess fission source convergence in addition to the k eff convergence by analyzing the trends of these quantities over multiple generations.

AZURE↗
DOE OSTI2022

Challenges in Markov Chain Monte Carlo for Bayesian Neural Networks

Markov chain Monte Carlo (MCMC) methods have not been broadly adopted in Bayesian neural networks (BNNs). This paper initially reviews the main challenges in sampling from the parameter posterior of a neural network via MCMC. Such challenges culminate to lack of convergence to the parameter posterior. Nevertheless, this paper shows that a nonconverged Markov chain, generated via MCMC sampling from the parameter space of a neural network, can yield via Bayesian marginalization a valuable posterior predictive distribution of the output of the neural network. Further, classification examples based on multilayer perceptrons showcase highly accurate posterior predictive distributions. The postulate of limited scope for MCMC developments in BNNs is partially valid; an asymptotically exact parameter posterior seems less plausible, yet an accurate posterior predictive distribution is a tenable research avenue.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

A test suite for magnetohydrodynamical simulations

A collection is presented of MHD problems that will provide researchers who have interests in modeling MHD phenomena with a battery of tests (a 'test suite') for calibrating their numerical algorithms. Use of these tests will provide a common reference for comparison of different numerical MHD algorithms. The test suite includes both 1D and 2D problems. Taken together, these problems test the abilities of a numerical method to propagate accurately all of the MHD wave families in moving and stationary media, to capture MHD shocks and contact discontinuities, and to model accurately Lorentz force terms in multiple dimensions. An example solution of each test problem is presented. Diagnostic convergence-testing procedures that provide a quantitative evaluation of MHD algorithms are demonstrated.

Stone, James M.↗
DOE OSTI2020

Hydro-scaling of direct-drive cylindrical implosions at the OMEGA and the National Ignition Facility

Deceleration-phase Rayleigh–Taylor instability (RTI) growth during inertial confinement fusion capsule implosions significantly affects the performance as it mixes cold ablator material into the fuel. Precise measurements of such instability growth are essential for both validating the existing simulation codes and improving our predictive capability. RTI measurements on the inner surface of a spherical shell are limited and are often inferred indirectly at limited convergence. In contrast, cylindrical implosions allow for direct diagnostic access to the converging interface by imaging down the cylinder axis while retaining the effects of convergence. We have performed direct-drive cylindrical implosion experiments at both the OMEGA and the NIF laser facilities using scaled targets. RTI growth is demonstrated to be scale-invariant between the cylindrical targets at OMEGA and similar targets at the NIF, which are scaled up by a factor of three in the radial dimension. Single-mode (m = 0) instability growth factors of ~14 are measured at a convergence ratio (CR) ~ 2.5 with nearly identical mode growth at both scales. The measurements are in agreement with xRAGE radiation-hydrodynamics simulations. In addition, we have developed the Bayesian-inference-engine method to account for the variations in the target alignment, magnification, and the parallax effect in the measurement, allowing a more precise comparison between the experimental data and the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Atmospheric teleconnection processes linking winter air stagnation and haze extremes in China with regional Arctic sea ice decline

Abstract. Recent studies suggested significant impacts of boreal cryosphere changes on wintertime air stagnation and haze pollution extremes in China. However, the underlying mechanisms of such a teleconnection relationship remains unclear. Here we use the Whole Atmosphere Community Climate Model (WACCM) to investigate dynamic processes leading to atmospheric circulation and air stagnation responses to Arctic sea ice changes. We conduct four climate sensitivity experiments by perturbing sea ice concentrations (SIC) and corresponding sea surface temperature (SST) in autumn and early winter over the whole Arctic and three subregions in the climate model. The results indicate distinct responses in circulation patterns and regional ventilation to the region-specific Arctic changes, with the largest increase of both the probability (by 132 %) and the intensity (by 30 %) of monthly air stagnation extremes being found in the experiment driven by SIC and SST changes over the Pacific sector of the Arctic (the East Siberian and Chukchi seas). The increased air stagnation extremes are mainly driven by an amplified planetary-scale atmospheric teleconnection pattern that resembles the negative phase of the Eurasian (EU) pattern. Dynamical diagnostics suggest that convergence of transient eddy forcing in the vicinity of Scandinavia in winter is largely responsible for the amplification of the teleconnection pattern. Transient eddy vorticity fluxes dominate the transient eddy forcing and produce a barotropic anticyclonic anomaly near Scandinavia and wave train propagation across Eurasia to the downstream regions in East Asia. The piecewise potential vorticity inversion analysis reveals that this long-range atmospheric teleconnection of Arctic origin takes place primarily via the middle and upper troposphere. The anomalous ridge over East Asia in the middle and upper troposphere worsens regional ventilation conditions by weakening monsoon northwesterlies and enhancing temperature inversions near the surface, leading to more and stronger air stagnation and pollution extremes over eastern China in winter. Ensemble projections based on state-of-the-art climate models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) corroborate this teleconnection relationship between high-latitude environmental changes and midlatitude weather extremes, though the tendency and magnitude vary considerably among each participating model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Santa Barbara Binary-disk Code Comparison

We have performed numerical calculations of a binary interacting with a gas disk, using 11 different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, 0.2(a 2 Ω B /v) is the kinematic viscosity. The largest discrepancies between codes resulted from the dimensionality of the setup (3D vs. 2D disks). We find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Water Mass Transformation Budgets in Finite‐Volume Generalized Vertical Coordinate Ocean Models

Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Double cylinder target design for study of hydrodynamic instabilities in multi-shell ICF

Cylindrical implosions are used to study hydrodynamic instability growth for inertial confinement fusion (ICF) applications, as the cylindrical geometry allows for easier diagnostic access while retaining convergence effects. In this work, we use the established cylindrical implosion platform [Palaniyappan et al., Phys. Plasmas 27, 042708 (2020)] to inform the double shell ICF campaign [Montgomery et al., Phys. Plasmas 25, 092706 (2018)]. We present a design for a double cylindrical target as an analogue to the double shell ICF capsule in order to study hydrodynamic instability growth on the high-Z inner shell. Our design work is done with two-dimensional (2D) Eulerian radiation-hydrodynamics simulations, considering the axial uniformity of the implosion and feasibility of measuring the instability growth of pre-seeded single mode sinusoidal perturbations. We discuss in depth the design for a target to be directly driven at the OMEGA laser facility [Boehly et al., Opt. Commun. 133, 495 (1997)]. We evaluate the design for axial implosion symmetry and visibility of instability growth using synthetic radiographs constructed from the simulations, as the instability growth on the inner cylinder is experimentally measured using x-ray radiography of the implosion. We find that the seeded perturbation growth on the inner cylinder should be visible in an experiment, even with axial implosion asymmetry and preheat. We compare our 2D simulations with linear theory predictions for perturbation growth and show that a cylinder with lower azimuthal mode number (mode-20) perturbations compares more favorably with linear theory, while a cylinder with higher azimuthal mode number (mode-40) perturbations at the same starting amplitude saturates and is over-predicted by linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

SE-FIT

The mathematical theory of capillary surfaces has developed steadily over the centuries, but it was not until the last few decades that new technologies have put a more urgent demand on a substantially more qualitative and quantitative understanding of phenomena relating to capillarity in general. So far, the new theory development successfully predicts the behavior of capillary surfaces for special cases. However, an efficient quantitative mathematical prediction of capillary phenomena related to the shape and stability of geometrically complex equilibrium capillary surfaces remains a significant challenge. As one of many numerical tools, the open-source Surface Evolver (SE) algorithm has played an important role over the last two decades. The current effort was undertaken to provide a front-end to enhance the accessibility of SE for the purposes of design and analysis. Like SE, the new code is open-source and will remain under development for the foreseeable future. The ultimate goal of the current Surface Evolver Fluid Interface Tool (SEFIT) development is to build a fully integrated front-end with a set of graphical user interface (GUI) elements. Such a front-end enables the access to functionalities that are developed along with the GUIs to deal with pre-processing, convergence computation operation, and post-processing. In other words, SE-FIT is not just a GUI front-end, but an integrated environment that can perform sophisticated computational tasks, e.g. importing industry standard file formats and employing parameter sweep functions, which are both lacking in SE, and require minimal interaction by the user. These functions are created using a mixture of Visual Basic and the SE script language. These form the foundation for a high-performance front-end that substantially simplifies use without sacrificing the proven capabilities of SE. The real power of SE-FIT lies in its automated pre-processing, pre-defined geometries, convergence computation operation, computational diagnostic tools, and crash-handling capabilities to sustain extensive computations. SE-FIT performance is enabled by its so-called file-layer mechanism. During the early stages of SE-FIT development, it became necessary to modify the original SE code to enable capabilities required for an enhanced and synchronized communication. To this end, a file-layer was created that serves as a command buffer to ensure a continuous and sequential execution of commands sent from the front-end to SE. It also establishes a proper means for handling crashes. The file layer logs input commands and SE output; it also supports user interruption requests, back and forward operation (i.e. undo and redo), and others. It especially enables the batch mode computation of a series of equilibrium surfaces and the searching of critical parameter values in studying the stability of capillary surfaces. In this way, the modified SE significantly extends the capabilities of the original SE.

Chen, Yongkang↗
DOE OSTI2022

Ab initio nucleon-nucleus elastic scattering with chiral effective field theory uncertainties

Effective interactions for nucleon-nucleus elastic scattering from first principles require the use of the same nucleon-nucleon interaction in the structure and reaction calculations and a consistent treatment of the relevant operators at each order. Systematic investigations of the effect of truncation uncertainties of chiral nucleon-nucleon ( N N ) forces have been carried out for scattering observables in the two- and three-nucleons system as well as for bound-state properties of light nuclei. Here we extend this type of study to proton and neutron elastic scattering for 16 O and 12 C. Using the frameworks of the spectator expansion of multiple scattering theory as well as the no-core shell model, we employ one specific chiral interaction from the LENPIC collaboration and consistently calculate the leading-order effective nucleon-nucleus interaction up to the third chiral order, from which we extract elastic scattering observables. Then we apply pointwise as well as correlated uncertainty quantification for the estimation of the chiral truncation error. We calculate and analyze proton elastic scattering observables for 16 O and 12 C as well as neutron elastic scattering observables for 12 C between 65 and 185 MeV projectile kinetic energy. We find in all cases qualitatively similar results for the chiral truncation uncertainties as in few-body systems and assess them using similar diagnostic tools. The order-by-order convergence of the scattering observables for 16 O and 12 C is very reasonable around 100 MeV, while for higher energies the chiral expansion parameter becomes too large for convergence. Comparing proton and neutron differential cross sections reveals that their description is comparable up to around 100 MeV. Here, we also find a nearly perfect correlation between the differential cross section for neutron scattering and the N N Wolfenstein amplitudes for small momentum transfers. The diagnostic tools for studying order-by-order convergence in observables in few-body systems can be employed for observables in nucleon-nucleus scattering with only minor modifications provided the momentum scale in the problems is not too large. We also find that the chiral N N interaction on which our study is based gives a very good description of differential cross sections for 16 O and 12 C as low as 65-MeV projectile energy, particularly in the forward direction. In addition, the very forward direction of the neutron differential cross section mirrors the behavior of the N N interaction amazingly well.

6 ≤ A ≤ 19↗
NASA NTRS1989

Global atmospheric moisture variability

Research efforts during FY-88 have focused on completion of several projects relating to analysis of FGGE data during SOP-1 and on expanded studies of global atmospheric moisture. In particular, a revised paper on the relationship between diabatic heating and baroclinicity in the South Pacific Convergence Zone (SPCZ) was submitted. A summary of completed studies on diagnostic convective parameterization was presented at the Satellite Meteorology and Oceanography Convergence last February. These investigations of diabatic heating in the SPCZ have demonstrated the requirement for a more quantitative description of atmospheric moisture. As a result, efforts were directed toward use of passive remote microwave measurements from the Nimbus-7 SMMR and the DOD's Special Sensor Microwave Imager (SSMI/I) as critical sources of moisture data. Activities this year are summarized.

Robertson, Franklin R.↗
DOE OSTI2020

Azimuthal Drive Asymmetry in Inertial Confinement Fusion Implosions on the National Ignition Facility

Data from nuclear diagnostics present correlated signatures of azimuthal implosion asymmetry in recent indirect-drive inertial confinement fusion (ICF) implosion campaigns performed at the National Ignition Facility (NIF). The mean hot-spot velocity, inferred from the Doppler shift of 14 MeV neutrons produced by deuterium-tritium (DT) fusion, is systematically directed toward one azimuthal half of the NIF target chamber, centered on φ ≈ 7 0 ° . Areal density ( ρ R ) asymmetry of the converged DT fuel, inferred from nuclear activation diagnostics, presents a minimum ρ R in the same direction as the hot-spot velocity and with Δ ρ R amplitude correlated with velocity magnitude. These two correlated observations, which are seen in all recent campaigns with cryogenic layers of DT fuel, are an understood signature of asymmetry in the fuel convergence, implying a systematic azimuthal drive asymmetry across a wide range of shot and target configurations. The direction of the implied radiation asymmetry is observed to cluster toward the hohlraum diagnostic windows. Furthermore, this low-mode asymmetry degrades hot-spot conditions at peak convergence and limits implosion performance and yield

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Merging flux-variance with surface renewal methods in the roughness sublayer and the atmospheric surface layer

In this article, Two micrometeorological methods utilizing high-frequency sampled air temperature were tested against eddy covariance (EC) sensible heat flux (𝐻) measurements at three sites representing agricultural, agro-forestry, and forestry systems. The two methods cover conventional and newly proposed forms of the flux-variance (FV) and surface renewal (SR) schemes of differing complexities. The sites represent measurements in surface, roughness, and roughness to surface transitional layers. Regression analyzes against EC show that the most reliable FV and SR forms estimate 𝐻 with slopes within ±10% from unity and coefficient of determination 𝑅2 > 0.9 across all the three sites. The best performance of both FV and SR was found at the agricultural site with measurements well within the surface layer, while the worst was found for the tall forest with measurements within the roughness sublayer where its thickness needed to be additionally estimated. The main variable driving 𝐻 in FV is the temperature variance, whereas in SR, it is the geometry of ramp-like structures. Since these structures are also responsible for most of the temperature variance, a novel FV-SR approach emerging from combining the methods is proposed and evaluated against EC measurements and conventional FV and SR schemes. The proposed FV-SR approach requiring only a single fast response thermocouple is potentially independent of calibration and ameliorates some of the theoretical objections that arise when combining ramp statistics with similarity arguments. The combination of methods also provides new insights into the contribution of coherent structures to the temperature variance and its dependence on atmospheric stratification. Other potential utility of the new method is to include it in multi-tool assessments of surface energy fluxes, since a convergence or divergence of the results has a high diagnostic value.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Assessing correlated truncation errors in modern nucleon-nucleon potentials

We test the BUQEYE model of correlated effective field theory (EFT) truncation errors on Reinert, Krebs, and Epelbaum's semilocal momentum-space implementation of the chiral EFT (𝜒⁢EFT ) expansion of the nucleon-nucleon (NN) potential. This Bayesian model hypothesizes that dimensionless coefficient functions extracted from the order-by-order corrections to NN observables can be treated as draws from a Gaussian process (GP). We combine a variety of graphical and statistical diagnostics to assess when predicted observables have a 𝜒⁢EFT convergence pattern consistent with the hypothesized GP statistical model. Our conclusions are that, first, the BUQEYE model is generally applicable to the potential investigated here, which enables statistically principled estimates of the impact of higher EFT orders on observables. Second, parameters defining the extracted coefficients such as the expansion parameter 𝑄 must be well chosen for the coefficients to exhibit a regular convergence pattern—a property we exploit to obtain posterior distributions for such quantities. Third, the assumption of GP stationarity across lab energy and scattering angle is not generally met; this necessitates adjustments in future work. We provide a workflow and interpretive guide for our analysis framework, and show what can be inferred about probability distributions for 𝑄, the EFT breakdown scale Λ 𝑏 , the scale associated with soft physics in the 𝜒⁢EFT potential 𝑚 eff , and the GP hyperparameters. All our results can be reproduced using a publicly available Jupyter notebook, which can be straightforwardly modified to analyze other 𝜒⁢EFT NN potentials.

Bayesian methods↗
DOE OSTI2021

Dissecting Anvil Cloud Response to Sea Surface Warming

Abstract We derive an anvil cloud diagnostic from the continuity equation of cloud ice and apply it to the output of convection‐permitting Energy Exascale Earth System Model (E3SM) simulations run in radiative‐convective equilibrium mode. This diagnostic shows that anvil cloud fraction can be reliably diagnosed as a product of cloud detrainment and lifetime. Detrainment is found to be approximated well by a product of clear sky convergence and cloud ice mixing ratio, while cloud lifetime is dominated by sedimentation. Taken together, this diagnostic expresses anvil cloud fraction as a function of five physically measurable quantities. Of these, clear‐sky convergence changes drive the anvil cloud reduction with warming while an increase in cloud ice mixing ratio buffers the decrease. Accordingly, this study provides a theoretical foundation upon which the Stability‐Iris hypothesis can be tested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗