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 1,693 records · Page 94

Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Iron Dissolution from Goethite (alpha-FeOOH) Surfaces in Water by Ab Initio Enhanced Free-Energy Simulations

Dissolution of redox-active metal oxides plays a key role in a variety of phenomena, including (photo)electrocatalysis, degradation of battery materials, corrosion of metal oxides, and biogeochemical cycling of metals in natural environments. Despite its widespread significance, mechanisms of metal-oxide dissolution remain poorly understood on the atomistic level. This study is aimed at elucidating the long-standing problem of iron dissolution from Fe(III)-oxide, a complex process involving coupled hydrolysis, surface protonation, electron transfer, and metal oxygen bond cleavage. We examine the case of goethite (alpha-FeOOH), a representative phase, bearing structural similarities with many other metal (hydr)oxides. By employing quantum molecular dynamics simulations (metadynamics combined with the Blue Moon ensemble approach), we unveil the mechanistic pathways and rates of both nonreductive and reductive dissolution of iron from the (110) and (021) goethite facets in aqueous solutions at room temperature. Our simulations reveal the interplay between concerted internal (structural) and external (from solution) protonations as essential for breaking Fe-O bonds as well as for stabilizing intermediate configurations of dissolving Fe. We demonstrate specifically how Fe(III) reduction to Fe(II) yields higher dissolution rates than the proton-mediated pathway, whereas the most rapid dissolution is expected for these two processes combined, in agreement with experiments.

Klyukin, Konstantin↗

Atomic scale understanding of organic anion separations using ion-exchange resins

A combination of ab initio and classical molecular dynamic simulations was used to explore the adsorption/desorption and diffusion characteristics of ion-exchange resins for extraction of organic anions in electrically-driven separation processes. We considered two classes of carboxylate mixtures that are commonly encountered in bioprocessing separations: a short chain fatty acid mixture (acetate/butyrate) and an aromatic mixture (ferulate/coumarate). The suitability of several resin materials including PFC100E, IRC86, PFA444 and IRA67 was interrogated. The decomposition of the interaction energies by the symmetry-adapted perturbation theory and the classical molecular dynamic simulations of organic diffusion together reveal that the geometries of the organic anions and the functional groups of the resins, as well as their Columbic interactions, are the controlling factors in the diffusion process of the organic compounds in these resins. Additionally, classical simulations also show that modifying the functionality of resin beads, the magnitude of electric fields, and the ratio of organic mixtures may provide an effective way to control the diffusion rate and obtain selective separation of these organic mixtures. Finally, a general suggestion for favorable separation conditions is summarized.

36 MATERIALS SCIENCE↗

Bastnaesite flotation chemistry issues associated with alkyl phosphate collectors

The hydrophobicity, selectivity, and adsorption state in bastnaesite flotation using alkyl phosphates were examined by contact angle and zeta potential measurements, microflotation and microcalorimetry experiments, as well as density functional theory calculations and molecular dynamics simulations. Here the contact angle and microflotation results show that bastnaesite flotation with alkyl phosphates (lauryl phosphate and 2-ethylhexyl phosphate) is accomplished at a remarkably lower collector concentration. Selective flotation of bastnaesite from calcite with 2-ethylhexyl phosphate increased the bastnaesite grade from 52% to 95% without sacrificing recovery when compared to lauryl phosphate. As expected, the results from density functional theory calculations, together with the microcalorimetry and molecular dynamics simulation results, suggest that the hydrophobic surface state is due to a chemisorbed monolayer of the alkyl phosphate and that the lower reactivity of 2-ethylhexyl phosphate accounts for higher selectivity when compared to results for lauryl phosphate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dehydration-Driven Ion Aggregation and the Onset of Gelation in ZnCl2 Solution

A minimal model of ionic aggregation in concentrated ZnCl2 is developed, guided by molecular dynamics simulations with a machine-learned potential. It explicitly incorporates solvent-site depletion, correlated chloride binding, and allows for loops within Zn–Cl clusters. Dehydration drives two coordination-controlled structural changes: a crossover at Z=2 from predominantly isolated Zn-containing species to Cl-bridged clusters and the onset of gelation near Z≈3. The model reproduces the concentration-dependent trends observed in the molecular dynamics simulations with two fitted parameters, while the cluster-size distribution at the highest concentration is consistent with three-dimensional percolation scaling over the accessible range of cluster sizes.

25 ENERGY STORAGE↗

Coarse-Grained Hydrodynamics from Correlation Functions

This paper will describe a formalism for using correlation functions between different grid cells as the basis for determining coarse-grained hydrodynamic equations for modeling the behavior of mesoscopic fluid systems. Configurations from a molecular dynamics simulation or other atomistic simulation are projected onto basis functions representing grid cells in a continuum hydrodynamic simulation. Equilibrium correlation functions between different grid cells are evaluated from the molecular simulation and used to determine the evolution operator for the coarse-grained hydrodynamic system. The formalism is demonstrated on a discrete particle simulation of diffusion with a spatially dependent diffusion coefficient. Correlation functions are calculated from the particle simulation and the spatially varying diffusion coefficient is recovered using a fitting procedure.

HPC↗

A numerical study of bubble and spike velocities in shock-driven liquid metals

In this study, we use detailed continuum hydrodynamics and molecular dynamics simulations to investigate the dynamics of ejecta that are initialized with large amplitude perturbations and non-sinusoidal shapes. Insights from the simulations are used to suggest a modified expression for the velocity associated with ejected spike structures, whereas a recently suggested model explains the observed bubble velocities. Specifically, we find the asymptotic bubble velocity prediction given by Mikaelian is in excellent agreement with the simulations, when a nonlinear correction for finite amplitudes is used in that model. In contrast, existing models can overpredict observed spike velocities if they do not include the modification of the initial spike growth rates due to nonlinearities. Instead, we find that when potential flow models are corrected with a suitable nonlinear prefactor, this leads to predictions in close agreement with our simulation data. We also propose a simple empirical expression for the nonlinear correction for spike velocities which is able to reproduce results from our simulations and published experimental and simulation data over a wide range of initial conditions and Mach numbers. We discuss extensions of these models to initial interfaces with arbitrary shapes. In particular, for non-sinusoidal shapes, the bubble and spike velocities are still predicted by these models provided we use an effective wavelength λ eff which is the wavelength of an equivalent sinusoid that has the same missing area. Lastly, the issues of nonlinearity, non-standard shapes and shock Mach number addressed in this work are relevant to recent experimental campaigns involving twice-shocked targets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving the reliability of machine learned potentials for modeling inhomogeneous liquids

The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Material Process Modeling at the Y-12 National Security Complex

Dynamic simulation modeling is used at Y-12 to evaluate and forecast nuclear material inventories and production capacities to ensure that future supply can meet mission demand. Model outputs are analyzed by numerous Y-12 organizations and programs and coordinated with NNSA’s Office of Secondary Stage Production Modernization. Data-driven decisions for both short-term and long-term strategic planning for Y-12 mission execution are informed by the model. Dynamic simulation modeling capabilities for Y-12 nuclear material production continue to be expanded and refined.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Structure of cyano-anion ionic liquids: X-ray scattering and simulations

Ionic liquids with cyano anions have long been used because of their unique combination of low-melting temperatures, reduced viscosities, and increased conductivities. Recently we showed that cyano anions in ionic liquids are particularly interesting for their potential use as electron donors to excited state photo-acceptors [B. Wu et al., J. Phys. Chem. B 119, 14790–14799 (2015)]. Here we report on bulk structural and quantum mechanical results for a series of ionic liquids based on the 1-ethyl-3-methylimidazolium cation, paired with the following five cyano anions: SeCN - , SCN - , N(CN) $-\atop{2}$ , C(CN) $-\atop{3}$ , and B(CN) $-\atop{4}$ . By combining molecular dynamics simulations, high-energy X-ray scattering measurements, and periodic boundary condition DFT calculations, we are able to obtain a comprehensive description of the liquid landscape as well as the nature of the HOMO-LUMO states for these ionic liquids in the condensed phase. Features in the structure functions for these ionic liquids are somewhat different than the commonly observed adjacency, charge-charge, and polarity peaks, especially for the bulkiest B(CN) $-\atop{4}$ anion. While the other four cyano-anion ionic liquids present an anionic HOMO, the one for Im $+\atop{2,1}$ /B(CN) $-\atop{4}$ is cationic.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Free-Flight CFD Simulations and Dynamic StabilityAnalysis of the Orion Crew Module

Dynamic stability analysis of the Multi-Purpose Crew Vehicle (MPCV) is performed us-ing the US3D flow solver and a plugin library to solve rigid-body dynamics, Free-Flight CFD(FF-CFD). Previous efforts investigated the free-flight behavior of higher speed, open-backaeroshell shapes. Three primary experimental sources are used to evaluate the predictivecapability of the FF-CFD solver in the low supersonic range (Mach≈1). First, the ballisticrange results obtained at the HFFAF facility at NASA Ames is used to verify the six degree-of-freedom (6-DoF) dynamic capability of the FF-CFD solver. Next, FF-CFD simulationsare preformed using restricted motion and the resultant trajectories are post-processed toobtain pitch damping coefficient as a function of angle-of-attack. The comparison of thepitch damping results obtained with one degree-of-freedom (1-DoF) FF-CFD compare wellwith experimental fits derived from ballistic range data. Finally, the atmospheric flight ca-pability of FF-CFD is compared to data from the Ascent Abort-2 (AA-2) flight experiment.Two simulations were performed using forced and forced-free flight to investigate surfacepressure predictive capability and free-flight aerodynamic performance through a varyingatmosphere at low speed (Mach 0.6-0.2). Surface pressure predicted with FF-CFD agreeswell with experimental trends, with slight over-prediction near the end of the trajectory.The total angle-of-attack for the free-flight portion agrees well with experimental data.

FFCFD↗

Real-time Simulation of Coordinated Sub-Transmission Volt-Var Control Tool under High Distributed PV Penetration

With the rapid penetration of intermittent solar photovoltaic (PV) and other distributed energy resources (DER) into the grid, and subsequent changes in power flow patterns in both distribution and sub-transmission, voltage regulation is becoming a major challenge. It is prudent to leverage PV and DER to provide ancillary services to the grid, such as voltage regulation. Recently, a quasi-static, coordinated real-time sub transmission volt-var control algorithm (CReST-VCT) was developed for voltage regulation under high PV penetration by dispatching the reactive power settings of the shunt devices and PV inverters. The algorithm was validated offline, on a quasistatic study. Such quasi-static dispatching algorithms cannot guarantee its performance for highly nonlinear dynamical power systems. In this paper the performance of the algorithm is validated with a real time nonlinear dynamic simulation of a modified IEEE 118 bus system in Opal-RT. The real time simulation is used to emulate the actual system operation, providing a more realistic testing environment for CReST-VCT. The dispatch and control signals are communicated between the power system (Opal solver) and the control center (GAMS solver) in real-time with a MODBUS bridge. The results demonstrate i) system is stable with the new dispatch points ii) significant improvement in system-wide voltage profiles compared to an uncontrolled scenario. Another significant contribution of the work is developing a framework for dispatchable, long duration dynamic simulations that can be leveraged for market/dispatch studies.

CReST-VCT, Opal-RT, Volt/Var Control, renewable in↗

Simulating bubble dynamics in a buoyant system

Summary Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two‐dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three‐dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.

Arai, Erin↗

Effect of Accurate Modeling of Converter Interfaced Generation (CIG) on a Practical Bulk Power System

With increased penetration of converter interfaced generation (CIG) in power systems, it becomes imperative to adequately model them in both steady-state and dynamic simulations. There are several modeling approaches and simulation platforms for CIG. In addition, it is noted that suitability of modeling approaches and models is also dependent on the system conditions. However, their suitability has yet to be systematically analyzed, and this becomes a challenge in the industry without a systemic method. To tackle these challenges, in this paper, we first develop a methodology for identifying the maximum penetration of aggregated rooftop sources that can be represented as negative loads in both the steady state and dynamic analysis. Then we verify the suitability of positive-sequence performance-based model and detailed electromagnetic transient (EMT) model for utility-scale PV plants under different fault conditions. Based on the results, we perform dynamic simulations on a large practical system with both positive-sequence transient stability(TS) and EMT-TS hybrid simulation. Our study shows that positive-sequence performance based PV plant model is more suitable for study cases with faults relatively far away from the point of connection, and that positive-sequence models and TS simulation approach can produce comparable results with reference to EMT-TS hybrid simulation for normally-cleared single-line-to-ground faults, but the inadequacy of positive-sequence based model and simulation approach becomes obvious for prolonged-tripping of unbalanced faults in proximity of the PV plants.

Huang, Qiuhua↗

Discovery of a metastable Al{sub 20}Sm{sub 4} phase

We present an efficient genetic algorithm, integrated with experimental diffraction data, to solve a nanoscale metastable Al{sub 20}Sm{sub 4} phase that evolves during crystallization of an amorphous magnetron sputtered Al{sub 90}Sm{sub 10} alloy. The excellent match between calculated and experimental X-ray diffraction patterns confirms an accurate description of this metastable phase. Molecular dynamic simulations of crystal growth from the liquid phase predict the formation of disordered defects in the devitrified crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Effect of precipitates in oxide film on hydrogen pick-up of zirconium alloys

In this study, we examined the effect of alloying elements on hydrogen pick-up. The alloying elements of Fe, Cr, and Ni form precipitates in Zircaloy because of their low dissolvability in Zirconium. The precipitates in oxide film are presumed to play an important role in the hydrogen pick-up process because oxide film is a pathway of hydrogen entry. A corrosion test was performed in pure water at 561 K, and the relation between the additive elements and the hydrogen pick-up was estimated by using Zircaloy-2 and Zircaloy-4. The behavior of precipitates in oxide film was estimated using TEM observation and molecular dynamic simulation. The results show that the un-oxidized precipitates had a tendency to substitute for the oxidized precipitates with distance from the metal-oxide film boundary. Additionally, precipitates containing a high concentration of oxygen were observed. The oxidation of the precipitates was confirmed. It was also shown that the Fe/Cr fraction of precipitates had a tendency to decrease with distance from the metal-oxide film boundary. It was confirmed that the diffusion coefficients of H{sup +} in the iron- or chromium-containing ZrO{sub 2} are smaller than those of the pure ZrO{sub 2} by molecular dynamics simulation. These results lead to the conclusion that the diffusion of Cr and Fe in oxide film contributes to the control of hydrogen pick-up on zirconium alloys as does the oxidation of precipitates.

36 MATERIALS SCIENCE↗

Importance of corners in fracture of highly crosslinked polymeric adhesives

In this study, very large molecular dynamics simulations with open ends between two solid adherends have been performed treating tensile deformation of coarse-grained, highly crosslinked polymer networks modeling epoxy systems. The open boundary and the presence of corners dramatically alter the fracture behavior. In contrast to systems with periodic boundaries, the failure strain decreases with increasing system size until a critical size is reached. This decrease greatly reduces the difference in the crack initiation strains between simulation and experiment. In the open geometry, the sides of the polymer network contract inward forming wedge shaped corners. The stress and strain are concentrated in the corners where the shear component is present and large. The nonuniformity of the strain results in accumulation of bond breaking in the corners and crack initiation there. Lastly, t the corner strain is system size dependent, which results in a system size dependence of the failure strain.

36 MATERIALS SCIENCE↗