Sparse finite element approximation of high-dimensional transport-dominated diffusion problems

Issue		ESAIM: M2AN Volume 42, Number 5, September-October 2008


Page(s)		777 - 819
DOI		10.1051/m2an:2008027
Published online		30 July 2008

ESAIM: M2AN 42 (2008) 777-819
DOI: 10.1051/m2an:2008027

Sparse finite element approximation of high-dimensional transport-dominated diffusion problems

Christoph Schwab¹, Endre Süli² and Radu Alexandru Todor¹

¹ Seminar für Angewandte Mathematik, Eidgenössische Technische Hochschule, 8092 Zürich, Switzerland. schwab@sam.math.ethz.ch; todor@math.ethz.ch
² University of Oxford, Computing Laboratory, Wolfson Building, Parks Road, Oxford OX1 3QD, UK. endre.suli@comlab.ox.ac.uk

Received March 7, 2007. Received February 11, 2008. Published online July 30, 2008.

Abstract
We develop the analysis of stabilized sparse tensor-product finite element methods for high-dimensional, non-self-adjoint and possibly degenerate second-order partial differential equations of the form $-a:\nabla\nabla u + b \cdot \nabla u + cu = f(x)$ , $x \in \Omega = (0,1)^d \subset \mathbb{R} ^d$ , where $a \in \mathbb{R} ^{d\times d}$ is a symmetric positive semidefinite matrix, using piecewise polynomials of degree $p\geq 1$ . Our convergence analysis is based on new high-dimensional approximation results in sparse tensor-product spaces. We show that the error between the analytical solution u and its stabilized sparse finite element approximation u_h on a partition of $\Omega$ of mesh size h=h_L = 2^-L satisfies the following bound in the streamline-diffusion norm $\vert\vert\vert\cdot\vert\vert\vert _{\rm SD}$ , provided u belongs to the space $\mathcal{H}^{k+1}(\Omega)$ of functions with square-integrable mixed (k+1)st derivatives:

$\begin{displaymath}\vert\vert\vert u-u_h\vert\vert\vert _{\rm SD}\leq C_{p,t} d^... ...vert u\vert _{\mathcal{H}^{t+1}(\Omega)}, \qquad \qquad \qquad \end{displaymath}$

where $\kappa_i=\kappa_i(p,t,L)$ , i=0,1, and $1 \leq t \leq \min(k,p)$ . We show, under various mild conditions relating L to p, L to d, or p to d, that in the case of elliptic transport-dominated diffusion problems $\kappa_0, \kappa_1 \in (0,1)$ , and hence for $p \geq 2$ the `error constant' $C_{p,t} d^2 \max\{(2-p)_+,\kappa_0^{d-1},\kappa_1^d\}$ exhibits exponential decay as $d \rightarrow \infty$ ; in the case of a general symmetric positive semidefinite matrix a, the error constant is shown to grow no faster than $\mathcal{O}(d^2)$ . In any case, in the absence of assumptions that relate L, p and d, the error $\vert\vert\vert u - u_h\vert\vert\vert _{\rm SD}$ is still bounded by $\kappa_\ast^{d-1} \vert\log_2 h_L\vert^{d-1}\mathcal{O}(\vert\sqrt{a}\vert h_L^t + \vert b\vert^{\frac{1}{2}} h_L^{t+\frac{1}{2}} + c^{\frac{1}{2}} h_L^{t+1})$ , where $\kappa_\ast \in (0,1)$ for all $L, p, d \geq 2$ .

Mathematics Subject Classification. 65N30.

Key words: High-dimensional Fokker-Planck equations, partial differential equations with nonnegative characteristic form, sparse finite element method.

© EDP Sciences, SMAI 2008

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.