Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gradient flow”

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 OSTI2026

High-precision scale setting with the Ω-baryon mass and gradient flow

The gradient-flow scale 𝑤 0 in lattice QCD is determined using the mass of the Ω − baryon to set the physical scale. Nine ensembles using the highly improved staggered quark (HISQ) action with lattice spacings of 0.15 fm down to 0.04 fm are used, seven of which have nearly physical light-quark masses. Electromagnetic corrections to the Ω − mass are defined in order to compute a pure-QCD Ω mass. The final result is 𝑤 0 = 0.17187⁢(68) fm, corresponding to a relative uncertainty of 0.40% and a central value in good agreement with previous calculations in the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Gradient flow based phase-field modeling using separable neural networks

Allen–Cahn equation is a reaction–diffusion equation and is widely used for modeling phase separation. Machine learning methods for solving the Allen–Cahn equation in its strong form suffer from inaccuracies in collocation techniques, errors in computing higher-order spatial derivatives, and the large system size required by the space–time approach. To overcome these challenges, we propose solving the gradient flow of the Ginzburg–Landau free energy functional, which is equivalent to the Allen–Cahn equation, thereby avoiding the second-order spatial derivatives associated with the Allen–Cahn equation. A minimizing movement scheme is employed to solve the gradient flow problem, eliminating the complexities of a space–time approach. We utilize a separable neural network that efficiently represents the phase field through low-rank tensor decomposition. As we use the minimizing movement scheme to numerically solve the gradient flow problem, we thus, refer to the proposed method as the Separable Deep Minimizing Movement (SDMM) method. The evaluation of the functional in the minimizing movement scheme using the Gauss quadrature technique bypasses the inaccuracies associated with collocation techniques traditionally used to solve partial differential equations. A hyperbolic tangent transformation is introduced on the phase field prior to the evaluation of the functional to ensure that it remains strictly bounded within the values of the two phases. For this transformation, theoretical guarantee for energy stability of the minimizing movement scheme is established. Our results suggest that this transformation helps to improve the accuracy and efficiency significantly. The proposed method resolves the challenges faced by state-of-the-art machine learning techniques, outperforming them in both accuracy and efficiency. It is also the first machine learning method to achieve an order of magnitude speed improvement over the finite element method. In addition to its formulation and computational implementation, several case studies illustrate the applicability of the proposed method.

42 ENGINEERING↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
DOE OSTI2026

Accessing the gluon momentum fraction of nucleons through the gradient flow

We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 ⁢(𝜇 =2 GeV) =0.482⁢(35), where we quote only statistical uncertainties.

Lattice QCD↗
DOE OSTI2026

Renormalized quark masses using gradient flow

We propose a new and simple method for determining the renormalized quark masses from lattice simulations. Renormalized quark masses are an important input to many phenomenological applications, including searching and modeling physics beyond the Standard Model. The nonperturbative renormalization is performed using gradient flow combined with the short-flow-time expansion that is improved by renormalization group (RG) running to match to the $\overline{MS}$ scheme. Implementing the RG running perturbatively, we demonstrate this method works reliably at least up to the charm-quark mass and exhibits an easily attainable “windowing condition.” Using RBC/UKQCD’s (2+1)-flavor Shamir domain-wall fermion ensembles with Iwasaki gauge action, we find 𝑚 $\overline{MS}$ 𝑠⁡ (𝜇 = 2 GeV) = 89⁢(3) MeV and 𝑚 $\overline{MS}$ 𝑐 ⁡(𝜇 = 3 GeV) = 972⁢(16) MeV. These results predict the scale-independent ratio 𝑚 𝑐 /𝑚 𝑠 = 12.1⁢(4). Generalization to other observables is possible, providing an efficient approach to determine nonperturbatively renormalized fermionic observables like form factors or bag parameters from lattice simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Gradient Flow for Parton Distribution Functions: First Application to the Pion

Parton distribution functions (PDFs) are central to precision QCD phenomenology. Their Mellin moments can be computed on the lattice, but direct determinations using local operators, besides ⟨𝑥⟩, face severe challenges from reduced hypercubic symmetry, limiting results to the lowest moments. A recently proposed method resolves these issues using gradient flow. We demonstrate the efficacy of this method by computing ratios of flavor nonsinglet pion PDF moments up to ⟨𝑥 5 ⟩ on four lattice spacings at 𝑚 𝜋 ≃ 411 MeV. The moments and reconstructed PDF agree quantitatively with recent phenomenological extractions.

Francis, Anthony [National Yang Ming Chiao Tung Un↗
DOE OSTI2026

Scale setting of SU⁡(𝑁) Yang–Mills theory, topology and large-𝑁 volume independence

We set the scale of SU⁡(𝑁) Yang-Mills theories for 𝑁 =3, 5, 8 and in the large-𝑁 limit via gradient flow, as a first step towards the computation of the large-𝑁 Λ-parameter using step scaling. We adopt twisted boundary conditions to achieve large-𝑁 volume reduction and the Parallel Tempering on Boundary Conditions algorithm to tame topological freezing. This setup allows accurate determinations of the gradient-flow scales down to lattice spacings as fine as ∼0.025 fm for all the explored values of 𝑁, a regime that has never been reached with ergodic algorithms. Moreover, we are able to precisely estimate the finite-size systematics related to topological freezing, and to show the suppression of finite-volume effects expected by virtue of large-𝑁 twisted volume reduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2025

Effects of an External Magnetic Field on Keyhole Mode Laser Melting of 316 Stainless Steel

Keyhole-mode laser melting is an efficient method for joining or cutting large, thick components, but controlling keyhole depth and fluctuations has remained challenging. Applying an external magnetic field can control melt pool flows and indirectly influence keyhole morphology and dynamics. The induced Lorentz force, comprising Seebeck and damping components, plays a crucial role in the melt pool dynamics, depending on temperature gradient, flow rates, and magnetic field orientation and magnitude. This research investigates the effects of an external magnetic field on keyhole behavior during laser spot melting of 316 stainless steel using synchronized high-speed synchrotron X-ray and thermal imaging. Findings revealed that a longitudinal magnetic field (120 mT) increased keyhole depth but exacerbated lateral fluctuations, resulted in a 20% increase in the melt pool temperature gradient and a 27% decrease in cooling rate. Conversely, a transverse magnetic field (760 mT) reduced keyhole depth and improved porosity formation. The findings suggest that a decrease in keyhole depth correlates with a decrease in fluctuations, and vice versa. These insights enhance understanding of external magnetic fields’ impact on laser melting, with implications for improving part quality.

Alamdari, Aslan Bafahm [The Ohio State Univ., Colu↗
DOE OSTI2024

Vortex rings in event-by-event relativistic heavy-ion collisions

We present event-by-event simulations for central asymmetric light+heavy and Au +Au collisions to investigate the formation and evolution of vortex-ring structures in the longitudinal flow velocity profile. The production-plane polarization of Λ hyperons, defined with respect to the Λ momentum and the beam, can track the “vortex-ring” feature in the event, a characteristic vortical structure generated by longitudinal flow gradients. We make comprehensive model predictions for the rapidity-dependent vortex-ring observables for different collision system sizes at √s NN = 200 and 72 GeV. Furthermore, our predictions at the latter energy can be explored in the future LHCb fixed-target experiment at the Large Hadron Collider.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Lattice study of correlators of chromoelectric fields for heavy quarkonium dynamics in the quark-gluon plasma

We perform a lattice calculation of the correlators of two chromoelectric fields in the adjoint representation connected by adjoint Wilson lines at nonzero temperature. These correlators arise in the study of quarkonium dynamics and of adjoint heavy quark diffusion in deconfined matter. We work in SU(3) gauge theory using either gradient flow or multilevel algorithms for noise reduction, and discuss the renormalization of the correlators on the lattice. We find that a Casimir factor rescaling relates the adjoint correlators corresponding to the diffusion of an adjoint heavy quark and the octet-octet quarkonium transitions to the chromoelectric correlator in the fundamental representation describing the diffusion of a heavy quark.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Enhancement to Fusion Reactivity in Sheared Flows

Sheared flow increases the reactivity of fusion plasma. In unmagnetized plasma with flow gradients comparable to the mean free path of reacting ions, fusion reactivity can be more than doubled. The effect is of particular relevance to inertial confinement fusion (ICF), where it allows implosion kinetic energy to contribute to the fusion burn even before thermalizing. Finally, because colder fuel stops alpha particles more quickly, ignition is possible in a smaller volume, substantially reducing energy requirements in fast-ignition designs.

fast ignition↗
DOE OSTI2026

Improved energies and local energies with weighted variational Monte Carlo

Neural network parametrizations have increasingly been used to represent the ground and excited states in variational Monte Carlo (VMC) with promising results. However, traditional VMC methods only optimize the wave function in regions of peak probability. The wave function is uncontrolled in the tails of the probability distribution, which can limit the accuracy of the trained wave function. To improve the approximation accuracy in the probability tails, this paper interprets VMC as a gradient flow in the space of wave functions, followed by a projection step. From this perspective, arbitrary probability distributions can be used in the projection step, allowing the user to prioritize accuracy in different regions of state space. Motivated by this theoretical perspective, the paper tests a weighted VMC method on the antiferromagnetic Heisenberg model for a periodic spin chain. Compared to traditional VMC, weighted VMC reduces the error in the ground state energy by a factor of 2, and it reduces the errors in the local energies away from the mode by large factors of 10 2 –10 4 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Scale Setting with the Omega Baryon

This work will discuss a recent calculation of the gradient flow scale $w_0$ in lattice QCD using the omega baryon mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Towards determination of the strong coupling $α_s(m_Z)$ from four-flavor lattice QCD using the continuous $β$-function method

The precise value of the strong coupling $α_s(m_{Z})$ at the $Z$-boson mass $m_{Z}$ is essential for high-energy phenomenology and precision tests of quantum chromodynamics (QCD). We present the status of a program targeting a $\sim 0.3\%$ determination of $α_s(m_{Z})$ using the renormalization group $β$-function in the infinite volume gradient flow scheme based on lattice QCD simulations of degenerate four-flavor highly improved staggered quark (HISQ) ensembles. In particular, we analyze both tree-level cutoff effects and finite-mass effects. We also outline the next steps of the analysis, including the infinite-volume and continuum extrapolations required for a precise determination of $α_s(m_Z)$.

Mandlecha, Yash (ORCID:000000020587962X)↗
DOE OSTI2026

Towards determination of the strong coupling $α_s(m_Z)$ from four-flavor lattice QCD using the continuous $β$-function method

The precise value of the strong coupling $α_s(m_{Z})$ at the $Z$-boson mass $m_{Z}$ is essential for high-energy phenomenology and precision tests of quantum chromodynamics (QCD). We present the status of a program targeting a $\sim 0.3\%$ determination of $α_s(m_{Z})$ using the renormalization group $β$-function in the infinite volume gradient flow scheme based on lattice QCD simulations of degenerate four-flavor highly improved staggered quark (HISQ) ensembles. In particular, we analyze both tree-level cutoff effects and finite-mass effects. We also outline the next steps of the analysis, including the infinite-volume and continuum extrapolations required for a precise determination of $α_s(m_Z)$.

Mandlecha, Yash [Michigan State U.; Michigan State↗
DOE OSTI2025

Strong gradient effects on neoclassical electron transport and the bootstrap current in large aspect ratio tokamaks

Standard approaches to neoclassical theory do not extend into regions of strong gradients in tokamaks such as the pedestal and internal transport barriers. Here, we calculate the modifications to neoclassical electron physics inside strong gradient regions of large aspect ratio tokamaks in the banana regime. We show that these modifications are due to the different ion flow and the strong poloidal variation of the potential. We also provide a physical interpretation of the mechanisms that drive poloidal asymmetries and hence a poloidal electric field. We apply our model to two specific example cases of pedestal profiles, calculating the neoclassical electron flux and the bootstrap current. We find that, depending on the ion flow, weak gradient neoclassical theory overestimates or underestimates the neoclassical electron transport and the bootstrap current in regions with strong gradients. We show that the determination of the mean parallel flow is more complex than in weak gradient neoclassical theory. For vanishing turbulence, we can determine the radial electric field for a given flow profile in the pedestal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗