Abstract:The continuous sedimentation process in a clarifier-thickener can be described by a scalar nonlinear conservation law for the local solids volume fraction. The flux density function is discontinuous with respect to spatial position due to feed and discharge mechanisms. Typically, the feed flow cannot be given deterministically and efficient numerical simulation requires a concept for quantifying uncertainty. In this paper uncertainty quantification is expressed by a new hybrid stochastic Galerkin (HSG) method … Show more
“…This is quite different from the approach taken in previous work, e.g. [23,22], and therefore their hyperbolicity analysis cannot be extended to this case.…”
Section: Contributions Of This Workmentioning
confidence: 62%
“…Sampling based generalized polynomial chaos methods, such as stochastic collocation [21] suffer from the curse of dimensionality, and they become infeasible for large problems due to the prohibitive computational cost. With the continuous growth of computer power, stochastic Galerkin methods including the efficient adaptive and parallelized hybrid stochastic Galerkin solver in [22], have been gaining popularity as a powerful alternative to sampling based methods.…”
The generalized polynomial chaos method is applied to the Buckley-Leverett equation. We consider a spatially homogeneous domain modeled as a random field. The problem is projected onto stochastic basis functions which yields an extended system of partial differential equations. Analysis and numerical methods leading to reduced computational cost are presented for the extended system of equations.The accurate representation of the evolution of a discontinuous stochastic solution over time requires a large number of stochastic basis functions. Adaptivity of the stochastic basis to reduce computational cost is challenging in the stochastic Galerkin setting since the change of basis affects the system matrix itself. To achieve adaptivity without adding overhead by rewriting the entire system of equations for every grid cell, we devise a basis reduction method that distinguishes between locally significant and insignificant modes without changing the actual system matrices.Results are presented for problems in one and two spatial dimensions, with varying number of stochastic dimensions. We show how to obtain stochastic velocity fields from realistic permeability fields and demonstrate the performance of the stochastic Galerkin method with local basis reduction. The system of conservation laws is discretized with a finite volume method and we demonstrate numerical convergence to the reference solution obtained *
“…This is quite different from the approach taken in previous work, e.g. [23,22], and therefore their hyperbolicity analysis cannot be extended to this case.…”
Section: Contributions Of This Workmentioning
confidence: 62%
“…Sampling based generalized polynomial chaos methods, such as stochastic collocation [21] suffer from the curse of dimensionality, and they become infeasible for large problems due to the prohibitive computational cost. With the continuous growth of computer power, stochastic Galerkin methods including the efficient adaptive and parallelized hybrid stochastic Galerkin solver in [22], have been gaining popularity as a powerful alternative to sampling based methods.…”
The generalized polynomial chaos method is applied to the Buckley-Leverett equation. We consider a spatially homogeneous domain modeled as a random field. The problem is projected onto stochastic basis functions which yields an extended system of partial differential equations. Analysis and numerical methods leading to reduced computational cost are presented for the extended system of equations.The accurate representation of the evolution of a discontinuous stochastic solution over time requires a large number of stochastic basis functions. Adaptivity of the stochastic basis to reduce computational cost is challenging in the stochastic Galerkin setting since the change of basis affects the system matrix itself. To achieve adaptivity without adding overhead by rewriting the entire system of equations for every grid cell, we devise a basis reduction method that distinguishes between locally significant and insignificant modes without changing the actual system matrices.Results are presented for problems in one and two spatial dimensions, with varying number of stochastic dimensions. We show how to obtain stochastic velocity fields from realistic permeability fields and demonstrate the performance of the stochastic Galerkin method with local basis reduction. The system of conservation laws is discretized with a finite volume method and we demonstrate numerical convergence to the reference solution obtained *
“…Due to the relative entropy framework we expect our theory to be extendable to this case. For further applications of our method, the construction of space-stochastic adaptive schemes using the Hybrid Stochastic Galerkin method (Bürger et al (2014)) and the residuals as local indicators, will be considered.…”
Section: Discussionmentioning
confidence: 99%
“…For an overview on recent work on Uncertainty Quantification for hyperbolic equations see Bijl et al (2013); Le Maître & Knio (2010); Pettersson et al (2015). New numerical schemes for the SG system can be found in Jin & Ma (2017); Bürger et al (2014); Wan & Karniadakis (2006) and for the convergence analysis of approximate solutions of the SG system see Gottlieb & Xiu (2008); Hu et al (2015); Zhou & Tang (2012). However, in all these works the dimension of the discrete stochastic space is chosen in an ad hoc way, in particular independent of the spatial-temporal resolution.…”
In this article we present an a posteriori error estimator for the spatial-stochastic error of a Galerkin-type discretization of an initial value problem for a random hyperbolic conservation law. For the stochastic discretization we use the Stochastic Galerkin method and for the spatial-temporal discretization of the Stochastic Galerkin system a Runge-Kutta Discontinuous Galerkin method. The estimator is obtained using smooth reconstructions of the discrete solution. Combined with the relative entropy stability framework of Dafermos (2016), this leads to computable error bounds for the space-stochastic discretization error. Moreover, it turns out that the error estimator admits a splitting into one part representing the spatial error, and a remaining term, which can be interpreted as the stochastic error. This decomposition allows us to balance the errors arising from spatial and stochastic discretization. We conclude with some numerical examples confirming the theoretical findings.
“…Further developments of the Multi-Element approach encompass h-and hp-adaptive refinements in the stochastic space ( [33,34,36]) or a multi-resolution discretization using wavelets instead of gPC, cf. [6,24]. Another approach which ensures hyperbolicity of the resulting nonlinear SG system is the intrusive polynomial moment method.…”
Intrusive Uncertainty Quantification methods such as stochastic Galerkin are gaining popularity, whereas the classical stochastic Galerkin approach is not ensured to preserve hyperbolicity of the underlying hyperbolic system. We apply a modification of this method that uses a slope limiter to retain admissible solutions of the system, while providing high-order approximations in the physical and stochastic space. This is done using a spatial discontinuous Galerkin scheme and a Multi-Element stochastic Galerkin ansatz in the random space. We analyze the convergence of the resulting scheme and apply it to the compressible Euler equations with various uncertain initial states in one and two spatial domains with up to three uncertainties. The performance in multiple stochastic dimensions is compared to the non-intrusive Stochastic Collocation method. The numerical results underline the strength of our method, especially if discontinuities are present in the uncertainty of the solution.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.