CORE
πΊπ¦Β
Β make metadata, not war
Services
Services overview
Explore all CORE services
Access to raw data
API
Dataset
FastSync
Content discovery
Recommender
Discovery
OAI identifiers
OAI Resolver
Managing content
Dashboard
Bespoke contracts
Consultancy services
Support us
Support us
Membership
Sponsorship
Community governance
Advisory Board
Board of supporters
Research network
About
About us
Our mission
Team
Blog
FAQs
Contact us
Sixth-Order Hybrid FDMs and/or the M-Matrix Property for Elliptic Interface Problems with Mixed Boundary Conditions
Authors
Qiwei Feng
Bin Han
Peter Minev
Publication date
22 June 2023
Publisher
View
on
arXiv
Abstract
In this paper, we develop sixth-order hybrid finite difference methods (FDMs) for the elliptic interface problem
β
β
β
(
a
β
u
)
=
f
-\nabla \cdot( a\nabla u)=f
β
β
β
(
a
β
u
)
=
f
in
Ξ©
\
Ξ
\Omega\backslash \Gamma
Ξ©\Ξ
, where
Ξ
\Gamma
Ξ
is a smooth interface inside
Ξ©
\Omega
Ξ©
. The variable scalar coefficient
a
>
0
a>0
a
>
0
and source
f
f
f
are possibly discontinuous across
Ξ
\Gamma
Ξ
. The hybrid FDMs utilize a 9-point compact stencil at any interior regular point of the grid and a 13-point stencil at irregular points near
Ξ
\Gamma
Ξ
. For interior regular points away from
Ξ
\Gamma
Ξ
, we obtain a sixth-order 9-point compact FDM satisfying the M-matrix property. Consequently, for the elliptic problem without interface (i.e.,
Ξ
\Gamma
Ξ
is empty), our compact FDM satisfies the discrete maximum principle, which guarantees the theoretical sixth-order convergence. We also derive sixth-order compact (4-point for corners and 6-point for edges) FDMs having the M-matrix property at any boundary point subject to (mixed) Dirichlet/Neumann/Robin boundary conditions. For irregular points near
Ξ
\Gamma
Ξ
, we propose fifth-order 13-point FDMs, whose stencil coefficients can be effectively calculated by recursively solving several small linear systems. Theoretically, the proposed high order FDMs use high order (partial) derivatives of the coefficient
a
a
a
, the source term
f
f
f
, the interface curve
Ξ
\Gamma
Ξ
, the two jump functions along
Ξ
\Gamma
Ξ
, and the functions on
β
Ξ©
\partial \Omega
β
Ξ©
. Numerically, we always use function values to approximate all required high order (partial) derivatives in our hybrid FDMs without losing accuracy. Our proposed FDMs are independent of the choice representing
Ξ
\Gamma
Ξ
and are also applicable if the jump conditions on
Ξ
\Gamma
Ξ
only depend on the geometry (e.g., curvature) of the curve
Ξ
\Gamma
Ξ
. Our numerical experiments confirm the sixth-order convergence in the
l
β
l_{\infty}
l
β
β
norm of the proposed hybrid FDMs for the elliptic interface problem
Similar works
Full text
Available Versions
arXiv.org e-Print Archive
See this paper in CORE
Go to the repository landing page
Download from data provider
oai:arXiv.org:2306.13001
Last time updated on 26/06/2023