edoc

H-matrix based second moment analysis for rough random fields and finite element discretizations

Dölz, Jürgen and Harbrecht, Helmut and Peters, Michael. (2017) H-matrix based second moment analysis for rough random fields and finite element discretizations. SIAM Journal on Scientific Computing, 39 (4). B618-B639.

[img]
Preview
PDF - Published Version
3910Kb

Official URL: http://edoc.unibas.ch/55822/

Downloads: Statistics Overview

Abstract

We consider the efficient solution of strongly elliptic partial differential equations with random load based on the finite element method. The solution's two-point correlation can efficiently be approximated by means of an H- matrix, in particular if the correlation length is rather short or the correlation kernel is nonsmooth. Since the inverses of the finite element matrices which correspond to the differential operator under consideration can likewise efficiently be approximated in the H- matrix format, we can solve the correspondent H- matrix equation in essentially linear time by using the H -matrix arithmetic. Numerical experiments for three-dimensional finite element discretizations for several correlation lengths and different smoothness are provided. They validate the presented method and demonstrate that the computation times do not increase for nonsmooth or shortly correlated data.
Faculties and Departments:05 Faculty of Science > Departement Mathematik und Informatik > Mathematik > Computational Mathematics (Harbrecht)
UniBasel Contributors:Harbrecht, Helmut and Dölz, Jürgen and Peters, Michael
Item Type:Article, refereed
Article Subtype:Research Article
Publisher:Society for Industrial and Applied Mathematics
ISSN:1064-8275
e-ISSN:1095-7197
Note:Publication type according to Uni Basel Research Database: Journal article
Language:English
Last Modified:02 Oct 2017 14:07
Deposited On:02 Oct 2017 14:07

Repository Staff Only: item control page