Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Integration”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Inference of bipolar neutrino flavor oscillations near a core-collapse supernova based on multiple measurements at Earth

Neutrinos in compact-object environments, such as core-collapse supernovae, can experience various kinds of collective effects in flavor space, engendered by neutrino-neutrino interactions. These include “bipolar” collective oscillations, which are exhibited by neutrino ensembles where different flavors dominate at different energies. Considering the importance of neutrinos in the dynamics and nucleosynthesis in these environments, it is desirable to ascertain whether an Earth-based detection could contain signatures of bipolar oscillations that occurred within a supernova envelope. To that end, we, in this study, continue examining a cost-function formulation of statistical data assimilation (SDA) to infer solutions to a small-scale model of neutrino flavor transformation. SDA is an inference paradigm designed to optimize a model with sparse data. Our model consists of two monoenergetic neutrino beams with different energies emanating from a source and coherently interacting with each other and with a matter background, with radially varying interaction strengths. We attempt to infer flavor transformation histories of these beams using simulated measurements of the flavor content at locations “in vacuum” (that is, far from the source), which could in principle correspond to Earth-based detectors. Within the scope of this small-scale model, we found that: (i) based on such measurements, the SDA procedure is able to infer whether bipolar oscillations had occurred within the protoneutron star envelope, and (ii) if the measurements sample the full amplitude of the neutrino oscillations in vacuum, then the amplitude of the prior bipolar oscillations is well predicted. This result intimates that the inference paradigm can well complement numerical integration codes, via its ability to infer flavor evolution at physically inaccessible locations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Impact of LHC precision measurements of inclusive jet and dijet production on the CTEQ-TEA global PDF fit

In this study, we investigate the impact of new Large Hadron Collider (LHC) inclusive jet and dijet measurements on parton distribution functions (PDFs) that describe the proton structure, with a particular focus on the gluon distribution at large momentum fraction, x , and the corresponding partonic luminosities. We assess constraints from these datasets using next-to-next-to-leading-order (NNLO) theoretical predictions, accounting for a range of uncertainties from scale dependence and numerical integration. From the scale choices available for the calculations, our analysis shows that the central predictions for inclusive jet production show a smaller scale dependence than dijet production. We examine the relative constraints on the gluon distribution provided by the inclusive jet and dijet distributions, and also explore the phenomenological implications for inclusive H , t t ¯ , and t t ¯ H production at the LHC at 14 TeV. Published by the American Physical Society 2025

Ablat, Alim (ORCID:0009000519877755)↗
DOE OSTI2020

Particle-number distribution in large fluctuations at the tip of branching random walks

Here, we investigate properties of the particle distribution near the tip of one-dimensional branching random walks at large times t , focusing on unusual realizations in which the rightmost lead particle is very far ahead of its expected position, but still within a distance smaller than the diffusion radius ~$\sqrt{t}$. Our approach consists in a study of the generating function $G_{Δx}(λ) = Σ_n$ ${λ^n}p_n(Δx)$ for the probabilities $p_n(Δx)$ of observing $\textit{n}$ particles in an interval of given size $Δ\textit{x}$ from the lead particle to its left, fixing the position of the latter. This generating function can be expressed with the help of functions solving the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) equation with suitable initial conditions. In the infinite-time and large-$Δ\textit{x}$ limits, we find that the mean number of particles in the interval grows exponentially with $Δ\textit{x}$, and that the generating function obeys a nontrivial scaling law, depending on $Δ\textit{x}$ and λ through the combined variable $[Δx — f(λ)]^3 / Δx^2$, where $\textit{f}$(λ) ≡ – ln(1 – λ) – ln [– ln(1 – λ)]. From this property, one may conjecture that the growth of the typical particle number with the size of the interval is slower than exponential, but, surprisingly enough, only by a subleading factor at large Δ$\textit{x}$. The scaling we argue is consistent with results from a numerical integration of the FKPP equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Discrete spherical harmonic functions for texture representation and analysis

A basis of discrete harmonic functions for efficient representation and analysis of crystallographic texture is presented. Discrete harmonics are a numerical representation of the harmonics on the sphere. A finite element formulation is utilized to calculate these orthonormal basis functions, which provides several advantageous features for quantitative texture analysis. These include high-precision numerical integration, a simple implementation of the non-negativity constraint and computational efficiency. Simple examples of pole figure and texture interpolation and of Fourier filtering using these basis sets are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Quantum Approach for Implementing Fixed-Point Arithmetic in Solving Ordinary Differential Equations

Differential equations (DEs) serve as fundamental tools in mathematical modeling across scientific disciplines, yet classical numerical solvers face limitations with large-scale or computationally intensive problems. This study explores a quantum-inspired approach to solving DEs, combining quantum-inspired techniques with classical methods. It focuses on fixed-point arithmetic on quantum circuits, utilizing basic quantum gates to manipulate DE solutions. We expand upon the techniques introduced by Zanger et al. [Quantum, 5, 502 (2021)] by offering a precise computation for a fixed-point signed multiplication scheme, while also presenting a quantum circuit capable of executing the fixed-point division algorithm. We demonstrate the feasibility of our approach through the simulation of a linear Ordinary Differential Equation (ODE), where initial conditions and parameters are encoded into quantum circuits using fixed-point representation. By executing sequences of quantum gates mimicking numerical integration steps, we obtain approximate solutions to the ODE with specified fixed-point precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Provably-Stable Overload Ride-Through Control for Grid-Forming Inverters Using System-Wide Lyapunov Function Analysis

A key challenge associated with a grid-forming (GFM) inverter based resource (IBR) is its behavior during severe grid disturbances: since a GFM inverter regulates voltage in the fast timescale instead of current or power, it may experience a transient overload of current, power and/or energy during a severe grid disturbance. While many promising control strategies for overload ride-through have been proposed over the past two decades, transient stability of the system during and after the transition to an overload ride-through control mode remains difficult to guarantee. In this work, a novel overload ride-through control strategy is proposed for a system of grid-forming inverters that takes both self-protection and system-wide transient stability into account. A proposed system-level supervisory control uses slow communication to pre-emptively assign a set of local ride through control parameters to individual GFM IBR, including a current-limiting virtual reactance, that guarantees that synchronism is still preserved for any set of anticipated grid disturbances. At the core of the supervisory control lies a Lyapunov-function-based routine capable of establishing a strong, albeit conservative, transient stability guarantee for the system. Here, the proposed overload ride-through control strategy is validated via numerical integration of a reduced-order model, as well as through detailed electromagnetic transient (EMT) simulation.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Fiber Uncertainty Visualization for Bivariate Data With Parametric and Nonparametric Noise Models

Visualization and analysis of multivariate data and their uncertainty are top research challenges in data visualization. Constructing fiber surfaces is a popular technique for multivariate data visualization that generalizes the idea of level-set visualization for univariate data to multivariate data. Here, in this paper, we present a statistical framework to quantify positional probabilities of fibers extracted from uncertain bivariate fields. Specifically, we extend the state-of-the-art Gaussian models of uncertainty for bivariate data to other parametric distributions (e.g., uniform and Epanechnikov) and more general nonparametric probability distributions (e.g., histograms and kernel density estimation) and derive corresponding spatial probabilities of fibers. In our proposed framework, we leverage Green's theorem for closed-form computation of fiber probabilities when bivariate data are assumed to have independent parametric and nonparametric noise. Additionally, we present a nonparametric approach combined with numerical integration to study the positional probability of fibers when bivariate data are assumed to have correlated noise. For uncertainty analysis, we visualize the derived probability volumes for fibers via volume rendering and extracting level sets based on probability thresholds. We present the utility of our proposed techniques via experiments on synthetic and simulation datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Development and Validation of Passive Yaw in the Open-Source WEC-Sim Code

A passive yaw implementation is developed, validated, and explored for the WEC-Sim, an open-source wave energy converter modeling tool that works within MATLAB/Simulink. The Reference Model 5 (RM5) is selected for this investigation, and a WEC-Sim model of the device is modified to allow yaw motion. A boundary element method (BEM) code was used to calculate the excitation force coefficients for a range of wave headings. An algorithm was implemented in WEC-Sim to determine the equivalent wave heading from a body’s instantaneous yaw angle and interpolate the appropriate excitation coefficients to ensure the correct time-domain excitation force. This approach is able to determine excitation force for a body undergoing large yaw displacement. For the mathematically simple case of regular wave excitation, the dynamic equation was integrated numerically and found to closely approximate the results from this implementation in WEC-Sim. A case study is presented for the same device in irregular waves. In this case, computation time is increased by 32x when this interpolation is performed at every time step. To reduce this expense, a threshold yaw displacement can be set to reduce the number of interpolations performed. A threshold of 0.01° was found to increase computation time by only 22x without significantly affecting time domain results. Similar amplitude spectra for yaw force and displacements are observed for all threshold values less than 1°, for which computation time is only increased by 2.2x.

50 EE - Wind and Water Power Program - Water (EE-4↗
DOE OSTI2019

The Bloch equation for spin dynamics in electron storage rings: Computational and theoretical aspects

In this paper, we describe our work on spin polarization in high-energy electron storage rings which we base on the Full Bloch equation (FBE) for the polarization density and which aims towards the [Formula: see text] option of the proposed Future Circular Collider (FCC-ee) and the proposed Circular Electron Positron Collider (CEPC). The FBE takes into account non spin-flip and spin-flip effects due to synchrotron radiation including the spin-diffusion effects and the Sokolov–Ternov effect with its Baier–Katkov generalization as well as the kinetic-polarization effect. This mathematical model is an alternative to the standard mathematical model based on the Derbenev–Kondratenko formulas. For our numerical and analytical studies of the FBE, we develop an approximation to the latter to obtain an effective FBE. This is accomplished by finding a third mathematical model based on a system of stochastic differential equations (SDEs) underlying the FBE and by approximating that system via the method of averaging from perturbative ODE theory. We also give an overview of our algorithm for numerically integrating the effective FBE. This discretizes the phase space using spectral methods and discretizes time via the additive Runge–Kutta (ARK) method which is a high-order semi-implicit method. We also discuss the relevance of the third mathematical model for spin tracking.

Physics↗
DOE OSTI2020

Analytic stability boundaries for compressional and global Alfven eigenmodes driven by fast ions. II. Interaction via Landau resonance.

Conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfven eigenmodes driven by the Landau resonance with super-Alfvenic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite \omega/\omega_{ci} and k_\parallel/k_\perp is retained and found to be responsible for the existence of the GAE instability via this resonance. Encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Bug fixes / enhancements to Parallel tempering Markov chain Monte Carlo sampler (PTMCMCSampler) v3

Markov chain Monte Carlo sampler useful for a variety of problems for optimization, numerical integration, and generating draws from a probability distribution. They are often used in physical and mathematical problems, especially when it is difficult or impossible to use other approaches. This disclosure is for new bug fixes and enhancements to an existing open source library that is already in the public domain under the MIT license. https://github.com/jellis18/PTMCMCSampler

Forrer, Mark↗
DOE OSTI2026

Kinetic Model of HoxEFU reduction by NADH [SWR-26-087]

This repository is used to release code generated for manuscripts on the Photosynthetic Energy Transduction core program. This code simulates the reduction of HoxEFU by NADH. The electro transfer rate constants for the simulation are specified in the .csv files. The two .csv files correspond tot he two models described in Dawson et al. Cell. Rep. Phys. Sci. 2026. The code utilizes a chemical master equation, a set of differential equations, defining the time evolution of the oxidation and reduction kinetics of NAD+, NADH, a FMN flavin, and a set of iron sulfur clusters. The kinetics of HoxEFU reduction by NADH are evaluated by numerical integration of the chemical master equation using a variable-time-step Runge-Kutta algorithm.

Dahl, Peter [National Laboratory of the Rockies (N↗
DOE OSTI2022

Theory of In-Cloud Activation of Aerosols and Microphysical Quasi-Equilibrium in a Deep Updraft

The microphysical quasi-equilibrium in an ascending adiabatic parcel is elucidated by an analytical theory in 0D with drastic simplifications. The theory predicts how in-cloud activation is most likely to be triggered by the onset of precipitation during sufficient ascent, with the ascent only needing to approach almost twice the cloud-base updraft speed aloft. The precipitation and cloud condensate mass fields are coupled in a closed system in a simplified version of the model. The initial state of no precipitation is unstable with respect to a perturbation. In the 2D phase space of both mass fields, there is a neutral line. Unstable growth of precipitation mass drives the system to cross the line into a regime of stability. An attractor is then approached consisting of precipitation mass balanced by accretion of cloud mass and fallout. The cloud-particle number concentration also approaches a stable equilibrium governed by the balance between in-cloud activation aloft from the inexorably increasing supersaturation during vertical acceleration and accretion of cloud droplets by precipitation. Dimensionless numbers characterizing the microphysical equilibria and their stability are derived mathematically, including a condensation–precipitation efficiency and an in-cloud activation efficiency. The theory explains common observations of the orders of magnitude of liquid water content in convective and stratiform clouds. Finally, sensitivity tests of the numerically integrated theoretical equations are documented with respect to variations in cloud condensation nucleus (CCN) aerosols and updraft speed. In conclusion, it is shown that this theory of in-cloud activation, with an increasing supersaturation during ascent, applies to both ice-only and liquid-only cloud.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Analytic phase solutions of three-wave interactions

Closed-form analytic phase solutions of three-wave interactions are presented for the first time. The cases from simple second harmonic generation to most general three wave interactions without any constraints are considered. The phase amplification or deamplification behavior in the phase-sensitive parametric process is illustrated using the results obtained here. The analytic solutions agree with the results from direct numerical integration. The validity range of an approximate phase solution is discussed.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Utah FORGE: Optimization of a Plug-and-Perf Stimulation (Fervo Energy)

Information around the plug-and-perf treatment design at Utah FORGE by Fervo Energy. Objective and Purpose: - Develop a multistage hydraulic stimulation approach designed specifically to target the top three factors that control the technical and commercial viability of an EGS system: i) Achieving sufficient injectivity to support high cross-well flow rates ii) Distributing flow evenly across the wellbore and reservoir to maximize heat mining efficiency, ensure sustained heat transfer, and mitigate thermal breakthrough iii) Overcoming the effects of stress heterogeneity, stress shadowing, and variations in natural fracture properties during the stimulation treatment, leading to a more predictable stimulated reservoir volume and offset well placement - The following activities will be performed: i) Design, plan, and execute a multistage plug-and-perf stimulation treatment at a Fervo site with data acquisition and well testing activities aimed at addressing key technical aspects of the issues above ii) Perform data processing and interpretation of field results to translate the results form the Fervo site to a site-specific design at the Utah FORGE site iii) Design, plan, and execute a multistage plug-and-perf stimulation treatment design at the Utah FORGE site Methods and Approach: - Design a detailed data acquisition plan to maximize learning around: i) DFIT testing ii) Petrophysical logging, image logging iii) Permanent DAS/DTS fiber optic monitoring iv) Deep borehole microseismic monitoring v) Shallow borehole induced seismicity monitoring vi) Injection/production testing (RTA analysis, tracer testing) vii) Integrated numerical modeling and production forecasting

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Machine learning for postprocessing ensemble streamflow forecasts

Skillful streamflow forecasts can inform decisions in various areas of water policy and management. We integrate numerical weather prediction ensembles, distributed hydrological model, and machine learning to generate ensemble streamflow forecasts at medium-range lead times (1–7 days). We demonstrate the application of machine learning as postprocessor for improving the quality of ensemble streamflow forecasts. Our results show that the machine learning postprocessor can improve streamflow forecasts relative to low-complexity forecasts (e.g., climatological and temporal persistence) as well as standalone hydrometeorological modeling and neural network. The relative gain in forecast skill from postprocessor is generally higher at medium-range timescales compared to shorter lead times; high flows compared to low–moderate flows, and the warm season compared to the cool ones. Overall, our results highlight the benefits of machine learning in many aspects for improving both the skill and reliability of streamflow forecasts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

SCEPTRE 2.2 Angular Quadrature Sets

This document includes details of the angular quadrature sets available in SCEPTRE for performing numerical integrations in the angular phase space. The angular dependence of the boundary and fixed-source terms an d initial angular flux are specified by angular index rather than by direction. It is, therefore, necessary to know the mapping from a specific direction to a direction index. This document includes angular quadrature weights and direction cosines for most of the quadrature sets available in SCEPTRE.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗