Stochastic magnetic pumping.
Collisionless plasma heating by random magnetic fields using simple betatron model
SEARCH · Engineering Papers
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.
Collisionless plasma heating by random magnetic fields using simple betatron model
We present a new analysis of the fundamental physics of charged-particle motion in a turbulent magnetic field using a numerical simulation. The magnetic field fluctuations are taken to be static and to have a power spectrum which is Kolmogorov. The charged particles are treated as test particles. It is shown that when the field turbulence is independent of one coordinate (i.e., k lies in a plane), the motion of these particles across the magnetic field is essentially zero, as required by theory. Consequently, the only motion across the average magnetic field direction that is allowed is that due to field-line random walk. On the other hand, when a fully three-dimensional realization of the turbulence is considered, the particles readily cross the field. Transport coefficients both along and across the ambient magnetic field are computed. This scheme provides a direct computation of the Fokker-Planck coefficients based on the motions of individual particles, and allows for comparison with analytic theory.
The basic theoretical work required to understand sound transmission into an enclosed space (that is, one closed by the transmitting structure) is developed for random pressure fields and for harmonic (tonal) excitation. The analysis is used to predict the noise reducton of unpressurized unstiffened cylinder, and also the interior response of the cylinder given a tonal (plane wave) excitation. Predictions and measurements are compared and the transmission is analyzed. In addition, results for tonal (harmonic) mechanical excitation are considered.
The performance of coherent Doppler lidar is determined by the statistics of the coherent Doppler signal. The derivation and calculation of the covariance of the Doppler lidar signal is presented for random atmospheric wind fields with wind shear. The random component is described by a Kolmogorov turbulence spectrum. The signal parameters are clarified for a general coherent Doppler lidar system. There are two distinct physical regimes: one where the transmitted pulse determines the signal statistics and the other where the wind field dominates the signal statistics. The Doppler shift of the signal is identified in terms of the wind field and system parameters.
Radar backscattering coefficient measurements by spaceborne scatterometers are presently simulated for the case of nonuniform wind fields, by means of a detailed numerical integration of the radar equation. The winds thus estimated are then compared with a nominal field which is defined as the average wind vector over the wind cell. The simulation results obtained for the NASA scatterometer are presented for cases of random wind fields whose spectra are consistent with the Seasat scatterometer sea surface wind spectrum. When the nonuniformity is small, system noise dominates the wind error; wind error degradation is therefore small for both perfect and imperfect coregistration cases. When it is relatively large, however, the wind error degradation persistently increases for both perfect and imperfect coregistrations.
Most crystalline materials present a highly heterogeneous response at the microscale, which can be affected by both internal factors (such as microstructural parameters) and external factors (such as loading). Relating microscale inhomogeneities to the macroscale response of a material requires the use of homogenisation techniques, usually based on the concept of a representative volume element (RVE)—the smallest volume of material that represents the global average response. Here, in this work, we present a new and robust experimental method of measuring the size of a strain-based RVE from high-resolution grain-scale strain fields obtained using digital image correlation (DIC). The proposed method is based on the statistical (stereological) nature of the RVE, which has been widely adopted in numerical studies, and involves dividing a strain field into randomly selected regions of varying sizes and statistically analysing the distributions of average strains within them. To validate the new method, we generate a large number of synthetic strain fields from a fractional Gaussian noise algorithm. The proposed stereological method is shown to be capable of producing reliable RVE measurements from a very large range of possible microscale strain fields while at the same time being robust in that it can produce RVE measurement results even in cases where other existing methods may be unable to do so. The proposed method has a low field-of-view requirement, only needing a field-of-view about 1.2 times as large as the RVE to produce reliable measurements. In addition, the stereological method offers significant flexibility since its statistical nature allows for control over how strict the RVE measurement should be in each case.
Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.
The methods used at the Goddard Laboratory for Atmospheres to retrieve sea surface temperatures from HIRS2/MSU data for the four months of the NASA sea surface temperature intercomparison workshop are described. Results are shown comparing anomaly fields produced using data from ships, AVHRR, HIRS2/MSU, SMMR AND VAS for the last three of these months. Fields from AVHRR and HIRS2 show the highest accuracy compared to ship fields. Errors in the HIRS2 fields appear more random while AVHRR data shows large area, spatially coherent errors. The random errors in the HIRS2 fields can be further reduced by performing the retrievals at a higher spatial resolution.
The methods used at the Goddard Laboratory for Atmospheres (GLA) to retrieve sea surface temperatures from HIRS2/MSU data for the last 2 months of the NASA sea surface temperature intercomparison workshop are described. Results are shown that compare anomaly fields produced by using data from ships, AVHRR, and HIRS2/MSU for these months. Monthly mean fields derived from AVHRR and HIRS2 data show roughly comparable statistical accuracy compared to ship fields. Errors in the HIRS2 fields appear more random, while AVHRR data have larger regionally dependent biases and show large area spatially coherent errors. The random errors in the HIRS2 fields should be further reduced by performing the retrievals at a higher spatial resolution.
Based on the complete Stokes wave model with the bias term and using a simple mapping approach and an iteration solution method, we established a formula for the joint probability density function of the surface slope elevation of a nonlinear random wave field. The formula requires three parameters to define the whole density function: the rms surface elevation and slope values and the significant slope. This model represents the dynamics of the wave in a more direct way than the Gram-Charlier approximation. Based on this new statistical model and laboratory experiments, formula and numerical values of EM bias and dynamics bias are derived. The results indicate that various biases should be considered seriously if accuracy of the altimeter measurement is required in centimeter range.
A popular prescription for the galaxy formation, i.e., identifying peaks in the linear density field as sites of galaxy formation, is numerically justified under the Cold Dark Matter cosmogony. In particular, the peak-background scheme, originally developed to simplify analytic calculations of peak properties in a Gaussian random density field, is proven to be an excellent prescription for assigning peak tracers in low-resolution N-body simulations through many statistical comparisons. It is found that statistical properties of peak tracers allocated by the peak-background scheme closely resemble those of true peaks when the background smoothing length is about three times the galaxy scale and the correspondence improves as they evolve gravitationally. It is also shown that these peak particles are indeed strongly associated with nonlinear objects that are presently collapsed.
A micromechanical method is employed for the prediction of unidirectional composites in which the fiber orientation can possess various statistical misalignment distributions. The method relies on the probability-weighted averaging of the appropriate concentration tensor, which is established by the micromechanical procedure. This approach provides access to the local field quantities throughout the constituents, from which initiation of damage in the composite can be predicted. In contrast, a typical macromechanical procedure can determine the effective composite elastic properties in the presence of statistical fiber misalignment, but cannot provide the local fields. Fully random fiber distribution is presented as a special case using the proposed micromechanical method. Results are given that illustrate the effects of various amounts of fiber misalignment in terms of the standard deviations of in-plane and out-of-plane misalignment angles, where normal distributions have been employed. Damage initiation envelopes, local fields, effective moduli, and strengths are predicted for polymer and ceramic matrix composites with given normal distributions of misalignment angles, as well as fully random fiber orientation.
The Earth's atmosphere is continually being bombarded by primary cosmic ray particles which are generally believed to be high-energy nuclei. The fact that the majority of cosmic ray primaries are charged particles and that space is permeated with random magnetic fields, means that the particles do not travel in straight lines. The arrival time distribution of EAS may also transfer some information about the primary particles. Actually, if the particles come to our Earth in a completely random process, the arrival time distribution of pairs of successive particles should fit an exponential law. The work reported here was arried out at Sydney University from May 1982 to January 1983. All the data are used to plot the arrival-time distribution of the events, that is, the distribution of time-separation between consecutive events on a 1 minute bin size. During this period more than 2300 showers were recorded. The results are discussed and compared with that of some other experiments.
An avctive flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A flat plate towed in a water channel is used as a test bed. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time during its 'on' period. When the jet is turned 'off', all of the turbulent fluid is released instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, the cyclic jet produces periodic structures that are qualitatively similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.
An active flow control device to generate large-scale, periodic structures in a turbulent shear flow is developed. Together with adaptive optics, the device may be used on airborne laser platforms to reduce or eliminate optical distortion caused by the turbulence in the aircraft's boundary layer. A cyclic jet issuing from a spanwise slot is used to collect the turbulent boundary layer for a finite time and then release all of the flow instantaneously in one large eddy that convects downstream. Flow visualization and hot-film probe measurements are used together with pattern recognition algorithms to demonstrate the viability of the flow control method. A flat plate towed in a water channel is used as a test bed. The instantaneous velocity signal is used to compute important statistical quantities of the random velocity field, such as the mean, the root-mean-square, the spectral distribution, and the probability density function. When optimized for a given boundary layer, it is shown that the cyclic jet will produce periodic structures that are similar to the random, naturally occurring ones. These structures seem to trigger the onset of bursting events near the wall of the plate. Thus, the present device generates periodic structures in both the outer and inner regions of a turbulent boundary layer.
The motion of charged particles in a stochastic magnetic field with nonzero mean is considered via a generalized quasi-linear expansion of Liouville's equation. The general result is an equation relating cosmic ray scintillations to magnetic fluctuations and to cosmic ray gradients. The resonant interaction between particles and the random magnetic field is considered in detail, and the effect of nonlinear terms in the equations is considered. The nonlinear terms are important in damping out initial conditions and in determining conditions near cyclotron resonances. The application of the theory to the propagation of cosmic rays during quiet times in interplanetary space is considered. It is concluded that cosmic ray scintillations in interplanetary space may provide useful information about interplanetary particles and fields and also about nonlinear plasma interactions.
Recent high-energy gamma-ray observations (E(sub gamma) greater than 100 MeV) of blazar Active Galactic Nuclei (AGNs) show emission spectra with no clear upper energy cutoff. AGNs, considered to be possible sources for the observed flux of cosmic rays beyond 10(exp 19) eV, may well have emission extending into the very high energy (VHE), (E(sub gamma) greater than 100 GeV) domain. Because VHE gamma-rays are absorbed by pair production on the intergalactic background radiation fields, much of this emission may not be directly visible. The electromagnetic cascades initiated by absorbed VHE gamma-rays, however, may be observable. Since, most probably, the velocities of (e(+), e(-)) pairs produced in a cascade are quickly isotropized by an ambient random magnetic field, extended 'halos' (R greater than 1 Mpc) of pairs will be formed around AGNs with VHE emission. The cascade radiation from these pair halos is emitted isotropically and should be observable at energies below a few TeV. The halo radiation can be distinguished by its characteristic variation in spectrum and intensity with angular distance from the central source. This variation depends weakly on the details of the central source model, e.g., the orientation and beaming/opening angle of an emitting jet. Limiting or determining the intensity of the pair halo can thus serve as a model-independent bound on or measure of the VHE power of AGNs. Next-generation Cherenkov telescopes may be able to image a pair halo.
Numerical linear algebra (NLA) underpins huge swaths of computational science and engineering. For scientists and engineers to make the most of the DOE’s computing resources, it is essential that they have access to high-performance implementations of algorithms with best-in-class scalability and reliability. Despite this, prevailing NLA libraries have little to no support for breakthrough algorithms from the field of randomized numerical linear algebra (RandNLA) that have been developed over the past twenty years. The goal of this LDRD was to break a log-jam that had prevented broad adoption of RandNLA. Our work had two thrusts. The first was to develop RandBLAS: a trustworthy and high-performance C++ library for randomized dimension reduction (an operation widely known as sketching). The second was the development of a novel randomized algorithm for computing a challenging type of matrix decomposition known as Householder QR with column pivoting (Householder QRCP). In this one-year late-start LDRD we successfully delivered RandBLAS 1.0 and new CPU and GPU codes for Householder QRCP. RandBLAS has extensive documentation at https://randblas.readthedocs.io/en/stable/. Papers on RandBLAS and and our high-performance QRCP codes are forthcoming.