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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Simulating Meson Scattering on Spin Quantum Simulators

Studying high-energy collisions of composite particles, such as hadrons and nuclei, is an outstanding goal for quantum simulators. However, preparation of hadronic wave packets has posed a significant challenge, due to the complexity of hadrons and the precise structure of wave packets. This has limited demonstrations of hadron scattering on quantum simulators to date. Observations of confinement and composite excitations in quantum spin systems have opened up the possibility to explore scattering dynamics in spin models. In this article, we develop two methods to create entangled spin states corresponding to wave packets of composite particles in analog quantum simulators of Ising spin Hamiltonians. One wave-packet preparation method uses the blockade effect enabled by beyond-nearest-neighbor Ising spin interactions. The other method utilizes a quantum-bus-mediated exchange, such as the native spin-phonon coupling in trapped-ion arrays. With a focus on trapped-ion simulators, we numerically benchmark both methods and show that high-fidelity wave packets can be achieved in near-term experiments. We numerically study scattering of wave packets for experimentally realizable parameters in the Ising model and find inelastic-scattering regimes, corresponding to particle production in the scattering event, with prominent and distinct experimental signals. Our proposal, therefore, demonstrates the potential of observing inelastic scattering in near-term quantum simulators.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Scalable multiscale modeling of platelets with 100 million particles

Here, we developed the core components of the AI-aided multiple time stepping algorithm for multiscale modeling of cell dynamics. This algorithm was implemented and analyzed on two supercomputer architectures with an application of simulating the aggregation of 250 platelets, or 102 million particles. To scale on these computers with complex memory and network architectures with GPUs, we devised a biomechanics-informed task mapping scheme to optimize load imbalance, communications, and memory utilization. Our simulations, scaling well up to 192 nodes on a Summit-like supercomputer with a peak speed of 11 petaflops, achieved a rate of 423 μs/day which is 500 times faster than the conventional algorithm using static time step and this has enabled studies of record size blood clots at record spatial–temporal resolutions. Additionally, we discovered the sensitive dependence of the scalability and execution time on the methods of decomposition, CPU–GPU coupling, and task mapping.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1999

MRI Measurements and Granular Dynamics Simulation of Segregation of Granular Mixture

A counter intuitive axial segregation phenomenon in a rotating horizontal cylinder has recently captured attention of many researchers in different disciplines. There is a growing consensus that the interplay between the particle dynamics and the evolution of the internal structure during the segregation process must be carefully investigated. Magnetic resonance imaging (MRI) has been used to non-invasively obtain much needed dynamic/static information such as velocity and concentration profiles, and it has proven to be capable of depicting the evolution of segregation processes. Segregation in a rotating cylinder involves two processes: the first is to transport small particles in the radial direction to form a radial core, and the second is to transform the radial core into axially segregated bands. Percolation and/or "stopping" have been proposed as mechanisms for the radial segregation. As to mechanisms for axial band formation, much less is known. The difference in the dynamic angle of repose has been proposed to segregate different components in the axial direction. Recently, Hill and Kakalios have reported that particles mix or demix depending upon the competition between diffusion and preferential drift whose order can be determined by the dynamic angle of repose through the adjustment of the rotation rate. We claim that the dynamic angle of repose could be one of the causes, however, it fails to offer reasonable explanations for certain aspects of the axial migration. For example, we always observe that the radial segregation precedes the axial segregation and small particles migrate in the radial direction to form an axially extended radial core. It then transforms into axially segregated bands. By definition, the effects of the dynamic angle of repose are restricted near the free surface where the flowing layer is present. However, during the process of transforming from the radially segregated core to axially segregated bands, small particles located in the deep core region, which is untouched by the flowing layer, also completely disappear. Usually, the dynamics angle of repose are uniquely defined for individual species to characterize particle properties, and the dynamic angle of repose thus defined provides little information for the dynamic angle of repose of the mixture since the concentration ratio and the internal packing structure do not remain the same during the segregation processes. Under microgravity environment, the dynamics angle of repose argument does not hold since there is simply no flowing layer to influence/determine the preferred directions of segregation. We have thus designed an experiment so that the effects of the dynamic angle of repose can be minimized by filling the cylinder almost completely full. Small particles still formed a radial core and also migrated to form axial bands. As ground based experiments we have designed and conducted both 2D and 3D segregation experiments. The 2D experiments are performed using a thin cylinder (the gap between two end caps is about 5 mm) filled with different combinations of particles. The 3D experiments are conducted with a long cylinder of its length and diameter of 27cm and 7cm, respectively. Results of 2D experiments indicate that different mechanisms govern particle motion in regions near and far from the axis of rotation. Results of 3D experiments indicate that a series of collapses of microstructures of particle packing (micro-collapses) may be responsible for the creation of voids for small particles to migrate through in the axial direction. We have successfully eliminated the dynamic angle of repose as a cause for segregation, however, by almost completely filling the cylinder with the particles, we have lost an opportunity to investigate a possibility of particle "mobility" being a cause for segregation which requires a flowing surface but not the difference in the angle of repose. This is currently being investigated.

Nakagawa, M.↗

DOE OSTI2021

Contribution of Electron Pressure to Ring Current and Ground Magnetic Depression Using RAM-SCB Simulations and Arase Observations During 7–8 November 2017 Magnetic Storm

Understanding the physical processes that control the dynamics of energetic particles in the inner magnetosphere is important for both space-borne and ground-based assets essential to the modern society. The storm time distribution of ring current particles in the inner magnetosphere depends strongly on their transport in the evolving electric and magnetic fields along with particle acceleration and loss. Here, we investigated the ring current particle variations using observations and simulations. We compared the ion (H + , He + , and O + ) and electron flux and plasma pressure variations from Arase observations with the self-consistent inner magnetosphere model: Ring current Atmosphere interactions Model with Self Consistent magnetic field (RAM-SCB) during the 7–8 November 2017 geomagnetic storm. We investigated the contribution of the different species (ions and electrons) to the magnetic field deformation observed at ground magnetic stations (09°–45° MLat) using RAM-SCB simulations. The results show that the ions are the major contributor with ~88% and electrons contribute ~12% to the total ring current pressure. It is also found that the electron contribution is non-negligible (~18%) to the ring current in dawn-side during the main phase of the storm. Thus, the electron contribution to the storm time ring current is important and should not be neglected.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Improving the electrostatic design of the Jefferson Lab 300 kV DC photogun

The 300 kV DC high voltage photogun at Jefferson Lab was redesigned to deliver electron beams with a much higher bunch charge and improved beam properties. The original design provided only a modest longitudinal electric field ( E z ) at the photocathode, which limited the achievable extracted bunch charge. To reach the bunch charge goal of approximately few nC with 75 ps full-width at half-maximum Gaussian laser pulse width, the existing DC high voltage photogun electrodes and anode–cathode gap were modified to increase E z at the photocathode. In addition, the anode aperture was spatially shifted with respect to the beamline longitudinal axis to minimize the beam deflection introduced by the non-symmetric nature of the inverted insulator photogun design. We present the electrostatic design of the original photogun and the modified photogun and beam dynamics simulations that predict vastly improved performance. We also quantify the impact of the photocathode recess on beam quality, where recess describes the actual location of the photocathode inside the photogun cathode electrode relative to the intended location. A photocathode unintentionally recessed/misplaced by sub-millimeter distance can significantly impact the downstream beam size.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

What Composition of High-Energy Physics Collaborations is Epistemically Optimal?

Computer simulations have recently come to the fore as a crucial tool for studying epistemic dynamics of scientific collectives. One pressing issue in big-ticket, long-term, and crowded particle physics research is the optimal organization of collaborations. This research developed a computation model of collaborations in high-energy physics and performed simulations to investigate the epistemic efficiency of groups and its dependence on the size of such groups and their composition: namely, percentages of pure, partially theoretically competent, and fully theoretically competent experimentalists. The present study reveals that in both small (100-member) and large (3000-member) groups, epistemic payoff as the measure of epistemic efficiency of collaborations is, generally, positively correlated with an increase in both the number of theoretically fully competent and pure (incompetent) experimentalists. Here, we are less certain of such conclusion in the case of the small collaborations. Although the subcommunity of experimentalists who are partially expert is not found to be immediately epistemically beneficial for collaborations, nevertheless they also crucially serve as a transitional community between the theorists and the experimentalists. Predicated on the toy model simulations, I suggest that institutions should provide measures to assist members of the latter subgroup in progressing toward developing a full theoretical expertise.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRSDate not supplied

Advances in Plume-Surface Interaction Simulation Capability Under Lunar Landing Conditions

Plume-Surface Interaction (PSI) between lander engine plumes and native regolith soil presents primary hazards in obscuration and contamination by particle clouds, high-energy ejecta streams, and landing area cratering damage. The MSFC Fluid Dynamics Branch has assembled a portfolio of simulation tools to develop a predictive PSI capability for NASA customers such as the Human Lander System (HLS) and Commercial Lunar Payload Services (CLPS). These tools are matured via funding by NASA’s STMD Game Changing Development Program and the NASA SBIR/STTR Program.

Peter A Liever↗

DOE OSTI2022

Charged particle motion and radiation in strong electromagnetic fields

The dynamics of charged particles in electromagnetic fields is an essential component of understanding the most extreme environments in our Universe. In electromagnetic fields of sufficient magnitude, radiation emission dominates the particle motion and effects of quantum electrodynamics (QED) in strong fields are crucial, which triggers electron-positron pair cascades and counterintuitive particle-Trapping phenomena. As a result of recent progress in laser technology, high-power lasers provide a platform to create and probe such fields in the laboratory. With new large-scale laser facilities on the horizon and the prospect of investigating these hitherto unexplored regimes, this review explores the basic physical processes of radiation reaction and QED in strong fields, how they are treated theoretically and in simulation, the new collective dynamics they unlock, recent experimental progress and plans, and possible applications for high-flux particle and radiation sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-Landau closure model *

Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

DOE OSTI2020

Machine Learning-Based Upscaling of Finite-Size Molecular Dynamics Diffusion Simulations for Binary Fluids

Molecular diffusion coefficients calculated using molecular dynamics (MD) simulations suffer from finite-size (i.e., finite box size and finite particle number) effects. Results from finite-sized MD simulations can be upscaled to infinite simulation size by applying a correction factor. For self-diffusion of single-component fluids, this correction has been well-studied by many researchers including Yeh and Hummer (YH); for binary fluid mixtures, a modified YH correction was recently proposed for correcting MD-predicted Maxwell–Stephan (MS) diffusion rates. In this study we use both empirical and machine learning methods to identify improvements to the finite-size correction factors for both self-diffusion and MS diffusion of binary Lennard-Jones (LJ) fluid mixtures. Using artificial neural networks (ANNs), the error in the corrected LJ fluid diffusion is reduced by an order of magnitude versus existing YH corrections, and the ANN models perform well for mixtures with large dissimilarities in size and interaction energies where the YH correction proves insufficient.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Hyperdiffusion of dust particles in a turbulent tokamak plasma

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\propto 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 $\Delta r^2$ is shown to depend on the ratio $Ku/St$, with $Ku\propto u_{rms}$ the Kubo number and $u_{rms}$ the fluctuation level of the plasma flow. These results are compared with a simple analytical model, predicting $\Delta r^2\propto (Ku/St)^2t^3$, or $\Delta r^2\propto (u_{rms}n_0/a_d)^2t^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 towards a ballistic one.

dust↗

DOE OSTI2020

ddcMD.os

ddcMD is a general purpose molecular dynamics (MD) code that supports MPI parallelism. MD codes are used for simulation of particles systems and capture all the many-body effects of the underlying particle potential that defines the physical system. Though MD can be used to model systems from the subatomic to astrological length scales ddcMD is mainly focused on the atomic scale length scale, In this release of ddcMD the support will be mainly for systems using the coarse-grain Martini potential, a particle potential for biological systems.

Glosli, JamesN↗

DOE OSTI2022

Understanding Phase and Interfacial Effects of Spall Fracture in Additively Manufactured Ti-5Al-5V-5Mo-3Cr

Additive manufactured Ti-5Al-5V-5Mo-3Cr (Ti-5553) is being considered as an AM repair material for engineering applications because of its superior strength properties compared to other titanium alloys. Here, we describe the failure mechanisms observed through computed tomography, electron backscatter diffraction (EBSD), and scanning electron microscopy (SEM) of spall damage as a result of tensile failure in as-built and annealed Ti-5553. We also investigate the phase stability in native powder, as-built and annealed Ti-5553 through diamond anvil cell (DAC) and ramp compression experiments. We then explore the effect of tensile loading on a sample containing an interface between a Ti-6Al-V4 (Ti-64) baseplate and additively manufactured Ti-5553 layer. Post-mortem materials characterization showed spallation occurred in regions of initial porosity and the interface provides a nucleation site for spall damage below the spall strength of Ti-5553. Preliminary peridynamics modeling of the dynamic experiments is described. Finally, we discuss further development of Stochastic Parallel PARticle Kinteic Simulator (SPPARKS) Monte Carlo (MC) capabilities to include the integration of alpha (α)-phase and microstructural simulations for this multiphase titanium alloy.

36 MATERIALS SCIENCE↗

NASA NTRS1978

The dynamics and control of large flexible space structures. Part A: Discrete model and modal control

Attitude control techniques for the pointing and stabilization of very large, inherently flexible spacecraft systems were investigated. The attitude dynamics and control of a long, homogeneous flexible beam whose center of mass is assumed to follow a circular orbit was analyzed. First order effects of gravity gradient were included. A mathematical model which describes the system rotations and deflections within the orbital plane was developed by treating the beam as a number of discretized mass particles connected by massless, elastic structural elements. The uncontrolled dynamics of the system are simulated and, in addition, the effects of the control devices were considered. The concept of distributed modal control, which provides a means for controlling a system mode independently of all other modes, was examined. The effect of varying the number of modes in the model as well as the number and location of the control devices were also considered.

Bainum, P. M.↗