[논문 리뷰] A convergent linear finite element scheme for the Maxwell-Landau-Lifshitz-Gilbert equation

Ľubomír Baňas, Marcus Page외 1명|arXiv (Cornell University)|2013. 03. 16.

Advanced Numerical Methods in Computational Mathematics참고 문헌 16인용 수 8

한 줄 요약

이 논문은 최대 한 번의 선형 시스템 해법을 필요로 하는 두 개의 조건부 수렴성 없는 선형 유한요소 스킴을 제안한다. 이 스킴들은 자기장과 자화의 진화를 분리하여 처리하면서도, 효과적 필드에 대한 약한 가정 하에 약한 해로의 수렴성을 유지한다.

ABSTRACT

We consider a lowest-order finite element discretization of the nonlinear system of Maxwell's and Landau-Lifshitz-Gilbert equations (MLLG). Two algorithms are proposed to numerically solve this problem, both of which only require the solution of at most two linear systems per timestep. One of the algorithms is fully decoupled in the sense that each timestep consists of the sequential computation of the magnetization and afterwards the magnetic and electric field. Under some mild assumptions on the effective field, we show that both algorithms converge towards weak solutions of the MLLG system. Numerical experiments for a micromagnetic benchmark problem demonstrate the performance of the proposed algorithms.

연구 동기 및 목표

• 맥스웰-랜도-라이프시츠-기르버 시스템을 위한 수치적으로 효율적이고 조건부 수렴성 없는 시간 스텝 스킴을 개발한다.
• 랜도-라이프시츠-기르버 방정식을 위한 기존의 선형 적분기법을 자발적 자기장의 정적 근사가 아닌 전체 시간의존 맥스웰 시스템으로 확장한다.
• 각 시간단계에서 비선형 시스템 해법이 필요로 하지 않도록 보장함으로써, 은폐된 시간단계-메쉬 크기 결합 현상이 나타나는 은닉 중간점 또는 고정점 스킴에서 관찰되는 문제를 피한다.
• 자기장과 필드의 해를 분리한 전략(자기장 먼저, 그 다음 필드)을 수렴성 분석에 엄밀하게 포함시켜 안정성과 정확성을 보장한다.
• 비선형 반복을 최소화하면서도 조건부 수렴성을 유지함으로써, 더 큰 시간단계를 사용한 실용적인 시뮬레이션을 가능하게 한다.

제안 방법

• 두 가지 알고리즘을 제안한다: 알고리즘 2는 각 시간단계에서 하나의 선형 시스템만을 풀며, 완전히 분리된 구조를 가진다. 알고리즘 3은 각 시간단계에서 두 개의 선형 시스템을 풀며, 더 정확한 필드 갱신을 위해 반복 보정을 허용한다.
• 공간 근사에 최소 차수의 유한요소를 사용하며, 자화에는 노드 기반 라그랑주 요소를, 전자기장에는 네드레크 요소를 사용한다.
• 맥스웰-LLG 시스템에 크랭크-니콜슨 유사 시간 이산화를 적용하고, 자화 진화에 대해 중점 규칙을 기반으로 한 선형화 전략을 사용한다.
• 자기장의 단위 길이 제약 조건을 변분 형식 내에서 라그랑주 승수를 통해 도입하여 자화가 단위 구면 위에 유지되도록 보장한다.
• 유도된 선형 시스템은 자화에 대해 직접 해법기법을, 전자기장에 대해 다중 격자 전처리를 적용한 우자와 알고리즘을 사용하여 풀었다.
• 시간 스텝에서 정밀도와 안정성을 향상시키기 위해 예측-수정 전략을 도입하였으며, 알고리즘 3에서는 필드 갱신을 위해 고정점 반복을 사용한다.

실험 결과

연구 질문

• RQ1 전체 맥스웰-LLG 시스템에 대해 비선형 시스템을 각 시간단계에서 풀지 않고도 선형 유한요소 스킴을 구성할 수 있는가?
• RQ2 자기장을 먼저 풀고 전자기장을 이후에 풀는 완전히 분리된 알고리즘이 여전히 약한 해로의 조건부 수렴성을 보장하는가?
• RQ3 시간에 따라 변화하는 효과적 필드가 정적 또는 단순화된 모델이 아닌 전체 맥스웰 시스템에서 유래하는 경우, 수렴성 분석을 엄밀하게 확장할 수 있는가?
• RQ4 제안된 선형 스킴의 성능은 정밀도와 계산 비용 측면에서 기존의 비선형 스킴(예: 은닉 중간점 규칙)과 비교해 볼 때 어떻게 되는가?
• RQ5 제안된 스킴은 비선형성 및 다스케일 항을 포함한 일반적인 효과적 필드 기여를 처리하면서도 수렴성과 안정성을 유지할 수 있는가?

주요 결과

• 제안된 두 알고리즘은 효과적 필드에 대한 약한 가정 하에 맥스웰-LLG 시스템의 약한 해로의 조건부 수렴성을 모두 확보한다.
• 알고리즘 2는 각 시간단계에서 하나의 선형 시스템만을 풀며, 완전히 분리된 구조를 가지며, 자화를 먼저 풀고 나서 전자기장을 푼다.
• 알고리즘 3는 각 시간단계에서 두 개의 선형 시스템을 풀며, 반복 보정을 통해 더 정확한 필드 갱신을 가능하게 하면서도 수렴성을 유지한다.
• 수치 실험 결과, 중간점 규칙(기준 비선형 방법)과 매우 유사한 결과를 보였으며, 중간점 규칙이 요구하는 것보다 더 큰 시간단계(k=0.05)를 사용함에도 불구하고 정확도가 높았다.
• 해석 오차(예: 자화장의 비정확한 계산 등)에 대해 강건하며, 이러한 오차는 수렴성 분석에 포함될 수 있다.
• 기존의 고정점 반복이나 명시적 시간 스텝 기반 스킴들과 달리, 안정성과 수렴성에 손상 없이 최대 k=1까지의 큰 시간단계를 사용할 수 있다.