Schwarz-Christoffel Mapping
本笔记不全,跳过了部分内容,并持续更新~
Chapter 1: Introduction
1.1 The Schwarz-Christoffel Idea
Schwarz-Christoffel变换背后的想法是,任何共形映射(conformal mapping)fff都有导数:f′=∏fkf'=\prod f_kf′=∏fk,其中fkf_kfk是一系列典型(canonical)函数。这个公式的重要性在于辐角argf′=∑argfk\arg f'=\sum\arg f_kargf′=∑argfk。在classical变换中,所有argfk\arg f_kargfk都设计为一个阶跃函数,所以argf′\arg f'argf′是一个带有特定阶跃的分段常数(如fff将实轴映射到多边形上)。设PPP为多边形Γ\GammaΓ包围的复平面,顶点为w1,⋯ ,wnw_1,\cdots,w_nw1,⋯,wn,内角为α1π,⋯ ,αnπ\alpha_1\pi, \cdots, \alpha_n\piα1π,⋯,αnπ,其中αk∈(0,2)\alpha_k\in(0,2)αk∈(0,2)。设共形映射fff将上半平面H+H^+H+映射到PPP上,并设zk=f−1(wk)z_k=f^{-1}(w_k)zk=f−1(wk)是第kkk个prevertex。我们可以不失一般性地假设zn=∞z_n=\infinzn=∞,即实轴上的无穷映射到多边形最后一个顶点上,因为即使∞\infin∞还不是prevertex,我们也可以将其image设为一个新的内角为π\piπ的顶点。其余的prevertices z1,⋯ ,zn−1z_1,\cdots,z_{n-1}z1,⋯,zn−1都是真的。
在线段(zk,zk+1)(z_k,z_{k+1})(zk,zk+1)之间,fff是分析连续的,即f′f'f′存在,且argf′\arg f'argf′是常数并在zkz_kzk处有阶跃,即[argf′(z)]zk−zk+=(1−αk)π=βkπ[\arg f'(z)]_{z_k^-}^{z_k^+}=(1-\alpha_k)\pi=\beta_k\pi[argf′(z)]zk−zk+=(1−αk)π=βkπ,其中βkπ\beta_k\piβkπ是在第k个顶点处的转角。我们可以根据以上条件设定函数fk=(z−zk)−βkf_k=(z-z_k)^{-\beta_k}fk=(z−zk)−βk。下图示意了实轴经过fkf_kfk映射的结果。
定理1.1: 设P是多边形Γ\GammaΓ的内部,顶点为w1,⋯ ,wnw_1,\cdots,w_nw1,⋯,wn(逆时针),内角为α1π,⋯ ,αnπ\alpha_1\pi, \cdots, \alpha_n\piα1π,⋯,αnπ,定义fff是将上半平面H+H_+H+变换到PPP的任意共形映射,且f(∞)=wnf(\infin)=w_nf(∞)=wn,则f(z)=A+C∫z∏k=1n−1(ζ−zk)ak−1dζf(z)=A+C\int^z\prod_{k=1}^{n-1}(\zeta-z_k)^{a_k-1} d \zetaf(z)=A+C∫zk=1∏n−1(ζ−zk)ak−1dζ其中A,CA,CA,C都是复常数,且对于k=1,⋯ ,n−1k=1, \cdots ,n-1k=1,⋯,n−1都有wk=f(zk)w_k=f(z_k)wk=f(zk)。
以上公式可以用于映射不同的区域(如单位圆)、有branch points的区域、多连通区域、边界为圆弧的区域、甚至是边界分段解析的区域,在第四章中说明。
问题来了:我们不知道prevertices zkz_kzk是多少,那也就没法用定理1.1计算映射。从定理中可以发现,无论zkz_kzk是多少,多边形的内角都是不变的,这由aka_kak决定。也就是说,如果zkz_kzk选的不好,映射出来的多边形只是内角与目标多边形相同,而边长不同。zkz_kzk的确定称为Schwarz-Christoffel parameter problem,并且需要在使用SC方程之前求解。具体求解方法在第二章中讲述。
在应用中,zkz_kzk非线性地取决于多边形的边长,没有解析解,因此经常使用数值计算来计算定理1.1中的积分。数值问题在第三章中讨论。
1.2 历史
共形映射起源于19世纪。黎曼在1851年提出了黎曼映射定理,说复平面上任何单连通区域,只要不是整个复平面,就可以互相进行共形映射。Schwarz-Christoffel公式随后提出。
Christoffel在1867年于ETH Zurich发表了第一篇关于SC公式的论文,发现了在多边形域内,格林方程可以通过对半平面进行共形映射获得。随后他扩展到多边形外部与曲线区域。
Schwarz与1869年独立提出相同的理论,但侧重于数值和特定情境,如三角形;他甚至发表了首个SC映射的plot。论文中还包含了他的著名的反射原理:如果一个解析函数fff被连续扩展到一个直线或圆弧边界,并将边界弧映射到另一个直线或者圆弧上,则fff可以被解析地continued across the arc by reflection。
在SC公式发现的130年后,它在理论复分析中产生了巨大影响,用于证明黎曼映射定理与其推论。在计算机发明以后,SC公式用于具体应用。
在20世纪下半夜,有很多关于算法和程序的工作,但质量不高、综合性不强;很多重要问题都被忽视,如SC积分的快速计算,还有在解参数问题时prevertices的顺序。最综合的程序应当是Trefethen在1980年写的SCPACK和Driscoll在1996年发布的SC Toolbox。
Chapter 2: Essentials of Schwarz-Christoffel mapping
2.1 多边形
定义多边形Γ\GammaΓ,顶点为w1,⋯ ,wnw_1,\cdots,w_nw1,⋯,wn(逆时针),内角为α1π,⋯ ,αnπ\alpha_1\pi, \cdots, \alpha_n\piα1π,⋯,αnπ,并扩展定义wn+1=w1w_{n+1}=w_1wn+1=w1与w0=wnw_0=w_nw0=wn。
内角αk\alpha_kαk的定义是从wkw_kwk的incoming边到outgoing边的扫过的角度。如果∣wk∣<∞|w_k|<\infin∣wk∣<∞,则αk∈(0,2]\alpha_k\in(0,2]αk∈(0,2]。如果αk=2\alpha_k=2αk=2,则两边重合,wkw_kwk是这个裂缝(slit)的尖尖。如果wk=∞w_k=\infinwk=∞,则αk∈[−2,0]\alpha_k\in[-2,0]αk∈[−2,0]。如果wkw_kwk与其相邻顶点都是有限的,则指定αk\alpha_kαk是冗余的,但如果不是,则就需要确定αk\alpha_kαk来保证多边形唯一。
同时要求外角和为2π2\pi2π,即∑k=1nαk=n−2\sum_{k=1}^n\alpha_k=n-2∑k=1nα