Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Particle Methods”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

When should PIC simulations be applied to atmospheric pressure plasmas? Impact of correlation heating

Abstract Molecular dynamics simulations are used to test when the particle-in-cell (PIC) method applies to atmospheric pressure plasmas. It is found that PIC applies only when the plasma density and macroparticle weight are sufficiently small because of two effects associated with correlation heating. The first is the physical effect of disorder-induced heating (DIH). This occurs if the plasma density is large enough that a species (typically ions) is strongly correlated in the sense that the Coulomb coupling parameter exceeds one. In this situation, DIH causes ions to rapidly heat following ionization. PIC is not well suited to capture DIH because doing so requires using a macroparticle weight of one and a grid that well resolves the physical interparticle spacing. These criteria render PIC intractable for macroscale domains. The second effect is a numerical error due to artificial correlation heating (ACH). ACH is like DIH in that it is caused by the Coulomb repulsion between particles, but differs in that it is a numerical effect caused by a macroparticle weight larger than one. Like DIH, it is associated with strong correlations. However, here the macroparticle coupling strength is found to scale as Γ w 2 / 3 , where Γ is the physical coupling strength and w is the macroparticle weight. So even if the physical coupling strength of a species is small, as is expected for electrons in atmospheric pressure plasmas, a sufficiently large macroparticle weight can cause the macroparticles to be strongly coupled and therefore heat due to ACH. Furthermore, it is shown that simulations in reduced dimensions exacerbate these issues.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1987

Fermi acceleration at shocks with arbitrary velocity profiles

Diffusive particle acceleration in modified shock fronts is considered. A method by means of which to obtain the resulting particle spectrum for an arbitrary flow profile is outlined; this method is based on a reformulation of the transport equation as a nonlinear first-order equation. Besides the shape of the flow profile, the spectral index of particles depends on the Peclet number of the shock. Results are given for some sample profiles.

Schneider, Peter↗
NASA NTRSDate not supplied

Development of Challenge Aerosols for Testing Filters in Spacecraft Air Revitalization Systems

The common means for reducing particle concentrations in air in enclosed spaces, including in space habitats, are source prevention and particle removal by air filters. While air filtration and testing is a well-established discipline and industry, testing and classifying filters according to commonly used standards rely on a test aerosol that is often arbitrary and chosen for the convenience of the test method. In space habitats, the particle size distributions are expected to be quite different than the particle size distributions prescribed in test standards, due to the partial or low gravity environment affecting sedimentation of large particulates like hair or cloth fibers, or the introduction of planetary dust to the pressurized volume. This means that the efficacy of the filter will be quite different in the space habitat than specified according to a prevailing standard. This paper will present a means to specify and generate a “composite” test aerosol that is similar to the measured and reported particle sizes utilizing the International Space Station (ISS) as an example. Application of this test aerosol is expected to yield filter efficiencies and loading effects closer to what one can expect on the ISS and to be useful in determining filter lifetime and replacement cycles. The method can be tuned, within reason, to match other particle size distributions one may encounter, such as the intrusion of Lunar dust for a Lunar lander or habitat.

Filtration↗
DOE OSTI2022

Collective plasma effects of electron–positron pairs in beam-driven QED cascades

Understanding the interplay of strong-field QED and collective plasma effects is important for explaining extreme astrophysical environments like magnetars. It has been shown that QED pair plasma can be produced and observed by passing a relativistic electron beam through an intense laser field. This work presents in detail multiple sets of 3D QED-particle-in-cell simulations to show the creation of pair plasma in the QED cascade. The beam driven method enables a high pair particle density and also a low particle Lorentz factor, which both play equal roles on exhibiting large collective plasma effects. Finite laser frequency upshift is observed with both ideal parameters (24 PW laser colliding with a 300 GeV electron beam) and with existing technologies (3 PW laser colliding with a 30 GeV electron beam).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

The fission track record of Apennine Front KREEP basalts

Whitlockite grains in two Apennine Front KREEP soil fragments contain extremely high densities of fission tracks. A new method is described in which particle tracks are examined in an actinide-poor phase bordering an actinide-rich whitlockite phase. This method clearly distinguishes the fission track contributions from other track sources; e.g., Fe-group cosmic rays, spallation recoils, and dislocations. Observed track excesses, when corrected for U- and Th-related fission sources, are 4-20 times greater than the contribution from spontaneous fission of U-238 (based on an age of 4 b.y.), and probably represent a large and variable contribution from the fission of Pu-244. These fragments may have pre-Imbrium ages, underscoring the importance of KREEP as a constituent of the pre-mare lunar crust. Track excesses are not correlated with the U or Th contents of the whitlockite grains. This behavior suggests that Pu fractionates differently from U and Th in lunar igneous KREEP.

Haines, E. L.↗
NASA NTRS1991

Summary of completed project

The primary goal of the work in this proposal has been the determination of the requirement for acceleration of particles, especially electrons, during solar flares. This is one of the important and yet unsolved problems in solar flares and other astronomical objects. This is done by detailed theoretical investigation of interaction of accelerated particles with ambient plasma including interaction with particles, waves, and large scale magnetic fields. Comparison with observed impulsive phase radiation such as X-rays, gamma-rays, microwaves, and type III radiation is used for constraining the models and for determination of the characteristics of the accelerated particles. Steady state and time-dependent numerical codes have been developed based on the Fokker-Planck method for study of the particle acceleration and transport during a flare. Using these, the correlation between characteristics of radiation at various wavelengths have been investigated and some constraints have been set on the model parameters. During the last three year period of this grant, dissertation research of three graduate students have been partially supported by this grant. Each of these theses has dealt with various observed properties of flare emission with particular emphasis on the interrelationship between radiation in at least two energy bands: X-rays and gamma-rays, X-rays and microwaves, and X-rays and type III radiation. This has allowed considerable tightening of constraints on models, and inclusion of wave particle interaction has allowed direct calculations of acceleration processes.

Petrosian, Vahe↗
DOE OSTI2022

Inorganic ternary halide semiconductors for hard radiation detection

Methods and devices for detecting incident radiation, such as incident X-rays, gamma-rays, and/or alpha particle radiation are provided. The methods and devices use high purity, high quality single-crystals of inorganic semiconductor compounds, including solid solutions, having the formula AB2X5, where A represents Tl or In, B represents Sn or Pb, and X represents Br or I, as photoelectric materials.

Kanatzidis, Mercouri G.↗
DOE OSTI2023

Phase diagram of twisted bilayer graphene at filling factor $v = ±3$

Here, we study the correlated insulating phases of twisted bilayer graphene (TBG) in the absence of lattice strain at integer filling $v = ±3$. Using the self-consistent Hartree-Fock method on a particle-hole symmetric model and allowing translation symmetry breaking terms, we obtain the phase diagram with respect to the ratio of AA interlayer hopping (w 0 ) and AB interlayer hopping (w 1 ). When the interlayer hopping ratio is close to the chiral limit (w 0 /w 1 ≲ 0.5), a quantum anomalous Hall state with Chern number v c = ±1 can be observed consistent with previous studies. Around the realistic value w 0 /w 1 ≈ 0.8, we find a spin and valley polarized, translation symmetry breaking, state with C 2⁢z T symmetry, a charge gap and a doubling of the moiré unit cell, dubbed theC 2⁢z T stripe phase. The real-space total charge distribution of this C 2⁢z T stripe phase in the flat band limit does not have modulation between different moiré unit cells, although the charge density in each layer is modulated, and the translation symmetry is strongly broken. Other symmetries, including C 2⁢z , C 2⁢x , particle-hole symmetry P, and the topology of the C 2⁢z T stripe phase, are also discussed in detail. We observed braiding and annihilation of the Dirac nodes by continuously turning on the order parameter to its fully self-consistent value, and provide a detailed explanation of the mechanism for the charge gap opening despite preserving C 2⁢z T and valley U⁡(1) symmetries. In the transition region between the quantum anomalous Hall phase and the C 2⁢z T stripe phase, we find an additional competing state with comparable energy corresponding to a phase with a tripling of the moiré unit cell.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Structure of the Human BK Ion Channel in Lipid Environment

Voltage-gated and ligand-modulated ion channels play critical roles in excitable cells. To understand the interplay among voltage sensing, ligand binding, and channel opening, the structures of ion channels in various functional states and in lipid membrane environments need to be determined. Here, the random spherically constrained (RSC) single-particle cryo-EM method was employed to study human large conductance voltage- and calcium-activated potassium (hBK or hSlo1) channels reconstituted into liposomes. The hBK structure was determined at 3.5 Å resolution in the absence of Ca 2+ . Instead of the common fourfold symmetry observed in ligand-modulated ion channels, a twofold symmetry was observed in hBK in liposomes. Compared with the structure of isolated hSlo1 Ca 2+ sensing gating rings, two opposing subunits in hBK unfurled, resulting in a wider opening towards the transmembrane region of hBK. In the pore gate domain, two opposing subunits also moved downwards relative to the two other subunits.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Design of a Combined Proportional Integral Derivative Controller to Regulate the Temperature Inside a High-Temperature Tubular Solar Reactor

Solar fuels are proven to be promising candidates for thermochemical energy storage. However, the transient nature of solar radiation is an obstacle to maintaining a stable operational temperature inside a solar reactor. To overcome this challenge, the temperature of a solar reactor can be regulated by controlling the incoming solar radiation or the feedstock flowrate inside the reactor. In this work, a combined proportional integral derivative (PID) controller is implemented to regulate the temperature inside a high-temperature tubular solar reactor with counter-current flowing gas/particles. The control model incorporates two control systems to regulate incoming solar radiation and gas flow simultaneously. The design of the controller is based on a reduced-order numerical model of a high-temperature tubular solar reactor that is vertically oriented with an upward gas flow and downward particle flow. The reactor receives heat circumferentially through its wall over a finite segment of its length. Formulation of the heat transfer model is presented by applying the energy balance for the reactor tube and considering heat and mass transfer inside. A set of governing differential equations are solved numerically by using the finite volume method to obtain reactor wall, particles, and gas temperatures along the reactor length with various boundary conditions. Simulation results are used to tune the PID controller parameters by utilizing the Ziegler–Nichols tuning method. Both the simulation results and the controller performance are visualized on the LABVIEW platform. The controller is challenged to track different temperature setpoints with different scenarios of transient solar radiation. Here, the performance of the PID controller was compared to experimental results obtained from an industrial PID controller embedded in a 7 kW electric furnace. Results show that the combined PID controller is successful in maintaining a stable temperature inside the reactor by regulating the incoming solar radiation and the flowrate via small steady-state error and reasonable settling time and overshoot.

14 SOLAR ENERGY↗
DOE OSTI2025

Enhanced Collisional Losses from a Magnetic Mirror Using the Lenard-Bernstein Collision Operator

Collisions are crucial in governing particle and energy transport in plasmas confined in a magnetic mirror trap. Modern gyrokinetic codes model transport in magnetic mirrors, but some utilize approximate model collision operators. This study focuses on a Pastukhov-style method of images calculation of particle and energy confinement times using a Lenard-Bernstein model collision operator. Prior work on parallel particle and energy balances used a different Fokker-Planck plasma collision operator. The method must be extended in non-trivial ways to study the Lenard-Bernstein operator. To assess the effectiveness of our approach, we compare our results with a modern finite element solver. Our findings reveal that the particle confinement time scales like a exp( a 2 ) using the Lenard-Bernstein operator, in contrast to the more accurate scaling that the Coulomb collision operator would yield a 2 exp( a 2 ), where a 2 is approximately proportional to the ambipolar potential. We propose that codes solving for collisional losses in magnetic mirrors utilizing the Lenard-Bernstein or Dougherty collision operator scale their collision frequency of any electrostatically confined species. This study illuminates the collision operator’s intricate role in the Pastukhov-style method of images calculation of collisional confinement.

fusion plasma↗
DOE OSTI2025

Muon tagging with flash ADC waveform baselines

Here, this manuscript describes an innovative method to tag muons using the baseline information of the Flash ADC (FADC) waveform of PMTs in the JSNS 2 (J-PARC Sterile Neutrino Search at J-PARC Spallation Neutron Source) experiment. The experiment is designed to search for evidence of sterile neutrinos, and a reliable method for muon tagging is an essential component for background rejection because the detector is located above ground, on the 3rd floor of the J-PARC Material and Life Science Experimental Facility (MLF). Cosmogenic muons that stop within the detector volume and produce a Michel electron are a particularly important background that must be rejected for our sterile neutrino search. Utilizing this innovative method, more than 99.8 % of Michel electrons can be rejected even without using information from the detector’s veto region PMTs. This technique can be employed by any experiments which uses a similar detector configuration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

New Probe of Cosmic Birefringence Using Galaxy Polarization and Shapes

We propose a novel statistical method to measure cosmic birefringence and demonstrate its power in probing parity violation due to axions. Exploiting an empirical correlation between the integrated radio polarization direction of a spiral galaxy and its apparent shape, we devise an unbiased minimum-variance estimator for the rotation angle, which should achieve an uncertainty of 5°–15° per galaxy. In conclusion, large galaxy samples from the forthcoming SKA continuum surveys, together with optical shape catalogs, promise a comparable or even lower noise power spectrum for the rotation angle than in the CMB Stage-IV (CMB-S4) experiment, with different systematics.

Axion-like particles↗
DOE OSTI2026

Reconstruction of 6D Phase Space

A ‘Phase Space’ is a mathematical representation of all possible physical states of matter of a physical volume or system. Particle beams at Fermilab exist in a six-dimensional phase space: three position dimensions, and three momentum dimensions. To ‘reconstruct’ that space means to take measurement data from virtual detectors along those beams and map out where particles are and how they’re moving, using computational methods. In this experiment, we simulated a particle beam (Fig 1) and used the position distribution of each particle at each detector to map out the magnitude of the beam’s initial momentum.

Shirlee, Jermaine [DuPage Coll.]↗
NASA NTRS1974

The measurement of the size distribution of artificial fogs

The size-distribution of the fog droplets at various fog particle concentrations in fog chamber was determined by two methods: (1) the Stokes' velocity photographic method and (2) using the active scattering particle spectrometer. It is shown that the two techniques are accurate in two different ranges of particle size - the former in the radii range (0.1 micrometers to 10.0 micrometers), and the latter for radii greater than 10.0 micrometers. This was particularly true for high particle concentration, low visibility fogs.

Deepak, A.↗
NASA NTRS2005

Flexible Multi-Shock Shield

Flexible multi-shock shield system and method are disclosed for defending against hypervelocity particles. The flexible multi-shock shield system and method may include a number of flexible bumpers or shield layers spaced apart by one or more resilient support layers, all of which may be encapsulated in a protective cover. Fasteners associated with the protective cover allow the flexible multi-shock shield to be secured to the surface of a structure to be protected.

Christiansen, Eric L.↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Solid state electrolytes

A product includes a solid state electrolyte particle coated with a coating consisting essentially of metal oxide. A method includes depositing a coating on a solid state electrolyte particle by atomic layer deposition. A method includes fabricating a product using a plurality of solid state electrolyte particles coated with a coating consisting essentially of metal oxide.

Ye, Jianchao↗