Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamical 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 973 records · Page 54

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↗

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↗

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↗

Thermal transport in UO 2 with defects and fission products by molecular dynamics simulations

The importance of the thermal transport in nuclear fuel has motivated a wide range of experimental and modelling studies. In this report, the reduction of thermal transport in UO 2 due to defects and fission products has been investigated using non-equilibrium MD simulations, with two sets of empirical potentials for studying the degregation of UO 2 thermal conductivity including a Buckingham type interatomic potential and a recently developed EAM type interatomic potential. Additional parameters for U 5+ and Zr 4+ in UO 2 have been developed for the EAM potential. The thermal conductivity results from MD simulations are then corrected for the spin-phonon scattering through Callaway model formulations. To validate the modelling results, comparison was made with experimental measurements on single crystal hyper-stoichiometric UO 2+x samples.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

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↗

Accelerated Molecular Dynamics Simulations of Reactive Hydrocarbon Systems

The research activities in this project consisted of four different sub-projects. Three different accelerated dynamics techniques (parallel replica dynamics, hyperdynamics, and temperature-accelerated dynamics) were applied to the modeling of pyrolysis of hydrocarbons. In addition, parallel replica dynamics was applied to modeling of polymerization.

02 PETROLEUM↗

Computational Fluid Dynamics Simulation of Lignocellulosic Biomass Transport in a Compression-Screw Feeder

There is continued interest in lignocellulosic biomass feedstocks for producing power, heat, and transportation liquid fuels in order to increase energy diversification and security. However, there are a few problems associated with the biomass production that must be resolved before biorefineries are considered economical. One of these problematic operations is the process of feeding biomass to reactor systems. Compression-screw feeders have been widely used as a way to fulfill this purpose. However, due to variability of feedstock properties and its complex challenges, production plants operate at less than 50% of their design capacity due to process upsets along with mechanical wear and failure. The objective of this study is to numerically investigate the screw-feeder operation and performance for various biomass physical properties and under different operating conditions. Therefore, the findings from this study will reduce the number of physical experiments needed to improve the plant design.

09 BIOMASS FUELS↗

A complete description of tunnelling using direct quantum dynamics simulation: Salicylaldimine proton transfer

We demonstrate here conclusively that the variational multiconfiguration Gaussian (vMCG) method converges to the grid based full quantum dynamics multiconfiguration time-dependent Hartree result for a tunnelling problem in many dimensions, using the intramolecular proton transfer in salicylaldimine as a model system. The 13-dimensional model potential energy surface was obtained from Hartree Fock energies with the 6-31G* basis set and the expectation value of the flux operator along the transition mode was used as a benchmark characteristic. As well as showing excellent convergence of the vMCG method on the model surface using a local harmonic approximation and a moderate number of basis functions, we show that the direct dynamics version of the vMCG also performs very well, usually needs the same number of Gaussians to converge, and converges to exact results if those are obtained on an accurately fitted surface. Finally, we make an important observation that the width of the Gaussian basis functions must be chosen very carefully to obtain accurate results with the use of the frozen-width approximation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic simulation and techno-economic analysis of a concentrated solar power (CSP) plant hybridized with both thermal energy storage and natural gas

The addition of thermal energy storage and natural gas as a complementary energy source improves the flexibility, reliability, and value of concentrated solar power (CSP) plants. Nevertheless, due to the transient nature of solar energy, transitions from solar-only mode and natural-gas mode to hybrid solar-natural gas mode is quite challenging, especially when the plant is equipped with thermal storage. Thus, it is important to develop proper dynamic modeling and control schemes to accurately simulate such transitions. The objective of this study is to address this subject by demonstrating a dynamic model with reliable control schemes for a highly integrated hybrid parabolic trough-natural gas plant equipped with thermal energy storage. The specific goal is to study the dynamics of adding thermal storage to the hybrid plant. It is found that the developed control schemes assist smooth transitions between different operational modes and effective utilization of thermal storage and natural gas to maintain steady power production and steam mass flow rates under different solar conditions. The results demonstrate that the integration of storage regulates power production by solar energy and natural gas during the day time. It also enables an increase in the solar fraction of the hybrid plant while it causes a small decrease in thermodynamic efficiency. The analysis shows that the hybrid plant with the storage has a substantially lower specific CO 2 emission (0.320 tonne/MWh) than single natural gas plant (0.413 tonne/MWh) although it has a higher levelized cost of electricity (86.32/MWh against 74.92/MWh). The hybrid plant with storage demonstrates a promising potential for reliable and clean production of electricity, although research and development should be conducted to lower its cost.

03 NATURAL GAS↗

Aladyn - Adaptive Neural Network Molecular Dynamics Simulation Code: Computational Materials Mini-Application

This report provides an overview and commands description of the Computational Materials mini-application, Aladyn. Aladyn is a simple molecular dynamics code written in FORTRAN 2008, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of quantum mechanics-based energy calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. The Aladyn code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing many core architectures based on graphic processing units (GPUs). The effort is related to the High Performance Computing Incubator (HPCI) project at NASA Langley Research Center.

Yamakov, Vesselin I.↗

Interactions of aqueous amino acids and proteins with the (110) surface of ZnS in molecular dynamics simulations

The growing usage of nanoparticles of zinc sulfide as quantum dots and biosensors calls for a theoretical assessment of interactions of ZnS with biomolecules. We employ the molecular-dynamics-based umbrella sampling method to determine potentials of mean force for 20 single amino acids near the ZnS (110) surface in aqueous solutions. We find that five amino acids do not bind at all and the binding energy of the remaining amino acids does not exceed 4.3 kJ/mol. Such energies are comparable to those found for ZnO (and to hydrogen bonds in proteins) but the nature of the specificity is different. Cysteine can bind with ZnS in a covalent way, e.g., by forming the disulfide bond with S in the solid. If this effect is included within a model incorporating the Morse potential, then the potential well becomes much deeper—the binding energy is close to 98 kJ/mol. We then consider tryptophan cage, a protein of 20 residues, and characterize its events of adsorption to ZnS. We demonstrate the relevance of interactions between the amino acids in the selection of optimal adsorbed conformations and recognize the key role of cysteine in generation of lasting adsorption. We show that ZnS is more hydrophobic than ZnO and that the density profile of water is quite different than that forming near ZnO—it has only a minor articulation into layers. Furthermore, the first layer of water is disordered and mobile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of a thermal neutron scattering law for polyethylene using molecular dynamics simulations

Polyethylene is an organic polymer that has been used as a neutron moderator and reflector material in various assemblies designed to perform critical and subcritical experiments. Neutrons that are produced within an assembly containing polyethylene will interact with its various atomic components (especially hydrogen) to slow down and eventually reach thermal energies. In this work, the thermal neutron scattering cross sections are generated in support of criticality safety analysis of neutronic assemblies that utilize polyethylene. Polyethylene is composed of ethylene (C{sub 2}H{sub 4}) monomers that are attached in long molecular chains producing the polymer structure. The principle neutron scatterer in this material is hydrogen. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of structural, thermal, and mechanical properties of portlandite through molecular dynamics simulations

The structural, thermal, and mechanical properties of portlandite, the primary solid phase of ordinary hydrated cement paste, are investigated using the molecular dynamics method. To understand the effects of temperature on the structural properties of portlandite, the coefficients of thermal expansion of portlandite are determined in the current study and validated with what reported from the experimental tests. The atomic structure of portlandite equilibrated at various temperatures is then subjected to uniaxial tensile strains in the three orthogonal directions and the stress-strain curves are developed. Based on the obtained results, the effect of the direction of straining on the mechanical properties of portlandite is investigated in detail. Structural damage analysis is performed to reveal the failure mechanisms in different directions. The energies of the fractured surfaces are calculated in different directions and compared to those of the ideal surfaces available in the literature. The key mechanical properties, including tensile strength, Young's modulus, and fracture strain, are extracted from the stress-strain curves. The sensitivity of the obtained mechanical properties to temperature and strain rate is then explored in a systematic way. This leads to valuable information on how the structural and mechanical properties of portlandite are affected under various exposure conditions and loading rates. - Graphical abstract: Fracture mechanism of portlandite under uniaxial strain in the z-direction. - Highlights: • The structural, thermal, and mechanical properties of portlandite are investigated. • The coefficients of thermal expansion are determined. • The stress-strain relationships are studied in three orthogonal directions. • The effects of temperature and strain rate on mechanical properties are examined. • The plastic energy required for fracture in the crystalline structure is reported.

36 MATERIALS SCIENCE↗

Computational fluid dynamics simulations of a binary particle bed in a riser-based carbon stripper for chemical looping combustion

Chemical-looping combustion (CLC) is a next generation combustion technology that shows great promise as a solution for the need of high-efficiency low-cost carbon capture from fossil fueled power plants. In this paper, numerical simulations are conducted of a binary particle bed associated with a coal-direct CLC system consisting of coal (represented by plastic beads) and oxygen carrier particles and validated against an experimental riser-based carbon stripper. The detailed particle dynamics and solid-gas and solid-solid interactions are investigated using the Lagrangian particle-tracking approach known as the discrete element method coupled with the computational fluid dynamics solution for the flow field. The simulation results of the fluidization behavior and the separation ratio of the particles are in excellent agreement with the experiment. In conclusion, a credible simulation of a binary particle bed is of particular importance for understanding the details of the fluidization behavior; the baseline simulation established in this work can be used as a tool for designing and optimizing the performance of such systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Particlization in fluid dynamical simulations of heavy-ion collisions: The IS3D module

The IS3D particlization module simulates the emission of hadrons from heavy-ion collisions via Monte-Carlo sampling of the Cooper–Frye formula which converts fluid dynamical information into local phase-space distributions for hadrons. The code package includes multiple choices for the non-equilibrium correction to these distribution functions: the 14-moment approximation, first-order Chapman–Enskog expansion, and two types of modified equilibrium distributions. This makes it possible to explore to what extent heavy-ion experimental data are sensitive to different choices for δf n , presently the main source of theoretical uncertainty in the particlization stage. Here, we validate our particle sampler with a high degree of precision by generating several million hadron emission events from a longitudinally boost-invariant hypersurface and comparing the event-averaged particle spectra and space–time distributions to the Cooper–Frye formula.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Application and use of Multi Body Dynamic Simulations for In Space Mechanisms

The Thermal and Mechanical Analysis Branch of the Space System Department (ES22) is actively growing and expanding a new mechanical analysis capability. In recent years, there has been an increase in the demand for Multi Body Dynamics (MBD) analysis to support small ESPA ring size spacecraft, cubesats, lunar surface operations and large spacecraft with complex operations or deployables. Requirements can be balanced for both needs to ensure mechanism is verified for launch. Multi-body dynamics is the study of multiple rigid or flexible bodies that are linked together or are in contact and the associated kinematics and dynamics of the system. The typical outputs include performance information, loads, and deflections of mechanical systems. They routinely support the development of robotic mechanisms, deployment mechanisms, landing simulations, control systems, and conops planning.

Multi Body Dynamics Analysis↗

Determining Atomistic SAXS Models of Tri-Ubiquitin Chains from Bayesian Analysis of Accelerated Molecular Dynamics Simulations

Small-angle X-ray scattering (SAXS) has become an increasingly popular technique for characterizing the solution ensemble of flexible biomolecules. However, data resulting from SAXS is typically low-dimensional and is therefore difficult to interpret without additional structural knowledge. In theory, molecular dynamics (MD) trajectories can provide this information, but conventional simulations rarely sample the complete ensemble. Here, we demonstrate that accelerated MD simulations can be used to produce higher quality models in shorter time scales than standard simulations, and we present an iterative Bayesian Monte Carlo method that is able to identify multistate ensembles without overfitting. This methodology is applied to several ubiquitin trimers to demonstrate the effect of linkage type on the solution states of the signaling protein. Here, we observe that the linkage site directly affects the solution flexibility of the trimer and theorize that this difference in plasticity contributes to their disparate roles in vivo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗