Search NASA⌕ Search

SEARCH · Search NASA

Results for “Radiation damping”

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 OSTI2025

Radiative damping of toroidal Alfvén eigenmode in low-shear plasmas

Instabilities of Alfvén eigenmodes (AEs) are of significant concern because they can enhance the cross-field transport of fusion-born alpha particles beyond the neoclassical level in magnetic fusion plasmas. The threshold value of alpha-particle pressure for exciting AEs depends critically on the damping rate of AEs. The damping mechanisms include kinetic damping due to interactions with thermal particles, continuum damping due to AE frequency crossing Alfvén continuum, and radiative damping due to emitting kinetic Alfvén waves (KAWs). The radiative damping is substantial and can even prevail in high-temperature burning plasmas [1]. We revisit the radiative damping analytic theory for TAE in plasmas with low positive magnetic shear, considering TAE with an eigenfrequency near the bottom of TAE-gap and with poloidal harmonics of the same sign (even TAE). In contrast to earlier papers, we provide the damping calculations in real space rather than Fourier space. This approach is straightforward technically and more enlightening from a physics standpoint for benchmarking numerical calculations of radiative damping. The parametric dependence of the resulting damping rate agrees with that of Refs. [2-5], but it has a smaller numerical factor in front of it.

Alpha-particle driven instability↗
DOE OSTI2024

Linear gyrokinetic simulations of toroidal Alfvén eigenmodes in the Mega-Amp Spherical Tokamak

Linear gyrokinetic (GK) simulations using the Gyrokinetic Toroidal Code (GTC) have been performed to investigate Toroidicity-driven Alfvén Eigenmodes (TAEs) driven by the neutral beam injection (NBI) induced fast ions in the Mega-Amp Spherical Tokamak (MAST) to identify the non-perturbative and kinetic effects of thermal plasma. A specific TAE in MAST discharge 26887, with an on-axis NBI power of approximately 1.5 MW and plasma current around 800 kA, exhibited frequency chirping, and the tangential soft x-ray camera array resolved the radial mode structure peaked near |q|=1.5. Various excitation methods were used in the GTC linear simulations, illustrating this code's capability to realistically represent the mechanisms and behaviors of fast ion-driven TAEs in spherical tokamaks. The radial structures from these GK simulations closely match measurements and calculations performed using the NOVA ideal MHD code, though with the frequencies approximately 10 kHz lower, likely due to various kinetic and non-perturbative effects. The simulations measured the damping rates due to continuum damping, radiative damping, and ion Landau damping, revealing that ion Landau damping has the most significant contribution to the total damping rate of the TAE. A comparison of growth rates of TAEs excited by fast ion Maxwellian and slowing-down distributions shows that the TAEs excited by a fast ion anisotropic pitch distribution (as part of the slowing-down distributions) are more unstable compared to those excited by a Maxwellian distribution with an equivalent fast ion beta. This shows that the use of fast ion anisotropy alters the number of fast ions to be in shear Alfvén resonance, and hence, it can greatly affect the stability of TAEs. These tests can be performed with the GTC but impossible with ideal MHD simulations, highlighting the necessity of kinetic simulations such as the GTC for a precise prediction of the TAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Dual Energy Electron Storage Ring Cooler Design for Relativistic Ion Beams

Collider experiments demand small beam emittances in order to achieve high luminosity. For light particles such as electrons, there exists a natural synchrotron radiation damping resulting in low emittance beams at equilibrium. In the case of heavy particle beams such as proton or ion beams, there is no significant synchrotron radiation damping effect and some cooling mechanism is needed to get to low emittance beams. A dual energy storage ring cooler is a novel concept proposed to cool hadron beams at higher energies. The design consists of two rings: a low energy ring and a high energy ring connected by the energy recovering linac (ERL) that provides the necessary energy difference. The low energy ring has an electron cooler where the cooling interaction takes place between the electron and hadron beams whereas electron beam going through the high energy ring undergoes synchrotron radiation damping. In this document, we present a possible design of a dual energy storage ring-based electron cooler for high energy proton beam cooling suitable for the Electron Ion Collider (EIC) to be built at Brookhaven National laboratory (BNL). A special feature of a dual energy storage ring cooler design is that the electron beam energy in the low energy section must be tuned to match the hadron beam velocity in the cooling section, and the electron beam energy in the high energy section is chosen to provide an adequate synchrotron radiation damping. In addition, the ring optics design is carried out considering chromaticity correction, dynamic aperture, momentum aperture, beam lifetime, radiation damping and the intrabeam scattering effect. Finally, the cooling performance is simulated for 275 GeV proton beams at the top energy of the EIC. In such a collider it is desirable to cool the hadron beams to balance the emittance growth rates due to intra-beam scattering (IBS) to maintain a high luminosity during the collision runs. In such a situation, a dual energy storage ring based cooler could be a practical approach to be used for hadron beam cooling.

Dhital, Bhawin↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Stabilizing dynamic subsea power cables using Bi-stable nonlinear energy sinks

This study investigates vibration mitigation of dynamic subsea cables through passive bi-stable nonlinear energy sinks (B-NESs). These devices suppress vibration energy in a broadband way, and can be regarded as extensions of classical linear tuned mass dampers (TMDs) which are narrowband devices. Through the open-source MoorDyn library, we simulated the vibrations of a vertical subsea cable equipped with a set of B-NESs. Multi-objective optimization was performed to detect the B-NES parameters for optimal mitigation of the cable vibrations. Advanced signal processing verified the efficacy of the optimized B-NESs not only to rapidly absorb and locally dissipate vibration energy, but also to nonlinearly scatter vibration energy from low to high frequencies within the cable itself. This last feature is especially beneficial for vibration mitigation of the undersea cable, as at higher frequencies the cable vibrations exhibit drastically reduced amplitudes and are more effectively dissipated by inherent structural damping and hydrodynamic radiation damping. This contrasts with traditional TMDs which can mitigate vibration only at a single frequency. Furthermore, our robustness study confirms the B-NES's effectiveness under even varying environmental conditions. Overall, the B-NES's capacity for broadband vibration mitigation renders it a promising retrofit solution for improving the performance and operational safety of dynamic power cables in offshore wind farms and other marine applications.

17 WIND ENERGY↗
DOE OSTI2024

Isotope effects and Alfvén eigenmode stability in JET H, D, T, DT, and He plasmas

While much about Alfvén eigenmode (AE) stability has been explored in previous and current tokamaks, open questions remain for future burning plasma experiments, especially regarding exact stability threshold conditions and related isotope effects; the latter, of course, requiring good knowledge of the plasma ion composition. In the JET tokamak, eight in-vessel antennas actively excite stable AEs, from which their frequencies, toroidal mode numbers, and net damping rates are assessed. The effective ion mass can also be inferred using measurements of the plasma density and magnetic geometry. Thousands of AE stability measurements have been collected by the Alfvén Eigenmode Active Diagnostic in hundreds of JET plasmas during the recent Hydrogen, Deuterium, Tritium, DT, and Helium-4 campaigns. In this novel AE stability database, spanning all four main ion species, damping is observed to decrease with increasing Hydrogenic mass, but increase for Helium, a trend consistent with radiative damping as the dominant damping mechanism. These data are important for confident predictions of AE stability in both non-nuclear (H/He) and nuclear (D/T) operations in future devices. In particular, if radiative damping plays a significant role in overall stability, some AEs could be more easily destabilized in D/T plasmas than their H/He reference pulses, even before considering fast ion and alpha particle drive. Active MHD spectroscopy is also employed on select HD, HT, and DT plasmas to infer the effective ion mass, thereby closing the loop on isotope analysis and demonstrating a complementary method to typical diagnosis of the isotope ratio.

Alfvén eigenmodes↗
DOE OSTI2020

Ring-based electron cooler for EIC

An Electron Ion Collider (EIC) at Brookhaven National Laboratory (BNL) is being proposed as a new discovery machine for the nuclear physics and quantum chromodynamics. The hadron beam cooling plays an important role in the EIC machine to achieve its physics goals. The most challenging is cooling of protons at the highest energy in the EIC. In this paper, we present a possible design of a ring-based electron cooler for the high energy hadron beam cooling. In proposed approach, the electrons will cool the hadrons while being cooled themselves by radiation damping in the storage ring. For design of the cooler using storage ring approach several aspects becomes very important, including electron ring optics design, chromaticity correction, calculate the dynamic aperture, radiation damping, quantum excitation, and intrabeam scattering (IBS). In addition, such effects as beam-beam scattering (BBS) due to interaction of electrons with hadrons become of special concern. In this paper, we discuss and take into design all the effects above, as well as introduce generalized treatment of the BBS and 3-D quantum lifetime. A special feature of our design is an effective use of dispersion in the cooling section, both for the ions and electrons, to redistribute the cooling rate between the longitudinal and horizontal planes. Finally, the cooling performance is simulated for proton beam at the top energy of the EIC. Our conclusion is that such ring-based cooler could be a feasible approach to provide required parameters of hadron beam at the top energy of 275 GeV for the EIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Investigation of Theoretical Solutions to a Bottom-Raised Oscillating Surge Wave Energy Converter (OSWEC) Through Experimental and Parametric Studies

Experiments were conducted on a wave tank model of a bottom raised oscillating surge wave energy converter (OSWEC) model in regular waves. The OSWEC model shape was a thin rectangular flap, which was allowed to pitch in response to incident waves about a hinge located at the intersection of the flap and the top of the supporting foundation. Torsion springs were added to the hinge in order to position the pitch natural frequency at the center of the wave frequency range of the wave maker. The flap motion as well as the loads at the base of the foundation were measured. The OSWEC was modeled analytically using elliptic functions in order to obtain closed form expressions for added mass and radiation damping coefficients, along with the excitation force and torque. These formulations were derived and reported in a previous publication by the authors. While analytical predictions of the foundation loads agree very well with experiments, large discrepancies are seen in the pitch response close to resonance. These differences are analyzed by conducting a sensitivity study, in which system parameters, including damping and added mass values, are varied. The likely contributors to the differences between predictions and experiments are attributed to tank reflections, standing waves that can occur in long, narrow wave tanks, as well as the thin plate assumption employed in the analytical approach.

analytical↗
DOE OSTI2026

Adiabatic fast passage spin manipulation measurements in solid polarized targets

Adiabatic fast passage (AFP) is a rapid method for reversing nuclear polarization and manipulating spin populations in polarized solid targets, avoiding the long repolarization times associated with dynamic nuclear polarization (DNP). We report AFP measurements in a 5 T, 1 K polarized-target system for irradiated 15 NH 3 , irradiated 14 ND 3 , and butanol-based materials prepared either with TEMPO doping or by irradiation. We also present a joint manipulated-lineshape analysis for spin-1 targets and demonstrate that vector and tensor polarizations can be extracted from AFP-manipulated deuteron NMR spectra even when the populations are not described by a single Boltzmann spin temperature. In conclusion, we report a reproducible polarization- and direction-dependent AFP response in a large irradiated 15 NH 3 sample. These ammonia results are presented as empirical observations under the specific sample–coil conditions of the experiment, with possible circuit-mediated mechanisms such as radiation damping or superradiant behavior discussed but not assigned as a definitive cause.

Adiabatic fast passage↗
DOE OSTI2021

Optimizing Power Generation of a Bottomraised Oscillating Surge Wave Energy Converter Using a Theoretical Model

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2021

Optimizing Power Generation of a Bottom-Raised Oscillating Surge Wave Energy Converter Using a Theoretical Model: Preprint

Preliminary sizing of an oscillating surge wave energy converter (OSWEC) is an iterative process that relies on knowledge of the relevant hydrodynamic coefficients for a given geometry. Often, the complex definition of the device geometry requires coefficients to be obtained through experiments or numerical boundary element solvers such as WAMIT. These techniques demand significant user and computational effort, therefore inhibiting the fine scale parametrization of object dimensions. In this study, a theoretical model, originally presented in Michele et al. (2016), is developed and demonstrated to efficiently optimize the power production for an OSWEC device (subjected to certain environmental conditions) with variations in device widths, heights, and distances from the seabed. Assuming negligible device thickness, the OSWEC motions are modelled as a bottom-raised two-dimensional flap in regular waves using potential flow theory formulated in elliptical coordinates. The solutions to this diffraction-radiation problem are obtained using Mathieu functions with appropriate boundary conditions. The resulting potentials are then used to derive frequency-dependent expressions for the added mass and radiation damping coefficients, along with wave excitation magnitude in surge, pitch, and coupled surge-pitch motions. Good agreement in hydrodynamic coefficient curves is shown between the theoretical model and the numerical results obtained from the boundary element-based program WAMIT. The theoretical model is then employed to maximize the time-averaged output power while maintaining or reducing the hinge reaction force, with variations in device dimensions, wave frequency and amplitude.

ACE↗
DOE OSTI2022

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗
DOE OSTI2022

Simultaneous measurements of unstable and stable Alfvén eigenmodes in JET

Here, we report the novel experimental observation of both unstable and stable toroidicity-induced Alfvén eigenmodes (TAEs) measured simultaneously in a JET tokamak plasma. The three-ion-heating scheme (D-DNBI-3He) is employed to accelerate deuterons to MeV energies, thereby destabilizing TAEs with toroidal mode numbers n = 3–5, each decreasing in mode amplitude. At the same time, the Alfvén eigenmode active diagnostic resonantly excites a stable n = 6 TAE with total normalized damping rate –γ/ω 0 ≈ 1%–4%. Hybrid kinetic-MHD modeling with codes NOVA-K and MEGA both find eigenmodes with similar frequencies, mode structures, and radial locations as in experiment. NOVA-K demonstrates good agreement with the n = 3, 4, and 6 TAEs, matching the damping rate of the n = 6 mode within uncertainties and identifying radiative damping as the dominant contribution. Improved agreement is found with MEGA for all modes: the unstable n = 3–5 and stable n = 2, 6 modes, with the latter two stabilized by higher intrinsic damping and lower fast ion drive, respectively. While some discrepancies remain to be resolved, this unique validation effort gives us confidence in TAE stability predictions for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A novel measurement of marginal Alfvén eigenmode stability during high power auxiliary heating in JET

The interaction of Alfvén eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25–250 kHz and toroidal mode numbers |n| < 20. During the 2019–2020 deuterium campaign, almost 7500 resonances and their frequencies f 0 , net damping rates γ < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient’s effect on accessibility and ELM-related noise’s impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ion cyclotron resonance heating and neutral beam injection) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f 0 , γ, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A novel measurement of marginal Alfven Eigenmode stability during high power auxiliary heating in JET

The interaction of Alfven Eigenmodes (AEs) and energetic particles is one of many important factors determining the success of future tokamaks. In JET, eight in-vessel antennas were installed to actively probe stable AEs with frequencies ranging 25-250 kHz and toroidal mode numbers |n| < 20. During the 2019-2020 deuterium campaign, almost 7500 resonances and their frequencies f0, net damping rates \gamma < 0, and toroidal mode numbers were measured in almost 800 plasma discharges. From a statistical analysis of this database, continuum and radiative damping are inferred to increase with edge safety factor, edge magnetic shear, and when including non-ideal effects. Both stable AE observations and their associated damping rates are found to decrease with |n|. Active antenna excitation is also found to be ineffective in H-mode as opposed to L-mode; this is likely due to the increased edge density gradient's effect on accessibility and ELM-related noise's impact on mode identification. A novel measurement is reported of a marginally stable, edge-localized Ellipticity-induced AE probed by the antennas during high-power auxiliary heating (ICRH and NBI) up to 25 MW. NOVA-K kinetic-MHD simulations show good agreement with experimental measurements of f0, \gamma, and n, indicating the dominance of continuum and electron Landau damping in this case. Similar experimental and computational studies are planned for the recent hydrogen and ongoing tritium campaigns, in preparation for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Performance characterization and modeling of an oscillating surge wave energy converter

Abstract Testing wave energy converters in the ocean could be expensive and complex, which necessitates the use of numerical modeling. However, accurately modeling the response of wave energy converters with high-fidelity simulations can be computationally intensive in the design stage where different configurations must be considered. Reduced-order models based on simplified equations of motion can be very useful in the design, optimization, or control of wave energy converters. Given the complex dynamics of wave energy converters, accurate representation, and evaluation of relative contributions by different forces are required. This effort is concerned with a performance characterization of the hydrodynamic response of an oscillating surge wave energy converter that is based on a reduced-order model. A state-space model is used to represent the radiation damping term. Morison’s representation of unsteady forces is used to account for the nonlinear damping. Wave tank tests are performed to validate simulations. A free response simulation is used to determine the coefficients of the state-space model. Torque-forced simulations are used to identify the coefficients of the nonlinear damping term for different amplitudes and wave frequencies. The impact of varying these coefficients on the response is investigated. An assessment of the capability of the model in predicting the hydrodynamic response under irregular forcing is performed. The results show that the maximum error is 3% when compared with high-fidelity simulations. It is determined that the nonlinear damping is proportional to the torque amplitude and its effects are more pronounced as the amplitude of the flap oscillations increases.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗
DOE OSTI2021

A Feynman-Kac based numerical method for the exit time probability of a class of transport problems

The exit time probability, which gives the likelihood that an initial condition leaves a prescribed region of the phase space of a dynamical system at, or before, a given time, is arguably one of the most natural and important transport problems. In this work, we present an accurate and efficient numerical method for computing this probability for systems described by non-autonomous (time-dependent) stochastic differential equations (SDEs) or their equivalent Fokker-Planck partial differential equations. The method is based on the direct approximation of the Feynman-Kac formula that establishes a link between the adjoint Fokker-Planck equation and the forward SDE. The Feynman-Kac formula is approximated using the Gauss-Hermite quadrature rules and piecewise cubic Hermite interpolating polynomials, and a GPU accelerated matrix representation is used to compute the entire time evolution of the exit time probability using a single pass of the algorithm. The method is unconditionally stable, exhibits second order convergence in space, first order convergence in time, and it is straightforward to parallelize. Applications are presented to the advection diffusion of a passive tracer in a fluid flow exhibiting chaotic advection, and to the runaway acceleration of electrons in a plasma in the presence of an electric field, collisions, and radiation damping. Benchmarks against analytical solutions as well as comparisons with explicit and implicit finite difference standard methods for the adjoint Fokker-Planck equation are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗