Search NASA⌕ Search

SEARCH · Search NASA

Results for “differential expansion”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Transverse momentum dependent PDFs at N3LO

We compute the quark and gluon transverse momentum dependent parton distribution functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Drell-Yan and gluon fusion Higgs production cross sections about their collinear limit. This method allows us to employ cutting edge multiloop techniques for the computation of cross sections to extract these universal building blocks of the collinear limit of QCD. The corresponding perturbative matching kernels for all channels are expressed in terms of simple harmonic polylogarithms up to weight five. As a byproduct, we confirm a previous computation of the soft function for transverse momentum factorization at N 3 LO. Our results are the last missing ingredient to extend the q T subtraction methods to N 3 LO and to obtain resummed q T spectra at N 3 LL' accuracy both for gluon as well as for quark initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Multi-messenger dynamic imaging of laser-driven shocks in water using a plasma wakefield accelerator

Understanding dense matter hydrodynamics is critical for predicting plasma behavior in environments relevant to laser-driven inertial confinement fusion. Traditional diagnostic sources face limitations in brightness, spatiotemporal resolution, and in their ability to detect relevant electromagnetic fields. In this work, we present a dual-probe, multi-messenger laser wakefield accelerator platform combining ultrafast X-rays and relativistic electron beams at 1 Hz, to interrogate a free-flowing water target in vacuum, heated by an intense 200 ps laser pulse. This scheme enables high-repetition-rate tracking the evolution of the interaction using both particle types. Betatron X-rays reveal a cylindrically symmetric shock compression morphology assisted by low-density vapor, resembling foam-layer-assisted fusion targets. The synchronized electron beam detects time-evolving electromagnetic fields, uncovering charge separation and ion species differentiation during plasma expansion – phenomena not captured by photons or hydrodynamic simulations. We show that combining both probes provides complementary insights spanning kinetic to hydrodynamic regimes, highlighting the need for hybrid physics models to accurately predict fusion-relevant plasma behavior.

imaging techniques↗
DOE OSTI2024

A systems genetic analysis identifies putative mechanisms and candidate genes regulating vessel traits in poplar wood

Wood is the water conducting tissue of tree stems. Like most angiosperm trees, poplar wood contains water-conducting vessel elements whose functional properties affect water transport and growth rates, as well as susceptibility to embolism and hydraulic failure during water stress and drought. Here we used a unique hybrid poplar pedigree carrying genomically characterized chromosomal insertions and deletions to undertake a systems genomics analysis of vessel traits. We assayed gene expression in wood forming tissues from clonal replicates of genotypes covering dosage quantitative trait loci with insertions and deletions, genotypes with extreme vessel trait phenotypes, and control genotypes. A gene co-expression analysis was used to assign genes to modules, which were then used in integrative analyses to identify modules associated with traits, to identify putative molecular and cellular processes associated with each module, and finally to identify candidate genes using multiple criteria including dosage responsiveness. These analyses identified known processes associated with vessel traits including stress response, abscisic acid and cell wall biosynthesis, and in addition identified previously unexplored processes including cell cycle and protein ubiquitination. We discuss our findings relative to component processes contributing to vessel trait variation including signaling, cell cycle, cell expansion, and cell differentiation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Conditional Karhunen-Loève expansion for uncertainty quantification and active learning in partial differential equation models

We use a conditional Karhunen-Lo` eve (KL) model to quantify and reduce uncertainty in a stochastic partial differential equation (SPDE) problem with partially-known space-dependent coefficient, Y (x). We assume that a small number of Y (x) measurements are available and model Y (x) with a KL expansion. We achieve reduction in uncertainty by conditioning the KL expansion coefficients on measurements. We consider two approaches for conditioning the KL expansion: In Approach 1, we condition the KL model first and then truncate it. In Approach 2, we first truncate the KL expansion and then condition it. We employ the conditional KL expansion together with Monte Carlo and sparse grid collocation methods to compute the moments of the solution of the SPDE problem. Uncertainty of the problem is further reduced by adaptively selecting additional observation locations using two active learning methods. Method 1 minimizes the variance of the PDE coefficient, while Method 2 minimizes the variance of the solution of the PDE. We demonstrate that conditioning leads to dimension reduction of the KL representation of Y (x). For a linear diffusion SPDE with uncertain log-normal coefficient, we show that Approach 1 provides a more accurate approximation of the conditional log-normal coefficient and solution of the SPDE than Approach 2 for the same number of random dimensions in a conditional KL expansion. Furthermore, Approach 2 provides a good estimate for the number of terms of the truncated KL expansion of the conditional field of Approach 1. Finally, we demonstrate that active learning based on Method 2 is more efficient for uncertainty reduction in the SPDE’s states (i.e., it leads to a larger reduction of the variance) than active learning using Method 2.

Conditioned Karhunen-Lo` eve expanion, machine lea↗
DOE OSTI2024

Demonstration and Verification of Thermo-Mechanical Bowing in a Limited Free-Bow SFR Concept Using MOOSE

Core bowing due to thermal gradients and irradiation induced swelling and creep introduces significant reactivity feedback effects in liquid metal-cooled fast spectrum reactors. During startup, normal operations, and transient events, differential thermal and flux gradients cause expansion, creep, and swelling, which results in bowing in preferential directions depending on local material properties and load-pad and restraint-ring design. The bowing phenomenon produces negative reactivity during accident conditions provided the restraint system has been properly designed to optimally guide the deformation outwards in the active core region. Under the Department of Energy Nuclear Energy Advanced Modeling and Simulation program, a Multiphysics Object Oriented Simulation Environment (MOOSE)-based multiphysics approach to model core bowing is being developed. The present work expands on previous modeling of simpler International Atomic Energy Agency (IAEA) verification problems with these tools by modeling IAEA Verification Problem 4, which involves a symmetric sector of a reactor core with ducted assemblies undergoing differential thermal expansion due to thermal gradients bowing outward with duct-to-duct contact. This model is available on the National Reactor Innovation Center Virtual Test Bed repository. This example was verified against benchmark participant results, which includes bowing deformation evaluation and duct-to-duct interactions at load pads with mechanical contact.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A scalable and tunable thermoreversible polymer for 3D human pluripotent stem cell biomanufacturing

Human pluripotent stem cells (hPSCs) are an exciting and promising source to enable cell replacement therapies for a variety of unmet medical needs. Though hPSCs can be successfully derived into numerous physiologically relevant cell types, effective translation to the clinic is limited by challenges in scalable production of high-quality cells, cellular immaturity following the differentiation process, and the use of animal-derived components in culture. To address these limitations, we have developed a fully defined, reproducible, and tunable thermoreversible polymer for high-quality, scalable 3D cell production. Our reproducible synthesis method enables precise control of gelation temperature (24°C–32°C), hydrogel stiffness (100–4000 Pa), and the prevention of any unintended covalent crosslinking. After material optimization, we demonstrated hPSC expansion, pluripotency maintenance, and differentiation into numerous lineages within the hydrogel. Overall, this 3D thermoreversible hydrogel platform has broad applications in scalable, high-quality cell production to overcome the biomanufacturing burden of stem cell therapy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Engineering Privacy at the Edge: A Practical Guide to Differential Privacy in System Architectures

The rapid expansion of distributed and edge computing platforms—spanning autonomous vehicles, IoT sensors, and healthcare monitors—has heightened concerns about data privacy. Differential Privacy (DP) offers a rigorous mathematical framework to protect sensitive information while retaining analytical utility. This tutorial introduces the foundations of DP for both numerical and categorical datasets and extends the discussion to correlation-aware techniques tailored for structured and high-dimensional data. Hands-on demonstrations will begin with the PETINA (Privacy prEservaTIoN Algorithms) package for numerical data and continue with MIC-DP (Maximum Information Correlated Differential Privacy) for tabular data. Designed for researchers and practitioners in secure systems, embedded architectures, and AI accelerators, the tutorial emphasizes practical and scalable methods for integrating DP into real-world system designs.

Kotevska, Olivera [ORNL] (ORCID:0000000316772243)↗
DOE OSTI2024

Calculation of machine precision second order derivatives using dual-complex numbers

It is well known that both complex and dual numbers can be employed to obtain machine precision first-order derivatives; however, neither, on their own, can compute machine precision 2nd order derivatives. To address this limitation, it is demonstrated in this paper that combined dual-complex numbers can be used to compute machine precision 1st and 2nd order derivatives. The dual-complex approach is simpler than utilizing multicomplex or hyper-dual numbers as existing dual libraries can be used as is or easily augmented to accept complex numbers, and the complexity of developing, integrating, and deploying multicomplex or hyper-dual libraries is avoided. The efficacy of this approach is demonstrated for both univariate and multivariate functions. Finally, source code examples using the Python, Julia, and Mathematica languages are provided as supplemental material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

TANTE: Time-adaptive operator learning via neural Taylor expansion

Operator learning for time-dependent partial differential equations (PDEs) has seen rapid progress in recent years, enabling efficient approximation of complex spatiotemporal dynamics. However, most existing methods rely on fixed time step sizes during rollout, which limits their ability to adapt to varying temporal complexity and often leads to error accumulation. In this work, we propose the Time-Adaptive Transformer with Neural Taylor Expansion (TANTE), a novel operator-learning framework that produces continuous-time predictions with adaptive step sizes. TANTE predicts future states by performing a Taylor expansion at the current state, where neural networks learn both the higher-order temporal derivatives and the local radius of convergence. This allows the model to dynamically adjust its rollout based on the local behavior of the solution, thereby reducing cumulative error and improving computational efficiency. We demonstrate the effectiveness of TANTE across a wide range of PDE benchmarks, achieving superior accuracy and adaptability compared to fixed-step baselines, delivering accuracy gains of 60-80 % and speed-ups of 30-40 % at inference time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗
DOE OSTI2023

Mef2c regulates bone mass through Sost -dependent and -independent mechanisms

We report Mef2c is a transcription factor that mediates key cellular behaviors that promote endochondral ossification and bone formation. Previously, Mef2c has been shown to regulate Sost transcription via its osteocyte-specific enhancer, ECR5, and conditional deletions of Mef2 cfl/fl with either Col1-Cre or Dmp1-Cre produced generalized high bone mass (HBM) consistent with Van Buchem Disease phenotypes. However, Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre mice produced a significantly higher bone mass phenotype that Sost-/- alone suggesting that Mef2c modulates bone mass through additional mechanisms, independent of Sost. To identify new Mef2c transcriptional targets important in bone metabolism, we profiled gene expression by single-cell RNA sequencing in subpopulations of cells isolated from Mef2 cfl/fl ; Dmp1-Cre and Mef2 cfl/fl ; Bglap-Cre femurs, both strains exhibiting similar high bone mass phenotypes. However, we found Mef2cfl/fl; Bglap-Cre to also display a growth plate defect characterized by an expansion of several osteoprogenitor subpopulations. Differential gene expression analysis identified a total of 96 up- and 2434 down- regulated genes in Mef2 cfl/fl ; Bglap-Cre and 176 up- and 1041 down- regulated genes in Mef2 cfl/fl ; Dmp1-Cre bone cell subpopulations compared to wildtype mice. Mef2c deletion affected the transcriptomes across several cell types including mesenchymal progenitors (MP), osteoprogenitors (OSP), osteoblast (OB), and osteocyte (OCY) subpopulations. Several energy metabolism genes such as Uqcrb, Ndufv2, Ndufs3, Ndufa13, Ndufb9, Ndufb5, Cox6a1, Cox5a, Atp5o, Atp5g2, Atp5b, Atp5 were significantly down regulated in Mef2c-deficient OBs and OCYs, in both strains. Binding motif analysis of promoter regions of differentially expressed genes identified Mef2c binding in Bone Sialoprotein (BSP/Ibsp), a gene known to cause increased trabecular BV/TV in the femurs of Ibsp -/- mice. Immunohistochemical analysis confirmed the absence of Ibsp protein in OBs and OCYs. These findings suggests that the HBM in Sost -/- ; Mef2 cfl/fl ; Dmp1-Cre is caused by a multitude of transcriptional changes in genes that regulate bone formation, two of which are Sost and Ibsp.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Evaluating the benefits of kinetic Monte Carlo and microkinetic modeling for catalyst design studies in the presence of lateral interactions

Popular computational catalyst design strategies rely on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends. The main benefit of this approach is related to the inherent dimensionality reduction of the large material space to just a few catalyst descriptors. Conversely, it is well documented that a small set of descriptors is insufficient to capture the intricacies and complexities of a real catalytic system. The inclusion of coverage effects through lateral adsorbate-adsorbate interactions can narrow the gap between simplified descriptor predictions and real systems, but mean-field MKMs cannot properly account for local coverage effects. This shortcoming of the mean-field approximation can be rectified by switching to a lattice-based kinetic Monte Carlo (kMC) method using cluster expansion representation of adsorbate–adsorbate lateral interactions. Using the prototypical CO oxidation reaction as an example, we critically evaluate the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. Here, the mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, the MKM approach does provide a better entry point for computational catalyst design.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unifying principles for catalytic hydrotreating processes (Final Technical Report)

This project builds on the hypothesis that the hydrotreating processes for the removal of oxygen and sulfur are fundamentally similar at the atomic-scale and existing knowledge from the treatment of petroleum derived feedstock can be leveraged for the design of novel catalysts for the upgrade of bio-oil. We tested this hypothesis by comparing computed potential energy diagrams for hydrodesulfurization (HDS) of thiophene over MoS 2 with hydrodeoxygenation (HDO) of furan over MoO 3 and concluded that certain aspects, such as catalyst promotion with transition metals, are valid strategies for both reactions. On the other hand, we also noticed significant differences in the mechanism for hydrogen (H 2 ) activation, which requires sites with metallic character. While MoS 2 is known to have metallic edge states that can catalyze H 2 dissociation, this elementary step is prohibitively slow on defect-free oxides. Only in the presence of vacancies or by creating metal/oxide interfaces can efficient H 2 activation sites during HDO be formed. The need for bifunctional catalyst when it comes to efficient and selective HDO or dehydrogenation reactions was further corroborated in joint experimental and theoretical studies of the Guerbet reaction for the coupling of biomass derived oxygenates over PdCu alloys, nitrate reduction over In-promoted Pd nanoparticles, and ethylene dehydroaromatization over Ga-exchanged ZSM-5 zeolites. All of these catalytic systems have in common that catalytic sites with distinct functional requirements are needed to create a working catalyst. Detailed computational studies were carried out for HDO of m -cresol and phenol on Ru-modified TiO 2 surfaces, which allowed us to attribute catalytic activity to the metal/oxide interface. A surprising finding was that heterolytic cleavage of the H-H bond across the Ru/TiO 2 interface was critically important, despite lower barriers for homolytic H 2 activation on Ru metal. The explanation lies in the high barriers for hydrogen spillover from Ru to TiO 2 , which becomes unnecessary in the heterolytic activation pathway. Moreover, we also reported that proton transfer steps between metal and oxide sites are mediated by weakly adsorbed surface water. During attempts to develop and validate a kinetic Monte Carlo (kMC) model for HDO reactions at the Ru/TiO 2 interface, it became clear that lateral interactions are paramount to describe realistic surface chemistry and without these interactions, the reduction and hydroxylation behavior from our simulations was inconsistent with reported experiments. To assess the importance of lateral interactions in popular computational catalyst design strategies relying on the identification of reactivity descriptors, which can be used along with Brønsted–Evans–Polanyi (BEP) and scaling relations as input to a microkinetic model (MKM) to make predictions for activity or selectivity trends, we compared predicted trends with those obtained from descriptor-based kMC models. We critically evaluated the benefits of kMC over MKM in terms of trend predictions and computational cost when using only a small set of input parameters. After confirming that in the absence of lateral interactions the kMC and MKM approaches yield identical trends and mechanistic information, we observed substantial differences between the two kinetic models when lateral interactions were introduced. The mean-field implementation applies coverage corrections directly to the descriptors, causing an artificial overprediction of the activity of strongly binding metals. In contrast, the cluster expansion in kMC implementation can differentiate among the highly active metals but it is very sensitive to the set of included interaction parameters. Considering that computational screening relies on a minimal set of descriptors, for which MKM makes reasonable trend predictions at a ca. three orders of magnitude lower computational cost than kMC, we concluded that the MKM approach does provide an overall better entry point for computational catalyst design. Overall, this project has led to 11 peer-reviewed publications, and their scientific is impact is well illustrated by their combined 702 citations.

08 HYDROGEN↗
DOE OSTI2023

Turboexpander for Direct Cooling in Hydrogen Vehicle Fueling Infrastructure

Hydrogen fuel cell electric vehicles (FCEVs) have been identified as one of a few options for zero carbon emissions transportation. A major advantage of FCEVs is that they can fuel quickly and follow a familiar fueling behavior to hydrocarbon-fueled vehicles. Whether light duty or heavy duty, the goal for a hydrogen dispenser is to fuel a vehicle in the same amount of time as the fossil fuel equivalent. When hydrogen is dispensed into the vehicle storage system, however, the temperature rises due to the Joule-Thomson effect and the heat of compression. Typically, vehicles store the compressed hydrogen in composite overwrapped pressure vessels that have a polymer liner with an operational temperature limit of 85°C. This temperature limit can be exceeded during fast fueling if hydrogen is not precooled. Precooling allows for a dispenser to fuel a vehicle at a faster flow rate by preventing the storage tank on the vehicle from overheating. Fueling protocols and requirements are presented in SAE J2601 Fueling Protocols for Light Duty Gaseous Hydrogen Surface Vehicles [1]. A heavy-duty equivalent is under development with similar requirements for precooling. Currently, conventional precooling for light-duty vehicle refueling uses a heat exchanger and chiller to cool the hydrogen gas to -40°C before entering the vehicle. The precooling system represents a significant part of the station capital and operating costs, so if the cost of the precooling system can be reduced by improving its efficiency, the overall station capital and operating cost can be reduced. In this project, National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) researchers teamed up to investigate the turboexpander precooling application. A turboexpander is a device that places a turbine in a flow path where a pressure differential can be attained. This expansion device will extract work and lower the temperature of the fluid as the pressure reduces. While initial calculations based on established principles showed potential for a turboexpander to generate cooled gas, much work needs to be done to prove the concept. Turboexpanders typically work best under steady state conditions, while the dispenser is a very dynamic flow system. Dynamic turboexpander systems have been proven, such as a turbocharger on a gasoline vehicle. The inlet pressure at a dispenser is also much higher than any other known turboexpander system but should behave similarly to higher density fluids at lower pressures. Having both performed initial calculations, NREL and SNL researchers teamed up to investigate the turboexpander precooling application further. A project was soon built around the idea with SNL performing system modeling using previously proven capabilities and NREL performing hardware characterization with established station capabilities. Creare LLC was contracted as the turbomachinery expert to design and build the concept device. Part way through the project, however, contracting issues with the funding partner caused the project to terminate early before building and characterizing the concept device. While the project could not continue, many key findings were already learned. This paper is a summary of those findings.

08 HYDROGEN↗
DOE OSTI2024

A Portable Miniature Cryogenic Environment for In Situ Neutron Diffraction

Neutron diffraction instruments offer a platform for materials science and engineering studies at extended temperature ranges far from ambient. As one of the widely used neutron sample environment types, cryogenic furnaces are usually bulky and complex, and they may need hours of beamtime overhead for installation, configuration, cooling, and sample change, etc. To reduce the overhead time and expedite experiments at the state-of-the-art high-flux neutron source, we developed a low-cost, miniature, and easy-to-use cryogenic environment (77–473 K) for in situ neutron diffraction. A travel-size mug serves for the environment where the samples sit inside. Immediate cooling and an isothermal dwell at 77 K are realized on the sample by direct contact with liquid N 2 in the mug. The designed Al inserts serve as the holder of samples and heating elements, alleviate the thermal gradient, and clear neutron pathways. Both a single-sample continuous measurement and multi-sample high-throughput measurements are demonstrated in this environment. High-quality and refinable in situ neutron diffraction patterns are acquired on model materials. The results quantify the orthorhombic-to-cubic phase transformation process in LiMn 2 O 4 and differentiate the anisotropic lattice thermal expansions and bond length evolutions between rhombohedral perovskite oxides with composition variation.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Analytical solutions of the Arrhenius-Semenov problem for constant volume burn

Analytical solutions to the Semenov thermal ignition problem for constant volume burn governed by Arrhenius reaction kinetics are derived. Specifically, an approximate analytical solution technique for the Arrhenius-Semenov differential equation is derived for reaction orders n ϵ R> 0 and exact solutions are also constructed for reaction orders n ϵ N : n ≤ 3. The approximation technique relies on expansion of the respective nondominant terms in the differential equation at the lower and upper bounds of the reaction progress variable in order to create a pair of integrable series. The two integrated series are then connected to create a single continuous analytical solution. Excellent agreement is observed between the analytical approximation and solutions obtained numerically. The presented approximation constitutes a simple and robust strategy for solving the Arrhenius-Semenov problem analytically.

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

Sparse-grid discontinuous Galerkin methods for the Vlasov–Poisson–Lenard–Bernstein model

Sparse-grid methods have recently gained interest in reducing the computational cost of solving high-dimensional kinetic equations. In this paper, we construct adaptive and hybrid sparse-grid methods for the Vlasov–Poisson–Lenard–Bernstein (VPLB) model. This model has applications to plasma physics and is simulated in two reduced geometries: a 0x3v space homogeneous geometry and a 1x3v slab geometry. Here we use the discontinuous Galerkin (DG) method as a base discretization due to its high-order accuracy and ability to preserve important structural properties of partial differential equations. We utilize a multiwavelet basis expansion to determine the sparse-grid basis and the adaptive mesh criteria. We analyze the proposed sparse-grid methods on a suite of three test problems by computing the savings afforded by sparse-grids in comparison to standard solutions of the DG method. The results are obtained using the adaptive sparse-grid discretization library ASGarD.

97 MATHEMATICS AND COMPUTING↗