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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Mechanism of Water Intrusion into Flexible ZIF-8: Liquid Is Not Vapor

Zeolitic Imidazolate Frameworks (ZIF) find applicationin storageand dissipation of mechanical energy. Their distinctive propertieslinked to their (sub)nanometer size and hydrophobicity allow for waterintrusion only under high hydrostatic pressure. Here we focus on thepopular ZIF-8 material investigating the intrusion mechanism in itsnanoscale cages, which is the key to its rational exploitation intarget applications. In this work, we used a joint experimental/theoreticalapproach combining in operando synchrotron experiments during high-pressureintrusion experiments, molecular dynamics simulations, and stochasticmodels to reveal that water intrusion into ZIF-8 occurs by a cascadefilling of connected cages rather than a condensation process as previouslyassumed. The reported results allowed us to establish structure/functionrelations in this prototypical microporous material, representingan important step to devise design rules to synthesize porous media.

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

Local well-posedness and singularity formation in non-Newtonian compressible fluids

We investigate the initial value problem of a very general class of 3 + 1 non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier–Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel–Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Identification of a quasi-liquid phase at solid–liquid interface

Abstract An understanding of solid–liquid interfaces is of great importance for fundamental research as well as industrial applications. However, it has been very challenging to directly image solid–liquid interfaces with high resolution, thus their structure and properties are often unknown. Here, we report a quasi-liquid phase between metal (In, Sn) nanoparticle surfaces and an aqueous solution observed using liquid cell transmission electron microscopy. Our real-time high-resolution imaging reveals a thin layer of liquid-like materials at the interfaces with the frequent appearance of small In nanoclusters. Such a quasi-liquid phase serves as an intermediate for the mass transport from the metal nanoparticle to the liquid. Density functional theory-molecular dynamics simulations demonstrate that the positive charges of In ions greatly contribute to the stabilization of the quasi-liquid phase on the metal surface.

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

Successive Procedure for Solution Verification Based on User Needs

This paper discusses a revised solution verification procedure for computational fluid dynamics simulations to estimate the uncertainties in the quantities of interest based on discretization error models. This proposed procedure builds upon current procedures described in ASME V&V 20 but provides more guidance in determining the necessary number of mesh levels to build reliable discretization error models. Such guidance is particularly useful for practicing engineers without prior experience in solution verification. The key features of this proposed solution verification procedure are the ability to determine the need for additional mesh levels iteratively and the seamless treatment for underdetermined, exact, and overdetermined solutions of the power series approximation to the discretization error models. This study applies the proposed procedure to a set of synthetic examples to demonstrate the revised procedure’s clarity in determining the number of mesh solutions required for a reliable estimate of the discretization error in computational fluid dynamics settings. Additionally, this proposed procedure prevents a potential pathway in the current procedure in ASME V&V 20 that may lead to unreasonably small discretization errors.

Weinmeister, Justin↗
DOE OSTI2020

Revisiting the π → π* transition of the nitrite ion at the air/water interface: A combined experimental and theoretical study

Broadband deep ultraviolet electronic sum frequency generation spectroscopy allows measurement of electronic |χ (2) | 2 -spectra at aqueous interfaces, providing considerable improvement over deep ultraviolet electronic second harmonic generation spectroscopy that yields error-prone, pointwise spectra. Reexamination of the π → π* transition of nitrite at the air/water interface reveals that the interfacial |χ (2) | 2 -spectrum is strikingly similar to the bulk absorption spectrum. Molecular dynamics simulations and SOS-CIS(D 0 ) electronic structure calculations provide no evidence for the previously reported contact ion pair-induced red-shift of the π → π* transition at the air/water interface. Finally, our results thus revise the previous description for nitrite interfacial adsorption as a contact ion pair.

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

Manipulating meso-scale solvent structure from Pd nanoparticle deposits in deep eutectic solvents

Deep Eutectic Solvents (DESs) are complex solutions that present unique challenges compared to traditional electrolytes. Unlike most electrolytes and ionic liquids, DESs have delicate hydrogen bond networks that are responsible for their highly sensitive compositional dependence on the melting point. Prior work has demonstrated unique nanoscale structure both experimentally and theoretically that brings both challenges and opportunities to their adoption in traditional electrochemical processes. In this study, we use in-situ sample-rotated ultra-small angle X-ray scattering to resolve the near-interface solvent structure after electrodepositing Pd nanoparticles onto a glassy carbon electrode in choline chloride:urea and choline chloride:ethylene glycol DESs. Our results indicate that solvent structure can be observed on the meso-scale in the choline chloride:urea and choline chloride:ethylene glycol systems. Importantly, this extended solvent structure peaks at -0.5 V (vs. Ag/AgCl) and decreases slightly with higher overpotentials during Pd electrodeposition. Experimentally, the nature of this structure is more pronounced in the ethylene glycol system, as evidenced by both the X-ray scattering as well as the electrochemical impedance spectroscopy (EIS). Here, molecular dynamics simulations and dipolar orientation analysis reveal that chloride delocalization near the Pd interface and long-range interactions between the choline and each hydrogen bond donor are very different and qualitatively consistent with the experimental data. These results show how the long-range solvent-deposit interactions can be tuned by changing the hydrogen bond donor in the DES and the applied potential.

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

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Trust-Region Approximation of Extreme Trajectories in Power System Dynamics

In this work we present a novel technique, based on a trust-region optimization algorithm and second-order trajectory sensitivities, to compute the extreme trajectories of power system dynamic simulations given a bounded set that represents parametric uncertainty. Furthermore, we show how this method, while remaining computationally efficient compared with sampling-based techniques, overcomes the limitations of previous sensitivity-based techniques to approximate the bounds of the trajectories when the local approximation loses validity because of the nonlinearity. We present several numerical experiments that showcase the accuracy and scalability of the technique, including a demonstration on the IEEE New England test system.

42 ENGINEERING↗
DOE OSTI2020

Quantifying the Influence of Jupiter on the Earth’s Orbital Cycles

A wealth of Earth-sized exoplanets will be discovered in the coming years, providing a large pool of candidates from which the targets for the search for life beyond the solar system will be chosen. The target selection process will require the leveraging of all available information in order to maximize the robustness of the target list and make the most productive use of follow-up resources. Here, we present the results of a suite of n-body simulations that demonstrate the degree to which the orbital architecture of the solar system impacts the variability of Earth’s orbital elements. By varying the orbit of Jupiter and keeping the initial orbits of the other planets constant, we demonstrate how subtle changes in solar system architecture could alter the Earth’s orbital evolution—a key factor in the Milankovitch cycles that alter the amount and distribution of solar insolation, thereby driving periodic climate change on our planet. The amplitudes and frequencies of Earth’s modern orbital cycles fall in the middle of the range seen in our runs for all parameters considered—neither unusually fast nor slow, neither large nor small. This finding runs counter to the “Rare Earth” hypothesis, which suggests that conditions on Earth are so unusual that life elsewhere is essentially impossible. Our results highlight how dynamical simulations of newly discovered exoplanetary systems could be used as an additional means to assess the potential targets of biosignature searches, and thereby help focus the search for life to the most promising targets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The effect of protein mutations on drug binding suggests ensuing personalised drug selection

The advent of personalised medicine promises a deeper understanding of mechanisms and therefore therapies. However, the connection between genomic sequences and clinical treatments is often unclear. We studied 50 breast cancer patients belonging to a population-cohort in the state of Qatar. From Sanger sequencing, we identified several new deleterious mutations in the estrogen receptor 1 gene (ESR1). The effect of these mutations on drug treatment in the protein target encoded by ESR1, namely the estrogen receptor, was achieved via rapid and accurate protein–ligand binding affinity interaction studies which were performed for the selected drugs and the natural ligand estrogen. Four nonsynonymous mutations in the ligand-binding domain were subjected to molecular dynamics simulation using absolute and relative binding free energy methods, leading to the ranking of the efficacy of six selected drugs for patients with the mutations. Our study shows that a personalised clinical decision system can be created by integrating an individual patient’s genomic data at the molecular level within a computational pipeline which ranks the efficacy of binding of particular drugs to variant proteins.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Computational approaches to model X-ray photon correlation spectroscopy from molecular dynamics

X-ray photon correlation spectroscopy (XPCS) allows for the resolution of dynamic processes within a material across a wide range of length and time scales. X-ray speckle visibility spectroscopy (XSVS) is a related method that uses a single diffraction pattern to probe ultrafast dynamics. Interpretation of the XPCS and XSVS data in terms of underlying physical processes is necessary to establish the connection between the macroscopic responses and the microstructural dynamics. To aid the interpretation of the XPCS and XSVS data, we present a computational framework to model these experiments by computing the X-ray scattering intensity directly from the atomic positions obtained from molecular dynamics simulations. We compare the efficiency and accuracy of two alternative computational methods: the direct method computing the intensity at each diffraction vector separately, and a method based on fast Fourier transform that computes the intensities at all diffraction vectors at once. Furthermore, the computed X-ray speckle patterns capture the density fluctuations over a range of length and time scales and are shown to reproduce the known properties and relations of experimental XPCS and XSVS for liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Broadening of the drumhead-mode spectrum due to in-plane thermal fluctuations of two-dimensional trapped ion crystals in a Penning trap

Two-dimensional crystals of ions stored in Penning traps are a leading platform for quantum simulation and sensing experiments. For small amplitudes, the out-of-plane motion of such crystals can be described by a discrete set of normal modes called the drumhead modes, which can be used to implement a range of quantum information protocols. However, experimental observations of crystals with Doppler-cooled and even nearground- state-cooled drumhead modes reveal an unresolved drumhead-mode spectrum. In this work, we establish in-plane thermal fluctuations in ion positions as a major contributor to the broadening of the drumhead-mode spectrum. In the process, we demonstrate how the confining magnetic field leads to unconventional in-plane normal modes, whose average potential and kinetic energies are not equal. This property, in turn, has implications for the sampling procedure required to choose the in-plane initial conditions for molecular-dynamics simulations. Further, for current operating conditions of the NIST Penning trap, our study suggests that the two-dimensional crystals produced in this trap undergo in-plane potential-energy fluctuations of the order of 10mK. Our study therefore motivates the need for designing improved techniques to cool the in-plane degrees of freedom.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Polarizable Water Potential Derived from a Model Electron Density

A new empirical potential for efficient, large scale molecular dynamics simulation of water is presented. The HIPPO (Hydrogen-like Intermolecular Polarizable POtential) force field is based upon the model electron density of a hydrogen-like atom. This framework is used to derive and parametrize individual terms describing charge penetration damped permanent electrostatics, damped polarization, charge transfer, anisotropic Pauli repulsion, and damped dispersion interactions. Initial parameter values were fit to Symmetry Adapted Perturbation Theory (SAPT) energy components for ten water dimer configurations, as well as the radial and angular dependence of the canonical dimer. The SAPT-based parameters were then systematically refined to extend the treatment to water bulk phases. The final HIPPO water model provides a balanced representation of a wide variety of properties of gas phase clusters, liquid water, and ice polymorphs, across a range of temperatures and pressures. This water potential yields a rationalization of water structure, dynamics, and thermodynamics explicitly correlated with an ab initio energy decomposition, while providing a level of accuracy comparable or superior to previous polarizable atomic multipole force fields. The HIPPO water model serves as a cornerstone around which similarly detailed physics-based models can be developed for additional molecular species.

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

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗
DOE OSTI2024

Neutrons in Structural Biology: Challenges and Opportunities (Workshop Report)

Gaining a thorough understanding of biological systems requires building our knowledge about biological processes from the level of atoms and electrons, and up to whole organisms. Such comprehensive knowledge will allow for a predictive understanding of complex biological systems behavior. It will guide us in the design and development of novel therapeutics and vaccines to tackle existing health threats and to prepare for future pandemics, and it will provide information necessary to create new biomaterials and bio-inspired technologies through manipulation of biological macromolecules, their assemblies, single cells and even microorganisms. Reaching these goals will require a synergistic combination of multiple experimental techniques with molecular calculations and predictive simulations, and the design and development of new techniques and capabilities that bridge current knowledge and technology gaps. Neutron scattering provides unique information about the biomacromolecular structure and function and can play a major role in achieving these goals. A workshop was held to engage the scientific community in identifying pressing challenges in biochemistry, structural biology, enzymology and structure-guided drug design not solved with the current neutron scattering technologies or utilizing other structural biology techniques such as X-ray crystallography, NMR, and cryo-EM. The workshop brought together structural biology, biochemistry and computational experts, as well as early career researchers and students, creating a forum for discussing scientific advancement and collaboration. The workshop included a one-day satellite training workshop where graduate students and postdoctoral researchers were educated in the application of neutron crystallography and small-angle scattering in structural biology. Furthermore, the Instrument Scientific Advisory Board (ISAB) for the development of a macromolecular neutron diffractometer at ORNL’s Second Target Station was introduced at the workshop. The major outcome was that neutrons can provide atomic-level understanding of biomacromolecular structure, function and dynamics which is of paramount importance for addressing the identified challenges. Neutron crystallography, in particular, can resolve long-standing biochemical issues regarding enzyme function by delineating the underlying chemistry and can have a major impact on the design of small-molecule therapeutics, especially in combination with molecular computation (quantum chemistry and molecular dynamics simulations) and the emerging artificial intelligence (AI)-assisted drug design technologies. The unique properties of neutrons, including their high sensitivity to hydrogen and their non-destructive nature, make them ideal probes of biological matter. There is a palpable need in the scientific community to expand and enhance the impact of neutron sciences on biology. Neutron crystallography is the only structural biology method capable of determining positions of all hydrogen atoms in proteins, nucleic acids and their complexes at near-physiological temperatures and of unstable species at cryogenic temperatures. Moreover, neutron analysis is non-ionizing, non-destructive and does not perturb the structure or redox chemistry of active site metal centers and clusters in proteins, which can be invaluable for studying radiation-sensitive metalloprotein complexes. Further, neutron energies used in scattering applications are similar to atomic motions, permitting neutron spectroscopies to characterize the dynamics of biomacromolecules on the picosecond to microsecond timescales. The different sensitivities of neutrons to protium (H) and deuterium (D) isotopes of hydrogen allow enhanced visibility of specific parts of biological complexes through isotopic labeling. The impact of neutrons will be most powerful when neutron scattering is combined with complementary experimental techniques that use photons and electrons, and with high-performance computing. The interconnection and mutuality of the experimental and theoretical capabilities will drive discoveries in biological and health sciences to generate more complete picture of complex biological systems. The major limitation in the field of biological neutron crystallography has been signal-to-noise, demanding large samples that are difficult to produce for the majority of biomacromolecules and limiting the applicability of this technique in biological sciences. A neutron crystallography instrument at the Second Target Station will revolutionize biological science with neutrons by engaging a large scientific community of structural biologists, enabling successful neutron diffraction experiments from radically smaller biomacromolecular crystals, resolving unanswered biochemical questions, and meaningfully contributing to rational drug design. The meeting highlighted 10 grand challenges that will be addressed with this advanced capability over the next decade and beyond, and the recommendations required to help address them are given below.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

A molecular dynamics study of collisional heat transfer to nanoclusters in the gas phase

Nanoclusters in the gas phase grow by condensation and coagulation, which heat up nanoclusters, creating thermal non-equilibrium with the surrounding gas. The extent of non-equilibrium for nanoclusters is dependent upon the rate of heating relative to the rate of cooling brought about by nanocluster-gas molecule collisions, hence the latter process plays an important role in nanocluster energy evolution during growth. We apply molecular dynamics simulations to investigate heat transfer between metal nanoclusters and gas molecules relevant to high temperature aerosol systems. In analysis, we first define a thermal reemission coefficient (TRC), which is a correction factor for the portion of energy possessed by reemitted gas molecules after collision. Here, the TRC differs from the thermal accommodation coefficient (TAC), which is a correction factor for the maximum possible energy transfer rate from nanoclusters to gas molecules based on the diffusive collision model. We find the TRC is insensitive to nanocluster size but is affected by the nanocluster to gas temperature ratio and the nanocluster to gas atomic mass ratio. Comparison to the traditionally-defined TAC shows that the TAC is insensitive to temperature, a phenomenon that has been widely indicated in literature but not quantitatively explained. We show how the temperature insensitivity of the TAC arises because of the manner in which TRCs vary with both nanocluster and gas temperature, i.e. the manner in which the reemitted gas molecule energy from a surface varies with the temperature of that surface relative to that of the gas. In addition, while calculations ultimately reveal temperature insensitive TACs and justify their continued use in modeling gas-phase heat transfer, we argue that the TRC is more appropriate to calculate a priori in analyzing simulations; it is physically more realistic to first introduce a correction factor solely for the energy of reemitted gas molecules after a collision (as TRCs) instead of a correction factor applying to both the reemission energy and the initial energy (as TACs), as the latter is exactly calculable and independent of nanocluster properties. In line with prior work, for monoatomic gases we find the TAC decreases as the nanocluster to gas atomic mass increases, approaching 0.08 for gold nanoclusters in helium.

77 NANOSCIENCE AND NANOTECHNOLOGY↗