Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle simulations”

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

What Controls Crystal Diversity and Microphysical Variability in Cirrus Clouds?

Abstract Variability of ice microphysical properties like crystal size and density in cirrus clouds is important for climate through its impact on radiative forcing, but challenging to represent in models. For the first time, recent laboratory experiments of particle growth (tied to crystal morphology via deposition density) are combined with a state‐of‐the‐art Lagrangian particle‐based microphysics model in large‐eddy simulations to examine sources of microphysical variability in cirrus. Simulated particle size distributions compare well against balloon‐borne observations. Overall, microphysical variability is dominated by variability in the particles' thermodynamic histories. However, diversity in crystal morphology notably increases spatial variability of mean particle size and density, especially at mid‐levels in the cloud. Little correlation between instantaneous crystal properties and supersaturation occurs even though the modeled particle morphology is directly tied to supersaturation based on laboratory measurements. Thus, the individual thermodynamic paths of each particle, not the instantaneous conditions, control the evolution of particle properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Anisotropic interactions for continuum modeling of protein–membrane systems

In this work, a model for anisotropic interactions between proteins and cellular membranes is proposed for large-scale continuum simulations. The framework of the model is based on dynamic density functional theory, which provides a formalism to describe the lipid densities within the membrane as continuum fields while still maintaining the fidelity of the underlying molecular interactions. Within this framework, we extend recent results to include the anisotropic effects of protein–lipid interactions. As applications, we consider two membrane proteins of biological interest: a RAS–RAF complex tethered to the membrane and a membrane embedded G protein-coupled receptor. A strong qualitative and quantitative agreement is found between the numerical results and the corresponding molecular dynamics simulations. Combining the scope of continuum level simulations with the details from molecular level particle simulations enables research into protein–membrane behaviors at a more biologically relevant scale, which crucially can also be accessed via experiment.

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

Generation models of electron conics

Electron distribution functions (EDFs) with a peak oblique to the magnetic field, adjacent to but distinct from loss-cone features, have been observed by the DE 1, Viking, and S3-3 spacecraft in passes through the nightside auroral zone, polar cap, dayside cusp, and extended dayside auroral oval. Using particle simulations, two types of wave excitation and particle acceleration mechanisms which may contribute to producing these electron conic distributions are investigated. The first involves excitation of upper-hybrid waves by the electron loss cone, and the subsequent perpendicular heating of the background and thermal electrons. The second involves excitation of downward propagating parallel modes by an auroral electron beam which frequently accompanies the upflowing electron conics. These modes provide parallel acceleration, which modifies the GDF. Those electrons which are not lost to the atmosphere and mirror back up the magnetic field line give rise to enhancements in the GDF at the edge of the loss cone.

Roth, I.↗
NASA NTRS1993

An overview of the cosmic dust analogue material production in reduced gravity: the STARDUST experience

The formation, properties and chemical dynamics of microparticles are important in a wide variety of technical and scientific fields including synthesis of semiconductor crystals from the vapour, heterogeneous chemistry in the stratosphere and the formation of cosmic dust surrounding the stars. Gravitational effects on particle formation from vapors include gas convection and buoyancy and particle sedimentation. These processes can be significantly reduced by studying condensation and agglomeration of particles in microgravity. In addition, to accurately simulate particle formation near stars, which takes place under low gravity conditions, studies in microgravity are desired. We report here the STARDUST experience, a recent collaborative effort that brings together a successful American program of microgravity experiments on particle formation aboard NASA KC-135 Reduced Gravity Research Aircraft and several Italian research groups with expertise in microgravity research and astrophysical dust formation. The program goal is to study the formation and properties of high temperature particles and gases that are of interest in astrophysics and planetary science. To do so we are developing techniques that are generally applicable to study particle formation and properties, taking advantage of the microgravity environment to allow accurate control of system parameters.

NASA Discipline Exobiology↗
NASA NTRS2016

Probing Acceleration and Turbulence at Relativistic Shocks in Blazar Jets

Diffusive shock acceleration (DSA) at relativistic shocks is widely thought to be an important acceleration mechanism in various astrophysical jet sources, including radio-loud active galactic nuclei such as blazars. Such acceleration can produce the non-thermal particles that emit the broad-band continuum radiation that is detected from extragalactic jets. An important recent development for blazar science is the ability of Fermi-Large Area Telescope spectroscopy to pin down the shape of the distribution of the underlying non-thermal particle population. This paper highlights how multiwavelength spectra spanning optical to X-ray to gamma-ray bands can be used to probe diffusive acceleration in relativistic, oblique, magnetohydrodynamic (MHD) shocks in blazar jets. Diagnostics on the MHD turbulence near such shocks are obtained using thermal and non-thermal particle distributions resulting from detailed Monte Carlo simulations of DSA. These probes are afforded by the characteristic property that the synchrotron vF(sub v), peak energy does not appear in the gamma-ray band above 100 MeV. We investigate self-consistently the radiative synchrotron and inverse Compton signatures of the simulated particle distributions. Important constraints on the diffusive mean free paths of electrons, and the level of electromagnetic field turbulence are identified for three different case study blazars, Mrk 501, BL Lacertae and AO 0235+164. The X-ray excess of AO 0235+164 in a flare state can be modelled as the signature of bulk Compton scattering of external radiation fields, thereby tightly constraining the energy-dependence of the diffusion coefficient for electrons. The concomitant interpretations that turbulence levels decline with remoteness from jet shocks, and the probable significant role for non-gyroresonant diffusion, are posited.

Baring, Matthew G.↗
DOE OSTI2024

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1991

Simulations of pulsed electron beam injection during active experiments

Two-dimensional electromagnetic particle simulations are used to investigate the characteristics of the spacecraft charging, particle acceleration, and wave emissions during beam injection, turnoff and subsequent pulsing of the beam. It is shown that, during beam injection, the beam current is neutralized by a spatially separate return current region extending several tens of meters from the beam region, with the currents being closed across the field lines by the perpendicular acceleration of ambient plasma ions into the beam region. After beam turnoff, this current system reverses after a time lag of about an ion plasma period. The current reversal is accompanied by prolonged electron collection by the spacecraft near the beam region, preferential ion collection by sections of the spacecraft magnetically connected to the initial return current regions, and the creation of hot plasma extending well into the return current regions. Because of the time lag, the currents induced in the plasma during periodic beam injection can be strongly modified from the imposed beam current.

Winglee, R. M.↗
NASA NTRS2000

Pulmonary Toxicity Study of Lunar and Martian Dust Simulants Intratracheally Instilled in Mice

NASA is contemplating sending humans to Mars and the Moon for further exploration. The properties of Hawaiian and Californian volcanic ashes allow them to be used to simulate Martian and lunar dusts, respectively. NASA laboratories use these dust simulants to test performance of hardware destined for Martian or lunar environments. Workers in these test facilities are exposed to low levels of these dusts. The present study was conducted to investigate the toxicity of these dust simulants. Particles of respirable-size ranges of lunar simulant (LS), Martian simulant (MS), TiO2 (negative control) and quartz (positive control) were each intratracheally instilled (saline as vehicle) to groups of 4 mice (C57BL, male, 2-3 month old) at a single treatment of 1 (Hi dose) or 0.1 (Lo dose) mg/mouse. The lungs were harvested at the end of 7 days or 90 days for histopathological examination. Lungs of the LS-Lo groups had no evidence of inflammation, edema or fibrosis. The LS-Hi-7d group had mild to moderate acute inflammation, and neutrophilic and lymphocytic infiltration; the LS-Hi-90d group showed signs of chronic inflammation and some fibrosis. Lungs of the MS-Lo-7d group revealed mild inflammation and neutrophilic and lymphocytic infiltration; the MS-Lo-90d group showed mild fibrosis and particle-laden macrophages (PLM). Lungs of the MS-Hi-7d group demonstrated mild to moderate inflammation and large foci of PLM; the MS-Hi-90d group showed chronic mild to moderate inflammation and fibrosis. To mimic the effects of the oxidative and reactive properties of Martian soil surface, groups of mice were exposed to ozone (3 hour at 0.5 ppm) prior to MS dust instillation. Lung lesions in the MS group were more severe with the pretreatment. The results for the negative and positive controls were consistent with the known pulmonary toxicity of these compounds. The overall severity of toxic insults to the lungs were TiO2<LS<MS<Quartz. For the mice in the 90-d study, blood samples were taken for immunotoxicity study. Antinuclear antibodies (such as those against Golgi, kinetichore, and centromere) were detected in serum in a greater extent in mice treated with LS than those with quartz, suggesting that LS, like quartz, has the potential of inducing antoimmune disease.

Lam, Chiu-Wing↗
DOE OSTI2021

Simulation of Electron-Proton Scattering Events by a Feature-Augmented and Transformed Generative Adversarial Network (FAT-GAN)

We apply generative adversarial network (GAN) technology to build an event generator that simulates particle production in electron-proton scattering that is free of theoretical assumptions about underlying particle dynamics. The difficulty of efficiently training a GAN event simulator lies in learning the complicated patterns of the distributions of the particles physical properties. We develop a GAN that selects a set of transformed features from particle momenta that can be generated easily by the generator, and uses these to produce a set of augmented features that improve the sensitivity of the discriminator. The new Feature-Augmented and Transformed GAN (FAT-GAN) is able to faithfully reproduce the distribution of final state electron momenta in inclusive electron scattering, without the need for input derived from domain-based theoretical assumptions. The developed technology can play a significant role in boosting the science of existing and future accelerator facilities, such as the Electron-Ion Collider.

Alanazi, Yasir↗
NASA NTRS1992

Numerical simulations of bursty radio emissions from planetary magnetospheres

One-dimensional electromagnetic particle simulations are used to investigate the characteristics of radiation from electron beams with a large temperature anisotropy, in order to identify the origin of the smooth and bursty radio emissions from Uranus and Neptune observed by the Voyager spacecraft. It is shown that these electron beams, which would typically originate from the sporadic or impulsive injections of energetic electrons, can generate electromagnetic radiation which should be able to escape into the solar wind despite the growth of the electrostatic instability. The amount of radiation with frequencies above the local x mode cutoff increases with the beam speed. It is also proposed that some of the radiation generated below the local x mode cutoff may also be able to escape the plasma and be detected remotely via mode conversion between regions where field-aligned currents produce local perturbations in the magnetic field.

Winglee, R. M.↗
NASA NTRS1993

Three-dimensional simulations of the solar wind interaction with Mars

Three-dimensional hybrid particle simulations of the dayside portion of the Mars/solar wind interaction have been performed. These simulations are compared with the in situ measurements taken during the elliptical orbits of Phobos 2. The comparisons show considerable agreement between the magnetic field data and the simulations results. The results of the simulations bring into question the type of structures created in the subsolar region of Mars. It appears that the large larmor radius of the solar wind ions prevents the formation of a traditional collisionless shock in the subsolar region of the interaction.

Brecht, Stephen H.↗
NASA NTRS2022

Asymmetrically Varying Guide Field During Magnetic Reconnection: Particle-In-Cell Simulations

Using fully kinetic particle-in-cell modeling, we investigate how magnetic reconnection responds to a varying guide field in one of the inflow regions. We find that the reconnection rate varies significantly when the orientation of the magnetic field changes between being strictly antiparallel and having a guide field. These variations are fairly consistent with the scaling relation for asymmetric reconnection developed by Cassak and Shay (2007). However, the rate is also found to be nonlinearly modulated by changes in the ion inflow velocity. The spatio-temporal change in the inflow velocity arises as the magnetic forces reconfigure to regions of different magnetic field strengths. The variations in the inflow magnetic field configuration allow for different gradients in the magnetic field, leading to asymmetries in the magnetic tension force. By momentum conservation, this facilitates asymmetries in the inflow velocity, which in turn affects the flux transport into the reconnection site. The outflow is found to be less laminar when the inflow varies, and various signatures of the inflow variations are identified in the outflow.

Susanne F Spinnangr↗
DOE OSTI2021

Deterministic scattering of relativistic electron beams by off-resonant circularly polarized electromagnetic waves

A mechanism is presented whereby relativistic electron beams localized in phase space are deterministically scattered by coherent circularly polarized electromagnetic waves without stochastic processes. It is shown via an exact single-particle analysis that the condition for maximal scattering is an off-resonant condition, contrary to previous kinetic analyses that predict maximal diffusion or interaction at exact resonance or its harmonics. The mechanism, verified by single-particle simulations, enables a fast, nonlinear redistribution of the beam particles. A possible application of this mechanism to runaway electron suppression is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Validation of the NOνA experiment 2023-tuning on simulated neutrino-matter interactions

NOνA is a long-baseline neutrino oscillation experiment that utilizes a two-detector design to study the oscillations of muon neutrinos into electron neutrinos over a baseline of 810 km. The Near Detector (ND) measures the neutrino beam spectrum and composition before oscillation, which is then compared to the oscillated neutrino energy spectrum observed in the Far Detector (FD). In the ND, the neutrinos are detected through their interactions with the heavy target nuclei within the detector. NOνA employs the GENIE neutrino event generator for simulating these neutrino-nucleus interactions. However, the default GENIE prediction does not adequately reproduce the ND data. To address this, NOνA developed a tune of the neutrino interaction models within GENIE version 3:0:6 to minimize discrepancies between the simulated predictions and the observed data in the ND. This dissertation tests the NOνA’s 2023 tune of the GENIE neutrino cross-section simulations by performing a data/simulations comparison for the ND. The analysis employed datasets comprising $2.55\times10^{21}$ protons-on-target (POT) in neutrino beam mode and $1.14\times10^{21}$ POT in antineutrino beam mode. The NOνA tuning of neutrino-matter interaction simulations matches with ND data within the $1\sigma$ error band, overestimating muon neutrino and antineutrino charged current interactions by approximately 6 % and 9 %, respectively. Discrepancies were observed in the energy region dominated by Quasi-Elastic-like interactions. Systematic uncertainties associated with the modeling of the neutrino cross-section, particularly those pertaining Quasi-Elastic like interactions, contributed considerably to the overall error in the simulations. Furthermore, the reconstruction algorithm used in NOνA for particle classification demonstrated significant misidentifications between charged pions and protons, as well as a tendency to overlook additional pions or protons in multi-particle simulated events.

Cortés Parra, Camilo Andrés↗
NASA NTRS1992

Two dimensional PIC simulations of plasma heating by the dissipation of Alfven waves

Two dimensional plasma particle simulations of the evolution of large amplitude circularly polarized Alfven waves propagating parallel to the magnetic field show that the waves decay via both one- and two- dimensional parametric decay instabilities. For parameters studied, one-dimensional processes dominate the simulations, but two-dimensional decay processes, including the recently predicted filamentation instability are also observed. The daughter waves generated by the parametric decay are primarily damped by the ions, leading to ion heating. The parametric decay processes efficiently convert the ordered fluid ion motion in the Alfven wave into ion thermal energy. These processes may be important for the dissipation of Alfven waves in the solar wind, the corona and other space plasma environments. The computations were performed on the Intel Touchstone parallel supercomputer.

Liewer, P. C.↗
NASA NTRS1994

Electron collection by a highly positive satellite in the ionosphere: Test aprticle simulation

Collection of electrons by a satellite at a relatively high potential in low Earth orbit, including the effects of the satellite's orbital motion, remains a poorly understood phenomenon. Using a test particle simulation in which charged particle motion is tracked in prescribed electric fields, we calculate here the current collection and the nature of the energy distribution function of the electrons collected by the satellite, including the effects of the satellite orbital motion. Calculations of the collected current without the orbital motion show an excellent agreement with the current from the Parker-Murphy model, but with the orbital motion of the satellite the current collection is enhanced, the degree of enhancement depending on the size of the satellite sheath extending along the magnetic field line. In the latter case, the flow of electrons around the satellite shows some interesting behaviors including azimuthal E x B drift around the satellite, axial trapping along the magnetic field, and formation of field-aligned flow of electrons in the wake region. The total energy of the collected electrons is approximately equal e phi(sub 0), where phi(sub 0) is the satellite potential, but the partition of the energy into components parallel (W(sub parallel)) and perpendicular (W(sub perpendicular)) to the magnetic field shows interesting features. For the magnetic field along the polar axis Z, the energy distributions near the poles are perfectly field aligned, that is, the parallel energy (W(sub parallel) = e phi(sub 0). The perpendicular component (W(sub perpendicular)) progressively increases toward the 'equator' of the spherical satellite. For a sufficiently large bias potential phi(sub 0), the equatorial energy distributions f(W(sub parallel)) and f(W(sub perpendicular)) show fine structures having multiple peaks. These features of the distribution functions result from the chaotic motion of electrons and the associated coupling between the electron motions parallel and perpendicular to the magnetic field in strongly nonuniform electric fields near the satellite.

Singh, Nagendra↗
DOE OSTI2020

Testing dark matter halo properties using self-similarity

ABSTRACT We use self-similarity in N-body simulations of scale-free models to test for resolution dependence in the mass function and two-point correlation functions of dark matter haloes. We use 10243 particle simulations performed with abacus , and compare results obtained with two halo finders: friends-of-friends (fof ), and rockstar . The fof mass functions show a systematic deviation from self-similarity which is explained by resolution dependence of the fof mass assignment previously reported in the literature. Weak evidence for convergence is observed only starting from haloes of several thousand particles, and mass functions are overestimated by at least as much as $20-25{{\ \rm per\ cent}}$ for haloes of 50 particles. The mass function of the default rockstar halo catalogue (with bound virial spherical overdensity mass), on the other hand, shows good convergence of the order of 50 to 100 particles per halo, with no detectable evidence at the few percent level of any systematic dependence for larger particle number. Tests show that the mass unbinding procedure in rockstar is the key factor in obtaining this much improved resolution. Applying the same analysis to the halo–halo two point correlation function, we find again strong evidence for convergence only for rockstar haloes, at separations sufficiently large so that haloes do not overlap. At these separations, we can exclude dependence on resolution at the $5-10{{\ \rm per\ cent}}$ level once haloes have of the order of 50 to 100 particles. At smaller separations results are not converged even at significantly larger particle number, and bigger simulations would be required to establish the resolution required for convergence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Advancing set-conditional set generation: Diffusion models for fast simulation of reconstructed particles

The computational intensity of detector simulation and event reconstruction poses a significant difficulty for data analysis in collider experiments. This challenge inspires the continued development of machine learning techniques to serve as efficient surrogate models. We propose a fast emulation approach that combines simulation and reconstruction. In other words, a neural network generates a set of reconstructed objects conditioned on input particle sets. To make this possible, we advance set-conditional set generation with diffusion models. Using a realistic, generic, and public detector simulation and reconstruction package (COCOA), we show how diffusion models can accurately model the complex spectrum of reconstructed particles inside jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗