Diffuse reflection and transmission by cloud and dust layers.
Radiative transfer in anisotropic phase function medium using Neumann solution, Romanova method and doubling method, noting probability distribution of photon optical paths
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.
Radiative transfer in anisotropic phase function medium using Neumann solution, Romanova method and doubling method, noting probability distribution of photon optical paths
Multiple light scattering in planetary atmospheres, discussing diffuse reflection and transmission by atmosphere of particles with anisotropic scattering pattern
Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.
When charged particles spiral along a large constant magnetic field, their trajectories are scattered by any random field components that are superposed on the guiding field. If the random field configuration embodies helicity, the scattering is asymmetrical with respect to a plane perpendicular to the guiding field, for particles moving into the forward hemisphere are scattered at different rates from those moving into the backward hemisphere. This asymmetry gives rise to new terms in the transport equations that describe propagation of charged particles. Helicity has virtually no impact on qualitative features of the diffusive mode of propagation. However, characteristic velocities of the coherent modes that appear after a highly anisotropic injection exhibit an asymmetry related to helicity. Explicit formulas, which embody the effects of helicity, are given for the anisotropies, the coefficient diffusion, and the coherent velocities. Predictions derived from these expressions are in good agreement with Monte Carlo simulations of particle transport, but the simulations reveal certain phenomena whose explanation calls for further analytical work.
For ferroelectric ionic-conductors, polarization switching is complicated by the interplay between ion mobility and charge screening effects. When the ionic charge carriers also play a key role in the domain reversal, such as in Rb-doped KTiOPO 4 (RKTP), a higher level of complexity is introduced. RKTP provides an ideal platform for investigating the relationship between ionic conductivity and polarization reversal because its highly anisotropic crystal properties allow selective modification of material characteristics through diffusive cation doping. Here, we use indiffused Ba/K doping to create a significant increase in the ionic conductivity. Time-of-flight secondary ion mass spectrometry is employed to map Ba/K doping within the RKTP crystal and correlate it to changes in ionic mobility and polarization switching characteristics under an external field applied to the nonpolar face. Using band-excitation piezoresponse force microscopy, we demonstrate a selective switching-inhibition mechanism driven by the enhanced charge screening.
The properties of the suprathermal particle distributions observed upstream of interplanetary shocks depend not only on the properties of the shocks but also on the transport conditions encountered by the particles as they propagate away from the shocks. The confinement of particles in close proximity to the shocks, as well as particle scattering processes during propagation to the spacecraft, lead to the common observation of upstream diffuse particle distributions. We present observations of a rare extended anisotropic low-energy(≤30 keV)proton beam together with a trapped ≥500 keV proton population observed in association with the arrival of an oblique interplanetary shock at the Advanced Composition Explorer, the Interplanetary Monitoring Platform-8, and the Wind spacecraft on 2001 January 31. Continuous injection of particles by the traveling shock into a smooth radial magnetic field region formed in the tail of a modest high-speed solar wind stream produced an extended foreshock region of energetic particles. The absence of enhanced magnetic field fluctuations upstream of the shock results in the observation of a prolonged anisotropic field-aligned beam of ≤30 keV protons as well as a population of higher-energy(≥500 keV)protons with small pitch-angle cosine(μ∼0)extending far from the shock.
The effects are studied of assumed random velocity fields on diffusion in a binary fluid. Random velocity fields can result, for example, from the high-frequency components of residual accelerations onboard spacecraft (often called g-jitter). An effective diffusion equation is derived for an average concentration which includes spatial and temporal correlations induced by the fluctuating velocity fields assumed to be Gaussianly distributed. The resulting equation becomes nonlocal, and if correlations between different components of the velocity field exist, it is also anisotropic. The simple limiting case of short correlation times is discussed and an effective diffusivity is obtained which reflects the enhanced mixing caused by the velocity fields. The results obtained in the limit of short correlation times are valid even if the probability distribution of the velocity field is not Gaussian.
Detailed analysis of electrons equal to or greater than 3 MeV and of protons 0.5 to 1.8 MeV and equal to or greater than 35 MeV for both the inbound and the outbound passes of the Pioneer 10 spacecraft. Conclusive evidence is obtained that the trapped radiation in Jupiter's inner magnetosphere is maintained and supplied by inward diffusion from the outer regions of the trapped radiation zone. It is shown that the time required for isotropization of an anisotropic flux by pitch angle scattering inside L approximately equal to 6 is long in comparison with the time required for particles to diffuse inward from L approximately equal to 6 to L approximately equal to 3, that the high-energy protons were not injected at high energies by the Crand (cosmic ray albedo neutron decay) process but were accelerated in the magnetosphere of Jupiter, and that the main conclusions of this analysis are unaffected by use of either the D sub 1 or the D sub 2 magnetic field models. Theoretical studies of the capture of trapped electrons and protons by Io have been carried out, and it is found that the probability of capture by Io depends strongly upon the particle species and kinetic energy.
Graphite intercalation compounds may provide an excellent material for the fabrication of electro-magnetic railgun armatures. As a pulse of power is fed into the armature the intercalate could be excited into the plasma state around the edges of the armature, while the bulk of the current would be carried through the graphite block. Such an armature would have the desirable characteristics of both diffuse plasma armatures and bulk conduction armatures. In addition, the highly anisotropic nature of these materials could enable the electrical and thermal conductivity to be tailored to meet the specific requirements of electromagnetic railgun armatures. Preliminary investigations were performed in an attempt to determine the feasibility of using graphite intercalation compounds as railgun armatures. Issues of fabrication, resistivity, stability, and electrical current spreading are addressed for the case of highly oriented pyrolytic graphite.
Graphite intercalation compounds may provide an excellent material for the fabrication of electromagnetic railgun armatures. As a pulse of power is fed into the armature the intercalate could be excited into the plasma state around the edges of the armature, while the bulk of the current would be carried through the graphite block. Such an armature would have both diffuse plasma armatures and bulk conduction armatures. In addition, the highly anisotropic nature of these materials could enable the electrical and thermal conductivity to be tailored to meet the specific requirements of electromagnetic railgun armatures. Preliminary investigations have been performed in an attempt to determine the feasibility of using graphite intercalation compounds as railgun armatures. Issues of fabrication, resistivity, stability, and electrical current spreading have been addressed for the case of highly oriented pyrolytic graphite.
Anisotropic unstructured mesh adaption is developed for a truly multi-dimensional upwind fluctuation splitting scheme, as applied to scalar advection-diffusion. The adaption is performed locally using edge swapping, point insertion/deletion, and nodal displacements. Comparisons are made versus the current state of the art for aggressive anisotropic unstructured adaption, which is based on a posteriori error estimates. Demonstration of both schemes to model problems, with features representative of compressible gas dynamics, show the present method to be superior to the a posteriori adaption for linear advection. The performance of the two methods is more similar when applied to nonlinear advection, with a difference in the treatment of shocks. The a posteriori adaption can excessively cluster points to a shock, while the present multi-dimensional scheme tends to merely align with a shock, using fewer nodes. As a consequence of this alignment tendency, an implementation of eigenvalue limiting for the suppression of expansion shocks is developed for the multi-dimensional distribution scheme. The differences in the treatment of shocks by the adaption schemes, along with the inherently low levels of artificial dissipation in the fluctuation splitting solver, suggest the present method is a strong candidate for applications to compressible gas dynamics.
The pressure recovery coefficient of a two-dimensional diffuser (Cp) was found to be significantly affected by the intensity and structure of the inlet free-stream turbulence. Significantly larger values of Cp occurred for the case of highly anisotropic inlet free-stream turbulence with the plane of rotation of the eddies parallel to the free-stream flow and perpendicular to the diverging walls of the diffuser, as compared to the values of Cp for the case of grid-generated near-isotropic inlet conditions. The relatively high values of Cp for the case of the highly anisotropic inlet free-stream turbulence conditions can be attributed to highly effective mixing, where the plane of rotation of the eddies facilitates a direct transport of momentum from the free-stream to the boundary layers along the diverging walls of the diffuser, as compared to the mixing that occurs with grid-generated highly unorganized eddies.
The electron diffusion region forms the inner core of the magnetic reconnection machine. In this region, physical processes act to sustain a reconnection electric field, they accelerate to provide current density, and the heat to provide particle pressure in the current layer. The relative roles of these processes depend on whether the overall geometry is anti-parallel or asymmetric, and whether a guide field is present of not, but the overall concept applies in each case. The EDR is embedded in successively larger regions, starting with the ion diffusion region to even larger scales, where the plasma increasingly behaves like an anisotropic MHD plasma. With our new understanding of how the EDR works, it is of great interest to study how EDR processes couple to the larger scale environment. The ultimate efficacy of the magnetic reconnection is shaped by this interesting example of micro-macro coupling. In this presentation, we present some thoughts pertaining to this topic. In particular, we will discuss how information transport between scales is mediated, and how the very small scales can reach balance with the larger, overall system. We will review any pertinent research and present a set of questions inviting future research.
PLoM (Probabilistic Learning on Manifolds) is a method introduced in 2016 for handling small training datasets by projecting an Itô equation from a stochastic dissipative Hamiltonian dynamical system, acting as the MCMC generator, for which the KDE-estimated probability measure with the training dataset is the invariant measure. PLoM performs a projection on a reduced-order vector basis related to the training dataset, using the diffusion maps (DMAPS) basis constructed with a time-independent isotropic kernel. In this paper, we propose a new ISDE projection vector basis built from a transient anisotropic kernel, providing an alternative to the DMAPS basis to improve statistical surrogates for stochastic manifolds with heterogeneous data. The construction ensures that for times near the initial time, the DMAPS basis coincides with the transient basis. For larger times, the differences between the two bases are characterized by the angle of their spanned vector subspaces. The optimal instant yielding the optimal transient basis is determined using an estimation of mutual information from Information Theory, which is normalized by the entropy estimation to account for the effects of the number of realizations used in the estimations. Consequently, this new vector basis better represents statistical dependencies in the learned probability measure for any dimension. Three applications with varying levels of statistical complexity and data heterogeneity validate the proposed theory, showing that the transient anisotropic kernel improves the learned probability measure.
Abstract While cosmic rays ( E ≳ 1 GeV) are well coupled to a galaxy’s interstellar medium (ISM) at scales of L > 100 pc, adjusting stratification and driving outflows, their impact on small scales is less clear. Based on calculations of the cosmic-ray diffusion coefficient from observations of the grammage in the Milky Way, cosmic rays have little time to dynamically impact the ISM on those small scales. Using numerical simulations, we explore how more complex cosmic-ray transport could allow cosmic rays to couple to the ISM on small scales. We create a two-zone model of cosmic-ray transport, with the cosmic-ray diffusion coefficient set at the estimated Milky Way value in cold gas but smaller in warm gas. We compare this model to simulations with a constant diffusion coefficient. Quicker diffusion through cold gas allows more cold gas to form compared to a simulation with a constant, small diffusion coefficient. However, slower diffusion in warm gas allows cosmic rays to take energy from the turbulent cascade anisotropically. This cosmic-ray energization comes at the expense of turbulent energy which would otherwise be lost during radiative cooling. Finally, we show our two-zone model is capable of matching observational estimates of the grammage for some transport paths through the simulation.
We have used a three-dimensional model of recovery phase storm hot plasmas to explore the signatures of pitch angle distributions (PADS) in global fast atom imagery of the magnetosphere. The model computes mass, energy, and position-dependent PADs based on drift effects, charge exchange losses, and Coulomb drag. The hot plasma PAD strongly influences both the storm current system carried by the hot plasma and its time evolution. In turn, the PAD is strongly influenced by plasma waves through pitch angle diffusion, a microscale effect. We report the first simulated neutral atom images that account for anisotropic PADs within the hot plasma. They exhibit spatial distribution features that correspond directly to the PADs along the lines of sight. We investigate the use of image brightness distributions along tangent-shell field lines to infer equatorial PADS. In tangent-shell regions with minimal spatial gradients, reasonably accurate PADs are inferred from simulated images. They demonstrate the importance of modeling PADs for image inversion and show that comparisons of models with real storm plasma images will reveal the global effects of these microscale processes.
The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).
Theoretical and experimental work on the fluid dynamics of physical vapor transport is reported. It is shown that diffusion in viscous interaction with container walls leads to concentration gradients normal to the main transport direction. Consequently any convection threshold is removed. This coupling with convective instabilities makes the theoretical interpretation of earthbound (anisotropic) morphological stability studies on vapor-solid interfaces intractable. Hence low-gravity experiments are suggested that will allow for the establishment of morphological stability criteria under diffusion controlled conditions.