SEARCH · Engineering Papers

Results for “parametric model dependence”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Linear gyrokinetic stability of a high β non-inductive spherical tokamak

Spherical tokamaks (STs) have been shown to possess properties desirable for a fusion power plant such as achieving high plasma β and having increased vertical stability. To understand the confinement properties that might be expected in the conceptual design for a high β ST fusion reactor, a 1 GW ST plasma equilibrium was analyzed using local linear gyrokinetics to determine the type of micro-instabilities that arise. Kinetic ballooning modes and micro-tearing modes are found to be the dominant instabilities. The parametric dependence of these linear modes was determined and, from the insights gained, the equilibrium was tuned to find a regime marginally stable to all micro-instabilities at θ 0 = 0.0. This work identifies the most important micro-instabilities expected to generate turbulent transport in high β STs. The impact of such modes must be faithfully captured in first-principles-based reduced models of anomalous transport that are needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Parametrization of Generalized Parton Distributions from 𝑡-Channel String Exchange in AdS Spaces

We introduce a string-based parametrization for nucleon quark and gluon generalized parton distributions (GPDs) that is valid for all skewness. Our approach leverages conformal moments, representing them as the sum of spin-𝑗 nucleon 𝐴-form factor and skewness-dependent spin-𝑗 nucleon 𝐷-form factor, derived from 𝑡-channel string exchange in AdS spaces consistent with Lorentz invariance and unitarity. This model-independent framework, satisfying the polynomiality condition due to Lorentz invariance, uses Mellin moments from empirical data to estimate these form factors. With just five Regge slope parameters, our method accurately produces various nucleon quark GPD types and symmetric nucleon gluon GPDs through pertinent Mellin-Barnes integrals. Our isovector nucleon quark GPD is in agreement with existing lattice data, promising to improve the empirical extraction and global analysis of nucleon GPDs in exclusive processes, by avoiding the deconvolution problem at any skewness, for the first time.

QCD phenomenology↗

DOE OSTI2025

Five Year Comparison of Mixing Height Determinations at the Savannah River Site

Air quality dispersion modeling is performed for the Savannah River Site (SRS) to demonstrate compliance with applicable regulations. The AMS/EPA Regulatory Model (AERMOD) modeling system is an EPA recommended model for air quality applications with a data preprocessor (AERMET) to incorporate meteorological data collected on site. AERMET parameterizes or calculates meteorological variables that are not directly measured onsite. One of the parameters estimated by AERMET is the atmospheric mixing height. While the mixing height is not currently a measurement input into AERMET, SRS has the capability to measure the local mixing height. The Savannah River National Laboratory (SRNL) operates a Vaisala CL31 Lidar Ceilometer which estimates mixing height from aerosol backscatter. This study compares the parameterized mixing height from AERMET to the ceilometer estimated mixing height for the current regulatory period at SRS incorporating data from 2015-2019. Results from this study showed the average daily minimum values (morning) from AERMET were an order of magnitude lower than the commonly used Holzworth (1972) method and the ceilometer estimated mixing heights. Additionally, on average, the ceilometer exhibited a daily maximum mixing height value that occurred 1-3 hours later than the AERMET estimated maximum. This difference is likely due to the nighttime atmospheric mixing height assumptions and calculations used by AERMET. The AERMET algorithm cuts off mixing height growth at sunset while the ceilometer data show ongoing evening convection typical of the southeastern United States. These results suggest that the AERMET parametrization scheme assumptions may not be representative of a forested landscape and evening convection which could account for more mixing overnight. The results obtained in this study are significant for air dispersion modeling applications for regulatory purposes and worker safety. Mixing height can impact model estimated pollutant concentrations. A greater mixing height will provide more volume for pollutant dispersion. This report documents efforts to quantify the dependence of mixing height inputs toward a conservative estimated pollutant concentration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2022

Self-organized patterns at the plasma–liquid anode interface in a helium glow discharge: temporal development and mechanisms

Self-organization at the plasma–liquid anode interface is a commonly observed phenomenon for atmospheric pressure glow plasmas, resulting in patterns with distinctive shapes such as circular ring, star-shaped, and gear-like structures, depending primarily on the discharge current and solution conductivity. Recent studies have shown that the electrode gap distance, solute used for liquid anode solution, and gas composition can also significantly impact pattern formation. Nonetheless, an overarching model or explanation of the key underlying mechanisms consistent with all experimentally observed trends is not yet reported. We propose a key underlying mechanism enabling pattern formation motivated by a detailed parametric study of pattern formation complemented by the temporal development of patterns and consistent with all observed trends. Pattern formation was observed to be on a time scale of 100 μs, similar to the time scales of gas heating and evaporation. It was found that a minimum water evaporation rate of (3.5 ± 0.5) × 10 –6 kg s –1 and reduced electric field in the positive column of 16.6 ± 0.4 Td is required for pattern formation in the investigated cases irrespective of solution conductivity and gas composition for NaCl solutions. Nonetheless, the presence of cations for which the corresponding metal atom has a low ionization energy was identified as a necessary condition for pattern formation. Furthermore, the reported results suggest that the presence of a small amount of metal atoms in the gas phase with low ionization energy enhances the overall ionization rate in the near anode region which triggers pattern formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

DESI DR2 results. II. Measurements of baryon acoustic oscillations and cosmological constraints

We present baryon acoustic oscillation (BAO) measurements from more than 14 million galaxies and quasars drawn from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2 (DR2), based on three years of operation. For cosmology inference, these galaxy measurements are combined with DESI Lyman-𝛼 forest BAO results presented in a companion paper (M. Abdul-Karim et al., companion paper, Phys. Rev. D 112, 083514 2025.). The DR2 BAO results are consistent with DESI DR1 and the Sloan Digital Sky Survey, and their distance-redshift relationship matches those from recent compilations of supernovae (SNe) over the same redshift range. The results are well described by a flat Λ cold dark matter (Λ⁢CDM) model, but the parameters preferred by BAO are in mild, 2.3⁢𝜎 tension with those determined from the cosmic microwave background (CMB), although the DESI results are consistent with the acoustic angular scale 𝜃 * that is well measured by Planck. This tension is alleviated by dark energy with a time-evolving equation of state parametrized by 𝑤0 and 𝑤𝑎, which provides a better fit to the data, with a favored solution in the quadrant with 𝑤 0 >−1 and 𝑤 𝑎 <0. This solution is preferred over Λ ⁢CDM at 3.1⁢𝜎 for the combination of DESI BAO and CMB data. When also including SNe, the preference for a dynamical dark energy model over Λ⁢ CDM ranges from 2.8 − 4.2⁢𝜎 depending on which SNe sample is used. We present evidence from other data combinations which also favor the same behavior at high significance. From the combination of DESI and CMB we derive 95% upper limits on the sum of neutrino masses, finding ∑𝑚 𝜈 < 0.064 eV assuming Λ ⁢CDM and ∑𝑚 𝜈 < 0.16 eV in the 𝑤0⁢𝑤𝑎 model. Unless there is an unknown systematic error associated with one or more datasets, it is clear that Λ⁢ CDM is being challenged by the combination of DESI BAO with other measurements and that dynamical dark energy offers a possible solution.

Baryon acoustic oscillations↗

DOE OSTI2025

Evolving Dark Sector and the Dark Dimension Scenario

String theory naturally leads to the expectation that dark energy is not stable, and may be evolving as captured by the Swampland de Sitter conjectures. Moreover, motivated by the distance conjecture, a unification of dark sector has been proposed, where the smallness of dark energy leads to one extra dimension of micron size with dark matter being the Kaluza--Klein graviton excitations in this extra dimension. We consider the natural possibility that the radius of the dark dimension varies as the dark energy decreases, leading to the variation of the dark matter mass. This correlates the decrease of the dark energy with the variation of the dark matter mass as they depend on the variations of a scalar field $ϕ$ controlling the radius of the extra dimension. A simple realization of this idea for small range of $ϕ$ is adequately captured by choosing a potential which is locally of the form $V=V_0\ {\rm exp}(-cϕ)$ and dark matter mass $m_{\rm DM}=m_0\ {\rm exp}(-c' ϕ)$ where the sign of $ϕ$ is chosen such that $c'\geq 0$ while we have two choices for the sign of $c$ depending on whether the dark dimension expands or shrinks when the dark energy dominates. We find excellent agreement with recent experimental data from DESI DR2 combined with SN measurements (from DES, Union3 or Pantheon+) and reproduces the same significance as CPL parametrization with the added benefit of providing a natural explanation for the apparent phantom behavior ($w<-1$) reported by DESI and DES based on a physical model. DESI and SN datasets independently favor non-zero values of $c'$ and $c$, respectively, both lying within the expected $\mathcal{O}(1)$ range suggested by the Swampland criteria. Moreover, our best fit value $c'\simeq 0.05 \pm 0.01$ is remarkably consistent with the experimental upper bound of $c'\lesssim 0.2$ demanded by the lack of detection of fifth force in the dark sector.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗

DOE OSTI2021

Experiment and non-local crystal plasticity finite element study of nanoindentation on Al-8Ce-10Mg alloy

Cerium and magnesium strengthened aluminum alloys, or Al-Ce-Mg, is a recently developed alloy family that exhibits good mechanical properties at elevated temperatures (~300 °C). To examine the single-crystal properties of Al-Ce-Mg alloys, nanoindentation experiments are conducted in this study. A crystal plasticity finite element model (CPFEM) with the evolution of geometrically necessary dislocations (GNDs) is applied to simulate the indentation in individual grains. A parametric study is carried out to investigate the sensitivity of each crystal plasticity model parameter to the indentation behavior. The highly sensitive parameters are calibrated by matching the indentation load-depth curves, while the rest parameters are obtained from bulk polycrystal uniaxial tension tests. Overall, satisfactory matching between experiment and simulation is obtained for each individual grain. The calculated hardness, as determined from the experiment, shows the dependence on indent depth, which is captured by the GND model. Furthermore, the effect of grain orientation and neighboring grains to nanoindentation behavior have been discussed with the comparison between the simulation and experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Modeling Sea Ice Effects for Wave Energy Resource Assessments

Wave-generated power has potential as a valuable coastal resource, but the wave climate needs to be mapped for feasibility before wave energy converters are installed. Numerical models are used for wave resource assessments to quantify the amount of available power and its seasonality. Alaska is the U.S. state with the longest coastline and has extensive wave resources, but it is affected by seasonal sea ice that dampens the wave energy and the full extent of this dampening is unknown. To accurately characterize the wave resource in regions that experience seasonal sea ice, coastal wave models must account for these effects. The aim of this study is to determine how the dampening effects of sea ice change wave energy resource assessments in the nearshore. Here, we show that by combining high-resolution sea ice imagery with a sea ice/wave dampening parameterization in an unstructured grid, the Simulating Waves Nearshore (SWAN) model improves wave height predictions and demonstrates the extent to which wave power decreases when sea ice is present. The sea ice parametrization decreases the bias and root mean square errors of wave height comparisons with two wave buoys and predicts a decrease in the wave power of up to 100 kW/m in areas around Prince William Sound, Alaska. The magnitude of the improvement of the model/buoy comparison depends on the coefficients used to parameterize the wave–ice interaction.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

Multiscale simulations of electron and ion dynamics in self-irradiated silicon

The interaction of energetic ions with the electronic and ionic system of target materials is an interesting but challenging multiscale problem, and understanding of the early stages after impact of heavy, initially charged ions is particularly poor. At the same time, energy deposition during these early stages determines later formation of damage cascades. In this work, we address the multiscale character by combining real-time time-dependent density functional theory for electron dynamics with molecular dynamics simulations of damage cascades. Our first-principles simulations prove that core electrons affect electronic stopping and have an unexpected influence on the charge state of the projectile. We show that this effect is absent for light projectiles, but dominates the stopping physics for heavy projectiles. By parametrizing an inelastic energy loss friction term in the molecular dynamics simulations using our first-principles results, we also show a qualitative influence of electronic stopping physics on radiation-damage cascades.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Estimating rates and patterns of diversification with incomplete sampling: a case study in the rosids

Premise Recent advances in generating large‐scale phylogenies enable broad‐scale estimation of species diversification. These now common approaches typically are characterized by (1) incomplete species coverage without explicit sampling methodologies and/or (2) sparse backbone representation, and usually rely on presumed phylogenetic placements to account for species without molecular data. We used empirical examples to examine the effects of incomplete sampling on diversification estimation and provide constructive suggestions to ecologists and evolutionary biologists based on those results. Methods We used a supermatrix for rosids and one well‐sampled subclade (Cucurbitaceae) as empirical case studies. We compared results using these large phylogenies with those based on a previously inferred, smaller supermatrix and on a synthetic tree resource with complete taxonomic coverage. Finally, we simulated random and representative taxon sampling and explored the impact of sampling on three commonly used methods, both parametric (RPANDA and BAMM) and semiparametric (DR). Results We found that the impact of sampling on diversification estimates was idiosyncratic and often strong. Compared to full empirical sampling, representative and random sampling schemes either depressed or inflated speciation rates, depending on methods and sampling schemes. No method was entirely robust to poor sampling, but BAMM was least sensitive to moderate levels of missing taxa. Conclusions We suggest caution against uncritical modeling of missing taxa using taxonomic data for poorly sampled trees and in the use of summary backbone trees and other data sets with high representative bias, and we stress the importance of explicit sampling methodologies in macroevolutionary studies.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Ripening of Rh Nanoparticle Catalysts in Reverse Water–Gas Shift via a Data-Driven Model Combining Physics, Theory, and Experiment

Degradation via sintering is an ongoing challenge that impedes the broad commercial success of supported metallic nanoparticle catalysts. To mitigate degradation via informed catalyst design and process operations, here we aim to disambiguate the underlying mechanisms of sintering by combining theory and experiment in a quantitative framework. While mechanistic sintering models exist, they only model a single sintering pathway, even though multiple sintering mechanisms can occur simultaneously or dominate at different stages of the process. Data-driven machine learning models have emerged as a means to represent complex processes through data regression. However, machine learning models have very large data needs and lack mechanistic insights due to their black-box encoding. To develop an interpretive model of catalyst degradation via sintering, we constructed a hybrid model combining mechanistic “physics-based” models and data-driven methods to obtain both reliable predictions and mechanistic insights regarding experimentally observed sintering phenomena. Focusing on nanoparticle sintering in the Rh–TiO 2 catalyst for the reverse water–gas shift (RWGS) reaction, the hybrid model couples a mechanistic term for Ostwald ripening with energy values calculated via density functional theory (DFT) with a parametric, data-driven discrepancy function term for unmodeled mechanisms. The hybrid model is trained using Bayesian inference with data collected from small-angle X-ray scattering (SAXS) in situ experiments wherein average nanoparticle diameter versus time was measured at three relevant operating temperatures. The calibrated hybrid model results show that an Ostwald ripening-only model parameterized with fixed DFT energies does not fully capture the time and temperature dependence of the SAXS-observed sintering kinetics, and that an additional functional contribution, or DFT energy calibration, is required to reconcile simulation and experiment. Analysis of the hybrid-model error confirms that the hybrid model outperforms both the purely mechanistic and purely data-driven alternatives in terms of expected predictive accuracy for time-evolving average particle sizes. Furthermore, the results support the hypothesis that the Ostwald ripening mechanism is less important for explaining the sintering phenomena as operating temperature increases under an assumed fixed DFT parameterization. This could be explained in one of two ways: either latent, unmodeled sintering mechanisms dominate at higher temperatures, or the DFT uncertainty increases with temperature. The proposed modeling approach directly links theory to experiments and simulations via a statistical hybrid modeling framework and can be extended to other catalytic systems to improve predictive models and mechanistic understanding.

Bayesian hybrid modeling↗

DOE OSTI2020

Effects of fuel composition and octane sensitivity on polycyclic aromatic hydrocarbon and soot emissions of gasoline–ethanol blend surrogates

The sooting propensity of a fuel is closely coupled with the fuel composition and chemistry. A detailed understanding of their effects is, therefore, needed to develop next-generation fuels which can minimize particulate emissions. With this overarching goal, the present work numerically investigates the effects of fuel composition and octane sensitivity (S) on polycyclic aromatic hydrocarbons (PAH) and soot emissions, for four-component gasoline-ethanol blend surrogates comprising isooctane, n-heptane, toluene, and ethanol. A partially-premixed counterflow flame is chosen as the canonical configuration for this study and simulations are performed using CHEMKIN-Pro-employing a kinetic mechanism developed by Park et al. (2017). In addition, a, detailed soot model based on the sectional method is used to capture the spatial characteristics of soot emissions. The kinetic mechanism and soot model are validated using available experimental data for various targets. A total of 86 TPRF-E mixtures, spanning a wide range of concentration of each component, and a wide range of S are analyzed. The effect of each non-paraffinic fuel component on the resultant PAH and soot emissions is investigated. PAH and soot emissions are found to vary significantly depending upon the blend composition. Additionally, based on the parametric sweeps, a regression analysis is carried out to identify global parameters that govern the formation of PAHs and soot. The analysis shows that both toluene content and S have a prominent effect on the formation of PAHs and soot, with toluene content having a stronger impact. Moreover, larger PAHs have higher dependency on toluene content and S. Furthermore, a detailed analysis is carried out to understand the physical and chemical phenomena associated with the observed trends of PAH and soot emissions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

From disorganized data to emergent dynamic models: Questionnaires to partial differential equations

Starting with sets of disorganized observations of spatially varying and temporally evolving systems, obtained at different (also disorganized) sets of parameters, we demonstrate the data-driven derivation of parameter dependent, evolutionary partial differential equation (PDE) models capable of generating the data. This tensor type of data is reminiscent of shuffled (multidimensional) puzzle tiles. The independent variables for the evolution equations (their “space” and “time”) as well as their effective parameters are all emergent , i.e. determined in a data-driven way from our disorganized observations of behavior in them. We use a diffusion map based questionnaire approach to build a smooth parametrization of our emergent space/time/parameter space for the data. This approach iteratively processes the data by successively observing them on the “space,” the “time” and the “parameter” axes of a tensor. Once the data become organized, we use machine learning (here, neural networks) to approximate the operators governing the evolution equations in this emergent space. Our illustrative examples are based (i) on a simple advection–diffusion model; (ii) on a previously developed vertex-plus-signaling model of Drosophila embryonic development; and (iii) on two complex dynamic network models (one neuronal and one coupled oscillator model) for which no obvious smooth embedding geometry is known a priori. This allows us to discuss features of the process like symmetry breaking, translational invariance, and autonomousness of the emergent PDE model, as well as its interpretability.

generative models↗

DOE OSTI2025

Using Calibrated Sodium Data for Preliminary Validation of the SRT Code for Advanced Reactors

Various types of non-light water reactors are currently engaged in the U.S. licensing process. Because of inherent differences compared with well-established large light water reactors, appropriate assessment tools are needed. Specifically, source term analysis, which determines environmental dose impacts from potential accident scenarios, is a crucial part of design and licensing. The U.S. Nuclear Regulatory Commission has emphasized the importance of mechanistic source term analysis for advanced reactor deployments. To align with these needs, Argonne National Laboratory has developed the Simplified Radionuclide Transport (SRT) source term analysis code for metal fuel Sodium-cooled Fast Reactors (SFRs) and microreactors. SRT conducts time-dependent radionuclide transport and retention in SFRs for core and ex-core radionuclide source accident sequences. The main objective of SRT is to provide rapid sensitivity and uncertainty analyses, incorporating parametric uncertainties and summarizing probabilistic results. As part of the code validation process, a study focused on the bubble scrubbing module was performed using an experiment recently carried out by the University of Wisconsin-Madison. Based on the analysis, the modeling approach in SRT provides accurate results for small and large aerosols, while slight underprediction of radionuclide aerosol removal are observed for medium sized aerosols. However, the deviation is minor, considering the highly uncertain phenomenon and range of results, and is in the conservative direction. In addition, uncertainty information derived from the experiments is further implemented, reflecting the actual span of parameters, which leads to enhanced agreement with code predictions. The results demonstrate that SRT provides reasonable predictions for the bubble scrubbing process in sodium pool.

Kam, Dong Hoon↗

DOE OSTI2026

Fantômas unconfined: global QCD fits with Bézier parameterizations

Fantômas is a C++ toolkit for exploring the parametrization dependence of parton distribution functions (PDFs) and other correlator functions in quantum chromodynamics (QCD). Fantômas facilitates the generation of adaptable polynomial parametrizations for PDFs, called metamorphs, to find best-fit PDF solutions and quantify the epistemic uncertainty associated with the parametrizations during their fitting. The method employs Bézier curves as universal approximators for a variety of PDF shapes. Integrated into the xFitter framework for the global QCD analysis, Fantômas provides a foundation for general models of PDFs, while reducing the computational time compared to the approaches utilizing traditional polynomial parametrizations as well as providing an interpretable alternative to neural-network-based models. This paper outlines the structure and practical usage of the Fantômas toolkit, including its inputs, outputs, and implementation within xFitter. It also provides a practical example of using Fantômas for uncertainty quantification as well as the combination of PDF fits into a single ensemble.

Bézier curves↗

DOE OSTI2025

Optimization problems governed by systems of PDEs with uncertainties

This paper reviews current theoretical and numerical approaches to optimization problems governed by partial differential equations (PDEs) that depend on random variables or random fields. Such problems arise in many engineering, science, economics and societal decision-making tasks. This paper focuses on problems in which the governing PDEs are parametrized by the random variables/fields, and the decisions are made at the beginning and are not revised once uncertainty is revealed. Examples of such problems are presented to motivate the topic of this paper, and to illustrate the impact of different ways to model uncertainty in the formulations of the optimization problem and their impact on the solution. A linear–quadratic elliptic optimal control problem is used to provide a detailed discussion of the set-up for the risk-neutral optimization problem formulation, study the existence and characterization of its solution, and survey numerical methods for computing it. Different ways to model uncertainty in the PDE-constrained optimization problem are surveyed in an abstract setting, including risk measures, distributionally robust optimization formulations, probabilistic functions and chance constraints, and stochastic orders. Furthermore, approximation-based optimization approaches and stochastic methods for the solution of the large-scale PDE-constrained optimization problems under uncertainty are described. Some possible future research directions are outlined.

Heinkenschloss, Matthias [Rice Univ., Houston, TX ↗