BDF2 기반 LLG 적분기: 무조건 약수렴과 2차 시간 정확도

초록

본 논문은 마이크로자기학에서 사용되는 Landau‑Lifshitz‑Gilbert 방정식(LLG)을 1차 유한요소와 2단계 BDF2 시간 적분법으로 전산화한다. 제시된 알고리즘은 각 시간 단계마다 선형 시스템 하나만 풀면 되며, CFL 조건 없이도 에너지 안정성을 보장한다. 재구성된 공간‑시간 자기화는 $H^1$ 약수렴을 통해 LLG의 약해 해로 수렴함을 증명하고, 수치 실험을 통해 공간 1차, 시간 2차 정확도를 확인한다.

상세 분석

이 연구는 LLG 방정식의 비선형성, 단위벡터 제약(|m|=1) 및 마그네틱 포텐셜의 복합 구조를 고려하면서, 계산 효율성을 극대화하는 새로운 BDF2‑type 적분기를 제안한다. 공간은 $S^1(\mathcal T_h)$, 즉 각 요소가 선형인 연속 유한요소를 사용하고, 노드 수준에서 $|m_h(z)|=1$을 강제함으로써 전통적인 Lagrange multiplier 없이도 제약을 만족한다. 시간 이산화는 두 단계 뒤로 차분(BDF2)을 적용하되, 첫 번째 단계는 기존의 1차 접선 평면(tangent‑plane) 방법으로 초기값 $m^1_h$을 얻는다. 핵심은 각 단계에서 이산 접선 공간 $T_h(\varphi_h)={\psi_h\in S^1(\mathcal T_h):\varphi_h(z)\cdot\psi_h(z)=0;\forall z\in\mathcal N_h}$에 제한된 선형 시스템을 푸는 것이다. 이 접근법은 비선형 연산자를 선형화하면서도 물리적 제약을 보존한다는 장점을 가진다.

수렴 분석은 에너지 방법을 기반으로 한다. 먼저 알고리즘이 Lax‑Milgram 정리를 만족해 유일한 해를 갖는 것을 보이고, 연속적인 시간 보간 $m_{h\tau}$와 전진/후진 보간 $m_h^\pm$를 정의한다. 이들 보간 함수는 $H^1(\Omega_T)$에서 균등 유계성을 가지며, 따라서 Banach‑Alaoglu 정리와 Rellich‑Kondrachov 정리를 이용해 약수렴(또는 강수렴) 부분수열을 추출한다. 약해 해 정의에 필요한 네 가지 조건(i)–(iv)를 차례로 검증한다. 특히 (iv) 에너지 부등식은 $\eta_0+\eta_n\to0$이라는 추가 조건을 통해 확보되는데, 이는 $\tau=o(h^2)$인 CFL‑type 조건이 첫·마지막 단계에만 필요함을 의미한다.

이 논문은 기존의 약수렴 결과를 제공한