Method of lines transpose: an efficient unconditionally stable solver for wave propagation

DOI10.1007/S10915-016-0268-8MaRDI QIDQ520216zbMATH OpenOpenAlexFDO

Authors Andrew Christlieb, Matthew F. Causley, Eric M. Wolf

Publication date 3 April 2017

Published in Journal of Scientific Computing (Search for Journal in Brave)

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

zbMATH Keywords

algorithmnumerical exampleperiodic boundary conditionswave equationinitial value problemalternating direction implicit methodsboundary integral methodsimplicit methodsDirichlet and Neumann boundary conditionsconsistent boundary conditionsmethod of lines transposetransverse method of lines

Mathematics Subject Classification ID

Wave equation (35L05) Method of lines for initial value and initial-boundary value problems involving PDEs (65M20) Boundary element methods for initial value and initial-boundary value problems involving PDEs (65M38)

Abstract: Building upon recent results obtained in [7,8,9], we describe an efficient second order, A-stable scheme for solving the wave equation, based on the method of lines transpose (MOL${}^T$), and the resulting semi-discrete (i.e. continuous in space) boundary value problem. In [7], A-stable schemes of high order were derived, and in [9] a high order, fast $mathcalO(N)$ spatial solver was derived, which is matrix-free and is based on dimensional-splitting. In this work, are interested in building a wave solver, and our main concern is the development of boundary conditions. We demonstrate all desired boundary conditions for a wave solver, including outflow boundary conditions, in 1D and 2D. The scheme works in a logically Cartesian fashion, and the boundary points are embedded into the regular mesh, without incurring stability restrictions, so that boundary conditions are imposed without any reduction in the order of accuracy. We demonstrate how the embedded boundary approach works in the cases of Dirichlet and Neumann boundary conditions. Further, we develop outflow and periodic boundary conditions for the MOL${}^T$ formulation. Our solver is designed to couple with particle codes, and so special attention is also paid to the implementation of point sources, and soft sources which can be used to launch waves into waveguides.

Recommendations

• Method of lines transpose: an implicit solution to the wave equation
• An efficient technique for solving linear conservative wave propagation problem
• scientific article; zbMATH DE number 4045058
• A linearized implicit finite-difference method for solving the equal width wave equation
• High-Accuracy Finite-Difference Schemes for Linear Wave Propagation
• A method of linearizing the inverse problem for the wave equation
• An effective numerical method for linear mode conversion problems
• An algorithm for direct simulation of linear wave propagation in irregular regions
• Method of lines transpose: high order L-stable \(\mathcal O(N)\) schemes for parabolic equations using successive convolution
• An efficient boundary element method for two-dimensional transient wave propagation problems

Cites work

• A Cartesian treecode for screened Coulomb interactions
• A Generalized Fast Multipole Method for Nonoscillatory Kernels
• A fast adaptive multipole algorithm in three dimensions
• A fast algorithm for particle simulations
• A particle method and adaptive treecode for vortex sheet motion in three-dimensional flow
• A particle-in-cell method for the simulation of plasmas based on an unconditionally stable field solver
• A treecode algorithm for simulating electron dynamics in a Penning-Malmberg trap
• An adaptive fast solver for the modified Helmholtz equation in two dimensions
• An integral evolution formula for the wave equation
• Difference Approximations of the Neumann Problem for the Second Order Wave Equation
• Divergence preservation in the ADI algorithms for electromagnetics
• High-order unconditionally stable FC-AD solvers for general smooth domains. I: Basic elements
• High-order unconditionally stable FC-AD solvers for general smooth domains. II: Elliptic, parabolic and hyperbolic PDEs; theoretical considerations
• Higher order A-stable schemes for the wave equation using a successive convolution approach
• Low order approximation of the spherical nonreflecting boundary kernel for the wave equation
• Method of lines transpose: an implicit solution to the wave equation
• Method of lines transpose: high order L-stable \(\mathcal O(N)\) schemes for parabolic equations using successive convolution
• Nonreflecting boundary conditions for the time-dependent wave equation
• Stability and accuracy of time-extrapolated ADI-FDTD methods for solving wave equations

Cited in

(6)

• Higher order A-stable schemes for the wave equation using a successive convolution approach
• Method of lines transpose: energy gradient flows using direct operator inversion for phase-field models
• A high order finite difference weighted essentially nonoscillatory schemes with a kernel-based constrained transport method for ideal magnetohydrodynamics
• On the operator splitting and integral equation preconditioned deferred correction methods for the ``good Boussinesq equation
• A particle-in-cell method for plasmas with a generalized momentum formulation. II: Enforcing the Lorenz gauge condition
• Method of lines transpose: high order L-stable \(\mathcal O(N)\) schemes for parabolic equations using successive convolution

This page was built for publication: Method of lines transpose: an efficient unconditionally stable solver for wave propagation

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