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 541 records · Page 30

Optimal Reward Functions in Distributed Reinforcement Learning

We consider the design of multi-agent systems so as to optimize an overall world utility function when (1) those systems lack centralized communication and control, and (2) each agents runs a distinct Reinforcement Learning (RL) algorithm. A crucial issue in such design problems is to initialize/update each agent's private utility function, so as to induce best possible world utility. Traditional 'team game' solutions to this problem sidestep this issue and simply assign to each agent the world utility as its private utility function. In previous work we used the 'Collective Intelligence' framework to derive a better choice of private utility functions, one that results in world utility performance up to orders of magnitude superior to that ensuing from use of the team game utility. In this paper we extend these results. We derive the general class of private utility functions that both are easy for the individual agents to learn and that, if learned well, result in high world utility. We demonstrate experimentally that using these new utility functions can result in significantly improved performance over that of our previously proposed utility, over and above that previous utility's superiority to the conventional team game utility.

Wolpert, David H.↗

The Distribution and Annihilation of Dark Matter Around Black Holes

We use a Monte Carlo code to calculate the geodesic orbits of test particles around Kerr black holes, generating a distribution function of both bound and unbound populations of dark matter (DM) particles. From this distribution function, we calculate annihilation rates and observable gamma-ray spectra for a few simple DM models. The features of these spectra are sensitive to the black hole spin, observer inclination, and detailed properties of the DM annihilation cross-section and density profile. Confirming earlier analytic work, we find that for rapidly spinning black holes, the collisional Penrose process can reach efficiencies exceeding 600%, leading to a high-energy tail in the annihilation spectrum. The high particle density and large proper volume of the region immediately surrounding the horizon ensures that the observed flux from these extreme events is non-negligible.

astroparticle physics↗

S 8 ring content in liquid sulfur determined from x-ray diffraction

The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .

lambda transition↗

Thermionic noise measurements for on-line dispenser cathode diagnostics for linear beam microwave tubes

A test stand has been set up to measure the current fluctuation noise properties of B- and M-type dispenser cathodes in a typical TWT gun structure. Noise techniques were used to determine the work function distribution on the cathode surfaces. Significant differences between the B and M types and significant changes in the work function distribution during activation and life are found. In turn, knowledge of the expected work function can be used to accurately determine the cathode-operating temperatures in a TWT structure. Noise measurements also demonstrate more sensitivity to space charge effects than the Miram method. Full automation of the measurements and computations is now required to speed up data acquisition and reduction. The complete set of equations for the space charge limited diode were programmed so that given four of the five measurable variables (J, J sub O, T, D, and V) the fifth could be computed. Using this program, we estimated that an rms fluctuation in the diode spacing d in the frequency range of 145 Hz about 20 kHz of only about 10 to the -5 power A would account for the observed noise in a space charge limited diode with 1 mm spacing.

Holland, C.↗

Swarm Intelligence Based Optimal Design of Local Volt/Var Control Function for Distributed Energy Resources

The increasing penetration of renewable based distributed energy resources (DERs) in distribution network (DN) leads to larger and more frequent voltage variation in distributions network (DN), thus posing challenges on voltage control. Real-time local voltage control method is a promising solution for the above issue. However, the local voltage control function needs to be customized and optimized according to real distribution system condition. In this paper, a swarm intelligence based Volt/Var control optimal design method (SO-VVC) is proposed to optimize the control function. Compared with existing approaches, the proposed method can not only represent the nonlinear behaviour of power flow but is also computation efficient. The performance of the proposed SO-VVC is demonstrated by case studies on a modified IEEE-123 bus system.

Zhang, Zhengfa [University of Tennessee, Knoxville↗

Evidence of short chains in liquid sulfur

High energy x-ray pair distribution function measurements show the average coordination number of the first shell in liquid sulfur is 1.86 ± 0.04 across the λ-transition, not precisely 2.0 as widely accepted. This indicates that upon melting, liquid sulfur does not comprise solely of S 8 rings but also possesses a significant number of short chains. Intensities of the pre-peak and first diffraction peak of the x-ray structure factor and third peak height of the pair distribution function all show deviations at the λ-transition temperature T λ , associated with the break-up of S 8 rings and the start of oligomer polymerization. A significant number of non-bonded or loosely bonded “interstitial atoms,” with an average coordination number of 0.20 ± 0.005, are also observed in the so-called “forbidden zone” between the first and second shells upon melting. The number of interstitial atoms is found to decrease to a minimum at the λ-transition, but the majority persist into the high temperature polymerized liquid. Furthermore, the existence of short chains and nearby interstitial atoms represent the two main factors required to initiate the S 8 -ring to chain transition, as proposed by recent molecular dynamics simulations.

Chemical bonding↗

X-ray and molecular dynamics study of the temperature-dependent structure of molten NaF-Zr⁢F 4

The local atomic structure of NaF-Zr⁢F 4 (53–47 mol%) molten system and its evolution with temperature are examined with x-ray scattering measurements which are then used to validate the quality of ab initio and neural network-based molecular dynamics (NNMD) calculations in the temperature range 515–700°⁢C. The machine-learning enhanced NNMD calculations offer improved efficiency while maintaining accuracy at higher distances compared to ab initio calculations. Looking at the evolution of the pair distribution function with increasing temperature, a fundamental change in the liquid structure within the selected temperature range, accompanied by a slight decrease in overall correlation is revealed. NNMD calculations indicate the coexistence of three different fluorozirconate complexes: [Zr⁢F 6 ] 2– , [Zr⁢F 7 ] 3– , and [Zr⁢F 8 ] 4– , with a shift in the dominant coordination state from the 7-coordinated Zr cation toward a 6-coordinated cation with increasing temperature. The study also highlights the metastability of different local coordination structures, with frequent interconversions between the 6- and 7-coordinate states. Analysis of the Zr-F-Zr angular distribution function reveals the presence of both “edge-sharing” and “corner-sharing” fluorozirconate complexes with specific bond angles and distances in accord with previous studies, while the next-nearest-neighbor cation-cation correlations demonstrate a clear preference for unlike cations as nearest-neighbor pairs, emphasizing nonrandom arrangement. Finally, these findings contribute to a comprehensive understanding of the complex local structure of the molten salt, providing insights into temperature-dependent preferences and correlations within the molten system.

36 MATERIALS SCIENCE↗

Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics

We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.

First-principles calculations↗

Reconciling calculations and measurements of inverse bremsstrahlung absorption

It was recently shown that the use of Coulomb logarithms appropriate for bremsstrahlung radiation (rather than transport processes) along with corrections for the Langdon effect and ion screening reproduced measurements of collisional absorption in well-characterized underdense plasmas [D. Turnbull et al., Phys. Rev. Lett. 130, 145103 (2023)]. However, it was recognized at the time that the use of the standard absorption-reduction factor from Langdon's seminal paper was inconsistent with the use of Coulomb logarithms that are thermally averaged over a Maxwellian electron-velocity distribution function. A more accurate approach would be to average over the expected super-Gaussian distribution function while accounting for the Gaunt factor's velocity dependence, which somewhat mitigates the Langdon effect; however, at that time, this theory matched the data less well. This conflict is now eliminated with the additional insight that the ionization state of our mid-Z ion species (when present) was lower than had been assumed, as evidenced by the Thomson-scattering data and time-dependent Cretin simulations. We are now able to show that an improved treatment of the Langdon effect provides the best match to data. Otherwise, the prior conclusions remain unchanged. We also show an example of the substantial expected impact to the absorption rate in calculations of indirect-drive hohlraums.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A kinetic line-driven radiation operator and its application to Gyrokinetics

A velocity dependent, kinetic model for line radiation is developed for continuum kinetic codes. It has been implemented in the full-f gyrokinetic code Gkeyll. The total radiation for a charge state is modeled as an advection in velocity space with a form of $\nabla_v \cdot(v\nu(v)f(v))$, guaranteeing particle conservation. The velocity dependence (in the form of an effective frequency $\nu(v)$) is found through fitting the energy loss of the operator, i.e. the second velocity moment, to the radiation data in the OpenADAS database. Therefore, each individual transition does not need to be evaluated every time step, significantly reducing the computational cost of including line radiation in a kinetic model. The dependence on velocity instead of the usual, temperature, allows the radiation to be computed from non-Maxwellian electron distribution functions: We benchmark the model against a collisional radiative model using isotropic non-Maxwellian distribution functions. A velocity dependent model of radiation can more accurately describe the radiation in the more kinetic regimes expected in reactor-scale devices. The velocity dependence qualitatively captures the quantum mechanical need for a minimum velocity before any radiation occurs.

kinetic↗

Solution of the Boltzmann kinetic equation for the relaxation of a gas mixture

The temporal behavior is considered of a quiescent mixture of gases of different temperatures with spatially uniform distribution. The process of heating a cold gas by a hot gas is treated on the basis of the Boltzmann kinetic equation. The mixture is assumed to be composed of absolutely hard smooth spheres, and the initial distribution functions for each gas is taken to the Maxwellian. With such a choice of initial distribution functions, it is shown that the solution of the Boltzmann kinetic equation depends only on the velocity modulus and the time.

Rykov, V. A.↗

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a five-dimensional continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular two-dimensional velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. Here we describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma using first- and second-order elements with the pseudo-inverse method and comparing with the bilinear mapping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC

The global total-f gyrokinetic particle-in-cell code XGC, used to study transport in magnetic fusion plasmas or to couple with a core gyrokinetic code while functioning as an edge gyrokinetic code, implements a 5-dimensional (5D) continuum grid to perform the dissipative operations, such as plasma collisions, or to exchange the particle distribution function information with a core code. To transfer the distribution function between marker particles and a rectangular 2D velocity-space grid, XGC employs a bilinear mapping. The conservation of particle density and momentum is accurate enough in this bilinear operation, but the error in the particle energy conservation can become undesirably large and cause non-negligible numerical heating in a steep edge pedestal. In the present work we update XGC to use a novel mapping technique, based on the calculation of a pseudo-inverse, to exactly preserve moments up to the order of the discretization space. We describe the details of the implementation and we demonstrate the reduced interpolation error for a tokamak test plasma by using 1st- and 2nd-order elements with the pseudo-inverse method and comparing to the bilinear mapping.

Fusion Plasma↗

Energetic electrons and plasma waves associated with a solar type III radio burst

Detailed in situ observations from the ISEE 3 spacecraft of energetic electrons, plasma waves, and radio emission for the type II solar radio burst of February 17, 1979, are presented. The reduced, one-dimensional electron distribution function is constructed as a function of time. Since the faster electrons arrive before the slower ones, a bump on tail distribution forms which is unstable to the growth of Langmuir waves. The plasma wave growth computed from the distribution function agrees well with the observed onset of the Langmuir waves, and there is qualitative agreement between variations in the plasma wave levels and in the development of regions of positive slope in the function. The evolution of the function, however, predicts far higher plasma wave levels than those observed. The maximum levels observed are approximately equal to the threshold for nonlinear wave processes, such as oscillation two-stream instability and soliton collapse.

Lin, R. P.↗

Interaction of field-aligned cold plasma flows with an equatorially-trapped hot plasma - Electrostatic shock formation

Effects of equatorially trapped hot plasma on the highly supersonic cold-plasma flow occurring during early stage plasmaspheric refilling are studied by means of numerical simulations. It is shown that the equatorially trapped hot ions set up a potential barrier for the cold ion beams and facilitate formation of electrostatic shocks by reflecting them from the equatorial region. Simulations with and without the hot plasma show different flow properties; the formation of electrostatic shocks occur only in the former case. The simulation with the hot plasma also reveals that the magnetic trapping in conjunction with the evolution of the electrostatic potential barrier produces ion velocity distribution functions consisting of a cold core and a hot ring in the perpendicular velocity. Such a distribution function provides a source of free energy for equatorial waves. The corresponding electron population is warm and field-aligned.

Singh, Nagendra↗

Research on Spectroscopy, Opacity, and Atmospheres

To make my calculations more readily accessible I have set up a web site cfaku5.harvard.edu that can also be accessed by FTP. it has 5 9GB disks that hold all of my atomic and diatomic molecular data, my tables of distribution function opacities, my grids of model atmospheres, colors, fluxes, etc, my program that are ready for distribution, most of my recent papers. Atlases and computed spectra will be added as they are completed. New atomic and molecular calculations will be added as they are completed. I got my atomic programs that had been running on a Cray at the San Diego Supercomputer Center to run on my Vaxes and Alpha. I started with Ni and Co because there were new laboratory analyses that included isotopic and hyperfine splitting. Those calculations are described in the appended abstract for the 6th Atomic Spectroscopy and oscillator Strengths meeting in Victoria last summer. A surprising finding is that quadrupole transitions have been grossly in error because mixing with higher levels has not been included. I now have enough memory in my Alpha to treat 3000 x 3000 matrices. I now include all levels up through n=9 for Fe I and 11, the spectra for which the most information is available. I am finishing those calculations right now. After Fe I and Fe 11, all other spectra are "easy", and I will be in mass production. ATL;LS12, my opacity sampling program for computing models with arbitrary abundances, has been put on the web server. I wrote a new distribution function opacity program for workstations that replaces the one I used on the Cray at the San Diego Supercomputer Center. Each set of abundances would take 100 Cray hours costing $100,000. 1 ran 25 cases. Each of my opacity CDs contains three abundances. I have a new program -iinning on the Alpha that takes about a week. I am going to have to get a faster processor or I will have to dedicate a whole workstation just to opacities.

Kurucz, Robert L.↗

Properties of Minor Ions in the Solar Wind and Implications for the Background Solar Wind Plasma

Ion charge states measured in situ in interplanetary space are formed in the inner coronal regions below 5 solar radii, hence they carry information on the properties of the solar wind plasma in that region. The plasma parameters that are important in the ion forming processes are the electron density, the electron temperature and the flow speeds of the individual ion species. In addition, if the electron distribution function deviates from a Maxwellian already in the inner corona, then the enhanced tail of that distribution function, also called halo, greatly effects the ion composition. The goal of the proposal is to make use of ion fractions observed in situ in the solar wind to learn about both, the plasma conditions in the inner corona and the expansion and ion formation itself. This study is carried out using solar wind models, coronal observations, and ion fraction calculations in conjunction with the in situ observations.

Esser, Ruth↗

Properties of Minor Ions in the Solar Wind and Implications for the Background Solar Wind Plasma

The scope of the investigation is to extract information on the properties of the bulk solar wind from the minor ion observations that are provided by instruments on board NASA space craft and theoretical model studies. Ion charge states measured in situ in interplanetary space are formed in the inner coronal regions below 5 solar radii, hence they carry information on the properties of the solar wind plasma in that region. The plasma parameters that are important in the ion forming processes are the electron density, the electron temperature and the flow speeds of the individual ion species. In addition, if the electron distribution function deviates from a Maxwellian already in the inner corona, then the enhanced tail of that distribution function, also called halo, greatly effects the ion composition. This study is carried out using solar wind models, coronal observations, and ion calculations in conjunction with the in situ observations.

Wagner, William↗