5.8 RWG基底関数を用いたポテンシャル積分について
RWG基底関数を用いたベクトルポテンシャルおよびスカラポテンシャルにおいて,観測点と波源が一致する特異点を含めた積分について説明する† .
† D. R. Wilton, S. M. Rao, A. W. Glisson, D. H. Schaubert, O. M. Al-Bundak, and C. M. Butler,
“Potential integrals for uniform and linear source distributions on polygonal and polyhedral domains,”
IEEE Trans. Antennas Propagat. , vol.32, no.3, pp.276-281 (1984).
三角形領域T q おける電流3成分( i = 1 , 2 , 3 ) から領域T p の中心点(位置ベクトルr p ( c ) )への寄与をまとめて考えると,
これらの積分は次のようになる.
(1) A p q , i = μ 4 π ∬ T q ( l q , i 2 A q ρ q , i ) e − j k R p R p d S ′ ( i = 1 , 2 , 3 ) (2) Φ p q , i = ∓ 1 4 π j ω ϵ ∬ T q ( l q , i A q ) e − j k R p R p d S ′ ( i = 1 , 2 , 3 )
ここで,
(3) R p = | r p ( c ) − r ′ | (4) ρ q , i = ± ( r ′ − v q , i ) ( i = 1 , 2 , 3 )
ただし,v q , i およびl q , i ( i = 1 , 2 , 3 ) は,三角形領域T q の頂点の位置ベクトルおよびその対辺の長さ,
A q はその面積を示す.また,上式の± の符号は,領域T q がT ( ± ) のどちらかで対応して決まる.
各辺i (i = 1 , 2 , 3 )に関する計算に用いる変数の定義
簡単のため,i = 1 を考え,観測点を一般化してr とおく.また,T q の符号は( + ) について考え,各変数の添字を略すとベクトルポテンシャルA は,
(5) A = μ 4 π l 2 A ∬ T ρ m ( + ) e − j k R R d S ′
上式の積分項 I を,次のような項に分けて計算していく.
I = ∬ T ρ m ( + ) e − j k R R d S ′ = ∬ T ( ρ ′ − ρ v ) e − j k R R d S ′ (6) = ∬ T ( ρ ′ − ρ ) e − j k R R d S ′ + ( ρ − ρ v ) ∬ T e − j k R R d S ′
ただし,ρ ,ρ ′ ,ρ v は,
観測点の位置ベクトルr ,電流源の位置ベクトルr ′ ,
共有する辺に対向する三角形の頂点v m ( + ) を各々面S に投影したベクトルを示す.
RWG基底関数に関わる変数の定義
さらに,
(7) e − j k R R = ( e − j k R R − 1 R ) + 1 R (8) lim R → 0 ( e − j k R R − 1 R ) = − j k
を考慮して,
I = ∬ T ( ρ ′ − ρ ) e − j k R − 1 R d S ′ + ∬ T ρ ′ − ρ R d S ′ (9) + ( ρ − ρ v ) ∬ T e − j k R − 1 R d S ′ + ( ρ − ρ v ) ∬ T d S ′ R
上式の第1項および第3項は特異点を持たないため,容易に数値積分できる.
3角形領域の積分について
まず,三角形領域T q における波源の位置ベクトルr ′ は,
(10) r ′ = v q , 1 + ξ q ( v q , 2 − v q , 1 ) + η q ( v q , 3 − v q , 1 ) ( 0 ≤ ξ q , η q ≤ 1 )
簡単のため,添字q は省略して,
r ′ = v 1 + ξ ( v 2 − v 1 ) + η ( v 3 − v 1 ) ( 0 ≤ ξ , η ≤ 1 ) (11) = ( 1 − ξ − η ) v 1 + ξ v 2 + η v 3
これより,面積要素d S ′ は,
d S ′ = | ∂ r ′ ∂ ξ × ∂ r ′ ∂ η | d ξ d η = | ( v 2 − v 1 ) × ( v 3 − v 1 ) | d ξ d η (12) = 2 A q d ξ d η
よって,三角形領域T q の積分は,次のような2重積分の形で表される.
(13) ∬ T q g ( r ′ ) d S ′ = 2 A q ∫ ξ = 0 1 ∫ η = 0 1 − ξ g ( ( 1 − ξ − η ) v 1 + ξ v 2 + η v 3 ) d η d ξ
3角形領域の面積分に関わる変数の定義
特異点を含む積分について
式(9) の第2項および第4項の積分I 2 ,I 4
(14) I 2 = ∬ T ρ ′ − ρ R d S ′ , I 4 = ∬ T d S ′ R
については,周回積分への変換と特異点のふるまいを考慮した積分を行う.まず,準備として,導体面に沿う2次元演算子
∇ s
を作用させたときの計算を行う.
(15) P = P u p ≡ ρ ′ − ρ (16) P = | ρ ′ − ρ | (17) R = | r − r ′ | = | d n + ρ − ρ ′ | = | d n − P u p | = d 2 + P 2 (18) ∇ s ′ R = ∂ R ∂ P u p = P R u p = P R = ρ ′ − ρ R
これより,積分I 2 に対して,特異点P = 0 の点を含む微小領域をs ϵ として,次のように領域s ϵ とそれ以外に分けて積分し,それから極限を求める.
I 2 = ∬ S ρ ′ − ρ R d S ′ (19) = lim s ϵ → 0 ∬ S − s ϵ ∇ s ′ R d S ′ + lim s ϵ → 0 ∬ s ϵ ρ ′ − ρ R d S ′
また,
(20) R P = d 2 + P 2 P = d 2 P 2 + 1
∂ ∂ P ( R P ) = ∂ ∂ P ( d 2 P 2 + 1 ) (21) = − d 2 P 3 P d 2 + P 2 = − d 2 P 2 1 R
∇ s ′ ⋅ ( R P u ρ ) = ∇ s ′ R P ⋅ u ρ + R P ∇ s ′ ⋅ u ρ = ∂ ∂ P ( R P ) + R P 1 P = − d 2 P 2 1 R + R P 2 (22) = − d 2 + R 2 P 2 R = P 2 P 2 R = 1 R
同様にして,積分I 4 も特異点P = 0 の点を含む微小領域をs ϵ として,
I 4 = ∬ T d S ′ R (23) = lim s ϵ → 0 ∬ S − s ϵ ∇ s ′ ⋅ ( R P P 2 ) d S ′ + lim s ϵ → 0 ∬ s ϵ d S ′ R
と分けて考えることにする.
周回積分への変換
式(23) の第1項はP = 0 の点を除いた積分であり,次のようにガウスの2次元発散定理により周回積分に変換できる.
lim ϵ → 0 ∬ S − s ϵ ∇ s ′ ⋅ ( R P P 2 ) d S ′ = lim ϵ → 0 ∮ C − c ϵ R P P 2 ⋅ d u ′ (24) = ∮ C R P 2 ( P ⋅ u ′ ) d l ′ + lim ϵ → 0 ∮ − c ϵ R P ( u p ⋅ u ′ ) d l ′
ただし,周回積分路C は面S の周辺,
c ϵ は面s ϵ の周辺の各々閉じた経路を示す.また,ベクトル線積分要素
d u ′ ( = u ′ d l ′ ) の方向は,面S ,s ϵ 上で積分路C ,c ϵ の法線方向にとる.上式の第1項は,積分路C が3角形領域の辺i に沿う直線の積分路L i となるので,
(25) ∮ C R P 2 ( P ⋅ u ′ ) d l ′ = ∑ i ( P ( 0 ) ⋅ u i ′ ) ∫ L i R P 2 d l ′
ここで,
(26) P 2 = ( P ( 0 ) ) 2 + l ′ 2 (27) R 2 = d 2 + ( P ( 0 ) ) 2 + l ′ 2 = d 2 + P 2
より,
R P 2 = R P 2 R R = P 2 + d 2 P 2 R = 1 R + d 2 P 2 R (28) = 1 d 2 + ( P ( 0 ) ) 2 + l ′ 2 + d 2 { ( P ( 0 ) ) 2 + l ′ 2 } 2 d 2 + ( P ( 0 ) ) 2 + l ′ 2
辺i (積分路L i )の両端の座標成分をl = l i ( − ) , l i ( + ) とすると,
∫ L i R P 2 d l ′ = ∫ L i ( 1 R + d 2 P 2 R ) d l ′ = ∫ l i ( − ) l i ( + ) d l ′ d 2 + ( P i ( 0 ) ) 2 + l ′ 2 (29) + d 2 ∫ l i ( − ) l i ( + ) d l ′ { ( P i ( 0 ) ) 2 + l ′ 2 } d 2 + ( P i ( 0 ) ) 2 + l ′ 2
辺i に沿う線積分に関わる変数の定義(i = 1 のとき)
上式は,置換積分(t ≡ x + x 2 + a 2 )より得られる不定積分公式(積分定数C は省略)
(30) ∫ d x x 2 + a 2 = ln | x + x 2 + a 2 |
を用い,
(31) R i ( ± ) = d 2 + ( P i ( 0 ) ) 2 + ( l i ( ± ) ) 2
とおいて,式(29) の第1項は,
∫ l i ( − ) l i ( + ) d l ′ d 2 + ( P i ( 0 ) ) 2 + l ′ 2 = [ ln | l ′ + d 2 + ( P i ( 0 ) ) 2 + l ′ 2 | ] l i ( − ) l i ( + ) = ln | l i ( + ) + d 2 + ( P i ( 0 ) ) 2 + ( l i ( + ) ) 2 | − ln | l i ( − ) + d 2 + ( P i ( 0 ) ) 2 + ( l i ( − ) ) 2 | = ln | l i ( + ) + R i ( + ) | − ln | l i ( − ) + R i ( − ) | (32) = ln | l i ( + ) + R i ( + ) l i ( − ) + R i ( − ) |
また,置換積分(| q | x ≡ p x 2 + p 2 + q 2 tan t )より得られる不定積分公式(積分定数C は省略)
(33) ∫ d x ( x 2 + p 2 ) x 2 + p 2 + q 2 = 1 | p q | tan − 1 ( | q | x p x 2 + p 2 + q 2 )
を用いると,式(29) の第2項は,
d 2 ∫ l ′ = l i ( − ) l i ( + ) d l ′ { ( P i ( 0 ) ) 2 + l ′ 2 } d 2 + ( P i ( 0 ) ) 2 + l ′ 2 = d 2 P i ( 0 ) | d | [ tan − 1 ( | d | l ′ P i ( 0 ) d 2 + ( P i ( 0 ) ) 2 + l ′ 2 ) ] l ′ = l i ( − ) l i ( + ) (34) = | d | P i ( 0 ) { tan − 1 ( | d | l i ( + ) P i ( 0 ) R i ( + ) ) − tan − 1 ( | d | l i ( − ) P i ( 0 ) R i ( − ) ) }
積分路− c ϵ においては,u p ⋅ u ′ = − 1 ,d l ′ = P d ϕ ,R = d 2 + ϵ 2 より,
lim ϵ → 0 ∮ − c ϵ R P ( u p ⋅ u ′ ) d l ′ = lim ϵ → 0 ∫ 0 α R P ( − 1 ) P d ϕ = lim ϵ → 0 ( − R ) ∫ 0 α d ϕ = lim ϵ → 0 − d 2 + ϵ 2 ⋅ α (35) = − α | d |
一方,式(23) のI 4 の第2項は,P = 0 の点を座標原点にしたローカルな極座標系を考え,積分範囲s ϵ を半径ϵ ,平面角α [rad]の円形領域にとり,
t ≡ d 2 + P 2 の置換積分を実行すると,次のようになる.
lim ϵ → 0 ∬ s ϵ d S ′ R = lim ϵ → 0 ∫ ϕ = 0 α ∫ P = 0 ϵ P d P d ϕ d 2 + P 2 = lim ϵ → 0 [ ϕ ] 0 α ∫ 1 t d t 2 = lim ϵ → 0 α [ t ] = lim ϵ → 0 α [ d 2 + P 2 ] 0 ϵ (36) = lim ϵ → 0 α ( d 2 + ϵ 2 − | d | ) = 0
よって,
I 4 = ∬ T d S ′ R = − α | d | + ∑ i ( u P 0 , i ⋅ u i ′ ) ⋅ [ P i ( 0 ) ln | l i ( + ) + R i ( + ) l i ( − ) + R i ( − ) | (37) + | d | { tan − 1 ( | d | l i ( + ) P i ( 0 ) R i ( + ) ) − tan − 1 ( | d | l i ( − ) P i ( 0 ) R i ( − ) ) } ]
ただし,
(38) α = ∑ i ( u P 0 , i ⋅ u i ′ ) { tan − 1 ( l i ( + ) P i ( 0 ) ) − tan − 1 ( l i ( − ) P i ( 0 ) ) }
また,
(39) β ≡ tan − 1 l i ( ± ) P i ( 0 ) (40) γ ≡ tan − 1 | d | l i ( ± ) P i ( 0 ) R i ( ± )
とおき,正接の加法定理
(41) tan ( β − γ ) = tan β − tan γ 1 + tan β tan γ
を用いると,
(42) ( l i ( ± ) ) 2 = ( P i ( ± ) ) 2 − ( P i ( 0 ) ) 2 (43) ( P i ( ± ) ) 2 = ( R i ( ± ) ) 2 − | d | 2 (44) ( P i ( 0 ) ) 2 = ( R i ( 0 ) ) 2 − | d | 2
より,次の関係が得られる(導出省略).
tan − 1 ( l i ( ± ) P i ( 0 ) ) − tan − 1 ( | d | l i ( ± ) P i ( 0 ) R i ( ± ) ) (45) = tan − 1 ( P i ( 0 ) l i ( ± ) ( R i ( 0 ) ) 2 + | d | R i ( ± ) )
これより,
I 4 = ∬ T d S ′ R = ∑ i ( u P 0 , i ⋅ u i ′ ) [ P i ( 0 ) ln | l i ( + ) + R i ( + ) l i ( − ) + R i ( − ) | − | d | { tan − 1 ( P i ( 0 ) l i ( + ) ( R i ( 0 ) ) 2 + | d | R i ( + ) ) (46) − tan − 1 ( P i ( 0 ) l i ( − ) ( R i ( 0 ) ) 2 + | d | R i ( − ) ) } ]
一方,積分I 2 についても,式(19) の第1項は2次元勾配定理より,周回積分(積分路C およびc ϵ )に変換でき,第2項はI 4 の計算と同様にしてゼロになる(導出省略).よって,
I 2 = lim s ϵ → 0 ∬ S − s ϵ ∇ s ′ R d S ′ = lim ϵ → 0 ∮ C − c ϵ R u ′ d l ′ = ∑ i u i ′ ∫ L i R d l ′ + lim ϵ → 0 ∫ − c ϵ R u ′ d l ′ (47) = ∑ i u i ′ ∫ l i ( − ) l i ( + ) d 2 + ( P i ( 0 ) ) 2 + l ′ 2 d l ′ + lim ϵ → 0 R ∫ 0 α u ′ ϵ d ϕ
上式の第2項はゼロ,
第1項は置換積分(t ≡ x + x 2 + a 2 )より得られる不定積分公式(積分定数C は省略)
(48) ∫ x 2 + a 2 d x = 1 2 ( x x 2 + a 2 + a 2 ln | x + x 2 + a 2 | )
より,
I 2 = ∑ i u i ′ 1 2 [ l ′ d 2 + ( P i ( 0 ) ) 2 + l ′ 2 + { d 2 + ( P i ( 0 ) ) 2 } ln | l ′ + d 2 + ( P i ( 0 ) ) 2 + l ′ 2 | ] l i ( − ) l i ( + ) = 1 2 ∑ i u i ′ [ { l ( + ) R i ( + ) + ( R i ( 0 ) ) 2 ln | l i ( + ) + R i ( + ) | } − { l ( − ) R i ( − ) + ( R i ( 0 ) ) 2 ln | l i ( − ) + R i ( − ) | } ] (49) = 1 2 ∑ i u i ′ { ( R i ( 0 ) ) 2 ln | l i ( + ) + R i ( + ) l i ( − ) + R i ( − ) | + l i ( + ) R i ( + ) − l i ( − ) R i ( − ) }
前のページに戻る
「電磁波工学特論」へ戻る