[논문 리뷰] A Two-dimensional eddy current model using thin inductors

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

Magnetic Properties and Applications참고 문헌 6인용 수 9

한 줄 요약

이 논문은 유도체 두께 $\epsilon \to 0$일 때의 점근적 분석을 이용하여 얇은 유도체를 위한 2차원 에디스 전류 모델을 개발한다. 전류를 지닌 얇은 도체선을 나타내는 딜라크 측도를 포함하는 특이 타원형 문제를 유도하며, $L^2$에서 수렴 속도 $O(\epsilon^\alpha)$($0 < \alpha < 1/2$)와 $W^{1,p}$에서 약한 수렴을 증명한다. 이는 전자기계 시스템 내 얇은 도체의 정확하고 잘 조절된 수치 시뮬레이션을 가능하게 한다.

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

연구 동기 및 목표

• 두 차원 기하학에서 얇은 유도체 내 에디스 전류에 대한 단순화되고 수학적으로 엄밀한 모델을 유도하는 것.
• 유도체 두께 $\epsilon \to 0$일 때의 점근적 분석을 이용하여 얇은 도체 메쉬 생성의 수치적 과제를 해결하는 것.
• 유도체 두께가 0으로 수렴할 때의 해가 전류원을 나타내는 딜라크 측도를 포함한 한계 문제로 수렴하는 것을 확립하는 것.
• 해의 수렴 속도를 $L^2$에서와 $W^{1,p}$ 공간에서 약한 수렴으로 증명하여 수치 근사의 안정성을 확보하는 것.

제안 방법

• 전류가 $x_3$-방향으로 흐르는 원통형 기하에서 3차원 에디스 전류 문제를 수립한다.
• 벡터 포텐셜과 컬 연산자를 이용하여 3차원 문제를 스칼라 포텐셜 형태로 2차원으로 축소한다.
• 지름이 $O(\epsilon)$인 얇은 유도체를 가정하고, $\epsilon \to 0$일 때의 한계 문제를 유도한다.
• 점근적 분석을 통해 한계 문제의 해가 유도체 위치에 딜라크 측도를 포함하는 특이 타원형 편미분방정식이 됨을 보인다.
• 에너지 추정과 소볼레프 포함을 통해 $L^2$에서 수렴 속도 $O(\epsilon^\alpha)$($0 < \alpha < 1/2$)를 확립한다.
• 일관된 유계성과 컴actness를 이용하여 $W^{1,p}$에서 $1 < p < 2$에 대해 약한 수렴을 증명한다.

실험 결과

연구 질문

• RQ1 두께가 0으로 수렴할 때 2차원 기하에서 에디스 전류 문제의 해는 어떻게 행동하는가?
• RQ2 고정된 전류를 지닌 얇은 유도체가 두께가 0이 되었을 때 스칼라 포텐셜의 한계 행동은 어떠한가?
• RQ3 해의 수렴 속도를 $\epsilon \to 0$일 때 $L^2$에서 정량화할 수 있는가?
• RQ4 한계 문제는 전류 보존과 자기장 연속성의 물리적 구조를 유지하는가?
• RQ5 소볼레프 공간 $W^{1,p}$에서 $1 < p < 2$에 대해 해의 정(regularity)과 수렴 행동은 어떠한가?

주요 결과

• 해는 $L^2$에서 수렴 속도 $O(\epsilon^\alpha)$($0 < \alpha < 1/2$)로 한계 문제로 수렴한다.
• 모든 유계 구역 $B$와 $1 < p < 2$에 대해 $W^{1,p}(B)$에서 약한 수렴이 성립하며, $L^q(B)$에서 $1 \leq q < 2p/(2-p)$에 대해 수렴한다.
• 한계 문제는 유도체 위치에 딜라크 측도를 포함하는 특이 타원형 편미분방정식이다: $-\Delta u + i\omega\mu\sigma u = \mu\sigma C_k$ in $\Omega_k$, $\Delta u = 0$ in $\Omega'$, $\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$로도 결과가 확장되며, 임의의 수의 도체로도 일반화된다.