Search NASA⌕ Search

SEARCH · Search NASA

Results for “denotational”

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

Insight into ideal shear strength of Ni-based dilute alloys using first-principles calculations and correlational analysis

Here the present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni11X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the calculated Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress–strain curves.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A better understanding of the mechanics of borehole breakout utilizing a finite strain gradient-enhanced micropolar continuum model

Borehole breakout denotes the failure in rock mass subjected to drilling, caused by stress concentrations exceeding the material strength. Depending on the material properties, the preexisting in situ stress state, and the borehole dimensions, different types of borehole breakout, such as spiral-shaped breakout or v-shaped breakout, are distinguished in the literature. In the present work, we address the influence of the material properties on the predicted borehole breakout mode in a comprehensive finite element study. To this end, herein we employ a gradient-enhanced micropolar damage-plasticity model based on the Mohr–Coulomb strength criterion, formulated in the finite strain regime, which is calibrated based on experimental results from plane strain compression tests on Red Wildmoor sandstone. In the numerical study, the influence of the in situ stress state, the material friction angle, plastic dilation, post peak residual strength, and the inherent material length scale parameters are investigated. Thereby, we demonstrate that depending on the material parameters, substantially different failure modes, characterized by strongly localized shear bands or diffuse failure zones, are predicted. It is shown that in particular the brittleness of the material in the post peak regime has a major influence on the predicted breakout type. Moreover, a statistical validation of the results is obtained by considering different random field distributions of the initial material strength.

58 GEOSCIENCES↗
DOE OSTI2023

The protist Aurantiochytrium has universal subtelomeric rDNAs and is a host for mirusviruses

Viruses are the most abundant biological entities in the world's oceans, where they play important ecological and biogeochemical roles. Metagenomics is revealing new groups of eukaryotic viruses, although disconnected from known hosts. Among these are the recently described mirusviruses, which share some similarities with herpesviruses.50 years ago, "herpes-type" viral particles were found in a thraustochytrid member of the labyrinthulomycetes, a diverse group of abundant and ecologically important marine eukaryotes, but could not be further characterized by methods then available. Long-read sequencing has allowed us to connect the biology of mirusviruses and thraustochytrids. We sequenced the genome of the genetically tractable model thraustochytrid Aurantiochytrium limacinum ATCC MYA-1381 and found that its 26 linear chromosomes have an extraordinary configuration. Subtelomeric ribosomal DNAs (rDNAs) found at all chromosome ends are interspersed with long repeated sequence elements denoted as long repeated-telomere and rDNA spacers (LORE-TEARS). We identified two genomic elements that are related to mirusvirus genomes. The first is a ~300-kbp episome (circular element 1 [CE1]) present at a high copy number. Strikingly, the second, distinct, mirusvirus-like element is integrated between two sets of rDNAs and LORE-TEARS at the left end of chromosome 15 (LE-Chr15). Similar to metagenomically derived mirusviruses, these putative A. limacinum mirusviruses have a virion module related to that of herpesviruses along with an informational module related to nucleocytoplasmic large DNA viruses (NCLDVs). CE1 and LE-Chr15 bear striking similarities to episomal and endogenous latent forms of herpesviruses, respectively, and open new avenues of research into marine virus-host interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Characterizing Long COVID: Deep Phenotype of a Complex Condition

Background: Numerous publications describe the clinical manifestations of post-acute sequelae of SARS-CoV-2 (PASC or "long COVID"), but they are difficult to integrate because of heterogeneous methods and the lack of a standard for denoting the many phenotypic manifestations. Patient-led studies are of particular importance for understanding the natural history of COVID-19, but integration is hampered because they often use different terms to describe the same symptom or condition. This significant disparity in patient versus clinical characterization motivated the proposed ontological approach to specifying manifestations, which will improve capture and integration of future long COVID studies. Methods: The Human Phenotype Ontology (HPO) is a widely used standard for exchange and analysis of phenotypic abnormalities in human disease but has not yet been applied to the analysis of COVID-19. Funding: We identified 303 articles published before April 29, 2021, curated 59 relevant manuscripts that described clinical manifestations in 81 cohorts three weeks or more following acute COVID-19, and mapped 287 unique clinical findings to HPO terms. We present layperson synonyms and definitions that can be used to link patient self-report questionnaires to standard medical terminology. Long COVID clinical manifestations are not assessed consistently across studies, and most manifestations have been reported with a wide range of synonyms by different authors. Across at least 10 cohorts, authors reported 31 unique clinical features corresponding to HPO terms; the most commonly reported feature was Fatigue (median 45.1%) and the least commonly reported was Nausea (median 3.9%), but the reported percentages varied widely between studies. Interpretation: Translating long COVID manifestations into computable HPO terms will improve analysis, data capture, and classification of long COVID patients. If researchers, clinicians, and patients share a common language, then studies can be compared/pooled more effectively. Furthermore, mapping lay terminology to HPO will help patients assist clinicians and researchers in creating phenotypic characterizations that are computationally accessible, thereby improving the stratification, diagnosis, and treatment of long COVID.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Computationally efficient models for aqueous organic redox flow batteries

The rising usage of intermittent energy has garnered the need for large scale energy storage systems. Redox flow batteries (RFB) based energy storage system shows promising potential. Numerical simulations and machine learning approaches have been widely used to study RFB performance. The development of autonomous material discovery framework and digital twin of energy storage system usually needs to query cell performance through fast response models. In this study, two computationally efficient models are introduced: a physics-based analytical flow battery model (EZBattery), and a machine learning operator model (Deep Operator Network, denoted by DeepONet). Both models can provide cell performance near instantly, and prediction accuracy was systematically examined on an application of evaluating the performances of a 780 cm 2 aqueous organic redox flow battery (AORFB), using potential anolyte candidates in dihydroxyphenazine (DHP)-based family of organic materials. A validated computationally expansive 3-dimensional multi-physics finite element model by COMSOL was used as the ground truth and provided the training data set for the DeepONet. 1280 samples were generated with 10 properties to mimic the different possible anolyte candidates, and the cell performances were evaluated under 10 different combined operating conditions. The accuracy comparisons for the two computationally efficient models show that both models can provide comparable accuracy in predicting cell charging/discharging voltage curves. DeepONet can provide slightly higher overall accuracy than EZBattery with faster calculation speed, but highly relies on the training dataset. EZBattery does not need a training dataset and can provide interpretable physics-based explanations of the results, while being more flexible to adjust to adapt any different cell designs, flow battery architectures, and electrolyte materials.

Analytical model↗
DOE OSTI2021

Search for meteoritic GEMS II: Comparison of inclusions in amorphous silicates from the Paris chondrite and from anhydrous chondritic interplanetary dust particles

Amorphous silicates containing abundant nano-inclusions have been reported in the Paris CM chondrite (Leroux et al., 2015). They have chemical and morphological similarities to glass with embedded metal and sulfides (GEMS) found in interplanetary dust particles (IDPs) and micrometeorites believed to originate from comets. We used scanning transmission electron microscopy (STEM) with energy-dispersive X-ray spectroscopy (EDS) and nanodiffraction to study the chemistry and mineralogy of these inclusions in order to understand the origin of the GEMS-like material in Paris and its possible relationships to other materials found in primitive chondritic materials including IDP GEMS. EDS and diffraction analyses indicate compositional and mineralogical differences between the nanophase inclusions in cometary GEMS and Paris GEMS-like material. Metal inclusions are notably absent within Paris amorphous silicate. Ni-rich sulfides, including pentlandite, are common in even the least altered matrix material of Paris, while they are absent in GEMS-bearing IDPs and Ultracarbonaceous Antarctic Micrometeorites (UCAMMs). From examination of the inclusions, we cannot yet confirm or refute the possibility that GEMS-like material in Paris is related to cometary GEMS. The distinct compositions and mineralogy of the Paris material may be due to aqueous alteration of cometary GEMS precursors, but they may also denote an independent origin for meteoritic GEMS-like assemblages.

58 GEOSCIENCES↗
DOE OSTI2023

Fast approximation to supershell partition functions: Explicit forms of the coefficients

In a previous work, we derived a formula for supershell partition functions, which are the cornerstone of the Super-Transition-Array approach to radiative-opacity calculations. The new expression takes the form of a functional of the distribution of energies within the supershell and allows for fast and accurate computations, truncating the number of terms in the expansion. The latter involves coefficients (denoted Γ $k$ ) for which we obtained a recursion relation. In the present short paper, as a complement of the previous one, we give an explicit formula for the coefficients Γ $k$ . Another recursion relation is also provided, as well as an alternative expression involving Bell polynomials. In conclusion, the connections with cycle indexes of permutation groups and with elementary symmetric polynomials are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Model Predictive Control of Discrete-Continuous Energy Systems via Generalized Disjunctive Programming

Generalized Disjunctive Programming (GDP) provides an alternative framework to model optimization problems with both discrete and continuous variables. The key idea behind GDP involves the use of logical disjunctions to represent discrete decisions in the continuous space, and logical propositions to denote algebraic constraints in the discrete space. Compared to traditional mixed-integer programming (MIP), the inherent logic structure in GDP yields tighter relaxations that are exploited by global branch and bound algorithms to improve solution quality. In this paper, we present a general GDP model for optimal control of hybrid systems that exhibit both discrete and continuous dynamics. Specifically, we use GDP to formulate a model predictive control (MPC) model for piecewise-affine systems with implicit switching logic. As an example, the GDP-based MPC approach is used as a supervisory control to improve energy efficiency in residential buildings with binary on/off, relay-based thermostats. A simulation study is used to demonstrate the validity of the proposed approach, and the improved solution quality compared to existing MIPbased control approaches.

Bhattacharya, Arnab↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2023

Accuracy and precision of ESI-Orbitrap-IRMS observations of hours to tens of hours via reservoir injection

Here, orbitrap isotope ratio mass spectrometry (Orbitrap-IRMS) has recently been applied to high-precision, natural-abundance isotope ratio measurements of a diverse range of compounds, including amino acids, oxyanions, fatty acids, and metals. These measurements can characterize many isotope ratios simultaneously at high (≈1.0‰) precision. In a successful experiment, observed precision will track the shot-noise limit and be limited by experimental time. Some isotope ratios, for example those involving 17 O in organic compounds or multiply-substituted (‘clumped’) isotopologues, require experimental times of hours to tens of hours to achieve desired precision, while current sample introduction techniques focus on observations on the order of seconds to tens of minutes. In this study, we characterize Orbitrap-IRMS performance for three long duration measurements (individual acquisitions ≥1 h and as long as 24 h) using an automated reservoir injection system coupled to a Q Exactive HF Orbitrap with an electrospray ionization (ESI) source. First, we characterize long-term intra-measurement stability through a 24-h long measurement of acetone. We report the following isotope ratios and precisions (as acquisition errors, errors on the observed ratio within this measurement (σ AE )): 13 C/ 12 C (σ AE = 0.07‰), 17 O/ 16 O (σ AE = 1.1‰), 18 O/ 16 O (σ AE = 0.3‰), and 13 C 13 C/ 12 C (σ AE = 0.65‰). The σ AE of each tracks the shot noise limit throughout and is limited by the challenging conditions (high resolution and low numbers of ions per scan) required for 17 O/ 16 O measurement in the presence of 13 C via Orbitrap. Second, we characterize inter-measurement stability via a sequence of seven 75-min analyses of perchlorate. We observe the following ratios and acquisition errors: 37 Cl/ 35 Cl (σ AE = 0.09‰); 17 O/ 16 O (σ AE = 1.6‰); 18 O/ 16 O (σ AE = 0.7‰), 37 Cl 17 O/ 35 Cl 16 O (σ AE = 2.7‰), and 37 Cl 18 O/ 35 Cl 16 O (σ AE = 1.2‰). However, we find that inter-measurement drift between acquisitions limits our accuracy and precision for standardized measurements (i.e., error on reported δ values) to ≈1‰ for the 37 Cl/ 35 Cl measurement. Hence, the benefits of low σ AE may not be fully realized. Third, we demonstrate accuracy via sample/standard comparisons of a methionine sample with 13 C enrichment of ≈20‰ relative to a known standard. Using a sequence of seven 60-min analyses, we recover the following isotope ratios and standardized precisions (i.e., error on reported δ values, denoted propagated acquisition errors, σ PAE ): 33 S/ 32 S (σ PAE = 1.0‰), 34 S/ 32 S (σ PAE = 0.7‰), 15 N/ 14 N (σ PAE = 2.1‰), 2 H/ 1 H (σ PAE = 3.2‰), 13 C/ 12 C (σ PAE = 0.4‰), 18 O/ 16 O (σ PAE = 1.6‰), & 13 C 13 C/ 12 C (σ PAE = 2.8‰) with confirmation of accurate results for the known 13 C/ 12 C and 13 C 13 C/ 12 C enrichments. Together, our results demonstrate the viability of Orbitrap-IRMS for long duration measurements of diverse sample types via an automated reservoir injection system. Inter-measurement stability remains a challenge; we expect our methods to be most applicable to extended measurements of hard-to-observe properties, such as 17 O in organics and clumped isotopologues.

58 GEOSCIENCES↗
DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗
DOE OSTI2022

Investigating the impact of micropore volume of aminosilica functionalized SBA-15 on catalytic activity for amine-catalyzed reactions

Mesoporous materials such as SBA-15 are used in many applications since the materials are robust and have tunable properties. Whereas these materials are thought to be ideal, the materials are complex with the common synthesis methods producing materials that contain micropores. This work will investigate the impact of micropores on catalytic activity of aminosilica materials. For comparison, materials are formed using a standard method that produces materials with micropores (REG-SBA-15) and using a modified method that produces materials with negligible micropore volume (NMP-SBA-15). After grafting the calcined support with an aminosilane (either (N-methyl amino propyl) trimethoxy silane (denoted as 2°Am) or (N,N-diethyl-3-amino propyl) trimethoxy silane (3°Am)), catalytic testing experiments reveal that NMP-SBA-15 functionalized materials are more active than REG-SBA-15 functionalized materials for a range of reactions. The reactions include (a) glucose isomerization where amines interacting with surface silanols limit catalytic activity and (b) aldol and Knoevenagel chemistry where amine-silanol interactions increase the reaction rate. For aldol chemistry, the observed difference in rate is further investigated using site quantification experiments. Interestingly, the results reveal that three types of sites exist for each material, including sites that are (1) highly active, (2) intermediate active, and (3) inactive. Additionally, NMR experiments indicate that amines immobilized on REG-SBA-15 supports are less mobile than amines on NMP-SBA-15. For REG-SBA-15, the materials are found to have aminosilanes that have limited mobility and the largest fraction of inactive sites, which suggests that the inactive catalytic sites are in the micropores. Overall, the design of highly active catalytic materials can be achieved using materials that have limited micropore volumes.

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

Amorphous 1-D nanowires of calcium phosphate/pyrophosphate: A demonstration of oriented self-growth of amorphous minerals

Amorphous inorganic solids are traditionally isotropic, thus, it is believed that they only grow in a non-preferential way without the assistance of regulators, leading to the morphologies of nanospheres or irregular aggregates of nanoparticles. However, in the presence of (ortho)phosphate (Pi) and pyrophosphate ions (PPi) which have synergistic roles in biomineralization, the highly elongated amorphous nanowires (denoted ACPPNs) form in a regulator-free aqueous solution (without templates, additives, organics, etc). Based on thorough characterization and tracking of the formation process (e.g., Cryo-TEM, spherical aberration correction high resolution TEM, solid state NMR, high energy resolution monochromated STEM-EELS), the microstructure and its preferential growth behavior are elucidated. In ACPPNs, amorphous calcium orthophosphate and amorphous calcium pyrophosphate are distributed at separated but close sites. The ACPPNs grow via either the preferential attachment of ~2 nm nanoclusters in a 1-dimension way, or the transformation of bigger nanoparticles, indicating an inherent driving force-governed process. Here, we propose that the anisotropy of ACPPNs microstructure, which is corroborated experimentally, causes their oriented growth. This study proves that, unlike the conventional view, amorphous minerals can form via oriented growth without external regulation, demonstrating a novel insight into the structures and growth behaviors of amorphous minerals.

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

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization↗
DOE OSTI2021

Dual-atom active sites embedded in two-dimensional C 2 N for efficient CO 2 electroreduction: A computational study

Double-atom catalysts (DACs) have emerged as an enhanced platform of single-atom catalyst for promoting electrocatalytic CO 2 reduction reaction (CO 2 RR). Herein, we present a density-functional theory study on CO 2 RR performance of seven C 2 N-supported homo- and heteronuclear DACs, denoted as M 2 @C 2 N. Our results demonstrate that there exists substantial synergistic effect of dual-metal-atom N 2 M 2 N 2 active site and C 2 N matrix on O&#x2550 C &#x2550O bond activation. The dual-atom M 2 sites are able to drive CO 2 RR beyond C 1 products with low limiting potential (U L ). Specifically, C 2 H 4 formation is preferred on FeM@C 2 N (M = Fe, Co, Ni, Cu) versus CH 4 formation on CuM@C 2 N (M = Co, Ni, Cu). Furthermore, *CO+*CO co-binding strength can serve as a descriptor for CO 2 RR activity for making C 2 products such that the moderate binding results in the lowest U L . Remarkably, C-affinity matters most to C—C bond coupling and C 2 H 4 formation while both C- and O-affinity control CH 4 formation. Furthermore, our results provide theoretical insight into rational design of DACs for efficient CO 2 RR.

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

Oxygen vacancy-rich amorphous FeNi hydroxide nanoclusters as an efficient electrocatalyst for water oxidation

Here, a one-pot strategy is presented to directly synthesize amorphous Fe x Ni y hydroxide nanoclusters (denoted as ANC-Fe x Ni y , <2 nm) with oxygen vacancies induced by ionic liquids. The ANC-Fe x Ni y catalyst presents abundant catalytic sites and high intrinsic conductivity. As such, the optimized ANC-Fe 1 Ni 2 exhibits high activity in oxygen evolution reaction (OER) with a Tafel slope of 39 mV dec –1 and an overpotential of 266 mV at 10 mA cm –2 . Notably, the optimized ANC-Fe 1 Ni 2 shows an extraordinarily large mass activity of 3028 A g FeNi –1 at the overpotential of 300 mV, which is ~24-fold of commercial RuO 2 catalyst. The superior activity of these Fe x Ni y hydroxide nanoclusters is ascribed to (i) the amorphous and distorted structure with abundant oxygen vacancies, and (ii) enhanced active site density by downsizing the ANC-Fe x Ni y clusters. This strategy provides a novel route for enhancing OER electrocatalytic performance and highly encouraging for the future application of amorphous metal hydroxides in catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗