電磁波伝搬
理論と物理
マクスウェル方程式と波動方程式
先生、電磁波の支配方程式を教えてください。
マクスウェル方程式から導かれる波動方程式：
位相速度$v_p = 1/\sqrt{\mu\varepsilon}$。真空中で光速$c = 3 \times 10^8$ m/s。
誘電体中では速度が遅くなるんですね。
比誘電率$\varepsilon_r$の媒質中では$v_p = c/\sqrt{\varepsilon_r}$。FR-4基板（$\varepsilon_r \approx 4.4$）では約$0.48c$。波長も短くなるため、高周波回路ではサイズが問題になる。
まとめ
マクスウェルが1865年に予言した電磁波——ヘルツが証明するまで23年
マクスウェルが方程式から電磁波の存在を予言したのは1865年。しかし「数式上の存在」に懐疑的な物理学者も多く、実際にハインリヒ・ヘルツが実験室で電磁波の送受信を実証したのは1888年、実に23年後のことです。ヘルツが使ったのはスパーク放電による数百MHzの電波。今日の高周波CAEが扱う数GHz〜数百GHzの電磁波理論は、この二人の巨人の仕事の上に成り立っています。マクスウェル方程式を眺めるたびに、23年という時間の重みを思い出してみてください。
各項の物理的意味
- 電場項 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：ファラデーの電磁誘導法則。時間変動する磁束密度が起電力を生じさせる。【日常の例】自転車のダイナモ（発電機）は、磁石を回転させることで近くのコイルに電圧が発生する——磁場が時間的に変化すると電場が誘起されるというこの法則の直接的応用。IHクッキングヒーターも同じ原理で、高周波磁場の変化が鍋底に渦電流を誘起し、ジュール熱で加熱する。
- 磁場項 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：アンペア-マクスウェルの法則。電流と変位電流が磁場を生成する。【日常の例】電線に電流を流すと周囲に磁場が生じる——これがアンペアの法則。電磁石はこの原理で動作し、コイルに電流を流して強力な磁場を作る。スマートフォンのスピーカーも、電流→磁場→振動板の力というこの法則の応用。高周波（GHz帯のアンテナ等）では変位電流 $\partial D/\partial t$ が無視できなくなり、電磁波の放射を記述する。
- ガウスの法則 $\nabla \cdot \mathbf{D} = \rho_v$：電荷が電束の発散源であることを示す。【日常の例】下敷きで髪の毛をこすると静電気で髪が逆立つ——帯電した下敷き（電荷）から電気力線が放射状に広がり、軽い髪の毛に力を及ぼす。コンデンサ（キャパシタ）の設計では、電極間の電場分布をこの法則で計算する。ESD（静電気放電）対策もガウスの法則に基づく電場解析が基盤。
- 磁束保存 $\nabla \cdot \mathbf{B} = 0$：磁気単極子が存在しないことを表す。【日常の例】棒磁石を半分に割っても、N極だけ・S極だけの磁石は作れない——必ずN極とS極がペアで存在する。これは磁力線が「始点も終点もない閉じたループ」を描くことを意味する。数値解析では、この条件を満たすためにベクトルポテンシャル $\mathbf{B} = \nabla \times \mathbf{A}$ という定式化を用い、磁束保存を自動的に保証する。
仮定条件と適用限界
- 線形材料仮定：透磁率・誘電率が磁場・電場強度に依存しない（飽和領域では非線形B-Hカーブが必要）
- 準静的近似（低周波）：変位電流項を無視可能（$\omega \varepsilon \ll \sigma$）。渦電流解析で一般的
- 2D仮定（断面解析）：電流方向が一様で、端部効果を無視できる場合に有効
- 等方性仮定：異方性材料（珪素鋼板の圧延方向等）では方向別の特性定義が必要
- 適用外ケース：プラズマ（電離気体）、超伝導体、非線形光学材料では追加の構成則が必要
数値解法と実装
数値解法
電磁波をどのような数値手法で解きますか？
| 手法 | 定式化 | 得意な問題 |
|---|---|---|
| FEM | 周波数領域 | 共振器、導波路、複雑形状 |
| FDTD | 時間領域 | 広帯域、過渡応答、大規模 |
| MoM | 積分方程式 | 開放空間、アンテナ |
| FIT | 積分形Maxwell | CST Studio Suiteの基盤 |
HFSSはFEM、CST Studio SuiteはFIT/FDTD、FEKO（Altair）はMoMが基盤。
メッシュサイズはどう決めますか？
波長$\lambda$の1/10以下が目安。FEMで2次要素なら$\lambda/5$でも良い精度。FDTDではCFL条件$\Delta t \leq \Delta x/(c\sqrt{3})$を満たす必要がある。
まとめ
レイトレーシング法と携帯基地局の電波設計
都市部の携帯電話基地局設計では、建物・地形による電磁波の反射・回折・散乱を考慮した電波伝搬シミュレーションが不可欠です。フルウェーブFDTDは精度は高いが数百m四方の都市モデルには非現実的なほど計算コストが大きい。そこで活躍するのがレイトレーシング（Ray Tracing）法で、光線追跡のアナロジーで各経路の伝搬損失・遅延を高速に計算します。5G基地局の最適配置を決める「エリア設計」に現在も使われており、建物モデルの精度がシミュレーション精度を大きく左右します。
辺要素（Nedelec要素）
電磁場解析に特化した要素。接線成分の連続性を自動的に保証し、スプリアスモードを排除。3D高周波解析の標準。
節点要素
スカラーポテンシャル定式化に使用。静磁場のスカラーポテンシャル法や静電場解析で有効。
FEM vs BEM（境界要素法）
FEM: 非線形材料・非均質媒質に対応。BEM: 無限領域（開領域問題）を自然に扱える。ハイブリッドFEM-BEMも有効。
非線形収束（磁気飽和）
B-Hカーブの非線形性をニュートン・ラフソン法で処理。残差基準: $||R||/||R_0|| < 10^{-4}$が一般的。
周波数領域解析
時間高調波仮定により定常問題に帰着。複素数演算が必要だが、広帯域特性は時間領域解析で取得。
時間領域の時間刻み
最高周波数成分の1/20以下の時間刻みが必要。暗黙的時間積分ではより大きな刻みも可能だが精度に注意。
周波数領域と時間領域の使い分け
周波数領域解析は「ラジオの特定の周波数に合わせる」ようなもの——1つの周波数での応答を効率的に計算できる。時間領域解析は「全チャンネルを同時に録画する」ようなもの——あらゆる周波数成分を含む過渡現象を再現できるが計算コストが高い。
実践ガイド
実務での適用
アンテナ設計、導波路設計、レーダ断面積（RCS）計算、電波伝搬予測が代表的。
実務チェックリスト
伝搬損失の測定——フリスの伝達方程式と現実のギャップ
自由空間でのアンテナ間の受信電力はフリス（Friis）の伝達方程式で計算できますが、実際の現場ではこの理論値より10〜40 dBも損失が大きくなるのが普通です。原因はマルチパス干渉・障害物による遮蔽・地表反射など。「現場での受信電力 vs シミュレーション値」を比較するキャリブレーション実測は避けられず、測定ポイントの取り方・アンテナ設置高度・周辺環境のモデル化が精度を左右します。実践ではまずフリス計算を上限として把握し、環境係数（パスロスエクスポーネント）で補正するのが定石です。
解析フローのたとえ
モータの電磁界解析は「ギターの調律」に近い感覚です。弦の太さ（コイル巻数）とブリッジの位置（磁石配置）を調整して、最も美しい音色（効率の良いトルク特性）を引き出す。1つのパラメータを変えると全体のバランスが変わる——だからパラメトリックスタディが重要なんです。
初心者が陥りやすい落とし穴
「空気領域？ なんで空気をメッシュで切るの？」——初めて電磁界解析に触れた人がほぼ全員抱く疑問です。答えは「磁力線は鉄心の外にも広がるから」。解析領域を鉄心ぎりぎりにすると、行き場を失った磁束が壁に「ぶつかって」反射し、実際にはありえない磁束集中が起きます。部屋が狭すぎてボールが壁に跳ね返りまくる状態を想像してみてください。
境界条件の考え方
遠方の境界条件って地味ですが超重要です。「ここから先は無限に広がる空間」ということを数値的に表現する必要がある。設定を間違えると、まるで「見えない壁」があるかのように磁束が跳ね返されてしまいます。
ソフトウェア比較
ツール
| ツール | 特徴 |
|---|---|
| Ansys HFSS | FEMベース。適応メッシュ。複雑形状に強い |
| CST Studio Suite | FIT/FDTD。広帯域。Sパラメータ抽出 |
| COMSOL RF | マルチフィジックス連成。熱-電磁連成 |
| Altair FEKO | MoM+MLFMM。大規模アンテナ・RCS |
電波伝搬シミュレーターの使い分け——フルウェーブ vs. 統計モデル
電磁波伝搬解析ツールは用途によって使い分けが大切です。ANSYS HFSSやCSTなどのフルウェーブシミュレータは数cm〜数mスケールのアンテナ近傍解析に向いています。一方、都市エリアスケール（数百m〜数km）の電波伝搬にはWinProp（Altair）やSiradel Volcano、あるいはITU-R推奨の統計的パスロスモデルが使われます。通信システム設計ではこの2段階を組み合わせることが多く、「近傍は電磁界シミュレーション、遠方は統計モデル」というワークフローが実務標準です。
選定で最も重要な3つの問い
- 「何を解くか」：電磁波伝搬に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術
大気の窓——電磁波伝搬と天文観測の関係
電磁波は大気中を伝搬する際、周波数によって吸収率が大きく異なります。水蒸気や酸素分子が特定の周波数帯（22 GHz、60 GHz、183 GHz付近など）を強く吸収するため、「大気の窓」と呼ばれる吸収の少ない帯域が存在します。電波天文学の望遠鏡はこの窓を利用して宇宙からの電波を観測し、5G mmWave通信の周波数選定もこの大気吸収特性を考慮しています。電磁波伝搬シミュレーションで大気の誘電体モデルを正確に組み込むことが、地上・衛星通信システム設計の精度に直結します。
トラブルシューティング
トラブル
シミュレーションが現場と合わない——吸収境界の設定ミスが元凶
電磁波伝搬シミュレーションで「計算領域の端で電波が反射して結果がおかしい」というトラブルは初心者に非常によくあります。原因は吸収境界条件（ABC）またはPML（完全整合層）の設定不足。計算ドメインの端が導体壁扱いになってしまい、本来は外に逃げるはずの電波がドメイン内で多重反射してしまう。対処は①PML層を厚くする（通常は波長の1/4以上）、②ドメインサイズをアンテナから十分に離す（最低でも半波長以上）、③反射係数を確認するテスト解析を実施する、の3ステップです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——電磁波伝搬の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告