Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of A Hazardous Gas Purge System in an EUS Aft Adapter Umbilical Plate Concept

The Fluid Dynamics Branch at the NASA Marshall Space Flight Center actively supports flight programs through a variety of in-line analyses. One area of expertise is the mitigation of hazardous gas mixtures that may form during nominal vehicle fueling operation. After the successful Artemis I launch, design cycles for future iterations of the vehicle are in progress. The subject of this work concerns the analysis of a hazardous gas mitigation system concept for a future upper stage.

Riley Bridges↗
DOE OSTI2023

$\mathrm{DRAGON}$: A multi-GPU orbital-free density functional theory molecular dynamics simulation package for modeling of warm dense matter

As progress in electronic structure theoretical methods is made, ab initio molecular dynamics (MD) based on orbital-free density functional theory (OF-DFT) is becoming increasingly more successful at substituting the traditional, very accurate but computationally costly Kohn–Sham (KS) approach for simulations of matter at the challenging warm dense matter (WDM) regime. However, despite the significant cost alleviation of eliminating the dependence on the KS orbitals, OF-DFT MD runs require ~10 2 to 10 3 CPU cores running for days, or even weeks, for simulations of systems comprised of 10 2 to 10 3 atoms, depending on thermodynamic conditions. Here we present DRAGON, a multi-GPU OF-DFT MD code for fast and efficient simulations of WDM. With a relatively small allocation of resources (4 to 8 GPU devices) it can provide an order of magnitude speedup for simulations containing $\mathscr{O}$(10 4 ) atoms and target systems composed of $\mathscr{O}$(10 5 ) atoms at conditions within the WDM regime, which is currently outside the capabilities of CPU codes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRSDate not supplied

Computational Fluid Dynamics Simulation of the Swordfish Oxygen/Methane Rotating Detonation Rocket Engine

Rotating Detonation Rocket Engines (RDREs) provide several theoretical performance benefits over traditional deflagration combustion rocket engines. RDREs can potentially generate higher specific impulse (ISP) for a given propellant combination enabling significant consumable mass savings for a spaceflight vehicle. Conversely, RDREs are a relatively new technology with a shallower analytical modeling base compared to traditional combustion engine architectures. In this work, we utilize time-accurate, Computational Fluid Dynamics (CFD) modeling to simulate the complex fluid environment inside an RDRE combustion chamber and propellant manifolds. The selected engine of study is the Swordfish RDRE which was developed and tested and at the National Aeronautics and Space Administration (NASA) Marshall Space Flight Center (MSFC). CFD simulations were conducted using the as-tested Swordfish geometry and conditions to facilitate comparison between the analytical predictions and experimental results.

Brian R. Richardson↗
DOE OSTI2016

Molecular dynamics simulations of cesium adsorption on illite nanoparticles

Here we present that the charged surfaces of micaceous minerals, especially illite, regulate the mobility of the major radioisotopes of Cs ( 134 Cs, 135 Cs, 137 Cs) in the geosphere. Despite the long history of Cs adsorption studies, the nature of the illite surface sites remains incompletely understood. To address this problem, we present atomistic simulations of Cs competition with Na for three candidate illite adsorption sites – edge, basal plane, and interlayer. Our simulation results are broadly consistent with affinities and selectivities that have been inferred from surface complexation models. Cation exchange on the basal planes is thermodynamically ideal, but exchange on edge surfaces and within interlayers shows complex, thermodynamically non-ideal behavior. The basal planes are weakly Cs-selective, while edges and interlayers have much higher affinity for Cs. The dynamics of Na-Cs exchange are rapid for both cations on the basal planes, but considerably slower for Cs localized on edge surfaces. Lastly, in addition to new insights into Cs adsorption and exchange with Na on illite, we report the development of a methodology capable of simulating fully-flexible clay mineral nanoparticles with stable edge surfaces using a well-tested interatomic potential model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2001

Textbook Multigrid Efficiency for Computational Fluid Dynamics Simulations

Considerable progress over the past thirty years has been made in the development of large-scale computational fluid dynamics (CFD) solvers for the Euler and Navier-Stokes equations. Computations are used routinely to design the cruise shapes of transport aircraft through complex-geometry simulations involving the solution of 25-100 million equations; in this arena the number of wind-tunnel tests for a new design has been substantially reduced. However, simulations of the entire flight envelope of the vehicle, including maximum lift, buffet onset, flutter, and control effectiveness have not been as successful in eliminating the reliance on wind-tunnel testing. These simulations involve unsteady flows with more separation and stronger shock waves than at cruise. The main reasons limiting further inroads of CFD into the design process are: (1) the reliability of turbulence models; and (2) the time and expense of the numerical simulation. Because of the prohibitive resolution requirements of direct simulations at high Reynolds numbers, transition and turbulence modeling is expected to remain an issue for the near term. The focus of this paper addresses the latter problem by attempting to attain optimal efficiencies in solving the governing equations. Typically current CFD codes based on the use of multigrid acceleration techniques and multistage Runge-Kutta time-stepping schemes are able to converge lift and drag values for cruise configurations within approximately 1000 residual evaluations. An optimally convergent method is defined as having textbook multigrid efficiency (TME), meaning the solutions to the governing system of equations are attained in a computational work which is a small (less than 10) multiple of the operation count in the discretized system of equations (residual equations). In this paper, a distributed relaxation approach to achieving TME for Reynolds-averaged Navier-Stokes (RNAS) equations are discussed along with the foundations that form the basis of this approach. Because the governing equations are a set of coupled nonlinear conservation equations with discontinuities (shocks, slip lines, etc.) and singularities (flow- or grid-induced), the difficulties are many. This paper summarizes recent progress towards the attainment of TME in basic CFD simulations.

Brandt, Achi↗
DOE OSTI2020

Protein Molecular Dynamics Simulations with Approximate QM: What Can We Learn?

Classical force fields are essential for computer simulations of proteins and are typically parameterized to reproduce secondary and tertiary structure of isolated proteins. However, while protein-protein interactions are ubiquitous in nature, they are not considered in parameterization efforts and are far less understood than isolated proteins. A better characterization of intermolecular interactions is widely recognized as a key to revolutionizing drug and therapeutic developments with high-throughput computational screening. Urgently needed is a critical assessment of the performance of modern protein force fields against first-principles electronic structure methods and experiments. In a daring step toward this goal, we here describe a comparison of peptide folding dynamics as predicted by a molecular mechanics force field on the one hand and by an approximate electronic structure quantum mechanical (QM) method based on density-functional tight-binding (DFTB) on the other. We further compare the dynamics from straightforward DFTB simulations with a near-linear scaling version of DFTB for massively parallel computation based on the fragment molecular orbital (FMO-DFTB) method. We illustrate differences between the phenomenology of the folding dynamics from these three methods for a small model peptide, as well as charge polarization and dynamic fluctuations, point out possible correlations and implications for force field developers, and discuss the lessons learned that might become applicable to future predictive high-throughput computer screening for personalized neoantigen cancer therapy.

Irle, Stephan↗
DOE OSTI2021

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Computational Fluid Dynamics Simulation of a Pilot-Scale Waste Glass Melter

The INL has been developing a high-fidelity computational fluid dynamics (CFD) model of a Joule-heated ceramic-lined waste glass melter. The purpose is to better understand the dynamics of heat transfer within the system, in order to provide support for the optimization of the Hanford Waste Treatment Plant (WTP) that will be used for vitrification of legacy tank waste. The simulation presented here is for the DM1200 pilot-scale system, which has been operated by the Vitreous State Laboratory (VSL) over the past several decades. This system is about one-third of the scale of the full size WTP melter. Experimental data from this system is being used to validate computational simulations. The examination of the melter with a fully coupled CFD model should give an accurate assessment of the heat transfer into the reactive cold cap region. Heat is provided to the cold cap region from conduction and convection from the molten glass pool as well as radiation from the hot plenum above the cold cap. The CFD model implements a mathematical model of the cold cap developed for the assessment of the melting rate during the vitrification of tank waste. This sub-model implemented directly into the CFD model combines the mass and enthalpy balance of the cold cap together with the conversion kinetics of the feed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Acceleration of diffusive molecular dynamics simulations through mean field approximation and subcycling time integration

Highlights: • Two numerical methods are presented to accelerate DMD simulations. • The methods greatly reduce the computation time without much loss of accuracy. • The accelerated DMD model is capable of simulating long-term mass transport at atomistic length scales. • Demonstrated for the propagation of an atomically sharp phase boundary over 30 seconds. - Abstract: Diffusive Molecular Dynamics (DMD) is a class of recently developed computational models for the simulation of long-term diffusive mass transport at atomistic length scales. Compared to previous atomistic models, e.g., transition state theory based accelerated molecular dynamics, DMD allows the use of larger time-step sizes, but has a higher computational complexity at each time-step due to the need to solve a nonlinear optimization problem at every time-step. This paper presents two numerical methods to accelerate DMD based simulations. First, we show that when a many-body potential function, e.g., embedded atom method (EAM), is employed, the cost of DMD is dominated by the computation of the mean of the potential function and its derivatives, which are high-dimensional random variables. To reduce the cost, we explore both first- and second-order mean field approximations. Specifically, we show that the first-order approximation, which uses a point estimate to calculate the mean, can reduce the cost by two to three orders of magnitude, but may introduce relatively large error in the solution. We show that adding an approximate second-order correction term can significantly reduce the error without much increase in computational cost. Second, we show that DMD can be significantly accelerated through subcycling time integration, as the cost of integrating the empirical diffusion equation is much lower than that of the optimization solver. To assess the DMD model and the numerical approximation methods, we present two groups of numerical experiments that simulate the diffusion of hydrogen in palladium nanoparticles. In particular, we show that the computational framework is capable of capturing the propagation of an atomically sharp phase boundary over a time window of more than 30 seconds. The effects of the proposed numerical methods on solution accuracy and computation time are also assessed quantitatively.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Using Dynamic Simulation to Evaluate Attemperator Operation in a Natural Gas Combined Cycle With Duct Burners in the Heat Recovery Steam Generator

A generic training simulator of a natural gas combined cycle was modified to match operations at a real plant. The objective was to use the simulator to analyze cycling operations of the plant. Initial operation of the simulator revealed the potential for saturation conditions in the final high pressure superheater as the attemperator tried to control temperature at the superheater outlet during gas turbine loading and unloading. Subsequent plant operational data confirmed simulation results. Multiple simulations were performed during loading and unloading of the gas turbine to determine operational strategies that prevented saturation and increased the approach to saturation temperature. The solutions included changes to the attemperator temperature control setpoints and strategic control of the steam turbine inlet pressure control valve.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2024

Elucidating and contrasting the mechanisms for Mg and Ca sulfate ion-pair formation with multi-level embedded quantum mechanics/molecular dynamics simulations

Here, solutions and minerals containing sulfate (SO 4 2- ), and Ca 2+ and Mg 2+ cations, are ubiquitous throughout the lithosphere and are significant components of seawater, thus presenting a prototypical system for the study of strong electrolytes and crystal nucleation mechanisms. However, despite their relative abundance, key questions remain unanswered about the most fundamental atomic-level steps of their mineralization pathways and aqueous dynamics. Here, we carry out enhanced sampling multi-level molecular dynamics (MD) embedded correlated wavefunction theory simulations to elucidate ion-pairing mechanisms for Mg–SO 4 and Ca–SO 4 in concentrated aqueous solution, accurately capturing effects arising from both structural dynamics and electron exchange–correlation. We predict contact-ion-pair formation to be barrierless and highly exoergic for Ca–SO 4 , in agreement with its minimal solubility, whereas for Mg–SO 4 , solvent-shared and contact ion pairs have similar free energies, qualitatively consistent with its higher solubility. Finally, we demonstrate that brief high-temperature pre-equilibration may be utilized to accelerate convergence of free energies in blue-moon-ensemble enhanced-sampling MD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Finite element for rotor/stator interactive forces in general engine dynamic simulation. Part 1: Development of bearing damper element

A general purpose squeeze-film damper interactive force element was developed, coded into a software package (module) and debugged. This software package was applied to nonliner dynamic analyses of some simple rotor systems. Results for pressure distributions show that the long bearing (end sealed) is a stronger bearing as compared to the short bearing as expected. Results of the nonlinear dynamic analysis, using a four degree of freedom simulation model, showed that the orbit of the rotating shaft increases nonlinearity to fill the bearing clearance as the unbalanced weight increases.

Adams, M. L.↗
NASA NTRS1995

A computational fluid dynamics simulation of the hypersonic flight of the Pegasus(TM) vehicle using an artificial viscosity model and a nonlinear filtering method

The computational fluid dynamics code, PARC3D, is tested to see if its use of non-physical artificial dissipation affects the accuracy of its results. This is accomplished by simulating a shock-laminar boundary layer interaction and several hypersonic flight conditions of the Pegasus(TM) launch vehicle using full artificial dissipation, low artificial dissipation, and the Engquist filter. Before the filter is applied to the PARC3D code, it is validated in one-dimensional and two-dimensional form in a MacCormack scheme against the Riemann and convergent duct problem. For this explicit scheme, the filter shows great improvements in accuracy and computational time as opposed to the nonfiltered solutions. However, for the implicit PARC3D code it is found that the best estimate of the Pegasus experimental heat fluxes and surface pressures is the simulation utilizing low artificial dissipation and no filter. The filter does improve accuracy over the artificially dissipative case but at a computational expense greater than that achieved by the low artificial dissipation case which has no computational time penalty and shows better results. For the shock-boundary layer simulation, the filter does well in terms of accuracy for a strong impingement shock but not as well for weaker shock strengths. Furthermore, for the latter problem the filter reduces the required computational time to convergence by 18.7 percent.

Mendoza, John Cadiz↗
DOE OSTI2018

Net-baryon diffusion in fluid-dynamic simulations of relativistic heavy-ion collisions

A hybrid (hydrodynamics + hadronic transport) theoretical framework is assembled to model the bulk dynamics of relativistic heavy-ion collisions at energies accessible in the Beam Energy Scan (BES) program at the Relativistic Heavy-Ion Collider (RHIC) and the NA61/SHINE experiment at CERN. The system's energy-momentum tensor and net baryon current are evolved according to relativistic hydrodynamics with finite shear viscosity and non-zero net baryon diffusion. Our hydrodynamic description is matched to a hadronic transport model in the dilute region. With this fully integrated theoretical framework, we present a pilot study of the hadronic chemistry, particle spectra, and anisotropic flow. Phenomenological effects of a non-zero net-baryon current and its diffusion on hadronic observables are presented for the first time. Furthermore, the importance of the hadronic transport phase is also investigated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Formation of Dislocations and Twins As a Result of Uniaxial Compression of Magnesium Single Crystals: Molecular Dynamics Simulation

An atomistic simulation of the deformation of ideal magnesium crystal along the [112̄0] crystallographic axis has been performed. The evolution of structural defects under load at T = 300–350 K is considered in detail. It is established that the nucleation of dislocations in an ideal crystal occurs when the stress reaches a level of 0.1G (G is the shear modulus). The acting deformation modes are found to be prismatic slip of a dislocations and {101̄3} twinning. The formation of dislocation networks and dislocation sites in the twinning plane is observed. Some reactions are proposed to describe the dislocation evolution in the (3̄034) plane.

36 MATERIALS SCIENCE↗
NASA NTRS1998

Quantum Molecular Dynamics Simulations of Nanotube Tip Assisted Reactions

In this report we detail the development and application of an efficient quantum molecular dynamics computational algorithm and its application to the nanotube-tip assisted reactions on silicon and diamond surfaces. The calculations shed interesting insights into the microscopic picture of tip surface interactions.

Menon, Madhu↗