A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion
- CETEC-UFRB, Centro (Brazil)
- UFPR, Departamento de Matemática (Brazil)
We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.
- OSTI ID:
- 22769330
- Journal Information:
- Computational and Applied Mathematics, Vol. 37, Issue 2; Other Information: Copyright (c) 2018 SBMAC - Sociedade Brasileira de Matemática Aplicada e Computacional; Country of input: International Atomic Energy Agency (IAEA); ISSN 0101-8205
- Country of Publication:
- United States
- Language:
- English
Similar Records
Characterization of a random anisotropic conductivity field with Karhunen-Loeve methods
Karhunen–Loève treatment to remove noise and facilitate data analysis in sensing, spectroscopy and other applications