The original problem in an unbounded domain is reduced to a linear parametric eigenvalue problem in a circle, which is convenient for numerical solution. The examination of the solvability of this problem is based on the spectral theory of compact self-adjoint operators. The existence of surface guided waves is proved, and properties of the dispersion curves are investigated. An algorithm for the numerical solution of the problem based on the flnite element method is proposed. The convergence of the numerical method is proved. Numerical results are discussed. 1. INTRODUCTION Optical waveguides are dielectric cylindrical structures that can conduct electromagnetic energy in the visible and infrared parts of the spectrum. The waveguides used in optical communication are ∞exible flbers made of transparent dielectrics. The cross section of a waveguide usually consists of three regions: the central region (core) is surrounded by a cladding which, in turn, is surrounded by a protective coating. The dielectric permittivity of the core can be constant or can vary over the cross section; the dielectric permittivity of the cladding is usually positive constant (denote it by 1). The coating is optically isolated from the core; for this reason, it is usually neglected in mathematical models, and it is assumed that the cladding is unbounded from the outside. We use the classical model (see (1)), in which the waveguide is assumed to be unbounded and linearly isotropic. A mathematical analysis of surface waves based on the theory of unbounded self- adjoint operators can be found in (2). In that paper, the original problem is considered as a problem of the form A(fl)H = k 2 H with respect to the spectral parameter k 2 , and the dependence k = k(fl) is studied (H is the magnetic vector amplitude, k is the wavenumber, fl is the propagation constant). In (3), a similar technique is used to extend the results obtained in (2) to the case of waveguides with a variable magnetic permeability. The results obtained in (2,3) give a complete understanding of the qualitative properties of the spectrum of surface guided waves; however, in order to calculate the spectral characteristics of waveguides, numerical methods are needed (see survey (4)). The formulations of the problems used in (2,3) are not quite convenient for obtaining numerical solutions. This is due to two speciflc features of those statements. 1. The problems are formulated for the entire plane R 2 . For a numerical solution, special measures must be taken to restrict the integration domain and to formulate additional boundary conditions. 2. Spectral problems (except for a point spectrum) have a continuous part of the spectrum. Although the location of this part is known exactly, a numerical solution requires that false ap- proximate solutions be detected and discarded. Statements of problems suggested in (5,6) are free of those drawbacks. In those papers, exact nonlocal boundary conditions (see (7,8)) are used to reduce the problems that were originally formulated for the entire planeR 2 to equivalent problems in a circle. In (5,6) the spectral problems are formulated in a circle › which includes waveguide's cross-section domain ›i (see Fig. 1); these problems have no continuous spectrum. Moreover, their spectrum is identical to the point part of the spectrum of the original problem. These statements are convenient for the flnite element method. The cost of this advantage is that the spectral parameter appears in the equation in a nonlinear fashion; more precisely, the problems have the form A(fl;‚)H = ‚H, where A is a compact self-adjoint operator. The solution of such problems requires the use of special iterative methods. In this paper, we use a new formulation of the problem proposed in (9). The original problem by exact nonlocal boundary conditions method is reduced to an equivalent linear self-adjoint eigenvalue problem A(p)H = fl 2 B(p)H in the circle ›. Here, the parameter p is the transverse wave number p =
Read full abstract