Phase reduction models for improving the accuracy of the finite element solution of time-harmonic scattering problems I: General approach and low-order models

Authors:
Xavier Antoine;Christophe Geuzaine
Affiliations:
Institut Elie Cartan Nancy (IECN), Nancy-Université, CNRS, INRIA Corida Team, Boulevard des Aiguillettes, B.P. 239 F-54506, Vandoeuvre-lès-Nancy Cedex, France and Institut National Polyt ...;University of Liège, Department of Electrical Engineering and Computer Science, Montefiore Institute, Sart Tilman, Bldg. B28, Parking P32, B-4000, Liège, Belgium
Venue:
Journal of Computational Physics
Year:
2009

Citing 5
Cited 2

A perfectly matched layer for the absorption of electromagnetic waves

Journal of Computational Physics
Finite Element Solution of the Helmholtz Equation with High Wave Number Part II: The h-p Version of the FEM

SIAM Journal on Numerical Analysis
Is the Pollution Effect of the FEM Avoidable for the Helmholtz Equation Considering High Wave Numbers?

SIAM Review
A hybrid numerical asymptotic method for scattering problems

Journal of Computational Physics
Iterative Methods for Sparse Linear Systems

Iterative Methods for Sparse Linear Systems

RETRACTED: Time domain computation and visualization of shock induced sound fields for a doubly fluid-loaded hollow cylinder

Computers and Structures
A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation

Journal of Computational Physics

Quantified Score

Hi-index	31.45

Visualization

Abstract

This paper introduces a new formulation of high frequency time-harmonic scattering problems in view of a numerical finite element solution. It is well-known that pollution error causes inaccuracies in the finite element solution of short-wave problems. To partially avoid this precision problem, the strategy proposed here consists in firstly numerically computing at a low cost an approximate phase of the exact solution through asymptotic propagative models. Secondly, using this approximate phase, a slowly varying unknown envelope is introduced and is computed using coarser mesh grids. The global procedure is called Phase Reduction. In this first paper, the general theoretical procedure is developed and low-order propagative models are numerically investigated in detail. Improved solutions based on higher order models are discussed showing the potential of the method for further developments.