We propose a high order numerical method for the first order Maxwell equations in the frequency domain, defined in media with arbitrary complex shape. Our approach is based on the combination of the least squares approach with the spectral element method. The former frees the solution from spurious modes. that can be found sometimes in classical finite element simulations. Many examples of such non-physical solutions exist in literature, and elimination of these spurious effects is a subject of great interest. Spectral elements are a numerical technique for solving partial differential equations which can be regarded as an extension of finite elements: they merge the flexibility of finite elements in dealing with complex geometries, and the better accuracy of spectral methods. Convergence to exact solution is improved by increasing (at run time) the polynomial degree, with no changes on the computational grid: this provides a significant advantage in respect to low order finite element, which necessarily have to resort to grid refinement. In the authors opinion this approach can be successfully used for the treatment of large scale electromagnetic problems or, alternatively, for applications where higher precision is required. We present a few numerical experiments which prove the capability of the method in object. Copyright (C) 2004 John Wiley Sons, Ltd.
Least squares spectral element method for 2D Maxwell equations in the frequency domain
MAZZARELLA, GIUSEPPE;
2004-01-01
Abstract
We propose a high order numerical method for the first order Maxwell equations in the frequency domain, defined in media with arbitrary complex shape. Our approach is based on the combination of the least squares approach with the spectral element method. The former frees the solution from spurious modes. that can be found sometimes in classical finite element simulations. Many examples of such non-physical solutions exist in literature, and elimination of these spurious effects is a subject of great interest. Spectral elements are a numerical technique for solving partial differential equations which can be regarded as an extension of finite elements: they merge the flexibility of finite elements in dealing with complex geometries, and the better accuracy of spectral methods. Convergence to exact solution is improved by increasing (at run time) the polynomial degree, with no changes on the computational grid: this provides a significant advantage in respect to low order finite element, which necessarily have to resort to grid refinement. In the authors opinion this approach can be successfully used for the treatment of large scale electromagnetic problems or, alternatively, for applications where higher precision is required. We present a few numerical experiments which prove the capability of the method in object. Copyright (C) 2004 John Wiley Sons, Ltd.
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.
