保角变换法求解定解问题

合集下载

第6章保角变换-数学物理方法

第6章保角变换-数学物理方法

f ( z0 ) 是经过映射Biblioteka f ( z ) 后通过点z0 的
的任何曲线C在 z0的伸缩率, 它与曲线C的形状及
方向无关. 所以这种映射又具有伸缩率的不变性.
5
2.共形映射(保角映射)
设函数w f ( z )在区域 D内解析, z0为 D内一点,
且 f ( z ) 0 , 那末映射w f ( z ) 在 z0 具有两个性 质: (1) 保角性; (2) 伸缩率不变性.
圆周的弧所围成的区域映射成一圆弧与一直线所 围成的区域. 3) 当二圆交点中的一个映射成无穷远点时, 这 二圆周的弧所围成的区域映成角形区域.
13
5. 几个初等函数所构成的映射
1) 幂函数 w z n ( n 2).
映射特点: 把以原点为顶点的角形域映射成以原
点为顶点的角形域, 但张角变成为原来的 n 倍.

b 1, a 1 i ,
(1 i ) z 1 ( i 1) z 1 所以 w 为所求. z (1 i ) z (1 i )
19
解3 利用典型区域映射公式
将所求映射设为 w e i
z z A , 1 z 1 z
保角变换
复变函数在几何意义上实际上相当于将平 面上的区域变成了平面上的另一个区域(简称 为映射). 应用:利用复变函数(特别是解析函数)所构 成的映射来实现复杂区域的简单化,这将给实 际问题的研究带来很大的方便.而利用保角变 换法求解数学物理方程边值问题.
1
本章内容: 1)保角射的概念; 2)分式线性映射和几个初等函数所构成的 映射; 3)典型实例描述保角映射的应用. 重点: 分式线性变换及其映射特点 难点:

保角变换法

保角变换法

R R i c ln 2 wm t 1 R 1 R
式中
1 wm w1 w2 , c 2
平板叶栅的一般绕流
为绕一个翼型的环量。
2.3.P9
(五)平板叶栅一般流动中环量的确定
环量的确定依据是弧立翼型绕流中的库塔 —— 恰普雷金条件。而栅中翼型尾缘点 B 必然 是后驻点,此外速度是一有限值。 经换算得
a) b)
2.3.P6
其复势为
t W 2 i R 1 R i ie ln ie ln R 1 R
流动奇点强度为
q1 q2 t sin
1 2 t cos
(三)平板叶栅纯环量绕流 b) 图示,栅前后只有 列线方向速度 w1、w2 。
可见 L L t b , ,具体 数值见图示曲线。 由上述已解得的平板叶栅 流动,可以求解由任何翼型组 成的等价平面直列叶栅流动。
平板叶栅环量修正曲线
三、平面环列叶栅流动的解法
2.3.P11
设图示环列叶栅由 n 个翼型组成,流动自中心 向外。可见,只要确定一个扇形区域内的流动即可。
平板叶栅无环量平行绕流
2
2.3.P4
q t cos
t sin
Z 平面复势
W z zei
表示速度为 1 的均匀流复势。 变换为 平面为 R 处相应放置点源、点汇
q 和点涡 的绕圆流动。
其复势
t W 2 i R i 1 R e ln e ln R 1 R
变换为 平面绕单位圆流动,且有
R R i W ln 1 1 4 R R

数学物理方法 保角变换法共40页

数学物理方法 保角变换法共40页
数学物理方法 保角变换法
1、战鼓一响,法律无声。——英国 2、任何法律的根本;不,不成文法本 身就是 讲道理 ……法 律,也 ----即 明示道 理。— —爱·科 克
3、法律是最保险的头盔。——爱·科 克 4、一个国家如果纲纪不正,其国风一 定颓败 。—— 塞内加 5、法律不能使人人平等,但是在法律 面前人 人是平 等的。 ——。 ——德 谟克利 特 67、今天应做的事没有做,明天再早也 是耽误 了。——裴斯 泰洛齐 68、决定一个人的一生,以及整个命运 的,只 是一瞬 之间。 ——歌 德 69、懒人无法享受休息之乐。——拉布 克 70、浪费时间是一桩大罪过。——卢梭

数值保角变换的新算法

数值保角变换的新算法
2数值保角变换的新算法通过上节讨论,根据模拟电荷法的原理,在求解模拟电荷点以及近似变换半径过程中,首先,要用围绕C的区域内部配置N个电荷点作为极的对数势场的1次结合的G(z)来高度近似g(z),而它的共轭调和函数h(z)则用H(z)来近似,再通过边界条件以及电荷点和约束点的配置构造出约束方程(5)。因此,数值保角变换的算法过程整理如下。
在图1中,C是z平面上任意的Jordan曲线,曲线C的区域外部作为D,灼j(j=1,2,…,N)是在区域内部配置的电荷点,zi(i=1,2,…,N)是边界C上的约束点,w=(f z)是从D軍(这里D軍=D胰C)到w平面上的单位圆外部w跃1(包括单位圆边界)的保角映射。在不失一般性的情况下,假定z=0在C的内部且(f 0)=0,保角变换函数w=(f z)满足正规化条件(f肄)=肄,f(忆肄)跃0时表示如下:f(z)=z酌exp(g(z)+ih(z))其中,酌是外部变换半径,g(z)是Dirichlet型场势问题:2g(z)=0,z沂Dg(z)=log酌-log z,z沂Cg(肄)=0扇墒设设设设缮设设设设的解。h(z)是g(z)的共轭调和函数,且h(肄)=0。以下文中均以G,H,祝表示g,h,酌的近似值。
数值保角变换的算法:淤给出模拟电荷法的模拟电荷和约束点数量N;于根据模拟电荷法的原理[11],给出电荷点灼1,…,灼N,约束点z1,…,zN;盂由模拟电荷点和约束点以及变换半径构造约束方程组(5);榆计算模拟电荷q1,…,qN和变换半径祝;虞根据(1)式和(2)式构造G(),H(z);愚构造近似保角变换函数F(z)=z祝exp(G(z)+iH(z))。
数值保角变换的新算法
摘要院保角变换理论在流体力学等许多领域中有着广泛的应用。但计算保角变换是个很困难的问题,因此寻求一种有效方法计算保角变换在实际应用中具有很大意义。本论文提出了一种数值保角变换的新算法,在这个新算法中我们在改进高斯消去法的基础上利用模拟电荷法计算新的电荷点,进而构造高精度的近似保角变换函数,并且通过典型图形的数值实验检验了新算法的有效性。

第十六保角变换法求解定解问题共37页文档

第十六保角变换法求解定解问题共37页文档
w f (z),可以将它转化为wuiv平面上
(u ,v)的拉普拉斯方程边值问题.
w = 同理可以证明,在单叶解析函数 f (z)
变换下,泊松方程
22(x,y)
x2 y2
(16.1.7a)
仍然变为泊松方程
2 u 2+ 2 v 2 |f(z)|2(x,y) (16.1.7b)
由上式可知,在保角变换下,泊松方程中的电荷密度
发生了变化.
同理可以证明,亥姆霍兹方程
2x2 2y2 k20 (16.1.8a)
经变换后仍然变为亥姆霍兹方程
22k2|f(z)|20 (16.1.8b) u2 v2
容易注意到方程要比原先复杂,且
能不是常系数.
前的系数可
下面将举例说明如何通过保角变换法来求解拉普拉斯方程.
保角变换法的优点不仅在于拉普拉斯方程、泊松方程 等方程的类型在保角变换下保持不变,更重要的是,能将 复杂边界问题变为简单边界问题,从而使问题得到解决.
(16.2.1)
作如下的保角变换.
(1)作分式线性变换
1
1
i1
za za
(16.2.2)
y
z 平面
1
1 平面
平面
πi
a
0
1
x
图图181.16.1
可以验证,考虑实轴 zx,(y0)的对应关系:
| (i)若 x | a ,则 axa,故
1
x x
a a
0 ,即有
1
0
(ii)若 | x | a 则 xa 或 xa
问题中的解析法――保角变换法,它是解决这类复杂边 界的最有效方法.它特别适合于分析平面场的问题,
例如静电场的问题,由于这种求解复杂边界的定解问 题具有较大的实用价值,所以有必要单独以一章的内 容进行介绍.复变函数论中已经系统介绍了保角变换

《流体力学》课件 3.9 保角变换

《流体力学》课件 3.9 保角变换

d w dW d d z d d z
在无穷远处,有:
d w d z
dW d
d dz
考虑到
dW d
kV

d dz
1 k
,有:
dw dz
V
三、环量的确定
1. 补充条件
dw 有限的常数
dz zB 2. 环量的确定
dz
d
E
0
dw 有限常数
dz zB
dw
d
E
w1 z
Q
2
lnz
i
h
Q
2
lnz
i
h
Q ln z2 h2 2
wz
w1 z
w1
a2 z
wz
Q
2
ln z2
a4 z2
h2
a4 h2
dz
d
k
;(其中:
k
是正的实数)
(根据黎曼定理这样的函数存在且是唯一的)
W
kV
w
z
kV R
kV
2
2 i
F z
ln
kV
F z
R
2
F z
ln
2 i
F
z
证明:1. 因W 是在 K D 上连续且在 D 内解析的函数, Fz是在 C D 上连续且在 D 内解析的函数。于是,根据复合函数的性质 wz W F z
一、保角变换的概念
w f z
V f lin w f ei Δz0 z
w wei f eiz f z ei
12
1 2
2 1` 2 1`
黎曼定理:
任何一个单连通区域必可通过某个保角变换 变为另一个任意给定的单连通区域。

144《高等渗流力学》—保角变换及应用

144《高等渗流力学》—保角变换及应用
四、保角变换及应用
定义:平面 z = x + iy 上给出某个流动其复势是F ( z ) ,引 入新复变函数 w = u + iv
两者间关系: z = z ( w) 或则 w = w( z )
z 实部和虚部关系: ( w) = z (u + iv ) = x (u , v) + iy (u , v ) x = x (u , v ), y = y (u , v ) …………….....(1)
⎛ R1 ⎞ Re > R1 , m ≥ 5 时 ⎜ ⎟ ⎝ Re ⎠
<< 1
2π ( Φ e − Φ w ) q= R R m ln e + ln 1 R1 mrw
14
保角变换及应用
例五:
取变换: w = 直线无限井列的变换。
ρe ⋅ e
iπ z a

w = ρe ⋅ e
θ= πx
a
iπ ( x + iy ) a
上式是长轴为 a ,短轴为b的椭圆方程,给定一个 ρ 值,z平面上 给定一条等势线(圆)。因此,上式为z平面等势线方程。另 外,由于 a 2 − b 2 = c 2 ,故Z平面上所有等式椭圆共焦,焦距为c。 x2 y2 a 2 − b2 − 2 2 = =1 2 2 2 c cos θ c sin θ c W平面上给定一v值,相当于给定一流线,故上式为Z平面的 流线方程。
dL
dφ vn = − —— 法线渗流速度; dn
dz dn = dv dw
5
保角变换及应用
dz dL = dλ 和 dw
例一:

dφ dL = dn

dφ dz dφ dλ = ∫ dλ dz dw dv dv dw

保角变换数学物理方法

保角变换数学物理方法
映射后跟 C1与C2对应的曲线 1与 2之间的夹角. 映射 w f (z) 具有保持两曲线间夹角的大小和 方向不变的性质, 此性质称为保角性.
2
第3页/共31页
4)伸缩率
极限
f
(z0 )
lim
zz0
s
(s表示C上点z0与z间的
弧长, 表示上C在z0的伸缩率.
16
第17页/共31页
解2 利用不变对称点
因 z 1 i 时, w , 所以 w az b , z (1 i)
又 z 1时, w 1, 故 i a b,
由对称点的不变性知, z 1 对应 w 0, 1 i
故 b 1, a 1 i,
所以 w (1 i)z 1 (i 1)z 1 为所求. z (1 i) z (1 i)
1 i 据分式线性映射不变对称点的性质知
15
第16页/共31页
w 0在z平面上的逆象为z 1 (z 1 i对应w ). 1 i
由交比不变性知
(1,0,
w,)
1,
1
1
i
,
z,1
i

w w
1
z
z
1 1
1i
1 i 1 1i 1
z1 , zi z i
1 i
所以 w (i 1)z 1 为所求. z (1 i)
定义 设 w f (z)在 z0 的邻域内是解析的,在 z0 具有保角性和伸缩率不变性,那末 w f (z) 在 z0 是共形的,或称w f (z) 在 z0 是共形映射. 也称为第一类共形映射.仅保持夹角的绝对值不 变而方向相反的映射, 称为第二类共形映射
4
第5页/共31页
3.分式线性映射
10
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

xa xa
0
x a (b)再考虑
的情况, 则
x a 0 ,x a 2 a 0 ,

1
xa xa
0
如图16.1所示,根据(16.2.1)式中的边界条件,对应于
|
x
|
a
u 处温度为
,故
0
1
平面的负实轴(即 1
0

u 温度保持为 0 ;而在| x | a 处有 1 0 ,故
1 平面的正实轴温度保持为零.
是下列定解问题的解
(3)可证 (x ,y )1 (x ,y )2 (x ,y )
是下列定解问题的解:
V0 (x1)
0(y0,x),(x,0) V1 (1x1)
V2 (1x)
其中
V 0 V 1 V 1 , V 1 V 2 V 1 , V 2 V 2 V 2
( x, y)又可改写成
u 处保持温度 = 0.求平板上的稳定温度分布.
【解】根据题意可得出定解问题
2u
x
2
2u y 2
0
u
y0
u 0 , ( x a ) 0 , ( x a )
(16.2.1)
作如下的保角变换.
(1)作分式线性变换
1
1
i1
za za
(16.2.2)
y
z 平面
1
1 平面
平面
πi
a
0
第十六章 保角变换法求解定解问题
在许多物理问题中(如电学、热学、光学、流体 力学和弹性力学等)经常会遇到解平面场的拉普拉斯 方程或泊松方程的问题.尽管可用前几章的理论方法 如:分离变量法或格林函数法等来解决,但当边值问 题中的边界形状变得十分复杂时,分离变量法和格林 函数法却显得十分困难,甚至不能解决.对于复杂的
(16.1.7b)
由上式可知,在保角变换下,泊松方程中的电荷密度
发生了变化.
同理可以证明,亥姆霍兹方程
2x2 2y2 k20 (16.1.8a)
经变换后仍然变为亥姆霍兹方程
22k2|f(z)|20 (16.1.8b) u2 v2
容易注意到方程要比原先复杂,且
能不是常系数.
前的系数可
下面将举例说明如何通过保角变换法来求解拉普拉斯方程.
u v, v u x y x y
(16.1.4)
以及解析函数的实部和虚部分别满足拉普拉斯方程的性质
x2u2 y2u2 0,
x2v2 y2v2 0(16.1.5)
将式(16.1.4)和式(16.1.5)代入到式(16.1.3)化简后得到
2 2 [ ( u ) 2 ( v ) 2 ] ( 2 + 2 ) |f(z )|2 ( 2 + 2 )
理论,本章主要介绍利用保角变换法求解定解问 题。
16.1 保角变换与拉普拉斯方程边值问题的关系
w 在复变函数论中我们已经知道,由解析函数 f (z)
实现的从z平面到 w 平面的变换在 f (z) 0 的点具有保
角性质,因此这种变换称为保角变换.下面我们主要讨论一一
对应的保角变换,即假定 w f (z)和它的反函数都是单值
Im π ,如图16.3
所以,在变换 ln a z i之下,定解问题变换为
az
v
2v 0 0 v 1
v π v 2
定解问题的解(仿上例16.2.1)为
vv2π v1v1v2π v1Im
z 将变量回到 平面,则
v
v1
v2
v1 π
Im[ln(azi)] az
v1
v2
v1 π
1
x
图图181.16.1
可以验证,考虑实轴 zx,(y0)的对应关系:
| (i)若 x | a 0
,即有
1
0
(ii)若 | x | a 则 xa 或 xa
xa (a)首先讨论
a 的情况,考虑到题给条件 0
则 x a 0 ,x a 2 a 0 ;

1
Im{ln[yi(ax)]ln[(xa)iy]
v1
v2
v1 π
[arctanayxarctanxya]v1
v2
v1 π
a2 arctan
x2 y2 2ay
v1
v2
v1 π
π
[ arctan
2
a2
2ay x2 y2
]v1
v2 2
v1
v2 π
arctan a2
2 ay x2 y2
化成极坐标形式,则上式又改写成
d2u
d
2
0
u 0 0, u π u 0
(16.2.6)
方程的解是 uAB
考虑边界条件即得到
z 回到 平面,则
u u 0 (16.2.7) π
u(x,y)u0Im[ln(za)]u0Im[ln(za)ln(za)] π za π
u0[arctg y arctg y ]u0arctg 2ay π xa xa π x2y2a2
边界形状,拉普拉斯方程定解问题常采用保角变换 法求解.
保角变换法解定解问题的基本思想是:通过解析
函数的变换(或映射,这部分知识在复变函数论中已经学
z w 习过)将
平面上具有复杂边界形状的边值问题变换为
平面上具有简单形状(通常是圆、上半平面或带形域)的
边值问题,而后一问题的解易于求得.于是再通过逆变换
(x ,y ) V 0 π V 1 a r g (z 1 ) V 1 π V 2a r g (z 1 ) V 2
(4)进一步推广
( x ,y ) V 0 π V 1 a r g ( z x 1 ) V 1 π V 2 a r g ( z x 2 ) L V n 1 π V n a r g ( z x n ) V n
2(v)2 2 2 u v (16.1.2)
v2 y
uv y y
两式相加得到
22[(u)2 (u)2]2+[(v)2 (v)2]2
x2 y2 x
y u2 x
y v2
2u +(x2
2u y2 )
u (x2v2
y2v2 )
v
+2( u v + u v ) 2 x x y y uv
(16.1.3)
利用解析函数 wf(z)uiv的C-R条件
(2)作变换 ln1 ln| 1| ia rg1
(16.2.3)
把 1 平面的上半平面变成 平面上平行于实轴,宽为
π 的一个带形区域, 1 平面的正实轴变换为
0 平面的实轴(正实轴辐角为零,故对应于
),
1 平面的负实轴变换为 平面的平行于实轴的直线
π (负实轴辐角为 ,故对应于 = π ).
函数;或者如果它们之中有多值函数就规定取它的黎曼面的一 叶.
定律16.1.1 如果将由z x iy 到 wuiv
的保角变换看成为二元(实变)函数 ( x, y ) 的变换由 x , y
u , v z w 到
的变量代换,则 平面上的边界变成了
平面上的边界.我们能证明,如果 ( x, y ) 满足拉普拉斯方
(u,v) V1 V2 v (16.2.10)
π
v w平面
y z平面
πi v1
v2
O
v1
v2
x
u
图图118.2 6.2
O

w u ivlnzlnz iargz
所以 vargz
于是,作反变换便可求得所求问题的解为
(x ,y ) V 1 V 2 a r g z V 1 V 2 V 1 V 2 a r c ta n y
z ln z (3)再作变换
2

平面的上半平面变成
2
π 平面上平行于实轴,宽为
的一个带形区域,其边界的
y
z 平面
0 x
y1 z1 平面
x1
y2
z2 平面 x2
图图118.63 .3
πi
v2
平面
v1
z 变换是将 2 平面的正半实轴变换为 平面的实轴,
z 2 平面的负半实轴变换为 平面的平行于实轴的直线
就求得了原始定解问题的解.
这就是本章将要介绍的一种解决数学物理方程定解
问题中的解析法――保角变换法,它是解决这类复杂边 界的最有效方法.它特别适合于分析平面场的问题,
例如静电场的问题,由于这种求解复杂边界的定解问 题具有较大的实用价值,所以有必要单独以一章的内 容进行介绍.复变函数论中已经系统介绍了保角变换
x 2 y 2 x x u 2 v 2
u 2 v 2
注意到上式已经使用了:
wf(z)uiv x x
对于保角变换 wf(z)0, 因而只要
( x, y ) 满足拉普拉斯方程,则 (u , v )也满足拉
普拉斯方程,即为
22
22
0 x2 y2
(u2+v2)0(16.1.6)
这样我们就有结论:如果在 z x iy 平面上给定了
程,则经过保角变换后得到的 (u,v )也满足拉普拉斯方程.
【证明】 利用复合函数求导法则有
u v x u x v x
2 2u 2 (u )2 2v x2 u x2 u2 x v x2
2 ( v )2 2 2 u v
v 2 x
uv x x
(16.1.1)
同理
2 2u 2(u)2 2v y2 u y2 u2 y v y2
内的电场分布情况.
【解】即求解定解问题
2v 0
v
a
v1 v2
( a) (0 π) (π 2π)
作如下的保角变换
(1) 作变换
z1
z a
相关文档
最新文档