Search NASA⌕ Search

SEARCH · Search NASA

Results for “particle dynamics simulation”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Isokinetic TWC Evaporator Probe: Calculations and Systemic Uncertainty Analysis

A new Isokinetic Total Water Content Evaporator (IKP2) was downsized from a prototype instrument, specifically to make airborne measurements of hydrometeor total water content (TWC) in deep tropical convective clouds to assess the new ice crystal Appendix D icing envelope. The probe underwent numerous laboratory and wind tunnel investigations to ensure reliable operation under the difficult high altitude/speed/TWC conditions under which other TWC instruments have been known to either fail, or have unknown performance characteristics and the results are presented in a companion paper. This paper presents the equations used to determine the total water content (TWC) of the sampled atmosphere from the values measured by the IKP2 or necessary ancillary data from other instruments. The uncertainty in the final TWC is determined by propagating the uncertainty in the measured values through the calculations to the final result. Two techniques were used and the results compared. The first is a typical analytical method of propagating uncertainty and the second performs a Monte Carlo simulation. The results are very similar with differences that are insignificant for practical purposes. The uncertainty is between 2 percent and 3 percent at most practical operating conditions. The capture efficiency of the IKP2 was also examined based on a computational fluid dynamic simulation of the original IKP and scaled down to the IKP2. Particles above 24 microns were found to have a capture efficiency greater than 99 percent at all operating conditions.

Aircraft Icing↗
NASA NTRS1985

Acoustic waves in the solar atmosphere. VII - Non-grey, non-LTE H(-) models

The propagation and shock formation of radiatively damped acoustic waves in the solar chromosphere are studied under the assumption that H(-) is the only absorber; the opacity is non-grey. Deviations from local thermodynamic equilibrium (LTE) are permitted. The results of numerical simulations show the depth dependence of the heating by the acoustic waves to be insensitive to the mean state of the atmosphere. After the waves have developed into shocks, their energy flux decays exponentially with a constant damping length of about 1.4 times the pressure scale height, independent of initial flux and wave period. Departures from LTE have a strong influence on the mean temperature structure in dynamical chromosphere models; this is even more pronounced in models with reduced particle density - simulating conditions in magnetic flux tubes - which show significantly increased temperatures in response to mechanical heating. When the energy dissipation of the waves is sufficiently large to dissociate most of the H(-) ions, a strong temperature rise is found that is reminiscent of the temperature structure in the transition zone between chromosphere and corona; the energy flux remaining in the waves then drives mass motions.

Schmitz, F.↗
DOE OSTI2021

Multi-fidelity kinetic theory-based approach for the prediction of particle attrition: Application to jet cup attrition system

The timescale difference between the “fast” flow dynamics of fluidized bed reactors and the relatively “slow” rate of particle degradation makes the direct computational prediction of attrition challenging. An approach to this challenge is a multi-fidelity strategy where a high fidelity model for the flow dynamics is coupled with a lower fidelity model for the long-time resolution of the bulk attrition of the reactor inventory. We implement this approach using high-fidelity kinetic theory simulations to calculate the flow dynamics which are post-processed to calculate the frequency and intensity of the particle-particle and particle-wall collisions (e.g. collision energy spectra). This is combined with the particle breakage properties to construct the coefficients for a low-fidelity model [e.g. Monazam et al., 2018, Powder Technology 340, p. 528-536]. Simulations are performed of a jet cup attrition system containing Canadian hematite (Monazam et al. 2018). Furthermore, these are first analyzed using the collision energy spectra. Quantitative predictions of the mass loss are made using a low-fidelity model derived from the collision-spectra and a calibrated material breakage coefficient. The results are found to compare favorably with the experimental measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Charging-driven coarsening and melting of a colloidal nanoparticle monolayer at an ionic liquid–vacuum interface

In this work, we induce and investigate the coarsening and melting dynamics of an initially static nanoparticle colloidal monolayer at an ionic liquid–vacuum interface, driven by a focused, scanning electron beam. Coarsening occurs through grain interface migration and larger-scale motions such as grain rotations, often facilitated by sliding dislocations. The progressive decrease in area fraction that drives melting of the monolayer is explained using an electrowetting model whereby particles at the interface are solvated once their accumulating charge recruits sufficient counterions to subsume the particle. Subject to stochastic particle removal from the monolayer, melting is recapitulated in simulations with a Lennard-Jones potential. This new driving mechanism for colloidal systems, whose dynamical timescales we show can be controlled with the accelerating voltage, opens the possibility to manipulate particle interactions dynamically without need to vary particle intrinsic properties or surface treatments. Furthermore, the decrease in particle size availed by electron imaging presents opportunities to observe force and time scales in a lesser-explored regime intermediate between typical colloidal and molecular systems.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2024

Dynamic population balance in molecular-level simulations of hypersonic flows

This report summarizes the work towards developing stochastic weighted particle methods (SWPM) for future application in hypersonic flows. Extensive changes to Sandia’s direct simulation Monte Carlo (DSMC) solver, SPARTA (Stochastic Particle Real Time Analyzer), were made to enable the necessary particle splitting and reduction capabilities for SWPM. The results from one-dimensional Couette and Fourier flows suggest that SWPM can reproduce the correct transport for a large range of Knudsen numbers with adequate accuracy. The associated velocity and temperature profiles are in good agreement with DSMC. An issue with particle placement during particle number reduction, is identified, to which, a simple but effective solution based on minimizing the center of mass error is proposed. High Mach wheel flows are simulated using the SWPM and DSMC methods. SWPM is capable of providing nearly an order of magnitude increase in efficiency over DSMC while retaining high accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Parametrizations of coupled betatron motion for strongly coupled lattices

The coupling of transverse motion is a natural occurrence in particle accelerators, either in the form of a residual coupling arising from imperfections or originating by design from strong systematic coupling fields. While the first can be treated perturbatively, the latter requires a robust approach adapted to strongly coupled optics, and a parametrization of the linear optics must be performed to explore beam dynamics in such peculiar lattices. This work highlights the key physical interpretations of the main parametrization formalisms to describe linear coupled optics, along with explicit links and comparisons of these parametrizations. Concepts rarely illustrated in other works, such as forced mode flips and local coupling, are explored in detail, clarifying some anomalies that can arise in lattice functions. The analytical methods have been implemented in a reference Python package and connected with ray-tracing and integration codes to explore examples of strongly coupled lattices, which are discussed in detail to highlight the key physical interpretations of the parametrizations and characteristics of the lattices.

47 OTHER INSTRUMENTATION↗
NASA NTRS2013

Representation of Nucleation Mode Microphysics in a Global Aerosol Model with Sectional Microphysics

In models, nucleation mode (1 nm<Dp <10 nm) particle microphysics can be represented explicitly with aerosol microphysical processes or can be parameterized to obtain the growth and survival of nuclei to the model's lower size boundary. This study investigates how the representation of nucleation mode microphysics impacts aerosol number predictions in the TwO-Moment Aerosol Sectional (TOMAS) aerosol microphysics model running with the GISS GCM II-prime by varying its lowest diameter boundary: 1 nm, 3 nm, and 10 nm. The model with the 1 nm boundary simulates the nucleation mode particles with fully resolved microphysical processes, while the model with the 10 nm and 3 nm boundaries uses a nucleation mode dynamics parameterization to account for the growth of nucleated particles to 10 nm and 3 nm, respectively.We also investigate the impact of the time step for aerosol microphysical processes (a 10 min versus a 1 h time step) to aerosol number predictions in the TOMAS models with explicit dynamics for the nucleation mode particles (i.e., 3 nm and 1 nm boundary). The model with the explicit microphysics (i.e., 1 nm boundary) with the 10 min time step is used as a numerical benchmark simulation to estimate biases caused by varying the lower size cutoff and the time step. Different representations of the nucleation mode have a significant effect on the formation rate of particles larger than 10 nm from nucleated particles (J10) and the burdens and lifetimes of ultrafinemode (10 nm≤Dp ≤70 nm) particles but have less impact on the burdens and lifetimes of CCN-sized particles. The models using parameterized microphysics (i.e., 10 nm and 3 nm boundaries) result in higher J10 and shorter coagulation lifetimes of ultrafine-mode particles than the model with explicit dynamics (i.e., 1 nm boundary). The spatial distributions of CN10 (Dp ≥10 nm) and CCN(0.2 %) (i.e., CCN concentrations at 0.2%supersaturation) are moderately affected, especially CN10 predictions above ∼700 hPa where nucleation contributes most strongly to CN10 concentrations. The lowermost-layer CN10 is substantially improved with the 3 nm boundary (compared to 10 nm) in most areas. The overprediction in CN10 with the 3 nm and 10 nm boundaries can be explained by the overprediction of J10 or J3 with the parameterized microphysics, possibly due to the instantaneous growth rate assumption in the survival and growth parameterization. The errors in CN10 predictions are sensitive to the choice of the lower size boundary but not to the choice of the time step applied to the microphysical processes. The spatial distribution of CCN(0.2 %) with the 3 nm boundary is almost identical to that with the 1 nm boundary, but that with the 10 nm boundary can differ more than 10-40% in some areas. We found that the deviation in the 10 nm simulations is partly due to the longer time step (i.e., 1 h time step used in the 10 nm simulations compared to 10 min time step used in the benchmark simulations), but, even with the same time step, the 10 nm cutoff showed noticeably higher errors than the 3 nm cutoff. In conclusion, we generally recommend using a lower diameter boundary of 3 nm for studies focused on aerosol indirect effects but down to 1 nm boundary for studies focused on CN10 predictions or nucleation.

microphysics↗
NASA NTRS1993

Design and testing of a unique randomized gravity, continuous flow bioreactor

A rotating, null gravity simulator, or Couette bioreactor was successfully used for the culture of mammalian cells in a simulated microgravity environment. Two limited studies using Lipomyces starkeyi and Streptomyces clavuligerus were also conducted under conditions of simulated weightlessness. Although these studies with microorganisms showed promising preliminary results, oxygen limitations presented significant limitations in studying the biochemical and cultural characteristics of these cell types. Microbial cell systems such as bacteria and yeast promise significant potential as investigative models to study the effects of microgravity on membrane transport, as well as substrate induction of inactive enzyme systems. Additionally, the smaller size of the microorganisms should further reduce the gravity induced oscillatory particle motion and thereby improve the microgravity simulation on earth. Focus is on the unique conceptual design, and subsequent development of a rotating bioreactor that is compatible with the culture and investigation of microgravity effects on microbial systems. The new reactor design will allow testing of highly aerobic cell types under simulated microgravity conditions. The described reactor affords a mechanism for investigating the long term effects of reduced gravity on cellular respiration, membrane transfer, ion exchange, and substrate conversions. It offers the capability of dynamically altering nutrients, oxygenation, pH, carbon dioxide, and substrate concentration without disturbing the microgravity simulation, or Couette flow, of the reactor. All progeny of the original cell inoculum may be acclimated to the simulated microgravity in the absence of a substrate or nutrient. The reactor has the promise of allowing scientists to probe the long term effects of weightlessness on cell interactions in plants, bacteria, yeast, and fungi. The reactor is designed to have a flow field growth chamber with uniform shear stress, yet transfer high concentrations of oxygen into the culture medium. The system described allows for continuous, on line sampling for production of product without disturbing fluid and particle dynamics in the reaction chamber. It provides for the introduction of substrate, or control substances after cell adaptation to simulated microgravity has been accomplished. The reactor system provides for the nondisruptive, continuous flow replacement of nutrient and removal of product. On line monitoring and control of growth conditions such as pH and nutrient status are provided. A rotating distribution valve allows cessation of growth chamber rotation, thereby preserving the simulated microgravity conditions over longer periods of time.

Lassiter, Carroll B.↗
NASA NTRS2024

Experimental Validation of Charged Lunar Dust Dynamics Simulants

The goal of this work is to experimentally verify the electrostatic interaction physics that have been incorporated into the discrete element method (DEM) modeling software package during the FY20 CIF project titled “Charged Particle Dynamics in the Lunar Environment.” Current state-of-the-art (SOA) granular gas dynamics models used to explain rocket plume impingement physics have not taken into account the natural or induced electrostatic environment of the lunar surface, nor the effect of charged regolith grains being present in the plume plasma. This work improved upon the SOA by adding long-range and contact inter-particle electrostatic interactions to a granular mechanics DEM modeling software package. Tribocharging mechanics between spherical particles has been experimentally examined along with the natural electrostatic plasma environments produced via emission from an electron gun and ultraviolet lamp.

Jay Phillips↗
DOE OSTI2024

Quantum Simulation of the First-Quantized Pauli-Fierz Hamiltonian

We provide an explicit recursive divide-and-conquer approach for simulating quantum dynamics and derive a discrete first-quantized nonrelativistic QED Hamiltonian based on the many-particle Pauli-Fierz Hamiltonian. We apply this recursive divide-and-conquer algorithm to this Hamiltonian and compare it to a concrete simulation algorithm that uses qubitization. Our divide-and-conquer algorithm, using lowest-order Trotterization, scales for fixed grid spacing as O ~ ( Λ N 2 η 2 t 2 / ϵ ) for grid size N , η particles, simulation time t , field cutoff Λ , and error ϵ . Our qubitization algorithm scales as O ~ ( N ( η + N ) ( η + Λ 2 ) t log ( 1 / ϵ ) ) . This shows that even a naive partitioning and low-order splitting formula can yield, through our divide-and-conquer formalism, superior scaling to qubitization for large Λ . We compare the relative costs of these two algorithms on systems that are relevant for applications such as the spontaneous emission of photons and the photoionization of electrons. We observe that for different parameter regimes, one method can be favored over the other. Finally, we give new algorithmic and circuit-level techniques for gate optimization, including a new way of implementing a group of multicontrolled- X gates that can be used for better analysis of circuit cost. Published by the American Physical Society 2024

Mukhopadhyay, Priyanka (ORCID:0000000164639100)↗
DOE OSTI2022

Sensitivity of puff dynamics and airborne droplet nuclei distribution to variations in violent expiration events

Here, we use large eddy simulations to investigate the puff and droplet dynamics from violent expiratory events such as coughs and sneezes in the first few seconds following an ejection. For each of the eleven simulations considered, over 60,000 droplets are ejected and individually tracked using the point-particle Euler–Lagrange approach. We test the sensitivity of the puff and droplet dynamics to various parameters including the ejection volume, momentum, and orientation. We also explore the effect of the mouth shape on the aforementioned dynamics by considering elliptical and circular inlet cross-sections. The results from the simulations compare favorably with a recent theoretical framework put forth by Balachandar et al. in terms of the puff size and propagation velocity. More importantly however, the theory is able to accurately predict the number and size spectra of the potentially virus-laden droplet nuclei that remain airborne within the puff. We observe that the ejection angle and mouth shape do not significantly affect the puff and droplet dynamics. Additionally, we quantify the carrying capacity of the detached puff portions in terms of the number and size spectra of droplets/droplet nuclei suspended within.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Roles of Hydrogen Bonds and Alignment in Oriented Attachment of Gibbsite Nanoparticles: Insights from Molecular Dynamics

Oriented attachment (OA) of nanoparticles is an important crystal growth pathway in the synthesis of hierarchical structures. Although a significant understanding of OA has been made, the effect of atomistic misalignment and the roles of solvent/particle and particle/particle interactions on the structure–energy relationship during an OA remain elusive. In this study, we perform molecular dynamics simulations to calculate the potential of mean force (PMF) profile for gibbsite particle translation on a gibbsite slab with 1 or 2 intervening water layers (1W or 2W). The structures of the gibbsite surfaces and the confined water are analyzed to determine how the number and type of hydrogen bonds (H-bonds) influence the free energy profile during the translation. The PMF profile exhibits a periodicity of length 5.078 Å, consistent with the gibbsite unit cell size along the translation direction. The changes in the surface–water and water–water hydrogen bond network and water and surface OH groups’ orientations during the translation are strongly coupled with the changes in the PMF profile in the 1W case. However, when increasing the number of intervening water layers from 1W to 2W, the particle/slab misalignment becomes a dominant factor controlling the behavior of the PMF profile. We also establish a method to quantify misalignment between the particle and the slab, which exhibits a strong correlation with the free energy for the 2W case. These results shed more light into the roles of particle/slab misalignment and hydrogen bond network in the OA of mineral particles in aqueous solution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

In this work, the effect of plasma turbulence on the trajectories of dust particles is investigated for the first time. The dynamics of dust particles is computed using the ad hoc developed Dust Injection Simulator code, using a 3D turbulent plasma background computed with the TOKAM3X code. As a result, the evolution of the particle trajectories is governed by the ion drag force, and the shape of the trajectory is set by the Stokes number St∝a d /n 0 , with a d the dust radius and n 0 the density at the separatrix. The plasma turbulence is observed to scatter the dust particles, exhibiting a hyperdiffusive regime in all cases. The amplitude of the turbulent spread of the trajectories Δr 2 is shown to depend on the ratio Ku/St, with Ku∝u rms the Kubo number and urms the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting Δr 2 ∝(Ku/St) 2 t 3 , or Δr 2 ∝(u rms n 0 /a d ) 2 t 3 . As the dust is heated by the plasma fluxes, thermionic emission sets the dust charge, originally negative, to slightly positive values. This results in a substantial reduction of the ion drag force through the suppression of its Coulomb scattering component. The dust grain inertia is then no longer negligible and drives the transition from a hyperdiffusive regime toward a ballistic one.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Editorial announcements and information — ICFA Beam Dynamics Newsletter No. 81

This is a general information article of the ICFA Beam Dynamics Newsletter No. 81. Thus, it contains two forewords from the Editor-in-Chief and the issue editor, a workshop and conference report section, a recent Ph.D. thesis section, a forthcoming beam dynamics events section, and a section of announcements from the beam dynamics panel.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Computational approaches to Coherent Synchrotron Radiation in two and three dimensions

Coherent Synchrotron Radiation (CSR) is an important and often detrimental effect in particle accelerators. While one-dimensional models have been successfully used to design and explain the behavior of modern machines, questions remain about their domain of validity. In recent years, two- and three-dimensional models have been developed that are amenable to efficient numerical computation. This article gives an overview of CSR computation from its discovery through the present state of the art.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Application of machine learning in the determination of impact parameter in the 132 Sn+ 124 Sn system

Here, 132 Sn + 124 Sn collisions at a beam energy of 270 MeV/nucleon were performed at the Radioactive Isotope Beam Factory (RIBF) in RIKEN to investigate the nuclear equation of state. Reconstructing the impact parameter is one of the important tasks in the experiment as it relates to many observable. In this work, we employ three commonly used algorithms in machine learning, the artificial neural network (ANN), the convolutional neural network (CNN), and the light gradient boosting machine (LightGBM), to determine the impact parameter by analyzing either the charged particle spectra or several features simulated with events from the ultrarelativistic quantum molecular dynamics (UrQMD) model. To closely imitate experimental data and investigate the generalizability of the trained machine learning algorithms, incompressibility of nuclear equation of state and the in-medium nucleon-nucleon cross sections are varied in the UrQMD model to generate the training data. The mean absolute error Δb between the true and the predicted impact parameter is smaller than 0.45 fm if training and testing sets are sampled from the UrQMD model with the same parameter set. However, if training and testing sets are sampled with different parameter sets, Δb would increase to 0.8 fm. The generalizability of the trained machine learning algorithms suggests that these machine learning algorithms can be used reliably to reconstruct the impact parameter in experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗