[論文レビュー] A Two-dimensional eddy current model using thin inductors

Youcef Amirat, Rachid Touzani|ArXiv.org|Jul 11, 2006

Magnetic Properties and Applications参考文献 6被引用数 9

ひとこと要約

本稿では、インダクタの厚さ $\epsilon \to 0$ の極限において、漸近解析を用いて薄いインダクタの2次元渦電流モデルを構築する。電流を運ぶ薄い導線を表すデルタ関数を含む特異な楕円型問題が導かれる。$L^2$ での収束率 $O(\epsilon^\alpha)$（$0 < \alpha < 1/2$）および $W^{1,p}$ での弱収束（$1 < p < 2$）が証明され、電磁気系における薄い導体の高精度で良好に条件付けられた数値シミュレーションが可能になる。

ABSTRACT

We derive a mathematical model for eddy currents in two dimensional geometries where the conductors are thin domains. We assume that the current flows in the $x\_3$-direction and the inductors are domains with small diameters of order $O(ε)$. The model is derived by taking the limit $ε o 0$. A convergence rate of $O(ε^α)$ with $0

研究の動機と目的

• 2次元幾何における薄いインダクタの渦電流に対して、簡略化され、数学的に厳密なモデルを導出すること。
• インダクタの厚さ $\epsilon \to 0$ の極限において、漸近解析を用いて、薄い導体のメッシュ生成に起因する数値的課題に対処すること。
• 薄いインダクタの問題の解が、電流源を表すデルタ関数を含む極限問題に収束することを確立すること。
• $L^2$ での収束率と $W^{1,p}$ 空間（$1 < p < 2$）での弱収束を証明し、数値近似の堅牢性を保証すること。

提案手法

• 電流が $x_3$ 方向に流れることを仮定した円柱座標系における3次元渦電流問題を定式化する。
• ベクトルポテンシャルと回転作用素を用いて、3次元問題を2次元スカラーポテンシャル形式に簡略化する。
• 直径 $O(\epsilon)$ の薄いインダクタを仮定し、$\epsilon \to 0$ の極限における問題を導出する。
• 漸近解析を用いて、インダクタ位置にデルタ関数を含む特異な楕円型偏微分方程式に極限問題が帰着することを示す。
• エネルギー推定とソボレフ埋没定理を用いて、$L^2$ での収束率 $O(\epsilon^\alpha)$（$0 < \alpha < 1/2$）を確立する。
• $W^{1,p}$ 空間におけるコンパクトネスと一様有界性を用いて、$1 < p < 2$ での弱収束を証明する。

実験結果

リサーチクエスチョン

• RQ1 2次元幾何において、薄いインダクタの厚さがゼロに近づくとき、渦電流問題の解はどのように振る舞うか？
• RQ2 固定された電流を有する薄いインダクタが厚さゼロに縮小されるとき、スカラーポテンシャルの極限的挙動は何か？
• RQ3 $\epsilon \to 0$ の極限において、$L^2$ での解の収束率を定量的に評価できるか？
• RQ4 極限問題は、電流保存則と磁場の連続性という物理的構造を保持しているか？
• RQ5 $1 < p < 2$ のソボレフ空間 $W^{1,p}$ における解の正則性と収束挙動は何か？

主な発見

• 解は $L^2$ で収束率 $O(\epsilon^\alpha)$（$0 < \alpha < 1/2$）で極限問題に収束する。
• $B$ を任意の有界球体として、すべての $B$ および $1 < p < 2$ に対して $W^{1,p}(B)$ での弱収束が成り立ち、$L^q(B)$（$1 \leq q < 2p/(2-p)$）でも収束する。
• 極限問題は、インダクタ位置にデルタ関数を含む特異な楕円型偏微分方程式である：$\Omega_k$ で $-\Delta u + i\omega\mu\sigma u = \mu\sigma C_k$、$\Omega'$ で $\Delta u = 0$、境界 $\Gamma$ でジャンプ条件を満たす。
• テスト関数 $\varphi \in W^{1,p'}(B)$ に対して、極限解は $\int_B \nabla u \cdot \nabla \bar\varphi \, dx + i\beta \int_{\Omega_0} u \bar\varphi \, dx = \mu I (\bar\varphi(z_1) - \bar\varphi(z_2))$ を満たす。
• 厚い導体が存在しない場合（$\Omega_0 = \emptyset$）、解は明示的に $u(x) = \frac{\mu I}{2\pi} \log \frac{|x - z_2|}{|x - z_1|}$ と与えられる。
• もし $\mu \in W^{1,\infty}$ であれば、$\mu$ や $\sigma$ が定数でない場合にも結果は拡張可能であり、任意の数の導体に対しても成立する。