Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic 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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Origin of increased helium density inside bubbles in Ni (1–$x$) Fe $x$ alloys

Due to virtually no solubility, He atoms implanted or created inside materials tend to form bubbles, which are known to damage material properties through embrittlement. Higher He density in nano-sized bub- bles was observed both experimentally and computationally in Ni (1–$x$) Fe $x$ -alloy samples compared to Ni. The bubbles in the Ni (1–$x$) Fe $x$ -alloys were observed to be faceted, whereas in elemental Ni they were more spherical. Molecular dynamics simulations showed that stacking fault structures formed around bubbles at maximum He density. Higher Fe concentrations stabilize stacking fault structures, suppress evolution of dislocation network around bubbles and suppress complete dislocation emission, leading to higher He density.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Enabling robust offline active learning for machine learning potentials using simple physics-based priors

Abstract Machine learning surrogate models for quantum mechanical simulations have enabled the field to efficiently and accurately study material and molecular systems. Developed models typically rely on a substantial amount of data to make reliable predictions of the potential energy landscape or careful active learning (AL) and uncertainty estimates. When starting with small datasets, convergence of AL approaches is a major outstanding challenge which has limited most demonstrations to online AL. In this work we demonstrate a Δ-machine learning (ML) approach that enables stable convergence in offline AL strategies by avoiding unphysical configurations with initial datasets as little as a single data point. We demonstrate our framework’s capabilities on a structural relaxation, transition state calculation, and molecular dynamics simulation, with the number of first principle calculations being cut down anywhere from 70%–90%. The approach is incorporated and developed alongside AMP torch , an open-source ML potential package, along with interactive Google Colab notebook examples.

Shuaibi, Muhammed↗
DOE OSTI2021

Energetics of polymeric carbon monoxide

The transformation of carbon monoxide (CO) from a molecular liquid to a polymeric solid under isothermal compression at room temperature is investigated using first principles theory. We report structural and thermodynamic properties from ambient density up to 2.45 g/cc obtained using density functional theory molecular dynamics simulations, including hybrid exchange corrections. The theoretical results are compared with newly obtained polymeric CO samples, synthesized in a large volume press. The explosive performance of polymeric CO is predicted and discussed. Under most favorable assumptions, it is found to be comparable to trinitrotoluene.

25 ENERGY STORAGE↗
DOE OSTI2022

Role of Local Structure in the Enhanced Dynamics of Deformed Glasses

External stress can accelerate molecular mobility of amorphous solids by several orders of magnitude. The changes in mobility are commonly interpreted through the Eyring model, which invokes an empirical activation volume. Here, we analyze constant-stress molecular dynamics simulations and propose a structure-dependent Eyring model, connecting activation volume to a machine-learned field, softness. We show that stress has a heterogeneous effect on the mobility that depends on local structure through softness. The barrier impeding relaxation reduces more for well-packed particles, which explains the narrower distribution of relaxation time observed under stress.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Functional lipid pairs as building blocks of phase-separated membranes

Biological membranes exhibit a great deal of compositional and phase heterogeneity due to hundreds of chemically distinct components. As a result, phase separation processes in cell membranes are extremely difficult to study, especially at the molecular level. It is currently believed that the lateral membrane heterogeneity and the formation of domains, or rafts, are driven by lipid–lipid and lipid–protein interactions. Nevertheless, the underlying mechanisms regulating membrane heterogeneity remain poorly understood. In the present work, we combine inelastic X-ray scattering with molecular dynamics simulations to provide direct evidence for the existence of strongly coupled transient lipid pairs. These lipid pairs manifest themselves experimentally through optical vibrational (a.k.a. phononic) modes observed in binary (1,2-dipalmitoyl- sn -glycero-3-phosphocholine [DPPC]–cholesterol) and ternary (DPPC–1,2-dioleoyl- sn -glycero-3-phosphocholine/1-palmitoyl-2-oleoyl-glycero-3-phosphocholine [DOPC/POPC]–cholesterol) systems. The existence of a phononic gap in these vibrational modes is a direct result of the finite size of patches formed by these lipid pairs. The observation of lipid pairs provides a spatial (subnanometer) and temporal (subnanosecond) window into the lipid–lipid interactions in complex mixtures of saturated/unsaturated lipids and cholesterol. Our findings represent a step toward understanding the lateral organization and dynamics of membrane domains using a well-validated probe with a high spatial and temporal resolution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Superlubric-pinned transition of a two-dimensional solid dusty plasma under a periodic triangular substrate

Here, the superlubric-pinned transition in the depinning dynamics of a two-dimensional (2D) solid dusty plasma modulated by 2D triangular periodic substrates is investigated using Langevin dynamical simulations. When the lattice structure of the 2D solid dusty plasma perfectly matches the triangular substrate, two distinctive pinned and moving ordered states are observed as the external uniform driving force gradually increases from zero. When there is a mismatch between the lattice structure and the triangular substrate, however, on shallow substrates, it is discovered that all of the particles can slide freely on the substrate even when the applied driving force is tiny. This is a typical example of superlubricity, which is caused by the competition between the substrate-particle and particle-particle interactions. If the substrate depth increases further, as the driving force increases from zero, there are three dynamical states consisting of the pinned state, the disordered plastic flow state, and the moving ordered state. In an underdense system, where there are fewer particles than potential well minima, it is found that the occurrence of the three different dynamical states is controlled by the depth of the substrate, which is quantitatively characterized using the average mobility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a machine learning model for polyethylene pyrolysis using a detailed reaction mechanism

Waste plastics have recently received significant attention as the issue of waste generation continues to increase. Thermal conversion processes, such as pyrolysis and gasification, are attractive potential technologies for utilizing waste plastics and reducing overall waste generation. Efficient utilization of plastics requires a detailed understanding of the conversion process such as pyrolysis and gasification. However, a mechanistic understanding of these processes lead to large and complex kinetic schemes that are not suited for large-scale and long-time simulation methods. Currently, most modeling approaches for pyrolysis and gasification rely on globally lumped, simplified kinetic schemes that provide results that are classified by their product type and not individual species, which limit the level of fidelity achieved via modeling. A machine learning (ML) model has been developed for the primary reactions of high-density polyethylene (HDPE) in an attempt to increase computational efficiency while still maintaining a high level of detail and accuracy. The ML model is trained on a detailed reaction mechanism containing 42 total species and 737 chemical reactions. A DeepONet branch and trunk architecture was adopted to train the model using time-steps relevant to computational fluid dynamics simulations. The ML used physics-informed loss functions to ensure mass conservation. The surrogate model has been deployed in simple MFiX CFD simulations, single particle and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗
DOE OSTI2024

Thermodynamic driving forces in contact electrification between polymeric materials

Abstract Contact electrification, or contact charging, refers to the process of static charge accumulation after rubbing, or even simple touching, of two materials. Despite its relevance in static electricity, various natural phenomena, and numerous technologies, contact charging remains poorly understood. For insulating materials, even the species of charge carrier may be unknown, and the direction of charge-transfer lacks firm molecular-level explanation. Here, we use all-atom molecular dynamics simulations to investigate whether thermodynamics can explain contact charging between insulating polymers. Based on prior work suggesting that water-ions, such as hydronium and hydroxide ions, are potential charge carriers, we predict preferred directions of charge-transfer between polymer surfaces according to the free energy of water-ions within water droplets on such surfaces. Broad agreement between our predictions and experimental triboelectric series indicate that thermodynamically driven ion-transfer likely influences contact charging of polymers. Furthermore, simulation analyses reveal how specific interactions of water and water-ions proximate to the polymer-water interface explain observed trends. This study establishes relevance of thermodynamic driving forces in contact charging of insulators with new evidence informed by molecular-level interactions. These insights have direct implications for future mechanistic studies and applications of contact charging involving polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of strong Coulomb coupling on diffusion in atmospheric pressure plasmas

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion–ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions (x i ≲ 10 —6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions (10 —6 ≲ x i ≲ 10 —2 ) there is a transition from ion-neutral to ion–ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion–ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion–ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions (x i ≳ 10 —2 ), strongly correlated ion–ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. Importantly, the results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

fast neutral gas heating↗
DOE OSTI2024

Impact of the Exciter and Governor Parameters on Forced Oscillations

In recent years, the frequency of forced oscillation events due to control system malfunctions or improper parameter settings has increased. Tuning the parameters of exciters and governor models is crucial for maintaining power system stability. Traditional simulation studies typically involve small transient disturbances or step changes to find optimal parameter sets, but existing optimization algorithms often fall short in fine-tuning for forced oscillations. Identifying the sensitive parameters within these control models is essential for ensuring stability during large, sustained disturbances. This study focuses on identifying these critical exciter and governor model parameters by analyzing their influence on sustained forced oscillations. Using Kundur’s two-area system, we analyze common exciter models such as SCRX, ESST1A, and AC7B, along with governor models like GAST, HYGOV, and GGOV1, utilizing PSS®E software version 34. Sustained forced oscillations are injected at generator-1 of area-1, with individual parameter changes dynamically simulated. By considering a local oscillation frequency of 1.4 Hz and an inter-area oscillation mode of 0.25 Hz, we analyze the impact of each parameter change on the magnitude and frequency of forced oscillations as well as on active and reactive power outputs. This novel approach highlights the most influential parameters of each tested model—such as exciter, governor, and turbine gains, as well as time constant parameters—on the impact of forced oscillations. Based on our findings, the sensitive parameters of each tested model are ranked. These would provide valuable insights for industry operators to fine-tune control settings during oscillation events, ultimately enhancing system stability.

42 ENGINEERING↗
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Defect formation and bending properties in graphite under He atom implantation investigated by molecular dynamics method

The formation and evolution of radiation defects and their effects on bending property of a single graphite under low energy helium atom implantation have been studied at atomic scale through molecular dynamics simulations. Various simple interstitial and vacancy defects are observed during the implantation process, which is a function of incident energy and cross section of helium-carbon interaction. Large defects are observed after continuous helium atom implantations. The effects of incident energy on the density, size and distribution of defects are analyzed in detail. Furthermore, the kinetic evolution process of radiation defects at different annealing temperatures have also been simulated, including the migration, aggregation, transformation, coalescence, and self-healing. Here, the large defect clusters, crossing several carbon layers, have been confirmed to be the primary reasons that affect the bending property of graphite. All these results provide a new understanding for further applications of graphite in nuclear reactors.

Atomic implantation↗
DOE OSTI2020

Entropy production in model colloidal suspensions under shear via the fluctuation theorem

Dissipative systems often exhibit novel and unexpected properties. This is, for instance, the case of simple liquids, which, when subjected to shear and after reaching a steady state, can exhibit a negative entropy production over finite length scales and timescales. This result, among others, is captured by nonequilibrium relations known as fluctuation theorems. Using nonequilibrium molecular dynamics simulations, we examine how, by fine-tuning the properties of the components of a complex fluid, we can steer the nonequilibrium response of the fluid. More specifically, we show how we control the nonequilibrium probability distribution for the shear stress and, in turn, how often states with a negative entropy production can occur. To achieve this, we start by characterizing how the size for the liquid matrix impacts the probability of observing negative entropy states, as well as the timescale over which these can be observed. We then measure how the addition of larger particles to this liquid matrix, i.e., simulating a model colloidal suspension, results in an increase in the occurrence of such states. As a result, this suggests how modifications in the composition of the mixture and in the properties of its components lead to an increase in the probability of observing states of negative entropy production and, thus, for the system to run in reverse.

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

Depletion-Driven Assembly of Polymer-Coated Nanocrystals

Depletion-driven assembly has been widely studied for micron-sized colloids, but questions remain at the nanoscale where the governing physics are impacted by stabilizing surface ligands or wrapping polymers, whose dimensions are similar to those of the colloidal core and the depletant. Here, we probe how wrapping colloidal tin-doped indium oxide nanocrystals with polymers affects their depletion-induced interactions and assembly in solutions of polyethylene glycol. Copolymers of polyacrylic acid grafted with polyethylene oxide provide nanocrystal wrappings with different effective polymer graft densities and molecular weights. (Ultra) small angle X-ray scattering, coarse-grained molecular dynamics simulation, and molecular thermodynamic theory were combined to analyze how depletant size and polymer wrapping characteristics affect depletion interactions, structure, and phase behavior. The results show how depletant molecular weight, as well as surface density and molecular weight of polymer grafts, set thresholds for assembly. These signatures are unique to depletion-driven assembly of nanoscale colloids and mirror phase behaviors of grafted nanoparticle–polymer composites. Optical and rheological responses of depletion-driven assemblies of nanocrystals with different polymer shell architectures were probed to learn how their structural differences impact properties. Finally, we discuss how these wrinkles for depletion-driven assembly at the nanoscale may provide fresh opportunities for designing responsive depletion interactions and dynamically reconfigurable materials.

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

Observations of grain-boundary phase transformations in an elemental metal

The theory of grain boundary (the interface between crystallites, GB) structure has a long history 1 and the concept of GBs undergoing phase transformations was proposed 50 years ago 2,3 . The underlying assumption was that multiple stable and metastable states exist for different GB orientations 4,5,6 . The terminology ‘complexion’ was recently proposed to distinguish between interfacial states that differ in any equilibrium thermodynamic property 7 . Different types of complexion and transitions between complexions have been characterized, mostly in binary or multicomponent systems 8,9,10,11,12,13,14,15,16,17,18,19 . Simulations have provided insight into the phase behaviour of interfaces and shown that GB transitions can occur in many material systems 20,21,22,23,24 . However, the direct experimental observation and transformation kinetics of GBs in an elemental metal have remained elusive. In this work we demonstrate atomic-scale GB phase coexistence and transformations at symmetric and asymmetric $[11\bar{1}]$ tilt GBs in elemental copper. Atomic-resolution imaging reveals the coexistence of two different structures at Σ19b GBs (where Σ19 is the density of coincident sites and b is a GB variant), in agreement with evolutionary GB structure search and clustering analysis 21,25,26 . We also use finite-temperature molecular dynamics simulations to explore the coexistence and transformation kinetics of these GB phases. Our results demonstrate how GB phases can be kinetically trapped, enabling atomic-scale room-temperature observations. Our work paves the way for atomic-scale in situ studies of metallic GB phase transformations, which were previously detected only indirectly 9,15,27,28,29 , through their influence on abnormal grain growth, non-Arrhenius-type diffusion or liquid metal embrittlement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Water under hydrophobic confinement: entropy and diffusion

The properties of liquid water are known to change drastically in confined geometries. A most interesting and intriguing phenomenon is that the diffusion of water is found to be strongly enhanced by the proximity of a hydrophobic confining wall relative to the bulk diffusion. We report a molecular dynamics simulation using a classical water model investigating the water diffusion near a non-interacting smooth confining wall, which is assumed to imitate a hydrophobic surface, revealing a pronounced diffusion enhancement within several water layers adjacent to the wall. We present evidence that the observed diffusion enhancement can be accounted for, with a quantitative accuracy, using the universal scaling law for liquid diffusion that relates the diffusion rate to the excess entropy. These results show that the scaling law, which has so far only been used for the description of the diffusion in simple liquids, can successfully describe the diffusion in water. It is shown that the law can be used for the analysis of water dynamics under nanoscale hydrophobic confinement, which is currently a subject of intense research activity.

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

Tuning conformational structures of imidazolium ionenes with 1-ethyl-3-methylimidazolium ionic liquid solvents

In this report the structures and conformations of imidazolium ionenes solvated in 1-ethyl-3-methylimidazolium ionic liquids (ILs) ([C 2 mim + ][X - ]) are investigated using molecular dynamics simulations. Four different ionenes poly(decylimidazolium) (PD 10 ), poly(tetraethyleneglycolimidazolium) (PE 10 ), alternating copolymer (P(ED) 5 ), and a block copolymer PE 5 D 5 are immersed in [C 2 mim + ][X - ], where [X - ] = thiocyanate [SCN - ], tetrafluoroborate [BF 4 - ], and trifluoromethanesulfonate [TfO - ]. The radius of gyration indicates that the ionene chains are more extended when immersed in [TfO - ] compared to [BF 4 - ] and [SCN - ], while the solvent accessible surface area shows that PE 10 is the most contracted ionene among all systems. The electrostatic interactions between the ionenes and the ILs are very consistent with the inherent electrostatic characteristics of the different anions, which can be quantified using the Ionic Polarity Index (IPI). The ionene configurations and dynamics strongly depend on the backbone functionality and architecture of the chain, due to the balance between the inter- and intra-molecular interactions in these systems.

42 ENGINEERING↗
DOE OSTI2022

Active learning of reactive Bayesian force fields applied to heterogeneous catalysis dynamics of H/Pt

Abstract Atomistic modeling of chemically reactive systems has so far relied on either expensive ab initio methods or bond-order force fields requiring arduous parametrization. Here, we describe a Bayesian active learning framework for autonomous “on-the-fly” training of fast and accurate reactive many-body force fields during molecular dynamics simulations. At each time-step, predictive uncertainties of a sparse Gaussian process are evaluated to automatically determine whether additional ab initio training data are needed. We introduce a general method for mapping trained kernel models onto equivalent polynomial models whose prediction cost is much lower and independent of the training set size. As a demonstration, we perform direct two-phase simulations of heterogeneous H 2 turnover on the Pt(111) catalyst surface at chemical accuracy. The model trains itself in three days and performs at twice the speed of a ReaxFF model, while maintaining much higher fidelity to DFT and excellent agreement with experiment.

42 ENGINEERING↗