Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphase chemistry”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Integrated low carbon H 2 conversion with in situ carbon mineralization from aqueous biomass oxygenate precursors by tuning reactive multiphase chemical interactions

Meeting our rising demand for clean energy carriers such as H 2 from renewable biomass resources is challenged by the co-emission of CO 2 and CH 4 . To address this challenge, we design novel reactive separation pathways that integrate multiphase chemical reactions by harnessing Ca and Mg bearing minerals as a sorbent to capture CO 2 released during the hydrothermal deconstruction of aqueous biomass oxygenates to produce H 2 and solid carbonates via low temperature aqueous phase reforming and thermodynamically downhill carbon mineralization. Earth abundant catalysts such as Ni/Al 2 O 3 are effective in producing H 2 yields as high as 79% and 74% using ethylene glycol and methanol in the presence of Ca(OH) 2 as an alkaline sorbent, without contaminating or deactivating the catalyst. H 2 yields with in situ carbon mineralization using a Ni or Pt/Al 2 O 3 catalyst are enhanced based on the following order of reactivity: acetate < glycerol < methanol < formate < ethylene glycol. These studies demonstrate that the multiphase chemical interactions can be successfully tuned to enhance H 2 yields through the selective cleavage of C–C bonds using Ni/Al 2 O 3 catalysts to deconstruct biomass oxygenates for producing H 2 and CO 2 , and in situ carbon mineralization by harnessing abundant alkaline materials, as demonstrated using ladle slag. This approach unlocks new scientific possibilities for harnessing multiple emissions including abundant organic-rich wastewater streams and alkaline industrial residues to co-produce low carbon H 2 and carbonate-bearing materials for use in construction by using renewable solar thermal energy resources.

09 BIOMASS FUELS↗
DOE OSTI2024

Revealing the Nature of Binary-Phase on Structural Stability of Sodium Layered Oxide Cathodes

The emergence of layered sodium transition metal oxides featuring a multiphase structure presents a promising approach for cathode materials in sodium-ion batteries, showcasing notably improved energy storage capacity. However, the advancement of cathodes with multiphase structures faces obstacles due to the limited understanding of the integrated structural effects. Herein, the integrated structural effects by an in-depth structure-chemistry analysis in the developed layered cathode system Na x Cu 0.1 Co 0.1 Ni 0.25 Mn 0.4 Ti 0.15 O 2 with purposely designed P2/O3 phase integration, are comprehended. The results affirm that integrated phase ratio plays a pivotal role in electrochemical/structural stability, particularly at high voltage and with the incorporation of anionic redox. In contrast to previous reports advocating solely for the enhanced electrochemical performance in biphasic structures, it is demonstrated that an inappropriate composite structure is more destructive than a single-phase design. The in situ X-ray diffraction results, coupled with density functional theory computations further confirm that the biphasic structure with P2:O3 = 4:6 shows suppressed irreversible phase transition at high desodiated states and thus exhibits optimized electrochemical performance. Finally, these fundamental discoveries provide clues to the design of high-performance layered oxide cathodes for next-generation SIBs.

36 MATERIALS SCIENCE↗
NASA NTRS2024

Critical Load Exceedances for North America and Europe using an Ensemble of Models and an Investigation of Causes for Environmental Impact Estimate Variability: An AQMEII4 Study

Exceedances of critical loads for deposition of sulphur (S) and nitrogen (N) to different ecosystems were estimated using European and North American ensembles of air quality models, under Phase 4 of the Air Quality Model Evaluation International Initiative (AQMEII4), to identify where risk of ecosystem harm is expected to occur based on model deposition estimates. The ensembles were driven by common emissions and lateral boundary condition inputs. Model output was regridded to common North American and Europe 0.125° resolution domains, which were then used to calculate critical load exceedances. New, targeted deposition diagnostics implemented in AQMEII4 allowed an unprecedented level of post-simulation analysis to be carried out and facilitated the identification of specific causes of model-to-model variability in critical load exceedance estimates. New datasets for North American critical loads for acidity for forest soil water and aquatic ecosystems were combined with the ensemble deposition predictions to show a substantial decrease in the area and number of locations in exceedance between 2010 and 2016 (forest soils: 13.2 % to 6.1 %; aquatic ecosystems: 21.2 % to 11.4 %). All models agreed in the direction of the ensemble exceedance change between 2010 and 2016. The North American ensemble also predicted a decrease in both severity and total area in exceedance between the years 2010 and 2016 for eutrophication-impacted ecosystems in the USA (sensitive epiphytic lichen: 81.5 % to 75.8 %). The exceedances for herbaceous community richness also decreased between 2010 and 2016, from 13.9 % to 3.9 %. The uncertainty associated with the North American eutrophication results is high; there were sharp differences between the models in both predictions of total N deposition and the change in N deposition, and hence in the predicted eutrophication exceedances between the two years. The European ensemble was used to predict relatively static exceedances of critical loads with respect to acidification (4.48 % to 4.32 % from 2009 to 2010) while eutrophication exceedance increased slightly (60.2 % to 62.2 %). While most models showed the same changes in critical load exceedances as the ensemble between the two years, the spatial extent and magnitude of exceedances varied significantly between the models. The reasons for this variation were examined in detail by first ranking the relative contribution of different sources of sulphur and nitrogen deposition in terms of deposited mass and model-to-model variability in that deposited mass, followed by their analysis using AQMEII4 diagnostics, along with evaluation of the most recent literature. All models in both the North American and European ensembles had net annual negative biases with respect to observed wet deposition of sulphate, nitrate and ammonium. Diagnostics and recent literature suggest that this bias may stem from insufficient cloud scavenging of aerosols and gases, and may be improved through the incorporation of multiphase hydrometeor scavenging within the modelling frameworks. The inability of North American models to predict the timing of the seasonal peak in wet ammonium ion deposition (observed maximum was in April, while all models predicted a June maximum) may also relate to the need for multiphase hydrometeor scavenging (absence of snow scavenging in all models employed here). High variability in the relative importance of particulate sulphate, nitrate and ammonium deposition fluxes between models was linked to the use of updated particle dry deposition parameterizations in some models. However, recent literature and further development of some of the models within the ensemble suggests these particulate biases may also be ameliorated via the incorporation of multiphase hydrometeor scavenging. Annual sulphur and nitrogen deposition prediction variability was linked to SO 2 and HNO 3 dry deposition parameterizations, and diagnostic analysis showed that the cuticle and soil deposition pathways dominate the deposition mass flux of these species. Further work improving parameterizations for these deposition pathways should reduce variability in model acidifying gas deposition estimates. The absence of base cation chemistry in some models was shown to be a major factor in positive biases in fine mode particulate ammonium and particle nitrate concentrations. Models employing ammonia bidirectional fluxes had both the largest and the smallest magnitude biases, depending on the model and bidirectional flux algorithm employed. A careful analysis of bidirectional flux models suggests that those with poor NH 3 performance may underestimate the extent of NH 3 emissions fluxes from forested areas. Based on these results, an increased process-research focus is therefore recommended for the following model processes and on observations which may assist in model evaluation and improvement: multiphase hydrometeor scavenging combined with updated particle dry deposition, cuticle and soil deposition pathway algorithms for acidifying gases, base cation chemistry and emissions, and NH 3 bidirectional fluxes. Comparisons with satellite observations suggest that oceanic NH 3 emissions sources should be included in regional chemical transport models. The choice of land use database employed within any given model was shown to significantly influence deposition totals in several instances, and employing a common land use database across chemical transport models and critical load calculations is recommended for future work.

critical loads↗
DOE OSTI2021

A New Approach to Fundamental Mechanism Discovery in Polymer Upcycling

We present a new experimental methodology for detailed experimental investigations of depolymerization reactions over solid catalysts. This project aims to address a critical need in fundamental research on chemical upcycling of polymers – the lack of rapid, sensitive, isomerselective probing techniques for the detection of reaction intermediates and products. Our method combines a heterogeneous catalysis reactor for the study of multiphase (gas/polymer melt/solid) systems, coupled to a vacuum UV photoionization time-of-flight mass spectrometer. This apparatus draws on our expertise in probing complex gas-phase chemistry and enables highthroughput, detailed chemical speciation measurements of the gas phase above the catalyst, providing valuable information on the heterogeneous catalytic reactions. Using this approach, we investigated the depolymerization of high-density polyethylene (HDPE) over Ir-doped zeolite catalysts. We showed that the product distribution was dominated by low-molecular weight alkenes with terminal C=C double bonds and revealed the presence of many methyl-substituted alkenes and alkanes, suggesting extensive methyl radical chemistry. In addition, we investigated the fundamental reactivity of model oligomer molecules n-butane and isobutane over ZSM-5 zeolites. We demonstrated the first direct detection of methyl radical intermediates, confirming the key role of methyl in zeolite-catalyzed activation of alkanes. Our results show the potential of this experimental method to achieve deep insight into the complex depolymerization reactions and pave the way for detailed mechanistic studies, leading to increased fundamental understanding of key processes in chemical upcycling of polymers.

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

CO and [C ii] line emission of molecular clouds: the impact of stellar feedback and non-equilibrium chemistry

We analyse synthetic 12 CO, 13 CO, and [C ii] emission maps of molecular cloud (MC) simulations from the SILCC-Zoom project. We present radiation, magnetohydrodynamic zoom-in simulations of individual clouds, both with and without radiative stellar feedback, forming in a turbulent multiphase interstellar medium following on-the-fly the evolution of e.g. H 2 , CO, and C + . We introduce a novel post-processing routine based on cloudy which accounts for higher ionization states of carbon due to stellar radiation in H ii regions. Synthetic emission maps of [C ii] in and around feedback bubbles show that the bubbles are largely devoid of [C ii], as recently found in observations, which we attribute to the further ionization of C+ into C 2+ . For both 12 CO and 13 CO, the cloud-averaged luminosity ratio, $L_\rm {CO}/L_\rm {[C\, \small {II}]}$, can neither be used as a reliable measure of the H 2 mass fraction nor of the evolutionary stage of the clouds. We note a relation between the $I_\rm {CO}/I_\rm {[C\, \small {II}]}$ intensity ratio and the H 2 mass fraction for individual pixels of our synthetic maps. The scatter, however, is too large to reliably infer the H 2 mass fraction. Finally, the assumption of chemical equilibrium overestimates H2 and CO masses by up to 150 and 50 per cent, respectively, and $L_\rm {CO}$ by up to 60 per cent. The masses of H and C + would be underestimated by 65 and 30 per cent, respectively, and $L_\rm {[C\, \small {II}]}$ by up to 35 per cent. Hence, the assumption of chemical equilibrium in MC simulations introduces intrinsic errors of a factor of 2 in chemical abundances, luminosities, and luminosity ratios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Iodide oxidation by ozone at the surface of aqueous microdroplets

The kinetics of iodide oxidation by ozone are studied in single microdroplets to explore chemical reactivity at the air–water interface. Molecular simulations and kinetic modeling are employed to understand the underlying multiphase mechanism.

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

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

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

Dynamic Binary Complexes (DBC) as Super-Adjustable Viscosity Modifiers for Hydraulic Fracturing Fluids

In the preceding project year two, we refined three DBC formulations from a selection of over 50 different chemistries. The optimization study primarily encompassed testing for (i) reversibility extent, (ii) performance in the presence of chemical additives, (iii) adhesion and friction behavior during displacement in wellbores and pipelines, (iv) corrosion protection performance, and (v) injection performance with model fracture systems at the laboratory scale. Highly promising results obtained from all these tests signify the significant potential of DBCs in enhancing hydrocarbon recovery from unconventional reservoirs. The primary activities in the third project year included publishing experimental findings across multiple articles and conducting outreach initiatives. Throughout the year, we undertook tasks such as replicating experimental results, further optimizing various formulations and their associated experimental sets, and conducting additional tests to address missing components based on reviewer feedback and suggestions. We also explored the surfactant and friction-reduction aspects of selected formulations through drag reduction tests. In addition, we constructed an improved fracturing performance setup and performed flow injection tests. The specific DBC formulations focused on during this project period were A8/B1, A12/B5, and A10/B12. We also obtained results for additional DBC formulations and a select few commercial fracturing fluids for the purpose of comparison. Within the project's scope, we aim to enhance the experimental findings with the development of various models. The first two years focused on two key aspects: (i) the creation of a high-fidelity hydraulic fracturing model for non-Newtonian fluids to gain insights into the implementation of DBC fluids in fracking environments, and (ii) the development of a multiphase flow simulator for estimating total production, fluid saturation in the reservoir, and the creation of a fracture propagation model and kinetic Monte Carlo (kMC) models for diverse applications. In the third year, we delved into the fundamental nanostructural properties of DBCs, exploring aspects such as material chemistry, pH tunability, and control of DBC formation and stability. Subsequently, in the extension year, we conducted a systematic investigation of various building blocks containing primary, secondary, and tertiary amine functional groups to understand their impact on rheological and viscoelastic properties. Furthermore, we explored a Dissipative Particle Dynamics (DPD) model to simulate self-assembly processes with precision, creating a high-fidelity representation of relevant nanostructures. The tasks performed this year with the significant results obtained have been discussed in Section 2.

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

Transient simulation of biomass combustion in a circulating fluidized bed riser

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

09 BIOMASS FUELS↗
DOE OSTI2025

A sticky situation: The heterogeneity and implications of magma mush unlocking

Magma, partially molten rock below Earth's surface, can be made up of both fluid (melt) and particles (crystals). The relative fraction of the fluid and particles matters for the magma's mobility and ability to reach Earth's surface. Here we conduct physical laboratory experiments to examine the mechanisms through which a magma becomes mobile and the speed at which pure melt and infiltrate a mixture of fluids and particles. The results have implications for understanding magma mobility and the mechanisms that move magma to Earth's surface. Poster presented at the 2025 International Association of Volcanology and Chemistry of Earth's Interior Meeting (IAVCEI).

multiphase flows↗
DOE OSTI2020

Inertial dynamics of an interface with interfacial mass flux: Stability and flow fields’ structure, inertial stabilization mechanism, degeneracy of Landau’s solution, effect of energy fluctuations, and chemistry-induced instabilities

This work focuses on the long-standing problem of inertial dynamics of an interface with interfacial mass flux and reports new mechanisms for the interface stabilization and destabilization. The interface is a phase boundary separating fluids of different densities and having interfacial mass flux. To analyze the interface dynamics from a far field, we develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations and analyze the interplay of interface stability with flow fields’ structure by directly linking rigorous mathematical attributes to physical observables. We find that the interface is stable when the dynamics conserves the fluxes of mass, momentum, and energy; the stabilization is due to an inertial mechanism causing small oscillations of the interface velocity. In the classic Landau’s dynamics, the postulate of perfect constancy of the interface velocity leads to the development of Landau–Darrieus instability. This destabilization is also linked to the imbalance of the perturbed energy at the interface. The classic Landau’s solution is found to have degeneracy; lifting of the degeneracy may lead to singularity and self-similar dynamics. Our results compare well with traditional theories of combustion and propose new experiments to study the dynamics of the interface and the flow fields in combustible systems. We further conduct reactive molecular dynamics simulations to elucidate the complexity of chemical processes, to study the destabilizing effect of energy fluctuations on the interface stability, and to illustrate the chemistry-induced instabilities. In summary, we identify the extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including the formal properties of fundamental solutions and the qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

42 ENGINEERING↗
DOE OSTI2024

The Scientific Case for Concurrent Neutron and X-ray Scattering and Spectroscopy

The interrogation of materials with X-rays or neutrons to determine the structure, energetics, and dynamics of materials is fundamental to advancing materials' physical and chemical science and developing innovative material technologies. A transcending challenge in developing novel materials is that progress hinges on understanding the structure and dynamics across multiple time and length scales in complex materials that feature multiple components, interfaces, and compositions. Despite the ever-growing demands on materials’ characterization, existing approaches are almost exclusively based on isolated X-ray or neutron scattering, i.e., an approach commensurate with the more narrowly defined needs of fifty years ago. A three-day workshop sponsored by the U.S. National Science Foundation (NSF) analyzed the demand for concurrent neutron and X-ray (NeX) experiments. It was held at the Spring Hill Suites, San Jose, California, from June 2 to 4, 2022. In this workshop, 70 national and international experts ascertained the crucial need to establish NeX capabilities to advance the science of complex materials and systems in the US. Here, we illustrate the need for NeX scattering and spectroscopy experiments by showcasing examples that span areas as diverse as biomaterials, energy science, soft matter, and nanomaterials. To provide NeX capability will require new instrumentation that enables concurrent experiments. Affected areas include chemistry, soft matter, quantum materials, pure and applied chemistry, bioscience, geoscience, and applied materials. NeX benefits research outcomes due to the complementarity of the two techniques, which is essential for better model refinement. While joint refinement of data from separate neutron and X-ray experiments is critical to avoid ambiguities, especially in multiphase-multicomponent materials, concurrent experiments overcome scientific and technical barriers associated with single measurements, separated by location and, thus, time. Among all the examples, these factors introduce uncertainties in the results that complicate data analysis. [1,2] [3] While models are strongly sample-dependent, the principles of joint refinement are generally applicable to these disciplines, including the development of advanced parameterization, modeling, and analysis techniques that also consider the temporal and spatial resolutions of the two methods, leading to unambiguous data interpretation. Solutions for technical barriers must be found to realize NeX experiments, including developing robust sample environments that meet the optical requirements of neutrons and X-rays.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular Simulation of Lithium Carbonate Reactive Vapor–Liquid Equilibria Using a Deep Potential Model

We developed a first-principles machine learning model for the reactive vapor–liquid phase behavior of molten Li 2 CO 3 . The model was trained on ab initio electronic density functional theory data using the Deep Potential (DP) methodology, and its accuracy was evaluated by comparing model predictions of density and viscosity to experimental measurements. Direct coexistence simulations with the DP model over time scales of tens of nanoseconds were used to observe equilibrium dissociation of Li 2 CO 3 into CO 2 residing primarily in the vapor phase and Li 2 O which remains dissolved in the liquid. The simulations covered a range of temperatures, overall system sizes, and vapor-to-liquid volume ratios. Results were analyzed in terms of the observed chemical composition of the liquid and vapor phases, product structure, and CO 2 partial pressures. In addition, we calculated equilibrium constants for the dissociation reaction by assuming ideal-solution behavior for the liquid. As expected on the basis of thermodynamic arguments and prior experiments for this system, the observed partial pressure of CO 2 in the gas phase depends on both the temperature and the ratio of vapor to liquid volumes, while the calculated equilibrium constants only depend on temperature. DP model predictions for the equilibrium constant of the reaction are generally consistent with the available experimental measurements. Furthermore, the present study establishes the validity of the DP methodology for the description of reactive, multiphase equilibria from first principles, with possible applications to many other systems of scientific and technological interest even in the absence of relevant experimental measurements.

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

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
DOE OSTI2023

Impact of organic solutes on capillary phenomena in water-CO2-quartz systems

Hypothesis: The migration of supercritical CO2 (scCO2) injected into underground reservoirs as part of carbon capture and storage is influenced by organic contamination affecting mineral wettability. Molecular dynamics (MD) simulations of relevant systems that incorporate representative organic solutes allow detailed investigation of changes in fundamental interfacial and capillary properties. Experiments: We use MD simulations to explore the effects of four organic solutes (quinoline, decanoic acid, coronene, sorgoleone) on the wettability of quartz by water in the presence of scCO2. We examine the impacts of polar, alkyl, and aromatic moieties as well as fluid flow velocity at elevated temperatures and pressures. Findings: Organic molecules accumulate at the water-CO2 interface, where they distribute according to their size and functional groups. Certain organics penetrate the adsorbed water film at the quartz-CO2 interface, revealing two modes of hydrogen bonding between polar organic functional group, water, and quartz surface –OH groups. Interfacial energies and contact angles are minimally impacted by organic adsorption at the water-CO2 interface, possibly due to simultaneous CO2 desorption. Non- equilibrium MD simulations reveal that flow-induced redistribution of organic compounds modulates the radii of curvature of the advancing and receding water-CO2 interfaces. Our results indicate that organic adsorption on water surfaces is likely ubiquitous during multi-phase flow in soils and sedimen- tary rocks, with implications for the mobilization and transport of organic compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗