無限線電荷・無限平面・球対称電荷分布の電場をガウスの法則で計算・可視化。ガウス面を動かして閉曲面内の電荷と電場の関係を体験。
線電荷: $E = \dfrac{\lambda}{2\pi\varepsilon_0 r}$
ガウスの法則の積分形です。任意の閉曲面（ガウス面）を考え、その面を垂直に貫く電場の総和（電束）は、面内に存在する全電荷に比例します。
$$\oint_S \mathbf{E}\cdot d\mathbf{A}= \frac{Q_{enc}}{\varepsilon_0}$$$\mathbf{E}$: 電場ベクトル [V/m], $d\mathbf{A}$: 面素ベクトル（大きさは微小面積、方向は面に垂直外向き）[m²], $Q_{enc}$: ガウス面内に閉じ込められた全電荷 [C], $\varepsilon_0$: 真空の誘電率 [F/m]
対称性を利用して電場を求めた具体例（本シミュレーターで可視化）。電荷分布の幾何学的対称性に合わせて適切なガウス面を選択することで、左辺の積分が簡単になります。
$$E(r) \cdot (ガウス面の面積) = \frac{Q_{enc}(r)}{\varepsilon_0}$$無限線電荷: ガウス面を半径$r$の円筒とすると、$E \cdot (2\pi r L) = \frac{\lambda L}{\varepsilon_0}$ ∴ $E = \frac{\lambda}{2\pi\varepsilon_0 r}$
無限平面: ガウス面を平面に垂直な円柱とすると、$E \cdot (2A) = \frac{\sigma A}{\varepsilon_0}$ ∴ $E = \frac{\sigma}{2\varepsilon_0}$
球対称（半径Rの球内一様帯電）: ガウス面を半径$r$の球面とすると、$r \ge R$で $E = \frac{Q}{4\pi\varepsilon_0 r^2}$, $r < R$で $E = \frac{Q r}{4\pi\varepsilon_0 R^3}$
電子部品・コンデンサ設計：平行板コンデンサの極板間電界計算は、無限平面モデルの近似が基礎となります。ガウスの法則から導かれる一定電界を用いて、静電容量や耐圧を設計します。
送電線・高電圧機器：送電線を無限長の線電荷とみなすことで、線路周辺の電界強度分布を推定できます。これにより、他の構造物や人体への影響評価、絶縁設計に活用されます。
半導体デバイス解析：MOSFETなどの半導体内部では、ドーパント分布に応じた電界分布が生じます。ポアソン方程式（ガウスの法則の微分形）を解くことで、デバイスの動作特性をシミュレーションします。
CAEによる静電場解析：ANSYS MaxwellやCOMSOL MultiphysicsなどのCAEソフトウェアは、数値計算（有限要素法）でガウスの法則を解き、複雑な形状のコンデンサ容量、PCBの寄生容量、絶縁破壊が起こりやすい高電界領域の可視化を行います。
このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず、「無限」という言葉の意味だ。「無限線電荷」や「無限平面」は、現実にはありえないけど、観測点から見て十分に長い、あるいは広いと近似できる場合に使える強力なモデルなんだ。例えば、長さ1mの細い電線の、その中心からわずか1cm離れた点の電場を求めるなら、ほぼ「無限線電荷」の式 $E = \frac{\lambda}{2\pi\varepsilon_0 r}$ で近似できる。でも、電線の端っこ近くではこの式は成り立たない。シミュレーターで「ガウス面の大きさ r」を思い切り大きくしてみると、電場が弱まる様子が視覚的にわかるでしょ。これが「無限」モデルの限界を理解する第一歩だ。
次に、ガウス面の「適切さ」。法則自体はどんなでたらめな閉曲面でも成り立つけど、電場Eを求めたいなら「対称性」が命だ。球対称電荷で立方体のガウス面を取っても計算が爆発的に難しくなる。シミュレーターで「球対称」を選び、わざとガウス面を球からずらして動かしてみて。面の各場所で電場ベクトルの向きと強さがバラバラで、とても「E×(面積)」とは計算できないことが体感できるはず。実務でCAEツールを使う時も、メッシュ（計算格子）の取り方はこの「対称性」を考慮して決めることが多いんだ。
最後に、「外に出ない」の意味。「無限平面」で電場が面から垂直に出ているのを見て、「平面の『外側』には電場がないんだ」と誤解する人がいる。違うんだ。ここで言う「外側」とは、ガウス面として設定した円柱の側面のこと。無限平面からは両側に垂直に一様な電場が出ているけど、円柱の側面とは平行だから、側面を貫く電束はゼロなんだ。コンデンサの極板のように2枚の平面があれば、極板間の電場は $E = \sigma / \varepsilon_0$ とさらに強まる。この違い、シミュレーターでσの値を変えながらしっかり確認しておこう。
このシミュレーターで遊んでいる計算は、実は様々な先端技術の根底にあるんだ。まず真っ先に挙がるのは静電容量（キャパシタンス）の設計と評価だ。スマホの基板に載っているチップコンデンサや、電力系統で使われる大型コンデンサの基本原理は、まさに「無限平面」のモデルが拡張されたもの。極板間の電界分布を理解しておくことは、容量値の計算や絶縁破壊（放電）のリスク評価に直結する。
もう一つの大きな分野は高電圧・絶縁工学。送電線（これは「無限線電荷」のモデルに近い！）からの電界計算は、他の構造物や人体への影響、コロナ放電（電暈）の発生予測に必須だ。例えば、500kVの送電線から地表までの電界強度が規制値以下かどうかを検討する時、この基本原理がベースになる。また、ガス絶縁開閉装置（GIS）内部の導体まわりの三次元的な電界集中を、CAE（例えばANSYS MaxwellやCOMSOL Multiphysics）で解析する際も、ガウスの法則は数値計算の核心にある方程式（ポアソン方程式）の源流だ。
さらに半導体デバイス工学にも深く関わる。MOSFETのゲート酸化膜中の電界は、チャネルの形成や信頼性に決定的な影響を与える。ここでは極薄の酸化膜を「無限平面」とみなした一次元の電界計算が基礎として活躍する。このシミュレーターで「面電荷密度σ」を変えた時の電界の変化は、ゲート電圧を変えた時の酸化膜中の電界変化と概念的には同じなんだ。
このシミュレーターでガウスの法則の「感覚」が掴めたら、次は数式とベクトル解析にもう一歩踏み込んでみよう。まずは、シミュレーターの画面に表示される積分形の式 $\oint_S \mathbf{E}\cdot d\mathbf{A}= Q_{enc}/\varepsilon_0$ の意味を、各記号レベルで理解する。$d\mathbf{A}$ が「面素ベクトル」で、$\mathbf{E}\cdot d\mathbf{A}$ がベクトルの内積（= $E dA \cos\theta$）であることを確認しよう。これが「面を垂直に貫く成分だけを足し上げる」操作なんだ。
その上で、微分形への展開を学ぶと世界が広がる。ガウスの法則の積分形を、任意の微小な領域で適用すると、微分形 $\nabla \cdot \mathbf{E} = \rho / \varepsilon_0$ が導かれる。ここで $\rho$ は電荷密度[C/m³]だ。この式は「空間の一点一点で、電場の発散（湧き出し）が、その点の電荷密度で決まる」ことを意味する。これこそが、CAEソフトウェアが静電場を計算する時に解いているポアソン方程式の核心部分なんだ。
具体的な次のステップとしては、「重ね合わせの原理」と組み合わせた応用問題に挑戦するのがおすすめだ。例えば、2本の平行な無限線電荷（一方は正、他方は負）が作る電界分布を考えてみる。シミュレーターでは1種類の電荷分布しか選べないが、現実は複数の分布が組み合わさっている。それぞれが作る電場をベクトル和で重ね合わせることで、より複雑で実用的な電界パターンが計算できるようになる。これが、実際のCAE解析でメッシュ全体で行われている計算の本質的なイメージだ。