Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surrogates”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Physics-Guided Continual Learning for Predicting Emerging Aqueous Organic Redox Flow Battery Material Performance

Aqueous organic redox flow batteries (AORFBs) have gained popularity in renewable energy storage due to their low cost, environmental friendliness and scalability. The rapid discovery of aqueous soluble organic (ASO) redox-active materials necessitates efficient machine learning surrogates for predicting battery performance. The physics-guided continual learning (PGCL) method proposed in this study can incrementally learn data from new ASO electrolytes while addressing catastrophic forgetting issues in conventional machine learning. Using a AORFB database with a thousand potential materials generated by a 780 $\text{cm}^2$ interdigitated cell model, PGCL incorporates AORFB physics to optimize the continual learning task formation and training strategies to retain previously learned battery material knowledge. Finally, the trained PGCL demonstrates its capability in assessing emerging ASO materials within the established parameter space when evaluated with the dihydroxyphenazine isomers.

25 ENERGY STORAGE↗
DOE OSTI2025

Learning to Simulate Aerosol Dynamics with Graph Neural Networks

Aerosol effects on climate, weather, and air quality depend on characteristics of individual particles, which are tremendously diverse and change in time. Particle-resolved models are the only models able to capture this diversity in particle physiochemical properties, and these models are computationally expensive. As a strategy for accelerating particle-resolved microphysics models, we introduce Graph-based Learning of Aerosol Dynamics (GLAD) and use this model to train a surrogate of the particle-resolved model PartMC-MOSAIC. GLAD implements a Graph Network-based Simulator (GNS), a machine learning framework that has been used to simulate particle-based fluid dynamics models. In GLAD, each particle is represented as a node in a graph, and the evolution of the particle population over time is simulated through learned message passing. Here, we demonstrate our GNS approach on a simple aerosol system that includes condensation of sulfuric acid onto particles composed of sulfate, black carbon, organic carbon, and water. A graph with particles as nodes is constructed, and a graph neural network (GNN) is then trained using the model output from PartMC-MOSAIC. The trained GNN can then be used for simulating and predicting aerosol dynamics over time. Results demonstrate the framework's ability to accurately learn chemical dynamics and generalize across different scenarios, achieving efficient training and prediction times. We evaluate the performance across four scenarios, highlighting the framework's robustness and adaptability in modeling aerosol microphysics and chemistry.

aerosol chemistry dynamics↗
DOE OSTI2022

Metabolomics Analysis of Bacterial Pathogen Burkholderia thailandensis and Mammalian Host Cells in Co-culture

The Tier 1 HHS/USDA Select Agent Burkholderia pseudomallei is a bacterial pathogen that is highly virulent when introduced into the respiratory tract and intrinsically resistant to many antibiotics. Transcriptomic- and proteomic-based methodologies have been used to investigate mechanisms of virulence employed by B. pseudomallei and Burkholderia thailandensis, a convenient surrogate; however, analysis of the pathogen and host metabolomes during infection is lacking. Changes in the metabolites produced can be a result of altered gene expression and/or post-transcriptional processes. Thus, metabolomics complements transcriptomics and proteomics by providing a chemical readout of a biological phenotype, which serves as a snapshot of an organism’s physiological state. However, the poor signal from bacterial metabolites in the context of infection poses a challenge in their detection and robust annotation. In this work, we coupled mammalian cell culture-based metabolomics with feature-based molecular networking of mono- and co-cultures to annotate the pathogen’s secondary metabolome during infection of mammalian cells. These methods enabled us to identify several key secondary metabolites produced by B. thailandensis during infection of airway epithelial and macrophage cell lines. Additionally, the use of in silico approaches provided insights into shifts in host biochemical pathways relevant to defense against infection. Using chemical class enrichment analysis, for example, we identified changes in a number of host-derived compounds including immune lipids such as prostaglandins, which were detected exclusively upon pathogen challenge. Taken together, our findings indicate that co-culture of B. thailandensis with mammalian cells alters the metabolome of both pathogen and host and provides a new dimension of information for in-depth analysis of the host–pathogen interactions underlying Burkholderia infection.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Tutorial: Machine-Learning-Based CREASE-2D Analysis of 2D SAXS Profiles to Characterize Anisotropic Nanostructures in Soft Materials

We present a tutorial to guide users on how to extend the Computational Reverse Engineering Analysis of Scattering Experiments-2D (CREASE-2D) framework to interpret their experimental two-dimensional small-angle scattering (SAS) data from soft materials (e.g., polymers, peptide amphiphiles, biomolecular fibrils). Unlike most traditional SAS analysis approaches, which typically rely on azimuthally averaged onedimensional (1D) profiles, CREASE-2D utilizes the complete 2D scattering profile to reveal information about anisotropy in the structure. In past applications, CREASE has provided insights into complex structural features, including the cross-sectional shapes of assembled nanostructures and dispersity in these features, which are difficult to discern with existing analytical models. While (1D- ) CREASE has been applied to SANS and SAXS data, this tutorial shares the steps for implementing CREASE-2D using an example of a dipeptide solution system, for which we have SAXS data. We present details for these steps involved in using CREASE-2D to interpret SAXS profiles: how to preprocess SAXS data, define relevant structural features, generate three-dimensional real-space structures for specific values of these features, train a machine learning (ML) surrogate model to predict scattering profiles for given structural features, and optimize these features using genetic algorithms (GA). Then, we use these steps to interpret complex 2DSAXS data collected from dipeptide solutions that, in microscopy images, exhibit nanoscale structures that could be elliptical tubes/ flat tapes/cylinders or a combination of these cross sections. Open-source codes, computational hardware, and software requirements, as well as the strengths and limitations of this protocol, are also presented. We expect researchers working with (soft) biomaterials, peptide amphiphiles, amphiphilic polymer solutions, polymer nanocomposites, and blends of particles/polymers will find this CREASE-2D method and this tutorial of use.

CREASE↗
DOE OSTI2020

Finite-Temperature Structures of Supported Subnanometer Catalysts Inferred via Statistical Learning and Genetic Algorithm-Based Optimization

Single-atom catalysts (SACs) minimize noble metal utilization and can alter the activity and selectivity of supported metal nanoparticles. However, the morphology of active centers, including single atoms and subnanometer clusters of a few atoms, remains elusive due to experimental challenges. The computational cost to describe numerous cluster shapes and sizes makes direct first-principles calculations impractical. We present a computational framework to enable structure determination for single-atom and subnanometer cluster catalysts. As a case study, we obtained the low energy structures of Pd n (n = 1-21) clusters supported on CeO 2 (111), which are critical components of automobile three-way catalysts. Trained on density functional theory data, a three-dimensional cluster expansion is established using statistical learning to describe the Hamiltonian and predict energies of supported Pdn clusters of any structure. Low energy stable and metastable structures are identified using a Metropolis Monte Carlo-based genetic algorithm in the canonical ensemble at 300 K. We observe that supported single atoms sinter to form bilayer clusters and large cluster isomers share similarities in both shape and energy, and elucidate the significance of the support and microstructure on cluster stability. We discovered a simple surrogate structure-energy model, where the energy per atom scales with the square root of the average first coordination number, which can be used to estimate energies and compare the stability of clusters. Our framework, applicable to any metal/support system, fills an important methodological gap to predict the stability of supported metal catalysts in the subnanometer regime.

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

Improved 140 Nd Production for the 140 Nd/ 140 Pr In Vivo Generator through Target Recycling and Radiochemical Optimization

Theranostic strategies that utilize f-block therapeutic radionuclides, including 161 Tb, 177 Lu, 225 Ac, and 227 Th, suffer from a shortage of positron emission tomography (PET) imaging counterparts in the same chemical space and often rely on 68 Ga as a surrogate. The 140 Nd/ 140 Pr in vivo PET generator, which belongs to the f-block, may address this issue and can be produced via the 141 Pr(p,2n) 140 Nd production route by using medium-energy cyclotrons. However, impurities in the target material, including stable Nd, and the inherent difficulty of adjacent lanthanide separations limit the achievable radionuclidic and chemical purity of 140 Nd. In this work, we address these challenges through the purification and recycling of praseodymium target material and optimization of Nd/Pr separation. The resulting purified 140Nd was evaluated using DOTA and Macropa chelators via radiolabeling and in vitro stability studies. A target material purification and recycling method was developed for the monoisotopic 141 Pr starting material to remove stable Nd impurities, yielding 90.3 ± 4.7% (n = 3) recovery. The purified 141 Pr was isolated as Pr 6 O 11 and irradiated with 24 MeV protons (20.07 MeV at the target surface) at 20 μA for 4 h, which produced 1417.0 ± 83.4 MBq (38.3 ± 2.2 mCi) of 140 Nd at the end of bombardment (EOB). The produced 140 Nd was purified through an optimized DGA normal method to recover 71.6 ± 6.3% pure 140 Nd. The amount of stable Nd reduced progressively in each target purification cycle from >340 ppm without purification to <250 ppb after three cycles, while other measured metallic impurities were below 30 ppb. This improvement in target purity was reflected in the direct increase of apparent molar activity (AMA), when purified 140 Nd was evaluated with DOTA and Macropa chelators. AMA of [ 140 Nd]Nd-DOTA and [ 140 Nd]Nd-Macropa increased from 70.3 MBq/μmol (1.9 mCi/μmol) and 74 MBq/μmol (2.0 mCi/μmol) to 8025.3 MBq/μmol (216.9 mCi/μmol) and 8473.0 MBq/μmol (229.0 mCi/μmol), respectively, after the third target purification cycle. Further evaluation of chelator-labeled 140 Nd showed that [ 140 Nd]Nd-DOTA was stable in phosphate-buffered saline (PBS), saline, human serum, and mouse serum, whereas [140Nd]Nd-Macropa was stable in all except human serum. This work established a practical methodological advance for the production of 140 Nd/ 140 Pr in vivo PET generators, combining optimized target recycling and radiochemical separation to enable scaled-up and high-molar activity 140 Nd suitable for preclinical imaging. These advances support broader development of 140 Nd/ 140 Pr as a robust PET analogue, especially for f-block therapeutics.

Irradiation↗
DOE OSTI2021

Optimal Design of Sustainable Ammonia-Based Food–Energy–Water Systems with Nitrogen Management

As the basis for virtually any form of nitrogen fertilizers, ammonia plays a vital role in agriculture; in addition, there has been an increased interest in its use as a carbon-free energy carrier. However, ammonia is also associated with two major environmental concerns: CO 2 emissions from the conventional production process and nitrogen pollution from the excessive use of ammonia-based fertilizers. To mitigate these environmental impacts, we develop an optimization framework for the design of a sustainable ammonia-based agricultural system that synergistically integrates the production of ammonia from renewable resources and effective measures for nitrogen management. The proposed model captures the effect of intermittency by incorporating both design and detailed operational decisions. Here, by applying a multiscale time representation that reduces the problem size and a tailored surrogate model that accurately approximates model nonlinearity, we are able to achieve optimal solutions within reasonable computation times. A computational case study is conducted using real-world data from a local farm in Morris, Minnesota, and the results indicate the trade-off between cost and nitrogen loss. Importantly, we show that practicing effective nitrogen management can significantly reduce the nitrogen loss with only a small increase in net present cost.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Accurate Thermochemistry of Complex Lignin Structures via Density Functional Theory, Group Additivity, and Machine Learning

A molecular-level understanding of lignin structures and bond dissociation energies could facilitate depolymerization technologies. Still, this information is currently limited due to the lack of databases and the simplification of surrogate models. Here, substitution effects on seven common linkages in lignin polymers are systematically investigated. An automated reaction network generator is employed to create a database of structures. A new group additivity (GA) model based on principal component analysis (PCA) descriptors is introduced and trained on gas-phase density functional theory data of 4100 species at the M06-2X/6-311++G(d,p) level. Hydrogen bonds, local steric, and nonaromatic ring contributions are also incorporated. Lastly, we improve the accuracy of the group additivity model to reach the G4 theory by computing a data set of 770 species at this level and using a data fusion approach.

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

Computational and Experimental Study for the Denitrification of Biomass-Derived Hydrothermal Liquefaction Oil

Hydrothermal liquefaction (HTL) is a promising method for processing wet biomass and waste feedstock to produce biofuels. During the HTL process, proteins and other biomolecules in certain feedstock get converted into nitrogenous compounds in produced biocrude, which represents a major challenge to further upgrading them into fuels. One promising approach is to separate nitrogenous compounds from the biocrude using polymeric resins. In this study, experiments were conducted to down-select sorbent and resin systems using a nitrogen-compound-containing surrogate biocrude. We model the binding interactions between an Amberlyst polymeric resin with various compounds present in the biocrude mixture such as nitrogenous compounds like pyrrole, pyridine, hexanamide, and representative co-existing compounds like phenol and dodecanoic acid. To ascertain the efficiency of various resins in the denitrogenation process, we have developed a quantitative structure–function model for the interacting components in the mixture. Our results suggest that the Amberlyst resin is a viable candidate for efficient removal of target nitrogen-containing compounds (such as pyridine) from the biocrude as a result of favorable interactions. The computational studies provide some insight into how and why the identified resin (Amberlyst) works in selective extraction of nitrogenous compound(s).

09 BIOMASS FUELS↗
DOE OSTI2026

Fe- and Ru-H-Mordenites for Polyethylene Upcycling: Insights from Thermo-Catalytic Pyrolysis and DFT Studies

Catalytic pyrolysis of polyethylene (PE) was performed using Fe- and Ru-impregnated H-mordenite (HM) catalysts to produce lighter hydrocarbons. Catalyst surface areas were analyzed by N 2 adsorption−desorption and acidity by NH 3 - TPD studies. Reducibility of metal from H 2 -TPR, experiments indicate multiple oxidation states of Fe due to higher H 2 consumption. The catalyst activity studies monitored by GC− MS analysis highlight the importance of temperature and PE: catalyst ratio in optimizing hydrocarbon selectivity to C 1 −C 4 alkanes and C 2 −C 5 olefins. Studies with polymer-to-catalyst ratios of 1:1 and 1:2 at higher temperature increased PE conversion with stronger acid sites favoring conversion and moderate acid sites improving selectivity toward lighter olefins. While propane dominated at lower temperatures, propene was the main product at 500 °C. Both catalysts exhibited overall similar conversion, while Ru-HM yielded >40% propene selectivity compared to <40% for Fe-HM. Ethylene selectivity exceeded 10% in both cases. Density functional theory simulations using a C 4 surrogate on Ru-HM confirmed that late-stage cracking dominates dehydrogenation, validating the proposed reaction mechanism.

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

PCET-Based Ligand Limits Charge Recombination with an Ir(III) Photoredox Catalyst

Upon photoinitiated electron transfer, charge recombination limits the quantum yield of photoredox reactions for which the rates for the forward reaction and back electron transfer are competitive. Taking inspiration from a proton-coupled electron transfer (PCET) process in Photosystem II, here a benzimidazole-phenol (BIP) has been covalently attached to the 2,2'-bipyridyl ligand of [Ir(dF(CF 3 )ppy) 2 (bpy)][PF 6 ] (dF(CF 3 )ppy = 2-(2,4-difluorophenyl)-5-(trifluoromethyl)pyridine; bpy = 2,2'-bipyridyl). Excitation of the [Ir(dF(CF 3 )ppy) 2 (BIP-bpy)][PF 6 ] photocatalyst results in intramolecular PCET to form a charge-separated state with oxidized BIP. Subsequent reduction of methyl viologen dication (MV 2+ ), a substrate surrogate, by the reducing moiety of the charge separated species demonstrates that the inclusion of BIP significantly slows the charge recombination rate. The effect of ~24-fold slower charge recombination in a photocatalytic phthalimide ester reduction resulted in a greater than 2-fold increase in reaction quantum efficiency.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Capture Instead of Release: Defect-Modulated Radionuclide Leaching Kinetics in Metal–Organic Frameworks

Comparison of defect-controlled leaching-kinetics modulation of metal–organic frameworks (MOFs) and porous functionalized silica-based materials was performed on the example of a radionuclide and radionuclide surrogate for the first time, revealing an unprecedented readsorption phenomenon. On a series of zirconium-based MOFs as model systems, we demonstrated the ability to capture and retain >99% of the transuranic 241 Am radionuclide after 1 week of storage. Here, we report the possibility of tailoring radionuclide release kinetics in MOFs through framework defects as a function of postsynthetically installed organic ligands including cation-chelating crown ether-based linkers. Based on comprehensive analysis using spectroscopy (EXAFS, UV–vis, FTIR, and NMR), X-ray crystallography (single crystal and powder), and theoretical calculations (nine kinetics models and structure simulations), we demonstrated the synergy of radionuclide integration methods, topological restrictions, postsynthetic scaffold modification, and defect engineering. This combination is inaccessible in any other material and highlights the advantages of using well-defined frameworks for gaining fundamental knowledge necessary for the advancement of actinide-based material development, providing a pathway for addressing upcoming challenges in the nuclear waste administration sector.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Adsorption of Natural Gas in Metal–Organic Frameworks: Selectivity, Cyclability, and Comparison to Methane Adsorption

Evaluation of metal–organic frameworks (MOFs) for adsorbed natural gas (ANG) technology employs pure methane as a surrogate for natural gas (NG). This approximation is problematic, as it ignores the impact of other heavier hydrocarbons present in NG, such as ethane and propane, which generally have more favorable adsorption interactions with MOFs compared to methane. Herein, using quantitative Raman spectroscopic analysis and Monte Carlo calculations, we demonstrate the adsorption selectivity of high-performing MOFs, such as MOF-5, MOF-177, and SNU-70, for a methane and ethane mixture (95:5) that mimics the composition of NG. The impact of selectivity on the storage and deliverable capacities of these adsorbents during successive cycles of adsorption and desorption, simulating the filling and emptying of an ANG tank, is also demonstrated. Finally, the study reveals a gradual reduction in the storage performance of MOFs, particularly with smaller pore volumes, due to ethane accumulation over long-term cycling, until a steady state is reached with substantially degraded storage performance.

03 NATURAL GAS↗
DOE OSTI2026

Hydrosilylation of a Molecular Molybdenum Nitride Provides Mechanistic Insights into Photodriven Ammonia Synthesis from N 2 and H 2

Addition of Ph 2 SiH 2 to [(depe) 2 Mo(N)][BAr F 4 ] (depe = 1,2-bis(diethylphosphino)ethane, BAr F 4 = B(3,5-(CF 3 ) 2 C 6 H 3 ) 4 ) at 60 °C generated the silyl imido molybdenum hydride complex, trans - [(depe) 2 Mo(NSiHPh 2 )H][BAr F 4 ], a surrogate for a proposed intermediate complex in the photodriven hydrogenation to free ammonia. Irradiation of a THF solution of trans -[(depe) 2 Mo(NSiHPh 2 )H]- [BAr F 4 ] with blue light under H 2 produced free amine along with [(depe) 2 MoH 5 ][BAr F 4 ] in 76% yield. This transformation occurred in the absence of a precious metal photocatalyst, suggesting that it was needed only for the initial addition of H 2 to the molybdenum nitride during the first N−H bond-forming step in the photodriven hydrogenation. Deuterium labeling and crossover studies support concerted Si−H bond addition across the Mo≡N bond, enabled by the nucleophilicity of the nitride. Subsequent hydrogenation involves an intramolecular H migration from Mo to the imido ligand, as supported by electronic absorption spectroscopy, transient absorption spectroscopy, initial rate measurements, and deuterium kinetic isotope effect measurements. These findings provide insights into the photodriven hydrogenation of [(depe) 2 Mo(N)][BAr F 4 ] to ammonia and the role of the photocatalyst in this transformation.

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

Development of Nontargeted Workflow of Occupational Exposure by Infrared Ion Spectroscopy and Silicone Wristbands’ Passive Sampling

Monitoring chemical exposure has well-established protocols and procedures for occupational health and industrial hygiene across various career fields and in a variety of occupational and environmental conditions. Iterative analytical development has led to validated methods to provide employers and employees with safe working conditions in terms of industrial hygiene. These validated methods are commonly supported by commercial laboratories. Here, this workflow is explicitly a targeted approach in which sampling media, sample preparation, chemical separation, quantitative target detection, and post processing are well established. This workflow also relies on commercially available chemical surrogates, isotopically labeled internal standards, and certified/standard reference materials from reputable vendors. An untargeted workflow lacks many key ingredients of validated analytical methods to include the known chemical target.

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

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗
DOE OSTI2020

Iron Speciation in PM 2.5 From Urban, Agriculture, and Mixed Environments in Colorado, USA

Atmospheric iron solubility varies depending on whether the particles are collected in rural or urban areas, with urban areas showing increased iron solubility. In this study, we investigate if the iron species present in different environments affects its ultimate solubility. Field data are presented from the Platte River Air Pollution and Photochemistry Experiment (PRAPPE), aimed at understanding the interactions between organic carbon and trace elements in atmospheric particulate matter (PM). 24-hr PM 2.5 samples were collected during the summer and winter (2016–2017), at three different sites on the Eastern Colorado plains: an urban, agricultural, and a mixed site. Downtown Denver had an average total and water-soluble iron air concentration of 181.2 and 7.7 ng m -3 , respectively. Platteville, the mixed site, had an average of total iron of 76.1 ng m -3 , with average water-soluble iron concentration of 9.1 ng m -3 . Jackson State Park (rural/agricultural) had the lowest total iron average of 31.5 ng m -3 and the lowest water-soluble iron average, 1.3 ng m -3 . The iron oxidation state and chemical speciation of 97 samples across all sites and seasons was probed by X-ray absorption near edge structure (XANES) spectroscopy. The most common iron phases observed were almandine (Fe 3 Al 2 Si 3 O 12 ) (Denver 21%, Platteville 16%, Jackson 24%), magnetite (Fe 3 O 4 ) (Denver 9%, Platteville 4%, Jackson 5%) and Fe (III)dextran (Denver 5%, Platteville 13%, Jackson 5%), a surrogate for Fe-organic complexes. Additionally, native iron [Fe(0)] was found in significant amounts at all sites. No correlation was observed between iron solubility and iron oxidation state or chemical speciation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Deep Learning Emulation of Subgrid-Scale Processes in Turbulent Shear Flows

Deep neural networks (DNNs) are developed from a data set obtained from the dynamic Smagorinsky model to emulate the subgrid-scale (SGS) viscosity (ν sgs ) and diffusivity (κ sgs ) for turbulent stratified shear flows encountered in the oceans and the atmosphere. These DNNs predict ν sgs and κ sgs from velocities, strain rates, and density gradients such that the evolution of the kinetic energy budget and density variance budget terms is similar to the corresponding values obtained from the original dynamic Smagorinsky model. These DNNs also compute ν sgs and κ sgs ~2–4 times quicker than the dynamic Smagorinsky model resulting in a ~2–2.5 times acceleration of the entire simulation. This study demonstrates the feasibility of deep learning in emulating the subgrid-scale (SGS) phenomenon in geophysical flows accurately in a cost-effective manner. In a broader perspective, deep learning-based surrogate models can present a promising alternative to the traditional parameterizations of the subgrid-scale processes in climate models.

58 GEOSCIENCES↗