Search NASA⌕ Search

SEARCH · Search NASA

Results for “Nonlinear diffusion acceleration”

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

Cluster for Research Excellence (CRE) for Accelerator Physics

The purpose of this Cooperative Research and Development Agreement (CRADA) was to establish a Cluster of Research Excellence (CRE) for Accelerator Physics and broaden the accelerator physics collaboration between Northern Illinois University (NIU) and Fermilab. The CRE was based in both Fermilab’s Accelerator Division and in NIU’s Department of Physics. A research agenda and strategy for the CRE for Accelerator Physics included specific research topics, experiments, and research facility upgrades. The scope of the research program covered two thrusts with subtopics in the Accelerator Science Research Program utilizing the infrastructure of FAST, IOTA, the Fermilab Accelerator Complex, and the Illinois Accelerator Research Center (IARC). Specific areas of proposed exploration included: 1.Tests of novel “nonlinear integrable” systems including beam self-fields; 2.Experimental tests of phase-space diffusion and chaos in intense beams; 3.Investigate technically feasible approaches towards next generation intenseneutrino facilities; 4.Understanding optical control of charged particle beams; 5.Exploitation of bright electron beams from FAST for novel applications; 6. Conduct research into accelerator-based ‘precision science’, such as Muon g-2, EDM, Mu2e, etc.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

A Monte Carlo Thermal Radiative Transfer Solver with Nonlinear Elimination

Here in this paper, we present a new Monte Carlo method for solving the thermal radiative transfer (TRT) equations via the method of nonlinear elimination (NLEM). This method is inspired by the previous application of NLEM to thermal radiation diffusion. Our approach, called diffusion accelerated Implicit Monte Carlo (DAIMC), is a hybrid technique which combines a Monte Carlo method for solving a purely-absorbing transport equation and a diffusion solution that accounts for effective scattering, or absorption–reemission. The method aims to improve the implicitness of the traditional implicit Monte Carlo (IMC) method. We derive DAIMC generally for 3D Cartesian geometries, but in this paper, we present results and analysis in 1D slab geometry. These preliminary results indicate that DAIMC implementations may provide more accurate and robust TRT solutions than IMC in certain test problems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

On the possibility of footprint compression with one lens in nonlinear accelerator lattice

Electromagnetic interaction of colliding beams along with other nonlinear fields often limits the beams' lifetimes and luminosities. Nonlinearities result in the spread of betatron frequencies (footprint) and, thus, may enhance dynamic diffusion of particles due to high order resonances. One of the possible ways to eliminate nonlinearities and overcome the corresponding difficulties is compensation of nonlinear forces, but, in practice, it is hardly possible to obtain exact linearity of the system. The compensation with a single nonlinear lens cannot cope with distributed nonlinearities, nonlinearities due to parasitic crossings, etc. Here, we present a method to compute parameters of nonlinear element (lens) that eliminates both the footprint and resonance strength without achieving full compensation.

43 PARTICLE ACCELERATORS↗
DOE OSTI2016

High-order solution methods for grey discrete ordinates thermal radiative transfer

This work presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation is accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Galactic cosmic-ray mediation of a spherical solar wind flow. 1: The steady state cold gas hydrodynamical approximation

Realistic models of the outer heliosphere should consider that the interstellar cosmic-ray pressure becomes comparable to pressures in the solar wind at distances more than 100 AU from the Sun. The cosmic-ray pressure dynamically affects solar wind flow through deceleration. This effect, which occurs over a scale length of the order of the effective diffusion length at large radial distances, has important implications for cosmic-ray modulation and acceleration. As a first step toward solution of this nonlinear problem, a steady state numerical model was developed for a relatively cold spherical solar wind flow which encounters the confining isotropic pressure of the surrounding Galactic medium. This pressure is assumed to be dominated by energetic particles (Galactic cosmic rays). The system of equations, which are solved self-consistently, includes the relevant hydrodynamical equations for the solar wind flow and the spherical cosmic-ray transport equation. To avoid the closure parameter problem of the two-fluid model, the latter equation is solved for the energy-dependent cosmic-ray distribution function.

Le Roux, J. A.↗
DOE OSTI2026

Searching for the Most Harmful Field Errors in the HSR IR Superconducting Magnets

In this project, we improve beam stability for the Electron-Ion Collider. Magnetic field errors can reduce beam stability, making it essential to identify the field errors that have the greatest impact on accelerator performance. However, this is particularly challenging because beam stability depends on the complex interactions of many magnetic field errors, resulting in a high-dimensional and nonlinear optimization problem. We determine which field errors are the most influential for the large physical aperture superconducting magnet B2PF, a critical magnet in the Interaction Region (IR) in the Hadron Storage Ring (HSR). We complete and analyze nearly 30,000 simulations on the Brookhaven National Laboratory Linux Cluster by varying 18 nonlinear magnetic field errors. We evaluate beam stability using the dynamic aperture and the tune diffusion. We identify the field errors that most strongly influence beam stability and establish quantitative field error tolerances that improve accelerator performance.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Rarefaction-driven Rayleigh–Taylor instability. Part 2. Experiments and simulations in the nonlinear regime

Experiments and large eddy simulation (LES) were performed to study the development of the Rayleigh–Taylor instability into the saturated, nonlinear regime, produced between two gases accelerated by a rarefaction wave. Single-mode two-dimensional, and single-mode three-dimensional initial perturbations were introduced on the diffuse interface between the two gases prior to acceleration. The rarefaction wave imparts a non-constant acceleration, and a time decreasing Atwood number, $A=(\unicode[STIX]{x1D70C}_{2}-\unicode[STIX]{x1D70C}_{1})/(\unicode[STIX]{x1D70C}_{2}+\unicode[STIX]{x1D70C}_{1})$ , where $\unicode[STIX]{x1D70C}_{2}$ and $\unicode[STIX]{x1D70C}_{1}$ are the densities of the heavy and light gas, respectively. Experiments and simulations are presented for initial Atwood numbers of $A=0.49$ , $A=0.63$ , $A=0.82$ and $A=0.94$ . Nominally two-dimensional (2-D) experiments (initiated with nearly 2-D perturbations) and 2-D simulations are observed to approach an intermediate-time velocity plateau that is in disagreement with the late-time velocity obtained from the incompressible model of Goncharov ( Phys. Rev. Lett. , vol. 88, 2002, 134502). Reacceleration from an intermediate velocity is observed for 2-D bubbles in large wavenumber, $k=2\unicode[STIX]{x03C0}/\unicode[STIX]{x1D706}=0.247~\text{mm}^{-1}$ , experiments and simulations, where $\unicode[STIX]{x1D706}$ is the wavelength of the initial perturbation. At moderate Atwood numbers, the bubble and spike velocities approach larger values than those predicted by Goncharov’s model. These late-time velocity trends are predicted well by numerical simulations using the LLNL Miranda code, and by the 2009 model of Mikaelian ( Phys. Fluids. , vol. 21, 2009, 024103) that extends Layzer type models to variable acceleration and density. Large Atwood number experiments show a delayed roll up, and exhibit a free-fall like behaviour. Finally, experiments initiated with three-dimensional perturbations tend to agree better with models and a simulation using the LLNL Ares code initiated with an axisymmetric rather than Cartesian symmetry.

42 ENGINEERING↗
DOE OSTI2024

Theory of Beam Echoes

We develop the theory of beam echoes in circular accelerators under several different conditions. We derive detailed expressions for the echo amplitude and pulse width with nonlinear quadrupole and dipole kicks, first without and then with momentum spread. We use the theory with the linearized dipole and quadrupole kicks to solve the diffusion equation for different dependencies of the diffusion coefficient on the action. We then consider the use of multiple quadrupole kicks to increase the maximum echo amplitude. We have extended these calculations partially to the 2D case and we also have partial results for longitudinal echoes.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Phase-space entropy cascade and irreversibility of stochastic heating in nearly collisionless plasma turbulence

We consider a nearly collisionless plasma consisting of a species of “test particles” in one spatial and one velocity dimension, stirred by an externally imposed stochastic electric field—a kinetic analog of the Kraichnan model of passive advection. The mean effect on the particle distribution function is turbulent diffusion in velocity space—known as stochastic heating. Accompanying this heating is the generation of fine-scale structure in the distribution function, which we characterize with the collisionless (Casimir) invariant C 2 ∝ ∫ ∫ d x d v 〈 f 2 〉 —a quantity that here plays the role of (negative) entropy of the distribution function. We find that C 2 is transferred from large scales to small scales in both position and velocity space via a phase-space cascade enabled by both particle streaming and nonlinear interactions between particles and the stochastic electric field. We compute the steady-state fluxes and spectrum of C 2 in Fourier space, with k and s denoting spatial and velocity wave numbers, respectively. In our model, the nonlinearity in the evolution equation for the spectrum turns into a fractional Laplacian operator in k space, leading to anomalous diffusion. Whereas even the linear phase mixing alone would lead to a constant flux of C 2 to high s (towards the collisional dissipation range) at every k , the nonlinearity accelerates this cascade by intertwining velocity and position space so that the flux of C 2 is to both high k and high s simultaneously. Integrating over velocity (spatial) wave numbers, the k -space ( s -space) flux of C 2 is constant down to a dissipation length (velocity) scale that tends to zero as the collision frequency does, even though the rate of collisional dissipation remains finite. The resulting spectrum in the inertial range is a self-similar function in the ( k , s ) plane, with power-law asymptotics at large k and s . Our model is fully analytically solvable, but the asymptotic scalings of the spectrum can also be found via a simple phenomenological theory whose key assumption is that the cascade is governed by a “critical balance” in phase space between the linear and nonlinear timescales. We argue that stochastic heating is made irreversible by this entropy cascade and that, while collisional dissipation accessed via phase mixing occurs only at small spatial scales rather than at every scale as it would in a linear system, the cascade makes phase mixing even more effective overall in the nonlinear regime than in the linear one. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Complete quasilinear model for the acceleration-driven lower hybrid drift instability and a computational assessment of its validity

A complete quasilinear model is derived for the electrostatic acceleration-driven lower hybrid drift instability in a uniform two-species low-beta plasma in which current is perpendicular to the background magnetic field. The model consists of coupled nonlinear velocity space diffusion equations for the volume-averaged ion and electron distribution functions. Each species' diffusion coefficient depends on a time-evolving spectral density of the electric-field energy per unit volume and a time-evolving dispersion relation. The dispersion relation is expressed analytically in integral form without the use of asymptotic limits and applies to arbitrary distribution functions, so long as they can be expressed as a function of one velocity coordinate, e.g., f⁡(vy) or f⁡(v⊥). The quasilinear model conserves energy and is complete in that it fully describes the evolution of the distribution functions, including resonant and nonresonant particle-wave interactions, while accounting for distribution-function-dependent mixed-complex frequencies. Further, the quasilinear diffusion model is solved numerically and self-consistently using a Crank-Nicolson temporal discretization and a second-order finite-volume velocity-space discretization. Numerical solutions are compared to nonlinear fourth-order accurate continuum kinetic Vlasov-Poisson simulations. Evolution of electric-field energy, growth rates, distribution functions, and diffusion coefficients are shown to be in agreement with Vlasov simulations. The quasilinear model is shown to predict anomalous transport terms, like resistivity and heating, to within a factor of order unity. Discrepancies between the quasilinear model and Vlasov simulations are assessed and attributed primarily to lack of damping in the quasilinear description and to the use of unperturbed-orbit susceptibilities in the linear theory dispersion relation. The results illuminate the predictive accuracy of the quasilinear model, place approximate bounds on its validity, and provide much needed vetting of quasilinear theory's ability to predict the nonlinear state of a microturbulent plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Investigation B: Laminar confined coaxial entrance flow with heat generation

The results of a parametric study on the entrance flow region in a gas core nuclear reactor are presented. The physical system is modeled as laminar confined, coaxial flow with heat generation in the inner fluid. The governing equations include the boundary layer approximations and the assumptions of only radial radiative transport of energy represented as an energy diffusion term. The Von Mises transformation and a zeta transformation are used to transform the equations into nonlinear nonhomogeneous convective-diffusion equations. A unique combination of forward and backward difference equations which yields accurate results at moderate computational times, is used in the numerical method. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Bobba, G. K. M.↗
NASA NTRS1979

Magnetospheric reconnection, substorms, and energetic particle acceleration

The steady state reconnection model of the terrestrial magnetosphere predicts a maximum potential drop of about 100 kV across the tail. During substorms particles are accelerated to energies above 1 MeV. At substorm onset, large inductive emfs may be generated by explosive tearing mode reconnection which is driven nonlinearly unstable by the solar wind and convection stresses on the tail plasma sheet. In the inner magnetosphere, energetic particles are also produced by stochastic wave turbulent acceleration and by convection driven inward radial diffusion.

Coroniti, F. V.↗
NASA NTRS1990

On neglect of nonlinear momentum terms in solar nebula accretion disk models

The standard model for accretion disks (e.g., Pringle, 1981; Lin and Papaloizou, 1985) parameterizes the radial eddy angular momentum flux which is quadratically nonlinear in velocity, solely by a viscous diffusion process which is linear in velocity. This paper shows that the eddy angular momentum flux, in addition to a gradient dissipative (viscous) flux, involves a countergradient accelerating flux. Therefore, the momentum parameterization in the standard model is in general unjustified. The formalism of the parameterization can be justified only by recognizing that the eddy viscosity used in the standard model is in fact a pseudoviscosity which includes the effects of countergradient fluxes, and which therefore may have a magnitude less than the eddy diffusivity K for mixing trace constituents in the disk.

Prinn, Ronald G.↗
NASA NTRS1989

Fokker-Planck electron diffusion caused by an obliquely propagating electromagnetic wave packet of narrow bandwidth

The relativistic motion of electrons in an intense electromagnetic wave packet propagating obliquely to a uniform magnetic field is analytically studied on the basis of the Fokker-Planck-Kolmogorov (FPK) approach. The wavepacket consists of circularly polarized electron-cyclotron waves. The dynamical system in question is shown to be reducible to one with three degrees of freedom. Within the framework of the Hamiltonian analysis the nonlinear diffusion tensor is derived, and it is shown that this tensor can be separated into zeroth-, first-, and second-order parts with respect to the relative bandwidth. The zeroth-order part describes diffusive acceleration along lines of constant unperturbed Hamiltonian. The second-order part, which corresponds to the longest time scale, describes diffusion across those lines. A possible transport theory is outlined on the basis of this separation of the time scales.

Hizanidis, Kyriakos↗
DOE OSTI2018

Suppression of Tearing Modes by Radio Frequency Current Condensation [Suppression of Tearing Modes by RF Current Condensation]

Currents driven by radio frequency (rf) waves in the interior of magnetic islands can stabilize deleterious tearing modes in tokamaks. Present analyses of stabilization assume that the local electron acceleration is unaffected by the presence of the island. However, the power deposition and electron acceleration are sensitive to the perturbation of the temperature. The nonlinear feedback on the power deposition in the island increases the temperature perturbation, and can lead to a bifurcation of the solution to the steady-state heat diffusion equation. The combination of the nonlinearly enhanced temperature perturbation with the rf current drive sensitivity to the temperature leads to an rf current condensation effect, which can increase the efficiency of rf current drive stabilization and reduce its sensitivity to radial misalignment of the ray trajectories. As a result, the threshold for the effect is in a regime that has been encountered in experiments, and will likely be encountered in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Gas Jet Profile Monitor for Use in IOTA Proton Beam

The Integrable Optics Test Accelerator (IOTA) at the Fermilab Accelerator Science and Technology (FAST) Facility will attempt to demonstrate novel techniques for high current accelerators. An electron beam of 150 MeV/c momentum and a proton beam of 70 MeV/c momentum will be used for a number of experiments, including nonlinear focusing integrable optics, space charge compensation, and optical stochastic cooling. The low energy proton beam will provide a space charge dominated regime with which to investigate Hamiltonian diffusion and chaos. A non-invasive beam instrumentation device is needed to study emittance evolution and halo formation without destroying the beam. A supersonic gas jet curtain can be injected perpendicular to the proton beam, which will preserve the integrity of the beam and provide a two-dimensional turn-by-turn profile measurement. Currently, IOTA and the gas jet monitor are being built and commissioned at Fermilab.

43 PARTICLE ACCELERATORS↗
NASA NTRSDate not supplied

Assessment of Edge-Based Viscous Method for Corner-Flow Solutions on Graphics Processing Units

A highly efficient, edge-based viscous (EBV) discretization method has been recently implemented in a practical, unstructured-grid, node-centered, finite-volume flow solver and evaluated for Reynolds-averaged Navier-Stokes (RANS) formulations. In comparison to a well-established cell-based viscous (CBV) method, the EBV method has demonstrated multifold acceleration of all viscous-kernel computations on general unstructured mixed-element grids. The viscous kernels include evaluation of viscous fluxes, diffusion terms in turbulence models, and the corresponding Jacobian terms. In this paper, an EBV implementation of a nonlinear extension of the Spalart-Allmaras turbulence model, SA-neg-QCR2000, is presented and verified. The SA-neg-QCR2000 model is used for simulating turbulent corner flows. Previously reported EBV computations have been conducted on traditional computing architectures based on central processing units (CPU). This paper assesses benefits of the EBV method on modern high-performance computing architectures based on graphics processing units (GPU). The GPU implementations of the CBV and EBV methods are verified by comparing solutions and iterative convergence with those observed in CPU computations on the same grids. A comprehensive assessment of the EBV speedup on CPU and GPU architectures is presented for established benchmark corner flows, namely, a supersonic flow through a long square duct and a subsonic flow around a NASA juncture flow model.

CFD↗