CWENO Interpolation for Non-Oscillatory Stochastic Collocation in Uncertainty Quantification Problems

Alina Chertock,Arsen S. Iskhakov,Anna Iskhakova,Alexander Kurganov
2024-11-20
Abstract:Uncertainty quantification (UQ) in mathematical models is essential for accurately predicting system behavior under variability. This study provides guidance on method selection for reliable UQ across varied functional behaviors in engineering applications. Specifically, we compare several interpolation and approximation methods within a stochastic collocation (SC) framework, namely: generalized polynomial chaos (gPC), B-splines, shape-preserving (SP) splines, and central weighted essentially nonoscillatory (CWENO) interpolation, to reconstruct probability density functions (PDFs) and estimate statistical moments. These methods are assessed for both smooth and discontinuous functions, as well as for the solution of the 1-D Euler and shallow water equations. While gPC and interpolation B-splines perform well with smooth data, they produce oscillations near discontinuities. Approximation B-splines and SP splines, while avoiding oscillations, converge more slowly. In contrast, CWENO interpolation demonstrates high robustness, effectively capturing sharp gradients without oscillations, making it suitable for complex, discontinuous data. Overall, CWENO interpolation emerges as a versatile and effective approach for SC, particularly in handling discontinuities in UQ.
Numerical Analysis
What problem does this paper attempt to address?
This paper attempts to solve the problems of how to accurately reconstruct the probability density function (PDF) and estimate statistical moments in the Uncertainty Quantification (UQ) problem. Specifically, the paper explores the performance of several interpolation and approximation methods in the Stochastic Collocation (SC) framework under different function behaviors (such as smooth and discontinuous functions). These methods include Generalized Polynomial Chaos Expansion (gPC), B - splines, Shape - Preserving Splines (SP splines), and Central - Weighted Essentially Non - Oscillatory (CWENO) interpolation. ### Main Problems and Background 1. **Uncertainty Quantification**: There is inherent uncertainty in many engineering problems, which may come from measurement errors or empirical approximations of input parameters or initial boundary conditions. Quantifying these uncertainties is crucial for conducting sensitivity analysis and improving model accuracy. 2. **Limitations of Existing Methods**: - **Monte Carlo Method**: Although reliable, it has a high computational cost and requires a large number of sample realizations. - **gPC Method**: It performs well for smooth data, but it will produce oscillation phenomena when dealing with discontinuities (such as shocks). - **B - splines**: Problems will occur at discontinuities whether in interpolation or approximation. Interpolating B - splines produce oscillations near discontinuities, while approximating B - splines will over - smooth the data. - **SP Splines**: It has a good local smoothing effect, but its convergence speed is slow, and its accuracy is limited in complex situations. 3. **Requirement for New Methods**: In order to effectively handle complex and discontinuous data, a method that can maintain high accuracy and avoid oscillation is needed. ### Research Objectives The main objective of this research is to find an interpolation method that can effectively handle complex and discontinuous data by comparing different SC methods, so as to achieve reliable UQ. Specifically, the research aims to evaluate the capabilities of these methods in the following aspects: - Reconstructing the probability density functions (PDF) of smooth and discontinuous functions. - Estimating statistical moments (such as mean and variance). - Solving problems in practical applications such as one - dimensional Euler equations and shallow water equations. ### Advantages of CWENO Interpolation The paper particularly emphasizes the advantages of the CWENO interpolation method: - **High Robustness**: It can maintain high accuracy in smooth regions and avoid oscillation near discontinuities. - **Suitable for Complex Data**: It is especially suitable for handling data with discontinuities and complex gradients, providing a stable PDF approximation. In conclusion, through a series of numerical experiments, this paper verifies the effectiveness and superiority of CWENO interpolation in handling complex and discontinuous data, providing new tools and ideas for future UQ research.