SEARCH · Search NASA

Results for “motion 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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dynamic response of a freely rotating butterfly valve in the advanced test reactor − dynamic coefficients modeling

Here, in evaluating the water hammer issue pertaining to the primary-coolant-regulating butterfly valve in the Advanced Test Reactor, the dynamic fluid body interaction (DFBI) approach was implemented in the analysis covered in Part I. Although DFBI modeling accurately and simultaneously solved the dynamic motion of the valve’s disk along with the flow field of the surrounding fluid, it shed little light on the reason behind such motion. For Part II, the reacting torque of the fluid on the disk was decomposed into representations of the dynamic coefficients in terms of stiffness, damping, and added mass. These were evaluated via simulations with steady-state static (stiffness), constant angular speed (damping), and variable angular speed (added mass) disks. Substituting the dynamic coefficients into Newton’s second law enabled the response trajectories to be obtained. Stable (by average) and unstable equilibrium positions and thrust tendencies of the valve were determined based on the stiffness coefficient (or static torque), the response amplitude was dampened or enlarged by the damping coefficient (minorly affected by added mass), and the response frequency was altered by the damping and added mass coefficients. Although the dynamic coefficient approach renders slightly different trajectories, due to the averaging effect of the torque in comparison to the DFBI method, the overall trend of the response aligns with the DFBI simulation, thus confirming the conclusion in Part I that a fix to the current butterfly valve is necessary.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

An unstructured body-of-revolution electromagnetic particle-in-cell algorithm with radial perfectly matched layers and dual polarizations

A novel electromagnetic particle-in-cell algorithm has been developed for fully kinetic plasma simulations on unstructured (irregular) meshes in complex body-of-revolution geometries. The algorithm, implemented in the BORPIC++ code, utilizes a set of field scalings and a coordinate mapping, reducing the Maxwell field problem in a cylindrical system to a Cartesian finite element Maxwell solver in the meridian plane. The latter obviates the cylindrical coordinate singularity in the symmetry axis. The choice of an unstructured finite element discretization enhances the geometrical flexibility of the BORPIC++ solver compared to the more traditional finite difference solvers. Symmetries in Maxwell’s equations are explored to decompose the problem into two dual polarization states with isomorphic representations that enable code reuse. The particle-in-cell scatter and gather steps preserve charge conservation at the discrete level. Our previous algorithm (BORPIC+) discretized the E and B field components of TE Φ and TM Φ polarizations on the finite element (primal) mesh. Here, we employ a new field-update scheme. Using the same finite element (primal) mesh, this scheme advances two sets of field components independently: (1) E and B of TE Φ polarized fields, (E z , E ρ , B Φ ) and (2) D and H of TM Φ polarized fields, (D Φ , H z , H ρ ). Since these field updates are not explicitly coupled, the new field solver obviates the coordinate singularity, which otherwise arises at the cylindrical symmetric axis, ρ = 0 when defining the discrete Hodge matrices (generalized finite element mass matrices). Here, a cylindrical perfectly matched layer is implemented as a boundary condition in the radial direction to simulate open space problems, with periodic boundary conditions in the axial direction. We investigate effects of charged particles moving next to the cylindrical perfectly matched layer. We model azimuthal currents arising from rotational motion of charged rings, which produce TMΦ polarized fields. Several numerical examples are provided to illustrate the first application of the algorithm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

The Effect of Core-Hole Shape on Attosecond Valence Electron Dynamics

Rapid X-ray ionization of a core electron is known to induce the attosecond motion of valence electrons; however, the effect of core-hole shape on the triggered dynamics remains relatively unknown. In this work, the sub-four fs response of prototypical functionalized/heterocyclic/polycyclic molecules was simulated using real-time time-dependent density functional theory (RT-TDDFT), a sudden approximation core-hole, and phenomenologically Auger–Meitner (AM) decay. These molecules included fluorobenzene, chlorobenzene, bromobenzene, phenol, thiophenol, pyridine, phosphorus, and azulene. It is observed that the valence electron dynamics are essentially independent of the core-hole created, provided that it is ionized from an inner-shell orbital and not an inner-valence orbital. This has broad implications for free-electron laser studies of X-ray pumped attosecond processes since the flexibility in edge allows for a wide range of experimental modalities, core-holes with longer AM lifetimes, and molecular targets.

aromatic compounds↗

DOE OSTI2024

Machine learning the electric field response of condensed phase systems using perturbed neural network potentials

Abstract The interaction of condensed phase systems with external electric fields is of major importance in a myriad of processes in nature and technology, ranging from the field-directed motion of cells (galvanotaxis), to geochemistry and the formation of ice phases on planets, to field-directed chemical catalysis and energy storage and conversion systems including supercapacitors, batteries and solar cells. Molecular simulation in the presence of electric fields would give important atomistic insight into these processes but applications of the most accurate methods such as ab-initio molecular dynamics (AIMD) are limited in scope by their computational expense. Here we introduce Perturbed Neural Network Potential Molecular Dynamics (PNNP MD) to push back the accessible time and length scales of such simulations. We demonstrate that important dielectric properties of liquid water including the field-induced relaxation dynamics, the dielectric constant and the field-dependent IR spectrum can be machine learned up to surprisingly high field strengths of about 0.2 V Å −1 without loss in accuracy when compared to ab-initio molecular dynamics. This is remarkable because, in contrast to most previous approaches, the two neural networks on which PNNP MD is based are exclusively trained on molecular configurations sampled from zero-field MD simulations, demonstrating that the networks not only interpolate but also reliably extrapolate the field response. PNNP MD is based on rigorous theory yet it is simple, general, modular, and systematically improvable allowing us to obtain atomistic insight into the interaction of a wide range of condensed phase systems with external electric fields.

Science & Technology - Other Topics↗

DOE OSTI2024

Filaments in and between galaxy clusters at low and mid-frequency with the SKA

Context. Understanding the magnetised Universe is a major challenge in modern astrophysics, and cosmic magnetism has been acknowledged as one of the key scientific drivers of the most ambitious radio instrument ever planned, the Square Kilometre Array. Aims. With this work, we aim to investigate the potential of the Square Kilometre Array and its precursors and pathfinders in the study of magnetic fields in galaxy clusters and filaments through diffuse synchrotron radio emission. Galaxy clusters and filaments of the cosmic web are indeed unique laboratories in which to investigate turbulent fluid motions and large-scale magnetic fields in action, and much of what is known about magnetic fields in galaxy clusters comes from sensitive radio observations. Methods. Based on cosmological magneto-hydrodynamic simulations, we predict radio properties (total intensity and polarisation) of a pair of galaxy clusters connected by a cosmic-web filament. Results. We use our theoretical expectations to explore the potential of polarimetric observations to study large-scale structure magnetic fields in the frequency ranges 50–350 MHz and 950–1760 MHz. We also present predictions for galaxy cluster polarimetric observations with the Square Kilometre Array precursors and pathfinders, such as the LOw frequency ARray 2.0 and the MeerKAT+ telescope. Conclusions. Our findings point out that polarisation observations are particularly powerful for the study of large-scale magnetic fields, since they are not significantly affected by confusion noise. The unprecedented sensitivity and spatial resolution of the intermediatefrequency radio telescopes make them the favourite instruments for the study of these sources through polarimetric data, potentially allowing us to understand if the energy density of relativistic electrons is in equipartition with the magnetic field or rather coupled with the thermal gas density. Our results show that low-frequency instruments also represent a precious tool to study diffuse synchrotron emission in total intensity and polarisation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S 2 ). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. Herein we use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S 2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Effect of artificial viscosity on shocked particle-laden flows for staggered grid Lagrangian methods

Abstract Shocked particle-laden flows are important to many natural and industrial processes. When simulating these systems, artificial viscosity is often required to prevent numerical artifacts, such as ringing, from arising in the pressure and density fields. The linear and quadratic coefficients of the artificial viscosity determine the amount of smoothing that occurs in these fields. For particle-laden flows, however, many of the fluid–particle interaction forces, for example, the pressure gradient force and unsteady forces, depend on gradients in the fluid fields. Furthermore, while the shock passes over a particle, these forces can be more dominant than drag. This means that the artificial viscosity coefficients affect how a particle and fluid interact when simulating shocked particle systems. Here this effect is investigated for isolated particles and for a particle curtain using a staggered grid Lagrangian approach. The artificial viscosity coefficients have a significant impact on the maximum force that a fluid imparts to a particle, which is important for determining whether a particle will break up in response to the shock. Furthermore, it is found that the density ratio between the particle and the fluid is important in determining whether the artificial viscosity coefficients have a significant impact on the particle’s motion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Nonreciprocal superfluidlike topological spin transport

We study superfluidlike spin transport facilitated by thermal diffusion of magnetic domain walls, where the positive and negative chiralities of domain walls act as opposite topological charges. The topological charge conservation leads to algebraic decay of the spin current carried by domain walls, allowing for the transport of spin over extended distances. We demonstrate that the presence of the Dzyaloshinskii-Moriya interaction can lead to nonreciprocity in spin flow, thus effectively realizing a spin ratchet. In one scenario, the nonreciprocity arises due to diode-like behavior where the nucleation of domain walls is governed by thermal activation for one direction of spin current and by viscous injection for the other direction of spin current. Here, we confirm our predictions by micromagnetic simulations of domain walls in TmIG nanowire.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

DOE OSTI2026

Three-Qubit Encoding in Ytterbium-171 Atoms for Simulating 1+1D Quantum Chromodynamics

Simulating nuclear matter described by quantum chromodynamics using quantum computers is notoriously inefficient because of the assortment of quark degrees of freedom such as matter/antimatter, flavor, color, and spin. Here, we propose to address this resource efficiency challenge by encoding three qubits within individual ytterbium-171 atoms of a neutral atom quantum processor. The three qubits are encoded in three distinct sectors: an electronic “clock” transition, the spin-1/2 nucleus, and the lowest two motional states in one radial direction of the harmonic trapping potential. We develop a family of composite sideband pulses and demonstrate a universal gate set and readout protocol for this three-qubit system. We then apply it to single-flavor quantum chromodynamics in 1+1D axial gauge for which the three qubits directly represent the occupancy of quarks in the three colors. We show that two atoms are sufficient to simulate both vacuum persistence oscillations and a screened hadron-number transition. We consider resource requirements and connections to error detection/correction. Our work is a step toward resource-efficient digital simulation of nuclear matter and opens new opportunities for versatile qubit encoding in neutral atom quantum processors.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Conceptual study on using Doppler backscattering to measure magnetic pitch angle in tokamak plasmas

We introduce a new approach to measure the magnetic pitch angle profile in tokamak plasmas with Doppler backscattering (DBS), a technique traditionally used for measuring flows and density fluctuations. The DBS signal is maximised when its probe beam's wavevector is perpendicular to the magnetic field at the cutoff location, independent of the density fluctuations [Hillesheim \emph{et al} 2022 \emph{Nucl. Fusion} \textbf{55} 073024]. Hence, if one could isolate this effect, DBS would then yield information about the magnetic pitch angle. By varying the toroidal launch angle, the DBS beam reaches cutoff with different angles with respect to the magnetic field, but with other properties remaining similar. Hence, the toroidal launch angle which gives maximum backscattered power is thus that which is matched to the pitch angle at the cutoff location, enabling inference of the magnetic pitch angle. We performed systematic scans of the DBS toroidal launch angle for repeated DIII-D tokamak discharges. Experimental DBS data from this scan were analysed and combined with Gaussian beam-tracing simulations using the Scotty code [Hall-Chen \emph{et al} 2022 \emph{Plasma Phys. Control. Fusion} \textbf{64} 095002]. The pitch-angle inferred from DBS is consistent with that from magnetics-only and motional-Stark-effect-constrained (MSE) equilibrium reconstruction in the edge. In the core, the pitch angles from DBS and magnetics-only reconstructions differ by one to two degrees, while simultaneous MSE measurements were not available. The uncertainty in these measurements was under a degree; we show that this uncertainty is primarily due to the error in toroidal steering, the number of toroidally separated measurements, and shot-to-shot repeatability. We find that the error of pitch-angle measurements can be reduced by optimising the poloidal launch angle and initial beam properties. Since DBS has high spatial and temporal resolutions, is non-perturbative, does not require neutral beams, and is likely robust to neutron damage of and debris on the first mirrors, using DBS to measure the pitch angle in future fusion energy systems is especially appealing.

beam tracing↗

DOE OSTI2024

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗

DOE OSTI2026

Curvilinear magnetic effects in helicoid nanotubes

We report the fabrication and characterization of ferromagnetic helicoid nanotubes with a tunable surface curvature and chirality. Through combining focused electron beam induced deposition and magnetron sputtering we realize these complex curvilinear magnetic architectures. Electron holography reveals two distinct states: a single domain remanent state, with a magnetic chirality governed by the geometric chirality, and a vortex-anti-vortex domain wall, preferentially forming at areas of lower curvature. Micromagnetic simulations reveal how the geometry, specifically pitch and local surface curvature, affects the stability and formation of magnetic domain walls. Furthermore, we show that the interplay between magnetic and geometric chirality has a strong impact on domain wall motion, whereby a right-handed geometric chirality leads to faster domain wall motion and a left-handed geometry hinders domain wall transport. The results shown here highlight the potential of 3D magnetic architectures in spintronic devices and offer insights into how geometry and curvature affect the magnetization, the distribution of magnetic solitons and the motion of domain walls.

3D Nanomagnetism↗

DOE OSTI2025

Galaxy formation with wave/fuzzy dark matter: The core-halo structure and the solitonic imprint

Dark matter-dominated cores have long been claimed for the well-studied local group dwarf galaxies. More recently, extended stellar halos have been uncovered around several of these dwarfs through deeper imaging and spectroscopy. Such core-halo structures (inner flat core and a characteristic r −3 asymptotic outer halo profile) are not a feature of conventional cold dark matter (CDM). In contrast, smooth and prominent dark matter cores are predicted for wave/fuzzy dark matter (ψDM). The question arises as to what extent the visible stellar profiles should reflect this dark matter core structure. Here we compare cosmological hydrodynamical simulations of CDM, “WDM” (model used as a proxy for ψDM) & ψDM, aiming to predict the stellar profiles for these three DM scenarios. We show that cores surrounded by extended halos are distinguishable for ψDM, where the stellar density is enhanced in the core due to the presence of the relatively dense soliton. Our analysis demonstrates that, in our simulations, a distinctive core-halo structure does not appear in the case of CDM in the DM, gas, or stars. Whereas we do find a core-halo transition for DM, gas, and stars for ψDM, and the scale of this transition is in line with the predicted core radius set by the soliton scale anticipated for the adopted boson mass of 2.5×10 −22 eV. The presence of a core-halo structure in the stellar profile for Galaxy 1 for ψDM is visible for the most massive and the first galaxy to form in the simulation. Clearly, further simulations are needed to establish how strict this possible relationship is between the DM and stellar core-halo profile as a potential observational discriminator. Furthermore, we observe the anticipated asymmetry for ψDM due to the soliton's motion (jumping and random walk), a distinctive characteristic not found in the symmetric distributions of stars in the warm and CDM models.

dark matter↗

DOE OSTI2025

Collective behavior of “flexicles”

In recent years the functionality of synthetic active microparticles has edged even closer to that of their biological counterparts. However, we still lack the understanding needed to recreate at the microscale key features of autonomous behavior exhibited by microorganisms or swarms of macroscopic robots. In this study, we propose a model for a three-dimensional deformable cellular composite particle consisting of self-propelled rod-shaped colloids confined within a flexible vesicle—representing a superstructure we call a “flexicle” that couples particle deformation to the internal dynamics of the internal active components. Using molecular dynamics simulations, we investigate the collective behavior of dense systems composed of many flexicles. We show that individual flexicles exhibit shape changes upon collisions with other flexicles that lead to rearrangements of the internal active rods, which slows flexicle motion. This shape deformability gives rise to a diverse set of motility-induced phase separation phenomena and the spontaneous flow of flexicles reminiscent of the migration of cells in dense tissues. Our findings establish a foundation for designing responsive, cell-like active particles and developing strategies for controlling swarm migration and other autonomous swarm behaviors at cellular and colloidal scales.

Science & Technology - Other Topics↗

DOE OSTI2024

Exploring Li-Ion Transport Properties of Li 3 TiCl 6 : A Machine Learning Molecular Dynamics Study

We performed large-scale molecular dynamics simulations based on a machine-learning force field (MLFF) to investigate the Li-ion transport mechanism in cation-disordered Li 3 TiCl 6 cathode at six different temperatures, ranging from 25°C to 100°C. In this work, deep neural network method and data generated by ab − initio molecular dynamics (AIMD) simulations were deployed to build a high-fidelity MLFF. Radial distribution functions, Li-ion mean square displacements (MSD), diffusion coefficients, ionic conductivity, activation energy, and crystallographic direction-dependent migration barriers were calculated and compared with corresponding AIMD and experimental data to benchmark the accuracy of the MLFF. From MSD analysis, we captured both the self and distinct parts of Li-ion dynamics. The latter reveals that the Li-ions are involved in anti-correlation motion that was rarely reported for solid-state materials. Similarly, the self and distinct parts of Li-ion dynamics were used to determine Haven’s ratio to describe the Li-ion transport mechanism in Li 3 TiCl 6 . Obtained trajectory from molecular dynamics infers that the Li-ion transportation is mainly through interstitial hopping which was confirmed by intra- and inter-layer Li-ion displacement with respect to simulation time. Ionic conductivity (1.06 mS/cm) and activation energy (0.29eV) calculated by our simulation are highly comparable with that of experimental values. Overall, the combination of machine-learning methods and AIMD simulations explains the intricate electrochemical properties of the Li 3 TiCl 6 cathode with remarkably reduced computational time. Thus, our work strongly suggests that the deep neural network-based MLFF could be a promising method for large-scale complex materials.

Selvaraj, Selva Chandrasekaran (ORCID:000000029023↗

DOE OSTI2026

Ultrafast Proton-Coupled Electron Transfer Reactions at Unit Quantum Yield

Proton-coupled electron transfer (PCET) reactions mediate essential biological energy-conversion processes, yet dynamical experiments that distinguish concerted electron/proton motion from sequential electron-transfer/proton-transfer pathways are unprecedented on ultrafast time scales. Here we report a tightly coupled donor–spacer–acceptor assembly, NDI-Sp-PhOH/BI, that enables direct interrogation of light-triggered and ultrafast PCET dynamics. Photoexcitation of the naphthalene diimide chromophore drives charge separation forming NDI radical anion and benzimidazolium phenoxy radical (PhO●/BIH+) products with a 440 fs time constant, followed by thermal charge recombination on a 5 ps time scale. Transient visible/near-infrared and mid-infrared spectroscopies independently track electron and proton motion and reveal coincident formation and decay of the charge-separated PCET products, excluding detectable phenol radical cation or phenoxide intermediates expected for sequential electron transfer-proton transfer (ET-PT) pathways. Surface-hopping quasiclassical trajectory simulations demonstrate rapid proton relocation coupled to evolution from locally excited to charge-transfer states, congruent with a concerted mechanism. These combined spectroscopic and computational results establish NDI-Sp-PhOH/BI as an archetype for ultrafast, concerted PCET occurring with unit quantum yield; further, because this platform provides distinct electronic and vibrational spectroscopic handles in its ground, electronically excited, and charge-separated states, it provides new opportunities to characterize mechanistic crossover in biomimetic PCET systems.

10 SYNTHETIC FUELS↗

DOE OSTI2025

Plasma instabilities dominate radioactive transients magnetic fields: the self-confinement of leptons in Type Ia and core-collapse supernovae, and kilonovae

The light curves of radioactive transients, such as supernovae and kilonovae, are powered by the decay of radioisotopes, which release high-energy leptons through $\beta ^+$ and $\beta ^-$ decays. These leptons deposit energy into the expanding ejecta. As the ejecta density decreases during expansion, the plasma becomes collisionless, with particle motion governed by electromagnetic forces. In such environments, strong or turbulent magnetic fields are thought to confine particles, though the origin of these fields and the confinement mechanism have remained unclear. Using fully kinetic particle-in-cell (PIC) simulations, we demonstrate that plasma instabilities can naturally confine high-energy leptons. These leptons generate magnetic fields through plasma streaming instabilities, even in the absence of pre-existing fields. The self-generated magnetic fields slow lepton diffusion, enabling confinement, and transferring energy to thermal electrons and ions. Our results naturally explain the positron trapping inferred from late-time observations of thermonuclear and core-collapse supernovae. Furthermore, they suggest potential implications for electron dynamics in the ejecta of kilonovae. We also estimate synchrotron radio luminosities from positrons for Type Ia supernovae and find that such emission could only be detectable with next-generation radio observatories from a Galactic or local-group supernova in an environment without any circumstellar material.

instabilities↗