Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamics Simulations”

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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

What do far-infrared spectra of solitary water in “water-in-solvent” systems reveal about water’s solvation and dynamics?

Classical molecular dynamics simulations of water in ionic and dipolar solvents were used to interpret the far-infrared (FIR) rotation/libration spectra of “solitary water” in terms of water’s rotational dynamics and interactions with solvents. Seven solvents represented by nonpolarizable all-atom force fields and a series of idealized variable-charge solvents were used to span the range of solvent polarities (hydrogen bonding) studied experimentally. Simulated spectra capture the solvent dependence observed, as well as the relationship between the frequencies of water libration (νL) and OH stretching bands (νOH). In more strongly interacting solvents, simulated νL are ∼20% higher than those of experiment. In all solvents, the simulated spectra are composites of rotational motions about the two axes perpendicular to water’s dipole moment, and the different frequencies of these two motions are responsible for the breadth of the libration band and the bimodal shape observed in halide ionic liquids. Simulations overestimate the separation of these two components in most solvents. The character of water rotational motions changes markedly with solvent polarity, from quasi-free rotation in nonpolar and weakly polar solvents to highly constrained libration in strongly hydrogen bonding environments. The changeover to librational motions dominating the spectrum occurs between solvents such as benzene (νL ∼ 250 cm−1) and acetonitrile (νL ∼ 400 cm−1). For solvents in the latter category, the mean frequency of the experimental FIR band provides a direct measure of mean-squared torques and, therefore, force constants associated with interactions constraining water’s librational motion.

Chemistry↗
DOE OSTI2020

Patchy particle model of hydrated amorphous silica

Glasses corroded in static aqueous conditions typically reach a slow, residual rate of dissolution as a result of the complex interplay between multiple processes. While intrinsically suited to address this problem, Monte Carlo (MC) models developed to date have relied on lattices to describe the amorphous structure of glass and are thus unable to reproduce the residual rate of glass corrosion. A recently developed MC approach based on amorphous structures derived from molecular dynamics simulations offers a solution [Kerisit and Du J. Non-Cryst. Solids 522 (2019) 119601], but it requires a simple model that can rapidly relax the glass-water interface upon each dissolution/condensation event to retain the ability of MC simulations to reach large spatial and temporal scales. To this end, a patchy particle model of hydrated amorphous silica was developed and evaluated against atomistic simulations in this work. Models of water and amorphous silica were first developed separately, and cross terms were then defined to enable simulations of hydrated amorphous silica structures with varying water content. While the nature of its interaction potential prevents the patchy particle model from reproducing simultaneously the structure and dynamics of the systems of interest with a single set of parameters, it described well the connectivity of hydrated amorphous silica structures and the local coordination geometry of individual species. Therefore, this model opens the door to MC simulations of the residual rate of glass corrosion and offers an alternative to atomistic models to investigate the evolution of silica gels over large spatial and temporal scales.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1981

UH-60A Black Hawk engineering simulation program. Volume 1: Mathematical model

A nonlinear mathematical model of the UR-60A Black Hawk helicopter was developed. This mathematical model, which was based on the Sikorsky General Helicopter (Gen Hel) Flight Dynamics Simulation, provides NASA with an engineering simulation for performance and handling qualities evaluations. This mathematical model is total systems definition of the Black Hawk helicopter represented at a uniform level of sophistication considered necessary for handling qualities evaluations. The model is a total force, large angle representation in six rigid body degrees of freedom. Rotor blade flapping, lagging, and hub rotational degrees of freedom are also represented. In addition to the basic helicopter modules, supportive modules were defined for the landing interface, power unit, ground effects, and gust penetration. Information defining the cockpit environment relevant to pilot in the loop simulation is presented.

Howlett, J. J.↗
DOE OSTI2023

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density–density and hydrostatic pressure–pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. Furthermore, these results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Unusual dynamics of tetrahedral liquids caused by the competition between dynamic heterogeneity and structural heterogeneity

Tetrahedral liquids exhibit intriguing thermodynamic and transport properties because of the various ways tetrahedra can be packed and connected. Recently, an unusual temperature dependence of the stretching exponent β in a model tetrahedral liquid ZnCl 2 from T m + 85 K to T m + 35 K has been reported using neutron-spin echo spectroscopy. This discovery stands in sharp contrast to other glass-forming liquids. In this study, we conducted neural network force field driven molecular dynamic simulations of ZnCl 2 . We found a non-monotonic temperature dependence of β from liquid to supercooled liquid temperatures. Further structural decomposition and dynamic analysis suggest that this unusual dynamic behavior is a result of the competition between the decrease in the diversity of tetrahedra motifs (structural heterogeneity) and the increase in glassy dynamic heterogeneity. Furthermore, this result may contribute to new understandings of the structural relaxation of other network liquids.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Orbiter altitude at ALT interface based on ALSES (Approach and Landing Shuttle Engineering Simulation) analyses. Mission planning, mission analysis and software formulation

The details of the orbiter altitude attainable at the approach and landing tests (ALT) interface determined by the approach and landing shuttle engineering simulation (manned) were documented. The analysis culminated in the verification of the trends observed in a similar analysis performed previously on the space vehicle dynamics simulation (unmanned). Altitude variations attributable to pilot steering variability ranged between 492 ft higher to 383 ft lower. The requirement for this parametric analysis is first elaborated. The specifications, assumptions, and analytical approach used to determine the orbiter altitude at the ALT interface are then presented, followed by the results of the analytical approach and the conclusions and recommendations.

Colwell, C. L.↗
NASA NTRS1992

Implementation of an intelligent control system

A laboratory testbed facility which was constructed at NASA LeRC for the development of an Intelligent Control System (ICS) for reusable rocket engines is described. The framework of the ICS consists of a hierarchy of various control and diagnostic functions. The traditional high speed, closed-loop controller resides at the lowest level of the ICS hierarchy. Above this level resides the diagnostic functions which identify engine faults. The ICS top level consists of the coordination function which manages the interaction between an expert system and a traditional control system. The purpose of the testbed is to demonstrate the feasibility of the OCS concept by implementing the ICS as the primary controller in a simulation of the Space Shuttle Main Engine (SSME). The functions of the ICS which are implemented in the testbed are as follows: an SSME dynamic simulation with selected fault mode models, a reconfigurable controller, a neural network for sensor validation, a model-based failure detection algorithm, a rule based failure detection algorithm, a diagnostic expert system, an intelligent coordinator, and a user interface which provides a graphical representation of the event occurring within the testbed. The diverse nature of the ICS has led to the development of a distributed architecture consisting of specialized hardware and software for the implementation of the various functions. This testbed is made up of five different computer systems. These individual computers are discussed along with the schemes used to implement the various ICS components. The communication between computers and the timing and synchronization between components are also addressed.

Simon, D. L.↗
NASA NTRS2013

NanoSail - D Orbital and Attitude Dynamics

NanoSail-D unfurled January 20th, 2011 and successfully demonstrated the deployment and deorbit capability of a solar sail in low Earth orbit. The orbit was strongly perturbed by solar radiation pressure, aerodynamic drag, and oblate gravity which were modeled using STK HPOP. A comparison of the ballistic coefficient history to the orbit parameters exhibits a strong relationship between orbital lighting, the decay rate of the mean semi-major axis and mean eccentricity. A similar comparison of mean solar area using the STK HPOP solar radiation pressure model exhibits a strong correlation of solar radiation pressure to mean eccentricity and mean argument of perigee. NanoSail-D was not actively controlled and had no capability on-board for attitude or orbit determination. To estimate attitude dynamics we created a 3-DOF attitude dynamics simulation that incorporated highly realistic estimates of perturbing forces into NanoSail-D torque models. By comparing the results of this simulation to the orbital behavior and ground observations of NanoSail-D, we conclude that there is a coupling between the orbit and attitude dynamics as well as establish approximate limits on the location of the NanoSail-D solar center of pressure. Both of these observations contribute valuable data for future solar sail designs and missions.

Heaton, Andrew F.↗
NASA NTRS2018

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗
DOE OSTI2025

Self-Assembly of Anisotropic Particles on Curved Surfaces

The surface curvature of membranes, interfaces, and substrates plays a crucial role in shaping the self-assembly of particles adsorbed on these surfaces. However, little is known about the interplay between particle anisotropy and surface curvature and how they couple to alter the free energy landscape of particle assemblies. Using molecular dynamics simulations, we investigate the effect of prescribed curvatures on a quasi-2D assembly of anisotropic patchy particles. By varying curvature and surface coverage, we uncover a rich geometric phase diagram, with curvature inducing ordered structures entirely absent on planar surfaces. Large spatial domains of ordered structures can contain hidden microdomains of orientational textures imprinted by the surface on the assembly. The dynamical landscape is also reshaped by surface curvature, with a glass-like state emerging at modest densities and high curvature. Changes to the symmetry of the surface curvature are found to result in distinct structures, including phases with mesoscale ordering. In conclusion, our findings show that the coupling between surface curvature and particle geometry opens an unexplored space of morphologies and structures that can be exploited for material design.

Fluids↗
DOE OSTI2021

Simulations of the IR and Raman spectra of water confined in amorphous silica slit pores

Water in nano-scale confining environments is a key element in many biological, material, and geological systems. The structure and dynamics of the liquid can be dramatically modified under these conditions. Probing these changes can be challenging, but vibrational spectroscopy has emerged as a powerful tool for investigating their behavior. A critical, evolving component of this approach is a detailed understanding of the connection between spectroscopic features and molecular-level details. In this paper, this issue is addressed by using molecular dynamics simulations to simulate the linear infrared (IR) and Raman spectra for isotopically dilute HOD in D 2 O confined in hydroxylated amorphous silica slit pores. The effect of slit-pore width and hydroxyl density on the silica surface on the vibrational spectra is also investigated. The primary effect of confinement is a blueshift in the frequency of OH groups donating a hydrogen bond to the silica surface. This appears as a slight shift in the total (measurable) spectra but is clearly seen in the distance-based IR and Raman spectra. Analysis indicates that these changes upon confinement are associated with the weaker hydrogen-bond accepting properties of silica oxygens compared to water molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Anomalous phase separation in a correlated electron system: Machine-learning–enabled large-scale kinetic Monte Carlo simulations

Significance Phase separation is crucial to the functionalities of many correlated electron materials with notable examples including colossal magnetoresistance in manganites and high- T c superconductivity in cuprates. However, the nonequilibrium phase-separation dynamics in such systems are poorly understood theoretically, partly because the required multiscale modeling is computationally very demanding. With the aid of machine-learning methods, we have achieved large-scale dynamical simulations in a representative correlated electron system. We observe an unusual relaxation process that is beyond the framework of classical phase-ordering theories. We also uncover a correlation-induced freezing behavior, which could be a generic feature of phase separation in correlated electron systems.

Zhang, Sheng↗
NASA NTRS1987

GTOSS: Generalized Tethered Object Simulation System

GTOSS represents a tether analysis complex which is described by addressing its family of modules. TOSS is a portable software subsystem specifically designed to be introduced into the environment of any existing vehicle dynamics simulation to add the capability of simulating multiple interacting objects (via multiple tethers). These objects may interact with each other as well as with the vehicle into whose environment TOSS is introduced. GTOSS is a stand alone tethered system analysis program, representing an example of TOSS having been married to a host simulation. RTOSS is the Results Data Base (RDB) subsystem designed to archive TOSS simulation results for future display processing. DTOSS is a display post processors designed to utilize the RDB. DTOSS extracts data from the RDB for multi-page printed time history displays. CTOSS is similar to DTOSS, but is designed to create ASCII plot files. The same time history data formats provided for DTOSS (for printing) are available via CTOSS for plotting. How these and other modules interact with each other is discussed.

Lang, David D.↗
DOE OSTI2023

Multiscale simulations reveal TDP-43 molecular-level interactions driving condensation

The RNA-binding protein TDP-43 is associated with mRNA processing and transport from the nucleus to the cytoplasm. TDP-43 localizes in the nucleus as well as accumulating in cytoplasmic condensates such as stress granules. Aggregation and formation of amyloid-like fibrils of cytoplasmic TDP-43 are hallmarks of numerous neurodegenerative diseases, most strikingly present in >90% of amyotrophic lateral sclerosis (ALS) patients. If excessive accumulation of cytoplasmic TDP-43 causes, or is caused by, neurodegeneration is presently not known. In this work, we use molecular dynamics simulations at multiple resolutions to explore TDP-43 self- and cross-interaction dynamics. A full-length molecular model of TDP-43, all 414 amino acids, was constructed from select structures of the protein functional domains (N-terminal domain, and two RNA recognition motifs, RRM1 and RRM2) and modeling of disordered connecting loops and the low complexity glycine-rich C-terminus domain. All-atom CHARMM36m simulations of single TDP-43 proteins served as guides to construct a coarse-grained Martini 3 model of TDP-43. The Martini model and a coarser implicit solvent Cα model, optimized for disordered proteins, were subsequently used to probe TDP-43 interactions; self-interactions from single-chain full-length TDP-43 simulations, cross-interactions from simulations with two proteins and simulations with assemblies of dozens to hundreds of proteins. Our findings illustrate the utility of different modeling scales for accessing TDP-43 molecular-level interactions and suggest that TDP-43 has numerous interaction preferences or patterns, exhibiting an overall strong, but dynamic, association and driving the formation of biomolecular condensates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite

Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.

Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗
DOE OSTI2025

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗
DOE OSTI2021

Quantification of energy and cost reduction from decreasing dissolved oxygen levels in full‐scale water resource recovery facilities

Abstract Aeration systems often lack the efficiency to maintain a desired residual dissolved oxygen (DO) concentration in the tank in part because little consideration is given to the dynamic daily and seasonal loading conditions. Although advanced aeration controllers exist, the majority of plants have DO set points typically based on common practice and literature values rather than site‐specific conditions, which can result in DO set points higher than those necessary to meet treatment objectives. DO set point reduction strategies have primarily been proposed through either static or dynamic simulations. In this study, the substantial improvements associated with DO set point reduction are demonstrated at full scale. A yearlong characterization of full‐scale aeration dynamics captured the effect of diurnal and seasonal fluctuations on oxygen transfer and energy demand and so facilitated the estimation of the potential savings of DO reduction strategies. Full‐scale validation provided direct evidence of DO reduction strategies inducing an overall enhancement of oxygen transfer efficiency along the different bioreactors, while confirming that energy savings as high as 20% were feasible. This study quantifies the influence of oxygen transfer efficiency on operating choices and site‐specific conditions (control strategy, loading conditions, and influent flow variability). Practitioner points We quantified the energy reduction and cost savings associated with a DO reduction in an aeration tank. For each 0.2 mg/L of DO decreased, the average power demand reduction per unit water treated exceeded 17%. Field measurements of dynamic alpha values eliminate the uncertainty in estimating aeration energy and cost savings from DO variations.

Pasini, Federico↗