Search NASA⌕ Search

SEARCH · Search NASA

Results for “adaptive robust 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics: A DOE-NIH Joint Workshop Executive Summary

The DOE-NIH Joint Workshop on Computational Modeling to Advance Novel Medical Isotopes for Radiotheranostics, held on September 27, 2024, brought together experts from government, academia, and industry to address critical challenges in radionuclide production and clinical translation. Here, the workshop emphasized interdisciplinary collaboration, particularly between the Department of Energy (DOE) and the National Institutes of Health (NIH), to strengthen the domestic isotope supply, streamline regulatory pathways, and further integrate computational tools into radiopharmaceutical therapy (RPT). Key discussions explored the role of AI-driven modeling, machine learning, and digital twin technologies in optimizing dosimetry, dynamically personalizing treatments, and reducing time to clinical adoption. Advances in predictive computational modeling were highlighted as essential for improving radionuclide yield, purity, and synthesis efficiency. Regulatory considerations and equitable access were central themes, with participants advocating for harmonized global standards, adaptive trial designs, and expanded infrastructure for clinical implementation. DOE computational and production infrastructure was emphasized. Future priorities identified include increased investment in radionuclide production infrastructure, expanded workforce development in radiopharmaceutical sciences and computational modeling, and the creation of robust public-private partnerships. The workshop concluded that continued strategic collaboration and sustained resources will be vital for advancing next-generation radiotheranostics, ensuring safe and effective therapies accessible to all patients.

digital twins↗
DOE OSTI2024

Online Monitoring of Catalytic Processes by Fiber-Enhanced Raman Spectroscopy

An innovative solution for real-time monitoring of reactions within confined spaces, optimized for Raman spectroscopy applications, is presented. This approach involves the utilization of a hollow-core waveguide configured as a compact flow cell, serving both as a conduit for Raman excitation and scattering and seamlessly integrating into the effluent stream of a cracking catalytic reactor. The analytical technique, encompassing device and optical design, ensures robustness, compactness, and cost-effectiveness for implementation into process facilities. Notably, the modularity of the approach empowers customization for diverse gas monitoring needs, as it readily adapts to the specific requirements of various sensing scenarios. As a proof of concept, the efficacy of a spectroscopic approach is shown by monitoring two catalytic processes: CO 2 methanation (CO 2 + 4H 2 → CH 4 + 2H 2 O) and ammonia cracking (2NH 3 → N 2 + 3H 2 ). Leveraging chemometric data processing techniques, spectral signatures of the individual components involved in these reactions are effectively disentangled and the results are compared to mass spectrometry data. This robust methodology underscores the versatility and reliability of this monitoring system in complex chemical environments.

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

Intraspecific Variation for Leaf Physiological and Root Morphological Adaptation to Drought Stress in Alfalfa ( Medicago sativa L.)

Drought stress reduces crop biomass yield and the profitability of rainfed agricultural systems. Evaluation of populations or accessions adapted to diverse geographical and agro-climatic environments sheds light on beneficial plant responses to enhance and optimize yield in resource-limited environments. This study used the morphological and physiological characteristics of leaves and roots from two different alfalfa subspecies during progressive drought stress imposed on controlled and field conditions. Two different soils (Experiments 1 and 2) imposed water stress at different stress intensities and crop stages in the controlled environment. Algorithm-based image analysis of leaves and root systems revealed key morphological and physiological traits associated with biomass yield under stress. The Medicago sativa subspecies (ssp.) sativa population, PI478573, had smaller leaves and maintained higher chlorophyll content (CC), leaf water potential, and osmotic potential under water stress. In contrast, M. sativa ssp. varia, PI502521, had larger leaves, a robust root system, and more biomass yield. In the field study, an unmanned aerial vehicle survey revealed PI502521 to have a higher normalized difference vegetation index (vegetation cover and plant health characteristics) throughout the cropping season, whereas PI478573 values were low during the hot summer and yielded low biomass in both irrigated and rainfed treatments. RhizoVision Explorer image analysis of excavated roots revealed a smaller diameter and a narrow root angle as target traits to increase alfalfa biomass yield irrespective of water availability. Root architectural traits such as network area, solidity, volume, surface area, and maximum radius exhibited significant variation at the genotype level only under limited water availability. Different drought-adaptive strategies identified across subspecies populations will benefit the plant under varying levels of water limitation and facilitate the development of alfalfa cultivars suitable across a broad range of growing conditions. The alleles from both subspecies will enable the development of drought-tolerant alfalfa with enhanced productivity under limited water availability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Hybrid Energy Management with Real-Time Control of a High-Power EV Charging Site

Decarbonization of transportation systems is driving higher capacity energy storage and faster charging power requirements in electric vehicles (EVs). Given the potential advantages - such as increased efficiency, reduced inverter capacity, and less total cable mass - there is a demand in the industry for more DC distribution for high-power charging (HPC) hubs. However, the cost-effective, adaptive, and robust operation of the DC-coupled HPC hub necessitates a robust site energy management system (SEMS). Validating SEMS operation using a digital twin of an HPC hub in a real-time simulator (RTS) platform is crucial before field deployment. In this study, we propose a hybrid energy management site controller designed to achieve high-level, long-term operational objectives while managing low-level power sharing control between hub assets. We develop a centralized model predictive controller (MPC) to optimize hub operating points and use these points to update the droop parameters of the site energy storage system (ESS). This approach ensures the hub follows an optimal operating point while maintaining the flexibility to respond to load surges. We tested and verified our proposed approach both offline and on a Controller Hardware-in-the-loop (C-HIL) simulation platform integrated with a SEMS framework, demonstrating real-time site operation and validating a cost-effective and robust site controller.

ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Analytical Identification Method of Generalized Short‐Circuit Ratio Using Phasor Measurement Units

This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.

Han, Zelei [Hohai University, Nanjing (China)] (OR↗
DOE OSTI2023

Coupling multi-fidelity xRAGE with machine learning for graded inner shell design optimization in double shell capsules

Bayesian optimization has shown promise for the design optimization of inertial confinement fusion targets. Specifically, in Vazirani et al. [Phys. Plasmas 28 , 122709 (2021)], optimal designs for double shell capsules with graded inner shells were identified using one-dimensional xRAGE simulation yield calculations. While the machine learning models were able to accurately learn and predict one-dimensional simulation target performance, using simulations with higher fidelity would improve design optimization and better match with the expected experimental performance. However, higher fidelity physics modeling, i.e., two-dimensional xRAGE simulations, requires significantly larger computational time/cost, usually at least an order of magnitude, in comparison with one-dimensional simulations. This study presents a multi-fidelity Bayesian optimization, in which the machine learning model leverages low-fidelity (one-dimensional xRAGE) and high-fidelity (two-dimensional xRAGE) simulations to more accurately predict “pre-shot” target performance with respect to the expected experimental performance. By building a multi-fidelity Bayesian optimization framework coupled with xRAGE, the low-fidelity and high-fidelity simulations are able to inform one another, such that we have: (1) improved physics modeling in comparison with using low-fidelity simulations alone, (2) reduced computational time/cost in comparison with using high-fidelity simulations alone, and (3) more confidence in the expected performance of optimized targets during real-world experiments. In the future, we plan to use this robust multi-fidelity Bayesian optimization methodology to expedite the design of graded inner shells further and eventually full capsules as a part of the current double shell campaign at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Voltage positioning using co-optimization of controllable grid assets in radial networks

With increasing penetration of solar PV, some distribution feeders are experiencing highly variable net-load flows and even reverse flows. To optimize distribution systems under such conditions, the scheduling of mechanical devices, such as OLTCs and capacitor banks, needs to take into account forecasted solar PV and actual grid conditions. However, these legacy switching assets are operated on a daily or hourly timescale, due to the wear and tear associated with mechanical switching, which makes them unsuitable for real-time control. Therefore, there is a natural timescale- separation between these slower mechanical assets and the responsive nature of inverter-based resources. In this paper, we present a network admissible convex formulation for holistically scheduling controllable grid assets to position voltage optimally against solar PV. An optimal hourly schedule is presented that utilizes mechanical resources to position the predicted voltages close to nominal values, while minimizing the use of inverter-based resources (i.e., DERs), making them available for control at a faster time-scale (after the uncertainty reveals itself). A convex, inner approximation of the OPF problem is adapted to a mixed-integer linear program that minimizes voltage deviations from nominal (i.e., maximizes voltage margins). Here, the resulting OPF solution respects all the network constraints and is, hence, robust against modeling simplifications. Simulation based analysis on IEEE distribution feeders validates the approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Learning to Solve Large-Scale Security-Constrained Unit Commitment Problems

Security-constrained unit commitment (SCUC) is a fundamental problem in power systems and electricity markets. In practical settings, SCUC is repeatedly solved via mixed-integer linear programming (MIP), sometimes multiple times per day, with only minor changes in input data. In this work, we propose a number of machine learning techniques to effectively extract information from previously solved instances in order to significantly improve the computational performance of MIP solvers when solving similar instances in the future. Based on statistical data, we predict redundant constraints in the formulation, good initial feasible solutions, and affine subspaces where the optimal solution is likely to lie, leading to a significant reduction in problem size. Computational results on a diverse set of realistic and large-scale instances show that using the proposed techniques, SCUC can be solved on average 4.3 times faster with optimality guarantees and 10.2 times faster without optimality guarantees, with no observed reduction in solution quality. Out-of-distribution experiments provide evidence that the method is somewhat robust against data-set shift. Summary of Contribution. The paper describes a novel computational method, based on a combination of mixed-integer linear programming (MILP) and machine learning (ML), to solve a challenging and fundamental optimization problem in the energy sector. The method advances the state-of-the-art, not only for this particular problem, but also, more generally, in solving discrete optimization problems via ML. We expect that the techniques presented can be readily used by practitioners in the energy sector and adapted, by researchers in other fields, to other challenging operations research problems that are solved routinely.

Machine Learning↗
DOE OSTI2026

Application of Peracetic Acid in Poultry Processing: Effects of Treatment Dynamics and Emerging Risk of Resistance Development in Salmonella spp.

Nontyphoidal Salmonella is a leading cause of foodborne illness, with poultry representing a major source. Peracetic acid (PAA), a widely adopted antimicrobial in poultry processing, offers advantages over traditional disinfectants but has sparked interest in its combined use with other antimicrobials and potential resistance. This review evaluates the efficacy of PAA in mitigating Salmonella in combination with other food-grade antimicrobials, explores possible synergism, efficacy under varying treatment parameters, resistance development against PAA, and its role in resistance evolution. While PAA demonstrates broad-spectrum efficacy, its performance varies with environmental parameters; higher temperatures generally enhance antimicrobial action but also accelerate PAA degradation. Organic matter diminishes PAA efficacy by reactive quenching. Variability in concentration and contact time further influences outcomes. Despite its oxidative mode of action and presumed low risk for resistance, emerging studies indicate that Salmonella can develop adaptive tolerance and potential cross/coresistance following repeated or sublethal exposure to PAA. These adaptations may involve genetic upregulation of oxidative stress response pathways, efflux systems, and modifications in cell membrane integrity, raising concerns about the long-term sustainability of PAA use. Additionally, combinatorial treatments (e.g., PAA with UV-C, enzymes, or organic acids) show promise in enhancing efficacy while mitigating resistance risks. Despite recognition of PAA's safety and effectiveness, knowledge gaps remain regarding standardized resistance definitions, serotype-specific tolerance, and optimal intervention strategies in commercial settings. Therefore, there is a need for standardized testing protocols, robust studies on potential resistance, and further exploration of synergistic PAA applications to ensure sustained poultry product safety and public health protection.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

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

Autonomous Cyber Defense Against Dynamic Multi-strategy Infrastructural DDoS Attacks

Dynamic Infrastructural Distributed Denial of Service (I-DDoS) attacks constantly change attack vectors to congest core backhaul links and disrupt critical network availability while evading end-system defenses. To effectively counter these highly dynamic attacks, defense mechanisms need to exhibit adaptive decision strategies for real-time mitigation. This paper presents a novel Autonomous DDoS Defense framework that employs model-based reinforcement agents. The framework continuously learns attack strategies, predicts attack actions, and dynamically determines the optimal composition of defense tactics such as filtering, limiting, and rerouting for flow diversion. Our contributions include extending the underlying formulation of the Markov Decision Process (MDP) to address simultaneous DDoS attack and defense behavior, and accounting for environmental uncertainties. We also propose a fine-grained action mitigation approach robust to classification inaccuracies in Intrusion Detection Systems (IDS). Additionally, our reinforcement learning model demonstrates resilience against evasion and deceptive attacks. Evaluation experiments using real-world and simulated DDoS traces demonstrate that our autonomous defense framework ensures the delivery of approximately 96 - 98% of benign traffic despite the diverse range of attack strategies.

Dutta, Ashutosh↗
DOE OSTI2025

Biochemical properties of glycerol kinase from the hypersaline-adapted archaeon Haloferax volcanii

ABSTRACT Extremophilic microorganisms are promising candidates for industrial and analytical biocatalysis.Haloferax volcanii, a halophilic archaeon that prefers glycerol over glucose, channels this substrate into central metabolism through glycerol kinase (GK). Here, we report the biochemical properties ofH. volcaniiGK and its potential for biotechnological applications. An N-terminal His-tagged GK was functionalin vivoand yielded 3 mg/L culture—4.5 times more enzyme than a C-terminal StrepII-tagged version. Size exclusion chromatography revealed a glycerol-induced oligomeric shift from homodimer to a dimer-dominant state with detectable tetramer. The purified enzyme showed robust activity across broad pH and salinity ranges, with optimal activity at 100 mM NaCl and 50°C–60°C. It retained catalytic activity in 5%–10% dimethyl sulfoxide (DMSO) and crude glycerol containing methanol. His-GK was freeze-thaw stable and thermotolerant in 2 M NaCl buffers. In the absence of ligands, the enzyme’s melting temperature (T m ) was 80°C. Glycerol increased the T m to 85°C, and combinations with MgCl₂ (84°C) or ATP (88°C) provided further stabilization. The highest T m (89°C) occurred with all three ligands, suggesting a cumulative stabilizing effect. Kinetic analyses revealed positive cooperativity for glycerol, ATP, and Mg² + ; Mn² + and Co² + also supported the activity.H. volcaniiGK is the first known GK to exhibit positive cooperativity with glycerol and ATP. Its high stability and substrate flexibility support its use in biodiesel waste valorization,in vitrobiocatalysis, and biosensor development—applications demanding robust, specific, and stable enzymes. IMPORTANCE This study reveals thatH. volcaniiGK exhibits positive cooperativity for glycerol, ATP, and Mg² + , a kinetic feature not previously reported for glycerol kinases. This behavior enables steep, switch-like responses to small substrate changes, offering unique advantages for biosensor design. Importantly,H. volcaniiGK also maintains high activity under extreme salinity, temperature, broad pH, and solvent conditions that typically limit enzyme use in industrial and environmental applications. These traits make this GK an ideal candidate for enzyme-based biosensors, which often suffer from poor tolerance to pH, solvent, and thermal stress. Its robustness supports its use in cross-linked enzyme crystals, an immobilization method that enhances enzyme stability and reusability under harsh conditions. Moreover, GKs are already employed in Mg² + detection kits; however,H. volcaniiGK’s ability to tolerate and respond to diverse divalent cations (e.g., Co² + , Mn² + ) broadens their potential for pollutant detection and environmental monitoring. These features collectively positionH. volcaniiGK as a valuable biocatalyst for biosensing,in vitrodiagnostics, and biotechnological applications requiring both precision and durability.

Biotechnology & Applied Microbiology↗
DOE OSTI2025

Dynamical Sketching for Enhanced Communication Efficiency in Federated Learning

Federated learning (FL) has revolutionized distributed machine learning by enabling collaborative model training without sharing local data. However, communication efficiency and privacy guarantees remain significant challenges. This paper introduces a dynamic sketching mechanism in FL, optimizing the trade-off between communication efficiency and model accuracy. By dynamically selecting the sketch matrix size, our approach adapts to the evolving characteristics of the data and the model, ensuring optimal performance across diverse scenarios. We leverage Bayesian optimization to systematically tune the sketch parameters, achieving an effective balance between resource efficiency and model performance. Experimental results on the MNIST dataset using a convolutional neural network (CNN) architecture validate the proposed method's efficiency and scalability. Our dynamic sketching approach significantly outperforms fixed-size sketching techniques, achieving higher compression ratios (up to 62x) and providing better privacy guarantees while maintaining high model accuracy. These findings highlight the robustness and versatility of our approach and make it a valuable solution for privacy-preserving, communication-efficient federated learning.

Afrose, Sharmin [ORNL]↗
DOE OSTI2020

Single Episode Policy Transfer in Reinforcement Learning

Transfer and adaptation to new unknown environmental dynamics is a key challenge for reinforcement learning (RL). An even greater challenge is performing near-optimally in a single attempt at test time, possibly without access to dense rewards, which is not addressed by current methods that require multiple experience rollouts for adaptation. To achieve single episode transfer in a family of environments with related dynamics, we propose a general algorithm that optimizes a probe and an inference model to rapidly estimate underlying latent variables of test dynamics, which are then immediately used as input to a universal control policy. This modular approach enables integration of state-of-the-art algorithms for variational inference or RL. Moreover, our approach does not require access to rewards at test time, allowing it to perform in settings where existing adaptive approaches cannot. In diverse experimental domains with a single episode test constraint, we report our method significantly outperforms existing adaptive approaches and shows favorable performance against baselines for robust transfer.

42 ENGINEERING↗
DOE OSTI2026

Surface-Functionalized, Two-Dimensional Polymer Electrochromic Layers as Ultrafast, Multi-State Infrared Optical Gates

Electrochromic devices have widespread application potential, but the currently available switching speeds limit broad real-world implementation of this technology. Here, we report surface-engineered two-dimensional polymers with ionophilic pores that offer unprecedented switching speeds in solid-state, two-terminal, electrochromic devices. In particular, we demonstrate that a crystalline donor–acceptor 2DP functionalized with ethylene glycol oligomers exhibits multistate infrared absorption that is 4× faster (t c = 320 ms) with 3× coloration efficiency (491 cm 2 C –1 ) compared to an alkyl functionalized 2DP constructed from the same chromophores. The functionalized nanoporous surfaces enable rapid switching in these materials under either oxidative or reductive conditions, allowing us to access a range of robust, stable optical responses in a single electrochromic layer. These attributes led us to leverage surface-functionalized 2DPs as multistate infrared logic gates. Collectively, this work demonstrates that surface engineering of nanoporous crystalline lattices is a promising approach to co-optimize the electronic and ionic conductivities required to achieve rapidly switchable electrochromic layers. Beyond speed and efficiency, the demonstration of multistate infrared characteristics shows that electrochromic frameworks are useful in integrated optoelectronic circuits. This positions surface-engineered 2DPs as improved electrochromic coatings and a new material platform for photonic information processing and adaptive devices.

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

FunDiff: diffusion models over function spaces for physics-informed generative modeling

Recent advances in generative modeling-particularly diffusion models and flow matching-have been widely used for synthesizing discrete data such as images and videos. However, adapting these models to physical applications remains challenging, as the quantities of interest are continuous functions governed by complex physical laws. To address this, we introduce FunDiff, an efficient and robust framework for generative modeling in function spaces. FunDiff combines a latent diffusion process with a function autoencoder architecture to handle input functions with varying discretizations, generates continuous functions that can be evaluated at arbitrary locations, and seamlessly incorporate physical priors. These priors are enforced through architectural constraints or physics-informed loss functions, ensuring that generated samples satisfy fundamental physical laws. We theoretically establish minimax optimality guarantees for density estimation in function spaces, demonstrating that diffusion-based estimators achieve optimal convergence rates under suitable regularity conditions. We further demonstrate the practical effectiveness of FunDiff across diverse applications in fluid dynamics and solid mechanics. Empirical results indicate that our method can generate physically consistent samples with high fidelity to the target distribution, and exhibit robustness to noisy and low-resolution data.

Wang, Sifan [Yale University, New Haven, CT (Unite↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗