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

DOI10.1016/J.JCP.2017.10.015zbMATH Open1380.65138arXiv1609.05379OpenAlexW2523108702MaRDI QIDQ1701282FDOQ1701282

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)

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.

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

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

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)

Cites Work

• Title not available (Why is that?)
• Numerical analysis of blood flow in the heart
• A non-oscillatory Eulerian approach to interfaces in multimaterial flows (the ghost fluid method)
• The Immersed Interface Method for Elliptic Equations with Discontinuous Coefficients and Singular Sources
• A Fast Iterative Algorithm for Elliptic Interface Problems
• The immersed interface method for acoustic wave equations with discontinuous coefficients
• 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
• 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 new interface method for hyperbolic problems with discontinuous coefficients: One-dimensional acoustic example
• 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
• Numerical treatment of two-dimensional interfaces for acoustic and elastic waves.
• Title not available (Why is that?)
• Edge element formulations of eddy current problems
• On the immersed interface method for solving time-domain Maxwell's equations in materials with curved dielectric interfaces
• High resolution sharp computational methods for elliptic and parabolic problems in complex geometries

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)