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 55 records · Page 3

Self-Assembly of Anisotropic Particles on Curved Surfaces

The surface curvature of membranes, interfaces, and substrates plays a crucial role in shaping the self-assembly of particles adsorbed on these surfaces. However, little is known about the interplay between particle anisotropy and surface curvature and how they couple to alter the free energy landscape of particle assemblies. Using molecular dynamics simulations, we investigate the effect of prescribed curvatures on a quasi-2D assembly of anisotropic patchy particles. By varying curvature and surface coverage, we uncover a rich geometric phase diagram, with curvature inducing ordered structures entirely absent on planar surfaces. Large spatial domains of ordered structures can contain hidden microdomains of orientational textures imprinted by the surface on the assembly. The dynamical landscape is also reshaped by surface curvature, with a glass-like state emerging at modest densities and high curvature. Changes to the symmetry of the surface curvature are found to result in distinct structures, including phases with mesoscale ordering. In conclusion, our findings show that the coupling between surface curvature and particle geometry opens an unexplored space of morphologies and structures that can be exploited for material design.

Fluids↗

Stokesian dynamics with odd viscosity

Stokesian dynamics is a well-established computational method for simulating dynamics of many particles suspended in a conventional passive fluid medium. Active fluids composed of self-driving particles with broken time-reversal symmetry permit the emergence of a so-called odd viscosity. Here, in this work, we extended the conventional Stokesian dynamics formalism to incorporate the additional hydrodynamic effects due to odd viscosity, which enables simulating collective behavior of many particles suspended in an active fluid medium with both even viscosity and odd viscosity.

42 ENGINEERING↗

Comparison of Eulerian Bin and Lagrangian Particle-Based Schemes in Simulations of Pi Chamber Dynamics and Microphysics

This paper presents a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as superdroplets. Each superdroplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh–Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

54 ENVIRONMENTAL SCIENCES↗

Investigation and validation of the dynamic response of an acoustically levitated particle using the lattice Boltzmann method

The stable levitation of an analyte sample in an acoustic levitator is a primary requirement for accurate x-ray characterization of its scientific structure. A rigid particle oscillates in an under-damped manner when introduced into the node of established standing acoustic waves. This investigation has employed the lattice Boltzmann method (LBM), a computational fluid dynamics technique, for the analysis of such rigid particle dynamics in acoustic levitation. The simulation uses the two dimensional and nine velocity (D2Q9) Bhatnagar-Gross-Krook formulation to levitate a rigid 1.6 mm diameter nylon ($\rho$ = 1150 kg/m 3 ) particle in the air at standard pressure and temperature conditions. The presented work is the first reported simulation of realistic acoustic levitator boundary conditions using the LBM. The simulation can capture the particle-fluid interactions that produce dynamic levitation at less than one-period timescale in the ultrasonic frequency regime. An experiment was conducted by levitating a 1.6mm-diameter nylon sphere to estimate the oscillations, and the oscillating frequency was found to be 50 Hz. The dynamic simulation results are consistent with experimental results for particle oscillations within the same order of magnitude, indicating that LBM formulation can be successfully used to study acoustic levitation to understand and mitigate particle jitter. The distortion of the acoustic field due to a levitating particle's presence was also analyzed to demonstrate how the presence of the particle can disrupt adjacent levitating nodes.

42 ENGINEERING↗

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Glassy interphases reinforce elastomeric nanocomposites by enhancing percolation-driven volume expansion under strain

For nearly a century, introduction of nanoparticles to elastomers has yielded extraordinarily tough nanocomposites that are critical to technologies from actuators to tires. The mechanisms by which this reinforcement occurs have nevertheless remained a central open question in material science. One widely debated hypothesis posits that strong interactions between polymer and particles induce "glassy bridges" that cement particles into a cohesive percolating network that resists elongation. Here, molecular dynamics simulations show that glassy particle shells do not primarily provide elongational cohesion. Instead, they amplify an underlying mechanism wherein competition between filler and elastomer networks causes the elastomer's volume to increase on deformation. This induces contributions from the elastomer's bulk modulus, which is of order 1000 times larger than its Young's modulus. These findings establish a unified understanding of low-strain reinforcement in filled elastomers as emanating from volumetric competition between coexisting particulate and elastomeric networks. This reframes and unifies our understanding of low-strain reinforcement, provides a clear-cut diagnostic for the presence of glassy bridging, and offers a new design principle for tough elastomeric nanocomposites.

Computational Physics (physics.comp-ph)↗

Investigation of inter-subject variation in ultrafine particle deposition across human nasal airways: A study involving children, adults, and the elderly

Occupational and environmental exposure to toxic nanoparticles, driven by the rapid expansion of nanotechnology, raises significant respiratory health concern. Numerous studies have explored airflow and particle dynamics in adult nasal airways, but understanding the impact of age-related anatomical changes in children and the elderly remains limited. This study systematically investigates age-related anatomical variations and associated influence on nasal airflow dynamics and ultrafine particle deposition characteristics. Using Computational Fluid-Particle Dynamics (CFPD) method, simulation was conducted under diverse inhalation conditions spanning a wide age range, including: two children (5 years old), two young adults (in their twenties), and two elderly (over 77 years old). Our results reveal distinctive variations across age groups in anatomical dimensions, which affect distribution of wall shear stress where the elderly and children display unique patterns distinct from the young adults. While total deposition efficiency differs significantly between children and adults, filtration efficiency in the subregion with most deposition, main respiratory, remains consistent. However, inter-subject differences are observed in the vestibular and olfactory regions,emphasizing nuanced impact of age-related anatomical variations. Overall and subregional empirical equations for deposition efficiency were developed by incorporating the combined diffusion parameter, Sc a Δ b , corroborating the use of geometrical characteristic parameters for each specific subject in predicting nasal deposition efficiency across age groups. Our findings contribute to predictive nanoparticle exposure analysis in nasal airways across different age groups, thereby enhancing respiratory healthcare for individuals across the life span.

60 APPLIED LIFE SCIENCES↗

Nonorthogonal Cascade Catalysis in Multicompartment Micelles

Multicompartment micelles (MCMs) containing acid and base sites in discrete domains are prepared from poly(norbornene)-based amphiphilic bottlebrush copolymers in aqueous media. Here, the acid and base sites are localized in different compartments of the micelle, enabling the nonorthogonal reaction sequence: deacetalization – Knoevenagel condensation – Michael addition of acetals to 2-amino chromene derivatives. Computational simulations using dissipative particle dynamics (DPD) elucidated the bottlebrush composition required to effectively site-isolate the nonorthogonal catalysts. This contribution presents MCMs as a new class of nanostructures for one-pot multistep nonorthogonal cascade catalysis, laying the groundwork for the isolation of three or more incompatible catalysts to synthesize value-added compounds in a single reaction vessel, in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Is there an exact magnetic moment for charged particle motion in a time-dependent, homogeneous magnetic field?

The non-perturbative guiding-centre model provides an exact alternative to full-orbit simulations of charged particle dynamics in situations where traditional guiding-centre theory may fail. We demonstrate that the charged particle motion in a homogeneous, time-varying magnetic field is a solvable example of the non-perturbative guiding-centre model. This entails showing that the exact magnetic moment of Qin and Davidson can be constructed to be asymptotic to the adiabatic invariant series of Kruskal. In contrast to the perturbative invariant, the exact invariant contains information about parametric resonances. These resonances destroy the conservation of the usual magnetic moment over very long times. This refutes some previous claims about the all-time invariance of the magnetic moment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Mesocrystal growth through oriented sliding and attachment of nanoplates

Oriented attachment is a critical, yet poorly understood, crystal growth pathway based on the self-assembly of nanocrystals. During oriented attachment, solvent-separated particles align and coalesce through forces that enable precise rotation and translation. While prior studies emphasized intragap forces driving crystallographic alignment, the forces enabling uniform stacking and superlattice formation remain unclear. Here, we demonstrate how macroscopic gibbsite mesocrystals emerge from nanoplates guided into staggered positions by directional sliding. Electron microscopy and X-ray scattering reveal the monoclinic superlattice structure, based on nanoplate stacking with a uniform ≈50° stagger along the gibbsite [010] direction. In situ liquid-cell TEM captures preferential sliding along the gibbsite [010] direction, decelerating with increasing particle overlap. Molecular dynamics simulations reveal that this staggered arrangement corresponds to a global free-energy minimum, rather than full alignment. The simulations also confirm that sliding along the [010] direction is energetically favored and provide insight into the role of interfacial water in achieving long-range ordered assemblies. These insights highlight the energy landscape’s role in oriented attachment, with implications for material synthesis and hierarchical structures in nature.

36 MATERIALS SCIENCE↗

Background studies and normalization of signal events in the Mu2e experiment

The Mu2e experiment is currently being constructed at Fermilab to search for the neutrino-less conversion of negative muons into electrons in the field of an aluminum nucleus. The experiment aims at a sensitivity of four orders of magnitude higher than previous related experiments, which implies highly demanding accuracy requirements both in the design and during the operation. To achieve such a goal, two important tasks should be accomplished. First, it is essential to estimate precisely the particle yields and all the backgrounds that could mimic the monoenergetic conversion electron signal. Second, it is necessary to normalize the signal events accurately. The normalization of the signal events is planned to be done using a detector system made of an HPGe detector and a Lanthanum Bromide detector, which will measure the rate of muons stopped on the aluminum target by looking at the emitted characteristic X-and γ-rays of energies up to 1809 keV. Therefore, it is essential to evaluate the detector system's performance before the start of the actual experiment. In this study, the first task was addressed by an extensive campaign of Monte Carlo simulations to investigate the relevant parameters and their impact on the experiment's sensitivity. The second task was handled by taking advantage of the Helmholtz-Zentrum Dresden-Rossendorf (HZDR) pulsed Bremsstrahlung photon beam at the ELBE facility. The detector system was tested at the ELBE facility under timing and background conditions similar to the ones expected at the Mu2e experiment. The study presents and discusses the simulation results and the detector system testing campaign.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Microphysics of liquid water in sub-10 nm ultrafine aerosol particles

Abstract. Ultrafine aerosol particles with sizes smaller than 50 nm have been shown in recent studies to serve as a large source of cloud condensation nuclei that can promote additional cloud droplet formation under supersaturation conditions. Knowledge of the microphysics of liquid water in these droplets remains limited, particularly in the sub-10 nm particle size range, due to experimental and theoretical challenges associated with the complexity of aerosol components and the small length scales of interest (e.g., difficulty of precisely sampling the liquid–air interface, questionable validity of mean-field theoretical representations). Here, we carried out molecular dynamics (MD) simulations of aerosol particles with diameters between 1 and 10 nm and characterized atomistic-level structure and water dynamics in well-mixed and phase-separated systems with different particle sizes, NaCl salinities, and pimelic acid (PML) organic surface loadings as a function of distance from the time-averaged Gibbs dividing interface or instantaneous water–air interface. We define a sphericity factor (ϕ) that can shed light on the phase-mixing state of nanodroplets, and we reveal an unexpected dependence of mixing state on droplet size. Our results also evidence an ion concentration enhancement in ultrafine aerosols, which should modulate salt nucleation kinetics in sub-10 nm droplets, and provide detailed characterization of the influence of droplet size on surface tension and on water self-diffusivity near the interface. Analysis of water evaporation free energy and water activity demonstrates the validity of the Kelvin equation and Köhler theory at droplet sizes larger than 4 nm under moderate salinities and organic loadings and the need for further extension to account for ion concentration enhancement in sub-10 nm aerosols, droplet-size-dependent phase separation effects, and a sharp decrease in the cohesiveness of liquid water in sub-4 nm droplets. Finally, we show that an idealized fractional surface coating factor (fs) can be used to categorize and reconcile water accommodation coefficients (α*) observed in MD simulations and experimental results in the presence of organic coatings, and we resolve the droplet size dependence of α*.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Skyrmion soliton motion on periodic substrates by atomistic and particle-based simulations

Abstract We compare the dynamical behavior of magnetic skyrmions interacting with square and triangular defect arrays just above commensuration using both atomistic and particle-based models. Under applied drives, the initial motion is a kink traveling through the pinned skyrmion lattice. For the square defect array, both models agree well and show a regime in which the soliton motion is locked along 45°. The atomistic model also produces locking of a soliton along 30°, which is absent in the particle-based model. For the triangular defect array, the atomistic model exhibits 30° soliton motion over a wide region of external current values. In contrast, the particle-based model gives 45° soliton motion over a small range of external driving force values. The difference arises because the nondeforming particle model facilitates meandering skyrmion orbits while the deformable atomistic model enables stronger skyrmion-skyrmion interactions that reduce the meandering. Our results indicate that soliton motion through pinned skyrmion lattices on a periodic substrate is a robust effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗

Modeling nanoparticle charge distribution in the afterglow of non-thermal plasmas and comparison with measurements

Particle charging in the afterglows of non-thermal plasmas typically take place in a non-neutral space charge environment. We model the same by incorporating particle-ion collision rate constant models, developed in prior work by analyzing particle-ion trajectories calculated using Langevin Dynamics simulations, into species transport equations for ions, electrons and charged particles in the afterglow. A scaling analysis of particle charging and additional Langevin Dynamics calculations of the particle-ion collision rate constant are presented to extend the range of applicability to ion electrostatic to thermal energy ratios of 300 and diffusive Knudsen number (that scales inversely with gas pressure) up to 2000. The developed collision rate constant models are first validated by comparing predictions of particle charge against measured values in a stationary, non-thermal DC plasma from past PK-4 campaigns published in Phys. Rev. Lett. 93(8): 085001 and Phys. Rev. E 72(1): 016406). The comparisons reveal excellent agreement within ±35% for particles of radius 0.6,1.0,1.3 μm in the gas pressure range of ~20-150 Pa. The experiments to probe particle charge distributions by Sharma et al. (J. Physics D: Appl. Phys. 53(24): 245204) are modeled using the validated particle-ion collision rate constant models and the calculated charge fractions are compared with measurements. The comparisons reveal that the ion/electron concentration and gas temperature in the afterglow critically influence the particle charge and the predictions are generally in qualitative agreement with the measurements. Along with critical assessment of the modeling assumptions, several recommendations are presented for future experimental design to probe charging in afterglows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Induced phase transformation in ionizable colloidal nanoparticles

Abstract Acid–base equilibria directly influence the functionality and behavior of particles in a system. Due to the ionizing effects of acid–base functional groups, particles will undergo charge exchange. The degree of ionization and their intermolecular and electrostatic interactions are controlled by varying the pH and salt concentration of the solution in a system. Although the pH can be tuned in experiments, it is hard to model this effect using simulations or theoretical approaches. This is due to the difficulty in treating charge regulation and capturing the cooperative effects in a colloidal suspension with Coulombic interaction. In this work, we analyze a suspension of ionizable colloidal particles via molecular dynamics (MD) simulations, along with Monte Carlo simulations for charge regulation (MC-CR) and derive a phase diagram of the system as a function of pH. It is observed that as pH increases, particles functionalized with acid groups change their arrangement from face-centered cubic (FCC) packing to a disordered state. We attribute these transitions to an increase in the degree of charge polydispersity arising from an increase in pH. Our work shows that charge regulation leads to amorphous solids in colloids when the mean nanoparticle charge is sufficiently high. Graphical abstract

36 MATERIALS SCIENCE↗