Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle dynamics”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Inertial Effects in Suspension Dynamics

The present work analyses the dynamics of a suspension of heavy particles in shear flow. The magnitude of the particle inertia is given by the Stokes number St = m(gamma/6(pi)a, which is the ratio of the viscous relaxation time of a particle tau(sub p) = m=6pi(eta)a to the flow time gamma(sup -1). Here, m is the mass of the particle, a is its size, eta is the viscosity of the suspending fluid and gamma is the shear rate. The ratio of the Stokes number to the Reynolds number, Re = (rho)f(gamma)a(exp 2)/eta, is the density ratio rho(sub p)/rho(sub f). Of interest is to understand the separate roles of particle (St) and fluid (Re) inertia in the dynamics of suspensions. In this study we focus on heavy particles, rho(sub p)/rho(sub f) much greater than 1, for which the Stokes number is finite, but the Reynolds number is sufficiently small for inertial forces in the fluid to be neglected; thus, the fluid motion is governed by the Stokes equations. On the other hand, the probability density governing the statistics of the suspended particles satisfies a Fokker-Planck equation that accounts for both configuration and momentum coordinates, the latter being essential for finite St. The solution of the Fokker-Planck equation is obtained to O(St) via a Chapman-Enskog type-procedure, and the conditional velocity distribution so obtained is used to derive a configuration-space Smoluchowski equation with inertial corrections. The inertial effects are responsible for asymmetry in the relative trajectories of two spheres in shear flow, in contrast to the well known symmetric structure in the absence of inertia. Finite St open trajectories in the plane of shear suffer a downward lateral displacement resulting from the inability of a particle of finite mass to follow the curvature of the zero-Stokes-number pathlines. In addition to the induced asymmetry, the O(St) inertial perturbation dramatically alters the nature of the near-field trajectories. The stable closed orbits (for St = 0) in the plane of shear now spiral in, approaching particle-particle contact in the limit. All trajectories starting from an initial offset of O(St(sup 1/2) or less (which remain open for St = 0) also spiral in. The asymmetry of the trajectories leads to a non-Newtonian rheology and diffusive behavior. The latter because a given particle (moving along a finite St open trajectory) suffers a net displacement in the transverse direction after a single interaction. A sequence of such uncorrelated displacements leads to the particle executing a random walk. The inertial diffusivity tensor is anisotropic on account of differing strengths of interaction in the gradient and vorticity directions. Since the entire region (constituting an in finite area) of closed orbits in the plane of shear spirals onto contact for #finite St, the latter represents a singular surface for the pair-distribution function. The exact form of the pair-distribution function at contact is still, however, indeterminate in the absence of non-hydrodynamic effects. It should also be noted that finite St non-rectilinear flows do not support a spatially uniform number density owing to the cross-streamline inertial migration of particles.

J. F. Brady↗
NASA NTRS1994

Coalescence of two current loops with a kink instability simulated by a three-dimensional electromagnetic particle code

We have studied the dynamics of a coalescence of current loops using three-dimensional electromagnetic (EM) particle simulation code. Our focus is the investigation of such kinetic processes as energy trasnfer, heating particles, and electromagnetic emissions associated with a current loop coalescence which cannot be studied by MHD simulations. First, the two loops undergo a pinching oscillation due to a pressure imbalance between the inside and outside of the current loop. During the pinching oscillation, a kinetic kink instability is excited and electrons in the loops are heated perpendicularly to an ambient magnetic field. Next, the two current loops collide and coalesce, while at the same time a helical structure grows further. Subsequently, the perturbed current, which is due to these helically bunched electrons, can drive a whistler instability. It should be noted in this case that the whistler wave is excited by the kinetic kink instability and not a beam instability. After the coalescence of two helical loops, tilting motions can be observed in the direction of left-hand rotation, and the helical structure will relax resulting in strong plasma heating mostly in the direction perpendicular to the ambient magnetic field. It is also shown that high-frequency electromagnetic waves can be emitted from the region where the two loops coalesce and propagate strongly in the direction of the electron drift velocity. These processes may be important in understanding heating mechansims for coronal loops as well as radio wave emission mechanisms from active regions of solar plasmas.

Nishikawa, K.-I.↗
DOE OSTI2026

Structure and Flow-Viscosity of Filled-Polymer-Based 3D Printing Ink: Exploration through Coarse-Grained Molecular Dynamics

The addition of nanofiller particles to a polymer matrix has long been known to enhance or modify the composite’s mechanical and rheological properties. However, quantitatively capturing such changes with molecular level simulations remains computationally challenging. Toward that goal, we performed coarse-grained molecular dynamics of a nanocomposite system at a fixed (25 vol %) filler loading under nonspecific, weak polymer–filler interactions representative of a broad class of technologically important materials. We report several interesting results, including: (1) the equilibrium chain-configuration remains Gaussian-like as in an unfilled melt; (2) smaller filler particles display a stronger tendency to cluster; (3) larger fillers act as plasticizers by reducing the entanglement density and accelerating the chain mobility; and (4) fillers enhance the tensile response modulus, with the effect being stronger for larger particles. We also simulate cluster breakup, yielding, and elongational flow under an applied time-linear tensile strain and study the flow viscosity as a function of filler-size and chain-length.

Materials science↗
NASA NTRS1983

High-energy particles

It is pointed out that the magnetosphere of Jupiter is in many respects quite different from that of the earth. The energy required to drive the Jovian magnetosphere is apparently extracted from Jupiter's rotational energy rather than from the solar wind. Jupiter is a strong source of energetic charged particles which can be detected as far away as the orbit of Mercury. The structure and dynamics of the energetic particle distribution in the inner magnetosphere is discussed, taking into account observations, transport and losses in the inner magnetosphere, satellite interactions, and electron synchrotron radiation. The subsolar hemisphere is considered, giving attention to particle fluxes in the subsolar magnetosphere, conditions in the middle magnetosphere, and the characteristics of the outer magnetosphere. A description of the predawn magnetosphere is also provided.

Schardt, A. W.↗
NASA NTRS1992

Statistical mechanics of violent relaxation

We propose a functional that is extremized through violent relaxation. It is based on the Ansatz that the wave-particle scattering during violent dynamical processes can be approximated as a sequence of discrete scattering events that occur near a particle's perigalacticon. This functional has an extremum whose structure closely resembles that of spheroidal stellar systems such as elliptical galaxies. The results described here, therefore, provide a simple framework for understanding the physical nature of violent relaxation and support the view that galaxies are structured in accord with fundamental statistical principles.

Spergel, David N.↗
NASA NTRSDate not supplied

Investigation of the Alfven Wave Turbulence on the Dynamics of SEPs

Solar energetic particles (SEPs) are high-energy particles originating from the Sun that pose potential risks to space missions, especially those outside Earth's protective magnetosphere. Understanding their behavior within the heliosphere and Earth's magnetosphere is vital for the safety and functionality of space exploration. The study focuses on how the pitch angle scattering, the angle change between a particle's velocity and the magnetic field due to interactions with magnetic irregularities, impacts the decay phase of SEP events. By solving the focused transport equation, the study models the transport of SEPs along magnetic field lines, extending these lines up to 5 AU (Astronomical Units) to understand their behavior even beyond the typical proximity to the Sun. Incorporating the effects of pitch angle scattering at such vast distances is crucial for understanding the decay dynamics of a SEP event. This modeling of SEPs is integrated with simulations of other vital space phenomena: the solar wind, interplanetary magnetic field, and Alfven wave turbulence. This presentation covers the modeling approach used in this research and the role of pitch angle scattering at various heliocentric distances on the dynamics of the SEP events' decay phase.

solar energetic particles↗
DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗
DOE OSTI2024

Investigation of thermal damage in explosive bridgewire detonators via discrete element method simulations

Exploding bridgewire (EBW) detonators are used to rapidly and reliably initiate energetic reactions by exploding a bridgewire via Joule heating. While the mechanisms of EBW detonators have been studied extensively in nominal conditions, comparatively few studies have addressed thermally damaged detonator operability. We present a mesoscale simulation study of thermal damage in a representative EBW detonator, using discrete element method (DEM) simulations that explicitly account for individual particles in the pressed explosive powder. We use a simplified model of melting, where solid spherical particles undergo uniform shrinking, and fluid dynamics are ignored. The subsequent settling of particles results in the formation of a gap between the solid powder and the bridgewire, which we study under different conditions. In particular, particle cohesion has a significant effect on gap formation and settling behavior, where sufficiently high cohesion leads to coalescence of particles into a free-standing pellet. This behavior is qualitatively compared to experimental visualization data, and simulations are shown to capture several key changes in pellet shape. We derive a minimum and maximum limit on gap formation during melting using simple geometric arguments. In the absence of cohesion, results agree with the maximum gap size. With increasing cohesion, the gap size decreases, eventually saturating at the minimum limit. In conclusion, we present results for different combinations of interparticle cohesion and detonator orientations with respect to gravity, demonstrating the complex behavior of these systems and the potential for DEM simulations to capture a range of scenarios.

cohesive powders↗
DOE OSTI2026

Rotational memory function of SPC/E water

Memory effects are essential for the dynamics of condensed materials and are responsible for non-exponential relaxation of correlation functions of dynamic variables through the memory function. Memory functions of dipole rotations for water have never been calculated directly from molecular dynamics simulations. We present here calculations of memory functions for single-dipole rotations and for the overall dipole moment of the sample for SPC/E water. The normalized memory functions for single-particle and collective dipole dynamics turn out to be nearly identical. This result validates theories of dielectric spectroscopy in terms of single-particle time correlation functions and the connection between the collective and single-particle relaxation times through the Kirkwood factor. The dielectric function in this formalism contains no new dynamic information that does not exist in the single-dipole correlation function. A short memory time, ≲1 fs, justifies the use of the mathematics of rotational diffusion to describe the dynamics of a single molecular dipole moment in bulk water. Here, an analytical equation for the rotational memory time is derived.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2026

Experimental Studies of Nonlinear Integrable Optics with Elliptic Potentials

The stable transport of charged particle beams is the core challenge in designing and operating accelerators. The current paradigm for transverse focusing is based on quadrupole and dipole elements, described as a linear integrable Hamiltonian by Courant and Snyder. The operational limits of high intensity accelerators are often determined by collective instabilities within the beam. Nonlinear elements may be added to suppress certain forms of these collective effects, but restrict the range of stable trajectories. These shortcomings motivate extending to a nonlinear integrable system which can offer the benefits of suppressing collective instabilities without limiting the stable trajectories. Danilov and Nagaitsev proposed a novel nonlinear integrable system with elliptical potentials, which could be implemented as magnetic elements in an accelerator. Such a system has been implemented for practical verification in the integrable optics test accelerator (IOTA), a small storage ring constructed for beam dynamics studies. Electron beam studies in IOTA treat the low-emittance beam as a macroparticle for detailed probing of the expected single particle dynamics. This dissertation details new measurements of the nonlinear integrable system relevant for practical implementation. Turn-by-turn measurements of the kicked beam responses are used for phase space reconstruction and analysis of the predicted nonlinear dynamics. The stability and aperture for various nonlinear configurations are measured with beam losses. Amplitude dependent detuning, a core figure of merit for suppressing instabilities, is measured and compared with high fidelity particle tracking simulations. The synchrotron radiation images of circulating beam allow direct measurements of the topology and lifetime in configurations where nonlinear focusing terms dominate.

Wieland, John [Michigan U.] (ORCID:000000028971852↗
DOE OSTI2025

Strain-associated nanoscale fluctuating lithium transport within single-crystalline LiNi 1/3 Mn 1/3 Co 1/3 O 2 cathode particles

Solid-state lithium diffusion dynamics are critical for the rate capability and longevity of Li-ion batteries. Conventionally, nanoscale lithium diffusion within individual battery particles has been simplified as being primarily driven by concentration gradients, despite the associated processes inducing local lattice expansion, contraction, and strain fields. Using operando scanning transmission soft X-ray microscopy with high spatial resolution and chemical sensitivity to track nanoscale intraparticle lithium transport, and post-cycling Bragg coherent diffraction X-ray imaging to directly reveal three-dimensional intraparticle strain fields, we uncover strain-associated lithium transport dynamics within single-crystalline LiNi 1/3 Mn 1/3 Co 1/3 O 2 (scNMC) particles during cycling. Contrary to the expected thermodynamic solid-solution behavior of scNMC, our observations reveal near-uniform but fluctuating regions of lithium-dense and lithium-dilute areas during cycling. These fluctuations suggest that nanoscale lithium diffusion can proceed counter to concentration gradients. Additionally, we demonstrate that an increased presence of lithium-dilute regions near the surface enhances lithium surface insertion kinetics, emphasizing the importance of controlling surface lithium distribution to improve rate performance. Our study provides insights into nanoscale solid-state ion transport, with potential applications in batteries, solid-state fuel cells, and memristors.

Lee, Danwon [Seoul National Univ. (Korea, Republic↗
DOE OSTI2025

Engineering spin coherence in core-shell diamond nanocrystals

Fluorescent diamond nanocrystals can host spin qubit sensors capable of probing the physical properties of biological systems with nanoscale spatial resolution. Sub-100 nm diamond nanosensors can readily be delivered into intact cells and even living organisms. However, applications beyond current proof-of-principle experiments require a substantial increase in sensitivity, which is limited by surface induced charge instability and electron-spin dephasing. In this work, we utilize engineered core-shell structures to achieve a drastic increase in qubit coherence times (T 2 ) from 1.1 to 35 μs in bare nanodiamonds to upward of 52 to 87 μs. We use electron-paramagnetic-resonance results to present a band bending model and connect silica encapsulation to the removal of deleterious mid-gap surface states that are negatively affecting the qubit’s spin properties. Combined with a 1.9-fold increase in particle luminescence these advances correspond to up to two-order-of-magnitude reduction in integration time. Probing qubit dynamics at a single particle level further reveals that the noise characteristics fundamentally change from a bath with spins that rearrange their spatial configuration during the course of an experiment to a more dilute static bath. The observed results shed light on the underlying mechanisms governing fluorescence and spin properties in diamond nanocrystals and offer an effective noise mitigation strategy based on engineered core-shell structures.

core-shell↗
DOE OSTI2025

Thermal Analysis of a Solid Particle Light-Trapping Planar Cavity Receiver Using Computational Fluid Dynamics

Concentrated solar power (CSP) is one of the most effective ways of harnessing solar power to create efficient, durable, and resilient energy systems. This study entails thermal modeling and analysis of a novel central tower receiver configuration. This receiver uses solid particles as the heat transfer fluid (HTF), a promising option for third-generation CSP systems. The configuration considered here is the light-trapping planar cavity receiver (LTPCR) introduced by the National Renewable Energy Laboratory. While heat transfer studies of various LTPCR subsystems have been done, system-level thermal analysis of the LTPCR receiver has not been attempted. This study also presents important sensitivity analyses of the operating parameters of the CSP system, which can help guide the design of future central tower receivers. This study employs Ansys Fluent as a computational fluid dynamics (CFD) tool to model fluid dynamics and heat transfer in the receiver, intending to quantify its thermal performance. The model seamlessly integrates Monte Carlo ray tracing data, which generates absorbed solar flux profiles from the heliostat field design, with the heat transfer characteristics of the fluidized particle bed. This unified model is designed to accurately predict the thermal behavior of the LTPCR. Analysis of preliminary results reveals that the primary loss mechanisms are radiative and natural convective losses, in that order. Based on observations from a baseline case, several strategies are suggested and numerically tested. These solutions include selective cooling of high-temperature regions and manipulation of particle bed parameters. Selective cooling of high-temperature regions reduced the peak temperature by 151 degrees C and decreased thermal losses by 0.9%. Improving the particle-wall heat transfer coefficient (P-W HTC) of the particle bed decreased the thermal losses by 1.7% and decreased the peak temperatures by 57 degrees C. Decreasing the particle inlet temperature (PIT) also reduced thermal losses by 3.5% and decreased peak temperatures by 29 degrees C. Compounding these strategies improved the thermal losses of the receiver from 13.5% in the baseline case to 7.5%. Additionally, the study explores the variation in thermal performance across different locations of the receiver, where a variation of thermal losses from 12.9% to 17.3% is found. This allows a comprehensive evaluation of potential improvements in efficiency and temperature management.

computational fluid dynamics↗
DOE OSTI2024

Towards an understanding of particle-scale flaws and microstructure evolution in cold spray via accumulation of single particle impacts

Cold spray coatings are the sum of countless individual bonding events between single particles impacting on top of one another at high velocities. Thus, the collective behavior of microparticles must be considered to elucidate the origins of coating flaws at the scale of the particles and larger, or the dynamic evolution of the overall coating microstructure. Laser-induced particle impact testing (LIPIT) has been extensively used to study single-particle impacts, and in this work is adapted to study the accumulation of numerous particles with knowledge of each individual particle's impact parameters (particle size, velocity). The method reproducibly deposits stacks of gold particles (>20 particles) with different characteristic spectra of impact velocity. The observation of impact-induced erosion lets us define a critical velocity for material-build-up that is higher than that for single-particle bonding. The quantitative single-particle data are analyzed in a correlative manner to the structure and flaws in the resulting stacks, providing some first statistical connections between, e.g., strain and recrystallization, or aberrant particle characteristics and defects. The results highlight opportunities for the study of many-particle phenomena in microparticle impact—from interaction of particles in cold spray to multi-step erosion processes—with a quantitative view of the behavior of single particles.

Erosion↗
NASA NTRS2000

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗