Magnetoplasma diffusion at F2 region altitudes
Magnetoplasma diffusion equation derivation for F-2 layer of ionosphere
F-2 LAYER↗
SEARCH · Engineering Papers
Results for “Diffusion equation”
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.
Energization in regions of CIRs unconnected to shocks are probably not the result of cross-field transport
Corotating energetic ion populations (CEIPs) associated with the forward and reverse shocks of corotating interaction regions (CIRs) are observed in CIRs at places where models say are magnetically unconnected to either shock. Such disconnections between CEIPs and shocks are common and have been documented with data from Pioneers 10 and 11 and confirmed with data from Ulysses. They pose a problem for models that account for these CEIPs in terms of ion energization at the shocks followed by ion propagation along field lines. Two possible resolutions to this problem have been suggested: diffusion of the ions across field lines and extension of the ion energization process to regions beyond the shock waves. Here we quantitatively examine the first of these possibilities. We give the Green's function solution to the convection-diffusion equation applied to idealized CIR geometry, with a source at the reverse shock -- the main producer of CEIPs. Two kinds of diffusion are considered: resonant diffusion and stochastic field line diffusion. We find that for resonant diffusion the computed ratio is many orders of magnitudes below the observed ratio. For stochastic field line diffusion, the computed ratio approximately equals the observed ratio if a diffusion coefficient appropriate to the free solar wind is used. It is several orders of magnitude below the observed ratio, however, if a diffusion coefficient appropriate to CIRs is used. We conclude that cross-field diffusion probably does not account for the presence of energetic ions in regions of CIRs that are magnetically unconnected to its shock waves. We suggest that the alternative possibility -- the energetic ions in regions magnetically unconnected to shocks result from an acceleration process that is independent of shocks -- be pursued to the point where quantitative tests can be performed.
Intriligator, Devrie S.↗
Random element method for numerical modeling of diffusional processes
The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.
Ghoniem, A. F.↗
Multigrid Techniques for Highly Indefinite Equations
A multigrid method for the solution of finite difference approximations of elliptic PDE's is introduced. A parallelizable version of it, suitable for two and multi level analysis, is also defined, and serves as a theoretical tool for deriving a suitable implementation for the main version. For indefinite Helmholtz equations, this analysis provides a suitable mesh size for the coarsest grid used. Numerical experiments show that the method is applicable to diffusion equations with discontinuous coefficients and highly indefinite Helmholtz equations.
Shapira, Yair↗
Cellular automaton formulation of passive scalar dynamics
Cellular automata modeling of the advection of a passive scalar in a two-dimensional flow is examined in the context of discrete lattice kinetic theory. It is shown that if the passive scalar is represented by tagging or 'coloring' automation particles a passive advection-diffusion equation emerges without use of perturbation expansions. For the specific case of the hydrodynamic lattice gas model of Frisch et al. (1986), the diffusion coefficient is calculated by perturbation.
Chen, Hudong↗
Radial diffusion in the Uranian radiatian belts - Inferences from satellite absorption loss models
Low-energy charged particle (LECP) phase space density profiles available from the Voyager/1986 Uranus encounter are analyzed, using solutions of the time-averaged radial diffusion equation for charged particle transport in a dipolar planetary magnetic field. Profiles for lower-energy protons and electrons are first analyzed to infer radial diffusion rate as a function of L, assuming that satellite absorption is the dominant loss process and local sources for these particles are negligible. Satellite macrosignatures present in the experimentally derived profiles are approximately reproduced in several cases, lending credence to the loss model and indicating that magnetospheric distributed losses are not as rapid as satellite absorption near the minimum satellite L shells for the particles. Diffusion rates and L dependences are found to be similar to those previously inferred in the inner Jovian magnetosphere (Thomsen et al., 1977) and for the inner Saturnian magnetosphere (Hood, 1985). Profiles for higher energy electrons and protons are also analyzed using solutions that allow for the existence of significant particle sources as well as sinks. Possible implications for radial diffusion mechanisms in the Uranian radiation belts are discussed.
Hood, L. L.↗
Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions
A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.
plasma↗
An Exponential Finite Difference Technique for Solving Partial Differential Equations
An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.
Handschuh, Robert F.↗
Seismic energy transmission in an intensively scattering environment
In order to account for some special features of lunar seismograms, namely, the gradual build-up of the signal, the extremely prolonged tail, and the lack of apparent coherence among three orthogonal components of ground motion, a statistical approach is proposed for describing transmission of seismic energy through a medium in which strong scattering takes place. A seismic diffusion theory is presented for a medium with randomly distributed scatterers of a given size distribution. A solution of the resulting diffusion equation for an impulsive energy source gives a curve which fairly closely reproduces the envelope of typical lunar impact seismograms. Since the model is based on constant diffusivity, long-range transmission will require a combination of diffusion and wave propagation treatments for accurate description.
Nakamura, Y.↗
The expected cosmic ray density and stream distributions at the heliolatitudinal asymmetry of solar wind
The results of the spatial distribution of cosmic ray density, gradients, and anisotropy obtained on the basis of the numerical solution of the anisotropic diffusion equation with an account of solar wind velocity change depending on the latitudinal angle theta of the form U=u sub oe sup alpha theta and the diffusion coefficient depending on the spatial coordinates and the particle rigidity are presented. It is shown that the increase of the solar wind velocity and the diffusion coefficient with heliolatitude leads to gradient distributions that are in accord with experimental data observed in space. The results of the energetic spectrum of 11 and 22-year cosmic ray variations obtained with an account of direction of the general magnetic field of the Sun are presented are given.
Alania, M. V.↗
Characteristics-based methods applied to infinite Prandtl number thermal convection in the hard turbulent regime
Characteristics-based methods for the advection-diffusion equation are presented and directly applied to study thermal convection with extremely large Rayleigh number (Ra). It is shown that the operator-splitting method for advection-diffusion problems is very accurate for determining the advected field at extremely high Peclet number (Pe). The technique presented is considered to have great potential for solving advection-dominated problems, while the Langrangian method is more accurate for lower Pe. It is noted that the accuracy of these characteristics-based methods strongly depends on the quality of interpolation. The computational time for the operator-splitting method grows with the number of time steps employed. The Langrangian method was used for simulations of convection at very high Ra, up to 3 x 10 to the 9th, and time-dependent, thermal convection solutions were obtained for infinite Prandtl number.
Malevsky, A. V.↗
Robust and Accurate Shock Capturing Method for High-Order Discontinuous Galerkin Methods
A simple yet robust and accurate approach for capturing shock waves using a high-order discontinuous Galerkin (DG) method is presented. The method uses the physical viscous terms of the Navier-Stokes equations as suggested by others; however, the proposed formulation of the numerical viscosity is continuous and compact by construction, and does not require the solution of an auxiliary diffusion equation. This work also presents two analyses that guided the formulation of the numerical viscosity and certain aspects of the DG implementation. A local eigenvalue analysis of the DG discretization applied to a shock containing element is used to evaluate the robustness of several Riemann flux functions, and to evaluate algorithm choices that exist within the underlying DG discretization. A second analysis examines exact solutions to the DG discretization in a shock containing element, and identifies a "model" instability that will inevitably arise when solving the Euler equations using the DG method. This analysis identifies the minimum viscosity required for stability. The shock capturing method is demonstrated for high-speed flow over an inviscid cylinder and for an unsteady disturbance in a hypersonic boundary layer. Numerical tests are presented that evaluate several aspects of the shock detection terms. The sensitivity of the results to model parameters is examined with grid and order refinement studies.
Atkins, Harold L.↗
Higher-order numerical methods derived from three-point polynomial interpolation
Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and Hermitian finite-difference discretization. The equations generally apply for both uniform and variable meshes. Hybrid schemes resulting from different polynomial approximations for first and second derivatives lead to the nonuniform mesh extension of the so-called compact or Pade difference techniques. A variety of fourth-order methods are described and this concept is extended to sixth-order. Solutions with these procedures are presented for the similar and non-similar boundary layer equations with and without mass transfer, the Burgers equation, and the incompressible viscous flow in a driven cavity. Finally, the interpolation procedure is used to derive higher-order temporal integration schemes and results are shown for the diffusion equation.
Rubin, S. G.↗
New approach to cosmic-ray diffusion theory.
We have investigated a new approach to deriving a diffusion equation for charged particles in a static, random magnetic field. Our method incorporates essential effects of the magnetic fluctuations in the lowest order particle orbits. Significant corrections to the usual quasilinear diffusion coefficient for cosmic rays with pitch angles near 90 deg are a consequence. Monte Carlo results bear out the validity of our theory.
Jones, F. C.↗
Numerical study of solar flare particle propagation in the heliosphere
Numerical solutions are presented for the propagation of solar cosmic rays in interplanetary space, including the effects of pitch-angle scattering and adiabatic focusing. The intensity-time profiles can be well fitted by a simple radial spatial diffusion equation. For low-rigidity particles the radial mean free path so obtained is significantly larger than the mean free path calculated from the scattering coefficient due to the inapplicability of the diffusive approximation early in the event. The well-known discrepancy between the scattering mean free path and the theoretical predictions may be resolved by these calculations.
Gombosi, T. I.↗
New developments in the method of space-time conservation element and solution element: Applications to the Euler and Navier-Stokes equations
A new numerical framework for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods--i.e., finite difference, finite volume, finite element, and spectral methods. It is conceptually simple and designed to avoid several key limitations to the above traditional methods. An explicit model scheme for solving a simple 1-D unsteady convection-diffusion equation is constructed and used to illuminate major differences between the current method and those mentioned above. Unexpectedly, its amplification factors for the pure convection and pure diffusion cases are identical to those of the Leapfrog and the DuFort-Frankel schemes, respectively. Also, this explicit scheme and its Navier-Stokes extension have the unusual property that their stabilities are limited only by the CFL condition. Moreover, despite the fact that it does not use any flux-limiter or slope-limiter, the Navier-Stokes solver is capable of generating highly accurate shock tube solutions with shock discontinuities being resolved within one mesh interval. An accurate Euler solver also is constructed through another extension. It has many unusual properties, e.g., numerical diffusion at all mesh points can be controlled by a set of local parameters.
Chang, Sin-Chung↗
Fermi-Compton scattering due to magnetopause surface fluctuations in Jupiter's magnetospheric cavity
The effects of boundary surface fluctuations on a spectrum of electromagnetic radiation trapped in a high Q (quality) cavity are considered. Undulating walls introduce small frequency shifts at reflection to the radiation, and it is argued that the process is entirely analogous to both Fermi (particle) acceleration and inverse Compton scattering. A Fokker-Planck formalism is pursued; it yields a diffusion equation in frequency for which the Green's function and steady-state solutions are found. Applying this analysis to the Jovian continuum radiation discovered by Voyager spacecraft, it is suggested that characteristic diffusion times are greater than 1 year, and that in order to account for the steep frequency spectra observed, an unidentified loss mechanism must operate in the cavity with a decay time constant approximately equal to the characteristic diffusion time divided by 28. A radiator-reactor model of the cavity is investigated to provide an estimate for the intrinsic luminosity of the low frequency (approximately 100 Hz) continuum source whose power is approximately 7 x 10 to the 6th W.
Barbosa, D. D.↗
Convergence of infinite dimensional sampled LQR problems - Theory and numerical results
A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.
Rosen, I. G.↗