Abstract
Summary Recently introduced non-reflecting boundary conditions are numerically exact: the solution on a given domain is the same as a subset of one on an enlarged domain where boundary reflections do not have time to reach the original domain. In 1D with second- or higher-order finite differences, a recurrence relation based on translation invariance provides the boundary conditions. In 2D or 3D, a recurrence relation was only found for a non-reflecting boundary on one or two opposing sides of the domain and zero Dirichlet or Neumann boundaries elsewhere. Otherwise, corners cause translation invariance to be lost. The proposed workaround restores translation invariance with classic, approximately non-reflecting boundary conditions on the other sides. As a proof of principle, the method is applied to the 2-D constant-density acoustic wave equation, discretized on a rectangular domain with a second-order finite-difference scheme, first-order Enquist-Majda boundary conditions as approximate ones, and numerically exact boundary conditions in the horizontal direction. The method is computationally costly but has the advantage that it can be reused on a sequence of problems as long as the time step and the sound speed values next to the boundary are kept fixed.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.