Research article

Numerical computation of preimage domains for spiral slit regions and simulation of flow around bodies

  • Received: 17 August 2022 Revised: 25 September 2022 Accepted: 07 October 2022 Published: 13 October 2022
  • In this paper, we propose the iterative numerical methods to calculate the conformal preimage domains for the specified logarithmic spiral slit regions and develop the applications of conformal mappings in the simulations of the flow around bodies. Firstly, we postulate that the boundaries of the preimage domains mapped onto logarithmic spiral slits are ellipses. The lengths of the long axes of ellipses and the coordinates of the centers are calculated using our iterative methods. Secondly, each type of the presented iterative method calculates numerical conformal mappings via solving the boundary integral equation with the generalized Neumann kernel. Finally, numerical examples show the convergence and availability of our iterative methods and display the simulations of the flow around the bodies as an application.

    Citation: Kang Wu, Yibin Lu. Numerical computation of preimage domains for spiral slit regions and simulation of flow around bodies[J]. Mathematical Biosciences and Engineering, 2023, 20(1): 720-736. doi: 10.3934/mbe.2023033

    Related Papers:

    [1] V. Volpert, B. Xu, A. Tchechmedjiev, S. Harispe, A. Aksenov, Q. Mesnildrey, A. Beuter . Characterization of spatiotemporal dynamics in EEG data during picture naming with optical flow patterns. Mathematical Biosciences and Engineering, 2023, 20(6): 11429-11463. doi: 10.3934/mbe.2023507
    [2] B. Wiwatanapataphee, D. Poltem, Yong Hong Wu, Y. Lenbury . Simulation of Pulsatile Flow of Blood in Stenosed Coronary Artery Bypass with Graft. Mathematical Biosciences and Engineering, 2006, 3(2): 371-383. doi: 10.3934/mbe.2006.3.371
    [3] Chiu-Yen Kao, Yuan Lou, Eiji Yanagida . Principal eigenvalue for an elliptic problem with indefinite weight on cylindrical domains. Mathematical Biosciences and Engineering, 2008, 5(2): 315-335. doi: 10.3934/mbe.2008.5.315
    [4] Guowei Wang, Yan Fu . Spatiotemporal patterns and collective dynamics of bi-layer coupled Izhikevich neural networks with multi-area channels. Mathematical Biosciences and Engineering, 2023, 20(2): 3944-3969. doi: 10.3934/mbe.2023184
    [5] Liqiong Pu, Zhigui Lin . A diffusive SIS epidemic model in a heterogeneous and periodically evolvingenvironment. Mathematical Biosciences and Engineering, 2019, 16(4): 3094-3110. doi: 10.3934/mbe.2019153
    [6] Alessandro Bertuzzi, Antonio Fasano, Alberto Gandolfi, Carmela Sinisgalli . Interstitial Pressure And Fluid Motion In Tumor Cords. Mathematical Biosciences and Engineering, 2005, 2(3): 445-460. doi: 10.3934/mbe.2005.2.445
    [7] Radu C. Cascaval, Ciro D'Apice, Maria Pia D'Arienzo, Rosanna Manzo . Flow optimization in vascular networks. Mathematical Biosciences and Engineering, 2017, 14(3): 607-624. doi: 10.3934/mbe.2017035
    [8] Pshtiwan Othman Mohammed, Donal O'Regan, Dumitru Baleanu, Y. S. Hamed, Ehab E. Elattar . Analytical results for positivity of discrete fractional operators with approximation of the domain of solutions. Mathematical Biosciences and Engineering, 2022, 19(7): 7272-7283. doi: 10.3934/mbe.2022343
    [9] Bestami TAŞAR, Fatih ÜNEŞ, Ercan GEMİCİ . Laboratory and numerical investigation of the 2-array submerged vanes in meandering open channel. Mathematical Biosciences and Engineering, 2023, 20(2): 3261-3281. doi: 10.3934/mbe.2023153
    [10] Tingxi Wen, Hanxiao Wu, Yu Du, Chuanbo Huang . Faster R-CNN with improved anchor box for cell recognition. Mathematical Biosciences and Engineering, 2020, 17(6): 7772-7786. doi: 10.3934/mbe.2020395
  • In this paper, we propose the iterative numerical methods to calculate the conformal preimage domains for the specified logarithmic spiral slit regions and develop the applications of conformal mappings in the simulations of the flow around bodies. Firstly, we postulate that the boundaries of the preimage domains mapped onto logarithmic spiral slits are ellipses. The lengths of the long axes of ellipses and the coordinates of the centers are calculated using our iterative methods. Secondly, each type of the presented iterative method calculates numerical conformal mappings via solving the boundary integral equation with the generalized Neumann kernel. Finally, numerical examples show the convergence and availability of our iterative methods and display the simulations of the flow around the bodies as an application.



    It is difficult to deal with field theory problems defined in the region W bounded by the non-Jordan curves. However, if a conformal mapping Υ exists that maps region D bounded by smooth Jordan curves onto domain W, then we can transform the problems specified in region W into domain D to simply the problems. Meanwhile, we can refer to region D as the preimage domain of the region W. Because the preimage domains have better properties, the method for calculating the preimage domain needs to be proposed and extended. Moreover, the preimage domain we discuss is given in the sense of conformal equivalence, which should be emphasized. Conformal mappings are an important theory of analytic function and potential theory, which find in the methods of solving problems in field theory. Conformal mappings can map the complex regions onto the simple domains, which simplifies the issues such as calculating the capacitance of irregular capacitors, harmonic measures [1], logarithmic capacity [2,3] in potential theory, and so on [4,5,6,7]. Many more practical applications of conformal mappings can be found in fluid mechanics [8,9] and optics [10]. As is known to all, the analytic expressions of conformal mappings are not available in many actual situations, only can be found in some specific regions, and conformal mappings have the superiority in dealing with the problems in field theory. Therefore, the numerical methods to compute conformal mappings have been extensively investigated in the past several decades [11,12,13,14]. In comparison, the numerical method presented by Nasser, based on the boundary integral equation with the generalized Neumann kernel, has the characteristic of a unified form and can use for thirty-nine canonical slit regions proposed by Koebe [15,16,17,18]. Using the method presented by Nasser to obtain the numerical solutions of conformal mappings for highly connected regions with piecewise smooth boundaries is quick and accurate [19,20].

    Several fluid mechanics phenomena are included in external flows past objects. Computing the preimage domains of specified slit regions is of great significance in explaining some fluid phenomena. For example, Crowdy finds several explicit formulas drawn up according to the Schottky-Klein prime function for the problems in fluid dynamics [21,22,23,24]. According to Ref.[24] Crowdy (2008), an explicit solution of complex potential for the irrotational incompressible potential fluid driven by an assembly of stirrers moving at specified velocities has been found. Examples presented in [24] show streamlines of the irrotational flow generated by circular disk stirrers and paddles, where the paddles can regard as rectilinear slits. The feasible methods to solve these problems are conformally mapping the specified slit region W onto the preimage domain D with the same connectivity and solving these problems in D. However, it is difficult to directly calculate conformal mapping that maps a specified slit region onto a preimage domain, but we can find many methods which compute conformal mapping from a preimage domain onto a slit region [11,16,17,25]. Nasser [26] reformulates the problems of stirring fluid as the classical Riemann-Hilbert (R-H) problems and obtains the numerical complex potentials for the flows. When the stirrers are finite rectilinear slits, Nasser [26] presents an iterative method, which calculates a preimage domain D bounded by ellipses and conformal mapping Υ that maps D onto the rectilinear slit region. In addition, conformal mappings are also a powerful tool for simulating the flow around bodies, i.e., the phenomenon of boundary layer separation caused by the fluid passing through the surface of objects [25,26]. The flow around bodies has been widely studied as a classical problem in fluid mechanics. According to Ref.[27] Chen et al. (2021), the analytical solution of the potential flow is derived in the region bounded by a circle, an ellipse, and a thin airfoil. They find that the singular equation or the hypersingular equation can be used to solve the symmetric or antisymmetric potential flow problem alone. According to Ref.[28] Chen et al. (2022), the analytical solution of the potential flow around two identical cylinders is derived by using the boundary integral equation with the bipolar degenerate-kernel. When considering the motion of spiral point vortices in a two-dimensional multiply connected domain with two parallel slits, Aoyama et al. [29] first propose an iterative method based on the particle charge simulation method, which can calculate the preimage domain of the slit region. According to Ref.[29] Aoyama et al. (2013), many equilibria are composed of two spiral point vortices that enhance the downward vertical force and the counterclockwise rotational force.

    As a result, computing the preimage domains of the specified slit regions and conformal mappings that map the preimage domains onto the slit regions are problems worthy of attention. According to Ref.[30] Nasser (2019), the iterative numeral method has successfully computed the preimage domains D and conformal mappings Υ that map D onto W, where W is two instances of bounded slit regions called an annulus with radial slit region and the unit disk with radial slit region. However, the thirty-nine canonical slit regions cataloged by Koebe have extensive application prospects. The logarithmic spiral slit regions deserve special attention because when the oblique angle is π2 or 0, the logarithmic spiral will degenerate to a circle with a center at the origin or a ray issuing from the origin, respectively [17,25,31,32]. Meanwhile, the logarithmic spiral slit regions are closely related to the spiral point vortices.

    However, the methods proposed in [29,30] can not directly calculate the preimage domains of the logarithmic spiral slit regions. Therefore, we present the improved iterative methods based on the method proposed in [30], which calculates the preimage domains D of the specified logarithmic spiral slit regions W and conformal mappings Υ that map D onto W. We postulate that the boundaries of the sequence {Dm} of the preimage domains are the ellipses. According to the distribution of logarithmic spiral slits in W, the lengths of the long axes and the coordinates of the centers of ellipses, which calculate during iteration, respectively converge to the lengths of the long axes and the coordinates of the centers of the boundaries D. In each iteration step, to calculate conformal mappings quickly and accurately, we use the method based on the boundary integral equation with the generalized Neumann kernel to calculate conformal mappings.

    The rest of this paper is organized as follows. In Sect.2, we introduce the related knowledge about the boundary integral equation with the generalized Neumann kernel. In Sect.3, we regard two instances of spiral slit regions called the unbounded logarithmic spiral slit region and a unit disk with logarithmic spiral slit region. We briefly review the process of formulating mapping functions as the R-H problems, and then the iterative methods are presented for calculating the preimage domains and conformal mappings that map the preimage domains onto the spiral slit regions. In Sect.4, numerical examples are proposed showing the availability of our iterative methods and simulating the flow around bodies. The results are summed up in the last section.

    In this Section, we introduce the related knowledge about the boundary integral equation with the generalized Neumann kernel. Further details can be found in [14,33].

    We suppose that D is a n-connected region in the extended complex plane C{} containing the point z=0, where the boundary of D is Γ:=D=nj=1Γj. Γ1,Γ2,,Γn are n closed smooth Jordan curves. Denote the orientation of Γ to ensure D where the left of Γ.

    Suppose that the curves Γ1,Γ2,,Γn can be parameterized by the 2π-periodic complex functions η1(t),η2(t),,ηn(t) with twice continuously derivative, respectively, in which the first derivative dηj(t)/dt=˙ηj(t)0 for tJj=[0,2π],j=1,2,,n. The parametrization of the boundary Γ can be defined by

    η(t):={η1(t),tJ1=[0,2π],η2(t),tJ2=[0,2π],ηn(t),tJn=[0,2π]. (2.1)

    Definition 2.1. A real function f is Hölder continuous, when there are nonnegative real constants C, α>0, such that

    |f(x)f(y)|Cxyα

    for all x and y in the domain of f.

    We denote by H the space of real functions on Γ whose are Hölder continuous. Thus, ϕH indicates that ϕ(t):=ϕ(η(t)) is a real Hölder continuous function with 2π-period.

    Definition 2.2. A piecewise constant function h(t) is defined on Jj such that

    h(t):={h1,tJ1=[0,2π],h2,tJ2=[0,2π],hn,tJn=[0,2π], (2.2)

    where hjR,j=1,2,,n. For convenience, we denote h(t)=(h0,h1,,hn).

    Let

    B(t):={ei(π2θ(t)),tJ,ifDisunbounded,ei(π2θ(t))(η(t)α),tJ,ifDisbounded (2.3)

    is a continuously differentiable complex function. αD is a fixed point, and θ(t)=(θ1,θ2,,θn) with constants θjR for j=1,2,,n representing the oblique angle of each slit is a specified piecewise function.

    The generalized Neumann kernel N(s,t) formed with B and η is defined as follows:

    N(s,t):=1πIm[B(s)B(t)˙η(t)η(t)η(s)],(s,t)J×J.

    The singular kernel M(s,t) is defined as follows:

    M(s,t):=1πRe[B(s)B(t)˙η(t)η(t)η(s)],(s,t)J×J.

    The Fredholm integral operator N and singular integral operator M on H can be defined as follows:

    Nμ(s):=JN(s,t)μ(t)dt,sJ, (2.4)

    and

    Mμ(s):=JM(s,t)μ(t)dt,sJ. (2.5)

    Theorem 2.1. ([17])For a specified function γH, there exists a unique real function μ and a unique piecewise constant real function h=(h1,h2,,hn) defined on Γ such that

    g(η(t))=γ(t)+h(t)+iμ(t)B(t),tJ, (2.6)

    are boundary values of an analytic function g in D, which g is a unique solution of R-H problem Re[Bg]=γ+h. Then the function μ is the solution of the integral equation

    (IN)μ=Mγ (2.7)

    and the function h is obtained by

    h=Mμ(IN)γ2. (2.8)

    In this Section, we briefly review the process of formulating mapping functions as the R-H problems, and the iterative methods are presented for calculating the preimage domains D and conformal mappings that map D onto W, which W is a specified unbounded or bounded logarithmic spiral slit region (For short: spiral slit region).

    In the w plane, the boundaries of the spiral slit regions are composed of the slits that have the origin as the asymptotic point and satisfy

    Im[eiθln(w)]=rwithθ,rR. (3.1)

    If any ray issuing from the origin intersects a slit, θ named the oblique angle of the slit is the angle between the ray with the tangent of the slit at the intersection. Especially, for the oblique angle θ=0, all points of the slit will fall on a ray, which the slit is a radial slit pointing to the origin; For the oblique angle θ=π2, the ray being always perpendicular to the tangent of any point in the slit, and the slit will be a circle slit centered at the origin[16,17,25,32].

    In this Subsection, we propose an iterative method based on the boundary integral equation, which calculates a preimage domain of a specified unbounded logarithmic spiral slit region.

    Theorem 3.1. ([25]) Let θ(t)=(θ1,θ2,,θn) and a fixed point α in the unbounded multiply connected region D if Υ that conformally maps D onto the unbounded logarithmic spiral slit region satisfies normalization conditions Υ(α)=0,Υ()=,limzΥ(z)z=1, and then the Laurent series

    Υ(z)=z+a0+a1z+a2z2+

    is uniquely determined.

    Suppose D is an unbounded n-connected region in the z plane. The unbounded spiral slit region W is the entire w plane with n spiral slits having the oblique angles of θ1,θ2,,θn, respectively. According to the normalization conditions, conformal mapping Υ(z) can be defined by (see [17,25,32])

    Υ(z)=(zα)eg(z),zDΓ, (3.2)

    which the analytic function g(z) defines on DΓ. Conformal mapping ω=Υ(z) maps Γj onto the spiral slit Lj for j=1,2,,n (see Figure 1). According to the boundary conditions (3.1), the boundary values of conformal mapping Υ(z) should satisfy

    Re[ei(π2θj)ln(Υ(z))]=Rj,zΓj,j=1,2,,n. (3.3)
    Figure 1.  Υ maps the unbounded region D onto the spiral slit region W, and ζ maps the spiral slits Lj onto the line segments Sj for j=1,2,,n..

    Plugging (3.2) into (3.3) leads to

    Re[ei(π2θj)g(z)]=RjRe[ei(π2θj)ln(zα)],zΓj,j=1,2,,n, (3.4)

    where RjR for j=1,2,,n are undetermined real constants obtained by solving the R-H problem with

    B(t)=ei(π2θ(t)),γ(t)=Re[ei(π2θ(t))ln(η(t)α)],tJ. (3.5)

    According to Theorem 2.1, by solving the integral equation (2.7) and calculating (2.8), we can obtain μ(t) and h(t)=(R1,R2,,Rn) such that

    ei(π2θ(t))g(η(t))=γ(t)+h(t)+iμ(t),tJ,

    which satisfie the boundary conditions (3.4). For zD, we use the following Cauchy integral to calculate the values of g(z)

    g(z)=12πiΓγ+h+iμξzdξ,zD. (3.6)

    Together with (3.2), we can obtain conformal mapping Υ(z).

    Note that the boundaries of a specified spiral slit region are fixed. However, Rj calculated by the above method are unable to specify in advance, which affect the location distributions of the boundaries Lj for j=1,2,,n, and conformal mapping cannot be computed directly if the preimage D is unknown. Thus, the problems will become finding a preimage domain D of a specified spiral slit region W and calculating conformal mapping Υ(z) that maps D onto W.

    In the ζ plane,

    Sj=ei(π2θj)ln(w),wLj

    are n line segments with lengths lj and upper endpoints cj for j=1,2,,n (see Figure 1). Any point on a line segment can be represented as

    qj=cjiλljwith 0λ1.

    Therefore, any point on a spiral slit can be written as

    Qj=exp{ei(π2θj)qj}.

    In the following, we will propose an iterative method based on the boundary integral equation, which calculates a preimage domain D of a specified unbounded spiral slit region and derives conformal mapping Υ(z) from D onto W at the same time. Each step iteration of our iterative method will generate a preimage domain denoted by Dm, and the sequence {Dm} will converge to our needed preimage domain D.

    Assume that the boundary of the domain Dm is composed of ellipses. In the iterating step m=0,1,2,, the ellipses Γmj can are parametrized by

    ηmj(t)=Qmj+amjeiarg(Qmj)(cos(t)+ibsin(t)),tJ,j=1,2,,n, (3.7)

    where

    Qmj denote the coordinates of the centers of ellipses Γmj for j=1,2,,n,

    amj represent half the lengths of the long axes of elliptical boundaries Γmj for j=1,2,,n,

    b(0,1) is specified uniformly, which denotes the ratio between the lengths of the short axes and the long axes of ellipses Γmj for j=1,2,,n.

    Qmj and amj can be derived by the following iterative method.

    Algorithm 1 Iterative method for calculating the preimage domain of the unbounded logarithmic spiral slit region
    Input: λ, Max, ε1, β, θ(t)=(θ1,θ2,,θn),a0j=0.5(10.5b)lj,Q0j=Qj, for j=1,2,,n.
    Iterate: For m=1,2,, Max
      Step 1 According to (3.7) computing ηm1j and solving the R-H problem constructed by (3.5), we obtain conformal mapping Υm1 from the unbounded Dm1 onto an unbounded domain Wm with logarithmic spiral slits Lmj for j=1,2,,n.
      Step 2 Compute the lengths lmj and upper endpoints cmj of Smj to get the points Qmj on Lmj, for j=1,2,,n.
      Step 3 If the following condition is satisfied, the iteration ends.
                Em=1nnj=1(|lmjlj|+|QmjQj|)<ε1.(3.8)
      If (3.8)is not satisfied, proceed to the next step.
      Step 4 Compute amj=am1j12β(10.5b)(lmjlj),Qmj=Qm1j1β(QmjQj)for j=1,2,,n. Go back to the step 1.
    Output: The preimage domain D=Dm1 and conformal mapping Υ=Υm1.

    Remark: The displacement of the ellipses centers and the incrementing of the half lengths of the long axes are controlled by giving different values of β in iterations. When β=1, the Algorithmic 1 is the iteration method proposed in [30]. Max and ε1 are the specified maximum number of iterations allowed and the tolerance, respectively.

    In this Subsection, we propose an iterative method based on the boundary integral equation, which calculates a preimage domain of a specified unit disk with logarithmic spiral slit region.

    Theorem 3.2. ([32]) For the specified oblique angles θ(t)=(π2,θ2,,θn) of spiral slits, suppose D is the bounded multiply connected region inside the unit disc. Γ1={|z|=1} is a boundary component of D and 0D. Then there exists a unique analytical function Υ(z) in DΓ, which satisfies the normalization conditions Υ(0)=0,Υ(0)>0 and conformally maps D onto the spiral slit region W bounded by L1={|w|=1}. The spiral slits L2,L3,,Ln are inside of L1.

    Conformal mapping conforming to Theorem 3.2 can be defined by(see [17,32])

    Υ(z)=czezg(z),zDΓ, (3.9)

    where the analytic function g(z) defines on DΓ, and c is an undetermined real constant. w=Υ(z) maps Γj onto the spiral slit Lj for j=1,2,,n (see Figure 2). According to the boundary conditions (3.1), the boundary values of Υ(z) should satisfy

    {Re[ln(Υ(z))]=0,zΓ1,Re[ei(π2θj)ln(Υ(z))]=Rj,zΓj,j=2,3,,n. (3.10)
    Figure 2.  Υ maps the bounded multiply connected region onto the unit disk with spiral slit domian, and ζ maps the spiral slits Lj onto the line segments Sj for j=2,,n..

    Plugging (3.9) into (3.10) leads to

    {Re[zg(z)]=ln(c)Re[ln(z)],zΓ1,Re[ei(π2θj)zg(z)]=hjRe[ei(π2θj)ln(z)],zΓj,j=2,3,,n, (3.11)

    where c and hj=RjRe[ei(π2θj)ln(c)] for j=2,3,,n are undetermined real constants obtained by solving the R-H problem with

    B(t)=ei(π2θ(t))η(t),γ(t)=Re[ei(π2θ(t))ln(η(t))],tJ. (3.12)

    According to Theorem 2.1, by solving the integral equation (2.7) and calculating (2.8), we can obtain μ(t) and h(t)=(ln(c),h2,,hn) such that

    ei(π2θ(t))η(t)g(η(t))=γ(t)+h(t)+iμ(t),tJ,

    which satisfie the boundary conditions (3.11). For zD, we also use the Cauchy integral (3.6) to calculate the values of g(z), and conformal mapping Υ(z) is obtained according to (3.9).

    The images of the spiral slits Lj under the mapping

    ζ=ei(π2θj)ln(w),wLj,j=2,3,,n

    are still the line segments Sj for j=2,3,,n in the ζ plane (see Figure 2). Thus we uniformly use lj and cj for j=2,3,,n to denote the lengths and upper endpoints of Sj. Any point on a line segment can be represented as

    qj=cjiλljwith 0λ1.

    Therefore, any point on a spiral slit still can be written as

    Qj=exp{ei(π2θj)qj}.

    As mentioned in Sect.3.1, for D is unknown, we will present an iterative method for computing a preimage domain D of a specified unit disk with spiral slit region W and conformal mapping ω=Υ(z) that maps D onto W. Each step iteration of our iterative method will also generate a preimage domain Dm, and the sequence {Dm} of the preimage domains will converge to our needed preimage domain D. For m=0,1,2,, we suppose the region Dm is inside the unit circle Γm1 that is parametrized by

    ηm1=cos(t)+isin(t),tJ, (3.13)

    and is outside the ellipses Γmj parametrized by

    ηmj(t)=Qmj+amjeiarg(Qmj)(cos(t)+ibsin(t)),tJ,j=2,3,,n, (3.14)

    where

    Qmj denote the coordinates of the centers of ellipses Γmj for j=2,3,,n,

    amj represent half the lengths of the long axes of elliptical boundaries Γmj for j=2,3,,n,

    b(0,1) is specified uniformly, which denotes the ratio between the lengths of the short axes and the long axes of ellipses Γmj for j=2,3,,n.

    Qmj and amj can be derived by the following iterative method.

    Algorithm 2 Iterative method for calculating the preimage domain of the unit disk with logarithmic spiral slit region
    Input: λ, Max, ε2, β, θ(t)=(π2,θ2,,θn),a0j=0.1(10.5b)lj,Q0j=Qj, for j=2,3,,n.
    Iterate: For m=1,2,, Max
      Step 1 According to (3.13) and (3.14) computing ηm1j and solving the R-H problem constructed by (3.12), we obtain conformal mapping Υm1 from Dm1 onto an unit disk domain Wm with logarithmic spiral slits Lmj for j=2,3,,n.
      Step 2 Compute the lengths lmj and upper endpoints cmj of Smj to get the points Qmj on Lmj, for j=2,3,,n.
      Step 3 If the following condition is satisfied, the iteration ends.
              Em=1n1nj=2(|lmjlj|+|QmjQj|)<ε2.(3.15)
      If (3.15)is not satisfied, proceed to the next step.
      Step 4 Compute amj=am1j110β(10.5b)(lmjlj),Qmj=Qm1j1β(QmjQj)for j=2,3,,n. Go back to the step 1.
    Output: The preimage domain D=Dm1 and conformal mapping Υ=Υm1.

    Remark: The displacement of the ellipses centers and the incrementing of the half lengths of the long axes are controlled by giving different values of β in iterations. Max and ε2 are the specified maximum number of iterations allowed and the tolerance, respectively.

    In this Section, we propose three numerical examples to test the availability of our iterative methods.

    It is of great significance to reach the flow around bodies in energy and power engineering, aerospace engineering, and ocean engineering. Conformal mappings have unique superiorities in simulating the flow around bodies. In three numerical examples, we will simulate the irrotational plane flow around bodies in region D. Let c(z)=ϕ(x,y)+iψ(x,y)and C(w) are the complex potentials of the flows in the fluid region D and W, respectively, where ϕ and ψ denote the velocity potential and the stream function. w=Υ(z) conformally maps D onto the slit region W. At the same time, the family of equipotential lines and streamlines in D will also be mapped by w=Υ(z) onto the equipotential lines and streamlines in W. Thus, when C(w) is known, the complex potential in D is given by

    C(w)=C(Υ(z))=c(z),zD.

    The boundary integral equation (2.7) and the function (2.8) can be calculated with precision by means of the MATLAB function fbie. We set the parameters restart=[], iprec=5, maxit=100, and gmrestol=0.5×1014, in fbie. N=2048 is the number of equidistant nodes in the discretization of each boundary component. Using the MATLAB function fcau, we can rapidly derive the values of g(z) for zD. For more details, see [19].

    In numerical examples, the figures illustrate the availability of our iterative methods and discuss the influences of β on the number of iterations needed for the iterative methods to converge:

    Figure 3, Figure 6 and Figure 9 show that both the specified slit region and the calculated preimage domain are composed of boundaries and orthogonal grids. Furthermore, using the mapping ζ(w), the boundaries of the given spiral slit domains are mapped onto the line segments in Figure 3(c), Figure 6(c) and Figure 9(c).

    Figure 3.  The computed preimage domain in Figure 3(a), the specified unbounded logarithmic spiral slit region in Figure 3(b) and the spiral slits are mapped onto the line segments in Figure 3(c) for Example 4.1 with b=0.5,λ=17,β=2..

    Figure 4(b) and Figure 8 show the influences generated by the different values of b and β on the number of iterations needed for the iterative methods to converge for Example 4.1 and Example 4.2.

    Figure 4.  Streamlines of the flow around bodies generated by a spiral point vortex in Figure 4(a) and the number of iterations needed for the iterative methods to converge with the different values of b and β in Figure 4(b) for Example 4.1.

    Figure 5 and Figure 10(b) show the influences of several values of the ratio b on the number of iterations needed for the iterative methods to converge and display the error variation in the iterative process, in which the value of β is optimal for Example 4.1, Example 4.2, and Example 4.3.

    Figure 5.  The number of iterations needed for the iterative method to converge as well as Ek variation in the iterative process with β=2 for the Example 4.1 in Figure 5(a) and Example 4.2 in Figure 5(b).

    Figure 4(a), Figure 7 and Figure 10(a) show the results of simulating the irrotational plane flow around bodies.

    Example 4.1 Regard the unbounded domain that is the entire w-plane with three logarithmic spiral slits for θ(t)=(π4,π4,π4) (see Figure 3(b)). When using the tolerance ε1=1×1010, the iterative method converges for 0.05b0.8 (see Figure 5(a)). Figure 4(a) shows streamlines of the flow around bodies generated by a spiral point vortex C(w)=55i2πiln(w).

    Example 4.2 Regard the unbounded domain that is the entire w-plane with five logarithmic spiral slits for θ(t)=(π4,π4,π4,π4,π4) (see Figure 6). When using tolerance ε1=1×1010 and 0.05b0.45, the iterative method should be convergent (see Figure 8). Figure 7 shows streamlines and equipotential lines of the flows around bodies generated by a uniform stream C(w)=w and a source point C(w)=ln(w), respectively.

    Figure 6.  The computed preimage domain in Figure 6(a), the specified unbounded logarithmic spiral slit region in Figure 6(b) and the spiral slits are mapped onto the line segments in Figure 6(c) for Example 4.2 with b=0.4,λ=19,β=2.
    Figure 7.  Streamlines of the flow around bodies generated by a uniform stream in Figure 7(a) and equipotential lines of the flow around bodies generated by a source point in Figure 7(b) for Example 4.2.
    Figure 8.  The number of iterations needed for the iterative methods to converge with the different values of b and β for Example 4.2.

    Example 4.3 Regard the unit disk with logarithmic spiral slit region that oblique angles are π4 (see Figure 9). When using tolerance ε2=1×104 and 0.05b0.25, the iterative method should be convergent. At the same time, equipotential lines of the flow around bodies generated by ten source points C(w)=nj=2ln(wQj) are shown in Figure 10.

    Figure 9.  The computed preimage domain in Figure 9(a), the unit disk with logarithmic spiral slit region in fFigure 9(b) and the spiral slits inside of circle are mapped onto the line segments in Figure 9(c) for Example 4.3 with b=0.2,λ=19,β=20.
    Figure 10.  Equipotential lines of the flow around bodies generated by ten source points in Figure 10(a) and the number of iterations needed for the iterative method to converge as well as Ek variation in the iterative process with β=20 in Figure 10(b) for Example 4.3.

    Figure 4(b), Figure 5, Figure 8 and Figure 10(b) illustrate the availability and convergence of the iterative methods proposed in Sect.3. Figure 4(b) and Figure 8 show that the values of β lead to a significant impact on the number of iterations needed for the iterative methods to converge. In Example 4.1 and Example 4.2, when β=1, the iterative methods are non-convergence for some ratio b. Therefore, it is necessary to choose an appropriate value of β, which can improve the convergence of our iterative methods and reduce the number of iterations needed for the iterative methods to converge.

    In this paper, for a specified multiply connected logarithmic spiral slit region, we propose an iterative method that successfully calculates a conformally equivalent preimage domain and conformal mapping that maps the preimage domain onto the specified spiral slit region. In the proposed iterative method, we introduce parameter β to control the changes to the boundaries in the preimage domain, ensuring the convergence of that approach. Through several numerical examples, we also explore how parameters β affects the proposed iterative method's convergence and show how many iterations are necessary when parameter β=2 and β=20 for unbounded and bounded logarithmic spiral slit regions, respectively. Moreover, in physical planes, using the above conformal mappings, we simulate the several flows around bodies, including the spiral point vortex, the uniform stream, and the source point, which develop the applications of conformal mappings in the simulations of the flow around bodies. Since our method can calculate the complex potentials of spiral point vortices and other potential flows at any point in the logarithmic spiral slit regions, it provides a possibility to explore the dynamic behavior of potential flows in that domain in the future.

    KW is grateful to Prof. Mohamed M. S. Nasser for his valuable suggestions. This work was supported by the National Natural Science Foundation of China (NSFC) (Grant No. 11461037) and Yunnan Fundamental Research Projects (Grant NO. 202101BE070001-050).

    The authors declare that they have no conflict of interest.



    [1] M. Nasser, M. Vuorinen, Computation of conformal invariants, Appl. Math. Comput., 389 (2021), 125617. https://doi.org/10.1016/j.amc.2020.125617 doi: 10.1016/j.amc.2020.125617
    [2] L. V. Ahlfors, Conformal Invariants: Topics in Geometric Function Theory, New York: McGraw-Hill, (1973). https://doi.org/10.1090/chel/371 doi: 10.1090/chel/371
    [3] C. Pommerenke, On the logarithmic capacity and conformal mapping, Duke Math. J., 35 (1968), 321–325. https://doi.org/10.1215/S0012-7094-68-03531-X doi: 10.1215/S0012-7094-68-03531-X
    [4] L. Bourchtein, Conformal mappings of multiply connected domains onto canonical domains using the Green and Neumann functions, Complex Var. Elliptic Equ., 58 (2013), 821–836. https://doi.org/10.1080/17476933.2011.622045 doi: 10.1080/17476933.2011.622045
    [5] D. Crowdy, Conformal slit maps in applied mathematics, Anziam J., 53 (2012), 171–189. https://doi.org/10.1017/S1446181112000119 doi: 10.1017/S1446181112000119
    [6] N. Hale, T. Wynn Tee, Conformal maps to multiply slit domains and applications, SIAM J. Sci. Comput., 31 (2009), 3195–3215. https://doi.org/10.1137/080738325 doi: 10.1137/080738325
    [7] X. Y. Liu, W. Li, M. Li, C. S. Chen, Circulant matrix and conformal mapping for solving partial differential equations, Comput. Math. Appl., 68 (2014), 67–76. https://doi.org/10.1016/j.camwa.2014.05.005 doi: 10.1016/j.camwa.2014.05.005
    [8] K. Amano, D. Okano, Numerical conformal mappings onto the canonical slit domains, Theor. Appl. Mechan. Japan, 60 (2012), 317–332. https://doi.org/10.11345/nctam.60.317 doi: 10.11345/nctam.60.317
    [9] W. Choi, Nonlinear surface waves interacting with a linear shear current, Math. Comput. Simul., 80 (2009), 29–36. https://doi.org/10.1016/j.matcom.2009.06.021 doi: 10.1016/j.matcom.2009.06.021
    [10] H. Li, Y. D. Xu, Q. N. Wu, H. Y. Chen, Carpet cloak from optical conformal mapping, Sci. China Inform. Sci., 56 (2013), 120411. https://doi.org/10.1007/s11432-013-5036-x doi: 10.1007/s11432-013-5036-x
    [11] K. Amano, A charge simulation method for numerical conformal mapping onto circular and radial slit domains, SIAM J. Sci. Comput., 19 (1998), 1169–1187. https://doi.org/10.1137/S1064827595294307 doi: 10.1137/S1064827595294307
    [12] K. Amano, Numerical conformal mapping onto the radial slit domains by the charge simulation method, Japan Soc. Industr. Appl. Math., 5 (1995), 267–280.
    [13] L. Trefethen, Numerical conformal mapping with rational functions, Comput. Methods Function Theory, 20 (2020), 369–387. https://doi.org/10.1007/s40315-020-00325-w doi: 10.1007/s40315-020-00325-w
    [14] M. Nasser, Numerical conformal mapping via a boundary integral equation with the generalized Neumann kernel, SIAM J. Sci. Comput., 31 (2009), 1695–1715. https://doi.org/10.1137/070711438 doi: 10.1137/070711438
    [15] M. Nasser, Fast computation of the circular map, Comput Methods Function Theory, 15 (2015), 187–223. https://doi.org/10.1007/s40315-014-0098-3 doi: 10.1007/s40315-014-0098-3
    [16] M. Nasser, Numerical conformal mapping of multiply connected regions onto the fifth category of Koebe's canonical slit regions, J. Math. Anal. Appl., 398 (2013), 729–743. https://doi.org/10.1016/j.jmaa.2012.09.020 doi: 10.1016/j.jmaa.2012.09.020
    [17] M. Nasser, Numerical conformal mapping of multiply connected regions onto the second, third and fourth categories of Koebe's canonical slit domains, J. Math. Anal. Appl., 382 (2011), 47–56. https://doi.org/10.1016/j.jmaa.2011.04.030 doi: 10.1016/j.jmaa.2011.04.030
    [18] M. Nasser, A. Murid, Z. Zamzamir, A boundary integral method for the Riemann–Hilbert problem in domains with corners, Complex Var. Elliptic Equ., 53 (2008), 989–1008. https://doi.org/10.1080/17476930802335080 doi: 10.1080/17476930802335080
    [19] M. Nasser, Fast solution of boundary integral equations with the generalized Neumann kernel, Electron. Transact. Numer. Ana., 44 (2015), 189–229.
    [20] M. Nasser, F. Al-Shihri, A fast boundary integral equation method for conformal mapping of multiply connected regions, SIAM J. Sci. Comput., 35 (2013), A1736–A1760. https://doi.org/10.1137/120901933 doi: 10.1137/120901933
    [21] D. Crowdy, A new calculus for two-dimensional vortex dynamics, Theor. Comput. Fluid Dynam., 24 (2010), 9–24. https://doi.org/10.1007/s00162-009-0098-5 doi: 10.1007/s00162-009-0098-5
    [22] D. Crowdy, Analytical solutions for uniform potential flow past multiple cylinders, European J. Mechan. B/Fluids, 25 (2006), 459–470. https://doi.org/10.1016/j.euromechflu.2005.11.005 doi: 10.1016/j.euromechflu.2005.11.005
    [23] D. Crowdy, Calculating the lift on a finite stack of cylindrical aerofoils, Proceed. Royal Soc. A Math. Phys. Eng. Sci., 462 (2006), 1387–1407. https://doi.org/10.1098/rspa.2005.1631 doi: 10.1098/rspa.2005.1631
    [24] D. Crowdy, Explicit solution for the potential flow due to an assembly of stirrers in an inviscid fluid, J. Eng. Math., 62 (2008), 333–344. https://doi.org/10.1007/s10665-008-9222-6 doi: 10.1007/s10665-008-9222-6
    [25] K. Amano, D. Okano, K. Endo, H. Ogata, Numerical conformal mappings onto the canonical slit domains of Koebe (1916) by the charge simulation method, Japan Soc. Indust. Appl. Math., 24 (2014), 157–183.
    [26] M. Nasser, C. C. Green, A fast numerical method for ideal fluid flow in domains with multiple stirrers, Nonlinearity, 31 (2018), 815–837. https://doi.org/10.1088/1361-6544/aa99a5 doi: 10.1088/1361-6544/aa99a5
    [27] J. T. Chen, J. H. Kao, Y. L. Huang, S. K. Kao, Applications of degenerate kernels to potential flow across circular, elliptical cylinders and a thin airfoil, European J. Mechan. B/Fluids, 90 (2021), 29–48. https://doi.org/10.1016/j.euromechflu.2021.07.012 doi: 10.1016/j.euromechflu.2021.07.012
    [28] J. T. Chen, Y. T. Chou, J. H. Kao, J. W. Lee, Analytical solution for potential flow across two circular cylinders using the BIE in conjunction with degenerate kernels of bipolar coordinates, Appl. Math. Letters, 132 (2022), 108137. https://doi.org/10.1016/j.aml.2022.108137 doi: 10.1016/j.aml.2022.108137
    [29] N. Aoyama, T. Sakajo, H. Tanaka, A computational theory for spiral point vortices in multiply connected domains with slit boundaries, Japan J. Industr. Appl. Math., 30 (2013), 485–509. https://doi.org/10.1007/s13160-013-0113-5 doi: 10.1007/s13160-013-0113-5
    [30] M. Nasser, Numerical computing of preimage domains for bounded multiply connected slit domains, J. Sci. Comput., 78 (2019), 582–606. https://doi.org/10.1007/s10915-018-0784-9 doi: 10.1007/s10915-018-0784-9
    [31] P. Koebe, Abhandlungen zur Theorie der konformen Abbildung, Ⅳ. Abbildung mehrfach zusammenhtingender schlichter Bereiche auf Schlitzbereiche, Acta Math., 41 (1916), 305–344. https://doi.org/10.1007/BF02422949 doi: 10.1007/BF02422949
    [32] G. C. Wen, Conformal Mappings and Boundary Value Problems, Providence: American Mathematical Society, (1992).
    [33] R. Wegmann, M. Nasser, The Riemann–Hilbert problem and the generalized Neumann kernel on multiply connected regions, J. Comput. Appl. Math., 214 (2008), 36–57. https://doi.org/10.1016/j.cam.2007.01.021 doi: 10.1016/j.cam.2007.01.021
  • This article has been cited by:

    1. Changqing Shen, Sihao Mao, Bingzhou Xu, Ziwei Wang, Xiaojian Zhang, Sijie Yan, Han Ding, Spiral complete coverage path planning based on conformal slit mapping in multi-connected domains, 2024, 0278-3649, 10.1177/02783649241251385
    2. Dan Luo, Yibin Lu, Xin Zhao, Computational method of conformal mapping from unbounded multi-connected regions onto annulus with spiral slit domains, 2025, 2964, 1742-6588, 012058, 10.1088/1742-6596/2964/1/012058
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1912) PDF downloads(98) Cited by(2)

Figures and Tables

Figures(10)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog