DOE PAGES 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. 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.

Authors:
 [1];  [2]
  1. Sandia National Lab. (SNL-CA), Livermore, CA (United States)
  2. Sapienza Univ. of Rome (Italy)
Publication Date:
Research Org.:
Sandia National Lab. (SNL-CA), Livermore, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Basic Energy Sciences (BES)
OSTI Identifier:
1106517
Report Number(s):
SAND2013-6540J
Journal ID: ISSN 0021-9991; PII: S0021999114002484
Grant/Contract Number:  
AC04-94AL85000
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 270; Related Information: Proposed for publication in Journal of Computational Physics.; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; Uncertainty quanti cation; polynomial chaos; intrusive; ordinary di erential equations; 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. https://doi.org/10.1016/j.jcp.2014.03.061
Najm, Habib N., and Valorani, Mauro. Sat . "Enforcing positivity in intrusive PC-UQ methods for reactive ODE systems". United States. https://doi.org/10.1016/j.jcp.2014.03.061. https://www.osti.gov/servlets/purl/1106517.
@article{osti_1106517,
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. 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.},
doi = {10.1016/j.jcp.2014.03.061},
journal = {Journal of Computational Physics},
number = ,
volume = 270,
place = {United States},
year = {Sat Apr 12 00:00:00 EDT 2014},
month = {Sat Apr 12 00:00:00 EDT 2014}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 9 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Making best use of model evaluations to compute sensitivity indices
journal, May 2002


Local and Global Uncertainty Analyses of a Methane Flame Model
journal, November 2005

  • Zádor, Judit; Zsély, István Gy.; Turányi, Tamás
  • The Journal of Physical Chemistry A, Vol. 109, Issue 43
  • DOI: 10.1021/jp053270i

The role of sensitivity and uncertainty analysis in combustion modelling
journal, January 2013


The Homogeneous Chaos
journal, October 1938

  • Wiener, Norbert
  • American Journal of Mathematics, Vol. 60, Issue 4
  • DOI: 10.2307/2371268

A Stochastic Projection Method for Fluid Flow
journal, November 2001

  • Le Maı̂ tre, Olivier P.; Knio, Omar M.; Najm, Habib N.
  • Journal of Computational Physics, Vol. 173, Issue 2
  • DOI: 10.1006/jcph.2001.6889

The Wiener--Askey Polynomial Chaos for Stochastic Differential Equations
journal, January 2002


Numerical Challenges in the Use of Polynomial Chaos Representations for Stochastic Processes
journal, January 2004

  • Debusschere, Bert J.; Najm, Habib N.; Pébay, Philippe P.
  • SIAM Journal on Scientific Computing, Vol. 26, Issue 2
  • DOI: 10.1137/S1064827503427741

Uncertainty Quantification and Polynomial Chaos Techniques in Computational Fluid Dynamics
journal, January 2009


Probabilistic characterization of transport in heterogeneous media
journal, June 1998


Stochastic Finite Element Analysis for Multiphase Flow in Heterogeneous Porous Media
journal, January 1998


Adaptive polynomial chaos expansions applied to statistics of extremes in nonlinear random vibration
journal, April 1998


Incorporation of parametric uncertainty into complex kinetic mechanisms: Application to hydrogen oxidation in supercritical water
journal, January 1998


A Stochastic Projection Method for Fluid Flow
journal, September 2002

  • Le Maı̂tre, Olivier P.; Reagan, Matthew T.; Najm, Habib N.
  • Journal of Computational Physics, Vol. 181, Issue 1
  • DOI: 10.1006/jcph.2002.7104

A multigrid solver for two-dimensional stochastic diffusion equations
journal, October 2003

  • Le Maı̂tre, O. P.; Knio, O. M.; Debusschere, B. J.
  • Computer Methods in Applied Mechanics and Engineering, Vol. 192, Issue 41-42
  • DOI: 10.1016/S0045-7825(03)00457-2

A new stochastic approach to transient heat conduction modeling with uncertainty
journal, November 2003


Modeling uncertainty in flow simulations via generalized polynomial chaos
journal, May 2003


Uncertainty quantification in reacting-flow simulations through non-intrusive spectral projection
journal, February 2003


A Two-Scale Nonperturbative Approach to Uncertainty Analysis of Diffusion in Random Composites
journal, January 2004

  • Xiu, Dongbin; Tartakovsky, Daniel M.
  • Multiscale Modeling & Simulation, Vol. 2, Issue 4
  • DOI: 10.1137/03060268X

Noisy Inflows Cause a Shedding-Mode Switching in Flow Past an Oscillating Cylinder
journal, April 2004


Spectral stochastic uncertainty quantification in chemical systems
journal, September 2004

  • Reagan, M. T.; Najm 4, H. N.; Debusschere, B. J.
  • Combustion Theory and Modelling, Vol. 8, Issue 3
  • DOI: 10.1088/1364-7830/8/3/010

Variational multiscale stabilized FEM formulations for transport equations: stochastic advection–diffusion and incompressible stochastic Navier–Stokes equations
journal, January 2005

  • Badri Narayanan, Velamur Asokan; Zabaras, Nicholas
  • Journal of Computational Physics, Vol. 202, Issue 1
  • DOI: 10.1016/j.jcp.2004.06.019

Using stochastic analysis to capture unstable equilibrium in natural convection
journal, September 2005

  • Asokan, Badrinarayanan Velamur; Zabaras, Nicholas
  • Journal of Computational Physics, Vol. 208, Issue 1
  • DOI: 10.1016/j.jcp.2005.02.007

A stochastic variational multiscale method for diffusion in heterogeneous random media
journal, November 2006

  • Velamur Asokan, Badrinarayanan; Zabaras, Nicholas
  • Journal of Computational Physics, Vol. 218, Issue 2
  • DOI: 10.1016/j.jcp.2006.02.026

Long-term behavior of polynomial chaos in stochastic flow simulations
journal, August 2006

  • Wan, Xiaoliang; Karniadakis, George Em
  • Computer Methods in Applied Mechanics and Engineering, Vol. 195, Issue 41-43
  • DOI: 10.1016/j.cma.2005.10.016

Stochastic analysis of transport in tubes with rough walls
journal, September 2006


Quantifying uncertainty in chemical systems modeling
journal, January 2005

  • Reagan, M. T.; Najm, H. N.; Pébay, P. P.
  • International Journal of Chemical Kinetics, Vol. 37, Issue 6
  • DOI: 10.1002/kin.20081

Stochastic Approaches to Uncertainty Quantification in CFD Simulations
journal, March 2005

  • Mathelin, Lionel; Hussaini, M. Yousuff; Zang, Thomas A.
  • Numerical Algorithms, Vol. 38, Issue 1
  • DOI: 10.1007/s11075-004-2866-z

Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations
journal, April 2005

  • Matthies, Hermann G.; Keese, Andreas
  • Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 12-16
  • DOI: 10.1016/j.cma.2004.05.027

A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data
journal, January 2007

  • Babuška, Ivo; Nobile, Fabio; Tempone, Raúl
  • SIAM Journal on Numerical Analysis, Vol. 45, Issue 3
  • DOI: 10.1137/050645142

Sparse grid collocation schemes for stochastic natural convection problems
journal, July 2007

  • Ganapathysubramanian, Baskar; Zabaras, Nicholas
  • Journal of Computational Physics, Vol. 225, Issue 1
  • DOI: 10.1016/j.jcp.2006.12.014

A Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data
journal, January 2008

  • Nobile, F.; Tempone, R.; Webster, C. G.
  • SIAM Journal on Numerical Analysis, Vol. 46, Issue 5
  • DOI: 10.1137/060663660

An Anisotropic Sparse Grid Stochastic Collocation Method for Partial Differential Equations with Random Input Data
journal, January 2008

  • Nobile, F.; Tempone, R.; Webster, C. G.
  • SIAM Journal on Numerical Analysis, Vol. 46, Issue 5
  • DOI: 10.1137/070680540

An adaptive hierarchical sparse grid collocation algorithm for the solution of stochastic differential equations
journal, May 2009


The multi-element probabilistic collocation method (ME-PCM): Error analysis and applications
journal, November 2008

  • Foo, Jasmine; Wan, Xiaoliang; Karniadakis, George Em
  • Journal of Computational Physics, Vol. 227, Issue 22
  • DOI: 10.1016/j.jcp.2008.07.009

Multi-element probabilistic collocation method in high dimensions
journal, March 2010


Hybrid Stochastic Finite Elements and Generalized Monte Carlo Simulation
journal, December 1998


Stochastic Finite-Element Analysis of Seismic Soil–Structure Interaction
journal, January 2002


Stochastic Response Surface Methods (SRSMs) for Uncertainty Propagation: Application to Environmental and Biological Systems
journal, June 1998


High-Order Collocation Methods for Differential Equations with Random Inputs
journal, January 2005

  • Xiu, Dongbin; Hesthaven, Jan S.
  • SIAM Journal on Scientific Computing, Vol. 27, Issue 3
  • DOI: 10.1137/040615201

Multiparameter Spectral Representation of Noise-Induced Competence in Bacillus Subtilis
journal, November 2012

  • Sargsyan, Khachik; Safta, Cosmin; Debusschere, Bert
  • IEEE/ACM Transactions on Computational Biology and Bioinformatics, Vol. 9, Issue 6
  • DOI: 10.1109/TCBB.2012.107

A Stochastic Multiscale Coupling Scheme to Account for Sampling Noise in Atomistic-to-Continuum Simulations
journal, January 2012

  • Salloum, Maher; Sargsyan, Khachik; Jones, Reese
  • Multiscale Modeling & Simulation, Vol. 10, Issue 2
  • DOI: 10.1137/110844404

A generalized spectral decomposition technique to solve a class of linear stochastic partial differential equations
journal, September 2007


Generalized spectral decomposition method for solving stochastic finite element equations: Invariant subspace problem and dedicated algorithms
journal, October 2008


A least-squares approximation of partial differential equations with high-dimensional random inputs
journal, July 2009


A non-adapted sparse approximation of PDEs with stochastic inputs
journal, April 2011


Finite elements for elliptic problems with stochastic coefficients
journal, February 2005

  • Frauenfelder, Philipp; Schwab, Christoph; Todor, Radu Alexandru
  • Computer Methods in Applied Mechanics and Engineering, Vol. 194, Issue 2-5
  • DOI: 10.1016/j.cma.2004.04.008

Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients
journal, April 2007

  • Todor, Radu Alexandru; Schwab, Christoph
  • IMA Journal of Numerical Analysis, Vol. 27, Issue 2
  • DOI: 10.1093/imanum/drl025

Uncertainty propagation using Wiener–Haar expansions
journal, June 2004

  • Le Maı̂tre, O. P.; Knio, O. M.; Najm, H. N.
  • Journal of Computational Physics, Vol. 197, Issue 1
  • DOI: 10.1016/j.jcp.2003.11.033

Multi-resolution analysis of Wiener-type uncertainty propagation schemes
journal, July 2004

  • Le Maı̂tre, O. P.; Najm, H. N.; Ghanem, R. G.
  • Journal of Computational Physics, Vol. 197, Issue 2
  • DOI: 10.1016/j.jcp.2003.12.020

Multi‐Resolution‐Analysis Scheme for Uncertainty Quantification in Chemical Systems
journal, January 2007

  • Le Maître, O. P.; Najm, H. N.; Pébay, P. P.
  • SIAM Journal on Scientific Computing, Vol. 29, Issue 2
  • DOI: 10.1137/050643118

An adaptive multi-element generalized polynomial chaos method for stochastic differential equations
journal, November 2005


Multi-Element Generalized Polynomial Chaos for Arbitrary Probability Measures
journal, January 2006

  • Wan, Xiaoliang; Karniadakis, George Em
  • SIAM Journal on Scientific Computing, Vol. 28, Issue 3
  • DOI: 10.1137/050627630

Uncertainty Quantification given Discontinuous Model Response and a Limited Number of Model Runs
journal, January 2012

  • Sargsyan, Khachik; Safta, Cosmin; Debusschere, Bert
  • SIAM Journal on Scientific Computing, Vol. 34, Issue 1
  • DOI: 10.1137/100817899

Spectral and multiresolution Wiener expansions of oscillatory stochastic processes
journal, July 2006


Uncertainty quantification of limit-cycle oscillations
journal, September 2006

  • Beran, Philip S.; Pettit, Chris L.; Millman, Daniel R.
  • Journal of Computational Physics, Vol. 217, Issue 1
  • DOI: 10.1016/j.jcp.2006.03.038

Wiener Chaos expansions and numerical solutions of randomly forced equations of fluid mechanics
journal, August 2006

  • Hou, Thomas Y.; Luo, Wuan; Rozovskii, Boris
  • Journal of Computational Physics, Vol. 216, Issue 2
  • DOI: 10.1016/j.jcp.2006.01.008

Modeling physical uncertainties in dynamic stall induced fluid–structure interaction of turbine blades using arbitrary polynomial chaos
journal, June 2007


An alternative unsteady adaptive stochastic finite elements formulation based on interpolation at constant phase
journal, December 2008

  • Witteveen, Jeroen A. S.; Bijl, Hester
  • Computer Methods in Applied Mechanics and Engineering, Vol. 198, Issue 3-4
  • DOI: 10.1016/j.cma.2008.09.005

Effect of randomness on multi-frequency aeroelastic responses resolved by Unsteady Adaptive Stochastic Finite Elements
journal, October 2009


Multiscale Stochastic Preconditioners in Non-intrusive Spectral Projection
journal, April 2011

  • Alexanderian, Alen; Le Maître, Oliver P.; Najm, Habib N.
  • Journal of Scientific Computing, Vol. 50, Issue 2
  • DOI: 10.1007/s10915-011-9486-2

Time-dependent generalized polynomial chaos
journal, November 2010

  • Gerritsma, Marc; van der Steen, Jan-Bart; Vos, Peter
  • Journal of Computational Physics, Vol. 229, Issue 22
  • DOI: 10.1016/j.jcp.2010.07.020

Stochastic Galerkin techniques for random ordinary differential equations
journal, April 2012


Eigenvalues of the Jacobian of a Galerkin-Projected Uncertain ODE System
journal, January 2011

  • Sonday, Benjamin E.; Berry, Robert D.; Najm, Habib N.
  • SIAM Journal on Scientific Computing, Vol. 33, Issue 3
  • DOI: 10.1137/100785922

On the convergence of generalized polynomial chaos expansions
journal, October 2011

  • Ernst, Oliver G.; Mugler, Antje; Starkloff, Hans-Jörg
  • ESAIM: Mathematical Modelling and Numerical Analysis, Vol. 46, Issue 2
  • DOI: 10.1051/m2an/2011045

Uncertainty quantification in chemical systems
journal, November 2009

  • Najm, H. N.; Debusschere, B. J.; Marzouk, Y. M.
  • International Journal for Numerical Methods in Engineering, Vol. 80, Issue 6‒7
  • DOI: 10.1002/nme.2551

Slow Manifold Structure in Explosive Kinetics. 1. Bifurcations of Points-at-Infinity in Prototypical Models
journal, November 2006

  • Creta, F.; Adrover, A.; Cerbelli, S.
  • The Journal of Physical Chemistry A, Vol. 110, Issue 50
  • DOI: 10.1021/jp0636064

Works referencing / citing this record:

On the Statistical Calibration of Physical Models: STATISTICAL CALIBRATION OF PHYSICAL MODELS
journal, February 2015

  • Sargsyan, K.; Najm, H. N.; Ghanem, R.
  • International Journal of Chemical Kinetics, Vol. 47, Issue 4
  • DOI: 10.1002/kin.20906