A correction function method for the wave equation with interface jump conditions

DOI10.1016/J.JCP.2017.10.015MaRDI QIDQ1701282zbMATH OpenOpenAlexFDO

Authors David S. Abraham, Alexandre Noll Marques, Jean-Christophe Nave

Publication date 22 February 2018

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1609.05379

zbMATH Keywords

Maxwell's equationswave equationhigh-order schemeimmersed methodcorrection function methodinterface jump

Mathematics Subject Classification ID

Finite difference methods applied to problems in optics and electromagnetic theory (78M20) Wave equation (35L05) PDEs with low regular coefficients and/or low regular data (35R05) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06)

Abstract: In this paper a novel method to solve the constant coefficient wave equation, subject to interface jump conditions, is presented. In general, such problems pose issues for standard finite difference solvers, as the inherent discontinuity in the solution results in erroneous derivative information wherever the stencils straddle the given interface. Here, however, the recently proposed Correction Function Method (CFM) is used, in which correction terms are computed from the interface conditions, and added to affected nodes to compensate for the discontinuity. In contrast to existing methods, these corrections are not simply defined at affected nodes, but rather generalized to a continuous function within a small region surrounding the interface. As a result, the correction function may be defined in terms of its own governing partial differential equation (PDE) which may be solved, in principle, to arbitrary order of accuracy. The resulting scheme is not only arbitrarily high order, but also robust, having already seen application to Poisson problems and the heat equation. By extending the CFM to this new class of PDEs, the treatment of wave interface discontinuities in homogeneous media becomes possible. This allows, for example, for the straightforward treatment of infinitesimal source terms and sharp boundaries, free of staircasing errors. Additionally, new modifications to the CFM are derived, allowing compatibility with explicit multi-step methods, such as Runge-Kutta (RK4), without a reduction in accuracy. These results are then verified through numerous numerical experiments in one and two spatial dimensions.

Recommendations

• A correction function method for Poisson problems with interface jump conditions
• A numerical method for solving the variable coefficient wave equation with interface jump conditions
• Imposing jump conditions on nonconforming interfaces for the correction function method: a least squares approach
• An Embedded Boundary Method for the Wave Equation with Discontinuous Coefficients
• High order immersed interface method for acoustic wave equation with discontinuous coefficients

Cites work

• scientific article; zbMATH DE number 1467363 (Why is no real title available?)
• scientific article; zbMATH DE number 3186688 (Why is no real title available?)
• A Fast Iterative Algorithm for Elliptic Interface Problems
• A correction function method for Poisson problems with interface jump conditions
• A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem
• A level set algorithm for tracking discontinuities in hyperbolic conservation laws. II: Systems of equations
• A level-set algorithm for tracking discontinuities in hyperbolic conservation laws. I: Scalar equations
• A new interface method for hyperbolic problems with discontinuous coefficients: One-dimensional acoustic example
• A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method)
• Edge element formulations of eddy current problems
• High resolution sharp computational methods for elliptic and parabolic problems in complex geometries
• Numerical analysis of blood flow in the heart
• Numerical treatment of two-dimensional interfaces for acoustic and elastic waves.
• On the immersed interface method for solving time-domain Maxwell's equations in materials with curved dielectric interfaces
• Stability and Instability in the Computation of Flows with Moving Immersed Boundaries: A Comparison of Three Methods
• The Explicit-Jump Immersed Interface Method: Finite Difference Methods for PDEs with Piecewise Smooth Solutions
• The Immersed Interface Method for Elliptic Equations with Discontinuous Coefficients and Singular Sources
• The immersed interface method for acoustic wave equations with discontinuous coefficients

Cited in

(8)

• Optimal order error estimates for discontinuous Galerkin methods for the wave equation
• High-order FDTD schemes for Maxwell's interface problems with discontinuous coefficients and complex interfaces based on the correction function method
• High-order grid-characteristic method for systems of hyperbolic equations with piecewise constant coefficients
• FDTD schemes for Maxwell's equations with embedded perfect electric conductors based on the correction function method
• A numerical method for solving the variable coefficient wave equation with interface jump conditions
• Imposing jump conditions on nonconforming interfaces for the correction function method: a least squares approach
• Treatment of complex interfaces for Maxwell's equations with continuous coefficients using the correction function method
• An immersed interface method for acoustic wave equations with discontinuous coefficients in complex geometries

This page was built for publication: A correction function method for the wave equation with interface jump conditions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1701282)