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

DOI10.1007/S10915-016-0268-8zbMATH Open1361.65070arXiv1511.01013OpenAlexW2509202927MaRDI QIDQ520216FDOQ520216

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)

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.

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

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

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)

Cites Work

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

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
• Method of lines transpose: high order L-stable \(\mathcal O(N)\) schemes for parabolic equations using successive convolution
• A particle-in-cell method for plasmas with a generalized momentum formulation. II: Enforcing the Lorenz gauge condition
• 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

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)