Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE VELOCITY”

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 1,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

The high explosives & affected targets (HEAT) dataset

Artificial Intelligence (AI) surrogate models offer a computationally efficient alternative to full-physics simulations, yet no existing datasets are publicly available for training, testing, and validation of machine learning models of the dynamics of high-explosive driven shocks through multiple materials. Shock propagation through materials is a computationally challenging problem because simulations must include material-specific equations of state (EOS) along with descriptions of other physical processes such as plastic deformation, phase change, damage processes, fluid instabilities, and multi-material interactions. Shocks are typically initiated by high-velocity impacts or explosive loading. The latter case necessitates the addition of models of reactive materials to represent high-explosive (HE) detonation. Here, to address the lack of an expansive dataset for multi-material shock propagation in the AI/ML community, we present the High-Explosives and Affected Targets (HEAT) Dataset. HEAT is a physics-rich collection of two-dimensional, cylindrically symmetric, simulations generated using an Eulerian, multi-material, shock-propagation code developed at Los Alamos National Laboratory. The dataset includes two partitions: (1) the expanding shock-cylinder (CYL) simulations, Figs. 1, and (2) the Perturbed Layered Interface (PLI) simulations, Fig. 2. Entries in both partitions consist of time series of arrays of thermodynamic fields (pressure, density, and temperature), kinematic fields (position and velocity), and additional fields that depend on thermodynamic and/or kinematic fields (e.g., material stress). Materials in the CYL partition include solids (aluminium, copper, depleted uranium, stainless steel, tantalum, and a generic polymer), a liquid (water), gases (air, nitrogen), and a generic detonating material (high explosive, HE). The PLI partition spans a highly varying geometry but consists of fixed materials across entries: Copper, aluminium, stainless steel, generic polymer, and generic HE. HEAT captures critical phenomena such as momentum transfer, shock propagation, plastic deformation, and thermal effects, making HEAT a valuable benchmark for development of AI/ML emulation of multi-material shock propagation.

36 MATERIALS SCIENCE↗
DOE OSTI2017

A History of H i Stripping in Virgo: A Phase-space View of VIVA Galaxies

We investigate the orbital histories of Virgo galaxies at various stages of H i gas stripping. In particular, we compare the location of galaxies with different H i morphology in phase space. This method is a great tool for tracing the gas stripping histories of galaxies as they fall into the cluster. Most galaxies at the early stage of H i stripping are found in the first infall region of Virgo, while galaxies undergoing active H i stripping mostly appear to be falling in or moving out near the cluster core for the first time. Galaxies with severely stripped, yet symmetric, H i disks are found in one of two locations. Some are deep inside the cluster, but others are found in the cluster outskirts with low orbital velocities. We suggest that the latter group of galaxies belong to a “backsplash” population. These present the clearest candidates for backsplashed galaxies observationally identified to date. We further investigate the distribution of a large sample of H i-detected galaxies toward Virgo in phase space, confirming that most galaxies are stripped of their gas as they settle into the gravitational potential of the cluster. In addition, we discuss the impact of tidal interactions between galaxies and group preprocessing on the H i properties of the cluster galaxies, and link the associated star formation evolution to the stripping sequence of cluster galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Reduced Order Description of Experimental Two-Phase Pipe Flows: Characterization of Flow Structures and Dynamics via Proper Orthogonal Decomposition

Multiphase pipe flow is investigated using proper orthogonal decomposition for tomographic X-ray data, where holdup, cross-sectional phase distributions and phase interface characteristics within the pipe are obtained. Six cases of stratified and mixed flow with water content of 10%, 30% and 80% are investigated to gain insight into effects of velocity and proportion of water on the flow fields. Dispersed and slug flows are separately analyzed to consider the added interface complexity of the flow fields. These regimes are also highly applicable to industry operational flows. Instantaneous and fluctuating phase fractions of the four flow regime are analyzed and reduced order dynamical descriptions are generated. Stratified flow cases display coherent structures that highlight the liquid-liquid interface location while the mixed flow cases show minimal coherence of the eigenmodes. The dispersed flow displays coherent structures for the first few modes near the horizontal center of the pipe, representing the liquid-liquid interface location while the slug flow case shows coherent structures that correspond to the cyclical formation and break up of the slug in the first 5 modes. The low order descriptions of the high water content, stratified flow field indicates that main characteristics can be captured with minimal degrees of freedom. Reconstructions of the dispersed flow and slug flow cases indicate that dominant features are observed in the low order dynamical description utilizing less than 1%of the full order model. POD temporal coefficients a1, a2 and a3 show a high level of interdependence for the slug flowcase. The coefficients also describe the phase fraction holdup as a function of time for both dispersed and slug flow. The second coefficient, a2, and the centerline holdup profile show a mean percent difference below 9% between the two curves. The mathematical description obtained from the decomposition will deepen the understanding of multiphase flow characteristics and is applicable to long distance multiphase transport pipelines, fluidized beds, hydroelectric power and nuclear processes to name a few.

Viggiano, Bianca↗
DOE OSTI2024

Sound Velocities of Stishovite at Simultaneous High Pressure and High Temperature Suggest an Eclogite‐Rich Layer Beneath the Hawaii Hotspot

Compressional and shear wave velocities of polycrystalline stishovite (SiO 2 ) have been measured at simultaneous high pressures and temperatures up to 14.5 GPa and 800°C. By fitting velocities to the finite strain equations, the elastic moduli and density were determined to be K S0 = 306.6(46) GPa, $K$$^{′}_{S}$ = 4.92(10), ∂K S /∂T = −0.024(1) GPa/K, G 0 = 229.0(34) GPa, G′ = 1.07(10), ∂G/∂T = −0.017(1) GPa/K, ρ 0 = 4.287(2) g/cm 3 . Our modeling suggested that, in the eclogite, coesite-stishovite transition can increase P and S wave velocities by 2.4% and 3.5%, respectively. A comparison between geophysical observations and our model shows that the coesite-stishovite phase transition in the eclogite can potentially be responsible for the occurrence of the X discontinuity beneath Hawaii. In addition, our current results suggest an eclogite-rich layer between 340 and 450 km depth beneath Hawaii. The eclogite concentration at the top and bottom of the layer is 41–55 vol% and >77 vol%, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2012

PHOSPHORUS CHEMISTRY IN THE SHOCKED REGION L1157 B1

We study the evolution of phosphorus-bearing species in one-dimensional C-shock models. We find that the abundances of P-bearing species depend sensitively on the elemental abundance of P in the gas phase and on the abundance of N atoms in the pre-shock gas. The observed abundance of PN and the non-detection of PO toward L1157 B1 are reproduced in C-shock models with shock velocity v = 20 km s{sup -1} and pre-shock density n(H{sub 2}) 10{sup 4}-10{sup 5} cm{sup -3}, if the elemental abundance of P in the gas phase is {approx}10{sup -9} and the N-atom abundance is n(N)/n{sub H} {approx}10{sup -5} in the pre-shock gas. We also find that P-chemistry is sensitive to O- and N-chemistry because N atoms are destroyed mainly by OH and NO. We identify the reactions of O-bearing and N-bearing species that significantly affect P-chemistry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Cholesteric–nematic transitions induced by a shear flow and a magnetic field

The untwisting of the helical structure of a cholesteric liquid crystal under the action of a magnetic field and a shear flow has been studied theoretically. Both factors can induce the cholesteric–nematic transition independently; however, the difference in the orienting actions of the magnetic field and the shear flow leads to competition between magnetic and hydrodynamic mechanisms of influence on the cholesteric liquid crystal. We have analyzed different orientations of the magnetic field relative to the direction of the flow in the shear plane. In a number of limiting cases, the analytic dependences are obtained for the pitch of the cholesteric helix deformed by the shear flow. The phase diagrams of the cholesteric–nematic transitions and the pitch of the cholesteric helix are calculated for different values of the magnetic field strength and the angle of orientation, the flow velocity gradient, and the reactive parameter. It is shown that the magnetic field stabilizes the orientation of the director in the shear flow and expands the boundaries of orientability of cholesterics. It has been established that the shear flow shifts the critical magnetic field strength of the transition. It is shown that a sequence of reentrant orientational cholesteric–nematic–cholesteric transitions can be induced by rotating the magnetic field in certain intervals of its strength and shear flow velocity gradients.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Dynamic density field measurements of an explosively driven α → ϵ phase transition in iron

We provide a unique set of observations of the behavior of the α→ϵ phase transition under a complex axially symmetric loading path created by sweeping a detonation wave along the end surface of a cylindrical sample. The primary data sets are the measured mass density distributions acquired at 5 independent times during the sweep of the detonation along the surface. Shocked regions and boundaries are measured, as well as regions and boundaries of elevated density (presumed to be the ϵ−phase iron). The formation and dynamics of these regions were captured and are available for comparisons to material descriptions. We also applied 16 Photon Doppler Velocimetry probes to capture the free surface velocity along a discrete set of radially distributed points in order to compare and correlate the density measurements with previous shock wave studies. The velocimetry data are in nearly exact agreement with previous shock wave studies of the α→ϵ phase transition, the density distributions, while generally in agreement with expectations evolved from the shock wave studies, show that the epsilon phase is generated in regions of high shear stress but at hydrostatic stresses below the typically quoted 13 GPa value. The density field measurements are particularly useful for observing the effects of the forward and reverse transformation kinetics, as well as the reverse transformation hysteresis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

The Earth's velocity for direct detection experiments

The Earth's velocity relative to the Sun in galactic coordinates is required in the rate calculation for direct detection experiments. We provide a rigorous derivation of this quantity to first order in the eccentricity of the Earth's orbit. We also discuss the effect of the precession of the equinoxes, which has hitherto received little explicit discussion. Comparing with other expressions in the literature, we confirm that the expression of Lee, Lisanti and Safdi is correct, while the expression of Lewin and Smith, the de facto standard expression, contains an error. For calculations of the absolute event rate, the leading order expression is sufficient while for modulation searches, an expression with the eccentricity is required for accurate predictions of the modulation phase.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A Herschel/HIFI legacy survey of HF and H{sub 2}O IN the galaxy: Probing diffuse molecular cloud chemistry

We combine Herschel observations for a total of 12 sources to construct the most uniform survey of HF and H{sub 2}O in our Galactic disk. Both molecules are detected in absorption along all sight lines. The high spectral resolution of the Heterodyne Instrument for the Far-infrared (HIFI) allows us to compare the HF and H{sub 2}O distributions in 47 diffuse cloud components sampling the disk. We find that the HF and H{sub 2}O velocity distributions follow each other almost perfectly and establish that HF and H{sub 2}O probe the same gas-phase volume. Our observations corroborate theoretical predictions that HF is a sensitive tracer of H{sub 2} in diffuse clouds, down to molecular fractions of only a few percent. Using HF to trace H{sub 2} in our sample, we find that the N(H{sub 2}O)-to-N(HF) ratio shows a narrow distribution with a median value of 1.51. Our results further suggest that H{sub 2}O might be used as a tracer of H{sub 2}—within a factor of 2.5—in the diffuse interstellar medium (ISM). We show that the measured factor of ∼2.5 variation around the median is driven by true local variations in the H{sub 2}O abundance relative to H{sub 2} throughout the disk. The latter variability allows us to test our theoretical understanding of the chemistry of oxygen-bearing molecules in the diffuse gas. We show that both gas-phase and grain-surface chemistry are required to reproduce our H{sub 2}O observations. This survey thus confirms that grain surface reactions can play a significant role in the chemistry occurring in the diffuse ISM (n{sub H} ⩽ 1000 cm{sup −3}).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Dynamical modeling of tidal streams

I present a new framework for modeling the dynamics of tidal streams. The framework consists of simple models for the initial action-angle distribution of tidal debris, which can be straightforwardly evolved forward in time. Taking advantage of the essentially one-dimensional nature of tidal streams, the transformation to position-velocity coordinates can be linearized and interpolated near a small number of points along the stream, thus allowing for efficient computations of a stream's properties in observable quantities. I illustrate how to calculate the stream's average location (its 'track') in different coordinate systems, how to quickly estimate the dispersion around its track, and how to draw mock stream data. As a generative model, this framework allows one to compute the full probability distribution function and marginalize over or condition it on certain phase-space dimensions as well as convolve it with observational uncertainties. This will be instrumental in proper data analysis of stream data. In addition to providing a computationally efficient practical tool for modeling the dynamics of tidal streams, the action-angle nature of the framework helps elucidate how the observed width of the stream relates to the velocity dispersion or mass of the progenitor, and how the progenitors of 'orphan' streams could be located. The practical usefulness of the proposed framework crucially depends on the ability to calculate action-angle variables for any orbit in any gravitational potential. A novel method for calculating actions, frequencies, and angles in any static potential using a single orbit integration is described in the Appendix.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Magnetic Properties of (Pr{sub 1−x}Sm{sub x}){sub 2}Fe{sub 14}B Melt-Spun Ribbons

The magnetic properties of (Pr{sub 1−x}Sm{sub x}){sub 2}Fe{sub 14}B melt-spun ribbons with varying x from 0 to 1.0 have been studied. The optimum magnetic property of Pr{sub 2}Fe{sub 14}B ribbon could be obtained with the surface velocity of copper wheel of 30 m/s. By melt spinning method, the formation of α-Fe impurity phase in the (Pr{sub 1−x}Sm{sub x}){sub 2}Fe{sub 14}B samples could be inhibited, and the grain could be aligned. Sm element could be incorporated into Pr{sub 2}Fe{sub 14}B main crystalline phase and formed (Pr,Sm){sub 2}Fe{sub 14}B phase with space group P4{sub 2}/mnm. The Curie temperature (T{sub c}) of the (Pr{sub 1−x}Sm{sub x}){sub 2}Fe{sub 14}B ribbons increases gradually from 555 to 617 K with the increase of Sm content x. Both the coercivity (H{sub c}) and remanence ratio (M{sub r}/M{sub s}) decrease firstly and then increase with the Sm content increasing. As long as the Sm content x does not exceed 0.1, the squareness of demagnetization curve remains good. The critical value of transition from easy c-axis to easy plane for (Pr{sub 1−x}Sm{sub x}){sub 2}Fe{sub 14}B ribbons is x= 0.6. When the Sm content x is up to 0.6, all the H{sub c} and M{sub r}/M{sub s} achieved the minimum value, i.e., 243 Oe and 0.346, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Study of solid molecular deuterium D 2 growth under gas pressure

The injection of high-speed cryogenic pellets made of frozen hydrogen-isotopes, represents to date the most effective method to fuel magnetically confined thermonuclear fusion plasmas. Additionally, the injection of very large pellets composed of cryogenic solid of some suitable impurity (typically a noble-gas such as H 2 , Ne, or H 2 /Ne, D 2 /Ne mixtures), shattered in relatively small fragments just before entering the plasma, seems to be the most promising method to reduce the damage risks for the plasma-facing components in case of a plasma disruption. This technology, known as "Shattered Pellet Injection" (SPI), allows to spread out the plasma energy and mitigate possible damage to the in-vessel components, as well as to densify the plasma to suppress the formation of runaway electrons, and/or dissipate their energy. Several techniques to produce and launch cryogenic pellets have been investigated in the past decades. "Pipe gun" injectors are reliable and relatively simple devices are still commonly used today. They make use of single- or two-stage pneumatic light-gas guns to accelerate the pellet at high speeds. In these injectors, the cryogenic pellets are formed “ in situ ” (i.e., inside the launching barrel), by de-sublimating them directly from the gas phase, i.e., at temperatures and pressures below those of the triple point. The simplest case is pure deuterium pellets ($\mathcal{T}$ < 18.7 K, $\mathcal{P}$ < 171.3 hPa). The production of good quality solid deuterium, capable of withstanding the mechanical stress during the acceleration of the pellets, is a key issue. To this end the phase transition of deuterium from gas to solid (and vice versa) is modeled with extensive molecular-dynamics (MD) simulations. Furthermore, the solid growth from the gas phase is simulated in an ample range of temperatures and pressures, to find the best compromise between growth velocity and mechanical properties of the resulting solid system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Temperature-driven massless Kane fermions in HgCdTe crystals

It has recently been shown that electronic states in bulk gapless HgCdTe offer another realization of pseudo-relativistic three-dimensional particles in condensed matter systems. These single valley relativistic states, massless Kane fermions, cannot be described by any other relativistic particles. Furthermore, the HgCdTe band structure can be continuously tailored by modifying cadmium content or temperature. At critical concentration or temperature, the bandgap collapses as the system undergoes a semimetal-to-semiconductor topological phase transition between the inverted and normal alignments. Here, using far-infrared magneto-spectroscopy we explore the continuous evolution of band structure of bulk HgCdTe as temperature is tuned across the topological phase transition. We demonstrate that the rest mass of Kane fermions changes sign at critical temperature, whereas their velocity remains constant. The velocity universal value of (1.07±0.05) × 106 m s -1 remains valid in a broad range of temperatures and Cd concentrations, indicating a striking universality of the pseudo-relativistic description of the Kane fermions in HgCdTe.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantitative spatiotemporal density evolution of aluminum heated purely by monochromatic electrons

A spatially resolved air-wedge shearing interferometer and shadowgraph diagnostic provides measurements of electron density with a resolution of ∼40 μm. A ∼100-ns-long, monoenergetic electron bunch at 19.8 MeV and a current of 1.4 kA (8.5×1014 e−) heats 100-μm-thick aluminum (Al) foils in a 1-mm-spot to Te∼1 eV. A 5-ns-long, ∼60 mJ, frequency doubled Nd:YAG laser probes the dense Al plasma. Electron densities up to 1020cm−3 are resolved; the maximum resolvable density is limited by opacity, transmission, and spatial fringe density achievable with the detector. This diagnostic provides measurements of the total phase shift, transmission, and electron density. Several measurements at different time slices provide the ability to determine the velocity of the leading edge of the shadowgraph and compare it to the motion of different density shells. These measurements are also compared to radiation hydrodynamics simulations. A rough quantitative agreement is shown between the hydro simulations and the measurements; there are differences in the exact density distributions.

Coleman, J. E. (ORCID:0000000308873757)↗
DOE OSTI2015

RADIATION FROM ELECTRON PHASE SPACE HOLES AS A POSSIBLE SOURCE OF JOVIAN S-BURSTS

Radio-frequency short burst emissions (10–40 MHz), known as Jovian S-bursts, have been observed from the Jovian aurora for over fifty years. These emissions, associated with Io’s motion, have a rapidly declining frequency and an exceptionally narrow bandwidth. While it is widely believed that S-bursts are generated by the electron cyclotron maser instability, the mechanism responsible for the rapidly declining frequency and narrow bandwidth currently is not well established. We explore a hypothesis that electron phase space holes radiate or stimulate radiation in the Jovian aurora plasma environment as a possible source of S-burst emissions. Electron phase-space holes (EHs) are ubiquitous in an auroral environment and travel at the implied speeds (∼20,000 km s{sup −1}) of the structures creating the Jovian S-bursts. Furthermore, EHs have the proper physical size to create the observed bandwidth, have sufficient energy content, and can create an environment whereby X mode emissions can be excited. If verified, these findings imply that EHs may be an important source of radiation from strongly magnetized or relativistic astrophysical plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

RESONANT ABSORPTION OF TRANSVERSE OSCILLATIONS AND ASSOCIATED HEATING IN A SOLAR PROMINENCE. I. OBSERVATIONAL ASPECTS

Transverse magnetohydrodynamic waves have been shown to be ubiquitous in the solar atmosphere and can, in principle, carry sufficient energy to generate and maintain the Sun’s million-degree outer atmosphere or corona. However, direct evidence of the dissipation process of these waves and subsequent heating has not yet been directly observed. Here we report on high spatial, temporal, and spectral resolution observations of a solar prominence that show a compelling signature of so-called resonant absorption, a long hypothesized mechanism to efficiently convert and dissipate transverse wave energy into heat. Aside from coherence in the transverse direction, our observations show telltale phase differences around 180° between transverse motions in the plane-of-sky and line-of-sight velocities of the oscillating fine structures or threads, and also suggest significant heating from chromospheric to higher temperatures. Comparison with advanced numerical simulations support a scenario in which transverse oscillations trigger a Kelvin–Helmholtz instability (KHI) at the boundaries of oscillating threads via resonant absorption. This instability leads to numerous thin current sheets in which wave energy is dissipated and plasma is heated. Our results provide direct evidence for wave-related heating in action, one of the candidate coronal heating mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

KINETIC ALFVÉN WAVE GENERATION BY LARGE-SCALE PHASE MIXING

One view of the solar wind turbulence is that the observed highly anisotropic fluctuations at spatial scales near the proton inertial length d{sub p} may be considered as kinetic Alfvén waves (KAWs). In the present paper, we show how phase mixing of large-scale parallel-propagating Alfvén waves is an efficient mechanism for the production of KAWs at wavelengths close to d{sub p} and at a large propagation angle with respect to the magnetic field. Magnetohydrodynamic (MHD), Hall magnetohydrodynamic (HMHD), and hybrid Vlasov–Maxwell (HVM) simulations modeling the propagation of Alfvén waves in inhomogeneous plasmas are performed. In the linear regime, the role of dispersive effects is singled out by comparing MHD and HMHD results. Fluctuations produced by phase mixing are identified as KAWs through a comparison of polarization of magnetic fluctuations and wave-group velocity with analytical linear predictions. In the nonlinear regime, a comparison of HMHD and HVM simulations allows us to point out the role of kinetic effects in shaping the proton-distribution function. We observe the generation of temperature anisotropy with respect to the local magnetic field and the production of field-aligned beams. The regions where the proton-distribution function highly departs from thermal equilibrium are located inside the shear layers, where the KAWs are excited, this suggesting that the distortions of the proton distribution are driven by a resonant interaction of protons with KAW fluctuations. Our results are relevant in configurations where magnetic-field inhomogeneities are present, as, for example, in the solar corona, where the presence of Alfvén waves has been ascertained.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗