Abstract
The Dirichlet problem for the Laplace equation in normalpolar annuli is addressed by using a suitable Fourierlike technique. Attention is in particular focused on the wide class of domains whose boundaries are defined by the socalled ‘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 [1]. In recent papers [29], the classical Fourier projection method [10,11] for solving boundaryvalue problems (BVPs) for the Laplace and Helmholtz equations in canonical domains has been extended in order to address similar differential problems in simply connected starlike domains, whose boundaries may be regarded as an anisotropically stretched unit circle centered at the origin.
In this contribution, a suitable technique useful to compute the coefficients of the Fourierlike expansion representing the solution of the Dirichlet boundaryvalue problem for the Laplace equation in complex annular domains is presented. In particular, the boundaries of the considered domains are supposed to be defined by the socalled Gielis formula [12]. Regular functions are assumed to describe the boundary values, but the proposed approach can be easily generalized in the case of weakened hypotheses. In order to verify and validate the developed methodology, a suitable numerical procedure based on the computer algebra system Mathematica^{©} has been adopted. By using such a procedure, a pointwise convergence of the Fourierlike series representation of the solution has been observed in the regular points of the boundaries, with Gibbslike phenomena potentially occurring in the quasicusped points. The obtained numerical results are in good agreement with theoretical findings by Carleson [13].
The Laplacian in stretched polar coordinates
Let us introduce in the real plane the usual polar coordinate system
and the polar equations
relevant to the boundaries of the supershaped annulus
with
In this way, upon introducing the stretched radius ϱ such that
with
Remark Note that in the stretched coordinate system
Let us consider a piecewise function
In the considered stretched coordinate system Δ can be represented by setting
In this way, by denoting
where the dot superscript denotes the differentiation with respect to the angle ϑ. Substituting equations (8)(10) into equation (6) finally yields
As it can be easily noticed, upon setting
The Dirichlet problem for the Laplace equation
Let us consider the interior Dirichlet problem for the Laplace equation in a starlike
annulus
Under the mentioned assumptions, one can prove the following theorem.
TheoremLet
where
For each indexm, define
and set, for shortness,
where
with
Proof Upon noting that in the stretched coordinate system
Substituting into the Laplace equation, one easily finds that the functions
respectively. The parameter μ is a separation constant whose choice is governed by the physical requirement that
at any fixed point in the real plane the scalar field
where
with
Enforcing the Dirichlet boundary condition readily yields
It is worth noting that the derived expressions still hold under the assumption that
Numerical procedure
In the following numerical examples, let us assume, for the boundaries
as introduced by Gielis in [12]. Very different characteristic geometries, including ellipses, Lamé curves, ovals,
and mfold symmetric figures are obtained by assuming suitable values of the parameters
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 as follows:
with
Remark It is to be noticed that where the boundary values exhibit a rapidly oscillating behavior, the order N of the expansion (29) approximating the solution of the problem should be increased accordingly in order to achieve the desired numerical accuracy.
First example
By assuming in (27)
Figure 1. Relative boundary error
Figure 2. Boundary behavior along
Figure 3. Spatial distribution of the Fourierlike series expansion
Figure 4. Magnitude (a), (b) and phase (c), (d) of the coefficients
Second example
In the second numerical example, we turn to the consideration of the class of annuli having one or both boundaries featuring a polygonal contour. In this respect, it is not difficult to show that the general ksided convex regular polygon can be readily described by the following specialized version of Gielis’ formula [14]:
In this way, the methodology detailed in the previous section can be used straightforwardly.
In particular, upon assuming in (27)
Figure 5. Relative boundary error
Figure 6. Boundary behavior along
Figure 7. Spatial distribution of the Fourierlike series expansion
Figure 8. Magnitude (a), (b) and phase (c), (d) of the coefficients
Remark It has been observed that an
Conclusion
A Fourierlike projection method, in combination with the adoption of a suitable stretched coordinate system, has been developed for solving the Dirichlet problem for the Laplace equation in supershaped annuli. In this way, analytically based expressions of the solution of the considered class of BVPs can be derived by using classical quadrature rules, thus overcoming the need for cumbersome numerical techniques such as finitedifference or finiteelement methods. The proposed approach has been successfully validated by means of a dedicated numerical procedure based on the computeraided algebra tool Mathematica^{©}. A pointwise convergence of the expansion series representing the solution seems to be verified with the only exception of a set of measure zero consisting of the quasicusped points along the boundary of the problem domain. In these points, Gibbslike oscillations may occur. The computed results are found to be in good agreement with the theoretical findings on Fourier series.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
DC proved the main theorem regarding the solution of the Laplace equation in supershaped annuli and drafted the paper. JG carried out the verification of the methodology and its application to Gielis domains. IT performed the numerical examples. PER derived the analytical expression of the Laplacian operator in stretched coordinates and helped to draft the manuscript. All authors read and approved the final manuscript.
Acknowledgements
Dedicated to Professor Hari M Srivastava.
This research has been carried out under the grant PRIN/2006 Cap. 7320.
References

Natalini, P, Patrizi, R, Ricci, PE: The Dirichlet problem for the Laplace equation in a starlike domain of a Riemann surface. Numer. Algorithms. 28, 215–227 (2001). Publisher Full Text

Caratelli, D, Ricci, PE: The Dirichlet problem for the Laplace equation in a starlike domain. Las Vegas, 1417 July 2008. (2008)

Caratelli, D, Germano, B, Gielis, J, He, MX, Natalini, P, Ricci, PE: Fourier Solution of the Dirichlet Problem for the Laplace and Helmholtz Equations in Starlike Domains, Tbilisi University Press, Tbilisi (2010)

Caratelli, D, Natalini, P, Ricci, PE, Yarovoy, A: The Neumann problem for the Helmholtz equation in a starlike planar domain. Appl. Math. Comput.. 216, 556–564 (2010). Publisher Full Text

Caratelli, D, Gielis, J, Natalini, P, Ricci, PE, Tavkelidze, I: The Robin problem for the Helmholtz equation in a starlike planar domain. Georgian Math. J.. 18, 465–480 (2011)

Caratelli, D, Gielis, J, Ricci, PE: Fourierlike solution of the Dirichlet problem for the Laplace equation in ktype Gielis domains. J. Pure Appl. Math., Adv. Appl.. 5, 99–111 (2011)

Caratelli, D, Ricci, PE, Gielis, J: The Robin problem for the Laplace equation in a threedimensional starlike domain. Appl. Math. Comput.. 218, 713–719 (2011)

Gielis, J, Caratelli, D, Fougerolle, Y, Ricci, PE, Gerats, T: Universal natural shapes from unifying shape description to simple methods for shape analysis and boundary value problems. PLoS ONE (2012) doi:10.1371/journal.pone.0029324

Lebedev, NN: Special Functions and Their Applications, Dover, New York (1972)

Gielis, J: A generic geometric transformation that unifies a wide range of natural and abstract shapes. Am. J. Bot.. 90, 333–338 (2003). PubMed Abstract  Publisher Full Text

Carleson, L: On convergence and growth of partial sums of Fourier series. Acta Math.. 116, 135–157 (1966). Publisher Full Text

Lenjou, K: Krommen en oppervlakken van Lamé and Gielis. MSc thesis, Catholic University of Leuven (2005)