Search NASA⌕ Search

SEARCH · Search NASA

Results for “Partial Optimization”

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 145 records · Page 8

Operation at Reduced Atmospheric Pressure and Concept of Reliability Redundancy for Optimized Design of Insulation Systems

Electrified transportation is calling for insulation design criteria that is adequate to provide elevated levels of power density, power dynamics and reliability. Increasing voltage levels are expected to cause accelerated intrinsic and extrinsic aging effects which will not be easily predictable at the design stage due to a lack of suitable modeling. Designing reliable insulation systems would require finding solutions able to control accelerated aging due to an unpredictable increase of intrinsic stresses and the onset of extrinsic stresses as partial discharges. This paper proposes the concept of reliability redundancy for the insulation design of aerospace electrical asset components, which is also validated at lower-than-standard atmospheric pressure. The principle is that extrinsic-aging-free design might be achieved upon determining the aging stress or abnormal service stresses distribution and being sure that aging will not generate conditions that can incept extrinsic aging (partial discharges) during operation life. However, such information is never, in practice, fully available to insulation system designers. Hence, especially in critical applications such as electrified aircraft, aerospace, and combat ships a further level of reliability should be added to a partial-discharge-free design, which can consist of the use of corona-resistant materials and/or of life models able to consider the accelerated aging effect of partial discharges (or any other type of extrinsic-accelerated aging factor). Innovative life modeling considering both extrinsic and intrinsic aging stresses, insulating material testing to estimate model parameters, and a metric for quantifying the extent of corona (or partial discharge) resistance can lead to establishing feasibility and limit conditions for optimized or fully reliability-redundant design. It is shown in the paper that if an extrinsic-aging-free design is not feasible, and it is therefore replaced by a redundant design, a further level of reliability redundancy can be provided by effective condition monitoring plans.

Montanari, Gian Carlo↗

Dyn$\mathrm{AMO}$: Multi-agent reinforcement learning for dynamic anticipatory mesh optimization with applications to hyperbolic conservation laws

Here we introduce DynAMO, a reinforcement learning paradigm for Dynamic Anticipatory Mesh Optimization. Adaptive mesh refinement is an effective tool for optimizing computational cost and solution accuracy in numerical methods for partial differential equations. However, traditional adaptive mesh refinement approaches for time-dependent problems typically rely only on instantaneous error indicators to guide adaptivity. As a result, standard strategies often require frequent remeshing to maintain accuracy. In the DynAMO approach, multi-agent reinforcement learning is used to discover new local refinement policies that can anticipate and respond to future solution states by producing meshes that deliver more accurate solutions for longer time intervals. By applying DynAMO to discontinuous Galerkin methods for the linear advection and compressible Euler equations in two dimensions, we demonstrate that this new mesh refinement paradigm can outperform conventional threshold-based strategies while also generalizing to different mesh sizes, remeshing and simulation times, and initial conditions.

97 MATHEMATICS AND COMPUTING↗

Crystallization and Phase Transformations of Aluminum (Oxy)hydroxide Polymorphs in Caustic Aqueous Solution

Gibbsite, bayerite, and boehmite are important aluminum (oxy)hydroxide minerals in nature and have been widely deployed in various industrial applications. They are also major components in caustic nuclear wastes stored at various U.S. locations. Knowledge of their crystallization and phase transformation processes contributes to understanding their occurrence and could help optimize waste treatment processes. While it has been reported that partial conversion of bayerite and gibbsite to boehmite occurs in basic solutions at elevated temperatures, systematic studies of factors affecting the phase transformation as well as the underlying reaction mechanisms are non-existent, particularly in highly alkaline solutions. We explored the effects of sodium hydroxide concentrations (0.1~3 M), reaction temperature (60~100 ?) and aluminum concentrations (0.1~1 M) on the crystallization and transformation of these aluminum (oxy)hydroxides. Detailed structural and morphological characterization by X-ray diffraction (XRD), scanning electron microscopy (SEM), and nuclear magnetic resonance (NMR) spectrometry revealed that these processes depend largely on the reaction temperature and the Al/OH- ratio. When 1 = Al/OH- = 2.5, the reactions favor formation of high crystallinity precipitates, whereas at Al/OH- ratio ? 2.5 precipitation ceases unless the Al concentration is higher than 1 M. We identified pseudoboehmite, bayerite and gibbsite as intermediate phases to bayerite, gibbsite and boehmite, respectively, all of which transform via dissolution-reprecipitation. Gibbsite transforms to boehmite in both acidic and weak caustic environments at temperatures above 80 oC. However, a ‘bar-shaped’ gibbsite morphology dominates in highly caustic environments (3 M NaOH). The findings enable a robust basis for selection of various solid phases by tuning the reaction conditions.

aluminum (oxy)hydroxide polymorphs, crystallizatio↗

Rational solvent molecule tuning for high-performance lithium metal battery electrolytes

Electrolyte engineering improved cycling of Li metal batteries and anode-free cells at low current densities; however, high-rate capability and tuning of ionic conduction in electrolytes are desirable yet less-studied. In this study, we design and synthesize a family of fluorinated-1,2-diethoxyethanes as electrolyte solvents. The position and amount of F atoms functionalized on 1,2-diethoxyethane were found to greatly affect electrolyte performance. Partially fluorinated, locally polar –CHF 2 is identified as the optimal group rather than fully fluorinated –CF 3 in common designs. Paired with 1.2 M lithium bis(fluorosulfonyl)imide, these developed single-salt-single-solvent electrolytes simultaneously enable high conductivity, low and stable overpotential, >99.5% Li||Cu half-cell efficiency (up to 99.9%, ±0.1% fluctuation) and fast activation (Li efficiency >99.3% within two cycles). Combined with high-voltage stability, these electrolytes achieve roughly 270 cycles in 50-μm-thin Li||high-loading-NMC811 full batteries and >140 cycles in fast-cycling Cu||microparticle-LiFePO 4 industrial pouch cells under realistic testing conditions. The correlation of Li + –solvent coordination, solvation environments and battery performance is investigated to understand structure–property relationships.

25 ENERGY STORAGE↗

EMSL Community Science Campaign Meeting: Critical Minerals and Materials - Rhizo Critical Campaign Breakout Session Report Summary

The “Critical Minerals Biogeochemistry in the Rhizosphere – Ultramafic Soils (Rhizo Critical)” campaign breakout (BO) session was organized to identify major knowledge gaps and fundamental research needs in rhizosphere microbiology and geochemistry that, if addressed, could transform our ability to recover critical minerals from ultramafic soil systems. We sought to identify significant challenges that must be surmounted in the pursuit of deeper science knowledge. Our ultimate goal is to understand this landscape well enough to identify and prioritize opportunities for EMSL to make the greatest impact with Environmental Transformations and Interactions (ETI) science area research campaigns focused on the biogeochemical processes controlling the behavior of critical minerals and materials in the rhizosphere. The increasing demand for critical materials and minerals (CMM) in the U.S. has heightened interest in low-grade ores with much attention on ultramafic soils, which contain valuable metals such as nickel (Ni), chromium (Cr), manganese, cobalt (Co), and copper (Lee et al., 2025; DOE CMM Report, 2023) used in advanced battery, magnet, wiring and wind turbines, and stainless steel technologies. Metal hyperaccumulating plants grown in ultramafic soils can extract economically valuable concentrations of CMMs through the process of phytomining. This technology has evolved from phytoremediation, which involves using plants to cleanse contaminated environments by removing, detoxifying, or stabilizing pollutants like metals and organic compounds. Hyperaccumulator plants are capable of storing metals in their living tissues at concentrations hundreds to thousands of times higher than those found in 'normal' plants. For instance, while the average concentration of Ni in the dry matter of plants growing in typical soils is usually less than 5 µg g?¹, Ni hyperaccumulation is defined by concentrations exceeding 1,000 µg g?¹ (Corzo Remigio et al., 2020; Reeves et al., 2018). Phytomining research has primarily focused on Ni (Rylott and van der Ent, 2025), for which the U.S. has very limited conventional mines in operation. Most soils typically contain Ni concentrations ranging from 7 to 50 mg kg-1, whereas serpentine soils exhibit significantly higher levels, with Ni content often ranging between 700 and 8,000 mg kg-1 (Sobczyk et al., 2017). While more than 500 plant species in over 50 different families have been identified as Ni hyperaccumulators (Kidd et al., 2018), Ni phytomining (and phytominng in general) remains largely untested because most studies are short-term, small-scale, and conducted under simplified or artificially enriched conditions, so they fail to capture the low metal concentrations, environmental variability, and management constraints that would be needed for a field-scale demonstration. Few hyperaccumulator species have been validated as true “metal crops,” and their biomass production, stress tolerance, and rooting characteristics are usually too poor to yield economically meaningful metal outputs. Critically, the basic mechanisms of metal uptake, transport, and sequestration, especially as shaped by belowground processes such as root exudation, rhizosphere chemistry, and root–microbe interactions that control metal mobility and bioavailability (Montreemuk et al., 2023; Kidd et al., 2018; Durand et al., 2023; Alford et al., 2010), are still only partially understood, and downstream metal recovery from biomass is rarely optimized. Because these limitations stem from gaps in fundamental knowledge rather than from a failure of the concept itself (Rylott and van der Ent, 2025; van der Ent et al., 2015), there is a strong need for basic science that dissects plant metal homeostasis, rhizosphere and microbial processes, and their integration with soil chemistry and process engineering to design more robust, scalable phytomining systems.

Ahkami, Amirhossein↗

Efficient Loads Surrogates for Waked Turbines in an Array

Accurately and efficiently predicting wind turbine structural loading is a crucial step in wind farm design. Without considering structural loading, wind farm optimization could negatively impact turbine fatigue and ultimate loads, especially for waked and partially waked turbines, which could result in higher maintenance costs and reduced turbine lifetime. However, predicting turbine loads throughout an array is a costly step, as these quantities require time-accurate results across long time histories, which is often intractable for large array optimization. Therefore, surrogate models that link array spacing to load outputs are often used, but the surrogates are then unique to the inflow conditions and array configurations in the training library. This work develops surrogate models for many wind turbine load outputs based solely on rotor plane velocity measurements, with no required input about array configuration or freestream inflow parameters. Surrogate models were constructed for many turbine quantities of interest (QoI), considering mean, standard deviation, ultimate, and fatigue loads. In general, most QoI statistics were accurately captured, as measured by predicted vs. actual correlation coefficient, confirming the suitability of the approach. Temporal mean values of the QoI required only temporal mean measurements of the rotor plane inflow velocity. However, accurate prediction of temporal standard deviation, ultimate, and fatigue values of QoI also required temporal standard deviations of the rotor plane velocity field. Poor surrogate performance was observed when too many correlated inputs were used, such as multiple velocity components. If the fewest inflow parameters are used to construct the surrogates, the average correlation coefficient value for all output QoI statistics is 0.89. Surrogates for standard deviations and damage equivalent loads (DELs) of turbine QoIs generally had lower accuracy and tower-base and shaft load channels posed the most difficult to capture accurately. The results suggest that these surrogates could be easily paired with analytic wake models, which are frequently used for pre-construction wind farm array optimization, to account for turbine loading in addition to power production. By including the optimal inflow conditions, the surrogate accuracy can improve to an average correlation coefficient value for all output QoI statistics of 0.92. This work has established the ability to build accurate surrogates for mean, standard deviation, ultimate load, and DEL turbine QoI values based on the rotor plane inflow velocity, and identified which inflow conditions lead to greater surrogate accuracy.

17 WIND ENERGY↗

Genetic variation in nitrogen‐use efficiency and its associated traits in dryland winter wheat ( Triticum aestivum L.) cultivars released from the 1940s to the 2010s in Shaanxi Province, China

Abstract BACKGROUND Improving the nitrogen‐use efficiency (NUE) of wheat can help mitigate the problems of poor soil fertility under dryland conditions. We conducted field experiments using three nitrogen (N) fertilization levels (0, 120, and 180 kg ha −1 ) applied to eight dryland wheat cultivars to assess NUE and its associated traits. RESULTS The grain yield significantly increased with the improvement in variety, mainly as a result of a substantial increase in 1000‐grain weight and harvest index. Modern wheat varieties have stabilized at an optimal plant height and exhibited improved performance in terms of NUE, partial N productivity, N harvest index, and grain protein content compared to older varieties. The NUE of wheat gradually increased with variety replacement. The net photosynthesis rate of the flag leaves in the filling stage improved with the year of cultivar release; Increasing soil–plant analysis development (SPAD) values of flag leaves in the flowering and filling stages were observed over time, with the flag leaves of modern varieties showing a high chlorophyll content in the filling stage. Additionally, the principal component analysis showed that the SPAD value, grain number per unit area, transpiration rate, leaf area, and grain protein content positively contributed to the clustering of the N180 and modern cultivars (from the 2000s to 2010s). CONCLUSION Overall, high levels of N application did not significantly improve the NUE of wheat. However, modern wheat varieties can optimize N distribution, increase flag leaf photosynthetic capacity, and improve photosynthesis ability, thus enhancing NUE to achieve high yields under a suitable level of N supply. © 2022 Society of Chemical Industry.

Lian, Huida↗

Empirically Optimized One-Electron Pseudopotential for the Hydrated Electron: A Proof-of-Concept Study

Mixed quantum-classical molecular dynamics simulations have been important tools for studying the hydrated electron. They generally use a one-electron pseudopotential to describe the interactions of an electron with the water molecules. Furthermore, this approximation shows both the strength and weakness of the approach. On the one hand, it enables extensive statistical sampling and large system sizes that are not possible with more accurate ab initio molecular dynamics methods. On the other hand, there has (justifiably) been much debate about the ability of pseudopotentials to accurately and quantitatively describe the hydrated electron properties. These pseudopotentials have largely been derived by fitting them to ab initio calculations of an electron interacting with a single water molecule. In this paper, we present a proof-of-concept demonstration of an alternative approach in which the pseudopotential parameters are determined by optimizing them to reproduce key experimental properties. Specifically, we develop a new pseudopotential, using the existing TBOpt model as a starting point, which correctly describes the hydrated electron vertical detachment energy and radius of gyration. In addition to these properties, this empirically optimized model displays a significantly modified solvation structure, which improves, for example, the prediction of the partial molar volume.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Random coordinate descent: A simple alternative for optimizing parameterized quantum circuits

Variational quantum algorithms rely on the optimization of parameterized quantum circuits in noisy settings. The commonly used back-propagation procedure in classical machine learning is not directly applicable in this setting due to the collapse of quantum states after measurements. Thus, gradient estimations constitute a significant overhead in a gradient-based optimization of such quantum circuits. This paper introduces a random coordinate descent algorithm as a practical and easy-to-implement alternative to the full gradient descent algorithm. This algorithm only requires one partial derivative at each iteration. Motivated by the behavior of measurement noise in the practical optimization of parameterized quantum circuits, this paper presents an optimization problem setting that is amenable to analysis. Under this setting, the random coordinate descent algorithm exhibits the same level of stochastic stability as the full gradient approach, making it as resilient to noise. The complexity of the random coordinate descent method is generally no worse than that of the gradient descent and can be much better for various quantum optimization problems with anisotropic Lipschitz constants. Theoretical analysis and extensive numerical experiments validate our findings. Published by the American Physical Society 2024

Ding, Zhiyan (ORCID:000000018863403X)↗

Where did the tumor start? An inverse solver with sparse localization for tumor growth models

In this work, we present a numerical scheme for solving an inverse problem for parameter estimation in tumor growth models for glioblastomas, a form of aggressive primary brain tumor. The growth model is a reaction–diffusion partial differential equation (PDE) for the tumor concentration. We use a PDE-constrained optimization formulation for the inverse problem. The unknown parameters are the reaction coefficient (proliferation), the diffusion coefficient (infiltration), and the initial condition field for the tumor PDE. Segmentation of magnetic resonance imaging (MRI) scans drive the inverse problem where segmented tumor regions serve as partial observations of the tumor concentration. Like most cases in clinical practice, we use data from a single time snapshot. Moreover, the precise time relative to the initiation of the tumor is unknown, which poses an additional difficulty for inversion. We perform a frozen-coefficient spectral analysis and show that the inverse problem is severely ill-posed. We introduce a biophysically motivated regularization on the structure and magnitude of the tumor initial condition. In particular, we assume that the tumor starts at a few locations (enforced with a sparsity constraint on the initial condition of the tumor) and that the initial condition magnitude in the maximum norm is equal to one. We solve the resulting optimization problem using an inexact quasi-Newton method combined with a compressive sampling algorithm for the sparsity constraint. Our implementation uses PETSc and AccFFT libraries. We conduct numerical experiments on synthetic and clinical images to highlight the improved performance of our solver over a previously existing solver that uses standard two-norm regularization for the calibration parameters. The existing solver is unable to localize the initial condition. Our new solver can localize the initial condition and recover infiltration and proliferation. In clinical datasets (for which the ground truth is unknown), our solver results in qualitatively different solutions compared to the two-norm regularized solver.

97 MATHEMATICS AND COMPUTING↗

Optimization of passive superconductors for shaping stellarator magnetic fields

Here, we consider the problem of optimizing a set of passive superconducting coils (PSCs) with currents induced by a background magnetic field rather than power supplies. In the nuclear fusion literature, such coils have been proposed to partially produce the 3D magnetic fields for stellarators and provide passive stabilization. We perform the first optimizations of PSC arrays with respect to the orientation, shape, and location of each coil, jointly minimized with the background fields. We conclude by generating passive coil array solutions for four stellarators.

coil optimization↗

Learning Optimal Multigrid Smoothers via Neural Networks

Multigrid methods are one of the most efficient techniques for solving large sparse linear systems arising from partial differential equations (PDEs) and graph Laplacians from machine learning applications. One of the key components of multigrid is smoothing, which aims at reducing high-frequency errors on each grid level. However, finding optimal smoothing algorithms is problem-dependent and can impose challenges for many problems. In this paper, we propose an efficient adaptive framework for learning optimized smoothers from operator stencils in the form of convolutional neural networks (CNNs). Here, the CNNs are trained on small-scale problems from a given type of PDEs based on a supervised loss function derived from multigrid convergence theories and can be applied to large-scale problems of the same class of PDEs. Numerical results on anisotropic rotated Laplacian problems and variable coefficient diffusion problems demonstrate improved convergence rates and solution time compared with classical hand-crafted relaxation methods.

97 MATHEMATICS AND COMPUTING↗

LANL Activities on Mechanistic Approach to Analyzing and Improving Unconventional Hydrocarbon Production

Hydrocarbon production from shale reservoirs is inherently inefficient and challenging since these are low permeability plays. In addition, there is a limited understanding of the fundamentals and the controlling mechanisms, further complicating how to optimize these plays. Herein, we summarize our experimental and computational efforts fully and partially supported by the fundamental shale portfolio to reveal unconventional shale fundamentals and devise development strategies to enhance extraction efficiency with a minimal environmental footprint. Integrating these fundamentals with machine learning, we outline a pathway to improve the predictive power of our models, which enhances the forecast quality of production, thereby improving the economics of operations in unconventional reservoirs. For instance, we have developed science informed workflows and platforms for optimizing pressure-drawdown at a site, which allow operators to make reservoir-management decisions that optimize recovery in consideration of future production. Recently, our work relies on the hybridization of physics-based prediction and machine learning, whereby accurate synthetic data (combined with available site data) can enable the application of machine learning methods for rapid forecasting and optimization. Consequently, the workflow and platform are readily extendable to operations at other sites, plays, and basins.

04 OIL SHALES AND TAR SANDS↗

Optimization of energy requirements for CO 2 post-combustion capture process through advanced thermal integration

The energy optimization modeling work described here was performed to determine efficiency improvements that could be achieved for existing coal-fired power plants to retrofit a partial CO 2 capture from the post-combustion flue gas for carbon sequestration through thermal integration. The work presented includes optimization of the mono-ethanol amine (MEA)-based post-combustion CO 2 capture to reduce energy requirements that could be achieved at existing power plants by thermal integration of the steam turbine cycle, boiler, CO 2 compression train and post-combustion CO 2 capture process to offset efficiency and capacity losses that would be incurred by retrofit or implementation of post-combustion CO 2 capture. Furthermore, partial CO 2 capture, involving treatment of less than 100% of the flue gas leaving the plant and modular design of the CO 2 scrubbing system, was also investigated. Thermal integration of the steam turbine cycle with boiler and CO 2 compression train improved cycle and plant performance and offset, in part, the negative effects of post-combustion CO 2 capture. The best-analyzed integration options improved gross power output by 5% and net unit efficiency by 1.57%, relative to the conventional MEA process. Operating with 40% CO 2 capture increased gross power output by 11.6–14% (depending on the MEA thermal integration option), relative to the conventional MEA integration and 90% CO 2 capture. The improvement in net unit performance is larger compared to the improvement in turbine cycle performance because of the CO 2 compression work, which is also reduced by partial CO 2 capture.

42 ENGINEERING↗

Disordered Rocksalts as High‐Energy and Earth‐Abundant Li‐Ion Cathodes

To address the growing demand for energy and support the shift toward transportation electrification and intermittent renewable energy, there is an urgent need for low‐cost, energy‐dense electrical storage. Research on Li‐ion electrode materials has predominantly focused on ordered materials with well‐defined lithium diffusion channels, limiting cathode design to resource‐constrained Ni‐ and Co‐based oxides and lower‐energy polyanion compounds. Recently, disordered rocksalts with lithium excess (DRX) have demonstrated high capacity and energy density when lithium excess and/or local ordering allow statistical percolation of lithium sites through the structure. This cation disorder can be induced by high temperature synthesis or mechanochemical synthesis methods for a broad range of compositions. DRX oxides and oxyfluorides containing Earth‐abundant transition metals have been prepared using various synthesis routes, including solid‐state, molten‐salt, and sol‐gel reactions. This review outlines DRX design principles and explains the effect of synthesis conditions on cation disorder and short‐range cation ordering (SRO), which determines the cycling stability and rate capability. In addition, strategies to enhance Li transport and capacity retention with Mn‐rich DRX possessing partial spinel‐like ordering are discussed. Finally, the review considers the optimization of carbon and electrolyte in DRX materials and addresses key challenges and opportunities for commercializing DRX cathodes.

Li-ion batteries↗

Toward real-time optimization through model reduction and model discrepancy sensitivities

Optimization problems arise in a range of scenarios, from optimal control to model parameter estimation. In many applications, such as the development of digital twins, it is essential to solve these optimization problems within wall-clock-time limitations. However, this is often unattainable for complex systems, such as those modeled by nonlinear partial differential equations. One strategy for mitigating this issue is to construct a reduced-order model (ROM) that enables more rapid optimization. In particular, the use of nonintrusive ROMs—those that do not require access to the full-order model at evaluation time—is popular because they facilitate the computation of optimization solutions within the wall-clock time requirements. However, the optimization solution will be unreliable if the iterates move outside the ROM training data. This article proposes the use of hyper-differential sensitivity analysis with respect to model discrepancy (HDSA-MD) as a computationally efficient tool to augment ROM-constrained optimization and improve its reliability. The proposed approach consists of two phases: (i) an offline phase where several full-order model evaluations are computed to train the ROM, and (ii) an online phase where a ROM-constrained optimization problem is solved, a limited number of full-order model evaluations are computed, and HDSA-MD is used to enhance the optimization solution. Numerical results are demonstrated for two examples, atmospheric contaminant control and wildfire ignition location estimation, in which a ROM is trained offline using inaccurate atmospheric data. In conclusion, the HDSA-MD update yields a significant improvement in the ROM-constrained optimization solution using only one full-order model evaluation online with corrected atmospheric data.

PDE-constrained optimization↗

The Baghdad Atlas: A relational database of inelastic neutron-scattering (n,n ' γ) data

A relational database has been developed based on the original (n,n'γ) work carried out by A. M. Demidov et al., at the Nuclear Research Institute in Baghdad, Iraq (Demidov et al., 1978) for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: γ-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted (n,n'γ) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial γ-ray cross section for the production of the 846.8-keV 21+→0gs+γ-ray transition in 56Fe, determined to be $\langle$σγ$\rangle$=143(29) mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas (n,n'γ) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Discovery of Probabilistic Dirichlet-to-Neumann Maps on Graphs

Dirichlet-to-Neumann maps enable the coupling of multiphysics simulations across computational subdomains by ensuring continuity of state variables and fluxes at artificial interfaces. We present a novel method for learning Dirichlet-to-Neumann maps on graphs using Gaussian processes, specifically for problems where the data obey a conservation law arising from an underlying partial differential equation. Our approach combines discrete exterior calculus and nonlinear optimal recovery to infer relationships between vertex and edge values. This framework yields data-driven predictions with uncertainty quantification across the entire graph, even when observations are limited to a subset of vertices and edges. By minimizing the reproducing kernel Hilbert space norm while penalizing kernel complexity through maximum likelihood estimation, our method ensures that the resulting surrogate strictly enforces conservation laws without overfitting. We demonstrate our method on two representative applications: subsurface flow in fracture networks and arterial blood flow. Finally, the results demonstrate that the method maintains high accuracy and well-calibrated uncertainty estimates even under severe data scarcity, highlighting its potential for scientific applications where limited data and reliable uncertainty quantification are critical.

Dirichlet-to-Neumann map↗