    # Method of fundamental solutions for the problem of doubly-periodic potential flow

In this paper, we propose a method of fundamental solutions for the problem of two-dimensional potential flow in a doubly-periodic domain. The solution involves a doubly-periodic function, to which it is difficult to give an approximation by the conventional method of fundamental solutions. We propose to approximate it by a linear combination of the periodic fundamental solutions, which has sources in a doubly-periodic array and contains the complex logarithmic functions and the theta functions. Numerical examples show the effectiveness of our method.

06/26/2020

### A method of fundamental solutions for doubly-periodic potential flow problems using the Weierstrass elliptic function

In this paper, we propose a method of fundamental solutions for the prob...
12/04/2021

### Fast Electromagnetic Validations of Large-Scale Digital Coding Metasurfaces Accelerated by Recurrence Rebuild and Retrieval Method

The recurrence rebuild and retrieval method (R3M) is proposed in this pa...
02/09/2021

06/28/2020

### A numerical study of variational discretizations of the Camassa-Holm equation

We present two semidiscretizations of the Camassa-Holm equation in perio...
10/27/2021

### Almost periodic functions and an analytical method of solving the number partitioning problem

In the present paper, we study the limit sets of the almost periodic fun...
03/02/2020

### Modified Bee Colony optimization algorithm for computational parameter identification for pore scale transport in periodic porous media

This paper discusses an optimization method called Modified Bee Colony a...
08/12/2020

### The AAAtrig algorithm for rational approximation of periodic functions

We present an extension of the AAA (adaptive Antoulas–Anderson) algorith...

## 1 Introduction

The method of fundamental solutions, or the charge simulation method, [6, 15] is a fast numerical solver for potential problems

 {−△u=0in  Du=fon  ∂D, (1)

where is a domain in the -dimensional Euclidean space , and is a function given on . In two dimensional problems , equalizing the Euclid plane with the complex plane , the method of fundamental solutions gives an approximate solution of (1) of the form 222The scheme shown here is the invariant scheme [16, 17] proposed by Murota.

 u(z)≃uN(z)=Q0−12πN∑j=1Qjlog|z−ζj|(z=x+iy), (2)

where are unknown real coefficients such that

 N∑j=1Qj=0, (3)

and are points given in . We call the “charges” and the “charge points”. We remark that the approximate solution exactly satisfies the Laplace equation in . Regarding the boundary condition, we pose the collocation condition on , namely, we assume that satisfies the equations

 uN(zi)=f(zi),i=1,…,N (4)

with given on the boundary , which are called the “collocation points”. The equations (4) are rewritten as

 Q0−12πN∑j=1Qjlog|zi−ζj|=f(zi),i=1,…,N, (5)

which, together with (3), form a system of linear equations with respect to . We determine the unknowns by solving the linear system (3) and (5) and obtain the approximate solution . The method of fundamental solutions has the advantages that it is easy to program, its computational cost is low, and it shows fast convergence such as exponential convergence  under some condition. It was first used for studies of electric field problems [25, 26], and now it is widely used in science and engineering.

The method of fundamental solutions is also used for the approximation of complex analytic functions. The idea of the method of fundamental solutions is applied to various problems other than potential problems, for example, problem of scattering of earthquake waves . Let be an analytic function in a domain . The real part of , which is a harmonic function in , can be approximated using the form (2), and the imaginary part of , which is the conjugate harmonic function of , is approximated using

 −12πN∑j=1Qjarg(z−ζj).

Then, the analytic function is approximated using a linear combination of the complex logarithmic functions

 Q0−12πN∑j=1Qjlog(z−ζj). (6)

From this point of view, Amano applied the method of fundamental solution to numerical conformal mappings [1, 2].

In this paper, we examine the problem of two-dimensional potential flow past an infinite doubly-periodic array of obstacles as shown in Figure 1. A two-dimensional potential flow is characterized by a complex velocity potential , which is an analytic function in the flow domain such that it gives the velocity field by , and its imaginary part satisfies the boundary condition

 Imf=constanton  ∂D. (7)

Physically, the boundary condition (7) means that the fluid flows along the boundary since the contour lines of are the streamlines. Therefore, in order to obtain the potential flow of our problem, we have to find an analytic function in satisfying the boundary condition (7).

We have, however, one problem in applying the method of fundamental solutions to our problem. The velocity field is obviously a doubly-periodic function due to the the double periodicity of the flow domain . Then, the complex velocity potential involves a doubly-periodic function, and it is difficult to approximate it using the form (6) of the conventional method of fundamental solutions. To overcome this challenge, we propose to solve our problem by a method of fundamental solutions using a doubly-periodic fundamental solution, that is, a logarithmic potential with a doubly-periodic array of sources. We construct a doubly-periodic fundamental solution using the theta functions and approximate the complex velocity potential by a linear combination of the doubly-periodic fundamental solutions.

The previous works related to this paper are as follows. As studies of problems of periodic flow, Zick and Homsy  proposed an integral equation method for three-dimensional Stokes flow problems with a three-dimensional periodic array of spheres, where the solution is given by an integral including the periodic fundamental solution. Greengard and Kropinski  proposed an integral equation method for two-dimensional Stokes flow problems in doubly-periodic domains, where an approximate solution is given as a complex variable formulation and the fast multipole method is used. Liron  studied Stokes flow due to infinite array of Stokeslets and applied it to the problems of fluid transport by cilia. As studies on methods of fundamental solutions applied to problems with periodicity, Ogata et al. proposed a method of numerical conformal mappings of complex domains with a single periodicity , where the mapping function, an analytic function involving a singly-periodic function, is approximated by the method of fundamental solution using the singly-periodic fundamental solutions. They proposed a method of fundamental solutions also to the problems of two-dimensional periodic Stokes flows [18, 19], where the solutions are approximated by the periodic fundamental solutions of the Stokes equation, that is, the Stokes flows induced by a periodic array of point forces. The author proposed a method of fundamental solutions for the two-dimensional elasticity problem with one-dimensional periodicity , where the solution is approximated using the periodic fundamental solutions of the elastostatic equation, that is, the displacements induced by concentrated forces in a periodic array.

The remainder of this paper is structured as follows. Section 2 proposes a method of fundamental solutions for our problem of potential flow with double periodicity. Section 3 shows some numerical examples which show the effectiveness of our method. In Section 4, we conclude this paper and present problems related to future studies.

## 2 Method of fundamental solutions

We consider a potential flow past the doubly-periodic array of obstacle. The flow domain is mathematically given by

 D=C∖⋃m,n∈Z¯¯¯¯¯¯¯¯¯¯Dmn,

where is one of the obstacles, a simply-connected domain and

with complex numbers and giving the periods of the obstacle array such that . We assume that the spatial average of the velocity field is a uniform flow, in whose direction we take the real axis, when the flow is observed on the direction of . Mathematically, this assumption is formulated as

 ∫Cv⋅dr=U(Reω1),∫Cv⋅nds=U(Imω1), (8)

where is the magnitude of the velocity of the unit flow and is the line segment connecting two points and

is the unit normal vector on

. Figure 1: Flow domain D with an infinite doubly-periodic array of obstacles Dmn,m,n∈Z.

As mentioned in the previous section, a two-dimensional potential flow is characterized by a complex velocity potential, , which is an analytic function in the flow domain such that it gives the velocity field by

 f′(z)=u−iv, (9)

and its imaginary part is constant on the boundary of the domain . The condition (8) that the average of the flow is uniform is summarized as

 ∫Cdf(z)=∫Cf′(z)dz=Uω1. (10)

Therefore, our potential flow problem is mathematically reduced to the problem to find an analytic function in satisfying

 Imf=constanton  ∂D (11)

and (10).

In solving our problem by the method of fundamental solution method, however, we have one problem. The complex velocity potential involves a periodic function since the velocity field given by (9) is obviously a periodic function, and it is impossible to approximate by the conventional method of fundamental solutions using the form (2). We propose to approximate the complex velocity potential by

 f(z)≃fN(z)=Uz−i2πN∑j=1Qjlogϑ1(z−ζjω1∣∣∣τ), (12)

where is the theta function  333 In , the theta function is defined by

 ϑ1(v|τ)= 2∞∑n=0q(n+1/2)2sin(2n+1)πv = 2q1/4sinπv∞∏n=1(1−q2n)(1−2q2ncos2πv+q4n)

with

 τ=ω2ω1,q=eiπτ,

are points given in and are unknown real coefficients such that

 N∑j=1Qj=0. (13)

We call the “charge” and the “charge points”. The theta function is an entire function which has simple zeros at , and satisfies the pseudo-periodicity

 ϑ1(v+1|τ)=−ϑ1(v|τ),ϑ1(v+τ|τ)=−q−1e−2πivϑ1(v|τ). (14)

Therefore, the functions appearing on the right-hand side of (12) can be regarded as the complex logarithmic potential with sources at the points and it is a periodic fundamental solution of the two-dimensional Poisson equation 444 Hasimoto  presented a periodic fundamental solution of the two-dimensional Poisson equation in terms of the Weierstrass elliptic functions .

The approximate velocity potential satisfies the condition (10). In fact,

 ∫CdfN= fN(z0+ω1)−fN(z0) = Uω1−i2πN∑j=1Qj{logϑ1(z−ζjω1+τ∣∣∣τ)−logϑ1(z−ζjω1∣∣∣τ)} = Uω1−i2πN∑j=1Qjlog(−1)=Uω1,

where we used the first formula of (14) on the third equality and (13) on the fourth equality.

The approximate velocity potential satisfies the pseudo-periodicity

 fN(z+ω1)=fN(z)+Uω1,fN(z+ω2)=fN(z)+Uω2−2πω1N∑j=1Qjζj. (15)

In fact, the first equation is proved above, and the second equation is proved similarly to the first equation using the second formula of (14). Therefore, the approximate velocity given by is a doubly-periodic function, namely,

 vN(z+ω1)=vN(z),vN(z+ω2)=vN(z).

It means that the approximate complex velocity is an elliptic function of periods .

The approximate velocity potential in (12) is an analytic function in such that it satisfies the condition (10) and gives velocity field which is a doubly-periodic function of periods . Regarding the boundary condition (11), we pose a collocation condition on , namely, we assume the equation

 ImfN(zi)=C,i=1,…,N, (16)

where are given points on called the “collocation points” and is an unknown real constant. The equations (16) are rewritten as

 −12πN∑j=1Qjlog∣∣∣ϑ1(zi−ζjω1∣∣∣τ)∣∣∣−C=−U(Imzi),i=1,…,N. (17)

The equations (13) and (17) form a system of linear equations with respect to the unknowns and . We determine the unknown charges by solving this linear system and obtain the approximate velocity potential . Owing to the pseudo-periodicity (15), the approximate potential approximately satisfies the boundary condition on the boundary of other obstacle , , that is,

 ImfN(zi+mω1+nω2)=constant,i=1,…,N.

## 3 Numerical examples

In this section, we show some numerical examples which show the effectiveness of our method. All the computations were performed using programs coded in C++ with double precision. Figure 2: The streamlines of potential flows past a doubly-periodic array of cylinders with periods ω1,ω2.

Figure 2 shows the examples of flows past a doubly-periodic array of cylinders, that is, flows in the domain

 D=C∖⋃m,n∈Z¯¯¯¯¯¯¯¯¯¯Dmn, where Dmn={z∈C||z−mω1−nω2|

with a positive constant . The charge points and the collocation points are respectively taken as

 ζj=qrexp(i2π(j−1)N),zj=rexp(i2π(j−1)N),j=1,…,N, (18)

where is a constant such that and was taken as in the examples shown in Figure 2

. To estimate the accuracy of our method, we evaluated the value

 ϵN=1Urmaxz∈∂D00|ImfN(z)−C|,

where is the constant obtained in solving the system of linear equations (13) and (16). The value shows how accurately the approximate potential satisfies the boundary condition (7). Figure 3 shows evaluated for the above numerical example with and for , which appears in giving the charge points by (18). The figure shows that the approximation of our method converges exponentially as the number of unknowns increases. Figure 3: The error estimate ϵN of the method of fundamental solutions.

## 4 Concluding Remarks

In this paper, we proposed a method of fundamental solutions for the problems of two-dimensional potential flow past a doubly-periodic array of obstacles. In terms of mathematics, our problem is to find the complex velocity potential, an analytic function in the flow domain with double periodicity. The solution obviously involves a doubly-periodic function, and it is difficult to approximate it by the conventional method of fundamental solutions. Then, we proposed a new method of fundamental solution for this problem using the periodic fundamental solutions, which is the logarithmic potential with a doubly-periodic array of sources and contains the theta functions. The numerical examples showed the effectiveness of our method.

We have two issues regarding this paper for future study. The first problem is a theoretical error estimate of our method. Theoretical studies on the accuracy of the method of fundamental solution have been presented for special cases such as two-dimensional Laplace equation and Helmholtz equation in a disk [11, 4, 20] and two-dimensional Laplace equation in a domain with an analytic boundary [12, 21]. However, theoretical error estimate still remains unknown for many types of problems and methods including our method. The author believes it one of the most important works on the method of fundamental solutions to give a theoretical error estimate.

The second problem is to extend our method to other problems with periodicity than two-dimensional potential problems. The author has already given methods of fundamental solutions for Stokes flow problems [19, 18] and elasticity problems  with periodicity. In the work on periodic Stokes flow , the author et al. proposed to use the periodic fundamental solutions, which was presented by Hasimoto  and is given by a Fourier series. It is expected to construct a method of fundamental solutions using periodic fundamental solutions given by the theta functions or elliptic functions [10, 5] as in this paper.

## References

•  K. Amano (1994) A charge simulation method for the numerical conformal mapping of interior, exterior and doubly-connected domains. J. Comput. Appl. Math. 53 (3), pp. 353–370. External Links: Document Cited by: §1.
•  K. Amano (1998) A charge simulation method for numerical conformal mapping onto circular and radial slit domains. SIAM J. Sci. Comput. 19 (4), pp. 1169–1187. External Links: Document Cited by: §1.
•  J. V. Armitage and W. F. Eberlein (2006) Elliptic functions. Cambridge University Press, Cambridge. External Links: Document Cited by: §2, footnote 3.
•  F. Chiba and T. Ushijima (2009) Exponential decay of errors of a fundamental solution method applied to a reduced wave problem in the exterior region of a disc. J. Comput. Appl. Math. 231, pp. 869–885. External Links: Document Cited by: §4.
•  D. Crowdy and E. Luca (2018) Fast evaluation of the fundamental singularities of two-dimensional doubly periodic Stokes flow. J. Eng. Math. 111, pp. 95–110. External Links: Document Cited by: §4.
•  G. Fairweather and A. Karageorghis (1998) The method of fundamental solutions for elliptic boundary value problems. Adv. Comp. Math. 9, pp. 69–95. External Links: Document Cited by: §1.
•  L. Greengard and M. C. Kropinski (2004) Integral equation methods for stokes flow in doubly-periodic domains. J. Eng. Math. 48, pp. 157–170. External Links: Document Cited by: §1.
•  H. Hasimoto (1959) On the periodic fundamental solutions of the stokes equations and their application to viscous flow past a cubic array of spheres. J. Fluid. Mech. 5 (2), pp. 317–328. External Links: Document Cited by: §4.
•  H. Hasimoto (2008) Periodic fundamental solution of a two-dimensional Poisson equation. J. Phys. Soc. Japan 77 (10), pp. 104601. External Links: Document Cited by: footnote 4.
•  H. Hasimoto (2009) Periodic fundamental solution of the two-dimensional Stokes equations. J. Phys. Soc. Japan 78 (7), pp. 074401. External Links: Document Cited by: §4.
•  M. Katsurada and H. Okamoto (1988) A mathematical study of the charge simulation method I. J. Fac. Sci. Univ. Tokyo, Sect IA 35 (3), pp. 507–518. Cited by: §1, §4.
•  M. Katsurada (1990) Asymptotic error analysis of the charge simulation method in a jordan region with an analytic boundary. J. Fac. Sci. Univ. Tokyo, Sect. IA Math. 37, pp. 635–657. Cited by: §4.
•  N. Liron (1978) Fluid transport by cilia between parallel plates. J. Fluid Mech. 86 (4), pp. 705–726. External Links: Document Cited by: §1.
•  L. M. Milne-Thomson (2011) Theoretical hydrodynamics. Dover, New York. Cited by: §1.
•  S. Murashima (1983) Charge simulation method and its applications. Morikita-Shuppan, Tokyo. Note: (in Japanese) Cited by: §1.
•  K. Murota (1993) On “invariance” of schemes in the fundamental solution method. Trans. IPS Japan 34 (3), pp. 533–535. Note: (in Japanese) Cited by: footnote 2.
•  K. Murota (1995) Comparison of conventional and “invariant” schemes of fundamental solutions method for annular domains. Japan J. Indust. Appl. Math. 12, pp. 61–85. Cited by: footnote 2.
•  H. Ogata, K. Amano, M. Sugihara, and D. Okano (2003) A fundamental solution method for viscous flow problems with obstacles in a periodic array. J. Comput. Appl. Math. 152 (1–2), pp. 411–425. External Links: Document Cited by: §1, §4.
•  H. Ogata and K. Amano (2010) Fundamental solution method for two-dimensional stokes flow problems with one-dimensional periodicity. Japan J. Indust. Appl. Math. 27, pp. 191–215. External Links: Document Cited by: §1, §4.
•  H. Ogata, F. Chiba, and T. Ushijima (2011) A new theoretical error estimate of the method of fundamental solutions applied to reduced wave problems in the exterior region of a disk. J. Comput. Appl. Math. 235 (12), pp. 3395–3412. External Links: Document Cited by: §4.
•  H. Ogata and M. Katsurada (2014) Convergence of the invariant scheme of the method of fundamental solutions for two-dimensional potential problems in a jordan region. Japan J. Indust. Appl. Math. 31, pp. 231–262. External Links: Document Cited by: §4.
•  H. Ogata, D. Okano, and K. Amano (2002) Numerical conformal mapping of periodic structure domains. Japan J. Indust. Appl. Math. 19, pp. 257–275. External Links: Document Cited by: §1.
•  H. Ogata (2008) Fundamental solution method for periodic plane elasticity. J. Numer. Anal. Indust. Appl. Math. (JNAIAM) 3 (3–4), pp. 249–267. Cited by: §1, §4.
•  F. J. Sanchez-Sezma and E. Rosenblueth (1979) Ground motion at canyons of arbitrary shape under incident sh waves. Int. J. Earthq. Eng. Struct. Dyn. 7, pp. 441–450. Cited by: §1.
•  H. Singer, H. Steinbigler, and P. Weiss (1974) A charge simulation method for the calculation of high voltage fields. IEEE Trans. Power Appar. Syst. PAS-93, pp. 1660–1668. External Links: Document Cited by: §1.
•  H. Steinbigler (1969) Note: dissertation, Tech. Univ. München Cited by: §1.
•  A. A. Zick and G. M. Homsy (1982) Stokes flow through periodic arrays of spheres. J. Fluid Mech. 115, pp. 13–26. External Links: Document Cited by: §1.