Higher-Order Finite Element Methods for the Nonlinear Helmholtz Equation

In this work, we analyze the finite element method with arbitrary but fixed polynomial degree for the nonlinear Helmholtz equation with impedance boundary conditions. We show well-posedness and error estimates of the finite element solution under a resolution condition between the wave number k, the mesh size h and the polynomial degree p of the form “k(kh)p+1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k(kh)^{p+1}$$\end{document} sufficiently small” and a so-called smallness of the data assumption. For the latter, we prove that the logarithmic dependence in h from the case p=1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p=1$$\end{document} in Wu and Zou (SIAM J Numer Anal 56(3):1338–1359, 2018) can be removed for p≥2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p\ge 2$$\end{document}. We show convergence of two different fixed-point iteration schemes. Numerical experiments illustrate our theoretical results and compare the robustness of the iteration schemes with respect to the size of the nonlinearity and the right-hand side data.