Results for “DISTRIBUTION FUNCTIONS” — Page 24

NASA NTRS1987

Comparison of the structure of harmonic aqueous glasses and liquid water

Glassy structures of water were generated by rapidly quenching configurations of 64 and 343 molecules of liquid water; and the potential energy was then expanded through quadratic order around local minima generated in this way, and the properties of the resulting harmonic system were calculated. The results obtained were used to test the degree to which the structure of liquid water is similar to that of a harmonic aqueous glass. The radial distribution functions for the glass were found to be remarkably similar to those of the liquid. Translational modes were found to be almost entirely responsible for the broadening of the oxygen-oxygen radial distribution function of the quenched configuration, and also primarily responsible for the broadening of other radial distribution functions.

Pohorille, Andrew↗

NASA NTRS1993

Compton scattering and the gamma-ray power-law spectrum in Markarian 421

The nearest BL Lac object, Mrk 421, has a gamma-ray spectrum which is approximately flat in EF-sub E from E less than about 50 MeV to E greater than about 1 TeV. Inverse Compton scattering can explain this smooth spectrum, despite the structure in the Klein-Nishina cross section, if the injected electron distribution function is proportional to gamma exp -2, where gamma is the electron Lorentz factor. When this is the case, the structure imprinted on the steady state electron distribution function by the structure in the Klein-Nishina cross section is almost exactly compensated in the radiated spectrum. Because particle acceleration in strong shocks injects particles with this distribution function, this shape injection function is in fact quite plausible. Other blazars may be explained by the same model if the cutoff below TeV energies observed in other objects is due to pair production on the cosmological IR background, as suggested by Stecker et al. (1992).

Zdziarski, Andrzej A.↗

DOE OSTI2022

Kinetic simulations of collision-less plasmas in open magnetic geometries

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (a) scrape-of-layer (SOL) and divertor regions in toroidal confinement fusion devices, (b) linear divertor simulators, (c) plasma-based thrusters and (d) magnetic mirrors etc. A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic particle-in-cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (a) the effect of non-uniform magnetic fields, (b) finite fully-absorbing boundaries for the particles and (c) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in a state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. The results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in a manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, the inclusion of a Monte Carlo based Fokker-Planck collision operator is discussed in the context of future work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1974

Cosmic-ray streaming perpendicular to the mean magnetic field

Starting from a quasi-linear approximation for the ensemble-averaged particle distribution function in a random magnetic field, the complete diffusion tensor is derived. This is done by assuming a simple form for the ensemble-averaged distribution function, explicitly retaining all components of the streaming flux. This derivation obtains the antisymmetric terms in a natural manner. The necessary dropping of higher-order terms gives a criterion for the lower-energy limit of validity of the perpendicular and antisymmetric diffusion coefficients. The limit for the assumed distribution function is about 0.8 GV rigidity in the interplanetary field near 1 AU.

Forman, M. A.↗

DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1986

High-frequency waves generated by auroral electrons

Measurements of marginally unstable electron distribution functions and high-frequency plasma waves were made on a sounding rocket flight through a quiet auroral arc. The waves appeared near the electron plasma frequency and had a large parallel electric field component such that k-parallel is greater than k-perpendicular. The appearance of these waves was correlated with the presence of marginally unstable parallel electron distributions. Analysis has shown that the waves were produced by parallel electron distribution function greater than 0 rather than the small perpendicular electron distribution function greater than 0 features. Wave levels and growth rates inside the arc were small, and nonlinear wave-wave and wave-particle interactions appear to have been minimal.

Mcfadden, J. P.↗

NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, r d a U production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembe1 (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and platelike), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and in the mid-latitude continent with different concentrations of CCN: a low "c1ean"concentration and a high "dirty" concentration. In addition, differences and similarities between bulk microphysics and spectral-bin microphysical schemes will be examined and discussed.

Tao, Wei-Kuo↗

NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensembel (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and in the mid-latitude continent with different concentrations of CCN: a low "c1ean"concentration and a high "dirty" concentration. In addition, differences and similarities between bulk microphysics and spectral-bin microphysical schemes will be examined and discussed.

Tao, Wei-Kuo↗

NASA NTRS2004

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, two detailed spectral-bin microphysical schemes were implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles (i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail). Each type is described by a special size distribution function containing many categories (i.e. 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. A spectral-bin microphysical model is very expensive from a computational point of view and has only been implemented into the 2D version of the GCE at the present time. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region, in the sub-tropics (Florida) and in the mid-latitude using identical thermodynamic conditions but with different concentrations of CCN: a low 'clean' concentration and a high 'dirty' concentration.

Tao, Wei-Kuo↗

NASA NTRS2005

The Impact of Aerosols on Cloud and Precipitation Processes: Cloud-Resolving Model Simulations

Cloud microphysics are inevitable affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds, Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effect of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, a detailed spectral-bin microphysical scheme was implemented into the the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bim microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.

Tao, Wei-Kuo↗

NASA NTRS2006

The Role of Aerosols on Precipitation Processes

Cloud physics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distribution below the clouds. Therefore, the size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, a detailed spectral--bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.

Tao, Wei-Kuo↗

NASA NTRS2012

The Role of Aerosols on Precipitation Processes: Cloud Resolving Model Simulations

Cloud microphysics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds. Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral bin microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e. pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions. The model is tested by studying the evolution of deep cloud systems in the west Pacific warm pool region, the sub-tropics (Florida) and midlatitudes using identical thermodynamic conditions but with different concentrations of CCN: a low "clean" concentration and a high "dirty" concentration. Results indicate that the low CCN concentration case produces rainfall at the surface sooner than the high CeN case but has less cloud water mass aloft. Because the spectral-bin model explicitly calculates and allows for the examination of both the mass and number concentration of species in each size category, a detailed analysis of the instantaneous size spectrum can be obtained for these cases. It is shown that since the low (CN case produces fewer droplets, larger sizes develop due to greater condensational and collection growth, leading to a broader size spectrum in comparison to the high CCN case. Sensitivity tests were performed to identify the impact of ice processes, radiation and large-scale influence on cloud-aerosol interactive processes, especially regarding surface rainfall amounts and characteristics (i.e., heavy or convective versus light or stratiform types). In addition, an inert tracer was included to follow the vertical redistribution of aerosols by cloud processes. We will also give a brief review from observational evidence on the role of aerosol on precipitation processes.

Tao, Wei-Kuo↗

NASA NTRS1984

Path-integrated growth of auroral kilometric radiation

Using Poeverlein's graphical method, three dimensional ray path calculations are performed to evaluate the path-integrated growth of auroral kilometric radiation (AKR). The ray tracing results indicate that waves whose initial wave vector lie in the local meridian plane continue to propagate in that plane and that among these waves, those with frequencies near the cutoff frequency (f sub R = 0) refract substantially, where as those with frequencies well above the cutoff frequency suffer little refraction. It is also shown that waves whose initial wave vector lie outside of the local meridian plane propagate in the longitudinal as well as the radial and the latitudinal directions. The refraction of these waves is also highly dependent upon the wave frequency, i.e., waves with frequencies near f sub R = 0 refract substantially, whereas waves with frequencies much above f sub R = 0 undergo little refraction. In order to test the electron cyclotron maser mechanisms as a method for generation of AKR, a typical electron distribution function measured in the auroral zone by the S3-3 satellite, is used to calculate path-integrated growths of representative rays. The results of this study indicate that electron distribution functions like those measured by the S3-3 satellite are not capable of amplifying cosmic noise background to the observed intensities of auroral kilometric radiation, and that much steeper slopes at the edges of the loss cone are required. The presence of such distribution functions in the auroral zone is plausible if one assumes that backscattered electrons in this region have energies less than a few hundred eV.

Omidi, N.↗

DOE OSTI2022

Mesh objective stochastic simulations of quasibrittle fracture

Continuum finite element (FE) modeling of damage and failure of quasibrittle structures suffers from the spurious mesh sensitivity due to strain localization. Here this issue has been addressed for deterministic analysis through the development of localization limiters. Here this study proposes a mechanism-based model to mitigate the mesh sensitivity in stochastic FE simulations of quasibrittle fracture. The interest is placed on the analysis of large-size structures, where the mesh size is conveniently chosen to be larger than the width of the fracture process zone as well as the correlation length of the random fields of constitutive properties. The present model is formulated within the framework of continuum damage mechanics. Two localization parameters are introduced to describe the evolution of the damage pattern of each finite element. These parameters are used to guide the energy regularization of the constitutive law, as well as to formulate the mesh-dependent probability distributions of constitutive properties. Depending on the prevailing damage pattern, different energy regularization schemes and mesh dependence of the probability distribution functions are used in the constitutive law. The model is applied to simulate the stochastic failure behavior of quasibrittle structures of different geometries featuring different failure processes including damage initiation, localization, and propagation. It is shown that using fixed probability distribution functions of constitutive properties could lead to strong mesh dependence of the prediction of the mean and variance of the structural load capacity. The probability distribution functions of constitutive properties must be linked to the damage pattern, which may evolve during the failure process. Such a mechanism-based modeling of the probability distributions of constitutive properties is essential for mitigating the spurious mesh sensitivity in stochastic FE analysis of quasibrittle fracture.

42 ENGINEERING↗

NASA NTRS1990

Observational signatures of nonlinear magnetotail particle dynamics

It has been predicted that the nonlinear particle dynamics in the magnetotail leads to a class of resonance structures in the quiet-time ion distribution functions in the central plasma sheet (CPS). These structures exhibit a scaling law of H to the 1/4 power. The first identification of such a scaling law in quiet-time CPS distribution functions obtained aboard ISEE 1 is reported. A method is proposed which uses the observed resonance structures to determine the quiet-time current sheet thickness based on a measurement of the distribution function and magnetic field obtained by one satellite.

Chen, James↗

NASA NTRS1992

Fuzzy operators and cyclic behavior in formal neuronal networks

Formal neuronal networks (FNN), which are comprised of threshold gates, make use of the unit step function. It is regarded as a degenerated distribution function (DDF) and will be referred to here as a non-fuzzy threshold operator (nFTO). Special networks of this kind generating long cycles of states are modified by introduction of fuzzy threshold operators (FTO), i.e., non-degenerated distribution functions (nDDF). The cyclic behavior of the new nets is compared with the original ones. The interconnection matrix and threshold values are not modified. It is concluded that the original long cycles change the fixed points and short cycles, and as the computer simulations demonstrate, the aperiodic motion that is associated with chaotic behavior appears. The emergence of the above changes depend on the steepness of the threshold operators.

Labos, E.↗

DOE OSTI2025

New Measurements of the Deuteron-to-Proton 𝐹 2 Structure-Function Ratio

Nucleon structure functions, as measured in lepton-nucleon scattering, have historically provided a critical observable in the study of partonic dynamics within the nucleon. However, at very large parton momenta, it is both experimentally and theoretically challenging to extract parton distributions due to the probable onset of nonperturbative contributions and the unavailability of high-precision data at critical kinematics. Extraction of the neutron structure and the d quark distribution have been further challenging because of the necessity of applying nuclear corrections when utilizing scattering data from a deuteron target to extract the free neutron structure. However, a program of experiments has been carried out recently at the energy-upgraded Jefferson Lab electron accelerator aimed at significantly reducing the nuclear correction uncertainties on the d quark distribution function at large partonic momentum. This allows leveraging the vast body of deuterium data covering a large kinematic range to be utilized for d quark parton distribution function extraction. In this Letter, we present new data from experiment E12-10-002, carried out in Jefferson Lab Experimental Hall C, on the deuteron to proton cross section ratio at large Bjorken 𝑥. These results significantly improve the precision of existing data and provide a first look at the expected impact on quark distributions extracted from parton distribution function fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2004

Representation of the Geosynchronous Plasma Environment for Spacecraft Charging Calculations

Historically, our ability to predict and postdict surface charging has suffered from both a lack of reliable secondary emission and backscattered electron yields and poor characterization of the plasma environment. One difficulty lies in the common practice of fitting the plasma data to a Maxwellian or Double Maxwellian distribution function, which may not represent the data well for charging purposes. For 13 years Los Alamos National Laboratory (LANL) has been accumulating measurements of electron and proton spectra from Magnetospheric Plasma Analyzer (MPA) instruments aboard a series of geosynchronous satellites. These data provide both a plasma characterization and the potential of the instrument ground. We use electron and ion flux spectra measured by the LANL MPA to examine how the use of different spectral representations of the charged particle environment in computations of spacecraft potentials during magnetospheric substorms affects the accuracy of the results. We calculate the spacecraft potential using both the measured fluxes and several different fits to these fluxes. These flux measurements and fits have been corrected for the difference between the measured and calculated potential. The potentials computed using the measured fluxes, the best available material properties of graphite carbon, and a secondary electron escape fraction of 81%, are within a factor of three of the measured potential for nearly all the data. Using a Kappa fit to the electron distribution function and a Maxwellian fit to the ion distribution function gives agreement similar to the calculations using the actual data. Alternative spectral representations, including Maxwellian and double Maxwellian for both species, lead to less satisfactory agreement between predicted and measured potentials.

Davis, V. A.↗