Abstract
In this work, we are concerned with neural network guided goal-oriented a posteriori error estimation and adaptivity using the dual weighted residual method. The primal problem is solved using classical Galerkin finite elements. The adjoint problem is solved in strong form with a feedforward neural network using two or three hidden layers. The main objective of our approach is to explore alternatives for solving the adjoint problem with greater potential of a numerical cost reduction. The proposed algorithm is based on the general goal-oriented error estimation theorem including both linear and nonlinear stationary partial differential equations and goal functionals. Our developments are substantiated with some numerical experiments that include comparisons of neural network computed adjoints and classical finite element solutions of the adjoints. In the programming software, the open-source library deal.II is successfully coupled with LibTorch, the PyTorch C++ application programming interface.Article HighlightsAdjoint approximation with feedforward neural network in dual-weighted residual error estimation.Side-by-side comparisons for accuracy and computational cost with classical finite element computations.Numerical experiments for linear and nonlinear problems yielding excellent effectivity indices.
Highlights
This work is devoted to an innovative solution of the adjoint equation in goal-oriented error estimation with the dual weighted residual (DWR) method [3,4,5]; we refer to [1, 7, 22, 41] for some important early work
This adjoint solution is usually obtained by global higher order finite element method (FEM) solutions or local higher order approximations [5]
Namely the necessity of working with strong formulations, the current study provides useful insights whether at all neural network guided adjoints can be an alternative concept for dual weighted residual error estimation
Summary
This work is devoted to an innovative solution of the adjoint equation in goal-oriented error estimation with the dual weighted residual (DWR) method [3,4,5] (based on former adjoint concepts [19]); we refer to [1, 7, 22, 41] for some important early work. This adjoint solution is usually obtained by global higher order finite element method (FEM) solutions or local higher order approximations [5]. The former is more stable, see e.g. As the adjoint solution is only required to evaluate the a posteriori error estimator, a cheap solution is of interest
Published Version (
Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have