skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems

Abstract

We explore the relation between the development of a non-negligible probability of negative states and the instability of numerical integration of the intrusive Galerkin ordinary differential equation system describing uncertain chemical ignition. To prevent this instability without resorting to either multi-element local polynomial chaos (PC) methods or increasing the order of the PC representation in time, we propose a procedure aimed at modifying the amplitude of the PC modes to bring the probability of negative state values below a user-defined threshold. This modification can be effectively described as a filtering procedure of the spectral PC coefficients, which is applied on-the-fly during the numerical integration when the current value of the probability of negative states exceeds the prescribed threshold. We demonstrate the filtering procedure using a simple model of an ignition process in a batch reactor. This is carried out by comparing different observables and error measures as obtained by non-intrusive Monte Carlo and Gauss-quadrature integration and the filtered intrusive procedure. The filtering procedure has been shown to effectively stabilize divergent intrusive solutions, and also to improve the accuracy of stable intrusive solutions which are close to the stability limits.

Authors:
 [1];  [2]
  1. Sandia National Laboratories, Livermore, CA 94551 (United States)
  2. Mechanical and Aerospace Engineering Dept., Sapienza University of Rome, Rome (Italy)
Publication Date:
OSTI Identifier:
22314884
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 270; Other Information: Copyright (c) 2014 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; CHAOS THEORY; COMPARATIVE EVALUATIONS; DIFFERENTIAL EQUATIONS; GAUSS FUNCTION; MATHEMATICAL SOLUTIONS; MONTE CARLO METHOD; POLYNOMIALS; PROBABILITY; QUADRATURES; STABILITY

Citation Formats

Najm, Habib N., and Valorani, Mauro. Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems. United States: N. p., 2014. Web. doi:10.1016/J.JCP.2014.03.061.
Najm, Habib N., & Valorani, Mauro. Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems. United States. doi:10.1016/J.JCP.2014.03.061.
Najm, Habib N., and Valorani, Mauro. Fri . "Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems". United States. doi:10.1016/J.JCP.2014.03.061.
@article{osti_22314884,
title = {Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems},
author = {Najm, Habib N. and Valorani, Mauro},
abstractNote = {We explore the relation between the development of a non-negligible probability of negative states and the instability of numerical integration of the intrusive Galerkin ordinary differential equation system describing uncertain chemical ignition. To prevent this instability without resorting to either multi-element local polynomial chaos (PC) methods or increasing the order of the PC representation in time, we propose a procedure aimed at modifying the amplitude of the PC modes to bring the probability of negative state values below a user-defined threshold. This modification can be effectively described as a filtering procedure of the spectral PC coefficients, which is applied on-the-fly during the numerical integration when the current value of the probability of negative states exceeds the prescribed threshold. We demonstrate the filtering procedure using a simple model of an ignition process in a batch reactor. This is carried out by comparing different observables and error measures as obtained by non-intrusive Monte Carlo and Gauss-quadrature integration and the filtered intrusive procedure. The filtering procedure has been shown to effectively stabilize divergent intrusive solutions, and also to improve the accuracy of stable intrusive solutions which are close to the stability limits.},
doi = {10.1016/J.JCP.2014.03.061},
journal = {Journal of Computational Physics},
number = ,
volume = 270,
place = {United States},
year = {Fri Aug 01 00:00:00 EDT 2014},
month = {Fri Aug 01 00:00:00 EDT 2014}
}
  • We explore the relation between the development of a non-negligible probability of negative states and the instability of numerical integration of the intrusive Galerkin ordinary differential equation system describing uncertain chemical ignition. To prevent this instability without resorting to either multi-element local polynomial chaos (PC) methods or increasing the order of the PC representation in time, we propose a procedure aimed at modifying the amplitude of the PC modes to bring the probability of negative state values below a user-defined threshold. This modification can be effectively described as a filtering procedure of the spectral PC coefficients, which is applied on-the-flymore » during the numerical integration when the current value of the probability of negative states exceeds the prescribed threshold. We demonstrate the filtering procedure using a simple model of an ignition process in a batch reactor. This is carried out by comparing different observables and error measures as obtained by non-intrusive Monte Carlo and Gauss-quadrature integration and the filtered intrusive procedure. Lastly, the filtering procedure has been shown to effectively stabilize divergent intrusive solutions, and also to improve the accuracy of stable intrusive solutions which are close to the stability limits.« less
  • This paper deals with stochastic spectral methods for uncertainty propagation and quantification in nonlinear hyperbolic systems of conservation laws. We consider problems with parametric uncertainty in initial conditions and model coefficients, whose solutions exhibit discontinuities in the spatial as well as in the stochastic variables. The stochastic spectral method relies on multi-resolution schemes where the stochastic domain is discretized using tensor-product stochastic elements supporting local polynomial bases. A Galerkin projection is used to derive a system of deterministic equations for the stochastic modes of the solution. Hyperbolicity of the resulting Galerkin system is analyzed. A finite volume scheme with amore » Roe-type solver is used for discretization of the spatial and time variables. An original technique is introduced for the fast evaluation of approximate upwind matrices, which is particularly well adapted to local polynomial bases. Efficiency and robustness of the overall method are assessed on the Burgers and Euler equations with shocks.« less
  • Software-based digital systems in nuclear applications offer many potential benefits in the fields of safety, functionality, flexibility, and control, but they also present substantial challenges in demonstrating software reliability. In at least one nuclear system, serious concerns over the protection-system software have been raised. Achieving the required high level of software dependability through techniques such as testing, inspections, or mathematical verification is difficult because of the quantity and complexity of the software. The goal of the research described here is to facilitate dependability analysis by using a novel kernel software architecture. The kernel encapsulates into a relatively small piece ofmore » software the implementation of a set of critical safety policies so that policy enforcement is isolated from the rest of the system. Provided the kernel operates correctly, safety policy compliance is assured irrespective of the actions of the majority of the software.« less
  • General-purpose numerical methods for ordinary initial value problems are usually found in the class of linear multistage multivalue methods, first formulated by J. C. Butcher. Among these the explicit methods are easiest to implement. For this reason there has been considerable research activity devoted to generating methods of this class which utilize independent function evaluations that can be performed in parallel. Each such group of concurrent function evaluations can be regarded as a stage of the method. It is desirable that methods have a large region of stability so that stability-dictated restrictions on the stepsize will occur less frequently. Inmore » order to do a fair comparison of stability regions, one should consider the scaled stability region. This is defined using not the nominal stepsize but rather the stepsize per unit of work. Each stage counts as one work unit; in the context of parallel computing each stage may involve any number of independent function evaluations. Very important results on scaled stability regions were obtained by Jeltsch and Nevanlinna, in the context of serial computation. They identify a type of method, actually a type of algebraic function, whose scaled stability region is optimal. This means that the region is not a proper subset of some other scaled stability region. The distinguished type of algebraic functions is one that satisfies Property C, which requires that the principal'' branch of the algebraic function be analytic outside the stability region except for a pole at infinity whose multiplicity equals the number of work units of the method. It is shown that by modifying Property C the result on optimal stability regions holds also in the context of parallel computation. An important consequence is that, although the use of parallelism enlarges the class of scaled stability regions, it does not result in scaled stability regions that properly contain those obtained without parallelism.« less