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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. Calculating the stochastic diffusion of electrons due to a spectrum of whistler modes is presented. The parametric dependence of the diffusion coefficients on the plasma particle density, magnitude of wave field, and the strength of magnetic field is studies. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
DOE OSTI2017

Two species drag/diffusion model for energetic particle driven modes

A nonlinear bump-on-tail model for the growth and saturation of energetic particle driven plasma waves has been extended to include two populations of fast particles—one dominated by dynamical friction at the resonance and the other by velocity space diffusion. The resulting temporal evolution of the wave amplitude and frequency depends on the relative weight of the two populations. The two species model is applied to burning plasma with drag-dominated alpha particles and diffusion-dominated ICRH accelerated minority ions, showing the stabilization of bursting modes. Further, the model also suggests an explanation for the recent observations on the TJ-II stellarator, where Alfvén Eigenmodes transition between steady state and bursting as the magnetic configuration varied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Two species drag/diffusion model for energetic particle driven modes

A nonlinear bump-on-tail model for the growth and saturation of energetic particle driven plasma waves has been extended to include two populations of fast particles—one dominated by dynamical friction at the resonance and the other by velocity space diffusion. The resulting temporal evolution of the wave amplitude and frequency depends on the relative weight of the two populations. The two species model is applied to burning plasma with drag-dominated alpha particles and diffusion-dominated ICRH accelerated minority ions, showing the stabilization of bursting modes. The model also suggests an explanation for the recent observations on the TJ-II stellarator, where Alfv'en Eigenmodes transition between steady state and bursting as the magnetic configuration varied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1988

Dynamics of two fluids under periodic acceleration

The evolution of the interface between two fluids confined in a rectangular cavity is investigated numerically to predict its transient behavior. These computations address mixing characteristics of fluids under microgravity conditions, particularly g-jitter conditions resulting from aircraft vibration or crew motion, and have applications in solution crystal growth. The two-dimensional formulation employs the Boussinesq approximation and treats the mixing of two fluids as an initial value problem with a prescribed concentration field. For fluid mixing of practical applications inside a cavity, it is shown that nonlinear convective transport can dominate over both viscous and molecular diffusion. However, viscous diffusion can become important for low Reynolds number or certain cavity sizes. A stacking phenomenon is shown to occur for aspect ratios (Ar) approximately between 5 to 10. For aspect ratios in the neighborhood of 0.1 to 0.2, inner cells evolve at the interface with increasing time. For a square cavity of Ar = 1, chaotic mixing of the fluid appears to occur for a Stokes-Reynolds number greater than about 5.

Duval, Walter M. B.↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Effects of ion two-stream instability on auroral ion heating

The ion two-stream instability develops between upflowing H(+) and O(+) in the auroral acceleration region, if the two species are accelerated to different velocities in a semistationary electrostatic potential. Parallel propagating modes are unstable only for small relative drifts, while oblique modes remain unstable at larger drifts and, consequently, higher altitudes. The effects of oblique modes on ion heating are investigated, and it is found that O(+) heats primarily by trapping in the direction of wave propagation during the early, coherent phase of wave growth in the temporal evolution problem. After undergoing trapping along the magnetic field direction, H(+) heats primarily by quasi-linear diffusion when the spectrum broadens nonlinearly. The parallel ion distributions evolve by formation of a high-energy (low-energy) tail on oxygen (hydrogen) along with bulk slowing of hydrogen, which results in more energetic oxygen than hydrogen. Total heating is greatest in the direction of wave propagation, and this effect is greater for hydrogen than for oxygen.

Roth, I.↗
DOE OSTI2021

Nonlinear dynamics of proton beams with hollow electron lens in the CERN high-luminosity LHC

The design stored beam energy in the CERN high-luminosity large hadron collider (HL-LHC) upgrade is about 700 MJ, with about 36 MJ in the beam tails, according to estimates based on scaling considerations from measurements at the LHC. Such a large amount of stored energy in the beam tails poses serious challenges on its control and safe disposal. In particular, orbit jitters can cause significant losses on primary collimators, which can lead to accidental beam dumps, magnet quenches, or even permanent damage to collimators and other accelerator elements. Thus, active control of the diffusion speed of halo particles is necessary and the use of hollow electron lenses (HELs) represents the most promising approach to handle overpopulated tails at the HL-LHC. HEL is a very powerful and advanced tool that can be used for controlled depletion of beam tails, thus enhancing the performance of beam halo collimation. For these reasons, HELs have been recently included in the HL-LHC baseline. In this paper, we present detailed beam dynamics calculations performed with the goal of defining HEL specifications and operational scenarios for HL-LHC. The prospects for effective halo control in HL-LHC are presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1985

Time-dependent nonlinear cosmic ray shocks confirming abstract

Numerical studies of time dependent cosmic ray shock structures in planar geometry are interesting because analytical time-independent solutions are available which include the non-linear reactions on the plasma flow. A feature of these time asymptotic solutions is that for higher Mach numbers (M approximately 5) and for a low cosmic ray upstream pressure the solution is not uniquely determined by the usual conservation laws of mass, momentum and energy. These numerical solutions clearly indicate that much work needs to be done before we understand shock acceleration as a time dependent process. The slowness of the process is possibly due to the fact that there is a diffusive flux into the downstream region in addition to the usual advective losses. Analytic investigations of thie phenomenon are required.

Dorfi, E. A.↗
NASA NTRS1988

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗
NASA NTRS1987

Weak double layers in the auroral ionosphere

Previous work on the evolution of weak double layers in a hydrogen plasma was extended to include H(+) and O(+) with relative drift. The relative drift between hydrogen and oxygen ions due to a quasi-static parallel electric field gives rise to a strong linear fluid instability which dominates the ion-acoustic mode at the bottom of the auroral acceleration region. This ion-ion instability can modify ion distributions at lower altitudes and the subsequent nonlinear evolution of weak double layers at higher altitudes in the ion-acoustic regime. Ion hole formation can occur for smaller relative electron-ion drifts than seen in previous simulations, due to the hydrogen-oxygen two-stream instability. This results in local modification of the ion distributions in phase space, and a partial filling of the valley between the hydrogen and oxygen peaks, which would be expected at higher altitudes on auroral field lines. The observed velocity diffusion does not necessarily preclude ion hole and double layer formation in hydrogen in the ion-acoustic regime. These simulation results are consistent with the experimentally measured persistence of separate hydrogen and oxygen peaks, and the observation of weak double layers above an altitude of 3000 km on auroral field lines.

Hudson, M. K.↗
DOE OSTI2021

Data-driven acceleration of thermal radiation transfer calculations with the dynamic mode decomposition and a sequential singular value decomposition

In this work, we present a method for accelerating discrete ordinates radiative transfer calculations for radiative transfer. Our method works with nonlinear positivity fixes, in contrast to most acceleration schemes. The method is based on the dynamic mode decomposition (DMD) and using a sequence of rank-one updates to compute the singular value decomposition needed for DMD. Using a sequential method allows us to automatically determine the number of solution vectors to include in the DMD acceleration. We present results for slab geometry discrete ordinates calculations with the standard temperature linearization. Compared with positive source iteration, our results demonstrate that our acceleration method reduces the number of transport sweeps required to solve the problem by a factor of about 3 on a standard diffusive Marshak wave problem, a factor of several thousand on a cooling problem where the effective scattering ratio approaches unity, and a factor of 20 improvement in a realistic, multimaterial radiating shock problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Design of double- and multi-bend achromat lattices with large dynamic aperture and approximate invariants

A numerical method to design nonlinear double- and multi-bend achromat (DBA and MBA) lattices with approximate invariants of motion is investigated. The search for such nonlinear lattices is motivated by Fermilab’s Integrable Optics Test Accelerator (IOTA), whose design is based on an integrable Hamiltonian system with two invariants of motion. While it may not be possible to design an achromatic lattice for a dedicated synchrotron light source storage ring with one or more exact invariants of motion, it is possible to tune the sextupoles and octupoles in existing DBA and MBA lattices to produce approximate invariants. In our procedure, the lattice is tuned while minimizing the turn-by-turn fluctuations of the Courant-Snyder actions J x and J y at several distinct amplitudes, while simultaneously minimizing diffusion of the on-energy betatron tunes. The resulting lattices share some important features with integrable ones, such as a large dynamic aperture, trajectories confined to invariant tori, robustness to resonances and errors, and a large amplitude-dependent tune-spread. Compared to the nominal NSLS-II lattice, the single- and multi-bunch instability thresholds are increased and the bunch-by-bunch feedback gain can be reduced.

43 PARTICLE ACCELERATORS↗
NASA NTRS2006

Dynamic Theory of Relativistic Electrons Stochastic Heating by Whistler Mode Waves with Application to the Earth Magnetosphere

In the Hamiltonian approach an electron motion in a coherent packet of the whistler mode waves propagating along the direction of an ambient magnetic field is studied. The physical processes by which these particles are accelerated to high energy are established. Equations governing a particle motion by group symmetries of the problem were transformed in to a closed pair of nonlinear difference equations. The solutions of these equations have shown there exists the energetic threshold below that the electron motion is regular, and when the initial energy is above the threshold an electron moves stochastically. It is proved that the upper boundary of particle stochastic heating is conditioned by intrinsic property of the particle chaotic motion. Particle energy spectra and pitch angle electron scattering are described by the Fokker-Planck-Kolmogorov equations. It is shown that significant pitch angle diffusion occurs for the Earth radiation belt electrons with energies from a few keV up to a few MeV.

Khazanov, G. V.↗
NASA NTRS2018

Precipitating Electron Energy Flux and Characteristic Energies in Jupiter's Main Auroral Region as Measured by Juno/JEDI

The relationship between electron energy flux and the characteristic energy of electron distributions in the main auroral loss cone bridges the gap between predictions made by theory and measurements just recently available from Juno. For decades such relationships have been inferred from remote sensing observations of the Jovian aurora, primarily from the Hubble Space Telescope, and also more recently from Hisaki. However, to infer these quantities, remote sensing techniques had to assume properties of the Jovian atmospheric structure - leading to uncertainties in their profile. Juno's arrival and subsequent auroral passes have allowed us to obtain these relationships unambiguously for the first time, when the spacecraft passes through the auroral acceleration region. Using Juno /Jupiter Energetic particle Detector Instrument (JEDI), an energetic particle instrument, we present these relationships for the 30-kiloelectronvolts to 1-megaelectronvolts electron population. Observations presented here show that the electron energy flux in the loss cone is a nonlinear function of the characteristic or mean electron energy and supports both the predictions from Knight (1973, https://doi.org/10.1016/0032-0633(73)90093-7) and magnetohydrodynamic turbulence acceleration theories (e.g., Saur et al., 2003, https://doi.org/10.1029/2002GL015761). Finally, we compare the in situ analyses of Juno with remote Hisaki observations and use them to help constrain Jupiter's atmospheric profile. We find a possible solution that provides the best agreement between these data sets is an atmospheric profile that more efficiently transports the hydrocarbons to higher altitudes. If this is correct, it supports the previously published idea (e.g., Parkinson et al., 2006, https://doi.org/10.1029/2005JE002539) that precipitating electrons increase the hydrocarbon eddy diffusion coefficients in the auroral regions.

Auroral↗
DOE OSTI2020

Understanding the eco-geomorphologic feedback of coastal marsh under sea level rise: vegetation dynamic representations, processes interaction, and parametric sensitivity

A growing number of coastal eco-geomorphologic modeling studies have been conducted to understand coastal marsh evolution under sea-level rise (SLR). Although these models quantify marsh topographic change as a function of sedimentation and erosion, their representations of vegetation dynamics that control organic sedimentation differ. How vegetation dynamic schemes contribute to simulation outcomes is not well quantified. Additionally, the sensitivity of modeling outcomes to parameter selection in the available formulations has not been rigorously tested to date, especially under the influence of an accelerating SLR. In this paper, we used a coastal eco-geomorphologic model with different vegetation dynamic schemes to investigate the eco-geomorphologic feedbacks of coastal marshes and parametric sensitivity under SLR scenarios. We found that marsh platform relief increased with sea level rise rate. The simulations with different vegetation schemes exhibited different spatial-temporal variations in elevation and biomass. The nonlinear Spartina scheme presented the most resilient prediction with generally the highest marsh accretion and vegetation biomass, and the least elevation relief under SLR. But the linear Spartina scheme predicts the lowest unvegetated-vegetated ratio. We also found that vegetation-related parameters and sediment diffusivity, which were not well measured or discussed in previous studies, were identified as some of the most critical parameters. Additionally, the model sensitivity to vegetation-related parameters increased with SLR rates. The identified most sensitive parameters may inform how to appropriately choose modeling representations of key processes and parameters for different coastal marsh landscapes under SLR, and demonstrate the importance of future field measurements of these key parameters.

58 GEOSCIENCES↗
NASA NTRS1999

Soot Aerosol Properties in Laminar Soot-Emitting Microgravity Nonpremixed Flames

The spatial distributions and morphological properties of the soot aerosol are examined experimentally in a series of 0-g laminar gas-jet nonpremixed flames. The methodology deploys round jet diffusion flames of nitrogen-diluted acetylene fuel burning in quiescent air at atmospheric pressure. Full-field laser-light extinction is utilized to determine transient soot spatial distributions within the flames. Thermophoretic sampling is employed in conjunction with transmission electron microscopy to define soot microstructure within the soot-emitting 0-g flames. The microgravity tests indicate that the 0-g flames attain a quasi-steady state roughly 0.7 s after ignition, and sustain their annular structure even beyond their luminous flame tip. The measured peak soot volume fractions show a complex dependence on burner exit conditions, and decrease in a nonlinear fashion with decreasing characteristic flow residence times. Fuel preheat by approximately 140 K appears to accelerate the formation of soot near the flame axis via enhanced fuel pyrolysis rates. The increased soot presence caused by the elevated fuel injection temperatures triggers higher flame radiative losses, which may account for the premature suppression of soot growth observed along the annular region of preheated-fuel flames. Electron micrographs of soot aggregates collected in 0-g reveal the presence of soot precursor particles near the symmetry axis at midflame height, The observations also verify that soot primary particle sizes are nearly uniform among aggregates present at the same flame location, but vary considerably with radius at a fixed distance from the burner. The maximum primary size in 0-g is found to be by 40% larger than in 1-g, under the same burner exit conditions. Estimates of the number concentration of primary particles and surface area of soot particulate phase per unit volume of the combustion gases are also made for selected in-flame locations.

Konsur, Bogdan↗
DOE OSTI2023

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Three-dimensional spontaneous magnetic reconnection

Magnetic reconnection is best known from observations of the Sun where it causes solar flares. Observations estimate the reconnection rate as a small, but non-negligible fraction of the Alfvén speed, so-called fast reconnection. Until recently, the prevailing pictures of reconnection were either of resistivity or plasma microscopic effects, which was contradictory to the observed rates. Alternative pictures were either of reconnection due to the stochasticity of magnetic field lines in turbulence or the tearing instability of the thin current sheet. In this paper we simulate long-term three-dimensional nonlinear evolution of a thin, planar current sheet subject to a fast oblique tearing instability using direct numerical simulations of resistive-viscous magnetohydrodynamics. The late-time evolution resembles generic turbulence with a −5/3 power spectrum and scale-dependent anisotropy, so we conclude that the tearing-driven reconnection becomes turbulent reconnection. The turbulence is local in scale, so microscopic diffusivity should not affect large-scale quantities. This is confirmed by convergence of the reconnection rate toward ∼0.015v{sub A} with increasing Lundquist number. In this spontaneous reconnection, with mean field and without driving, the dissipation rate per unit area also converges to ∼0.006ρv{sub A}{sup 3}, and the dimensionless constants 0.015 and 0.006 are governed only by self-driven nonlinear dynamics of the sheared magnetic field. Remarkably, this also means that a thin current sheet has a universal fluid resistance depending only on its length to width ratio and to v{sub A}/c.

79 ASTRONOMY AND ASTROPHYSICS↗