Imperial College London
Portrait Picture

Radu Cimpeanu

Faculty of Natural Sciences, Department of Mathematics

Honorary Lecturer
 
 
 
//

Contact

 

radu.cimpeanu11 CV

 
 
//

Location

 

742Huxley BuildingSouth Kensington Campus

//

Summary

 

Publications

Citation

BibTex format

@article{Cimpeanu:2015:10.1016/j.jcp.2015.05.006,
author = {Cimpeanu, R and Martinsson, A and Heil, M},
doi = {10.1016/j.jcp.2015.05.006},
journal = {Journal of Computational Physics},
pages = {329--347},
title = {A parameter-free perfectly matched layer formulation for the finite-element-based solution of the Helmholtz equation},
url = {http://dx.doi.org/10.1016/j.jcp.2015.05.006},
volume = {296},
year = {2015}
}
Download

RIS format (EndNote, RefMan)

TY  - JOUR
AB  - This paper presents a parameter-free perfectly matched layer (PML) method for the finite-element-based solution of the Helmholtz equation. We employ one of Bermúdez et al.'s unbounded absorbing functions for the complex coordinate mapping underlying the PML. With this choice, the only free parameter that controls the accuracy of the numerical solution for a fixed numerical cost (characterised by the number of elements in the bulk and the PML regions) is the thickness of the perfectly matched layer, δPML. We show that, for the case of planar waves, the absorbing function performs best for PMLs whose thickness is much smaller than the wavelength. We then perform extensive numerical experiments to explore its performance for non-planar waves, considering domain shapes with smooth and polygonal boundaries, different solution types (smooth and singular), and a wide range of wavenumbers, k  , to identify an optimal range for the normalised PML thickness, kδPML, such that, within this range, the error introduced by the presence of the PML is consistently small and insensitive to change. This implies that if the PML thickness is chosen from within this range no further PML optimisation is required, i.e. the method is parameter-free. We characterise the dependence of the error on the discretisation parameters and establish the conditions under which the convergence of the solution under mesh refinement is controlled exclusively by the discretisation of the bulk mesh.
AU  - Cimpeanu,R
AU  - Martinsson,A
AU  - Heil,M
DO  - 10.1016/j.jcp.2015.05.006
EP  - 347
PY  - 2015///
SN  - 1090-2716
SP  - 329
TI  - A parameter-free perfectly matched layer formulation for the finite-element-based solution of the Helmholtz equation
T2  - Journal of Computational Physics
UR  - http://dx.doi.org/10.1016/j.jcp.2015.05.006
UR  - http://hdl.handle.net/10044/1/40633
VL  - 296
ER  -
Download