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 181 records · Page 10

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗

Stochastic Fermi acceleration in solar flares

Proton spectra, valid from non- to ultra-relativistic energies, resulting from stochastic Fermi acceleration in solar flares are calculated. These spectra were obtained by numerically solving the Fokker-Planck equation, in which the escape of the particles from the acceleration region is characterized by an energy-independent escape time. In addition to equilibrium spectra, time-dependent energy spectra showing the approach to equilibrium are also presented. These numerical equilibrium spectra are compared with previous results which were obtained either by Monte Carlo simulations or approximate analytical treatments. There are no analytic solutions valid in the transrelativistic regime, which is very important for the production of pions and neutrons in solar flares. The acceleration efficiency is related to physical parameters, in particular the energy density in either magnetosonic or Alfven waves, and a lower limit is placed on either of these energy densities from acceleration times implied by gamma-ray observations. Also discussed is the physical interpretation of the escape time.

Miller, James A.↗

Non-Maxwellian electron distributions in models of the solar atmosphere

In this paper, the solar models of Fontenla et al. (1990) are extended to coronal temperatures using optically thin radiative losses and a semicircular loop geometry. The effects of a non-Maxwellian high-energy tail of the electron distribution function is tested by performing a Fokker-Planck calculation of the tail. The results show that the downward conductive heat flux is very clsoe to that given by Braginskii (1965) for small departures from a Maxwellian distribution. The effects of the high-energy tail of the electron distribution on element ionization are computed and found to be negligible for all ionization states of O, Ne, and Si. It is concluded that, for quiet sun loop models in stationary energy balance, the ionization state of all elements is not significantly affected by the non-Maxwellian tail of the electron velocity distribution, and that the heat flux can be computed by considering only small first-order departures from the Maxwellian distribution.

Macneice, P.↗

Performance evaluation of a Doppler radar system for wind shear detection

Nonlinear stochastic differential equations are used to model wind shear, and extended Kalman filters are used to generate state estimates from measurements received from a Doppler radar onboard an aircraft. Likelihood-ratio tests are then used to detect the presence of wind shear. The performance of the system is evaluated by deriving theoretical expressions for the false alarm and miss error probabiilties. The approach uses a Fokker-Planck equation. The overall methodology is general and should be of interest in other applications.

Khalaf, Camille S.↗

Focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection

Evidence is presented for focused transport of energetic particles along magnetic field lines draped around a coronal mass ejection. This evidence was obtained with the University of Maryland/Max-Planck-Institute experiment on the ISEE-3 spacecraft during the decay phase of the June 6, 1979, solar particle event. During the early portion of the decay phase of this event, interplanetary magnetic field lines were apparently draped around a coronal mass ejection, leading to a small focusing length on the western flank where ISEE 3 was located. A period of very slow decrease of particle intensity was observed, along with large sunward anisotropy in the solar wind frame, which is inconsistent with predictions of the standard Fokker-Planck equation models for diffusive transport. It was found possible to fit the observations, assuming that focused transport dominates and that the particle pitch angle scattering is isotropic.

Tan, L. C.↗

The Sobolev approximation for line formation with partial frequency redistribution

Attention is given to the formation of a spectral line in a uniformly expanding infinite medium in the Sobolev approximation, with emphasis on the various mechanisms for frequency redistribution. Numerical and analytic solutions of the transfer equation are presented of a number of redistribution functions and their approximations, including type I and type II partial redistribution, coherent scattering and complete redistribution, and the Fokker-Planck and uncorrelated approximation to the R sub II function. The solutions for the mean intensity are shown to depend very much on the type of redistribution mechanism, while for the frequency-weighted mean intensity, which enters the rate equations, this dependence is weak. It is inferred that use of Sobolev escape probabilities based on complete redistribution can be an adequate approximation for many calculations for which only the radiative excitation rates are needed.

Hummer, D. G.↗

Evolving, dynamical models for collapsed-core globular clusters - M15 and NGC 6624

The stellar populations of the collapsed-core globular clusters M15 and NGC 6624 are investigated by fitting observed surface-brightness and projected velocity-dispersion profiles. The present evolving cluster models were generated by the direct Fokker-Planck method and incorporate realistic stellar mass spectra and energy input from binaries formed by three-body interactions. An evolved power-law mass function with nonluminous remnants of maximum mass 1.0-1.4 solar mass is adopted. M15 is found to be best fitted by postcollapse evolving models with a mass-function slope of x = 0.9 (where x = 1.35 corresponds to the Salpeter mass function) and a maximum remnant mass of 1.3 solar mass. Postcollapse core oscillations, driven by energy input from hard binaries, can produce sufficiently rapid core expansion to explain the resolution of the core of M15 by HST observations reported by Lauer et al. (1991).

Grabhorn, Robert P.↗

Kinetic description of electron beams in the solar chromosphere

We formulate the relativistic Fokker-Plank equation for a beam of accelerated electrons interacting with a partially ionized plasma. In our derivation we conserved those terms contributing to velocity diffusion and found that this effect cannot be neglected a priori. We compute the terms accounting for elastic and inelastic collisions with neutral hydrogen and helium. Collisions with neutral hydrogen are found to be dominant throughout the chromosphere, except at the uppermost layers close to the transition region. As an application, we compute the loss of energy and momentum for a power-law beam impinging on the solar chromosphere, for a particular case in which the Fokker-Planck equation can be integrated analytically. We find that most of the beam energy is deposited in a relatively thin region of the chromosphere, a result which is largely insensitive to the theoretical method employed to compute the energy deposition rate.

Gomez, Daniel O.↗

Stochastic inflation lattice simulations - Ultra-large scale structure of the universe

Non-Gaussian fluctuations for structure formation may arise in inflation from the nonlinear interaction of long wavelength gravitational and scalar fields. Long wavelength fields have spatial gradients, a (exp -1), small compared to the Hubble radius, and they are described in terms of classical random fields that are fed by short wavelength quantum noise. Lattice Langevin calculations are given for a toy model with a scalar field interacting with an exponential potential where one can obtain exact analytic solutions of the Fokker-Planck equation. For single scalar field models that are consistent with current microwave background fluctuations, the fluctuations are Gaussian. However, for scales much larger than our observable Universe, one expects large metric fluctuations that are non-Gaussian. This example illuminates non-Gaussian models involving multiple scalar fields which are consistent with current microwave background limits.

Salopek, D. S.↗

U-band CCD surface brightness profiles for candidate collapsed-core globular clusters

We have determined U-band CCD surface brightness profiles for 15 of the 21 globular clusters that have been characterized as having 'collapsed' cores by Djorgovski and King (1986). We have fitted these profiles with two types of seeing-convolved models: single-mass King models and power laws. We find that in many cases the characterization of a profile as being either 'King' or 'cusp' is somewhat subtle. Of the clusters in our sample, 10 appear to have central power-law cusps with either unresolved or marginally resolved cores, three appear to have resolved cores with surrounding power-law structure that is unlike a King model, and two are approximately described by King model fits. The observed power-law cusps are consistent with the predictions of Fokker-Planck models for core-collapse.

Lugger, Phyllis M.↗

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗

Nonlinear damping model for flexible structures

The study of nonlinear damping problem of flexible structures is addressed. Both passive and active damping, both finite dimensional and infinite dimensional models are studied. In the first part, the spectral density and the correlation function of a single DOF nonlinear damping model is investigated. A formula for the spectral density is established with O(Gamma(sub 2)) accuracy based upon Fokker-Planck technique and perturbation. The spectral density depends upon certain first order statistics which could be obtained if the stationary density is known. A method is proposed to find the approximate stationary density explicitly. In the second part, the spectral density of a multi-DOF nonlinear damping model is investigated. In the third part, energy type nonlinear damping model in an infinite dimensional setting is studied.

Zang, Weijian↗

Stochastic modeling of turbulent reacting flows

Direct numerical simulations of a single-step irreversible chemical reaction with non-premixed reactants in forced isotropic turbulence at R(sub lambda) = 63, Da = 4.0, and Sc = 0.7 were made using 128 Fourier modes to obtain joint probability density functions (pdfs) and other statistical information to parameterize and test a Fokker-Planck turbulent mixing model. Preliminary results indicate that the modeled gradient stretching term for an inert scalar is independent of the initial conditions of the scalar field. The conditional pdf of scalar gradient magnitudes is found to be a function of the scalar until the reaction is largely completed. Alignment of concentration gradients with local strain rate and other features of the flow were also investigated.

Fox, R. O.↗

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.]↗

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↗