Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fokker-Planck”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Experiments on excitation of Alfvén eigenmodes by alpha-particles with bump-on-tail distribution in JET DTE2 plasmas

Dedicated experiments were performed in JET DTE2 plasmas for obtaining an α-particle bump-on-tail (BOT) distribution aiming at exciting Alfvén eigenmodes (AEs). Neutral beam injection-only heating with modulated power was used so that fusion-born α-particles were the only ions present in the MeV energy range in these DT plasmas. The beam power modulation on a time scale shorter than the α-particle slowing down time was chosen for modulating the α-particle source and thus sustaining a BOT in the α-particle distribution. High-frequency modes in the toroidicity-induced Alfven eigenmode (TAE) frequency range and multiple short-lived modes in a wider frequency range have been detected in these DT discharges with interferometry, soft x-ray cameras, and reflectometry. The modes observed were localised close to the magnetic axis and were not seen in the Mirnov coils. Analysis with the TRANSP and Fokker-Planck FIDIT codes confirms that α-particle distributions with BOT in energy were achieved during some time intervals in these discharges though no clear correlation was found between the times of the high-frequency mode excitation and the BOT time intervals. The combined magneto-hydrodynamic (MHD) and kinetic modelling studies show that the high-frequency mode in the TAE frequency range is best fitted with a TAE of toroidal mode number n= 9. This mode is driven mostly by the on-axis beam ions while the smaller drive due to the pressure gradient of α-particles allows overcoming the marginal stability and exciting the mode. The observed multiple short-lived modes in a wider frequency range are identified as the on-axis kinetic AEs predicted in Rosenbluth and Rutherford.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Relativistic diffusion model for hadron production in p -Pb collisions at energies available at the CERN Large Hadron Collider

We investigate charged-hadron production in relativistic heavy-ion collisions of asymmetric systems within a nonequilibrium-statistical framework. Calculated centrality-dependent pseudorapidity distributions for p-Pb collisions at s NN = 5.02 and 8.16 TeV are compared with data from the Large Hadron Collider. Our approach combines a relativistic diffusion model with formulations based on quantum chromodynamics while utilizing numerical solutions of a Fokker-Planck equation to account for the shift and broadening of the fragmentation sources for particle-production with respect to the stopping (net-baryon) rapidity distributions. To represent the centrality dependence of charged-hadron production in asymmetric systems over a broad region of pseudorapidities, the consideration and precise modeling of the fragmentation sources, along with the central gluon-gluon source, is found to be essential. Specifically, this results in an inversion of the particle-production amplitude from backward to forward dominance when transitioning from central to peripheral collisions, in agreement with recent ATLAS and ALICE p-Pb data at s NN = 5.02 T e V . Published by the American Physical Society 2024

Schulz, Philipp (ORCID:0009000266642148)↗
DOE OSTI2025

Relaxation times for disoriented isospin condensates in high energy heavy ion collisions

Fluctuations between charged and neutral kaons measured by the ALICE Collaboration in Pb-Pb collisions at the CERN Large Hadron Collider (LHC) exceed conventional explanations. Previously it was shown that if the scalar condensate is accompanied by an electrically neutral isospin-1 field then the combination can produce large equilibrium fluctuations where ⟨$\overline{𝑢}$ ⁢𝑢⟩≠⟨ $\overline{𝑑}$𝑑⟩. Hadronizing strange and antistrange quarks might then strongly fluctuate between charged ($\overline{𝑢}$⁢𝑠 or $\overline{𝑠}$⁢𝑢) and neutral (𝑑⁢$\overline{𝑠}$⁢ or 𝑠⁢$\overline{𝑑}$) kaons. Here, we estimate the times for the condensates to achieve their equilibrium probability distributions within causal volumes in high energy heavy ion collisions. This is achieved by modeling the temperature dependence of the condensates, mesonic collective excitations, decay rates of the associated fields, and employing the Langevin and Fokker-Planck equations. Within this model, we find that the equilibration times are short compared with the expansion time, making disoriented isospin condensates a viable explanation for the anomalous fluctuations observed at the LHC.

Quark-gluon plasma↗
DOE OSTI2022

Direct Measurement of the Return Current Instability in a Laser-Produced Plasma

Measurements were made of the return current instability growth rate, demonstrating its concurrence with nonlocal transport. Thomson scattering was used to measure a maximum growth rate of 5.1 × 10 9 Hz, which was 3 times less than classical Spitzer-Härm theory predicts. Here, the measured plasma conditions indicate the heat flux was nonlocal, and Vlasov–Fokker-Planck simulations that account for nonlocality reproduce the measured growth rates. Furthermore, the threshold for the return current instability was measured (δ T = 0.017 ± 0.002) to be in good agreement with previous theoretical models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Shifting and Splitting of Resonance Lines due to Dynamical Friction in Plasmas

We report a quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Correlated fission fragment spin dynamics

Here, this study explores the role of nucleon exchange for the generation of the fission fragment angular momenta. For a number of typical fission cases, samples of 10 4 shape evolutions are generated by Langevin simulation and, subsequently, for each such evolution, the nucleon exchange transport theory previously developed for damped nuclear reactions is used to obtain the development of the fragment spin-spin distribution within the Fokker-Planck transport framework. The characteristic evolution of both parallel and perpendicular spin components is discussed. A common feature is that the rotational modes fall out of equilibrium before scission when the temperature rises rapidly while the concurrent shrinking of the neck suppresses further exchange. A number of fission observables are extracted from the event ensembles: the distribution of the magnitude of the fragment spin and its orientation relative to the fission axis, as well as the correlation between the two spins and the distribution of their opening angle. The dependence of these observables on the mass asymmetry is also examined.

fission↗
DOE OSTI2025

Bias-Variance Trade-Off in Physics-Informed Neural Networks with Randomized Smoothing for High-Dimensional PDEs

Physics-Informed Neural Networks (PINNs) have triggered a paradigm shift in scientific computing, leveraging mesh-free properties and robust approximation capabilities. While proving effective for low-dimensional partial differential equations (PDEs), the computational cost of PINNs remains a hurdle in high-dimensional scenarios. This is particularly pronounced when computing high-order and high-dimensional derivatives in the physics-informed loss. Randomized Smoothing PINN (RS-PINN) introduces Gaussian noise for stochastic smoothing of the original neural net model, enabling the use of Monte Carlo methods for derivative approximation, which eliminates the need for costly automatic differentiation. Despite its computational efficiency, especially in the approximation of high-dimensional derivatives, RS-PINN introduces biases in both loss and gradients, negatively impacting convergence, especially when coupled with stochastic gradient descent (SGD) algorithms. We present a comprehensive analysis of biases in RS-PINN, attributing them to the nonlinearity of the Mean Squared Error (MSE) loss as well as the intrinsic nonlinearity of the PDE itself. We propose tailored bias correction techniques, delineating their application based on the order of PDE nonlinearity. The derivation of an unbiased RS-PINN allows for a detailed examination of its advantages and disadvantages compared to the biased version. Specifically, the biased version has a lower variance and runs faster than the unbiased version, but it is less accurate due to the bias. To optimize the bias-variance trade-off, we combine the two approaches in a hybrid method that balances the rapid convergence of the biased version with the high accuracy of the unbiased version. In addition to methodological contributions, we present an enhanced implementation of RS-PINN. Extensive experiments on diverse high-dimensional PDEs, including Fokker-Planck, Hamilton-Jacobi-Bellman (HJB), viscous Burgers’, Allen-Cahn, and Sine-Gordon equations, illustrate the bias-variance trade-off and highlight the effectiveness of the hybrid RS-PINN. Empirical guidelines are provided for selecting biased, unbiased, or hybrid versions, depending on the dimensionality and nonlinearity of the specific PDE problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Score-Based Physics-Informed Neural Networks for High-Dimensional Fokker–Planck Equations

The Fokker-Planck (FP) equation is a foundational partial differential equation (PDE) in stochastic processes involving Brownian motions. However, the curse of dimensionality (CoD) poses a formidable challenge when dealing with high-dimensional FP equations. Although Monte Carlo simulation and (vanilla) Physics-Informed Neural Networks (PINNs) have shown the potential to tackle CoD, both methods exhibit significant numerical errors in high dimensions when dealing with the probability density function (PDF) associated with Brownian motion. The point-wise PDF values tend to decrease exponentially as dimensionality increases, surpassing the precision of numerical simulations and resulting in substantial errors. In addition, due to its massive sampling, Monte Carlo fails to offer fast sampling. Modeling the logarithm likelihood (LL) via vanilla PINNs transforms the FP equation into a notoriously difficult Hamilton-Jacobi-Bellman (HJB) equation, which is impractical for PINN learning, whose error grows rapidly with dimension. To this end, we propose a novel approach utilizing a score-based solver to fit the score function in stochastic differential equations (SDEs). The score function, defined as the gradient of the LL, plays a fundamental role in inferring LL and PDF and enables fast SDE sampling, offering an effective means to overcome the CoD. Three fitting methods, Score Matching (SM), Sliced Score Matching (SSM), and Score-PINN, are introduced, each contributing unique advantages in computational complexity, accuracy, and generality. The proposed score-based SDE solver operates in two stages: first, employing score matching or Score-PINN to acquire the score function; and second, solving the LL via an ordinary differential equation (ODE) using the obtained score function. Comparative evaluations across these methods showcase varying trade-offs. The proposed methodology is evaluated across diverse SDEs, including anisotropic Ornstein-Uhlenbeck processes, geometric Brownian motion, and Brownian motion with varying eigenspace. We also test various distributions, including Gaussian, Log-normal, Laplace, and Cauchy distributions. The numerical results demonstrate the score-based SDE solver’s stability, speed, and performance across different experimental settings, solidifying its potential as a solution to CoD for high-dimensional FP equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dataset from: Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species primitive moments are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, avoids unphysical temperatures and satisfies the H-theorem unlike an analogous Bhatnagar-Gross-Krook operator. Formulas for both a Cartesian velocity-space and a gyroaveraged operator are provided for use in Vlasov as well as long-wavelength gyrokinetic models. We present an algorithm for the discontinuous Galerkin discretization of this operator, and provide results from relaxation and Landau damping benchmarks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Formation and Redshift Evolution of Dark Matter Spikes

Dark matter density spikes forming around adiabatically growing black holes can dramatically enhance indirect and direct detection signals. Canonical predictions, however, assume a zero-mass seed in a purely dark matter environment and do not track the long-term dynamical impact of surrounding stars. We present a semi-analytic framework that first generalizes adiabatic spike formation to include finite seed masses, stellar cusps, and non-circular orbits, and then studies the subsequent cosmic evolution by solving coupled Fokker-Planck equations for the dark matter and stellar phase-space distributions, with a heating rate modulated by the cosmic star formation rate. Starting conservatively from canonical Gondolo-Silk spikes and marginalizing over astrophysical uncertainties, we find that stellar gravitational heating drives the inner slope towards $γ_χ\simeq 1.5$ within a few Gyrs (e.g by $z \lesssim 2$ for spikes formed at $z\simeq 10$), yielding overdensities two to four orders of magnitude below canonical expectations but still well above an NFW-like cusp. We provide redshift-dependent benchmarks for the column density and $J$-factor relevant to scattering, decay and annihilation signatures. Any robust interpretation of indirect dark matter signals from galactic nuclei must account for this evolution.

Herrera, Gonzalo [MIT, LNS; MIT, MKI; Harvard U.]↗
DOE OSTI2024

Electron tail suppression and effective collisionality due to synchrotron emission and absorption in mildly relativistic plasmas

Synchrotron radiation losses are a significant cause of concern for high-temperature aneutronic fusion reactions such as proton–Boron 11. The fact that radiation losses occur primarily in the high-energy tail, where the radiation itself has a substantial impact on the electron distribution, necessitates a self-consistent approach to modeling the diffusion and drag induced by synchrotron absorption and emission. Furthermore, an accurate model must account for the fact that the radiation emission spectrum is momentum-dependent, and the plasma opacity is frequency-dependent. Here, we present a simple Fokker–Planck operator, built on a newly solved-for blackbody synchrotron diffusion operator, which captures all relevant features of the synchrotron radiation. Focusing on magnetic mirror fusion plasmas, we show that significant suppression of the electron distribution occurs for relativistic values of the perpendicular electron momentum, which therefore emit much less radiation than predicted under the assumption of a Maxwell–Jüttner distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Velocity-space compression from Fermi acceleration with Lorentz scattering

The Fermi acceleration model describes how cosmic ray particles accelerate to great speeds by interacting with moving magnetic fields. In this work, we identify a variation of the model where light ions interact with a moving wall while undergoing pitch angle scattering through Coulomb collisions due to the presence of a heavier ionic species. The collisions introduce a stochastic component which adds complexity to the particle acceleration profile and sets it apart from collisionless Fermi acceleration models. The unusual effect captured by this simplified variation of Fermi acceleration is the nonconservation of phase space, with the possibility for a distribution of particles initially monotonically decreasing in energy to exhibit an energy peak upon compression. A peaked energy distribution might have interesting applications, such as to optimize fusion reactivity or to characterize astrophysical phenomena that exhibit nonthermal features.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Thermal conductivity of a laser plasma

Here, we present a model of the electron thermal conductivity of a laser-produced plasma. The model, supported by Vlasov–Fokker–Planck simulations, predicts that laser absorption reduces conductivity by forcing electrons out of a Maxwell–Boltzmann equilibrium, which results in depletion of both low-velocity bulk electrons and high-velocity tail electrons. We show that both the bulk and tail electrons approximately follow super-Gaussian distributions, but with distinct exponents that each depend on the laser intensity and wavelength through the parameter a = Zv$^{2}_{E}$/v$^{2}_{T}$. For a value of a = 0.5, tail depletion reduces the thermal conductivity to half its zero-intensity value. We present our results as simple analytic fits that can be readily implemented in any radiation-hydrodynamics code or used to correct the local limit of non-local conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Regulation of Alfvén Eigenmodes by Microturbulence in Fusion Plasmas

Global gyrokinetic simulations of mesoscale reversed shear Alfven eigenmodes (RSAE) excited by energetic particles (EP) in fusion plasmas find that RSAE amplitude and EP transport are much higher than experimental levels at nonlinear saturation, but quickly diminish to very low levels after the saturation when background microturbulence is artificially suppressed. In contrast, in simulations coupling micro-meso scales, the RSAE amplitude and EP transport decrease drastically at the initial saturation but later increases to the experimental levels in the quasi-steady state with bursty dynamics due to regulation by thermal ion temperature gradient (ITG) microturbulence. The quasi-steady state EP transport is larger for a stronger microturbulence. Here, the RSAE amplitude in the quasi-steady state ITG-RSAE turbulence from gyrokinetic simulations, for the first time, agrees very well with experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Absorption and Transport Effects Induced in Plasmas by the Interaction of Electrons with Laser Speckles

Here, we show that the ponderomotive force associated with laser speckles can scatter electrons in a laser-produced plasma in a manner similar to Coulomb scattering. Analytic expressions for the effective collision rates are given. The electron-speckle collisions become important at high laser intensity or during filamentation, affecting both long- and short-pulse laser intensity regimes. As an example, we find that the effective collision rate in the laser-overlap region of hohlraums on the National Ignition Facility is expected to exceed the Coulomb collision rate by 1 order of magnitude, leading to a fundamental change to the electron transport properties. At the high intensities characteristic of short-pulse laser-plasma interactions (I ≳ 10 17 W cm –2 ), the scattering is strong enough to cause the direct absorption of laser energy, generating hot electrons with energy scaling as E ≈ 1.44 (I/10 18 W cm –2 ) 1/2 MeV , close to experimentally observed results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗