Computational uncertainty quantification for a clarifier-thickener model with several random perturbations: A hybrid stochastic Galerkin approach

Abstract Continuous sedimentation processes in a clarifier-thickener unit can be described by a scalar nonlinear conservation law whose flux density function is discontinuous with respect to the spatial position. In the applications of this model, which include mineral processing and wastewater treatment, the rate and composition of the feed flow cannot be given deterministically. Efficient numerical simulation is required to quantify the effect of uncertainty in these control parameters in terms of the response of the clarifier-thickener system. Thus, the problem at hand is one of uncertainty quantification for nonlinear hyperbolic problems with several random perturbations. The presented hybrid stochastic Galerkin method is devised so as to extend the polynomial chaos approximation by multiresolution discretization in the stochastic space. This approach leads to a deterministic hyperbolic system, which is partially decoupled and therefore suitable for efficient parallelisation. Stochastic adaptivity reduces the computational effort. Several numerical experiments are presented.

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

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

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

[4]  Aslak Tveito,et al.  The Solution of Nonstrictly Hyperbolic Conservation Laws May Be Hard to Compute , 1995, SIAM J. Sci. Comput..

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

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

[7]  Rémi Abgrall,et al.  A semi-intrusive deterministic approach to uncertainty quantification in non-linear fluid flow problems , 2013, J. Comput. Phys..

[8]  A. Vande Wouwer,et al.  Settler dynamic modeling and MATLAB simulation of the activated sludge process , 2009 .

[9]  Pak Ching Li,et al.  Construction of analytical solutions and numerical methods comparison of the ideal continuous settling model , 2015, Comput. Chem. Eng..

[10]  Ingmar Nopens,et al.  On reliable and unreliable numerical methods for the simulation of secondary settling tanks in wastewater treatment , 2012, Comput. Chem. Eng..

[11]  Kai Schneider,et al.  Adaptive Multiresolution Methods for the Simulation of Waves in Excitable Media , 2010, J. Sci. Comput..

[12]  Stefan Diehl,et al.  Operating Charts for Continuous Sedimentation II: Step Responses , 2005 .

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

[14]  W. T. Martin,et al.  The Orthogonal Development of Non-Linear Functionals in Series of Fourier-Hermite Functionals , 1947 .

[15]  A Vande Wouwer,et al.  Modeling and numerical simulation of secondary settlers: a method of lines strategy. , 2009, Water research.

[16]  Daniel G. Sbarbaro-Hofer,et al.  Simple mass balance controllers for continuous sedimentation , 2013, Comput. Chem. Eng..

[17]  D. Xiu,et al.  Modeling uncertainty in flow simulations via generalized polynomial chaos , 2003 .

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

[19]  John D. Towers,et al.  Well-posedness in BVt and convergence of a difference scheme for continuous sedimentation in ideal clarifier-thickener units , 2004, Numerische Mathematik.

[20]  Raimund Bürger,et al.  Discontinuous finite volume element discretization for coupled flow-transport problems arising in models of sedimentation , 2015, J. Comput. Phys..

[21]  P. Servais,et al.  Dynamic modeling of sludge compaction and consolidation processes in wastewater secondary settling tanks. , 2009, Water environment research : a research publication of the Water Environment Federation.

[22]  Alexandre Ern,et al.  Adaptive Anisotropic Spectral Stochastic Methods for Uncertain Scalar Conservation Laws , 2012, SIAM J. Sci. Comput..

[23]  G. J. Kynch.  A theory of sedimentation , 1952 .

[24]  G. Ekama,et al.  Assessing the applicability of the 1D flux theory to full-scale secondary settling tank design with a 2D hydrodynamic model. , 2004, Water research.

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

[26]  Siegfried Müller,et al.  Adaptive multiresolution discontinuous Galerkin schemes for conservation laws: multi-dimensional case , 2013, Math. Comput..

[27]  A. Harten.  Multiresolution algorithms for the numerical solution of hyperbolic conservation laws , 2010 .

[28]  Barna L. Bihari,et al.  Application of generalized wavelets: an adaptive multiresolution scheme , 1995 .

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

[30]  S. Osher,et al.  One-sided difference approximations for nonlinear conservation laws , 1981 .

[31]  J. Ph. Chancelier,et al.  Analysis of a Conservation PDE With Discontinuous Flux: A Model of Settler , 1994, SIAM J. Appl. Math..

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

[33]  Kai Schneider,et al.  A MULTIRESOLUTION METHOD FOR THE SIMULATION OF SEDIMENTATION IN INCLINED CHANNELS , 2012 .

[34]  Witold Nocon.  Mathematical modelling of distributed feed in continuous sedimentation , 2006, Simul. Model. Pract. Theory.

[35]  Stefan Diehl,et al.  Scalar conservation laws with discontinuous flux function: I. The viscous profile condition , 1996 .

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

[37]  Yabing Guo,et al.  Application Analysis of One-Dimensional Sedimentation Model , 2010, 2010 4th International Conference on Bioinformatics and Biomedical Engineering.

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

[39]  Benedek Gy Plósz,et al.  One-dimensional modelling of the secondary clarifier-factors affecting simulation in the clarification zone and the assessment of the thickening flow dependence. , 2007, Water research.

[40]  Monika Bargieł,et al.  Extension of the Richardson–Zaki equation to suspensions of multisized irregular particles , 2013 .

[41]  K. B. Oldham,et al.  An Atlas of Functions. , 1988 .

[42]  I Nopens,et al.  Impact on sludge inventory and control strategies using the benchmark simulation model no. 1 with the Bürger-Diehl settler model. , 2015, Water science and technology : a journal of the International Association on Water Pollution Research.

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

[44]  Stefan Diehl,et al.  Operating charts for continuous sedimentation IV: limitations for control of dynamic behaviour , 2008 .

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

[46]  Stefan Diehl,et al.  Operating Charts for Continuous Sedimentation III: Control of Step Inputs , 2006 .

[47]  Stefan Diehl,et al.  Operating charts for continuous sedimentation I: Control of steady states , 2001 .

[48]  Kai Schneider,et al.  Fully adaptive multiresolution schemes for strongly degenerate parabolic equations in one space dimension , 2008, 0807.0400.

[49]  Siegfried Müller,et al.  Adaptive Multiscale Schemes for Conservation Laws , 2002, Lecture Notes in Computational Science and Engineering.

[50]  Murat Kulahci,et al.  iCFD: Interpreted Computational Fluid Dynamics - Degeneration of CFD to one-dimensional advection-dispersion models using statistical experimental design - The secondary clarifier. , 2015, Water research.