A matrix-free isogeometric Galerkin method for Karhunen-Loève approximation of random fields using tensor product splines, tensor contraction and interpolation based quadrature

DOI10.1016/J.CMA.2021.113730MaRDI QIDQ2021879zbMATH OpenOpenAlexFDO

Authors Michal L. Mika, Dominik Schillinger, P. Wriggers, R. R. Hiemstra, T. J. R. Hughes

Publication date 27 April 2021

Published in Computer Methods in Applied Mechanics and Engineering (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2011.13861

Fredholm integral eigenvalue problem isogeometric analysis Kronecker products matrix-free solver random fields

Random fields (60G60) Numerical computation using splines (65D07) Finite element, Rayleigh-Ritz and Galerkin methods for boundary value problems involving PDEs (65N30) Probabilistic methods, particle methods, etc. for boundary value problems involving PDEs (65N75) Finite element methods applied to problems in solid mechanics (74S05) Stochastic and other probabilistic methods applied to problems in solid mechanics (74S60)

Abstract: The Karhunen-Lo`eve series expansion (KLE) decomposes a stochastic process into an infinite series of pairwise uncorrelated random variables and pairwise

L^{2}

-orthogonal functions. For any given truncation order of the infinite series the basis is optimal in the sense that the total mean squared error is minimized. The orthogonal basis functions are determined as the solution of an eigenvalue problem corresponding to the homogeneous Fredholm integral equation of the second kind, which is computationally challenging for several reasons. Firstly, a Galerkin discretization requires numerical integration over a

2 d

dimensional domain, where

d

, in this work, denotes the spatial dimension. Secondly, the main system matrix of the discretized weak-form is dense. Consequently, the computational complexity of classical finite element formation and assembly procedures as well as the memory requirements of direct solution techniques become quickly computationally intractable with increasing polynomial degree, number of elements and degrees of freedom. The objective of this work is to significantly reduce several of the computational bottlenecks associated with numerical solution of the KLE. We present a matrix-free solution strategy, which is embarrassingly parallel and scales favorably with problem size and polynomial degree. Our approach is based on (1) an interpolation based quadrature that minimizes the required number of quadrature points; (2) an inexpensive reformulation of the generalized eigenvalue problem into a standard eigenvalue problem; and (3) a matrix-free and parallel matrix-vector product for iterative eigenvalue solvers. Two higher-order three-dimensional benchmarks illustrate exceptional computational performance combined with high accuracy and robustness.

Recommendations

Cites work

Cited in

(24)

Describes a project that uses

Uses Software

This page was built for publication: A matrix-free isogeometric Galerkin method for Karhunen-Loève approximation of random fields using tensor product splines, tensor contraction and interpolation based quadrature

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2021879)