巻数・電流・長さ・比透磁率を変えて内部磁場・軸上磁場分布・インダクタンスをリアルタイム計算。
軸上の磁場分布。中央で均一、端部で約半値に低下する様子が確認できます。
各巻線（円電流）のビオ-サバール則から軸上磁場を積分します。1巻線の軸上磁場は $B = \frac{\mu_0 I r^2}{2(r^2+x^2)^{3/2}}$。有限長ソレノイド（N巻、長さL）の軸上x点での磁場は全巻線の積分：$B(x) = \frac{\mu_0 nI}{2}[\sin\alpha_2 - \sin\alpha_1]$（α は端部からの角度）。端部ではちょうど中心値の半分になります。
LCRメーターで周波数を変えたインピーダンス測定から求めるのが実用的です。理論値は $L = \mu_0 \mu_r N^2 A / l$（但し短いコイルは補正係数Kが必要：Nagaoka係数）。インダクタンスは巻数²に比例するので、巻数を2倍にすると4倍になります。
交流電流でソレノイドを励磁すると、鉄心の磁化はB-Hループ（ヒステリシスループ）を描きます。1サイクルあたりのループ面積に比例したエネルギーが熱として失われます（ヒステリシス損失）。電動機・変圧器の鉄損の主要成分で、珪素鋼板はこれを抑制するために設計されています。
通常の銅線コイルは電気抵抗 R > 0 のためジュール熱 I²R が発生します。超伝導体（液体ヘリウム冷却の NbTi, Nb₃Sn 等）は R=0 なので電気的損失なしに大電流を維持できます。一度励磁したら電源を切っても電流が流れ続ける「永久電流モード」が可能。MRI・核融合炉（ITER）・粒子加速器（LHC）の磁石に不可欠です。
ソレノイドに変動電流を流すと内部の磁束が変化し、ファラデーの法則（$\mathcal{E} = -d\Phi/dt$）から起電力が生じます。これが「インダクタンス L」の本質——電流変化を妨げる「慣性」として働きます。電源を切る瞬間に大電圧スパイクが発生するのはこの自己誘導起電力で、フライバックダイオードが必要になります。
ソレノイド磁場計算シミュレーターの物理モデルでは、有限長ソレノイドの中心軸上磁場をビオ・サバールの法則に基づき数値積分する。巻数 \(N\)、電流 \(I\)、長さ \(L\)、半径 \(a\)、比透磁率 \(\mu_r\) を入力とし、真空の透磁率 \(\mu_0\) を用いて、軸上位置 \(z\) における磁束密度 \(B(z)\) は次式で与えられる。 $ B(z) = \frac{\mu_0 \mu_r N I}{2L} \left( \frac{z + L/2}{\sqrt{a^2 + (z + L/2)^2}} - \frac{z - L/2}{\sqrt{a^2 + (z - L/2)^2}} \right) $ この式から、ソレノイド中心 (\(z=0\)) の内部磁場は \(B(0) = \frac{\mu_0 \mu_r N I}{\sqrt{L^2 + 4a^2}}\) と近似され、長さが半径より十分大きい理想ソレノイドでは \(B \approx \mu_0 \mu_r n I\)（\(n = N/L\)）に漸近する。インダクタンス \(L_{\text{ind}}\) は磁気エネルギー \(U = \frac{1}{2} L_{\text{ind}} I^2\) と体積積分から導出され、長岡係数を用いた補正式で計算する。これにより、巻数や比透磁率の変化が磁場分布や蓄積エネルギーに与える影響をリアルタイムに可視化できる。
産業での実際の使用例
自動車業界では、電磁弁（ソレノイドバルブ）の設計に活用されています。例えば、トヨタやデンソーが開発するハイブリッド車用の油圧制御ソレノイドでは、コイルの巻数やコア材の比透磁率を調整し、応答速度と保持力を最適化。また、半導体製造装置の精密ガス流量制御バルブでも、磁場分布の均一性を確認しながら製品設計が行われています。
研究・教育での活用
大学の電磁気学実験では、学生が巻数や電流を変えながら内部磁場とインダクタンスの変化をリアルタイムで観察し、理論式（ソレノイドの磁場公式）と比較する教材として利用。研究では、MRI用超伝導マグネットの予備設計や、核融合炉のプラズマ閉じ込め用コイルの磁場分布解析に、簡易シミュレーターとして導入されています。
CAE解析との連携や実務での位置付け
本シミュレーターは、本格的な3D電磁界CAE（例：ANSYS MaxwellやJMAG）の前段階として使用されます。設計初期にパラメトリックスタディを行い、大まかな仕様（巻数・寸法）を決定。その後、詳細な渦電流解析や熱連成解析へ移行することで、開発期間を短縮。実務では、試作回数の削減や設計変更の迅速な検討に貢献しています。
「巻数を増やせば内部磁場も比例して強くなる」と思いがちですが、実際にはソレノイドの長さが有限である場合、端部の磁場低下や漏れ磁束の影響を受けるため、無限長近似が成り立つ条件（長さ≫半径）から外れると計算値と実測値に乖離が生じます。特に比透磁率の高い磁性体コアを使用する際は、コアの飽和特性を考慮せずに線形計算を続けると、実際の磁場が予想より大幅に低くなる点に注意が必要です。また、インダクタンスは単に巻数と透磁率だけで決まると思われがちですが、実際にはコイルの形状や巻線の配置（分布容量や近接効果）にも依存するため、高周波領域では周波数依存性を無視できません。さらに、軸上磁場分布は中央部で最大になりますが、端部では急激に減少するため、均一磁場が必要な実験ではコイル長の中央付近のみを使用する設計が重要です。