Search NASA⌕ Search

SEARCH · Search NASA

Results for “function approximation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Analytic solutions of the DGLAP evolution and theoretical uncertainties

The energy dependence for the singlet sector of Parton Distributions Functions (PDFs) is described by an entangled pair of ordinary linear differential equations. Although there are no exact analytic solutions, it is possible to provide approximated results depending on the assumptions and the methodology adopted. These results differ in their sub-leading, neglected terms and ultimately they are associated with different treatments of the theoretical uncertainties. In this work, a novel analytic approach in Mellin space is presented and a new methodology for obtaining closed and exponentiated analytic solutions is devised. Different results for the DGLAP evolution at Next-Leading-Order are compared, discussing advantages and disadvantages for each solution. The generalizations to higher orders are addressed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Tethered Balloon Observations of Vertically Resolved Aerosol Size Distributions during the U.S. DOE ARM CoURAGE and BNF Campaigns

This dataset contains vertically resolved aerosol particle number size distributions measured using a miniature Scanning Electrical Mobility Spectrometer (mSEMS) deployed aboard the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) Tethered Balloon System (TBS) during the 2025 Coast-Urban-Rural Atmospheric Gradient Experiment (CoURAGE) and Bankhead National Forest (BNF) field campaigns. The mSEMS measured aerosol size distributions over approximately 10–300 nm during repeated TBS vertical profiles through the atmospheric boundary layer. The dataset provides measurements of particle number size distributions as a function of particle diameter, altitude, and time and supports investigation of aerosol vertical structure, elevated ultrafine and Aitken-mode particle layers, boundary-layer evolution, and residual-layer processes.

ARM↗
DOE OSTI2024

Measurement report: A comparison of ground-level ice-nucleating-particle abundance and aerosol properties during autumn at contrasting marine and terrestrial locations

Abstract. Ice-nucleating particles (INPs) are an essential class of aerosols found worldwide that have far-reaching but poorly quantified climate feedback mechanisms through interaction with clouds and impacts on precipitation. These particles can have highly variable physicochemical properties in the atmosphere, and it is crucial to continuously monitor their long-term concentration relative to total ambient aerosol populations at a wide variety of sites to comprehensively understand aerosol–cloud interactions in the atmosphere. Hence, our study applied an in situ forced expansion cooling device to measure ambient INP concentrations and test its automated continuous measurements at atmospheric observatories, where complementary aerosol instruments are heavily equipped. Using collocated aerosol size, number, and composition measurements from these sites, we analyzed the correlation between sources and abundance of INPs in different environments. Toward this aim, we have measured ground-level INP concentrations at two contrasting sites, one in the Southern Great Plains (SGP) region of the United States with a substantial terrestrially influenced aerosol population and one in the Eastern North Atlantic Ocean (ENA) region with a primarily marine-influenced aerosol population. These measurements examined INPs mainly formed through immersion freezing and were performed at a ≤ 12 min resolution and with a wide range of heterogeneous freezing temperatures (Ts above −31 °C) for at least 45 d at each site. The associated INP data analysis was conducted in a consistent manner. We also explored the additional offline characterization of ambient aerosol particle samples from both locations in comparison to in situ data. From our ENA data, on average, INP abundance ranges from ≈ 1 to ≈ 20 L−1 (−30 °C ≤ T ≤ −20 °C) during October–November 2020. Backward air mass trajectories reveal a strong marine influence at ENA with 75.7 % of air masses originating over the Atlantic Ocean and 96.6 % of air masses traveling over open water, but analysis of particle chemistry suggests an additional INP source besides maritime aerosols (e.g., sea spray aerosols) at ENA. In contrast, 90.8 % of air masses at the SGP location originated from the North American continent, and 96.1 % of the time, these air masses traveled over land. As a result, organic-rich SGP aerosols from terrestrial sources exhibited notably high INP abundance from ≈ 1 to ≈ 100 L−1 (−30 °C ≤ T ≤ −15 °C) during October–November 2019. The probability density function of aerosol surface area-scaled immersion freezing efficiency (ice nucleation active surface site density; ns) was assessed for selected freezing temperatures. While the INP concentrations measured at SGP are higher than those of ENA, the ns(T) values of SGP (≈ 105 to ≈ 107 m−2 for −30 °C ≤ T ≤ −15 °C) are reciprocally lower than ENA for approximately 2 orders of magnitude (≈ 107 to ≈ 109 m−2 for −30 °C ≤ T ≤ −15 °C). The observed difference in ns(T) mainly stems from varied available aerosol surface areas, Saer, from two sites (Saer,SGP > Saer,ENA). INP parameterizations were developed as a function of examined freezing temperatures from SGP and ENA for our study periods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Characterizing Artificial Viscosity Parameters with Approximate Symmetries

We are often faced with trying to capture the physics of compressible shocks, which are governed by the Euler equations. However, Euler shocks are formally discontinuous at the shock front (translating to a step-function behavior of rel evant flow variables). This poses a practical problem for codes with finite-sized grid elements. As a result, one must make a concession in simulating the be havior of shocks within a discrete framework. In particular, we must blur, or ‘regularize’ Euler shocks so that they may be captured on a finite grid.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Helicity evolution at small x : quark to gluon and gluon to quark transition operators

We include the quark to gluon and gluon to quark shock-wave transition operators into the small Bjorken-x evolution equations for helicity in the flavor-singlet channel derived earlier in [1,2,3]. While such transitions do not affect the large-N c version of the evolution equations for helicity, the large-N c & N f equations are affected. (N c and N f are the numbers of quark colors and flavors, respectively.) We derive the corresponding corrected large-N c & N f equations for the polarized dipole amplitudes contributing to the flavor-singlet quark and gluon helicity distributions in the double-logarithmic approximation (DLA), resumming powers of α s ln 2 (1/x) with αs the strong coupling constant. We solve these equations iteratively and extract the polarized splitting functions up to four loops. We show that our splitting functions agree with the fixed-order perturbative calculations up to and including the existing three-loops results [4,5,6,7]. Similar to the large-Nc helicity evolution in the shock-wave approach [8], our large-N c & N f small-x splitting functions agree with those obtained in the infrared evolution equations framework from [9, 10] up to three loops, but appear to slightly disagree at four loops.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Numerical Analysis of Offshore Wind-Farm-Induced Drag Effects on Coastal Upwelling Dynamics

Wind farms extract momentum from the atmospheric flow, generating wind-speed deficits both within the plant, and extending downstream. When located offshore, these deficits modulate air-sea coupling, potentially impacting coastal upwelling in sensitive regions. We investigate impacts of wind farms on coastal upwelling using kilometre-scale, three-way-coupled simulations with the coupled ocean-atmosphere-wave-sediment transport system for the US West Coast. Wind-farm effects are represented by a generalised turbine drag formulation, an idealised, height-dependent body force whose magnitude is systematically varied. This approach isolates the leading-order fluid-dynamical response in a realistic coastal configuration. The atmospheric adjustment exhibits an approximately linear relation between drag force and wind-speed deficit, with wakes that expand downstream and increase in magnitude as drag increases. An empirical orthogonal function analysis of sea-surface-temperature anomalies reveals the emergence of a canonical dipole pattern under strong drag forcing. Subsurface diagnostics show consistent shoaling of the mixed layer and suppressed upward velocities in areas near wind-farm region, accompanied by compensating enhancements of shoaling closer to the coast. These results identify turbine drag as a control parameter in assessing interactions between wind-farm wake and coastal upwelling and provide scaling relationships for understanding offshore wind-farm effects on the coastal circulation dynamics.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Understanding Majorana braiding in superconductors with a fixed total number of particles

Mean-field one-dimensional topological superconductors host edge Majorana zero modes (MZMs) that encode a topologically protected ground-state degeneracy and enable robust braiding operations within this subspace. Mean-field states lack definite particle number and thus cannot represent isolated systems. Nevertheless, projecting them onto fixed particle number can yield good approximations to an isolated superconductor ground state. In earlier work [Sajith et al., Phys. Rev. B 109, 184509 (2024)] we showed that the projected Kitaev wave function of a single wire preserves some important mean-field features, such as the zero-energy spectral peaks near the wire edges. However, uniqueness of the fixed-number ground state does not allow for any nontrivial operations in the ground-state subspace. To overcome this limitation, here we consider the case of multiple wires with a conserved total charge, using the same approach. In the limit of vanishing interwire coupling, this system has a macroscopic ground-state degeneracy. We show how this degeneracy is resolved by coherent single-particle tunneling, and identify special many-body states which play the same role as the MZM parity states in the mean field. In this work, we demonstrate how braiding operations can be implemented and discuss both intrinsic and extrinsic limits on their fidelity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

QCD splitting functions beyond kinematical limits

We present a systematic decomposition of QCD splitting functions into scalar dipole radiators and pure splitting remainders up to second order in the strong coupling. The individual components contain terms that are formally subleading in soft or collinear scaling parameters, but well understood and universal due to their origin in scalar QCD. The multipole radiator functions which we derive share essential features of the known double-soft and one-loop soft gluon currents, and are not based on kinematical approximations.

Campbell, John M. [Fermilab] (ORCID:00000002136193↗
DOE OSTI2025

Quantum Hamilton-Jacobi theory, spectral path integrals, and exact WKB analysis

We propose a new way to perform path integrals in quantum mechanics by using a quantum version of Hamilton-Jacobi (HJ) theory. In classical mechanics, Hamilton-Jacobi theory is a powerful formalism, however, its utility is not explored in quantum theory beyond approximation schemes. The canonical transformation enables one to set the new Hamiltonian to constant or zero, but keeps the information about solution in Hamilton’s characteristic function. To benefit from this in quantum theory, one must work with a formulation in which classical Hamiltonian is used. This uniquely points to phase space path integral. However, the main variable in HJ formalism is energy, not time. Thus, we are led to consider the Fourier transform of the path integral, the spectral path integral Z ˜ ( E ) . The evaluation of path integrals reduces to determining the quantum Hamilton characteristic functions (which can be achieved via an asymptotic analysis) and a discrete sum over the quantum period lattice, generalizing Gutzwiller’s sum. Published by the American Physical Society 2025

Türe, Mustafa (ORCID:0009000975968618)↗
DOE OSTI2025

Tracing the Cosmic Evolution of the Cool Circumgalactic Medium of Luminous Red Galaxies with DESI Year 1 Data

We investigate the properties of the cool circumgalactic medium (CGM) of massive galaxies and their cosmic evolution. By using the year 1 dataset of luminous red galaxies (LRGs) and QSOs from the Dark Energy Spectroscopic Instrument survey, we construct a sample of approximately 600,000 galaxy-quasar pairs and measure the radial distribution and kinematics of the cool gas traced by Mg II absorption lines as a function of galaxy properties from redshift 0.4 to redshift 1.2. Our results show that the covering fraction of the cool gas around LRGs increases with redshift, following a trend similar to the global evolution of galaxy star formation rate. At small radii (< 0.3rvir), the covering fraction anti-correlates with stellar mass, suggesting that mass-dependent processes suppress the cool gas content in the inner region. In addition, we measure the gas dispersion by modeling the velocity distribution of absorbers with a narrow and a broad components -- sigma_n ~ 160 and sigma_b ~ 380 km/s -- and quantify their relative contributions. The results show that the broad component becomes more prominent in the outer region, and its relative importance in the central region grows with increasing stellar mass. Finally, we discuss possible origins of the cool gas around massive galaxies, including the contribution of satellite galaxies and the precipitation scenario.

Chang, Yu-Ling [Taiwan, Natl. Taiwan U.] (ORCID:00↗
DOE OSTI2025

DIF3D-VARIANT 12.0: Updates and New Features

The DIF3D code has been a workhorse of fast reactor analysis work at Argonne National Laboratory for over 40 years. In 1995, a transport option called VARIANT was added to DIF3D to improve the flux solutions for fast reactor problems which we term DIF3D-VARIANT today. DIF3D-VARIANT performs nodal neutron transport calculations using P N or SP N theory in Cartesian and hexagonal two- and three-dimensional geometries. The limited computing capabilities of the time restricted DIF3D-VARIANT to use at most a 6 th order spatial approximation combined with a P3 flux approximation and P1 scattering kernel for a 33 group structure on most studied reactor problems. Computer capabilities have increased steadily since 1995 and today much larger space-angle-energy approximations are possible. This manuscript serves as an update to the theory section of the original DIF3D-VARIANT manual and details more than twenty years of changes made to DIF3D to make version 12 which was released on November 1 st , 2024. The primary focus of the initial work was to extend the space-angle approximations available in DIF3D-VARIANT such that the error due to transport approximations could be better understood. This work was started and completed in 2002 and marked the official version 10. Unfortunately, those higher order approximations could not be used at that time due to the memory constraints of the BPOINTER part of DIF3D (limited to 2 GB). In version 11, completed in 2012, BPOINTER was circumvented in DIF3D-VARIANT for the largest arrays by introducing a Fortran 90 module called LMA (Large Memory Array). This seamlessly replaces all of the functionality of the BPOINTER concept, but it allows 64 bit addressing for every array such that they can be larger than 2 GB. It is now common for DIF3D-VARIANT jobs to consume 50 GB of memory on modern workstations when using high order space-angle approximations and a large number of groups. Many improvements were made to version 11 from 2012 to 2022 when work to create version 12 started. For version 12, several parts of DIF3D were updated to improve performance and thread parallelism was introduced to further reduce the runtime. Numerous minor bugs were discovered in DIF3D-VARIANT as part of the process of creating the perturbation and sensitivity code PERSENT. All of these algorithmic problems were identified in the transition from version 10 to version 11 which prevented DIF3D-VARIANT from running efficiently and reliably. Firstly, the coarse mesh rebalance scheme would routinely diverge and a study detailed in this report demonstrates how it was also typically not effective. This is not a failure of the coarse mesh rebalance methodology, but a failure of its implementation in DIF3D-VARIANT for hexagonal geometries. The fission source extrapolation algorithm was also found to be unreliable on larger group structure problems, leading to divergence in some cases and a negligible improvement in performance overall. Finally, the “Omega” acceleration applied to the partial current solver routine of DIF3D-VARIANT was found to cause DIF3D-VARIANT to converge to the wrong answer. To resolve these issues, both the coarse mesh rebalance and fission source extrapolation were permanently disabled in version 11. The Tchebychev acceleration was put in as a temporary reliable alternative but it is generally inferior to coarse mesh rebalance or coarse mesh finite difference. For the Omega acceleration, the factor was restricted to guarantee that it would not cause follow-on errors in PERSENT. Due to limited funding to support maintenance and development of DIF3D in the last 10 years, no effort was spent since to resolve the outer iteration acceleration. Except for the threading work, all of the changes discussed in this manuscript refer to changes made between version 10 and version 11. Performance comparisons are done to demonstrate the improvements from version 9 to version 12. As will be demonstrated, the updated versi

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Exclusive J / ψ , ψ ( 2 s ) , and e + e - pair production in Au + Au ultraperipheral collisions at the BNL Relativistic Heavy Ion Collider

Measurements of exclusive J/ψ, ψ(2⁢s), and electron-positron (e + ⁢e - ) pair photoproduction in Au+Au ultraperipheral collisions are reported by the STAR experiment at $\sqrt{s_{NN}}$ = 200GeV. Here, we report several first measurements at the BNL Relativistic Heavy Ion Collider, which are (i) J/ψ photoproduction with large momentum transfer up to 2.2⁢(GeV/c) 2 , (ii) coherent J/ψ photoproduction associated with neutron emissions from nuclear breakup, (iii) the rapidity dependence of incoherent J/ψ photoproduction, (iv) the ψ⁡(2⁢s) photoproduction cross section at midrapidity, and (v) e + ⁢e - pair photoproduction up to high invariant mass of 6GeV/c 2 . For measurement (ii), the coherent J/ψ total cross section of γ + Au → J/ψ + Au as a function of the center-of-mass energy W γN has been obtained without photon energy ambiguities. The data are quantitatively compared with the Monte Carlo models STARlight, Sartre, BeAGLE, and theoretical calculations of gluon saturation with color glass condensate, nuclear shadowing with leading twist approximation, quantum electrodynamics, and the next-to-leading-order perturbative QCD. At the photon-nucleon center-of-mass energy of 25.0 GeV, the coherent and incoherent J/ψ cross sections of Au nuclei are found to be 71% ± 10% and 36% ± 7%, respectively, of that of free protons. These data provide an important experimental constraint for nuclear parton distribution functions and a unique opportunity to advance the understanding of the nuclear modification effect at the top RHIC energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

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

Resonant Raman in armchair graphene nanoribbons from first-principles

Resonant Raman spectra of armchair graphene nanoribbons (AGNRs) are computed using Density Functional Theory (DFT) and third-order perturbation theory. Results are benchmarked against available experimental data and compared to previously used theoretical approaches based on the Placzek approximation. Comparable agreement with experiments is found for both previously and presently used methods. In addition, a numerical analysis is carried out to provide a justification for the resonant modeling method based on the use of the frequency-dependent dielectric tensor in the Placzek approximation. Finally, this work also provides additional predictions and references for wide AGNRs that might be investigated with Raman scattering experiments in the future.

42 ENGINEERING↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Comparison between explicit and implicit discretization strategies for a dissipative thermal environment

We investigate strategies for simulating open quantum systems coupled to dissipative baths by comparing explicit wave function-based discretization [via multi-layer multi-configuration time-dependent Hartree (ML-MCTDH)] and the implicit density matrix-based master equation method [via tree tensor network hierarchical equations of motion (TTN-HEOM)]. For dissipative baths characterized by exponentially decaying bath correlation functions, the implicit discretization approach of HEOM—rooted in bath correlation function decompositions—proves significantly more efficient than explicit discretization of the bath into discrete harmonic modes. Explicit methods, like ML-MCTDH, require extensive mode discretization to approximate continuum baths, leading to computational bottlenecks. Case studies for two-level systems and a Fenna–Matthews–Olson complex model highlight TTN-HEOM’s superiority in capturing dissipative dynamics with relaxations with a minimal number of auxiliary modes, while the explicit methods are as exact as the HEOM in pure dephasing regimes. This comparison is enabled by the TENSO package, which has both ML-MCTDH and TTN-HEOM implemented using the same computational structure and propagation strategy.

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

Producing High-fidelity Synthetic Population Ensembles at Scale

Used within social simulations, synthetic population ensembles enable uncertainty quantification (UQ) methods for obtaining more robust model inference and prediction. A synthetic population ensemble is a series of plausible virtual reconstructions of an area’s population at the granularity of people and residences, generated stochastically to preserve privacy of the source population survey’s respondents. In this paper, we demonstrate the production of large synthetic population ensembles for the US via Oak Ridge National Laboratory’s UrbanPop framework to support modeling of high spatial resolution energy affordability metrics from nationwide social surveys in collaboration with the fusionACS project. Our initial task involves creating ensembles for 17 US metropolitan areas, each consisting of 41 population instances (a base realization and 40 replicates). To accomplish this task at scale, we configured an integrated system comprised of a research cloud, virtual containerization, GPU-enhanced functionality, and a dual API/CLI to interact with UrbanPop’s maturing Likeness Python ecosystem. We observe a reduction in theoretical execution time while maintaining high-fidelity approximations of residential totals by metropolitan area and the demographic characteristics of neighborhoods. We discuss expansion of our approach to produce synthetic population ensembles for the entire US, particularly plans to establish automated workflows for job orchestration to increase computational efficiency, as well as provide outlook for broadening applications of the ensembles.

Gaboardi, James [ORNL] (ORCID:0000000247766826)↗