ラプラス方程式(静電場)
ラプラス方程式(静電場)の理論基礎
ラプラス方程式
先生、ラプラス方程式はポアソン方程式とどう違いますか?
電荷がない領域($\rho_v = 0$)ではポアソン方程式がラプラス方程式になる:
電極間の絶縁体内部、誘電体内部など、自由電荷がない空間の電位分布はこれで決まる。
高電圧機器の絶縁設計は大体ラプラス方程式ですか?
そう。電極の電位(Dirichlet BC)だけで電位分布が決まる。絶縁体の電界分布、等電位線の計算はラプラス方程式の領域。
最大値原理
ラプラス方程式の重要な性質: 電位の最大値・最小値は必ず境界上にある(領域内部にはない)。
内部に極値がないということは、FEMの結果で領域内部に電位のピークがあったらおかしいんですね。
そう。結果の妥当性検証に使える。もし内部にピークがあれば、メッシュの問題か電荷源の設定ミス。
まとめ
ラプラス方程式の解には「最大値・最小値がない」——実務への影響
ラプラス方程式の解が持つ「最大値原理」は、CAE実務で地味に重要な性質です。電荷が存在しない空間( $\nabla^2\phi = 0$ が成立する領域)では、電位の最大値と最小値は必ず境界上に現れます。つまり内部空間に「電位の極大点」が突然現れることはないんです。これが絶縁設計に使えて、「最も電界が強い場所=境界(電極や誘電体界面)上」と確定できるので、内部全域を細かくチェックしなくても境界付近だけ精度を上げればいい。高圧機器設計では、この原理を活用してメッシュの集中箇所を絞り込む工夫が現場でよく使われます。
ラプラス方程式(静電場)の数値計算手法
FEM解法
ラプラス方程式のFEMはポアソン方程式の右辺をゼロにしたもの。$[K]\{\phi\} = \{0\}$(電荷なし)+ 境界条件。
解法は完全に同じ。違いは右辺ベクトルが電極のDirichlet条件由来の項のみになること。
解析解との比較
| 問題 | 解析解 | 用途 |
|---|---|---|
| 平行平板 | $\phi = V_0 (1 - x/d)$ | 最も簡単な検証 |
| 同軸円筒 | $\phi = V_0 \ln(r/b)/\ln(a/b)$ | ケーブルの検証 |
| 同心球 | $\phi = V_0 (1/r - 1/b)/(1/a - 1/b)$ | 3D球対称の検証 |
まとめ
「ラプラス方程式を流体でも解く」——調和関数の普遍性
ラプラス方程式 $\nabla^2\phi = 0$ は静電位だけでなく、非圧縮ポテンシャル流の速度ポテンシャル、定常熱伝導の温度分布、多孔質媒体中の圧力場など、驚くほど多様な物理現象を支配します。これは「境界条件が違っても、同じソルバーコードで解ける」ことを意味します。実際に静電界の解析コードを書いた人が「誘電率を熱伝導率に置き換えたら熱解析になった」という経験をすることがあります。この数学的な共通性を「アナロジー(類比)」と呼び、昔は電気的なアナログ回路を使って流体の圧力場を測定する実験が行われていました。ラプラス方程式は自然界の「共通言語」なんです。
ラプラス方程式(静電場)の実務適用
実務
チェックリスト
高圧変電所のシールドリング——ラプラス方程式設計の実例
高圧変電所の機器端子によく付いている「シールドリング(コロナリング)」は、ラプラス方程式の解を活用した設計の実例です。端子の角やエッジに電界が集中すると、コロナ放電が発生して絶縁劣化や電波障害の原因になります。そこで滑らかな金属リングを端子付近に設置し、等電位面を丸めることで電界の集中を緩和します。設計者はラプラス方程式の静電界解析でリングのサイズ・位置・形状を最適化します。「電界が均一になる形とは何か」——これを数値で見つけるのは、まさにCAEが最も得意とする仕事です。
ラプラス方程式(静電場)のソフトウェア比較
ツール
全静電場ソルバーで対応。ポアソン方程式と同じツール。
ラプラス方程式専用ツール vs 汎用EM解析ツールの選び方
純粋な静電界(ラプラス・ポアソン方程式)だけを解くなら、汎用電磁界解析ソフトの全機能は不要で、静電界専用の軽量ソルバーが使えます。SIMION(荷電粒子軌跡解析向け)やElmer(オープンソース)は静電界に特化しておりコストパフォーマンスが高い。一方で「高圧機器の絶縁設計→電界集中→コロナ放電→プラズマへの遷移」まで連続して追いたい場合は、CST StudioやANSYS Maxwellのような汎用ツールが必要です。「今の仕事にラプラス方程式だけで十分か」を最初に見極めることで、ツールのコストと計算時間を大幅に節約できます。
ラプラス方程式(静電場)の先端研究
先端
ラプラス方程式と「等角写像」——複素関数が絵を描く
2Dのラプラス方程式には「等角写像(共形写像)」という解析的な解法があります。複素関数 $w = f(z)$ を使って、複雑な形状の境界値問題を単純な形(例えば円の内側)に変換して解く手法です。コンピュータがなかった時代、変電所の碍子断面の電界計算はこの手法で手計算されていました。現代でもFEMの解析結果を検証したいとき、単純な形状なら等角写像の解析解と比較するのが手軽です。複素関数論が「電気工学の実務ツール」として使われていたという事実は、数学が苦手だったエンジニアには少し驚きかもしれません。
ラプラス方程式(静電場)のトラブル対応
トラブル
「等電位線が鋭角になっている」——メッシュ品質を疑え
ラプラス方程式の解析結果で電位コンター図を描いたとき、「等電位線が急に曲がる」「くびれがある」という箇所が出たら要注意です。通常、ラプラス方程式の解は滑らかな等電位面を持つはずなので、不自然な形状が見えたらメッシュの粗さや品質の問題であることが多い。特に三角形要素で「細長い要素(アスペクト比が大きい要素)」が集中している部分に誤差が乗りやすい。実務では等電位線の滑らかさを「目視でのメッシュ品質チェック」の代わりに使うエンジニアもいます。コンター図は計算結果を見る道具であり、同時に解析品質を確認する道具でもあります。
なった
詳しく
報告