Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sensitivity formulations”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Reductions in wind farm main bearing rating lives resulting from wake impingement

This paper studies the impacts of wake impingement on main bearing rating lives predicted during the wind turbine design stage. A computational tool chain was developed to explore and quantify these effects across a wind farm populated by 10 MW wind turbines. Wind field and turbine load modelling was undertaken using the Dynamiks Python package, including application of a dynamic wake meandering model. The ISO 281 basic bearing rating life formulation was subsequently applied in order to evaluate impacts from wake effects. Analyses included a two-turbine parametric analysis, followed by a full wind farm analysis undertaken for the TotalControl 32-turbine reference wind farm, including full wind rose simulations across all operational wind speeds. Site conditions were accounted for using a Weibull wind speed distribution and a range of parametric wind direction rose models. Results indicate that wind farm main bearing rating lives are negatively impacted by the effects of wake impingement, resulting in rating life reductions for the analysed wind farm of the order of 16 % on average and as much as 20 %–25 %, both for the locating main bearing. Despite these high sensitivities, it is important to note that these resultant rating lives (i.e. the predicted lives) still far exceed the standard wind turbine operational lifetimes of 20–30 years. Wake impacts were also found to be asymmetrically related to the side on which the rotor is impinged, suggesting that, for the main bearing, there may be a “better” side for wake impingement to occur. Rating life sensitivities to wind rose shape were also observed. While these findings must be interpreted with due consideration for the various methodological limitations present, they provide compelling evidence that wake effects at the wind farm level should necessarily be included when undertaking main bearing operational load modelling, rating life assessment, or other load-related analyses.

17 WIND ENERGY↗
DOE OSTI2026

Extrusion‐Spheronization of Energetic Materials

The prevailing method to produce plastic‐bonded explosive (PBX) molding powder, or “prills”, is a complex, multiphase, and bespoke process that was developed by the high explosives (HEs) manufacturing industry several decades ago. This work demonstrates the utility of a simpler, widely‐used mechanical process—extrusion‐spheronization—to produce PBX prills. We begin by detailing precautions taken to enable safe remote operation of extrusion‐spheronization equipment with HE. We then perform a study investigating the effect of lacquer solvent composition on the particle packing, pressed density, and compressive strength properties of a 95 wt.% TATB/5 wt.% polymer binder formulation akin to PBX 9502. It was found that increased composition of low vapor pressure solvents caused prolonged retention of the solvent, resulting in tackier materials that would agglomerate and form larger prills. The larger prills also led to lower poured density and tapped density of HE prills and compressive strength of pressed PBX articles. The samples prepared with a 75% propyl acetate/25% butyl acetate lacquer solvent composition exhibited the highest compressive strength. However, it is believed that the prill packing and compressive strength properties are primarily driven by the prill size rather than the chemical composition of the lacquer itself. Extrusion‐spheronization remains a promising method to reliably and repeatably produce HE prills that is less sensitive to feedstock or process variation than traditional methods.

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

Multicriterion benefit evaluation methodology for safety enhancements in nuclear power plants and application for FLEX strategies

Here we propose a multicriterion benefit evaluation (MCBE) methodology for evaluating costs and benefits of implementing safety enhancements in nuclear power plants. The MCBE methodology establishes a comprehensive evaluation scope using multiple evaluation criteria (including public risk, occupational risk, plant revenue, and plant cost) and multiple contexts (including plant normal operations, incidents, and accidents). The MCBE methodology formulates a benefit or cost element as the product of a monetary consequence and its occurrence-likelihood change. The MCBE methodology also incorporates a decision maker’ preferences toward occurrence likelihood and monetary consequence into cost and benefit estimations through a modified usage of cumulative prospect theory. This paper applies the MCBE methodology to the implementation of diverse and flexible coping strategy (FLEX) in a generic nuclear power plant as an illustrative example. The FLEX benefits are evaluated through three types of influence paths: plant safety enhancement, risk-informed activities, and potential additional credit of FLEX in supporting plant operation and maintenance-efficiency improvement. The benefit-to-cost ratio (BCR) of FLEX implementation is calculated. Sensitivity analyses are conducted to examine the impacts of modeling assumptions and input parameter values on the BCR.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Dynamic analysis of fully constrained Cable-Driven Parallel Robots for automated prefabricated component installation

This paper presents a dynamic analysis and validation framework to assess a fully constrained six-anchor Cable-Driven Parallel Robot (CDPR) for automated installation of prefabricated facade components. Compared with conventional eight-anchor systems, the six-anchor configuration simplifies setup and reduces cost, but it also reduces control authority, shrinks the wrench-feasible workspace, and tightens orientation limits. Consequently, it is unclear a priori whether dynamically feasible trajectories exist to move the end effector from pickup to the facade. A constrained trajectory optimization is formulated to enforce the system dynamics, cable-tension bounds, and pose/velocity limits, and the framework is evaluated in simulation at three levels: (i) an idealized reference model, (ii) a lab-scale prototype model incorporating measured anchor misalignments and identified damping, and (iii) a full-scale three-story building model with load decomposition for structural feasibility checks. Across these scenarios, the analysis shows that optimal, constraint-satisfying trajectories exist that move the end effector from pickup to installation while maintaining a near-plumb, level orientation at the final pose. Collectively, this multi-scale dynamic analysis and validation framework supports the deployment readiness of the six-anchor CDPR and provides a prototype-based sensitivity case study of how measured anchor placement deviations affect feasibility.

CDPR↗
DOE OSTI2026

In-Band Scattering and Absorption of Infrared Blocking Foam Filters for Millimeter-wave Cameras

Expanded closed-cell polymer foams are widely used as thermal infrared (IR) blocking filters in millimeter-wave cameras, particularly for Cosmic Microwave Background observations. Precise knowledge of their millimeter-wave properties is essential for optimizing sensitivity. We present broadband (150 GHz - 2 THz) transmittance spectroscopy of Styroace-II and several Zotefoam filters, fitting their spectra with a radiative transfer model incorporating dielectric absorption and Rayleigh, Mie, and higher-order scattering. For a typical 5~cm thick filter stack at 280~GHz, Styroace-II exhibits ${\sim}10\%$ scattering with absorption estimated as ${\lesssim}5\%$ by effective-medium theory, while Zotefoam HD30 offers superior performance at ${\sim}3\%$ scattering and absorption likewise bounded to ${{\lesssim}0.3\%}$. Each model component is constrained at the ${\sim}0.1\%$ transmittance level for millimeter wavelengths. We observe batch-to-batch scattering variability of up to 2 percentage points in foams with multiple tested batches. Less commonly used Zotefoam formulations (LD15 and LD24) can further reduce in-band scattering to ${<}1\%$ while maintaining negligible in-band absorption and likely comparable IR blocking due to shared polyethylene absorption features and similar cell sizes. Based on this work, a filter constructed from the best measured LD24 batch has replaced the Styroace-II filter in a Simons Observatory 220/280 GHz Small Aperture Telescope.

Thomas, Alex [Chicago U., Astron. Astrophys. Ctr.;↗
DOE OSTI2026

Stochastic tensor contraction for quantum chemistry

Many computational methods in ab initio quantum chemistry are formulated in terms of high-order tensor contractions, whose cost determines the size of system that can be studied. We introduce stochastic tensor contraction to perform such operations with greatly reduced cost, and present its application to the gold-standard quantum chemistry method, coupled cluster theory with up to perturbative triples. For total energy errors more stringent than chemical accuracy, we reduce the computational scaling to that of mean-field theory, while starting to approach the mean-field absolute cost, thereby challenging the existing cost-to-accuracy landscape. Benchmarks against state-of-the-art local correlation approximations further show that we achieve an order-of-magnitude improvement in both total computation time and error, with significantly reduced sensitivity to system dimensionality and electron delocalization. We conclude that stochastic tensor contraction is a powerful computational primitive to accelerate a wide range of quantum chemistry.

Chemical Physics (physics.chem-ph)↗
DOE OSTI2024

Extreme Risk Mitigation in Reinforcement Learning using Extreme Value Theory

Risk-sensitive reinforcement learning (RL) has garnered significant attention in recent years due to the growing interest in deploying RL agents in real-world scenarios. A critical aspect of risk awareness involves modelling highly rare risk events (rewards) that could potentially lead to catastrophic outcomes. These infrequent occurrences present a formidable challenge for data-driven methods aiming to capture such risky events accurately. While risk-aware RL techniques do exist, they suffer from high variance estimation due to the inherent data scarcity. Our work proposes to enhance the resilience of RL agents when faced with very rare and risky events by focusing on refining the predictions of the extreme values predicted by the state-action value distribution. To achieve this, we formulate the extreme values of the state-action value function distribution as parameterized distributions, drawing inspiration from the principles of extreme value theory (EVT). We propose an extreme value theory based actor-critic approach, namely, Extreme Valued Actor-Critic (EVAC) which effectively addresses the issue of infrequent occurrence by leveraging EVT-based parameterization. Importantly, we theoretically demonstrate the advantages of employing these parameterized distributions in contrast to other risk-averse algorithms. Our evaluations show that the proposed method outperforms other risk averse RL algorithms on a diverse range of benchmark tasks, each encompassing distinct risk scenarios.

Wang, Yu↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fuel Property-Informed Process Design for the Direct Catalytic Conversion of Cellulosics

The direct catalytic conversion of cellulosics (DC3) with supercritical methanol presents a promising pathway to completely solubilize woody biomass and produce >60% mass yield of C2-C6 oxygenates for light duty vehicle applications. Although the oxygenate composition can be tuned by modifying catalyst and process conditions, limited efforts to date have examined DC3 fuel properties for light duty fuel applications to iteratively inform conversion process design. In this work, we designed and evaluated multi-component surrogate bioblendstocks that represent major DC3 light oxygenates and refined catalyst formulations and distillation parameters to improve the resulting light duty fuel quality. Surrogate design for this novel fuel production pathway was based on a preliminary product slate that was systematically adjusted to determine impacts on fuel properties of interest including heating value, oxidative stability, and heat of vaporization, among others. The resulting predictions and measurements of these surrogates were then used to inform improvement of conversion process, and specifically, reformulation of the catalyst which resulted in production of a higher energy density fuel product primarily comprised of saturated alcohols. Fuel testing outcomes also informed downstream oxygenate separations requirements. This resulted in the production of a promising DC3 light-duty fuel blendstock with >15% greater energy density relative to ethanol, while maintaining high octane number and octane sensitivity.

BIOMASS FUELS↗
DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Securing Federated Learning Against Active Reconstruction Attacks

Federated Learning (FL) has amassed notable attention for its ability to preserve user privacy while emphasizing the retainment of model training efficiency. Due to this potential, FL has been integrated in many domains, such as healthcare, finance, law, and industrial engineering, where data cannot be easily exchanged due to sensitive information and strict privacy laws. However, current research has indicated that FL protocols are easily compromised by active data reconstruction attacks employed by actively dishonest servers. The malicious modification of global model parameters allows an actively dishonest server to obtain a direct copy of users’ private data via gradient inversion. Here, this class of attacks is highly underexplored and continues to be a major challenge due to the intense threat model. In this paper, we propose OASIS as a scalable and modality-agnostic defense based on data augmentation that counteracts active data reconstruction attacks while preserving model performance. To generalize our defense, we uncover the intuition behind gradient inversion that enables these attacks and theoretically establish the conditions by which the defense can be considered robust regardless of attack design. From this, we formulate our defense with data augmentation that illustrates its ability to undermine the attack principle. We evaluate OASIS on five real-world datasets–two image-based (ImageNet and CIFAR100) and three text-based (Wikitext, Stack Overflow, and Shakespeare)–which span diverse uses cases such as vision tasks and language modeling. Comprehensive evaluations on these datasets exhibit the efficacy of OASIS and highlight its feasibility as a solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Refueling infrastructure planning in intercity networks considering route choice and travel time delay for mixed fleet of electric and conventional vehicles

The range anxiety has been a major factor that affects the market acceptance of electric vehicles. Even with the recent development of battery technologies, a lack of charging stations and range anxiety are still significant concerns, specifically for intercity trips. This calls for more investments in building charging stations and advancing battery technologies to increase the market share of electric vehicles and improve sustainability. This study suggests a configuration for plug-in electric vehicle charging infrastructure to support long-distance intercity trips of electric vehicles at the network level. A model is proposed to minimize the total system cost including infrastructure investment (building charging stations/spots) and travel time delays (charging time, waiting time in the queue, and detour time to access charging stations). This study fills existing gaps in the literature by capturing realistic patterns of travel demand and considering flow-dependent charging delays at charging stations. Furthermore, the proposed model, which is formulated as a mixed-integer program with nonlinear constraints, solves the optimization problem at the network level. At the network level, impacts of charging station locations on the traffic assignment problem with a mixed fleet of electric and conventional vehicles need to be considered. To this end, a traffic assignment module is integrated with a simulated annealing algorithm. The numerical experiments show a satisfactory application of the model for a full-scale case study (intercity network in Michigan). The solution quality and efficiency of the proposed solution algorithm are evaluated against those of an enumeration approach for a small case study. The results suggest that even for the current market share and charging stations’ setting, a significant investment is needed to support intercity trips without range anxiety issues and with acceptable delays. Additionally, through sensitivity analyses, the required infrastructure and battery investments to support intercity trips with acceptable delays are established for hypothetical increased market shares and battery size in the future.

42 ENGINEERING↗
DOE OSTI2021

Gradient sensing via cell communication

Experimental evidence lends support to the conjecture that cell-to-cell communication plays a role in the gradient sensing of chemical species by certain chains of cells. Models have been formulated to explore this idea. For cells with no identifiable sensing structure, Mugler et al. [Proc. Natl. Acad. Sci. (U.S.A.) 113, E689 (2016)] have defined a particular local excitation, global inhibition (LEGI) model that pits nearest-neighbor communication against local reactions in a noisy environment to suggest how this sensing capability might arise in a physical system. In this study, we generalize the nearest-neighbor communication mechanism in the aforementioned LEGI model in order to explore the extent to which the gradient sensing characteristics depend on the parametrization of the communication itself, as well as on the cell size, the radius of influence of neighboring cells, and the influence of the background noise. Using our generalization and a collection of particular candidate communication models, we find that the precision of gradient sensing is indeed sensitive to the particular communication model, and we derive physical and analytic explanations for these results. The framework established and the associated results should prove useful in understanding the appropriateness of particular cell-to-cell communication models in gradient sensing studies.

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

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Defining Reactivity–Deconstructability Relationships for Copolymerizations Involving Cleavable Comonomer Additives

The incorporation of cleavable comonomers as additives into polymers can imbue traditional polymers with controlled deconstructability and expanded end-of-life options. The efficiency with which cleavable comonomer additives (CCAs) can enable deconstruction is sensitive to their local distribution within a copolymer backbone, which is dictated by their copolymerization behavior. While qualitative heuristics exist that describe deconstructability, comprehensive quantitative connections between CCA loadings, reactivity ratios, polymerization mechanisms, and deconstruction reactions on the deconstruction efficiency of copolymers containing CCAs have not been established. Here, in this work, we broadly define these relationships using stochastic simulations characterizing various polymerization mechanisms (e.g., controlled/living, free-radical, and reversible ring-opening polymerizations), reactivity ratio pairs (spanning 2 orders of magnitude between 0.01 and 100), CCA loadings (2.5% to 20%), and deconstruction reactions (e.g., comonomer sequence-dependent deconstruction behavior). We show general agreement between simulated and experimentally observed deconstruction fragment sizes from the literature, demonstrating the predictive power of the methods used herein. These results will guide the development of more efficient CCAs and inform the formulation of deconstructable materials.

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

Topology and spectral entanglement in cavity-mediated photon scattering

Here, we develop a microscopic diagrammatic theory for cavity-mediated photon scattering in a topological one-dimensional insulator described by the Su–Schrieffer–Heeger model. Within the velocity-gauge formulation, we derive the photon self-energy and vertex corrections arising from virtual electron–hole excitations coupled to a quantized cavity mode, and we evaluate the resulting polariton dispersion and two-photon correlation spectra. Our analysis shows that vacuum fluctuations of the cavity field induce a momentum-resolved self-energy that mixes conduction and valence bands through virtual photon exchange, producing interband hybridization and avoided crossings in the electronic dispersion. This “cavity dressing” is symmetry-dependent, vanishing at the Brillouin-zone edge where the dipole matrix element is zero, and its strength is controlled by the spatial coherence range ζ ≈ (l c /a) 2 of virtual excitations. We further examine how the cavity modifies nonlinear optical observables, including the Kerr nonlinearity and biphoton spectral entanglement, and identify the regimes where these effects become sensitive to the underlying topological phase. The theoretical framework established here provides a unified description of light–matter coupling in topological and polaritonic systems, bridging solid-state cavity QED with the emerging field of cavity-modified quantum materials. Our results suggest that engineered photonic environments can coherently reshape the electronic landscape of topological insulators, offering new routes to control collective electronic and optical phenomena through vacuum-field fluctuations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

A Review of Quantum Computing Technologies in Power System Optimization

As modern power grids increasingly integrate variable renewable generation, distributed energy resources, and energy storage systems, classical optimization techniques are facing unprecedented challenges. This review examines the emerging application of quantum computing to overcome these challenges in power system optimization, including optimal power flow (OPF), unit commitment (UC), economic dispatch (ED), and intelligent switching and topology optimization (IS-TO). Recent research has introduced various quantum methodologies—such as gate-based, annealing-based, variational algorithms, and quantum-inspired algorithms—to address the combinatorial complexity inherent in grid reconfiguration and energy management. The review summaries the quantum algorithms, quantum devices and the power system test cases, highlighting hybrid quantum–classical strategies that leverage the complementary strengths of both paradigms. Some quantum advantages have been observed, including theoretical speedup, accurate simulation results, scalable qubit usage, efficient QUBO mapping. In particular, the review emphasizes the importance of integrating quantum optimization techniques with classical control frameworks, these hybrid approaches demonstrate the potential to improve real-time grid management and operational reliability. A significant portion of the analysis is devoted to the practical limitations of current quantum devices. Present-day quantum hardware, operating in the noisy intermediate-scale quantum (NISQ) era, remains highly sensitive to noise and limited in qubit connectivity, which constrains the scale and accuracy of implemented algorithms. The review delves into specific challenges such as the need for qubit-efficient encoding techniques and error mitigation strategies that are critical for handling real-world grid optimization problems. In addition, the work draws attention to the performance discrepancies between theoretical quantum speedups and experimental validations, underscoring the importance of rigorous benchmark studies using representative power grid test cases. In summary, this review highlights both the promise and limitations of quantum computing for power system optimization. It provides a comprehensive overview of the state-of-the-art technologies, categorizes recent advancements in algorithm design, and discusses practical considerations for implementation, and serves as an informative resource on current research. Future research directions include developing robust hybrid frameworks, advancing qubit-efficient formulations, and scaling up experimental demonstrations to confirm the theoretical advantages of quantum methods in large-scale power system operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Stochastic Framework for Optimal Control of Planetary Reentry Trajectories Under Multilevel Uncertainties

We present a novel stochastic optimal control framework that accounts for various types of uncertainties, with application to reentry trajectory planning. The formulation of the optimal trajectory control problem is presented in the context of an indirect method where a functional objective associated with the terminal vehicle speed is to be minimized. Uncertain input parameters in the optimal trajectory control model, including aerodynamic parameters and initial and terminal conditions, are modeled as aleatory random variables, while the statistical parameters of these aleatory distributions are themselves random variables. The parametric and model uncertainties are simultaneously propagated through an extended polynomial chaos expansion (EPCE) formalism. Several metrics are described to evaluate response statistics and presented as insightful tools for robust decision making. Specifically, the response probability density function (PDF) reflecting influence of both epistemic and aleatory uncertainties is obtained. By sampling over the random variables representing model error, an ensemble of response PDFs is generated and the associated failure probability is estimated as a random variable with its own polynomial chaos expansion. Besides, the sensitivity index functions of response PDF with respect to the statistical parameters are evaluated. Coupling parametric and model uncertainties within the EPCE framework leads to a robust and efficient paradigm for multilevel uncertainty propagation and PDF characterization in general optimal control problems.

Engineering↗