Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Coupling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamics of homogeneous cavitation with pressure feedback

Theoretical studies of homogeneous cavitation have largely been based on the classical nucleation theory. However, existing cavitation models cannot adequately describe its dynamics at nanosecond timescale, which has been called for in other fields. We develop a model coupling nucleation kinetics with cavity growth and pressure feedback as saturation mechanisms. Our numerical studies reveal the exponential dependence of cavitation characteristics such as saturation cavity density and most probable cavity radius on model parameters: Tolman length and initial pressure. This work also sheds light on the possibility of accurately determining Tolman length, whose value has a large spread in the literature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

CompLaB v1.0: a scalable pore-scale model for flow, biogeochemistry, microbial metabolism, and biofilm dynamics

Abstract. Microbial activity and chemical reactions in porous media depend on the local conditions at the pore scale and can involve complex feedback with fluid flow and mass transport. We present a modeling framework that quantitatively accounts for the interactions between the bio(geo)chemical and physical processes and that can integrate genome-scale microbial metabolic information into a dynamically changing, spatially explicit representation of environmental conditions. The model couples a lattice Boltzmann implementation of Navier–Stokes (flow) and advection–diffusion-reaction (mass conservation) equations. Reaction formulations can include both kinetic rate expressions and flux balance analysis, thereby integrating reactive transport modeling and systems biology. We also show that the use of surrogate models such as neural network representations of in silico cell models can speed up computations significantly, facilitating applications to complex environmental systems. Parallelization enables simulations that resolve heterogeneity at multiple scales, and a cellular automaton module provides additional capabilities to simulate biofilm dynamics. The code thus constitutes a platform suitable for a range of environmental, engineering and – potentially – medical applications, in particular ones that involve the simulation of microbial dynamics.

58 GEOSCIENCES↗
DOE OSTI2014

Ensemble-Based Parameter Estimation in a Coupled General Circulation Model

Parameter estimation provides a potentially powerful approach to reduce model bias for complex climate models. Here, in a twin experiment framework, the authors perform the first parameter estimation in a fully coupled ocean–atmosphere general circulation model using an ensemble coupled data assimilation system facilitated with parameter estimation. The authors first perform single-parameter estimation and then multiple-parameter estimation. In the case of the single-parameter estimation, the error of the parameter [solar penetration depth (SPD)] is reduced by over 90% after ~40 years of assimilation of the conventional observations of monthly sea surface temperature (SST) and salinity (SSS). The results of multiple-parameter estimation are less reliable than those of single-parameter estimation when only the monthly SST and SSS are assimilated. Assimilating additional observations of atmospheric data of temperature and wind improves the reliability of multiple-parameter estimation. The errors of the parameters are reduced by 90% in ~8 years of assimilation. Finally, the improved parameters also improve the model climatology. With the optimized parameters, the bias of the climatology of SST is reduced by ~90%. Altogether, this study suggests the feasibility of ensemble-based parameter estimation in a fully coupled general circulation model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Equation-based languages – A new paradigm for building energy modeling, simulation and optimization

Most of the state-of-the-art building simulation programs implement models in imperative programming languages. This complicates modeling and excludes the use of certain efficient methods for simulation and optimization. In contrast, equation-based modeling languages declare relations among variables, thereby allowing the use of computer algebra to enable much simpler schematic modeling and to generate efficient code for simulation and optimization. We contrast the two approaches in this paper. We explain how such manipulations support new use cases. In the first of two examples, we couple models of the electrical grid, multiple buildings, HVAC systems and controllers to test a controller that adjusts building room temperatures and PV inverter reactive power to maintain power quality. In the second example, we contrast the computing time for solving an optimal control problem for a room-level model predictive controller with and without symbolic manipulations. As a result, exploiting the equation-based language led to 2, 200 times faster solution

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Increased Asian Sulfate Aerosol Emissions Remarkably Enhance Sahel Summer Precipitation

Observational evidence shows that Sahel summer precipitation has experienced a considerable increase since the 1980s, coinciding with significant diverging trends of increased sulfate emissions in Asia and decreased emissions in Europe (dipole pattern of aerosols between Asia and Europe). The decrease in European sulfate aerosols has substantial effects on the Sahel summer precipitation increase, but the corresponding effect of increased Asian sulfate is unknown. Multi-model simulations in the Precipitation Driver and Response Model Intercomparison Project (PDRMIP) show, compared to decreased European aerosols, that increased Asian aerosols similarly enhance the Sahel summer precipitation but with different large-scale atmospheric circulation changes. Further analysis of the Sixth Coupled Model Intercomparison Project (CMIP6) simulations under historical attribution and various emission scenarios reinforces the results about the climate impacts of anthropogenic aerosols and suggests that in future scenarios with strong international cooperation and rapid climate mitigations (SSP2-45), the Sahel drought will be intensified likely due to the decline in Asian aerosol emissions. Our results suggest that Asian anthropogenic aerosols are likely a non-negligible driver of the recent recovery in Sahel precipitation amounts

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Research Priorities to Incorporate Terrestrial-Aquatic Interfaces in Earth System Models. Workshop Report, September 7-9, 2016

The terrestrial-aquatic interface (TAI) is a highly dynamic component of the Earth system, developed from a near balance between terrestrial and aquatic conditions and forming unique processes and community assemblages. Furthermore, TAIs are known to play a critical role in carbon biogeochemical cycling and have the potential to provide major feedbacks to the Earth system (e.g., methane production). However, there is a lack of basic data and multiscale models to adequately describe how a changing climate can influence the key processes (also an unknown) related to Earth system–relevant feedbacks in these unique, ubiquitous ecosystems. In general, Earth system models (ESMs) have excluded TAI ecosystem processes, thus creating tremendous uncertainty as to how TAIs will influence climate feedbacks across a range of scales, spanning from local to global. For example, ESMs represent wetlands very simplistically at best (e.g., the Accelerated Climate Modeling for Energy project) and include only a static fraction of dry land or open water, entirely lacking key processes in these hybrid areas and critical climate and biological feedbacks between land and water. This new area of research, therefore, must integrate a variety of important research topics—plants, soil, hydrology, reactive transport, microbiology, genomics, and modeling—into a systems-level understanding that can be extended to improve predictive modeling capabilities. The goal of the September 2016 Research Priorities to Incorporate Terrestrial-Aquatic Interfaces in Earth System Models Workshop was to engage the scientific community in an open discussion about critical scientific gaps. These gaps and research topics demand immediate field investigations to gather data for representing these important, yet understudied and underrepresented, ecosystems in ESMs. Workshop results will inform the U.S. Department of Energy’s Office of Biological and Environmental Research (BER) as it plans and prepares for future research efforts that address the TAI gap through model-informed and model-inspired field studies. The resulting data will enable iterative refinement of high-resolution, next-generation ESMs. Specifically, the workshop (1) summarized past and current field, process, and modeling TAI research; (2) identified critical sensitivities and uncertainties in the systems; (3) identified key processes, traits, existing data, and environmental variables needed to adequately characterize these systems; and (4) discussed idealized strategies that couple models and experiments to advance the state of the science in TAI modeling, including potential experiments that would test and improve land model fidelity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Phase Field Modeling of Corrosion for Design of Next-Generation Magnesium-Aluminum Vehicle Joints

In this study, WPI and its partners sought to understand galvanic corrosion in Mg-Al friction stir weld (FSW) joints. Unlike fusion welds, FSW is done in the solid state. This limits (but does not eliminate) the formation of brittle Al-Mg intermetallics, resulting in better joint strength and corrosion resistance. This study aimed to understand corrosion behavior of Mg-Al FSW joints by using a Cahn-Hilliard phase field electrochemistry model based on the work of Pongsaksawad et al. Its context challenge problem was an ultra-light door designed by Magna with 6061-Al sheet outer and ZEK100-Mg sheet inner panels, joined by an adhesive and FSW. Unlike most FSW joints, which weld through the softer material into the harder one, the Al sheet’s hem joint around Mg required the FSW tool to through the harder Al into the softer Mg, which was itself a challenge. Despite this challenge, PNNL succeeded in consistently making very strong FSW lap joints, by using a triflute tool and power control. Characterization showed mostly separate Al and Mg regions in the weld, usually with a hook protruding from Al into Mg, which likely contributed to high joint strength. Nanohardness mapping showed higher hardness in thin Mg swirls into the Al nugget, likely due to the intermetallics. The study used the industry standard SAE J2334 Cyclic Corrosion Test (CCT) to simulate corrosion conditions over the life of a vehicle, as well as linear polarization testing and ASTM G71 pitting corrosion testing. Corrosion reactions were: H 2 O + ½O 2 + 2e - → 2OH - at the Al cathode, and either Mg → Mg 2+ + 2e - or Mg + 2OH - → Mg(OH) 2 + 2e - at the Mg anode, with electron transfer through the joint. Mass loss increased with Al section length, supporting the hypothesis that galvanic corrosion in Al-Mg joints is limited by the cathodic reaction. Both pitting corrosion (first anode reaction) and hydroxide film formation (second reaction) were observed. The J2334 CCT test method was slightly modified, adding a chromium solution cleaning step after each week of corrosion testing, in order to remove Mg oxides and hydroxides and accurately measure mass loss over time. That said, FSW joints treated in this way did not did not exhibit significant reduction in lap shear strength of the joints vs. newly welded samples which had not undergone corrosion. Corrosion appeared to consist mainly of pitting in Mg in a way which did not directly affect the bond between the materials where the joints failed. A Cahn-Hilliard phase field model coupled with electrical potential described the electronically mediated galvanic corrosion electrochemical reactions described above in the four-component Al-Mg-H 2 O-(H 2 O+½O 2 ) system. This 2-D model successfully predicted Mg(OH)₂ formation at the anode in some circumstances, and Mg pitting corrosion in others, and showed the correct electric field directions in both cases. The model predicted Mg pit corrosion depth within a factor of two of measured pits.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

The Carbon Dioxide Removal Model Intercomparison Project (CDRMIP): rationale and experimental protocol for CMIP6

The recent IPCC reports state that continued anthropogenic greenhouse gas emissions are changing the climate, threatening severe, pervasive and irreversible impacts. Slow progress in emissions reduction to mitigate climate change is resulting in increased attention to what is called geoengineering, climate engineering, or climate intervention – deliberate interventions to counter climate change that seek to either modify the Earth's radiation budget or remove greenhouse gases such as CO 2 from the atmosphere. When focused on CO 2 , the latter of these categories is called carbon dioxide removal (CDR). Future emission scenarios that stay well below 2 °C, and all emission scenarios that do not exceed 1.5 °C warming by the year 2100, require some form of CDR. At present, there is little consensus on the climate impacts and atmospheric CO 2 reduction efficacy of the different types of proposed CDR. To address this need, the Carbon Dioxide Removal Model Intercomparison Project (or CDRMIP) was initiated. This project brings together models of the Earth system in a common framework to explore the potential, impacts, and challenges of CDR. Here, we describe the first set of CDRMIP experiments, which are formally part of the 6th Coupled Model Intercomparison Project (CMIP6). These experiments are designed to address questions concerning CDR-induced climate reversibility , the response of the Earth system to direct atmospheric CO 2 removal (direct air capture and storage), and the CDR potential and impacts of afforestation and reforestation, as well as ocean alkalinization.>

Keller, David P.↗
DOE OSTI2025

Baseline Climate Variables for Earth System Modelling

The Baseline Climate Variables for Earth System Modelling (ESM-BCVs) are defined as a list of 135 variables which have high utility for the evaluation and exploitation of climate simulations. The list reflects the most frequently used elements of the Coupled Model Intercomparison Project Phase 6 (CMIP6) archive. Successive phases of CMIP have supported strong results in science and substantially influence international climate policy formulation. This paper responds to both interest in exploiting CMIP data standards in a broader range of climate modelling activities and a need to achieve greater clarity about the significance and intention of variables in the CMIP Data Request. As Earth system modelling archives grow in scale and complexity, there are emerging problems associated with weak standardisation at the variable collection level. That is, there are good standards covering how specific variables should be archived, but this paper fills a gap in the standardisation of which variables should be archived. The ESM-BCV list is intended as a resource for ESM intercomparison projects (MIPs) developing requests to enable greater consistency among MIPs and as a reference for modelling centres to enhance consistency within MIPs. Provisional planning for the CMIP7 Data Request exploits the ESM-BCVs as a core element. The baseline variable list includes 98 variables which have modest or minor data volume footprints and could be generated systematically when simulations are produced and archived for exploitation by the World Climate Research Programme (WCRP) community. A further 35 variables are classed as “high volume” and are only suitable for production when the resource implications are justified.

Juckes, Martin [University of Oxford (United Kingd↗
DOE OSTI2015

Problem-free time-dependent variational principle for open quantum systems

Methods of quantum nuclear wave-function dynamics have become very efficient in simulating large isolated systems using the time-dependent variational principle (TDVP). However, a straightforward extension of the TDVP to the density matrix framework gives rise to methods that do not conserve the energy in the isolated system limit and the total system population for open systems where only energy exchange with environment is allowed. These problems arise when the system density is in a mixed state and is simulated using an incomplete basis. Thus, the basis set incompleteness, which is inevitable in practical calculations, creates artificial channels for energy and population dissipation. To overcome this unphysical behavior, we have introduced a constrained Lagrangian formulation of TDVP applied to a non-stochastic open system Schrödinger equation [L. Joubert-Doriol, I. G. Ryabinkin, and A. F. Izmaylov, J. Chem. Phys. 141, 234112 (2014)]. While our formulation can be applied to any variational ansatz for the system density matrix, derivation of working equations and numerical assessment is done within the variational multiconfiguration Gaussian approach for a two-dimensional linear vibronic coupling model system interacting with a harmonic bath.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

DEVELOPING ATMOSPHERIC RETRIEVAL METHODS FOR DIRECT IMAGING SPECTROSCOPY OF GAS GIANTS IN REFLECTED LIGHT. I. METHANE ABUNDANCES AND BASIC CLOUD PROPERTIES

Upcoming space-based coronagraphic instruments in the next decade will perform reflected light spectroscopy and photometry of cool directly imaged extrasolar giant planets. We are developing a new atmospheric retrieval methodology to help assess the science return and inform the instrument design for such future missions, and ultimately interpret the resulting observations. Our retrieval technique employs a geometric albedo model coupled with both a Markov chain Monte Carlo Ensemble Sampler ( emcee ) and a multimodal nested sampling algorithm ( MultiNest ) to map the posterior distribution. This combination makes the global evidence calculation more robust for any given model and highlights possible discrepancies in the likelihood maps. As a proof of concept, our current atmospheric model contains one or two cloud layers, methane as a major absorber, and a H{sub 2}–He background gas. This 6-to-9 parameter model is appropriate for Jupiter-like planets and can be easily expanded in the future. In addition to deriving the marginal likelihood distribution and confidence intervals for the model parameters, we perform model selection to determine the significance of methane and cloud detection as a function of expected signal-to-noise ratio in the presence of spectral noise correlations. After internal validation, the method is applied to realistic spectra of Jupiter, Saturn, and HD 99492c, a model observing target. We find that the presence or absence of clouds and methane can be determined with high confidence, while parameter uncertainties are model dependent and correlated. Such general methods will also be applicable to the interpretation of direct imaging spectra of cloudy terrestrial planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution

Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Late Holocene Cliff Retreat in Del Mar, CA, Revealed From Shore Platform 10 Be Concentrations and Numerical Modeling

Rocky coast cliff retreat presents a hazard to coastal communities and infrastructure that is potentially amplified under rising sea level conditions, among other factors. Unfortunately, constraints on retreat rates are typically limited to those derived from imagery and maps spanning the last ~100 years. Here, we use a newly developed coupled model of shore platform profile evolution and cosmogenic radionuclide production that considers the influence of relative sea level (RSL) rise, weathering, material resistance, and wave height decay to model cliff retreat over millennial timescales and its potential drivers in Del Mar, California. We demonstrate the ability to use topographic and bathymetric measurements from a narrow shore platform along with a limited data set of nine cosmogenic 10Be concentrations extending ~125 m from the cliff base to obtain modeled cliff retreat rates that steadily range from 5.0 to 12.5 cm yr -1 over the last two millennia until 100 years before present. These rates are consistent with modern retreat rates of about 2–19 cm yr -1 here. RSL rise in Southern California remained relatively constant during the late Holocene, potentially explaining the relatively stable modeled cliff retreat rate over this time. We also explore the relative influence of weathering, material resistance, and wave erosion efficacy and find that both weathering and wave-driven erosion are necessary to replicate the measured data at this location, with the latter exerting a stronger control on model acceptance, suggesting that waves may provide a possible mechanism by which RSL rise may influence coastal cliff erosion in southern California.

58 GEOSCIENCES↗
DOE OSTI2019

Exergy surface shaping and thermodynamic flow control of electro-mechanical-thermal systems

Our work extends the concepts and tools of Hamiltonian Surface Shaping and Power Flow Control (HS SPFC) for electro-mechanical (EM) systems(i.e., adiabatic irreversible work processes and Hamiltonian natural systems)to Exergy Surface Shaping and Thermodynamic Flow Control (ESSTFC) for electro-mechanical-thermal (EMT) systems (i.e., irreversible work processes with heat and mass flows). The extension of HSSPFC requires the development of exergy potential functions, irreversible entropy production terms of the entropy balance equation to obtain the exergy destruction terms for inclusion in the exergy balance equation, and variational principles for producing consistent equations of motion for coupled EMT systems. The Hamiltonian for natural EM systems is an exergy potential function which leaves the development of exergy potential functions for the thermal part of the coupled models. This development is completed by integrating the exergy function over the control volume subject to the modeling assumptions. The irreversible entropy production terms are the exergy destruction terms of the exergy balance equation and the generalization of the mechanical dissipation and electrical resistance within EM systems. These generalized dissipation terms enable the derivation of a consistent set of coupled equations of motion for EMT systems. For this paper, Extended Irreversible Thermodynamics will be utilized to produce consistent thermal equations of motion that directly include the exergy destruction terms. There are several variational principles that are available for application to EMT systems. We focus on the variational principles developed by Biot and Fung [1, 2]. Furthermore, a simplified EMT system that models the EMT dynamics of a Navy ship equipped with a railgun is used to demonstrate the application of ESSTFC for designing high performance, stable nonlinear controllers for EMT systems.

42 ENGINEERING↗
DOE OSTI2022

Heat and Mass Transfer in Hygroscopic Hydrogels

Hygroscopic hydrogels have emerged as a scalable material capable of high-performance vapor sorption for atmospheric water harvesting, dehumidification, and passive cooling due to their fast kinetics and high water uptake. Despite extensive research interest, there is a lack of understanding of the governing energy and mass transport mechanisms within the hydrogels, which hinders their optimization. Here we present a model that explains sorption and desorption in hygroscopic hydrogels for the first time. Our model couples vapor transport in the hydrogel micropores, liquid transport in the polymer nanopores, and heat transport in the porous hydro-gel to accurately capture experimentally observed sorption and desorption data. The knowledge obtained from this model can therefore guide the design of hygroscopic hydrogels and enable improved performance in multi-ple sorption-based applications.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Modeling a Typical Winter-time Dust Event over the Arabian Peninsula and the Red Sea

We used WRF-Chem, a regional meteorological model coupled with an aerosol-chemistry component, to simulate various aspects of the dust phenomena over the Arabian Peninsula and Red Sea during a typical winter-time dust event that occurred in January 2009. The model predicted that the total amount of emitted dust was 18.3 Tg for the entire dust outburst period and that the two maximum daily rates were ~2.4 Tg/day and ~1.5 Tg/day, corresponding to two periods with the highest aerosol optical depth that were well captured by ground- and satellite-based observations. The model predicted that the dust plume was thick, extensive, and mixed in a deep boundary layer at an altitude of 3-4 km. Its spatial distribution was modeled to be consistent with typical spatial patterns of dust emissions. We utilized MODIS-Aqua and Solar Village AERONET measurements of the aerosol optical depth (AOD) to evaluate the radiative impact of aerosols. Our results clearly indicated that the presence of dust particles in the atmosphere caused a significant reduction in the amount of solar radiation reaching the surface during the dust event. We also found that dust aerosols have significant impact on the energy and nutrient balances of the Red Sea. Our results showed that the simulated cooling under the dust plume reached 100 W/m2, which could have profound effects on both the sea surface temperature and circulation. Further analysis of dust generation and its spatial and temporal variability is extremely important for future projections and for better understanding of the climate and ecological history of the Red Sea.

Arabian↗
DOE OSTI2021

Evaluation of CMIP6 models in simulating the statistics of extreme precipitation over Eastern Africa

We report the Eastern Africa region experiences frequent extreme precipitation events that can cause destruction of property and environment, and loss of lives. Thus, there is a need to understand how these events may change in the future and how well the global climate models that are used to make projections can simulate precipitation extremes in this region before they can be used in downscaling or flood and drought impact assessment studies. In this work, we evaluated the ability of sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models to simulate present-day precipitation extremes over the Eastern Africa region during the two rainy seasons (March–May and September–November). We used nine extreme precipitation indices (including seven (one) indices of wet (dry) extremes) defined by the Expert Team on Climate Change Detection and Indices. The CMIP6 models were evaluated against two gridded observation datasets: Global Precipitation Climatology Project One-Degree Daily Dataset and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis 3B42. Three model performance metrics (percentage bias, normalized root-mean-square error, and pattern correlation coefficient) were employed to further assess the strengths and weakness of the models. Our results show that the multi-model ensemble mean generally provides a better representation of observed precipitation and related extremes compared to individual models when considering all metrics and seasons. Several consistent biases are evident across CMIP6 models, which tend to overestimate the total-wet day precipitation and consecutive wet days, and underestimate very wet days and maximum 5-day precipitation in both seasons. Furthermore, no single model consistently performs best, model performance varies with the season and index under consideration and is generally independent of horizontal resolution.

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

Physics informed deep neural network embedded in a chemical transport model for the Amazon rainforest

Secondary organic aerosols (SOA) are fine particles in the atmosphere, which interact with clouds, radiation and affect the Earth’s energy budget. SOA formation involves chemistry in gas phase, aqueous aerosols, and clouds. Simulating these chemical processes involve solving a stiff set of differential equations, which are computationally expensive steps for three-dimensional chemical transport models. Deep neural networks (DNNs) are universal function approximators that could be used to represent the complex nonlinear changes in aerosol physical and chemical processes; however, key challenges such as generalizability to extended time periods, preservation of mass balance, simulating sparse model outputs, and maintaining physical constraints have limited their use in atmospheric chemistry. Here, we develop an approach of using a physics-informed DNN that overcomes previous such challenges and demonstrates its applicability for the chemical formation processes of isoprene epoxydiol SOA (IEPOX-SOA) over the Amazon rainforest. The DNN is trained with data generated by simulating IEPOX-SOA over the entire atmospheric column, using the Weather Research and Forecasting Model coupled with Chemistry (WRF-Chem). The trained DNN is then embedded within WRF-Chem to replace the computationally expensive default solver of IEPOX-SOA formation. The trained DNN predictions generalizes well with the default model simulation of the IEPOX-SOA mass concentrations and its size distribution (20 size bins) over several days of simulations in both dry and wet seasons. The embedded DNN reduces the computational expense of WRF-Chem by a factor of 2. Our approach shows promise in terms of application to other computationally expensive chemistry solvers in climate models.

54 ENVIRONMENTAL SCIENCES↗