Computational uncertainty quantification for some strongly degenerate parabolic convection-diffusion equations

Abstract Strongly degenerate parabolic convection–diffusion equations arise as governing equations in a number of applications such as traffic flow with driver reaction and anticipation distance and sedimentation of solid–liquid suspensions in mineral processing and wastewater treatment. In these applications several parameters that define the convective flux function and the degenerating diffusion coefficient are subject to stochastic variability. A method to evaluate the variability of the solution of the governing partial differential equation in response to that of the parameters is presented. To this end, a generalized polynomial chaos (gPC) expansion of the solution is approximated by its projection onto a finite-dimensional space of piecewise polynomial functions defined on a suitable discretization of the stochastic domain, according to the basic principle of the hybrid stochastic Galerkin (HSG) approach. This approach is combined with a finite volume (FV) method, resulting in a so-called FV–HSG method, to compute the sought deterministic coefficient functions of the truncated polynomial-chaos-based expansion of the solution. Since the stochastic parameter space is now spanned by piecewise polynomial functions, one may employ the numerical result to compute the reconstruction of the numerical solution for arbitrary values of the random variables. The expectation, the variance or other stochastic quantities of the solution (as functions of time and position) can also be computed from these coefficient functions. The method is illustrated by a number of numerical examples.

[1]  Raimund Bürger,et al.  A Model of Continuous Sedimentation of Flocculated Suspensions in Clarifier-Thickener Units , 2005, SIAM J. Appl. Math..

[2]  Hermann G. Matthies,et al.  Galerkin methods for linear and nonlinear elliptic stochastic partial differential equations , 2005 .

[3]  Stefan Diehl,et al.  A conservation Law with Point Source and Discontinuous Flux Function Modelling Continuous Sedimentation , 1996, SIAM J. Appl. Math..

[4]  E. Tadmor,et al.  New High-Resolution Central Schemes for Nonlinear Conservation Laws and Convection—Diffusion Equations , 2000 .

[5]  Raimund Bürger,et al.  Uncertainty Quantification for a Clarifier–Thickener Model with Random Feed , 2011 .

[6]  BabuskaIvo,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007 .

[7]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[8]  Raimund Bürger,et al.  Model equations for gravitational sedimentation-consolidation processes , 2000 .

[9]  B. Alpert A class of bases in L 2 for the sparse representations of integral operators , 1993 .

[10]  Ilse Smets,et al.  SENSITIVITY ANALYSIS OF A ONE-DIMENSIONAL CONVECTION-DIFFUSION MODEL FOR SECONDARY SETTLING TANKS , 2004 .

[11]  Christian Rohde,et al.  A stochastically and spatially adaptive parallel scheme for uncertain and nonlinear two-phase flow problems , 2015, Computational Geosciences.

[12]  Stefan Diehl,et al.  Numerical identification of constitutive functions in scalar nonlinear convection-diffusion equations with application to batch sedimentation , 2015 .

[13]  R. Ghanem,et al.  Multi-resolution analysis of wiener-type uncertainty propagation schemes , 2004 .

[14]  Christian Rohde,et al.  Stochastic Modeling for Heterogeneous Two-Phase Flow , 2014 .

[15]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[16]  Guodong Wang,et al.  Semidiscrete central-upwind scheme for conservation laws with a discontinuous flux function in space , 2011, Appl. Math. Comput..

[17]  Fabio Nobile,et al.  A Stochastic Collocation Method for Elliptic Partial Differential Equations with Random Input Data , 2007, SIAM Rev..

[18]  P. Nelson Traveling-wave solutions of the diffusively corrected kinematic-wave model , 2002 .

[19]  Ingmar Nopens,et al.  A consistent modelling methodology for secondary settling tanks: a reliable numerical method. , 2013, Water science and technology : a journal of the International Association on Water Pollution Research.

[20]  Andrea Barth,et al.  Multi-level Monte Carlo Finite Element method for elliptic PDEs with stochastic coefficients , 2011, Numerische Mathematik.

[21]  Raimund Bürger,et al.  Sedimentation and Thickening , 1999 .

[22]  P. I. Richards Shock Waves on the Highway , 1956 .

[23]  Bruno Després,et al.  Uncertainty quantification for systems of conservation laws , 2009, J. Comput. Phys..

[24]  Siddhartha Mishra,et al.  Sparse tensor multi-level Monte Carlo finite volume methods for hyperbolic conservation laws with random initial data , 2012, Math. Comput..

[25]  M. Ratto,et al.  Sensitivity and uncertainty quantification techniques applied to systems of conservation laws , 2010 .

[26]  B. Alpert Wavelets and other bases for fast numerical linear algebra , 1993 .

[27]  G. Karniadakis,et al.  An adaptive multi-element generalized polynomial chaos method for stochastic differential equations , 2005 .

[28]  Alexandre Ern,et al.  Intrusive Galerkin methods with upwinding for uncertain nonlinear hyperbolic systems , 2010, J. Comput. Phys..

[29]  Raimund Bürger,et al.  Convexity-preserving flux identification for scalar conservation laws modelling sedimentation , 2013 .

[30]  R. Bürger,et al.  Efficient parameter estimation in a macroscopic traffic flow model by discrete mollification , 2015 .

[31]  B D Greenshields,et al.  A study of traffic capacity , 1935 .

[32]  G. Karniadakis,et al.  Multi-Element Generalized Polynomial Chaos for Arbitrary Probability Measures , 2006, SIAM J. Sci. Comput..

[33]  Raimund Bürger,et al.  Mathematical model and numerical simulation of the dynamics of flocculated suspensions in clarifier–thickeners , 2005 .

[34]  Alexander Kurganov,et al.  Central‐upwind schemes on triangular grids for hyperbolic systems of conservation laws , 2005 .

[35]  Raimund Bürger,et al.  A hybrid stochastic Galerkin method for uncertainty quantification applied to a conservation law modelling a clarifier‐thickener unit , 2014 .

[36]  Ingmar Nopens,et al.  Extending and calibrating a mechanistic hindered and compression settling model for activated sludge using in-depth batch experiments. , 2008, Water research.

[37]  M J Lighthill,et al.  On kinematic waves II. A theory of traffic flow on long crowded roads , 1955, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[38]  Raimund Bürger,et al.  ON A DIFFUSIVELY CORRECTED KINEMATIC-WAVE TRAFFIC FLOW MODEL WITH CHANGING ROAD SURFACE CONDITIONS , 2003 .

[39]  Raimund Bürger,et al.  Computational uncertainty quantification for a clarifier-thickener model with several random perturbations: A hybrid stochastic Galerkin approach , 2016, Comput. Chem. Eng..

[40]  G. Gagneux,et al.  Solution forte entropique de lois scalaires hyperboliques-paraboliques dégénérées , 1999 .