Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion problems”

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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A space-time smooth artificial viscosity method with wavelet noise indicator and shock collision scheme, Part 2: The 2-D case

This is the second part to our companion paper. Herein, we generalize to two space dimensions the C-method developed in for adding localized, space-time smooth artificial viscosity to nonlinear systems of conservation laws that propagate shock waves, rarefaction waves, and contact discontinuities. For gas dynamics, the C-method couples the Euler equations to scalar reaction-diffusion equations, which we call C-equations, whose solutions serve as space-time smooth artificial viscosity indicators for shocks and contacts. We develop a high-order numerical algorithm for gas dynamics in 2-D which can accurately simulate the Rayleigh-Taylor (RT) instability with Kelvin-Helmholtz (KH) roll-up of the contact discontinuity, as well as shock collision and bounce-back. Solutions to our C-equations not only indicate the location of the shocks and contacts, but also track the geometry of the evolving fronts. This allows us to implement both directionally isotropic and anisotropic artificial viscosity schemes, the latter adding diffusion only in directions tangential to the evolving front. We additionally produce a novel shock collision indicator function, which naturally activates during shock collision, and then smoothly deactivates. Moreover, we implement a high-frequency 2-D wavelet-based noise detector together with an efficient and localized noise removal algorithm. To test the methodology, we use a highly simplified WENO-based discretization scheme. We provide numerical results for some classical 2-D test problems, including the RT problem, the Noh problem, a circular explosion problem from the Liska & Wendroff [13] review paper, the Sedov blast wave problem, the double Mach 10 reflection test, and a shock-wall collision problem. In particular, we show that our artificial viscosity method can eliminate the wall-heating phenomenon for the Noh problem, and thereby produce an accurate, non-oscillatory solution, even though our simplified WENO-type scheme fails to run for this problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Deep learning analysis of defect and phase evolution during electron beam-induced transformations in WS2

Abstract Recent advances in scanning transmission electron microscopy (STEM) allow the real-time visualization of solid-state transformations in materials, including those induced by an electron beam and temperature, with atomic resolution. However, despite the ever-expanding capabilities for high-resolution data acquisition, the inferred information about kinetics and thermodynamics of the process, and single defect dynamics and interactions is minimal. This is due to the inherent limitations of manual ex situ analysis of the collected volumes of data. To circumvent this problem, we developed a deep-learning framework for dynamic STEM imaging that is trained to find the lattice defects and apply it for mapping solid state reactions and transformations in layered WS 2 . The trained deep-learning model allows extracting thousands of lattice defects from raw STEM data in a matter of seconds, which are then classified into different categories using unsupervised clustering methods. We further expanded our framework to extract parameters of diffusion for sulfur vacancies and analyzed transition probabilities associated with switching between different configurations of defect complexes consisting of Mo dopant and sulfur vacancy, providing insight into point-defect dynamics and reactions. This approach is universal and its application to beam-induced reactions allows mapping chemical transformation pathways in solids at the atomic level.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Elucidation of the physical factors that control activated transport of penetrants in chemically complex glass-forming liquids

Understanding the activated transport of penetrant or tracer atoms and molecules in condensed phases is a challenging problem in chemistry, materials science, physics, and biophysics. Many angstrom- and nanometer-scale features enter due to the highly variable shape, size, interaction, and conformational flexibility of the penetrant and matrix species, leading to a dramatic diversity of penetrant dynamics. Based on a minimalist model of a spherical penetrant in equilibrated dense matrices of hard spheres, a recent microscopic theory that relates hopping transport to local structure has predicted a novel correlation between penetrant diffusivity and the matrix thermodynamic dimensionless compressibility, S 0 ( T ) (which also quantifies the amplitude of long wavelength density fluctuations), as a consequence of a fundamental statistical mechanical relationship between structure and thermodynamics. Moreover, the penetrant activation barrier is predicted to have a factorized/multiplicative form, scaling as the product of an inverse power law of S 0 ( T ) and a linear/logarithmic function of the penetrant-to-matrix size ratio. This implies an enormous reduction in chemical complexity that is verified based solely on experimental data for diverse classes of chemically complex penetrants dissolved in molecular and polymeric liquids over a wide range of temperatures down to the kinetic glass transition. The predicted corollary that the penetrant diffusion constant decreases exponentially with inverse temperature raised to an exponent determined solely by how S 0 ( T ) decreases with cooling is also verified experimentally. Our findings are relevant to fundamental questions in glassy dynamics, self-averaging of angstrom-scale chemical features, and applications such as membrane separations, barrier coatings, drug delivery, and self-healing.

Science & Technology - Other Topics↗
DOE OSTI2015

Compressed-sensing-based content-driven hierarchical reconstruction: Theory and application to C-arm cone-beam tomography

Purpose: This paper addresses the reconstruction of x-ray cone-beam computed tomography (CBCT) for interventional C-arm systems. Subsampling of CBCT is a significant issue with C-arms due to their slow rotation and to the low frame rate of their flat panel x-ray detectors. The aim of this work is to propose a novel method able to handle the subsampling artifacts generally observed with analytical reconstruction, through a content-driven hierarchical reconstruction based on compressed sensing. Methods: The central idea is to proceed with a hierarchical method where the most salient features (high intensities or gradients) are reconstructed first to reduce the artifacts these features induce. These artifacts are addressed first because their presence contaminates less salient features. Several hierarchical schemes aiming at streak artifacts reduction are introduced for C-arm CBCT: the empirical orthogonal matching pursuit approach with the ℓ{sub 0} pseudonorm for reconstructing sparse vessels; a convex variant using homotopy with the ℓ{sub 1}-norm constraint of compressed sensing, for reconstructing sparse vessels over a nonsparse background; homotopy with total variation (TV); and a novel empirical extension to nonlinear diffusion (NLD). Such principles are implemented with penalized iterative filtered backprojection algorithms. For soft-tissue imaging, the authors compare the use of TV and NLD filters as sparsity constraints, both optimized with the alternating direction method of multipliers, using a threshold for TV and a nonlinear weighting for NLD. Results: The authors show on simulated data that their approach provides fast convergence to good approximations of the solution of the TV-constrained minimization problem introduced by the compressed sensing theory. Using C-arm CBCT clinical data, the authors show that both TV and NLD can deliver improved image quality by reducing streaks. Conclusions: A flexible compressed-sensing-based algorithmic approach is proposed that is able to accommodate for a wide range of constraints. It is successfully applied to C-arm CBCT images that may not be so well approximated by piecewise constant functions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

HARD: A performance portable radiation hydrodynamics code based on FleCSI framework

Hydrodynamics And Radiation Diffusion (HARD) is an open-source application for high-performance simulations of compressible hydrodynamics with radiation-diffusion coupling. Built on the FleCSI (Bergen et al., 2021 [1]) (Flexible Computational Science Infrastructure) framework, HARD expresses its computational units as tasks whose execution can be orchestrated by multiple back-end runtimes, including Legion (Bauer et al., 2012 [2]), MPI (Forum, 1994 [3]), and HPX (Kaiser et al., 2020 [4]). Node-level parallelism is handled through Kokkos (Edwards et al., 2014 [5]), providing a single-source, portable code base that runs efficiently on laptops, small homogeneous clusters, and the largest heterogeneous supercomputers currently available. To ensure scientific reliability, HARD includes a regression test suite that automatically reproduces canonical verification problems such as the Sod and LeBlanc shock tubes, and the Sedov blast wave, comparing numerical solutions against known analytical results. The project is distributed under an OSI-approved license, hosted on GitHub, and accompanied by reproducible build scripts and continuous integration workflows. This combination of performance portability, verification infrastructure, and community-focused development makes HARD a sustainable platform for advancing radiation hydrodynamics research across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Crystal-growth-dominated fabrication of metal-organic frameworks with orderly distributed hierarchical porosity

Constructing architectures with hierarchical porosity has been widely considered as the most efficient way to bypass the problems related to slow mass transfer and inaccessibility of internal space in MOFs. Now, a crystal-growth-dominated strategy is proposed to fabricate hierarchically porous MOFs (HP-MOFs). When the crystal growth is dominated by the monomer attachment, the aggregation of nonionic surfactant or polymer can be easily captured and released during the crystal growth process, resulting in the formation and ordering hierarchical pores along the radial direction. Furthermore, owing to the accelerated mass diffusion and more exposed active sites of this design, HP-MOFs exhibited an enhanced catalytic efficiency in styrene oxidation.

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

Using Radiogenic Noble Gas Nuclides to Identify and Characterize Rock Fracturing

Abstract Fracture‐released radiogenic noble gas nuclides are used to identify locations and constrain the volume of new fracture creation during subsurface detonations. Real‐time, in situ noble gases and reactive gases were monitored using a field‐deployed mass spectrometer and automated sampling system in a multilevel borehole array. Released gases were measured after two different detonations having distinct energy, pressure, and gas volume characteristics. Explosive‐derived gases (N 2 O, CO 2 ) and excess radiogenic 4 He and 40 Ar above atmospheric background are used to identify locations of gas transport and new fracture creation after each detonation. Fracture‐released radiogenic 4 He is used to constrain the volume of newly created fractures with a model of helium release from fracturing. Explosive by‐product gas was observed in multiple locations both near and distal to the shot locations for both detonations. Radiogenic 4 He and 40 Ar release from rock damage was observed in locations near the detonation after the second, more powerful detonation. Observed 4 He response is consistent with a model of diffusive release from newly created fractures. Volume of new fractures estimated from the 4 He release ranges from 1 to 5 m 2 with apertures ranging from 0.1 to 1 m. Our results provide evidence that radiogenic noble gases released during fracture creation can be identified at the field scale in real time and used to identify timing and location of fracture creation during deformation events. This technique could be useful in subsurface science and engineering problems where the location and amount of newly created rock fracturing is of interest including fault rupture, mine safety, subsurface detonation monitoring and reservoir stimulation.

58 GEOSCIENCES↗
DOE OSTI2016

Comparison of the telegraph and hyperdiffusion approximations in cosmic-ray transport

The telegraph equation and its generalizations have been repeatedly considered in the models of diffusive cosmic-ray transport. Yet the telegraph model has well-known limitations, and analytical arguments suggest that a hyperdiffusion model should serve as a more accurate alternative to the telegraph model, especially on the timescale of a few scattering times. We present a detailed side-by-side comparison of an evolving particle density profile, predicted by the telegraph and hyperdiffusion models in the context of a simple but physically meaningful initial-value problem, compare the predictions with the solution based on the Fokker–Planck equation, and discuss the applicability of the telegraph and hyperdiffusion approximations to the description of strongly anisotropic particle distributions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Ratcheting Up The Search for Dark Matter

The last several years have included remarkable advances in two of the primary areas of fundamental particle physics: the search for dark matter and the discovery of the Higgs boson. This dissertation will highlight some contributions made on the forefront of these exciting fields. Although the circumstantial evidence supporting the dark matter hypothesis is now almost undeniably significant, indisputable direct proof is still lacking. As the direct searches for dark matter continue, we can maximize our prospects of discovery by using theoretical techniques complementary to the observational searches to rule out additional, otherwise accessible parameter space. In this dissertation, I report bounds on a wide range of dark matter theories. The models considered here cover the spectrum from the canonical case of self-conjugate dark matter with weak-scale interactions, to electrically charged dark matter, to non-annihilating, non-fermionic dark matter. These bounds are obtained from considerations of astrophysical and cosmological data, including, respectively: diffuse gamma ray photon observations; structure formation considerations, along with an explication of the novel local dark matter structure due to galactic astrophysics; and the existence of old pulsars in dark-matter-rich environments. I also consider the prospects for a model of neutrino dark matter which has been motivated by a wide set of seemingly contradictory experimental results. In addition, I include a study that provides the tools to begin solving the speculative ``inverse'' problem of extracting dark matter properties solely from hypothetical nuclear energy spectra, which we may face if dark matter is discovered with multiple direct detection experiments. In contrast to the null searches for dark matter, we have the example of the recent discovery of the Higgs boson. The Higgs boson is the first fundamental scalar particle ever observed, and precision measurements of the production and decay of the Higgs boson represent a unique entry p! oint to searches for new kinds of physics. Continuing to refine our understanding of the Higgs boson will also allow us to learn about a vast array of possible new physics. This dissertation includes work parameterizing some of the scenarios that are most likely to be discovered with future Higgs data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

TOGA: A TOUGH code for modeling three-phase, multi-component, and non-isothermal processes involved in CO 2 -based Enhanced Oil Recovery

TOGA is a numerical reservoir simulator for modeling non-isothermal flow and transport of water, CO 2 , multicomponent oil, and related gas components for applications including CO 2 -enhanced oil recovery (CO 2 -EOR) and geologic carbon sequestration in depleted oil and gas reservoirs. TOGA uses an approach based on the Peng-Robinson equation of state (PR-EOS) to calculate the thermophysical properties of the gas and oil phases including the gas/oil components dissolved in the aqueous phase, and uses a mixing model to estimate the thermophysical properties of the aqueous phase. The phase behavior (e.g., occurrence and disappearance of the three phases, gas + oil + aqueous) and the partitioning of non-aqueous components (e.g., CO 2 , CH 4 , and n-oil components) between coexisting phases are modeled using K-values derived from assumptions of equal-fugacity that have been demonstrated to be very accurate as shown by comparison to measured data. Models for saturated (water) vapor pressure and water solubility (in the oil phase) are used to calculate the partitioning of the water (H 2 O) component between the gas and oil phases. All components (e.g., CO 2 , H 2 O, and n hydrocarbon components) are allowed to be present in all phases (aqueous, gaseous, and oil). TOGA uses a multiphase version of Darcy’s Law to model flow and transport through porous media of mixtures with up to three phases over a range of pressures and temperatures appropriate to hydrocarbon recovery and geologic carbon sequestration systems. Transport of the gaseous and dissolved components is by advection and Fickian molecular diffusion. New methods for phase partitioning and thermophysical property modeling in TOGA have been validated against experimental data published in the literature for describing phase partitioning and phase behavior. Flow and transport has been verified by testing against related TOUGH2 EOS modules and CMG. The code has also been validated against a CO 2 -EOR experimental core flood involving flow of three phases and 12 components. Results of simulations of a hypothetical 3D CO 2 -EOR problem involving three phases and multiple components are presented to demonstrate the field-scale capabilities of the new code. This user guide provides instructions for use and sample problems for verification and demonstration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

A theory of self-organized zonal flow with fine radial structure in tokamak

The (low frequency) zonal flow-ion temperature gradient (ITG) wave system, constructed on Braginskii's fluid model in tokamak, is shown here to be a reaction-diffusion-advection system; it is derived by making use of a multiple spatiotemporal scale technique and two-dimensional (2D) ballooning theory. For real regular group velocities of ITG waves, two distinct temporal processes, sharing a very similar meso-scale radial structure, are identified in the nonlinear self-organized stage. The stationary and quasi-stationary structures reflect a particular feature of the poloidal group velocity. The equation set posed to be an initial value problem is numerically solved for JET low mode parameters; the results are presented in several figures and two movies that show the spatiotemporal evolutions as well as the spectrum analysis—frequency-wave number spectrum, auto power spectrum, and Lissajous diagram. This approach reveals that the zonal flow in tokamak is a local traveling wave. For the quasi-stationary process, the cycle of ITG wave energy is composed of two consecutive phases in distinct spatiotemporal structures: a pair of Cavitons growing and breathing slowly without long range propagation, followed by a sudden decay into many Instantons that carry negative wave energy rapidly into infinity. A spotlight onto the motion of Instantons for a given radial position reproduces a Blob-Hole temporal structure; the occurrence as well as the rapid decay of Caviton into Instantons is triggered by zero-crossing of radial group velocity. A sample of the radial profile of zonal flow contributed from 31 nonlinearly coupled rational surfaces near plasma edge is found to be very similar to that observed in the JET Ohmic phase [J. C. Hillesheim et al., Phys. Rev. Lett. 116, 165002 (2016)]. The theory predicts an interior asymmetric dipole structure associated with the zonal flow that is driven by the gradients of ITG turbulence intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

THE APPLICATION OF MAMMOTH FOR A DETAILED TIGHTLY COUPLED FUEL PIN SIMULATION WITH A STATION BLACKOUT

Accurate calculation of desired quantities to predict fuel behavior requires the solution of interlinked equations representing different physics. Traditional fuels performance codes often rely on internal empirical models for the pin power density and a simplified boundary condition on the cladding edge. These simplifications are performed because of the difficulty of coupling applications or codes on differing domains and mapping the required data. To demonstrate an approach closer to first principles, the neutronics application Rattlesnake and the thermal hydraulics application RELAP-7 were coupled to the fuels performance application BISON under the master application MAMMOTH. A single fuel pin was modeled based on the dimensions of a Westinghouse 17x17 fuel rod. The simulation consisted of a depletion period of 1343 days, roughly equal to three full operating cycles, followed by a station blackout (SBO) event. The fuel rod was depleted for 1343 days for a near constant total power loading of 65.81 kW. After 1343 days the fission power was reduced to zero (simulating a reactor shut-down). Decay heat calculations provided the time-varying energy source after this time. For this problem, Rattlesnake, BISON, and RELAP-7 are coupled under MAMMOTH in a split operator approach. Each system solves its physics on a separate mesh and, for RELAP-7 and BISON, on only a subset of the full problem domain. Rattlesnake solves the neutronics over the whole domain that includes the fuel, cladding, gaps, water, and top and bottom rod holders. Here BISON is applied to the fuel and cladding with a 2D axi-symmetric domain, and RELAP-7 is applied to the flow of the circular outer water channel with a set of 1D flow equations. The mesh on the Rattlesnake side can either be 3D (for low order transport) or 2D (for diffusion). BISON has a matching ring structure mesh for the fuel so both the power density and local burn up are copied accurately from Rattlesnake. At each depletion time step, Rattlesnake calculates a power density, fission density rate, burn-up distribution and fast flux based on the current water density and fuel temperature. These are then mapped to the BISON mesh for a fuels performance solve. BISON calculates the fuel temperature and cladding surface temperature based upon the current power density and bulk fluid temperature. RELAP-7 then calculates the fluid temperature, water density fraction and water phase velocity based upon the cladding surface temperature. The fuel temperature and the fluid density are then passed back to Rattlesnake for another neutronics calculation. Six Picard or fixed-point style iterations are preformed in this manner to obtain consistent tightly coupled and stable results. For this paper a set of results from the detailed calculation are provided for both during depletion and the SBO event. We demonstrate that a detailed calculation closer to first principles can be done under MAMMOTH between different applications on differing domains.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Two noncentrosymmetric polyphosphates featuring infinite one-dimensional (PO{sub 3}){sub ∞} chain, LiMP{sub 2}O{sub 6} (M = Rb, Cs): Synthesis, structure and optical properties

Highlights: • The crystal structure of two noncentrosymmetric LiMP{sub 2}O{sub 6} (M = Rb, Cs) is described in detail. • The linear and nonlinear optical properties for them were characterized firstly. • Band structures and density of states for them have been calculated firstly. - Abstract: Two noncentrosymmetric polyphosphates, namely LiCsP{sub 2}O{sub 6} and LiRbP{sub 2}O{sub 6} have been synthesized from the high temperature solution. The structural analysis suggests that both of them crystallize in orthorhombic space group Fdd2 (No.43), with lattice parameters: a = 18.517(4) Å, b = 19.037(4) Å, c = 13.011(3) Å, and Z = 32 for LiRbP{sub 2}O{sub 6,} while a = 19.009(13) Å, b = 19.402(13) Å, c = 13.180(9) Å, and Z = 32 for LiCsP{sub 2}O{sub 6}. In the structure, the PO{sub 4} polyhedra are connected through corner-sharing to form 1D infinite (PO{sub 3}) {sub ∞} chains stretching along the c axis, which are further linked by the LiO{sub 4} polyhedra to construct a 3D framework with two types of tunnels. And the Rb/Cs atoms are filled in or around the periphery of the tunnels to keep the charge equilibrium. The IR spectra verify their structural validity. In addition, other characterizations including the UV−vis−NIR diffuse reflectance spectra, as well as first-principles theoretical studies have been performed on the two compounds.

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

Effect of drop size on the impact thermodynamics for supercooled large droplet in aircraft icing

Supercooled large droplet (SLD), which can cause abnormal icing, is a well-known issue in aerospace engineering. Although efforts have been exerted to understand large droplet impact dynamics and the supercooled feature in the film/substrate interface, respectively, the thermodynamic effect during the SLD impact process has not received sufficient attention. This work conducts experimental studies to determine the effects of drop size on the thermodynamics for supercooled large droplet impingement. Through phenomenological reproduction, the rapid-freezing characteristics are observed in diameters of 400, 800, and 1300 μm. The experimental analysis provides information on the maximum spreading rate and the shrinkage rate of the drop, the supercooled diffusive rate, and the freezing time. A physical explanation of this unsteady heat transfer process is proposed theoretically, which indicates that the drop size is a critical factor influencing the supercooled heat exchange and effective heat transfer duration between the film/substrate interface. On the basis of the present experimental data and theoretical analysis, an impinging heating model is developed and applied to typical SLD cases. The model behaves as anticipated, which underlines the wide applicability to SLD icing problems in related fields.

42 ENGINEERING↗
DOE OSTI2015

Method for estimating optimal spectral and energy parameters of laser irradiation in photodynamic therapy of biological tissue

We have solved the problem of layer-by-layer laser-light dosimetry in biological tissues and of selecting an individual therapeutic dose in laser therapy. A method is proposed for real-time monitoring of the radiation density in tissue layers in vivo, concentrations of its endogenous (natural) and exogenous (specially administered) chromophores, as well as in-depth distributions of the spectrum of light action on these chromophores. As the background information use is made of the spectrum of diffuse light reflected from a patient's tissue, measured by a fibre-optic spectrophotometer. The measured spectrum is quantitatively analysed by the method of approximating functions for fluxes of light multiply scattered in tissue and by a semi-analytical method for calculating the in-depth distribution of the light flux in a multi-layered medium. We have shown the possibility of employing the developed method for monitoring photosensitizer and oxyhaemoglobin concentrations in tissue, light power absorbed by chromophores in tissue layers at different depths and laser-induced changes in the tissue morphology (vascular volume content and ratios of various forms of haemoglobin) during photodynamic therapy. (biophotonics)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

A Computational Study of the Thermodynamic Conditions Leading to Autoignition in Nanosecond Pulsed Discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results toward resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model (BRM). This can afford the use of very detailed chemical kinetic information. The full nonequilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

42 ENGINEERING↗
DOE OSTI2018

Advanced Modeling of Ions in Solutions, on Surfaces, and in Biological Environments (Final Scientific/Technical Report)

The research supported by this grant focused on modeling liquid water from basic microscopic theory and on developing advanced methodologies for molecular simulation and electronic structure calculation of large and disordered systems. Water is arguably the most important molecule on earth: without it there would be no life, as we know it. Water does not behave as a simple liquid. Understanding its unique and anomalous properties at the molecular level has been the focus of intense study during the entire last century and still is a major goal of research today. To address the problem we adopted a multipronged approach that included methodological and algorithmic developments as well as simulations of important model systems. We made progress on several key issues. (i) We implemented accurate electronic structure schemes for modeling water from quantum mechanical theory that make accessible in simulations processes driven by formation and breaking of chemical bonds. We studied in this way neat bulk water and simple solvated ions. We could understand, in particular, why hydronium, i.e. the positive water ion, diffuses faster than its negative counterpart, the hydroxyl, solving a 200-year-old ionic mystery of water, as reported in a news article featuring our study on Chemistry World, the magazine of the Royal Society of Chemistry (U.K.). (ii) Using advanced statistical mechanical techniques we studied water in the deeply undercooled regime, which may hold the clue to understanding many of the water anomalies. These conditions are currently only accessible to simulations based on empirical force fields. We found that deeply undercooled water, as described by the popular force field ST2, possesses a secondary critical point, below which it would spontaneously separate into two liquids, one of low-density (LDL) and one of high-density (HDL), much in the same way in which water and oil phase-separate at room temperature. The secondary critical point is supported indirectly by many experiments, but rapid ice nucleation has so far hindered direct observation. In related studies we gained new insight on the connection of the two metastable liquids with the two corresponding amorphous ice phases, LDA and HDA, which are observed experimentally. (iii) We developed and implemented several new algorithms for computing the electronic structure of large systems and to accelerate the calculation of electron exchange energies. We applied these methodologies to model nanomaterials, such as nanoribbons and nanoflakes, and to study water absorption and splitting on nanomaterials. (iv) We developed new local order descriptors. (v) In the last year of the project we developed a methodology based on deep neural networks to construct model interatomic potentials that can accurately reproduce ab-initio quantum mechanical simulations at the computational cost of empirical force fields. These potentials will make possible to access the size and time scales that are currently only accessible with empirical force fields without loosing the predictive power of ab-initio simulations.

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

A computational study of the thermodynamic conditions leading to autoignition in nanosecond pulsed discharges

Nanosecond pulsed discharges have attracted the attention of engine manufacturers due to the possibility of attaining distributed ignition sites that accelerate burn rates while resulting in very little electrode erosion. Multidimensional modeling tools currently capture the electrical structure of such discharges accurately, but resolving the chemical structure remains a challenging problem owing to the disparity of time-scales in streamer propagation (nanoseconds) and ignition phenomena (microseconds). The purpose of this study is to extend multidimensional results towards resolving the chemical structure in the wake of streamers (or the afterglow) by using a batch reactor model. This can afford the use of very detailed chemical kinetic information. The full non-equilibrium nature of the electrons is taken into account, along with fast gas heating, shock wave propagation, and thermal diffusion. The results shed light on ignition phenomena brought about by such discharges.

ignition, plasma, kinetics↗