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

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Energy Transfer Mechanisms in Large Low-Bandgap Polymers from Time-Resolved Experiments and Nonadiabatic Molecular Dynamics Calculations

Conjugated polymers offer unprecedented chemical tunability for modulating energy transfer in a multitude of infrared light applications. In this work, we use a combination of time-resolved spectroscopic experiments and nonadiabatic molecular dynamics calculations to probe the photochemistry and nonradiative transitions in a recently synthesized narrow bandgap donor–acceptor conjugated polymer based on alternating cyclopentadithiophene and electronegative benzothiadiazole heterocycles. Using large-scale semi-empirical nonadiabatic molecular dynamics, which can treat a large 260-atom hexamer, we calculate an S 5 → S 1 lifetime of 34.75 fs, which is consistent with our time-resolved spectroscopic data. Our simulations suggest that vibronic motions of the central carbons in the cyclopentadithiophene functional groups are predominantly involved in the nonradiative transitions, and the excitation becomes more localized on a monomer fragment over time. The combined use of time-resolved experiments and nonadiabatic molecular dynamics calculations in this work provides mechanistic insight into chemical functionalities that can be tuned to enhance energy transfer in other prospective low-bandgap polymer materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Nanoconfined Grain Boundaries Increase the Conductivity of Polycrystalline Molecular Crystals

Soft-solid molecular crystals consist of crystalline grains and fluid grain boundaries (GBs) that enhance the grain binding and transport of Li + ions between the grains. The total ionic conductivity consists of ion migration in both the grains and GBs. To unravel these contributions in adiponitrile (Adpn):LiPF 6 molecular crystals, the GB volume fraction was varied by changing the size of the crystals and the Adpn:LiPF 6 molar ratio. Molecular dynamics (MD) simulations indicate that ion motion was subdiffusive in the grains and “well-diffusive” in the GBs, with GBs characterized as disordered nanoconfined regions of higher charge carrier concentration (∼1 M) than in saturated Adpn:LiPF 6 solutions (0.04 M), and Li + ions predominantly solvated by cyano groups with few contact ion pairs. The diffusivity in the GBs is at least an order of magnitude higher than that in the crystalline grains. The emergent picture is the grains as a reservoir of ions that migrate to faster-conducting GBs.

Crystals↗

DOE OSTI2025

Autonomous Operations for Advanced Reactors Utilizing Supervisory Control

Automation is a critical tenet of reactor plant operations as reliance on nuclear energy increases. Nuclear power plants require a large workforce which does not scale with output; that is, the cost per megawatt increases as reactor output becomes smaller. The economic viability of advanced reactors, particularly small modular reactors (SMRs) and microreactors, requires a significantly reduced onsite workforce. The logical solution is establishing a systematic process of elimination of reliance on human operators, and to the extent possible, replacing these actions with automated functions. In this paper, we propose a method for such transformation to establish a robust technical basis to enable transition to autonomy. Our method is based on finite state automata (FSA)—also known as finite state machines (FSMs). Relying on this method allows us to exploit the rich set of mathematical proofs available in the field of regular languages. FSA are one of the mathematical tools to model discrete event systems (DES). These properties are applied to produce an automated startup controller for the Massachusetts Institute of Technology Research Reactor (MITR). The startup procedure is captured in terms of discrete changes from one state to another while an independent supervisory control system directs the sequence of states and alerts a human in the event of an abnormal operation. First, the design and behavior of the MITR rod control system were modeled in Simulink. Then, the startup procedure was applied to the rod control system and the DES performed a startup by procedurally withdrawing rods to the subcritical position. The simulation also stops rod motion in response to an uncontrollable event and restarts rod motion once the event has been cleared.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

DOE OSTI2025

Rapid wavefield forecasting for earthquake early warning via deep sequence to sequence learning

We propose a deep learning model, WaveCastNet, to forecast high-dimensional wavefields. WaveCastNet integrates a convolutional long expressive memory architecture into a sequence-to-sequence forecasting framework, enabling it to model long-term dependencies and multiscale patterns in both space and time. By sharing weights across spatial and temporal dimensions, WaveCastNet requires significantly fewer parameters than more resource-intensive models such as transformers, resulting in faster inference times. Crucially, WaveCastNet also generalizes better than transformers to rare and critical seismic scenarios, such as high-magnitude earthquakes. Here, we show the ability of the model to predict the intensity and timing of destructive ground motions in real time, using simulated data from the San Francisco Bay Area. Furthermore, we demonstrate its zero-shot capabilities by evaluating WaveCastNet on real earthquake data. Our approach does not require estimating earthquake magnitudes and epicenters, steps that are prone to error in conventional methods, nor does it rely on empirical ground-motion models, which often fail to capture strongly heterogeneous wave propagation effects.

Geophysics↗

DOE OSTI2026

Sub-Doppler Cooling of a Trapped Ion in a Phase-Stable Polarization Gradient

Trapped ions provide a highly controlled platform for quantum sensors, clocks, simulators, and computers, all of which depend on cooling ions close to their motional ground state. Existing methods like Doppler, resolved sideband, and dark resonance cooling balance trade-offs between the final temperature and cooling rate. A traveling polarization gradient has been shown to cool multiple modes quickly and in parallel, but utilizing a stable polarization gradient can achieve lower ion energies, while also allowing tailorable light-matter interactions in the sub-wavelength regime. In this Letter, we demonstrate cooling of a trapped ion below the Doppler limit using a phase-stable polarization gradient created using trap-integrated photonic devices. At an axial frequency of 2⁢𝜋×1.45 MHz we achieve ⟨𝑛⟩=1.56±0.07 in 150 μ⁢s and cooling rates of ∼0.3 quanta/μ⁢s. Here, we examine ion dynamics under different polarization gradient phases, detunings, and intensities, showing reasonable agreement between experimental results and a multilevel model. Cooling is fast and power efficient, with lower average motional Fock state occupation when compared to simulated operation under the corresponding running wave configuration. Our results demonstrate a well-controlled test bed for studying the dynamics of multilevel atomic systems in a phase-stable polarization gradient.

Atom & ion cooling↗

DOE OSTI2025

Moving lens effect: Simulations, forecasts, and foreground mitigation

The peculiar motion of massive objects across the line of sight imprints a dipolar temperature anisotropy pattern on the cosmic microwave background known as the moving lens effect. This effect provides a unique probe of the transverse components of the peculiar velocity field, but has not yet been detected due to its small size. We implement and validate a stacking estimator for the moving lens signal using a galaxy catalog as a tracer of massive haloes combined with reconstructed velocities from the galaxy number density field. Using simulations, we forecast detection prospects for the moving lens signal from current and upcoming microwave background and galaxy surveys. Here, we demonstrate a new foreground mitigation strategy likely sufficient for current datasets, and discuss various sources of systematic error and noise. Upcoming galaxy surveys will provide high-significance statistical detections of the moving lens effect.

Astrophysical & cosmological simulations↗

DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗

DOE OSTI2025

Respiration Signal Pattern Analysis for Doppler Radar Sensor with Passive Node and Its Application in Occupancy Sensing of a Stationary Subject

Doppler radar node occupancy sensors are promising for applications in smart buildings due to their simple circuits and price advantage compared to quadrature radar sensors. However, single-channel sensitivity limitations may result in low sensitivity and misinterpreted motion rates if the detected subject is at or close to “null” points. We designed and tested a novel method to eliminate such limits, demonstrating that passive nodes can be used to detect a sedentary person regardless of position. This method is based on characteristics of chest motion due to respiration, found via both simulations and experiments based on a sinusoidal model and a more realistic model of cardiorespiratory motion. In addition, respiratory rate variability is considered to distinguish a true human presence from a mechanical target. Sensor node data were collected simultaneously with an infrared camera system, which provided a respiration signal reference, to test the algorithm with 19 human subjects and a mechanical target. The results indicate that a human presence was detected with 100% accuracy and successfully differentiated from a mechanical target in a controlled environment. The developed method can greatly improve the occupancy detection accuracy of single-channel radar-based occupancy sensors and facilitate their adoption in smart building applications.

Song, Chenyan↗

DOE OSTI2025

Review of recent WEC-Sim (v6.1) advanced features

WEC-Sim (Wave Energy Converter SIMulator) is an opensource software for modeling the motions, loads and power generation of wave energy converters. WEC-Sim performs simulations in the time domain using hydrodynamic coefficients calculated by boundary element method (BEM) frequency-domain potential flow solvers such as WAMIT, NEMOH, Capytaine, or Ansys AQWA. WEC-Sim development is ongoing, including various features, applications, and example cases used to demonstrate potential use-cases to meet the needs of the growing marine energy industry. Through input from a broad user base and an extensive team of developers and collaborators, new features of WEC-Sim are developed to expand the software’s use cases and improve overall functionality. Three new WEC-Sim features highlighted in this paper include updating WEC-Sim to be compatible with MoorDyn Version 2, incorporation of second order excitation loads (quadratic transfer functions) and allowing for dynamically changing hydrodynamics.

13 HYDRO ENERGY↗

DOE OSTI2025

Kinetic Plasma Simulation Capabilities in the MOOSE Framework: Verification of Particle-Particle Collisions

High-fidelity simulations of complex plasma systems allow researchers to gain key insights into and understanding of these systems. To facilitate massively parallel high-fidelity plasma simulations, finite-element-based particle-in-cell capabilities are being developed within the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE) based framework called Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER). While SALAMANDER’s primary objective is modeling edge plasmas and plasma-facing components in fusion devices, the particle-in-cell capabilities being developed are general and will support modeling low-temperature plasmas as well. Previously, collisionless magnetostatic simulation capabilities have been verified with the two-stream and Dorey-Guest-Harris instabilities, and single particle motion. Collisions were implemented using the direct simulation Monte Carlo method, and verification of this capability will be presented here several verification problems: relaxation of a randomly initialized gas to a Maxwellian distribution, Fourier heat flow, and comparison of reaction rates to both analytic calculations and those calculated using a multi-term Boltzmann solver.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Parametric amplification of spin-motion coupling in three-dimensional trapped-ion crystals

Three-dimensional (3D) crystals offer a route to scaling up trapped-ion systems for quantum sensing and quantum simulation applications; however, engineering coherent spin-motion couplings and effective spin-spin interactions in large crystals poses technical challenges associated with decoherence and prolonged timescales to generate appreciable entanglement. Here, we explore the possibility of speeding up these interactions in 3D crystals via parametric amplification. For this purpose, we derive a general Hamiltonian for the parametric amplification of spin-motion coupling that is broadly applicable to normal modes with motion transverse to or along the spatial extent of the crystal. Unlike in lower-dimensional crystals, we find that the ability to faithfully (uniformly) amplify the spin-spin interactions in 3D crystals depends on the physical implementation of the spin-motion coupling. We consider the light-shift gate, and the so-called phase-insensitive and phase-sensitive Mølmer-Sørensen (MS) gates, and we find that only the phase-sensitive MS gate can be faithfully amplified in general 3D crystals. We discuss a situation where nonuniform amplification can be advantageous. We also reconsider the effect of counter-rotating terms on parametric amplification and find that they are not as detrimental as previous studies suggest.

Optical parametric oscillators & amplifiers↗

DOE OSTI2026

Molecular origin of negative lithium transference in electrolytes with star-shaped multivalent anions

Large-scale molecular dynamics simulations illustrate that highly correlated cation–anion motion leads to negative t 0+ on the order of 1 in lithium electrolytes with star-shaped multivalent anions. Large multivalent anions have gained increasing attention for their potential to improve lithium transference in electrolytes. We employ large-scale molecular dynamics simulations based on the Onsager transport framework to investigate ion transport in a lithium electrolyte with star-shaped multivalent anions. The simulations show that t 0+, the cation transference number with respect to solvent velocity, is negative over a wide range of concentration. This is consistent with experimental data reported previously. The simulation-based Onsager transport coefficients reveal that the magnitudes of the cation–cation, anion–anion, and cation–anion correlations are comparable, a signature of highly correlated motion in the electrolyte. Examination of the cation solvation environment indicates the presence of strong cation–anion association across the entire concentration range, which leads to negative t 0+ on the order of −1. Both simulation and experiment also show that the maximum value of t 0+ reaches 0 when the cation concentration is c + = 0.4 M. This is the concentration at which the anions begin to spatially overlap, and lithium ions serve as dynamic linkers to balance cation–cation and cation–anion correlations. Our results provide molecular-level insights into the origin of transference in multivalent electrolytes.

Fang, Chao↗

DOE OSTI2024

Emittance preservation for the electron arm in a single PWFA-LC stage using quasi-adiabatic plasma density ramp matching sections

Plasma-based acceleration (PBA) is being considered for a next generation linear collider (LC). In some PBA-LC designs for the electron arm, the extreme beam parameters are expected to trigger background ion motion within the witness beam, which can lead to longitudinally varying nonlinear focusing forces and result in an unacceptable emittance growth of the beam. To mitigate this, we propose to use quasi-adiabatic plasma density ramps as matching sections at the entrance and exit of each stage. We match the witness electron beam to the low density plasma entrance, where the beam initially has a large matched spot size so the ion motion effects are relatively small. As the beam propagates in the plasma density upramp, it is quasi-adiabatically focused, and its distribution maintains a non-Gaussian equilibrium distribution in each longitudinal slice throughout the process, even when severe ion collapse has occurred. This only causes small amounts of slice emittance growth. The phase mixing between slices with different betatron frequencies leads to additional projected emittance growth within the acceleration stage. A density downramp at the exit of an acceleration section can eliminate much of the slice and projected emittance growth as the beam and ion motion adiabatically defocuses and decreases, respectively. Simulation results from QuickPIC with Azimuthal Decomposition show that within a single acceleration stage with a 25 GeV energy gain, this concept can limit the projected emittance growth to only ∼2% for a 25 GeV, 100 nm emittance witness beam and ∼20% for a 100 GeV, 100 nm normalized emittance witness beam. The trade-off between the adiabaticity of the plasma density ramp and the initial ion motion at the entrance for a given length of the plasma density ramp is also discussed.

Physics↗

DOE OSTI2025

Comment on “Brownian motion of droplets induced by thermal noise”

We simulate phase separated fluids using the Cahn-Hillard fluctuating hydrodynamic (CH-FHD) model and measure the statistical properties of capillary waves generated by thermal fluctuations. Our measurements are in good agreement with stochastic lubrication theory and molecular dynamics simulations but differ significantly from recent CH-FHD results by Zhang et al. [Phys. Rev. E 109, 024208 (2024)2470-004510.1103/PhysRevE.109.024208]. Specifically, we find that capillary wave statistics at thermodynamic equilibrium are independent of transport properties, namely viscosity and species diffusion.

Bell, J B↗

DOE OSTI2025

OpenFAST simulation of floating wind turbines with large heading change

Previous versions of OpenFAST, the physics-based wind turbine engineering and design tool developed by the National Renewable Energy Laboratory, were limited to small rotations of the floating platform. This prevented OpenFAST from being used to simulate important events, such as the loss of a mooring line, or specific floater concepts that might experience large platform yaw motion. To overcome this limitation, we modify the structural dynamics and hydrodynamics modules of OpenFAST to allow unrestricted platform yaw motion. We apply the improved version of OpenFAST to simulate the drifting of a floating wind turbine system after the loss of a mooring line. The results, including both global motion and the internal structural loads at selected locations, appear credible and consistent with expectations.

17 WIND ENERGY↗

DOE OSTI2024

Nonlinear simulations of GAEs in NSTX-U

A set of nonlinear simulations has been performed in order to study the nonlinear evolution of unstable global Alfvén eigenmodes in the National Spherical Torus Experiment-Upgrade (NSTX-U). Results of the single toroidal mode number, n, simulations are compared with a full nonlinear simulation (all toroidal harmonics included). In single-n simulations, the conservation of two integrals of motion of a particle in a cyclotron resonance with a monochromatic wave is demonstrated, resulting in a one-dimensional evolution of the particle distribution in (E,μ,pϕ) phase-space. Nonlinear simulations (both single-n and full nonlinear) show a significant redistribution of the resonant fast ions, especially in the pitch parameter. Thus, the changes in the resonant particle's parallel and perpendicular energies can be several times larger than the total particle energy change, with only a small fraction transferring into the excitation of the mode itself. This implies that even a relatively small amplitude mode can significantly modify the beam distribution in the resonant region. For the NSTX-U case considered, the single-n simulation results are close to full nonlinear simulation only for the most unstable mode, in which case the saturation amplitudes and changes in the fast ion distribution are comparable. In contrast, peak amplitudes of subdominant modes in all-n simulations are smaller by a factor of 3–10 compared to single-n runs due to the flattening of the beam ion distribution by the fastest growing mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Mobility of twinning dislocations in copper up to supersonic speeds

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Observation and Numerical Simulation of Cold Ions Energized by EMIC Waves

This is the first report of significant energization (up to 7,000 eV) of low-energy He + ions, which occurred simultaneously with H-band electromagnetic ion cyclotron (EMIC) wave activity, in a direction mostly perpendicular to the ambient magnetic field. The event was detected by the Arase satellite in the dayside plasmatrough region off the magnetic equator on 15 May 2019. The peak energy of the He + flux enhancements is mostly above 1,000 eV. At some interval, the He + ions are energized up to ~7,000 eV. The H-band waves are excited in a frequency band between the local crossover and helium gyrofrequencies and are close to a linear polarization state with weakly left-handed or right-handed polarization. The normal angle of the waves exhibits significant variation between 0° and 80°, indicating a non-parallel propagation. Here, we run a hybrid code with parameters estimated from the Arase observations to examine the He + energization. The simulations show that cold He + ions are energized up to more than 1,000 eV, similar to the spacecraft observations. From the analysis of the simulated wave fields and cold plasma motions, we found that the ratio of the wave frequency to He + gyrofrequency is a primary factor for transverse energization of cold He + ions. As a consequence of the numerical analysis, we suggest that the significant transverse energization of He + ions observed by Arase is attributed to H-band EMIC waves excited near the local helium gyrofrequency.

79 ASTRONOMY AND ASTROPHYSICS↗