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 685 records · Page 38

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Multinucleon structure and dynamics via quantum computing

We propose a framework for computing the structure and dynamics for second-quantized many-nucleon Hamiltonians on quantum computers. We develop an oracle-based Hamiltonian input model that computes the many-nucleon states and nonzero Hamiltonian matrix elements of the many-nucleon system. With our Fock-state based input model, we show how to implement the sparse matrix simulation algorithms to calculate the dynamics of the second-quantized many-nucleon Hamiltonian. Based on the dynamics simulation methods, we also present the methodology for structure calculations of the many-nucleon system. In this work, we provide an explicit circuit design of our input model of the second-quantized Hamiltonian within a direct encoding scheme that maps the occupation of each available single-particle state in the many-nucleon state to the state of specific qubit in a quantum register. Here, we analyze our method and provide the asymptotic cost in computing resources for structure and dynamics calculations of many-nucleon systems. For pedagogical purposes, we demonstrate our input model with two model problems in restricted model spaces.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Performance of Coupled Physics Solvers for Multidisciplinary Hypersonic Flow Simulations on Several Classes of Computer Architectures

The application of hypersonic flow simulation tools to realistic flight scenarios will require the coupling of multiple physical effects to the baseline fluid dynamics. Such multiphysics effects can include the aerooelastic response of the airframe or engine components, dynamic transport of atmospheric particles, the deformation of solid-fluid interfaces that can ablate, pyrolyze, or erode, as well as a host of other processes, all of which are governed by unique sets of physical equations and models. Coupling multiple (and potentially disparate) physics solvers to a robust compressible flow solver poses additional challenges related to the stability, performance and scalability of the combined solver. The choices made during the software design process can therefore lead to a variation in simulation efficiency across different computer architectures. In this paper, we will consider two representative multiphysics hypersonic flow scenarios: the interaction of solid particulates with the flow field created by a hypersonic lifting body and the aerooelastic deformation of a model airframe under high-Mach-number flow conditions. For these simulations we explore the behavior of several hypersonic simulation tools, including Kestrel, FUN3D, US3D, and JENRE multiphysics framework, on several high performance computing systems containing various CPU and GPU architectures.

architecture↗

DOE OSTI2024

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Dynamical simulations of the Muon Campus at Fermilab

The Muon Campus at Fermilab is a system through which muons are delivered to the storage ring of the Muon g-2 Experiment (E989). It consists of a set of 1 km beamlines that transport and prepare a highly polarized muon beam out of secondaries produced downstream a target station. Realistic simulations of this beam delivery system (BDS) using COSY INFINITY, and presented here, contribute to the understanding and characterization of the muon beam production in relation to the statistical and systematic uncertainties of the E989 measurement, intended to be smaller than 0.14 parts per million to achieve the goals of the experiment. The impact of nonlinearities from fringe fields and high-order contributions on the BDS performance are presented, as well as detailed studies of the interactions between secondaries and the beamline elements apertures, particle decay channels, spin dynamics and beamline misalignments.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

Development of Composite Photocatalyst Materials that are Highly Selective for Solar Hydrogen Production and their Evaluation in Z-Scheme Reactor Designs

The key technology gap preventing a vertically stacked dual-bed particle suspension reactor from achieving the DOE MYRD&D ultimate cost target for H 2 production remains the lack of materials in particle form factor that exhibit ≥10% solar-to-H 2 energy conversion (STH) efficiency as a suspension. Therefore, our project goals centered around strategies to increase the STH efficiency by enhancing photophysical properties of perovskite oxide particles including increased visible-light absorption, increased selectivity for electrocatalysis of the H 2 evolution reaction (HER) and the O 2 evolution reaction (OER) through development of ultrathin oxide coatings, correlating composition and structure to function, and improving understanding of multiscale transport and kinetic processes.

08 HYDROGEN↗

NASA NTRS2013

Comparison of Computed and Measured Vortex Evolution for a UH-60A Rotor in Forward Flight

A Computational Fluid Dynamics (CFD) simulation using the Navier-Stokes equations was performed to determine the evolutionary and dynamical characteristics of the vortex flowfield for a highly flexible aeroelastic UH-60A rotor in forward flight. The experimental wake data were acquired using Particle Image Velocimetry (PIV) during a test of the fullscale UH-60A rotor in the National Full-Scale Aerodynamics Complex 40- by 80-Foot Wind Tunnel. The PIV measurements were made in a stationary cross-flow plane at 90 deg rotor azimuth. The CFD simulation was performed using the OVERFLOW CFD solver loosely coupled with the rotorcraft comprehensive code CAMRAD II. Characteristics of vortices captured in the PIV plane from different blades are compared with CFD calculations. The blade airloads were calculated using two different turbulence models. A limited spatial, temporal, and CFD/comprehensive-code coupling sensitivity analysis was performed in order to verify the unsteady helicopter simulations with a moving rotor grid system.

Evolution↗

DOE OSTI2024

Modeling transient edge plasma transport with dynamic recycling

The work presents numerical simulation studies of the role that dynamic plasma recycling on the main wall and divertor target surfaces plays in transient edge plasma transport phenomena, such as edge localized modes (ELMs). The studies are performed by coupling the edge plasma transport code UEDGE [Rognlien et al., J. Nucl. Mater. 196–198, 347 (1992)] and the wall reaction–diffusion transport code FACE [Smirnov et al., Fusion Sci. Technol. 71, 75 (2017)]. The two-dimensional, time-dependent, two-way coupling of the codes, in a realistic tokamak geometry, is accomplished using the Integrated Plasma Simulator framework [Elwasif et al., in 18th Euromicro Conference on Parallel, Distributed and Network-Based Processing (PDP 2010), Pisa, Italy (IEEE, 2010), pp. 419–427] for all modeled material plasma boundaries. The simulations show that dynamic plasma recycling has substantially different characteristics on the main wall and on the divertor plates. It is demonstrated that during an ELM cycle the outer wall can dynamically absorb and release a number of particles comparable to that expelled by the ELM from the core plasma, by far exceeding the dynamic retention capacity of the divertor surfaces. The resulting evolution of the edge and divertor plasma conditions during an ELM cycle is analyzed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Corrosion resistant hot melt adhesive to bind metals

Hot melt adhesives (HMAs) play an important role in many industries, and their demand is expected to grow. HMAs don't require any solvents, and their application results in the formation of strong bonds with the substrate upon cooling within seconds. These properties differentiate them from liquid glues and make them preferable for practical application. Currently, commercial HMAs are used in bonding lightweight materials such as paper, polymers, and cartons, and have limited usage in areas necessitating the bonding of heavier objects like metals. Here, in this study, we report a design and testing of versatile platform comprising an ion-coordinating polymer and ionic fillers for performance optimization and understanding of structure–property relationships, enabling the rational design of HMAs with improved adhesion to metal surfaces. All-atom Molecular dynamics (MD) simulations and various characterization methods are used to elucidate the adhesion mechanism in model composite system containing polyethylene oxide mixed with chemically diverse salt particles. The maximum adhesion strength is found in composites with Al(OH) 3 and FeCl 3 fillers. Interestingly, the presence of Al(OH) 3 also provides a multifunctional anticorrosion property as measured electrochemically using the Tafel method. The discovered path to formulations with improved adhesion to metal surfaces constitutes an important step toward advancing HMAs for use in the structural and semi-structural metal work domain.

Adhesive and cohesive forces↗

DOE OSTI2021

Shock-compressed silicon: Hugoniot and sound speed up to 2100 GPa

In this work, the high-pressure equation of state and isentropic sound speed of fluid silicon were studied using laser-driven shock waves. Principal Hugoniot measurements to 2100 GPa were performed using impedance matching to an α-quartz reference. Sound speed was determined by time correlating imposed shock-velocity perturbations in both the sample (Si) and reference material (α-quartz). A change in shock velocity versus particle velocity (u s –u p ) slope on the fluid silicon principal Hugoniot was detected at 200 GPa. Density functional theory-based quantum molecular dynamics simulations suggest that an increase in ionic coordination and 50% increase in average ionization is coincident with the experimentally observed change in slope.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1994

Momentum transport near a magnetic X line in collisionless reconnection

Plasma dynamics and momentum transport near an X line during time-dependent magnetic reconnection in a collisionless plasma are investigated based on two-dimensional particle simulations. We find that a weakly skewed velocity distribution is formed near the magnetic X line, leading to the presence of off-diagonal elements of the plasma pressure tensor. Let the reconnection electric field be in the y direction. The gradients of the off-diagonal elements of the pressure tensor can provide a transport of the y momentum. During the normal magnetic reconnection, the momentum transport associated with the off-diagonal terms of the pressure tensor mediates a transfer of the y momentum from the region near the X line to regions outside the X line. A period of 'reverse magnetic reconnection,' during which the plasma kinetic energy is converted into magnetic energy, is also observed in the simulation. When reverse reconnection occurs, the gradients of the off-diagonal pressure tensor elements can mediate a transfer of y momentum into the X line. It is found that the inertial term also plays a significant role in the force balance near the magnetic X line. An explanation for the origin of the off-diagonal pressure terms is also given in this paper.

Cai, H. J.↗

NASA NTRS1998

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗

NASA NTRSDate not supplied

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume-Surface Interaction (PSI) between lander engine plumes and landing area regolith poses risks to the landers and space exploration missions through view obscuration and high-energy ejecta impacts from eroded particle clouds and vehicle tilting from plume-induced craters. To address these risks, NASA MSFC Fluid Dynamics Branch (ER42) has developed a cascade of predictive simulation capabilities for PSI, including the recent development of a hybrid CFD/engineering model tool that allows for rapid PSI simulations with a descending/ascending vehicle and an eroding surface. The development of this fully-coupled, moving body/eroding surface simulation capability is detailed along with a demonstration of the tool’s capabilities on the Apollo 12 Lunar Module landing. Predictions of eroded mass flow rates as a function of time are verified against Apollo 12 flight data, from which the developed, viscous erosion model has been calibrated.

Plume Surface Interaction↗

NASA NTRSDate not supplied

A Hybrid CFD/Engineering Model Tool for Lunar Lander Surface Erosion Prediction

Plume Surface Interaction↗

DOE OSTI2020

High-Performance Simulation of Dynamic Hydrologic Exchange and Implications for Surrogate Flow and Reactive Transport Modeling in a Large River Corridor

Hydrologic exchange flows (HEFs) have environmental significance in riverine ecosystems. Key river channel factors that influence the spatial and temporal variations of HEFs include river stage, riverbed morphology, and riverbed hydraulic conductivity. However, their impacts on HEFs were often evaluated independently or on small scales. In this study, we numerically evaluated the combined interactions of these factors on HEFs using a high-performance simulator, PFLOTRAN, for subsurface flow and transport. The model covers 51 square kilometers of a selected river corridor with large sinuosity along the Hanford Reach of the Columbia River in Washington, US. Three years of spatially distributed hourly river stages were applied to the riverbed. Compared to the simulation when riverbed heterogeneity is not ignored, the simulation using homogeneous riverbed conductivity underestimated HEFs, especially upwelling from lateral features, and overestimated the mean residence times derived from particle tracking. To derive a surrogate model for the river corridor, we amended the widely used transient storage model (TSM) for riverine solute study at reach scale with reactions. By treating the whole river corridor as a batch reactor, the temporal changes in the exchange rate coefficient for the TSM were derived from the dynamic residence time estimated from the hourly PFLOTRAN results. The TSM results were evaluated against the effective concentrations in the hyporheic zone calculated from the PFLOTRAN simulations. Our results show that there is potential to parameterize surrogate models such as TSM amended with biogeochemical reactions while incorporating small-scale process understandings and the signature of time-varying streamflow to advance the mechanistic understanding of river corridor processes at reach to watershed scales. However, the assumption of a well-mixed storage zone for TSM should be revisited when redox-sensitive reactions in the storage zones play important roles in river corridor functioning.

Fang, Yilin↗

DOE OSTI2022

Characterizing particle-like charge-migration dynamics with high-order harmonic sideband spectroscopy

Here, we introduce high-order harmonic sideband spectroscopy and show that it can be a robust probe of charge migration (CM) in a halogenated carbon-chain molecule. We simulate both the CM and high-order harmonic generation (HHG) dynamics using ab initio time-dependent density-functional theory. We find that CM dynamics initiated along the molecular backbone induces sidebands in the HHG spectrum driven by a delayed laser pulse that is polarized perpendicular to the molecular axis. Monitoring the spectrum as either the HHG laser frequency or the relative delay is scanned allows for the extraction of detailed information about the time-domain characteristics of the CM process.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1987

Collisional dynamics of perturbed particle disks in the solar system

Investigations of the collisional evolution of particulate disks subject to the gravitational perturbation of a more massive particle orbiting within the disk are underway. Both numerical N-body simulations using a novel collision algorithm and analytical kinetic theory are being employed to extend our understanding of perturbed disks in planetary rings and during the formation of the solar system. Particular problems proposed for investigation are: (1) The development and testing of general criteria for a small moonlet to clear a gap and produce observable morphological features in planetary rings; (2) The development of detailed models of collisional damping of the wavy edges observed on the Encke division of Saturn's A ring; and (3) The determination of the extent of runaway growth of the few largest planetesimals during the early stages of planetary accretion.

Roberts, W. W.↗