Matrix-oriented FEM formulation for stationary and time-dependent PDEs on x-normal domains

09/02/2021
by   Massimo Frittelli, et al.
0

When numerical solution of elliptic and parabolic partial differential equations is required to be highly accurate in space, the discrete problem usually takes the form of large-scale and sparse linear systems. In this work, as an alternative, for spatial discretization we provide a Matrix-Oriented formulation of the classical Finite Element Method, called MO-FEM, of arbitrary order k∈ℕ. On structured 2D domains (e.g. squares or rectangles) the discrete problem is then reformulated as a Sylvester matrix equation, that we solve by the reduced approach in the associated spectral space. On a quite general class of domains, namely normal domains, and even on special surfaces, the MO-FEM yields a multiterm Sylvester matrix equation where the additional terms account for the geometric contribution of the domain shape. In particular, we obtain a sequence of these equations after time discretization of parabolic problems by the IMEX Euler method. We apply the matrix-oriented form of the Preconditioned Conjugate Gradient (MO-PCG) method to solve each multiterm Sylvester equation for MO-FEM of degree k=1,…,4 and for the lumped ℙ_1 case. We choose a matrix-oriented preconditioner with a single-term form that captures the spectral properties of the whole multiterm Sylvester operator. For several numerical examples, we show a gain in computational time and memory occupation wrt the classical vector approach solving large sparse linear systems by a direct method or by the vector PCG with same preconditioning. As an application, we show the advantages of the MO-FEM-PCG to approximate Turing patterns with high spatial resolution in a reaction-diffusion PDE system for battery modeling.

READ FULL TEXT

page 9

page 28

research
06/17/2022

An integral equation method for the advection-diffusion equation on time-dependent domains in the plane

Boundary integral methods are attractive for solving homogeneous linear ...
research
05/02/2019

Lifting Vectorial Variational Problems: A Natural Formulation based on Geometric Measure Theory and Discrete Exterior Calculus

Numerous tasks in imaging and vision can be formulated as variational pr...
research
05/17/2021

Full operator preconditioning and the accuracy of solving linear systems

Unless special conditions apply, the attempt to solve ill-conditioned sy...
research
06/20/2021

A PDE approach to centroidal tessellations of domains

We introduce a class of systems of Hamilton-Jacobi equations that charac...
research
06/10/2022

Preconditioned infinite GMRES for parameterized linear systems

We are interested in obtaining approximate solutions to parameterized li...
research
05/07/2020

A Reduced Basis Method For Fractional Diffusion Operators II

We present a novel numerical scheme to approximate the solution map s u(...
research
10/05/2022

Geometric discretization of diffeomorphisms

Many partial differential equations in mathematical physics describe the...

Please sign up or login with your details

Forgot password? Click here to reset