Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulator motion”

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 217 records · Page 12

Understanding and Controlling Photoexcited Molecules in Complex Environments (Final Technical Report)

Predicting the outcome of a photochemical reaction is complicated by factors outside the scope of tradition theoretical chemistry. On such small scales, fluctuations of the environment render outcomes stochastic, nonadiabatic dynamics blur the separation between nuclear and electronic motion, and ultrafast relaxation renders typical assumptions of equilibrium inappropriate. This project aimed to elucidate and control nonadiabatic molecular dynamics in condensed phases. We developed accurate, efficient simulation tools to study irreversible molecular processes computational across a range of scales. We applied these methodological advances to study electron transfer under exotic conditions where spin and topology opened new pathways for reactivity. We also brought these ideas to bare on light harvesting systems, where energy and charge transport are coupled, and where long lifetimes can be deleterious to function. These advances helped to establish general principles for chemical efficiency and guide the design of nanoscale energy materials, molecular machines, and related technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theoretical modeling of a bottom-raised oscillating surge wave energy converter structural loadings and power performances

Here, this study presents theoretical formulations to evaluate the fundamental parameters and performance characteristics of a bottom-raised oscillating surge wave energy converter (OSWEC) device. Employing a flat plate assumption and potential flow formulation in elliptical coordinates, closed-form equations for the added mass, radiation damping, and excitation forces/torques in the relevant pitch-pitch and surge-pitch directions of motion are developed and used to calculate the system's response amplitude operator and the forces and moments acting on the foundation. The model is benchmarked against numerical simulations using WAMIT and WEC-Sim, showcasing excellent agreement. The sensitivity of plate thickness on the analytical hydrodynamic solutions is investigated over several thickness-to-width ratios ranging from 1:80 to 1:10. The results show that as the thickness of the benchmark OSWEC increases, the deviation of the analytical hydrodynamic coefficients from the numerical solutions grows from 3% to 25%. Differences in the excitation forces and torques, however, are contained within 12%. While the flat plate assumption is a limitation of the proposed analytical model, the error is within a reasonable margin for use in the design space exploration phase before a higher-fidelity (and thus more computationally expensive) model is employed. A parametric study demonstrates the ability of the analytical model to quickly sweep over a domain of OSWEC dimensions, illustrating the analytical model's utility in the early phases of design.

13 HYDRO ENERGY↗

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

Shear Response of Ionizable Polymer Melts at the Crossover from Ionomers to Polyelectrolytes

Nonlinear shear response of polymers is affected by inherent barriers for diffusive motion, including entanglements and topology. In ionizable polymers, ionic clusters further constrain the intrinsic dynamics of the polymers, significantly enhancing their viscosity. Here, using fully atomistic molecular dynamics simulations, the nonlinear shear response of ionizable polymers is presented, across the transition from the ionomer regime where distinctive clusters dominate the structure to the polyelectrolyte regime where clusters percolate, in polystyrene randomly sulfonated with fractions of SO 3 − groups of f = 0.20 and 0.35, in pristine and tetrahydrofuran (THF) swollen polyelectrolyte melts. For f = 0.20, the ionic clusters first fracture into smaller clusters followed by splitting into individual ionic groups and eventually reform. At higher f, the clusters morph in shape but do not break under high shear. At very high shear rates, all of the chains stretch and recoil rapidly. As the shear rate is reduced, some chains stretch and recoil, while others remain largely unaffected by the shear. Macroscopically, for all systems, the shear viscosity displays initially an elastic response, followed by nonlinear shear stress overshoot and, eventually, a steady state. The evolution of viscosity with time and shear reflects that of the ionizable domains.

cluster chemistry↗

ZFP: A compressed array representation for numerical computations

HPC trends favor algorithms and implementations that reduce data motion relative to FLOPS. We investigate the use of lossy compressed data arrays in place of traditional IEEE floating point arrays to store the primary data of calculations. Simulation is fundamentally an exercise in controlled approximation, and error introduced by finite-precision arithmetic (or lossy compression) is just one of several sources of error that need to be managed to ensure sufficient accuracy in a computed result. We describe ZFP, a compressed numerical format designed for in-memory storage of multidimensional arrays, and summarize theoretical results that demonstrate that the error of repeated lossy compression can be bounded and controlled. Furthermore, we establish a relationship between grid resolution and compression-induced errors and show that, contrary to conventional floating point, ZFP reduces finite-difference errors with finer grids. We present example calculations that demonstrate data reduction by 4x or more with negligible impact on solution accuracy. Our results further demonstrate several orders-of-magnitude increase in accuracy using ZFP over IEEE floating point and Posits for the same storage budget.

Lindstrom, Peter↗

Long-Duration High-Resolution Large-Eddy Simulations of Hurricane Laura (2020) for Infrastructure Design

This dataset provides long-duration large-eddy simulations (LES) of Hurricane Laura (2020), capturing the turbulent boundary-layer wind field during the storm's passage at horizontal grid spacings of 55.55 m and 10 m. Hurricane Laura made landfall as a Category 4 storm in southwestern Louisiana at 06:00 UTC on 27 August 2020, producing widespread damage to coastal infrastructure. The simulations characterize the spatial and temporal variability of the turbulent wind at heights spanning the lowest 300 m of the boundary layer, throughout the approach and passage of the storm's eyewall over a fixed location. To resolve the spatial heterogeneity of wind and turbulence conditions across the storm, seven independent refined LES runs are performed at radial locations spanning the storm relative to its direction of motion. The refinement positions are defined on a normalized radial coordinate r ̂=r/R, where R≈21 km is Laura's radius of maximum wind, and lie between r ̂=-1.5 and r ̂=+1.5 in ∆r ̂=0.5 increments along a line perpendicular to the storm's track at the midpoint of the simulation. Negative r ̂ corresponds to locations south-southwest of the storm center and positive r ̂ to locations north-northeast. The dataset is segmented into seven radial refinement groups.

17 WIND ENERGY↗

Signatures of Antisymmetric Vibrations in the Ultrafast Dynamics of Quadrupolar Dyes

Antisymmetric molecular vibrations are central to ultrafast, nonadiabatic photophysical and photochemical processes such as conical intersection dynamics, Herzberg–Teller couplings and, potentially, singlet fission. Direct spectroscopic identification of such vibrations is, however, challenging, since they are typically Raman inactive and affect optical transitions only weakly. Here, we report experimental signatures of vibronic couplings to a high-frequency antisymmetric vibration in the excited state dynamics of a quasi-quadrupolar molecule by ultrafast two-dimensional electronic spectroscopy (2DES). The early time, sub-50 fs 2DES maps reveal an asymmetric peak pattern with characteristic low-energy cross-peaks. We show that these peaks arise from stimulated emission transitions from an anharmonic, double-minimum excited state potential energy surface formed by vibronic coupling to a high-frequency antisymmetric mode. Simulations based on a phenomenological essential state model support the results. Our findings offer a new approach for identifying antisymmetric vibrations in ultrafast 2DES and track excited state wavepacket motion before relaxation washes out the spectra.

36 MATERIALS SCIENCE↗

Perfect spinfluid: A divergence-type approach

We present a new formulation of nondissipative relativistic spin hydrodynamics that incorporates spin degrees of freedom into the divergence-type theory framework. Due to the divergence-type structure, it is straightforward to enforce nonlinear causality and symmetric hyperbolicity of the equations of motion, ensuring local well-posedness of the initial-value problem and stability of the theory. Furthermore, in a specific realization based on spin kinetic theory, we prove that the equations of motion remain nonlinearly causal and symmetric-hyperbolic to all orders in the spin potential, provided a specific thermodynamic constraint is satisfied. Here, this framework can be applied for numerical simulations to study the dynamics of spin-polarized fluids, such as the quark-gluon plasma in heavy-ion collisions.

Chirality↗

Sub‐5 Ångstrom Porosity Tuning in Calixarene‐Derived Porous Liquids via Supramolecular Complexation Construction

Abstract Sub‐Ångstrom‐level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub‐5 Ångstrom scale is created via sophisticated porosity tuning in calixarene‐derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy‐, X‐ray‐, and neutron‐scattering techniques. The presence of permanent porosity in calixarene‐derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub‐5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. The approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

Li, Errui [Department of Chemistry University of T↗

Sub-5 Ångstrom Porosity Tuning in Calixarene-Derived Porous Liquids via Supramolecular Complexation Construction

Sub-Ångstrom-level porosity engineering, which is appealing in gas separations, has been demonstrated in solid carbon, polymer, and framework materials but rarely achieved in the liquid phase. In this work, a gas molecular sieving effect in the liquid phase at sub-5 Ångstrom scale is created via sophisticated porosity tuning in calixarene-derived porous liquids (PLs). Type II PLs are constructed via supramolecular complexation between the sodium salts of calixarene derivatives and crown ether solvents. The chemical structure variation and assembly behavior of the porous host upon PL construction are monitored by spectroscopy-, X-ray-, and neutron-scattering techniques. The presence of permanent porosity in calixarene-derived PLs is verified by pressure swing gas uptake, altered CO 2 physisorption behavior, and molecular simulations. Sub-5 Ångstrom porosity tuning within the PL phase is achieved by introducing bulky substituted groups on the benzene ring of the calixarene host, which then greatly affects the dynamic motion and transport behavior of CO 2 molecules and the Xe uptake performance. Further, the approach being demonstrated in this work represents a promising pathway to tune and leverage the porosity effect for enhanced gas uptake capacity and selectivity in liquid sorbents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discrete-Element and Material-Point Method (DEM and MPM) Based Solvers for Sustainable Technologies

We present the use of discrete element method (DEM) and material point method (MPM) in three relevant green technology applications that include biomass feedstock handling, lithium-ion battery manufacturing, and high-pressure reverse osmosis. Our open-source DEM and MPM solvers are developed using performance portable grid and particle management library, AMReX, thus enabling superior performance on NVIDIA and AMD GPUs with > 100 million particles. Our DEM solver resolves the motion of individual particles in a granular system and includes a bonded sphere method for modeling non-spherical particles along with Hertzian and liquid bridge-based contact models. We simulate highly variable biomass feedstock flows in large-scale hoppers for biofuel production and electrode calendering in battery manufacturing using DEM. Our simulations predict flow blockage in large scale biomass hoppers and electrode microstructure variations, thus providing valuable information for biofuel and battery manufacturers, respectively. The second half of the talk will be on MPM and its application towards pore resolved simulations of reverse osmosis membranes under compressive loads. We present a validation study of our MPM simulations with membrane microscopy imaging thus providing useful insights on membrane stability under high pressure conditions. We also present a spectral stability analysis of using linear hat, quadratic and cubic spline basis in MPM indicating regions of numerical stability.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

Raman Sideband Cooling of Molecules in an Optical Tweezer Array to the 3D Motional Ground State

Ultracold polar molecules are promising for quantum information processing and searches for physics beyond the standard model. Laser cooling to ultracold temperatures is an established technique for trapped diatomic and triatomic molecules. Further cooling of the molecules to near the motional ground state is crucial for reducing various dephasings in quantum and precision applications. In this work, we demonstrate Raman sideband cooling (RSC) of CaF molecules in optical tweezers to near their motional ground state, with average motional occupation quantum numbers of n ¯ x = 0.16 ( 12 ) , n ¯ y = 0.17 ( 17 ) (radial directions), and n ¯ z = 0.22 ( 16 ) (axial direction), and a 3-D motional-ground-state probability of 54 ± 18 % of the molecules that survive the RSC. This process paves the way to increase molecular coherence times in optical tweezers for robust quantum computation and simulation applications. Published by the American Physical Society 2024

Physics↗

The Fuel Motion Monitoring System at TREAT - Current Status and Future Plans

An important component of the United States Nuclear Fuel Safety Transient Testing Program, the Fuel Motion Monitoring System (FMMS) at Idaho National Laboratory's Transient Reactor Test Facility (TREAT) is fast-neutron hodoscope capable of imaging the location, movement, and relocation of nuclear fuel experiments under simulated transient accident conditions. The FMMS was refurbished in parallel with the TREAT restart program starting in 2014, restoring 96-channels of fast-neutron detection. Since returning to operation in 2017 the FMMS has supported many fuel safety experiments supporting accident tolerant fuel development, light-water reactor safety, space thermal nuclear propulsion fuel development, and advanced reactor research and development. In Phase 2 of the FMMS restoration. work is now under way to expand the FMMS' field-of-view by adding an additional 96 channels of fast neutron detectors to the system's hodoscope, along with an expanded data acquisition system and associated transient timing electronics. An overview of the FMSS system and its fast-neutron detectors will be presented along with examples of current FMMS imaging performance and associated information.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Solvation governs cation transference in glyme-based lithium battery electrolytes

The efficacy of electrochemical systems is governed by the cation transference number, which represents the fraction of current carried by the working ion. Energy is wasted when field-induced motion also drives anions and solvent molecules, decreasing the transference number to near-zero. We present a systematic study of cation transference in a series of electrolytes: tetraglyme (TG), octaglyme (OG), and poly(ethylene oxide) (PEO) mixed with lithium bis(trifluoromethanesulfonyl)imide. In all three electrolytes, starting from the dilute salt concentration limit, the experimentally measured cation transference number decreases with increasing concentration, reaching a minimum between -0.1 and -0.2, before rising back to positive values. Explicit measurements of field-induced species' velocities by electrophoretic nuclear magnetic resonance indicate that negative cation transference numbers in TG and OG electrolytes are dictated by solvation interactions with minimal contribution from anion-cation interactions. Simulation-based solvation structures indicate that OG serves as a bridge between TG and PEO. Multi-charge positive clusters, which are negligible in TG, become increasingly important at higher chain lengths (OG and PEO). As migrating cations drag their solvation shells, this solvation-induced motion is amplified in glyme electrolytes because of covalent interactions between solvating glyme molecules and free glyme molecules.

Im, Julia↗

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↗

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↗

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↗

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↗