Abstract:A finite element method for approximating the solution of a mathematical model for the response of a penetrable, bounded object (obstacle) to the excitation by an external electromagnetic field is presented and investigated. The model consists of a nonlinear Helmholtz equation that is reduced to a spherical domain. The (exemplary) finite element method is formed by Courant-type elements with curved facets at the boundary of the spherical computational domain. This method is examined for its well-posedness, in particular the validity of a discrete inf-sup condition of the modified sesquilinear form uniformly with respect to both the truncation and the mesh parameters is shown. Under suitable assumptions to the nonlinearities, a quasi-optimal error estimate is obtained. Finally, the satisfiability of the approximation property of the finite element space required for the solvability of a class of adjoint linear problems is discussed.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to mathematically simulate the response of a penetrable, finite object (obstacle) to the excitation of an external electromagnetic field. Specifically, the paper focuses on how to approximately solve a nonlinear Helmholtz equation by the finite element method (FEM), which describes the scattering and propagation of electromagnetic waves in a nonlinear polarizable medium. To make the problem numerically solvable, the original problem is transformed into a boundary - value problem. A non - local Dirichlet - to - Neumann (DtN) operator is used to deal with the problem of the infinite domain and confine it within a spherical computational domain.
The main contributions of the paper are as follows:
1. **Establishment of a mathematical model**: A mathematical model for describing the scattering and propagation of electromagnetic waves in a nonlinear medium is proposed, which is based on a nonlinear Helmholtz equation.
2. **Transformation of the problem**: By introducing the DtN operator, the original full - space transmission problem is transformed into an equivalent boundary - value problem, so that it can be numerically solved within a finite spherical domain.
3. **Application of the finite element method**: A finite element method is designed to approximately solve the transformed nonlinear problem. In particular, Courant - type elements with curved edges are used to handle the boundary of the spherical computational domain.
4. **Theoretical analysis**: The well - posedness of the discrete problem is proved. In particular, it is verified that the truncated bilinear form satisfies the parameter - uniform discrete inf - sup condition. Moreover, under appropriate assumptions, a quasi - optimal error estimate is obtained.
5. **Error estimate**: The approximation properties of the finite element space for the solutions of a class of adjoint linear problems are discussed, and the error estimate of the finite element solution is given. Combined with the truncation error estimate obtained from the paper [Ang23], the total error estimate is obtained.
Overall, this paper aims to provide an effective numerical method to solve the scattering and propagation problems of electromagnetic waves in nonlinear media, especially in applications such as optical switches or frequency multipliers.