A Saddle Point Numerical Method for Helmholtz Equations

Abstract

In a previous work, the author and D.C. Dobson proposed a numerical method for solving the complex Helmholtz equation based on the minimization variational principles developed by Milton, Seppecher, and Bouchitte. This method results in a system of equations with a symmetric positive definite coefficient matrix, but at the same time requires solving simultaneously for the solution and its gradient. Herein is presented a method based on the saddle point variational principles of Milton, Seppecher, and Bouchitte, which produces symmetric positive definite systems of equations, but eliminates the necessity of solving for the gradient of the solution. The result is a method for a wide class of Helmholtz problems based completely on the Conjugate Gradient algorithm.