An official website of the United States government
Abstract The inverse problem for radiative transfer is important in many applications, such as optical tomography and remote sensing. Major challenges include large memory requirements and computational expense, which arise from high-dimensionality and the need for iterations in solving the inverse problem. Here, to alleviate these issues, we propose adaptive-mesh inversion: a goal-orientedhp-adaptive mesh refinement method for solving inverse radiative transfer problems. One novel aspect here is that the two optimizations (one for inversion, and one for mesh adaptivity) are treated simultaneously and blended together. By exploiting the connection between duality-based mesh adaptivity and adjoint-based inversion techniques, we propose a goal-oriented error estimator, which is cheap to compute, and can efficiently guide the mesh-refinement to numerically solve the inverse problem. We use discontinuous Galerkin spectral element methods to discretize the forward and the adjoint problems. Then, based on the goal-oriented error estimator, we propose anhp-adaptive algorithm to refine the meshes. Numerical experiments are presented at the end and show convergence speed-up and reduced memory occupation by the goal-oriented mesh adaptive method.
more »
« less
A unified hp-HDG framework for Friedrichs' PDE systems
This work proposes a unified hp-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs’ type. In particular, we present unified hp-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified hp-formulations and demonstrate the effectiveness of hp-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, and singularities.
more »
« less
- PAR ID:
- 10537016
- Publisher / Repository:
- CAMWA
- Date Published:
- Journal Name:
- Computers & Mathematics with Applications
- Volume:
- 154
- Issue:
- C
- ISSN:
- 0898-1221
- Page Range / eLocation ID:
- 236 to 266
- Subject(s) / Keyword(s):
- Hybridized Discontinuous Galerkin Friedrichs’ system Discontinuous Galerkin Hybridization hp-adaptation
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We present the lowest-order hybridizable discontinuous Galerkin schemes with numerical integration (quadrature), denoted as HDG-P0 for the reaction-diffusion equation and the generalized Stokes equations on conforming simplicial meshes in two- and three-dimensions. Here by lowest order, we mean that the (hybrid) finite element space for the global HDG facet degrees of freedom (DOFs) is the space of piecewise constants on the mesh skeleton. A discontinuous piecewise linear space is used for the approximation of the local primal unknowns. We give the optimal a priori error analysis of the proposed HDG-P0 schemes, which hasn’t appeared in the literature yet for HDG discretizations as far as numerical integration is concerned. Moreover, we propose optimal geometric multigrid preconditioners for the statically condensed HDG-P0 linear systems on conforming simplicial meshes. In both cases, we first establish the equivalence of the statically condensed HDG system with a (slightly modified) nonconforming Crouzeix–Raviart (CR) discretization, where the global (piecewise-constant) HDG finite element space on the mesh skeleton has a natural one-to-one correspondence to the nonconforming CR (piecewise-linear) finite element space that live on the whole mesh. This equivalence then allows us to use the well-established nonconforming geometry multigrid theory to precondition the condensed HDG system. Numerical results in two- and three-dimensions are presented to verify our theoretical findings.more » « less
-
In this article, we present a three-dimensional anisotropic hp-mesh refinement strategy for ultraweak discontinuous Petrov-Galerkin (DPG) formulations with optimal test functions. The refinement strategy utilizes the built-in residual-based error estimator accompanying the DPG discretization. The refinement strategy is a two-step process: (a) use the built-in error estimator to mark and isotropically hp-refine elements of the (coarse) mesh to generate a finer mesh; (b) use the reference solution on the finer mesh to compute optimal ℎ-and 𝑝-refinements of the selected elements in the coarse mesh. The process is repeated with coarse and fine mesh being generated in every adaptation cycle, until a prescribed error tolerance is achieved. We demonstrate the performance of the proposed refinement strategy using several numerical examples on hexahedral meshes.more » « less
-
Abstract We derive a residual-based a posteriori error estimator for the conforming hp -Adaptive Finite Element Method ( hp -AFEM) for the steady state Stokes problem describing the slow motion of an incompressible fluid. This error estimator is obtained by extending the idea of a posteriori error estimation for the classical h -version of AFEM. We also establish the reliability and efficiency of the error estimator. The proofs are based on the well-known Clément-type interpolation operator introduced in [28] in the context of the hp -AFEM. Numerical experiments show the performance of an adaptive hp-FEM algorithm using the proposed a posteriori error estimator.more » « less
-
Abstract In their article “Coupling at a distance HDG and BEM” , Cockburn, Sayas and Solano proposed an iterative coupling of the hybridizable discontinuous Galerkin method (HDG) and the boundary element method (BEM) to solve an exterior Dirichlet problem. The novelty of the numerical scheme consisted of using a computational domain for the HDG discretization whose boundary did not coincide with the coupling interface. In their article, the authors provided extensive numerical evidence for convergence, but the proof of convergence and the error analysis remained elusive at that time. In this article we fill the gap by proving the convergence of a relaxation of the algorithm and providing a priori error estimates for the numerical solution.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.camwa.2023.12.009
