Spherical Harmonic Solution of the Robin Problem for the Helmholtz Equation in a Supershaped Shell

The Robin problem for the Helmholtz equation in normal-polar shells is addressed by using a suitable spherical harmonic expansion technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by a generalized version of the so-called “superformula” introduced by Gielis. A dedicated numerical procedure based on the computer algebra system Mathematica is developed in order to validate the proposed methodology. In this way, highly accurate approximations of the solution, featuring properties similar to the classical ones, are obtained.


Introduction
Many problems of mathematical physics and electromagnetics are related to the Laplacian differential operator.Among them, it is worth mentioning those relevant to the Laplace and Helmholtz equations.However, most of the mentioned differential problems can be solved in explicit way only in canonical domains with special symmetries, such as intervals, cylinders or spheres [1].The solution in more general domains can be obtained by using the Riemann theorem on conformal mappings and the relevant invariance of the Laplacian [2].However, it is clear that conformal mapping techniques can not be used in the three-dimensional case where approaches based on suitable spatial discretization procedures, such as such as finite-difference or finite-element methods, are usually adopted [3].
Different techniques have been proposed in order to solve the mentioned class of differential problems both from a theoretical and numerical point of view (e.g., representing the solution by means of boundary layer techniques [4], solving the corresponding boundary integral equation by iterative methods [5], approximating the relevant Green function by means of the least squares fitting technique [6], solving the linear system relevant to an elliptic partial differential equation by means of relaxation methods [7]).However, none of the contribu-tions already available in the scientific literature deals with the classical Fourier projection method [8] which has been extended in recent papers [9][10][11][12][13][14][15][16] in order to address boundary-value problems (BVPs) in simply connected starlike domains whose boundaries may be regarded as an anisotropically stretched unit circle or sphere centered at the origin.
In this contribution, a suitable methodology, based on the theory of spherical harmonics [17], has been developed in order to compute the solution of the Robin problem for the Helmholtz equation in normal-polar shell-like domains.In particular, the boundaries of the considered domains are supposed to be defined by a generalized version of the so-called Gielis formula (also known as "superformula") [18].Regular functions are assumed to describe the boundary values, but the proposed approach can be easily generalized in case of weakened hypotheses.
In order to verify and validate the developed technique, a suitable numerical procedure based on the computer algebra system Mathematica © has been adopted.By using such procedure, a point-wise convergence of the spherical harmonic series representation of the solution has been observed in the regular points of the boundaries, with Gibbs-like phenomena potentially occurring in the quasi-cusped points.The obtained numerical results are in good agreement with theoretical findings by Carleson [19].

The Laplacian in Stretched Spherical Coordinates
Let us introduce in the real space the usual spherical coordinate system: sin cos , sin sin , cos , and the polar equations: relevant to the boundaries of the supershaped shell which is described by the following chain of inequalities:  are assumed to be piece-wise functions satisfying the condition: In this way, upon introducing the stretched radius  such that: with , the considered shell-like domain can be readily obtained by assuming Remark: Note that, in the stretched coordinate system , ,    the original domain is transformed into the spherical shell of radii a and b, respectively.Hence, in this system one can use classical techniques to solve the Helmholtz equation, including the eigenfunction method [1].


Let us consider a piece-wise function ν (x, y, z) the Laplace operator,  , in spherical coordinates: In the considered stretched coordinate system,  can be represented by setting: In this way, by denoting as R  for the sake of shortness, one can readily find: , and where the subscripts denote the partial differentiation with respect to the polar angles  and  .Substituting Equations ( 8)-( 12) into Equation ( 6) finally yields Equation (13). 1 sin As it can be easily noticed, upon setting , the classical expression of the Laplacian in spherical coordinates is recovered.

 


The Robin Problem for the Helmholtz Equation
Let us consider the interior Robin problem for the Helmholtz equation in a starlike shell , whose boundaries are described by the polar equations where denotes the propagation constant, are the outward-pointing normal unit vectors to the domain boundaries , respectively, and Under the mentioned assumptions, one can prove the following theorem. Theorem.Let: , sin cos , , sin sin , , cos , cos cos sin , where: such that the series expansion (20) holds true.
In Equation (20 denotes the spherical Hankel function of kind and order n. introduce the terms reported in Equation (21), where: and can be determined by solving the infinite linear system: , , Copyright © 2013 SciRes.AM D. CARATELLI ET AL. 266 1,1, , , , p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p q n m p , where:    introduced in the , , x y z space, the considered domain turns into the spherical shell of radii a and b, one can readily adopt the usual eigenfunction method [1] in combination with the separation of variables (with respect to  , r  and  ).As a consequence, elementary solutions of the problem can be searched in the form: Substituting into the Helmholtz equation, one easily finds that the functions     , P    , and     must satisfy the ordinary differential equations: respectively.The parameters  and  are separation constants, whose choice is governed by the physical requirement that at any fixed point in space the scalar field  , , u r    must be single-valued.So, by setting: one can easily find: In this way, it is readily shown that   r  satisfies: and, hence, is a cylinder function of half order that, without loss of generality, can be expressed as a linear combination of ordinary Hankel functions of first and second kind must, so that: .
Enforcing the Robin boundary condition yields: where: and: Hence, combining equations above and using a classical harmonic projection method, the Equations ( 21)-(26) follow after some algebraic manipulations.
It is worth noting that the derived expressions still hold under the assumption that   , R    are piecewise continuous functions, and the boundary values are described by square integrable, not necessarily continuous, functions, so that the relevant spherical harmonic coefficients , , , in Equation ( 19) are finite quantities.

Numerical Procedure
In the following numerical examples, let us assume, for the boundaries of the considered annulus, general polar equations of the type: which provides an extension, to the three-dimensional case, of the "superformula" introduced by Gielis in [18] .Very different characteristic geometries, including ellipsoids, ovaloids, and Lamé-type domains (also called "superellipsoids") can be obtained by assuming suitable values of the parameters , , , , , , , , ,  in (42).It is to emphasized that almost all three-dimensional normal-polar shell-like domains can be described, or closely approximated, by the considered formula.
In order to assess the performance of the proposed methodology in terms of numerical accuracy and convergence rate, the relative boundary error has been evaluated according to Equation (43), where  is the usual norm, and denotes the partial sum of order N relevant to the spherical harmonic expansion representing the solution of the boundary-value problem for the Helmholtz equation (see Equation (44)).
Remark: It is to be noted that, where the boundary values exhibit a rapidly oscillating behavior, the number N of terms in the spherical harmonic expansion approximating the solution of the problem should be increased accordingly in order to achieve the desired numerical accuracy.

Example
Upon assuming in (42) 2, 4,  as been observed that norm of the di nd 2 L fference between the exact solution a the relevant approximation is generally negligible.Point-wise convergence seems to be verified in the considered domains, with the only exception of a set of measure zero consisting of quasi-cusped points.In the neighbourhood of these points, oscillations of the truncated order solution, recalling the classical Gibbs phenomenon, usually take place., , , cos sin .

Conclusion
A harmonic projection meth ord m, obin for expressions of the solution of the con-Ps can be derived by using classical adrature rules, so overcoming the need for cum erh as finite-difference or finite-element methods.The proposed approach has been lidated by dica puter-aided algebra tool Mathematica © .A point-wise convergence of the expan- oints, Gibbs-like oscillations may occur.The computed results are fo to be in goo agreement with the theoretical findings on Fourier series.

m
being the usual Neumann's symbol and    m n P  the associated Legendre function of the first kind and orders n, m.Then, the boundary-value problem (14) for the Helmholtz equation admits a classical solution    denote arbitrary constants.In order to identify the radial function   P  introduced in (27), it is convenient to set: coefficients in the Robin condition, the relative boundary error N e as function of the number N of terms in the truncated ries expansion (44) exhibits the behavior shown in Figure1.As it appears from Figure2, the selection of the expan-

Figure 1 .
Figure 1.Relative boundary error e N as function of the order N of the truncated spherical harmonic expansion representing the solution of the considered Robin problem for the Helmholtz equation in the supershaped shell described by the generalized Gielis formula with parameters  , , , , ,

Figure 2 .
Figure 2. Boundary behavior along  (a) and  (b) of the partial sum Figure 3. Magnitude of the coefficients and , , n m A  the shell  turns to feature a cuboidal shap