Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Mixture modeling”

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 883 records · Page 49

Estimation of Daily PM(sub 10) Concentrations in Italy (2006-2012) Using Finely Resolved Satellite Data, Land Use Variables and Meteorology

Health effects of air pollution, especially particulate matter (PM), have been widely investigated. However, most of the studies rely on few monitors located in urban areas for short-term assessments, or land use/dispersion modelling for long-term evaluations, again mostly in cities. Recently, the availability of finely resolved satellite data provides an opportunity to estimate daily concentrations of air pollutants over wide spatio-temporal domains. Italy lacks a robust and validated high resolution spatio-temporally resolved model of particulate matter. The complex topography and the air mixture from both natural and anthropogenic sources are great challenges difficult to be addressed. We combined finely resolved data on Aerosol Optical Depth (AOD) from the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm, ground-level PM10measurements, land-use variables and meteorological parameters into a four-stage mixed model framework to derive estimates of daily PM10concentrations at 1-km2 grid over Italy, for the years 2006-2012. We checked performance of our models by applying 10-fold cross-validation (CV) for each year. Our models displayed good fitting, with mean CV-R2=0.65 and little bias (average slope of predicted VS observed PM10=0.99). Out-of-sample predictions were more accurate in Northern Italy (Po valley) and large conurbations (e.g. Rome), for background monitoring stations, and in the winter season. Resulting concentration maps showed highest average PM10levels in specific areas (Po river valley, main industrial and metropolitan areas) with decreasing trends over time. Our daily predictions of PM10concentrations across the whole Italy will allow, for the first time, estimation of long-term and short-term effects of air pollution nationwide, even in areas lacking monitoring data. Copyright © 2016 Elsevier Ltd. All rights reserved.

Particulate Matter/analysis↗

Experimental and computational investigation of the influence of iso -butanol on autoignition of n -decane and n -heptane in non-premixed flows

Experimental and computational investigation is carried out to elucidate the influence of iso-butanol on critical conditions of autoignition of n-decane and n-heptane. The counterflow configuration is employed. In this configuration an axisymmetric stream of air is directed over the surface of an evaporating pool of a liquid fuel. A stagnation plane is established. The temperature of the air is increased until autoignition is achieved in the mixing layer in the vicinity of the stagnation plane. The temperature of the air stream at autoignition, T ig , is measured at various values of the strain rate, which is defined as the axial gradient of the axial component of the flow velocity at the stagnation plane. Fuels tested are n-decane, n-heptane, iso-butanol and various mixtures of n-decane/iso-butanol and n-heptane/iso-butanol. Kinetic modeling is carried out using the recently updated version of the comprehensive CRECK chemical-kinetic mechanism. Critical conditions of autoignition are predicted for all fuels and fuel mixtures and the results are compared with the measurements. Computations show that low-temperature chemistry plays a significant role in promoting autoignition of n-decane and n-heptane, and the influence of low-temperature chemistry decreases with increasing strain rates because there is insufficient residence for the low-temperature reactions to take place. Experimental data and numerical simulations show that addition of even small amounts of iso-butanol to n-decane or n-heptane increases the value of T ig at low strain rates indicating that iso-butanol strongly inhibits the low-temperature chemistry of n-decane and n-heptane. Furthermore, the simulations show that when iso-butanol is present in the liquid fuel, only a small amount of n-decane is available in the gas phase and its low concentration cannot effectively sustain the low temperature degenerate branching oxidation process. Here, predicted flame structures show that the peak values of mole fraction of ketohydroperoxide are significantly reduced when iso-butanol is added to n-decane indicating that the kinetic pathway to low temperature ignition is blocked. This observation is confirmed from sensitivity analysis.

42 ENGINEERING↗

Understanding and controlling water-organic co-transport in amorphous microporous materials

The transport of molecules in microporous material is a significant and active area of research in separation applications. The movement of vapor/liquid molecules in a high-loading condition in microporous spaces is especially challenging to interrogate. This project aims to understand and control the transport of complex water and organic solvent mixtures in varying structures of microporous carbon molecular sieve (CMS) and activated carbon membranes. Modeling of transport behavior in such carbonaceous samples using computational modeling is difficult due to their amorphous structure. Therefore, this research intends to study the transport mechanism of water-organic mixtures by experimentally revealing fundamental transport properties, such as guest sorption amounts as well as diffusion and permeation rates. Considering the different size and guest-host affinity of water and organic solvent molecules, the objective of this research is to understand the structural conditions within the ultramicropores and micropores that generate different types of molecular transport and selection mechanisms within such complex systems. Carbonaceous materials are developed using tailored pyrolysis techniques to pyrolyze polymeric precursors, including polyvinylidene chloride (PVDC),polyvinylidene fluoride (PVDF), polymer of intrinsic microporosity (PIM) -1, and fully aromatic polyamide. The initial stage of the research will focus on the microscopic diffusion and sorption studies of pure component water, xylene isomers, and n,n- dimethylformamide within various carbons. Additional structural investigation on CMS, such as gas physisorption and neutron scattering studies, will be executed to gain deeper insight into these structure-transport relationships. Scanning electron microscopy and X-ray photoelectron microscopy are also used to further characterize these materials. The research will progress to probe mutual diffusion of complex water-organic mixtures within various CMS microstructures. Competitive sorption and permeation studies of the binary mixtures in different microstructural CMS will be followed. The obtained binary mixture transport parameters were used to predict binary water-organic mixtures transport.

36 MATERIALS SCIENCE↗

A Study of Cavitation-Ignition Bubble Combustion

We present the results of an experimental and computational study of the physics and chemistry of cavitation-ignition bubble combustion (CIBC), a process that occurs when combustible gaseous mixtures are ignited by the high temperatures found inside a rapidly collapsing bubble. The CIBC process was modeled using a time-dependent compressible fluid-dynamics code that includes finite-rate chemistry. The model predicts that gas-phase reactions within the bubble produce CO and other gaseous by-products of combustion. In addition, heat and mechanical energy release through a bubble volume-expansion phase are also predicted by the model. We experimentally demonstrate the CIBC process using an ultrasonically excited cavitation flow reactor with various hydrocarbon-air mixtures in liquid water. Low concentrations (< 160 ppm) of carbon monoxide (CO) emissions from the ultrasonic reactor were measured, and found to be proportional to the acoustic excitation power. The results of the model were consistent with the measured experimental results. Based on the experimental findings, the computational model, and previous reports of the "micro-diesel effect" in industrial hydraulic systems, we conclude that CIBC is indeed possible and exists in ultrasonically- and hydrodynamically-induced cavitation. Finally, estimates of the utility of CIBC process as a means of powering an idealized heat engine are also presented.

Nguyen, Quang-Viet↗

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-resolved measurements of HO 2 radical in a heated plasma flow reactor

Time-resolved, absolute HO 2 number density in diluted H 2 –O 2 –Ar, CH 4 –O 2 -Ar, and C 2 H 4 –O 2 –Ar mixtures excited by a repetitive ns pulse discharge in a heated plasma flow reactor is measured by Cavity Ringdown Spectroscopy (CRDS). The experimental results are obtained at $\textit{T}$ = 300-600 K and $\textit{P}$ = 130 Torr, both during the discharge pulse burst and in the afterglow. In this work, the HO 2 number density is inferred from the CRDS data using a spectral model exhibiting good agreement with previous measurements of absolute HO 2 absorption cross sections. In the room-temperature H 2 –O 2 mixture, as well as in CH 4 –O 2 and C 2 H 4 –O 2 mixtures over the entire temperature range studied, HO2 is generated only during the discharge burst and decays in the afterglow. However, in the H 2 –O 2 mixture at elevated temperatures, $\textit{T}$ = 400-600 K, HO 2 persists in the afterglow up to 10 ms after the discharge burst, comparable with the flow residence time in the reactor. Comparison with kinetic modeling shows that the sustained reactivity after the source of radicals is turned off is due to a chain propagation / hydrogen oxidation process, which dominates the radical recombination reactions. The kinetic modeling predictions are in good agreement with the relative HO 2 number density measured in all three mixtures, although the model underpredicts the absolute number densities in H 2 –O 2 at $\textit{T}$ = 400-600 K by up to a factor of two. Detection of the sustained low-temperature reactivity in H 2 –O 2 , initiated by the radical generation in the plasma, suggests that the plasma excitation may also affect kinetics of oxidation and reforming of fuels exhibiting low-temperature chemistry below hot ignition point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Redox equilibria of multivalent ions in silicate glasses

Experimental studies were made on the compositional dependence of the redox equilibrium of Eu in synthetic silicate liquids, together with an empirical model describing the observed compositional dependence. Electron paramagnetic resonance (EPR) was used to measure the concentration ratio of Eu(2+) to Eu(3+) in various glasses formed by rapidly quenching silicate liquids. The compositional field studied comprised mixtures of SiO2, TiO2, Al2O3, CaO, MgO, and Na2O. The proposed model describes the Eu(2+)/Eu(3+) ratio over the entire compositional field in terms of parameters easily related to each glass composition. The general applicability and utility of the model is further demonstrated by its application to the Fe(2+)-Fe(3+), Ce(3+)-Ce(4+), and Cr(3+)-Cr(6+) redox reactions in binary alkali oxide silicate glasses of Li, Na, and K.

Lauer, H. V., Jr.↗

Synchrophasor Measurements-based Events Detection Using Deep Learning

Deep learning algorithms have been developed for phasor measurement units (PMUs) analysis aiming at providing grid operators to observe and react to significant real-time changes in the grid associated with multiple factors (e.g., power generation and load variations, different type of faults, and equipment mailfunction), or for offline post-event system diagnostics. In this study, a Long Short-Term Memory (LSTM)-based deep neural network (DNN) is adopted and evaluated to identify the most appropriate model configurations for event detection and longer-term anomalous pattern extraction. The proposed DNN model shows the potential on long-term predictions with the ability to capture nonlinear and nonstationary mixture complex patterns in PMU datasets. Real-world PMU in the WECC system were used for model development and validation.

deep learning↗

Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"

Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.

EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗

Comparison of numerical simulation and experimental data for steam-in-place sterilization

A complex problem involving convective flow of a binary mixture containing a condensable vapor and noncondensable gas in a partially enclosed chamber was modelled and results compared to transient experimental values. The finite element model successfully predicted transport processes in dead-ended tubes with inside diameters of 0.4 to 1.0 cm. When buoyancy driven convective flow was dominant, temperature and mixture compositions agreed with experimental data. Data from 0.4 cm tubes indicate diffusion to be the primary air removal method in small diameter tubes and the diffusivity value in the model to be too large.

Young, Jack H.↗

Automation-Accelerated Electrolyte Design Mitigates Solubility Competition between Redox-Active Molecules and Supporting Salts

In nonaqueous redox-flow batteries (NRFBs), redox-active organic molecules (ROMs) and supporting salts compete for solvation sites, limiting achievable energy density. We combine automated high-throughput experimentation (HTE) with camera-based saturation monitoring and quantitative NMR to measure paired (ROM, salt) solubilities across single and mixed organic solvents. Using 2,1,3-benzothiadiazole (BTZ) with lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) as a model system, we find that a binary m-xylene/acetonitrile mixture dissolves ≈3 M of both BTZ and LiTFSI─surpassing the previously reported 2 M ceiling for neat acetonitrile─by leveraging complementary solvation (MX is BTZ-philic and salt-phobic; ACN stabilizes LiTFSI). A random-forest model (RMSE ≈ 0.24) trained on solvent descriptors highlights log P and salt concentration as dominant predictors and predicts MX/ACN ≈0.3/0.7 (v/v) to be near-optimal. These formulations retain practical viscosity and ∼5 mS·cm –1 conductivity at high loading. In conclusion, the workflow provides a reproducible, data-centric route to NRFB electrolyte design and motivates an open, standardized dual-solute solubility resource for accelerated electrolyte discovery.

Electrolytes↗

Water-Triggered Photoinduced Electron Transfer in Acetonitrile–Water Binary Solvent. Solvent Microstructure-Tuned Reactivity of Hydrophobic Solutes

The solvent-composition dependence of quenching triplet states of benzophenone (3 BP) by anisole in acetonitrile-water (ACN-H2O) mixtures was investigated by laser flash photolysis over the water mole fraction (xw) increasing from 0 to 0.92. Single exponential decay of 3 BP was observed over the whole composition range. The quenching rate constant consistently increased with the water content but increased far more rapidly with xw > 0.7. The water-triggered electron-transfer (ET) mechanism was confirmed by a steeply growing quantum yield of the benzophenone ketyl radical anion, escaping back-ET when the partial water volume exceeded the acetonitrile one. The water-content influence on the 3 BP quenching rate was described by a kinetic model accounting for the microheterogeneous structure of the ACN-H2O mixtures and the very different solubility of the reactants in the solvent components. According to the model, the ET mechanism occurs at a rate constant of 1.46 × 109 M-1 s -1 and is presumably assisted by the ACN-H2O hydrogen-bonding interaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Solvent Selection for the Separation of Lignin-Derived Monomers Using the Conductor-like Screening Model for Real Solvents

The separation of desired monomers from a liquid-phase mixture of lignin depolymerization products is necessary to facilitate their upscaling and upgrading for industrial applications. One effective method to separate multiple liquid-phase products is countercurrent chromatography (CCC), which is a common liquid chromatography technique that separates target solutes based on differences in their partitioning in a biphasic solvent system. Effective CCC separation requires the selection of solvent compositions to tune solute partition coefficients. Here, to alleviate the experimental burden of selecting optimal solvent systems, we apply the conductor-like screening model for real solvents (COSMO-RS) method to compute partition coefficients of representative lignin monomers from four currently used depolymerization strategies in standard solvent systems. We further design new ternary and quaternary solvent systems that are predicted to further improve separation efficacy. On the basis of these predicted partition coefficients and empirical solvent selection criteria for CCC measurements, we suggest a range of solvent systems that would be suitable for the effective separation of aromatic lignin-derived products via CCC or similar liquid–liquid extraction methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Photo-sensitized ignition of hydrogen/oxygen mixtures for hypersonic flight vehicles

An approach toward improving supersonic ignition and combustion efficiency was investigated under a feasibility program with NASA-Ames. The idea is to selectively irradiate targetted species in the combustion zone. The ensuing photodissociation reactions produce highly reactive radicals which modify the gas phase kinetics. The overall objective of this research was to demonstrate photosensitized ignition of H2/O2 mixtures, and to gain a preliminary understanding of the ignition mechanism, especially as it pertains to the concept of supersonic combustion. To this end, a series of experiments were performed with a model sensitizer, NO2. At subatmospheric pressures and moderate temperatures (above 100 F), mixtures of H2/O2/NO2 are ignitable by a pulsed eximer laser at 308 nm, 115 to 200 mJ. Ignitability was found to be a function of NO2 pressure, temperature, laser energy per pulse, and laser lens focal length (fluence at focal volume). Mixtures could not be ignited in the absence of the sensitizer. While increasing temperature improved sensitized ignition in the lower temperature range, a retarding effect was observed at the highest temperatures. To more thoroughly understand these thermal effects, the UV-Vis and IR spectra of NO2 and H2/O2/NO2 mixtures were studied over the temperature range 70 to 900 F. Initially, increasing temperature shifts the N2O4 in a reversible reaction to 2 NO2 equilibrium toward NO2, favoring photochemical ignition. However, at still higher temperature, NO2 converts to non-absorbing NO + 1/2 O2, thus preventing photosensitized ignition.

Lavid, Moshe↗

Fuel Spray Mixing and Wall Interaction

All future high-efficiency engines will have fuel directly sprayed into the engine cylinder. Engine developers agree that a major barrier to the rapid development and design of these high-efficiency, clean engines is the lack of accurate fuel spray computational fluid dynamics (CFD) models. The spray injection process largely determines the fuel–air mixture processes in the engine, which subsequently drive combustion and emissions in both direct-injection gasoline and diesel systems, particularly at cold-start conditions when aftertreatment is ineffective. Engines must be tolerant to a range of fuels, and there must be an understanding of how specific fuel properties affect the spray mixing and evaporation processes to intentionally create better fuels and better injectors. More predictive spray combustion models will enable rapid design and optimization of future high-efficiency engines, providing more affordable vehicles and saving fuel.

33 ADVANCED PROPULSION SYSTEMS↗

Gas-Liquid Flow Modeling for Renewable Fuels Production

Aerobic/anaerobic and gas fermentation pathways have emerged as promising new technologies for the generation of renewable fuels/chemicals from biomass derived sugars, and mixtures of greenhouse/energy rich gas streams (CO2/CH4/H2/CO) via microbial action. Example pathways include sugars-to-ethanol conversion, biomethanation (CO2/H2 to CH4), biogas upgrading, CO fermentation and wet-waste conversion. Gas and liquid phase transport, mass-transfer, and mixing physics at large length scales can significantly affect microbial conversion rates, particularly when the microbial reaction requires a narrow set of conditions. These phenomena are difficult to study in small-scale bench-top reactors that are typically well-mixed. Predictive computational fluid dynamics (CFD) based simulations can therefore aid in the scale-up, design and optimization of these reactors. This work presents multiphase Euler-Euler CFD simulations of at-scale (~500 m3) bioreactors. Our mathematical model treats the gas and liquid as interpenetrating phases. This approach reduces the computational complexity of tracking individual gas bubbles that are several orders of magnitude smaller than reactor dimensions. We solve the Reynolds averaged Navier-Stokes (RANS) multiphase equations that account for phase and chemical species transport, interphase mass and momentum transfer and uses a phenomenological model for gas uptake by microbes. We use a customized solver derived from open-source CFD toolbox, OpenFOAM [1], to perform these simulations, which has been validated against small-scale reactors in our previous work [2]. There is currently a knowledge-gap regarding bubble-size distributions when using gas mixtures with vastly different properties, which can have a significant impact overall mass-transfer. For example, hydrogen bubbles are more buoyant compared to other relatively heavier gases (CO2/CH4/CO), resulting in a large distribution of residence times and bubble sizes. This work therefore develops a deeper understanding of bubble dynamics and interphase mass transfer in such heterogenous gas mixtures through well-resolved computational models. We use a population balance model (PBM) for bubble-size-distribution modeling that is validated against small-scale experiments in our solver with an uncertainty quantification study for bubble coalescence and break-up model parameters. Results pertaining to multiple simulations of gas-fermentation reactors are presented where gas mixtures with varying compositions of CO2/CH4/CO/H2 are imposed at the sparger boundaries. The spatio-temporal variations in bubble-size distribution and mass transfer coefficient are analyzed for varying superficial velocities and gas-compositions for varying sizes of bubble-column and airlift reactors. This work will also examine the performance of different reactor designs, viz. bubble column reactor, airlift reactor with an internal draft tube, and a stirred-tank reactor with Rushton impellers. Reactor mass-transfer coefficient, gas hold-up, and dissolved gas distribution are critically analyzed among reactors, and sensitivity studies pertaining to gas flow rates and reactor geometry will be presented. [1] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998. [2] Rahimi, M., Sitaraman, H., Humbird, D. and Stickel, J., Computational fluid dynamics study of full-scale aerobic bioreactors: Evaluation of gas-liquid mass transfer, oxygen uptake, and dynamic oxygen distribution, Chemical Engineering Research and Design, 139: 283-295.

BIOMASS FUELS↗

Electrochemical Performance of Mixed Redox-Active Organic Molecules in Redox Flow Batteries

Designing electrolytes based on mixture of different organic redox active molecules brings the opportunity of enhancing the volumetric energy density of flow batteries and removes the requirement of high solubility for individual organic species in the mixture. Here, in the present work, we conduct computational and experimental analysis to investigate the electrochemical performance of mixed redox-active organic molecules. A zero-dimensional transient model is employed to investigate the changes in the half-cell potential and the concentrations and partial currents of individual redox reactions in a mixture of organic molecules over time. The model demonstrates the effects of individual properties of species such as kinetic rate constants, mass transfer coefficients, concentration ratios and standard redox potentials and reports the effect of energy-losing homogenous chemical redox reaction on the voltage efficiency and concentration ratios of the mixed species. Pairs of anthraquinone negolyte species were selected for an experimental case study. A mixture of 2,6-N-TSAQ and 2,6-DHAQ showed 40% increase in the volumetric energy density compared to the performance of 2,6-DHAQ alone. Based on the results of the experimental and computational analysis, we propose guidelines for the design of suitable mixed redox-active organic species.

25 ENERGY STORAGE↗

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗