5. PowerPointファイル「parabora」

9月27日
パラボラミラーによる
ミリ波ビーム絞り
1m先までのビームの様子
@30,80,140GHz
1
目次
1.パラボラミラーの焦点距離の違いによる反射ビームの絞られ方の違い
2.周波数の違いによる反射ビームの絞られ方の違い
2
1.パラボラミラーの焦点距離の違いによる
反射ビームの絞られ方の違い
概要
まず、以下の状況を考える。
状況
2次元FDTD法を用いて計算する。
まず、ソースがミリ波ビームを出す。
それを、パラボラミラー2枚を用いて絞る。
ソースはホーンアンテナを想定する。
ビームはZ字のように進ませる。
このとき、後段のパラボラミラーの形を焦点距離の異なる数種類に変える。
後段パラボラからのある一定距離離れた場所でのビーム断面を評価する。
3
計算水準の詳細
周波数:30GHz(一定)
金属製
パラボラミラー(上側)
金属製
パラボラミラー(下側)
ソース、ミラーの配置位置
図1が基準の配置位置、形状である。
図1において、ソースから出たビームの軸
は次の通り進む。まず、下側パラボラミラー
反射中心点まで真左に167mm進む。
そこで120度折れて、上側パラボラミラー
反射中心点まで458mm進む。さらにそこで
120度折れて真左に進む。
上側パラボラミラー反射点中心から
上側パラボラミラー焦点までの距離を
500mm、1000mm(図1)、1500mm、
∞mm(平面ミラー) と変化させる。
ソース中心位置は下側パラボラミラーの
焦点の位置に置く。ソースサイズ30mm。
パラボラミラーのサイズは上下とも
直線近似で約220mm。
ソース
図1 ミラーとソースの位置を示す図
白矢印がミリ波ビーム進行方向
計算領域は青色部分
4
計算水準の詳細(続き)
ソースとミラーの位置座標の式
10目盛り=1mmで書いた場合(1次変換して図1にする前)の位置座標の式は
以下の通り
○ソース
(0,1100)-(0,1400)の線分を、(0,1250)を中心に反時計回りに30度回転
させた線分(パラボラミラー(下側)の反射中心(-1443,416)に向いている)
○パラボラミラー(下側)
Y = X^2/5000 (-2443<=X<=-443)
○パラボラミラー(上側)
Y = -(X+5773)^2/15000 + 6250 (-2443<=X<=-443) @f=500mm
Y = -(X+10103)^2/30000 + 7500 (-2443<=X<=-443) @f=1000mm
Y = -(X+14433)^2/45000 + 8750 (-2443<=X<=-443) @f=1500mm
Y = -0.577(X+1443) + 5000 (-2443<=X<=-443)
@f=∞mm
但し、f は焦点距離を表す。
4つの式での反射中心の座標は一定(-1443,5000)となっている。
5
全水準で共通な条件
ソース(ホーンアンテナ)を表す線分上の電界強度は空間的にcosカーブを
維持しながら時間的に以下のようにcosカーブで変化するとした。
時間の流れ
電界強度
3[V/m]
位置
図2
ソース線分上での電界強度の変化
(青と赤で1周期)
6
計算結果の例 @焦点距離1000mm
焦点
1/e幅は128mm
最大強度は1.281
[V/m]
図3
EZ(紙面に垂直方向の
電界強度)片側振幅分布
(目盛りは計算させたときのメッシュ
の細かさになっている)
計算上のビーム軸中心はy=505
図4
図3の白縦線線分上(上側ミラー
から1000mm先)でのEZ片側振幅
プロファイル
7
各水準間でのEZ片側振幅分布(図3に相当)の比較
f:500mm
f:1000mm
[V/m]
f:1500mm
[V/m]
f:∞mm
[V/m]
[V/m]
8
図5
各水準間での1000mm先のビーム断面(図4に相当)の比較
f:500mm
f:1000mm
f:1500mm
f:∞mm
9
図6
最大電界振幅強度
[V/mm]
1000mm先におけるビーム軸中心での電界振
幅強度とミラー焦点距離との関係
注
0.0014
0.0012
0.001
0.0008
0.0006
0.0004
0.0002
0
500
1000
1500
∞
上側パラボラミラーの焦点距離 [mm]
注:軸中心で電界強度が最大ではない場合もありえるがほぼ一致していた
10
電界振幅強度が1/eに
なるまでのビーム幅
[mm]
図7
1000mm先におけるビーム幅と
ミラー焦点距離との関係
200
150
100
50
0
500
1000
1500
∞
上側パラボラミラーの焦点距離[mm]
11
2.周波数の違いによる反射ビームの絞られ方の違い
概要
1のときと同様、以下の状況を考える。
状況
2次元FDTD法を用いて計算する。
まず、ソースがミリ波ビームを出す。
それを、パラボラミラー2枚を用いて絞る。
ソースはホーンアンテナを想定する。
ビームはZ字のように進ませる。
また、後段のパラボラミラーから1000mm先でのビーム断面を見る。
今度は、後段のパラボラミラーの焦点距離を1000mmに固定する。
構成物体の位置は先ほどと同じである。
周波数を3種類振る。30GHz、80GHz、140GHzとする。
12
図8
焦点距離1000mmにおける
各周波数でのビームの軌跡
30GHz
[V/m]
焦点(着目点)
80GHz
140GHz
[V/m]
[V/m]
13
1000mm先におけるビーム軸中心での電界
振幅強度と周波数との関係
最大電界振幅強度
[V/mm]
図9
0.0014
0.0012
0.001
0.0008
0.0006
0.0004
0.0002
0
30GHz
80GHz
140GHz
ソースの送信周波数
14
図10
電界振幅強度が1/eに
なるまでのビーム幅
[mm]
1000mm先におけるビーム幅と
周波数との関係
200
150
100
50
0
30GHz
80GHz
140GHz
ソースの送信周波数
15
補足
FDTD法パラメータ
実際に使用したパラメータは以下のようなものになる
(140GHzでメッシュを比較的小さくとったものの場合)。
(input) X方向の分割数+2の値 =
4657
(input) Y方向の分割数+2の値 =
2532
(input) 全ステップ回数 =
12000
(input) セルのX方向サイズ [m]=
0.00025
(input) セルのY方向サイズ [m]=
0.00025
(input) 時間ステップ[s]=
5e-013
(input) 調べる電磁波の周波数[Hz]=
1.4e+011
1秒当たりのステップ回数 =
2e+012
1周期当たりのステップ回数 =
14.28571
(check) セルのX方向サイズ(=セルのY方向サイズ)/波長=
0.1167474
(check)クーラントパラメータ(<1)=
0.8479411
その他
誘電体:なし
導電体:金属製パラボラミラー、ソース以外はなし
解析領域の周辺部:2次のMurの吸収境界条件を用いた
16