An MP-DWR method for h-adaptive finite element methods

09/13/2022
by   Chengyu Liu, et al.
0

In a dual weighted residual based adaptive finite element method for solving partial differential equations, a new finite element space needs to be built for solving the dual problem, due to the Galerkin orthogonality. Two popular approaches in the literature are h-refinement and p-refinement for the purpose, respectively, which would cause nonignorable requirements on computational resources. In this paper, a novel approach is proposed for implementing the dual weighted residual method through a multiple precision technique, i.e., the primal and dual problems are solved in two finite element spaces, respectively, built with the same mesh, degrees of freedom, and basis functions, yet different precisions. The feasibility of such an approach is discussed in detail. Besides effectively avoiding the issue introduced by the Galerkin orthogonality, remarkable benefits can also be expected for the efficiency from this new approach, since i). operations on reorganizing mesh grids and/or degrees of freedom in building a finite element space in h- and p-refinement methods can be avoided, and ii). both CPU time and storage needed for solving the derived system of linear equations can be significantly reduced compared with h- and p-refinement methods. In coding the algorithm with a library AFEPack, no extra module is needed compared with the program for solving the primal problem with the employment of a C++ feature template, based on which the floating point precision is used as a parameter of the template. A number of numerical experiments show the effectiveness and improvement of the efficiency of our method. Further discussion towards the limitation of the method is also delivered.

READ FULL TEXT
research
06/25/2020

A simplified primal-dual weak Galerkin finite element method for Fokker-Planck type equations

A simplified primal-dual weak Galerkin (S-PDWG) finite element method is...
research
01/19/2020

New Primal-Dual Weak Galerkin Finite Element Methods for Convection-Diffusion Problems

This article devises a new primal-dual weak Galerkin finite element meth...
research
04/11/2021

A spatial local method for solving 2D and 3D advection-diffusion equations

In this study, an implicit-explicit local differential transform method ...
research
02/24/2021

Neural network guided adjoint computations in dual weighted residual error estimation

In this work, we are concerned with neural network guided goal-oriented ...
research
07/17/2023

Goal-Adaptive Meshing of Isogeometric Kirchhoff-Love Shells

Mesh adaptivity is a technique to provide detail in numerical solutions ...
research
08/14/2023

Towards the efficient calculation of quantity of interest from steady Euler equations II: a CNNs-based automatic implementation

In <cit.>, a dual-consistent dual-weighted residual-based h-adaptive met...
research
04/09/2018

Poly-Spline Finite Element Method

We introduce an integrated meshing and finite element method pipeline en...

Please sign up or login with your details

Forgot password? Click here to reset