Search NASA⌕ Search

SEARCH · Search NASA

Results for “motion simulators”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Dynamic Trajectory Fit to Multisensor Fireball Observations

Meteorites with known orbital origins are key to our understanding of solar system formation and the source of life on Earth. Fireball networks have been developed globally in a unified effort to record and ultimately retrieve these cosmic samples. However, the accuracy of the determined orbit and the likelihood of meteorite recovery depend directly on the accuracy of the chosen meteoroid triangulation method. There are three leading techniques for meteoroid triangulation discussed in the literature: the method of planes, the straight-line least-squares method, and the multiparameter fit method. Here we describe an alternative method to meteoroid triangulation, called the dynamic trajectory fit. This approach uses the meteoroid’s 3D dynamic equations of motion to fit a realistic trajectory directly to multisensor line-of-sight observations. This method has the ability to resolve fragmentation events, fit systematic observatory timing offsets, and determine mass estimates of the meteoroid along its observable trajectory. Through a comprehensive Monte Carlo analysis of over 100,000 trajectory simulations, we find this new method to more accurately estimate meteoroid trajectories of slow entry events (<25 km s{sup −1}) and events observed from low convergence angles (<10°) compared to existing meteoroid triangulation techniques. Additionally, we triangulate an observed fireball event with visible fragmentation using the various triangulation methods to show that the proposed dynamic trajectory fit implementing fragmentation to best match the captured multisensor line-of-sight data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Clarifying the definition of ‘transonic’ screw dislocations

We report a number of recent Molecular Dynamics (MD) simulations have demonstrated that screw dislocations in face centred cubic (fcc) metals can achieve stable steady state motion above the lowest shear wave speed (v shear ) which is parallel to their direction of motion (often referred to as transonic motion). This is in direct contrast to classical continuum analyses which predict a divergence in the elastic energy of the host material at a crystal geometry dependent ‘critical’ velocity v crit . Within this work, we first demonstrate through analytic analyses that the elastic energy of the host material diverges at a dislocation velocity (v crit ) which is greater than v shear , i.e. v crit >v shear . We argue that it is this latter derived velocity (v crit ) which separates ‘subsonic’ and ‘supersonic’ regimes of dislocation motion in the analytic solution. In addition to our analyses, we also present a comprehensive suite of MD simulation results of steady state screw dislocation motion for a range of stresses and several cubic metals at both cryogenic and room temperatures. At room temperature, both our independent MD simulations and the earlier works find stable screw dislocation motion only below our derived v crit . Nonetheless, in real-world polycrystalline materials v crit cannot be interpreted as a hard limit for subsonic dislocation motion. In fact, at very low temperatures our MD simulations of Cu at 10 Kelvin confirm a recent claim in the literature that true ‘supersonic’ screw dislocations with dislocation velocities v>v crit are possible at very low temperatures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Multiscale CFD simulation of biomass fast pyrolysis with a machine learning derived intra-particle model and detailed pyrolysis kinetics

Coupling particle and reactor scale models is as essential as reactor fluid dynamics and particle motion for accurate Computational Fluid Dynamic (CFD) simulations of biomass fast pyrolysis reactors due to intraparticle heat transfer and chemical reactions controlling conversion time and product distributions. Direct online coupling of a particle model with a reactor model is computationally expensive, while offline coupling is case-dependent. In this research, solutions from a series of particle pyrolysis simulations were regressed with Artificial Neural Network (ANN). Furthermoer, this machine learning-derived model predicted the same temperature and conversion profiles compared with particle resolved simulation while the isothermal approach overpredicted the temperature by 130 K and underpredicted the conversion time by 30 s. The ANN model was then integrated into CFD simulations of fluidized bed biomass fast pyrolysis with varied feedstocks via coupling PyTorch and MFiX. The averaged error of simulation predicted bio-oil yields with four feedstocks is 6.4%. This multi-scale approach provides an efficient tool for the coupled particle and reactor scale simulations of biomass pyrolysis.

09 BIOMASS FUELS↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
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 OSTI2022

Unified many-worlds browsing of arbitrary physics-based animations

Manually tuning physics-based animation parameters to explore a simulation outcome space or achieve desired motion outcomes can be notoriously tedious. This problem has motivated many sophisticated and specialized optimization-based methods for fine-grained (keyframe) control, each of which are typically limited to specific animation phenomena, usually complicated, and, unfortunately, not widely used. In this paper, we propose Unified Many-Worlds Browsing (UMWB), a practical method for sample-level control and exploration of physics-based animations. Our approach supports browsing of large simulation ensembles of arbitrary animation phenomena by using a unified volumetric WORLDPACK representation based on spatiotemporally compressed voxel data associated with geometric occupancy and other low-fidelity animation state. Beyond memory reduction, the WORLDPACK representation also enables unified query support for interactive browsing: it provides fast evaluation of approximate spatiotemporal queries, such as occupancy tests that find ensemble samples ("worlds") where material is either IN or NOT IN a user-specified spacetime region. WORLDPACKS also support real-time hardware-accelerated voxel rendering by exploiting the spatially hierarchical and temporal RLE raster data structure. Our UMWB implementation supports interactive browsing (and offline refinement) of ensembles containing thousands of simulation samples, and fast spatiotemporal queries and ranking. We show UMWB results using a wide variety of physics-based animation phenomena---not just JELL-O ® .

Computer Science↗
DOE OSTI2022

Directed motion of liquid crystal skyrmions with oscillating fields

Abstract Using continuum simulations, we show that under a sinusoidal electric field, liquid crystal skyrmions undergo periodic shape oscillations which produce controlled directed motion. The speed of the skyrmion is non-monotonic in the frequency of the applied field, and exhibits multiple reversals of the motion as a function of changing frequency. We map out the dynamical regime diagram of the forward and reverse motion for two superimposed ac driving frequencies, and show that the reversals and directed motion can occur even when only a single ac driving frequency is present. Using pulsed ac driving, we demonstrate that the motion arises due to an asymmetry in the relaxation times of the skyrmion shape. We discuss the connection between our results and ratchet effects observed in systems without asymmetric substrates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

Designing, synthesizing, and modeling active fluids

We review recent advances in the design, synthesis, and modeling of active fluids. Active fluids have been at the center of many technological innovations and theoretical advances over the past two decades. Research on this new class of fluids has been inspired by the fascinating and remarkably efficient strategies that biological systems employ, leading to the development of biomimetic nano- and micro-machines and swimmers. The review encompasses active fluids on both the nano- and micro-scale. We start with examples of biological active systems before we discuss how experimentalists leverage novel propulsion mechanisms to power nano- and micro-machines. Furthermore, we then examine how the study of these far-from-equilibrium systems has prompted the development of new simulation methods and theoretical models in nonequilibrium physics to account for their mechanical, thermodynamic, and emergent properties. Recent advances in the field have paved the way for the design, synthesis, and modeling of autonomous systems at the nano- and micro-scale and opened the door to the development of soft matter robotics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗
DOE OSTI2022

The driving mode of shock-driven turbulence

Turbulence in the interstellar medium (ISM) is crucial in the process of star formation. Shocks produced by supernova explosions, jets, radiation from massive stars, or galactic spiral-arm dynamics are amongst the most common drivers of turbulence in the ISM. However, it is not fully understood how shocks drive turbulence, in particular whether shock driving is a more solenoidal (rotational, divergence-free) or a more compressive (potential, curl-free) mode of driving turbulence. The mode of turbulence driving has profound consequences for star formation, with compressive driving producing three times larger density dispersion, and an order of magnitude higher star formation rate than solenoidal driving. Here, we use hydrodynamical simulations of a shock inducing turbulent motions in a structured, multiphase medium. This is done in the context of a laser-induced shock, propagating into a foam material, in preparation for an experiment to be performed at the National Ignition Facility (NIF). Specifically, we analyse the density and velocity distributions in the shocked turbulent medium, and measure the turbulence driving parameter $b=(\sigma _{\rho /\langle \rho \rangle }^{2\Gamma }-1)^{1/2}(1-\sigma _{\rho /\langle \rho \rangle }^{-2})^{-1/2}\mathcal {M}^{-1}\Gamma ^{-1/2}$, with the density dispersion σρ/<ρ>, the turbulent Mach number $\mathcal {M}$, and the polytropic exponent Γ. Purely solenoidal and purely compressive driving correspond to b ~ 1/3 and b ~ 1, respectively. As a result, using simulations in which a shock is driven into a multiphase medium with structures of different sizes and Γ < 1, we find b ~ 1 for all cases, showing that shock-driven turbulence is consistent with strongly compressive driving.

79 ASTRONOMY AND ASTROPHYSICS↗
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 OSTI2023

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Real-space inversion and super-resolution of ultrafast scattering

Ultrafast scattering using x-rays or electrons is an emerging method to obtain structure dynamics at the atomic lengthscales and timescales. However, directly resolving in real-space atomic motions is inherently limited by the finite detector range and the probe energy. So, as a result, the time-resolved signal interpretation is mostly done in reciprocal space and relies on modeling and simulations of specific structures and processes. Here, we introduce a model-free approach to directly resolve scattering signals in real space, surpassing the diffraction limit, using scattering kernels and signal priors that naturally arise from the measurement constraints. We demonstrate the approach on simulated and experimental data, recover multiple atomic motions at sub-angstrom resolutions, and discuss the recovery accuracy and resolution limits versus signal fidelity. The approach offers a robust path to obtain high-resolution real-space information of atomic-scale structure dynamics using current time-resolved x-ray or electron scattering sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Ion Conductivity and Correlations in Model Salt-Doped Polymers: Effects of Interaction Strength and Concentration

Correlated anion and cation motion can greatly reduce the overall ion conductivity in electrolytes versus the ideal conductivity calculated based on the diffusion constants alone. Using coarse-grained molecular dynamics simulations, we calculate conductivity and the degree of uncorrelated ion motion in salt-doped homopolymers and block copolymers as a function of concentration and interaction strengths. Calculating conductivity from ion mobility under an applied electric field increases accuracy versus the typical use of fluctuation dissipation relationships in equilibrium simulations. In typical electrolytes, correlation in cation-anion motion is often expected to be reduced at low ion concentrations. Yet, for these polymer electrolytes with strong ion-polymer and ion-ion interactions, we find cation-anion motion is more correlated at lower concentrations when other variables are held constant. We show this phenomenon is related to the slower ion cluster relaxation rate at low concentrations rather than the static spatial state of ion aggregation or the fraction of free ions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗