Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution function”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
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.↗
DOE OSTI2020

Machine learning approaches for structural and thermodynamic properties of a Lennard-Jones fluid

Predicting the functional properties of many molecular systems relies on understanding how atomistic interactions give rise to macroscale observables. However, current attempts to develop predictive models for the structural and thermodynamic properties of condensed-phase systems often rely on extensive parameter fitting to empirically selected functional forms whose effectiveness is limited to a narrow range of physical conditions. Here, we illustrate how these traditional fitting paradigms can be superseded using machine learning. Specifically, we use the results of molecular dynamics simulations to train machine learning protocols that are able to produce the radial distribution function, pressure, and internal energy of a Lennard-Jones fluid with increased accuracy in comparison to previous theoretical methods. The radial distribution function is determined using a variant of the segmented linear regression with the multivariate function decomposition approach developed by Craven et al. [J. Phys. Chem. Lett. 11, 4372 (2020)]. The pressure and internal energy are determined using expressions containing the learned radial distribution function and also a kernel ridge regression process that is trained directly on thermodynamic properties measured in simulation. The presented results suggest that the structural and thermodynamic properties of fluids may be determined more accurately through machine learning than through human-guided functional forms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2019

Behavior of Triple Langmuir Probes in Non-Equilibrium Plasmas

A model of current collection in Langmuir probes is used to investigate the effect representative non- equilibrium plasmas under various conditions have on the electron temperature and number density that would be calculated through analysis of the probe collection characteristics. The model uses the distribution function to calculate the charged particle flux to a probe and then, for fixed applied voltages, the current values that satisfy continuity in the probes are determined. The triple probe is not scanned in voltage, so there is no practical way using experimental triple probe data to determine if the plasma is in equilibrium. As a consequence, a triple probe analysis typically relies on the assumption that the plasma is in equilibrium. Proceeding from this point, the numerically-generated non-equilibrium triple probe data are analyzed assuming that the plasma is in equilibrium, with the data compared to the initial distribution function inputs of plasma temperature and number density to determine the effect the non-equilibrium distribution has on plasma measurements. The temperature and number density are both significantly affected when a fraction of the particles in the distribution are shifted from the equilibrium configuration into the non-equilibrium part of the distribution function. For all instances studied, the computed electron temperature and number density are extremely sensitive to small deviations from equilibrium ( ≤ 5% of the plasma shifted into the non-equilibrium function). Shifting more of the plasma into the non-equilibrium distribution beyond this initial level does not produce a significant additional shift in the computed plasma properties.

Polzin, Kurt A.↗
NASA NTRS1991

The sensitivity of fracture location distribution in brittle materials

The Weibull weak-link theory allows for the computation of distribution functions for both the fracture location and the applied far-field stress. Several authors have suggested using the fracture location information from tests to infer Weibull parameters, and others have used the predictive capabilities of the theory to calculate average fracture locations for brittle bodies. By a simple set of example calculations, it is shown that the fracture location distribution function is distinctly more sensitive to perturbations in the stress state than the fracture stress distribution function is. In general, the average fracture location is more subject to stress perturbations than the average fracture stress. The results indicate that care must be exercised in applying fracture location theory.

Wetherhold, Robert C.↗
NASA NTRS1982

Comparison of transport equations based on Maxwellian and bi-Maxwellian distributions for anisotropic plasmas

A wide variety of plasma flow conditions is found in aeronomy and space plasma physics. Transport equations based on an isotropic Maxwellian vilecity distribution function can be used to describe plasma flows which contain 'small' temperature anisotropies. However, for plasma flows characterized by large temperature anisotropies, transport equations based on an anisotropic bi-Maxwellian (or two-temperature) velocity distribution function are expected to provide a much better description of the plasma transport properties. The present investigation is concerned with the extent to which transport equations based on both Maxwellian and bi-Maxwellian series expansions can describe plasma flows characterized by non-Maxwellian velocity distributions, giving particular attention to a modelling of the anisotropic character of the distribution function. The obtained results should provide clues as to the extent to which a given series expansion can account for the anisotropic character of a plasma.

Barakat, A. R.↗
DOE OSTI2020

Heavy-light mesons on the light front

Abstract We study the heavy-light mesons within basis light-front quantization. The resulting mass spectra of D , $$D_s$$ D s , B , and $$B_s$$ B s agree reasonably well with experiments. We also predict states which could be measured in the near future. In the light-front formalism, we calculate the light-front wave functions and additional experimental observables, such as parton distribution functions, distribution amplitudes, and decay constants by means of integrations over light-front wave functions. We also provide ratios of decay constants for selected pseudoscalar meson decays ( $$D_s$$ D s to D and $$B_s$$ B s to B ) as they may prove to be theoretically more robust and more reliably determined in experiments. We find that our ratios are systematically smaller than existing experiments and other approaches by 5–18%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1983

Plasma electron analysis: Voyager plasma science experiment

The Plasma Science Experiment (PLS) on the Voyager spacecraft provide data on the plasma ions and electrons in the interplanetary medium and the magnetospheres of the giant planets Jupiter and Saturn. A description of the analysis used to obtain electron parameters (density, temperature, etc.) from the plasma science experiment PLS electron measurements which cover the energy range from 10 eV to 5950 eV is presented. The electron sensor (D cup) and its transmission characteristics are described. A derivation of the fundamental analytical expression of the reduced distribution function F(e) is given. The electron distribution function F(e), used in the moment integrations, can be derived from F(e). Positive ions produce a correction current (ion feedthrough) to the measured electron current, which can be important to the measurements of the suprathermal electron component. In the case of Saturn, this correction current, which can either add to or subtract from the measured electron current, is less than 20% of the measured signal at all times. Comments about the corrections introduced by spacecraft charging to the Saturn encounter data, which can be important in regions of high density and shadow when the spacecraft can become negatively charged are introduced.

Sittler, E. C., Jr.↗
NASA NTRS1985

The propagation and growth of whistler mode waves generated by electron beams in earth's bow shock

In this study, the propagation and growth of whistler mode waves generated by electron beams within earth's bow shock is investigated using a planar model for the bow shock and a model electron distribution function. Within the shock, the model electron distribution function possesses a field-aligned T greater than T beam that is directed toward the magnetosheath. Waves with frequencies between about 1 and 100 Hz with a wide range of wave normal angles are generated by the beam via Landau and anomalous cyclotron resonances. However, because the growth rate is small and because the wave packets traverse the shock quickly, these waves do not attain large amplitudes. Waves with frequencies between about 30 and 150 Hz with a wide range of wave normal angles are generated by the beam via the normal cyclotron resonance. The ray paths for most of these waves are directed toward the solar wind although some wave packets, due to plasma convection travel transverse to the shock normal. These wave packets grow to large amplitudes because they spend a long time in the growth region. The results suggest that whistler mode noise within the shock should increase in amplitude with increasing upstream theta sub Bn. The study provides an explanation for the origin of much of the whistler mode turbulence observed at the bow shock.

Tokar, R. L.↗
NASA NTRS1991

Chaotic jumps in the generalized first adiabatic invariant in current sheets

The present study examines how the changes in the generalized first adiabatic invariant J derived from the separatrix crossing theory can be incorporated into the drift variable approach to generating distribution functions. A method is proposed for determining distribution functions for an ensemble of particles following interaction with the tail current sheet by treating the interaction as a scattering problem characterized by changes in the invariant. Generalized drift velocities are obtained for a 1D tail configuration by using the generalized first invariant. The invariant remained constant except for the discrete changes caused by chaotic scattering as the particles cross the separatrix.

Brittnacher, M. J.↗
NASA NTRS1973

Airborne Laser Polar Nephelometer

A polar nephelometer has been developed at NCAR to measure the angular variation of the intensity of light scattered by air molecules and particles. The system has been designed for airborne measurements using outside air ducted through a 5-cm diameter airflow tube; the sample volume is that which is common to the intersection of a collimated source beam and the detector field of view within the airflow tube. The source is a linearly polarized helium-neon laser beam. The optical system defines a collimated field-of-view (0.5deg half-angle) through a series of diaphragms located behind a I72-mm focal length objective lens. A photomultiplier tube is located immediately behind an aperture in the focal plane of the objective lens. The laser beam is mechanically chopped (on-off) at a rate of 5 Hz; a two-channel pulse counter, synchronized to the laser output, measures the photomultiplier pulse rate with the light beam both on and off. The difference in these measured pulse rates is directly proportional to the intensity of the scattered light from the volume common to the intersection of the laser beam and the detector field-of-view. Measurements can be made at scattering angles from 15deg to 165deg with reference to the direction of propagation of the light beam. Intermediate angles are obtained by selecting the angular increments desired between these extreme angles (any multiple of 0.1deg can be selected for the angular increment; 5deg is used in normal operation). Pulses provided by digital circuits control a stepping motor which sequentially rotates the detector by pre-selected angular increments. The synchronous photon-counting system automatically begins measurement of the scattered-light intensity immediately after the rotation to a new angle has been completed. The instrument has been flown on the NASA Convair 990 airborne laboratory to obtain data on the complex index of refraction of atmospheric aerosols. A particle impaction device is operated simultaneously to collect particles from the same airflow tube used to make the scattered-light measurements. A size distribution function is obtained by analysis of the particles collected by the impaction device. Calculated values of the angular variation of the scattered-light intensity are obtained by applying Mie scattering theory to the observed size distribution function and assuming different values of the complex index of refraction of the particles. The calculated values are then compared with data on the actual variation of the scattered-light intensity obtained with the polar nephelometer. The most probable value of the complex refractive index is that which provides the best fit between the experimental light scattering data and the parameters calculated from the observed size distribution function.

Grams, Gerald W.↗
NASA NTRS2007

A Simple Stochastic Model for Generating Broken Cloud Optical Depth and Top Height Fields

A simple and fast algorithm for generating two correlated stochastic twodimensional (2D) cloud fields is described. The algorithm is illustrated with two broken cumulus cloud fields: cloud optical depth and cloud top height retrieved from Moderate Resolution Imaging Spectrometer (MODIS). Only two 2D fields are required as an input. The algorithm output is statistical realizations of these two fields with approximately the same correlation and joint distribution functions as the original ones. The major assumption of the algorithm is statistical isotropy of the fields. In contrast to fractals and the Fourier filtering methods frequently used for stochastic cloud modeling, the proposed method is based on spectral models of homogeneous random fields. For keeping the same probability density function as the (first) original field, the method of inverse distribution function is used. When the spatial distribution of the first field has been generated, a realization of the correlated second field is simulated using a conditional distribution matrix. This paper is served as a theoretical justification to the publicly available software that has been recently released by the authors and can be freely downloaded from http://i3rc.gsfc.nasa.gov/Public codes clouds.htm. Though 2D rather than full 3D, stochastic realizations of two correlated cloud fields that mimic statistics of given fields have proved to be very useful to study 3D radiative transfer features of broken cumulus clouds for better understanding of shortwave radiation and interpretation of the remote sensing retrievals.

Prigarin, Sergei M.↗