[논문 리뷰] A Four Function Variational Principle for Barotropic Magnetohydrodynamics

Asher Yahalom|arXiv (Cornell University)|2008. 11. 14.

Solar and Space Plasma Dynamics참고 문헌 24인용 수 12

한 줄 요약

이 논문은 비정상적인 등온성 자기유체역학에 대한 네 기능 유도형 변분 원리를 제안하며, 속도, 자기장, 밀도의 7개 필드를 χ, η, ν, ρ의 네 스칼라 함수로 감소시킨다. 이 형식은 속도 및 자기장의 잠재함수 표현을 기반으로 한 라그랑지안 밀도를 사용하여, 변분 도함수를 통해 모든 지배 방정식을 도출하며, 계산 복잡도가 감소함에 따라 수치 시뮬레이션과 안정성 분석이 가능해진다.

ABSTRACT

Variational principles for magnetohydrodynamics were introduced by previous authors both in Lagrangian and Eulerian form. In a previous work Yahalom & Lynden-Bell introduced a simpler Eulerian variational principles from which all the relevant equations of magnetohydrodynamics can be derived. The variational principle was given in terms of six independent functions for non-stationary flows and three independent functions for stationary flows. This is less then the seven variables which appear in the standard equations of magnetohydrodynamics which are the magnetic field $\vec B$ the velocity field $\vec v$ and the density $ρ$. In this work I will improve on the previous results showing that non-stationary magnetohydrodynamics should be described by four functions .

연구 동기 및 목표

• 비정상적인 등온성 자기유체역학에서 독립 변수의 수를 7개에서 4개로 감소시키기.
• 추가 제약 조건 없이 모든 MHD 방정식을 생성하는 변분 원리 개발하기.
• 야하롬과 린든벨의 이전 여섯 기능 형식에 비해 더 컴act하고 계산 효율성이 높은 형식으로 개선하기.
• 대칭에 적합한 라그랑지안 프레임워크를 통해 직접 수치 시뮬레이션과 안정성 분석을 가능하게 하기.

제안 방법

• 자기장 B = ∇χ × ∇η로 표현하여 ∇·B = 0이 자동으로 만족되도록 하기.
• 속도장 v = ∇ν + α∇χ + β∇η로 표현하여 클레브시 표현을 일반화하기.
• χ, η, ν, ρ에 의존하는 라그랑지안 밀도를 구성하며, 시간 도함수와 공간 기울기를 포함한 명시적 형태 제공하기.
• 행동의 변분 도함수를 0으로 설정하여 운동 방정식 유도하기, 이를 통해 보존 법칙과 동역학 방정식 도출하기.
• 자기행렬식이 1인 무한 대칭군을 활용하여 노이터 정리에 의해 두 개의 국소 보존 법칙 도출하기.
• 이중 차수 방정식을 두 개의 일차 및 두 개의 이차 PDE 혼합계로 변환하여 수치 해석 가능하게 하기.

실험 결과

연구 질문

• RQ1 비정상적인 설정에서 등온성 자기유체역학을 여섯 개 이하의 독립 함수로 기술할 수 있는가?
• RQ2 보조 제약 조건 없이 모든 MHD 방정식을 생성할 수 있는 변분 원리는 어떻게 구성할 수 있는가?
• RQ3 비정상적인 등온성 자기유체역학을 완전히 기술하기 위해 필요한 최소 스칼라 함수의 집합은 무엇인가?
• RQ4 라그랑지안의 대칭성은 교차 운동량과 같은 보존 법칙을 어떻게 도출하는가?
• RQ5 이 감소된 형식은 수치 시뮬레이션의 효율성과 안정성 분석을 향상시킬 수 있는가?

주요 결과

• 논문은 비정상적인 등온성 자기유체역학을 기술하는 데 χ, η, ν, ρ의 네 스칼라 함수로 충분하다고 규명한다.
• 변분 원리는 첫 번째 변분 도함수를 통해 연속 방정식, 운동량 방정식, α와 β의 진화 방정식을 도출한다.
• 라그랑지안 밀도는 자동행렬식이 1인 무한 대칭군을 갖으며, 노이터 정리에 의해 두 개의 국소 보존 법칙을 도출한다.
• 시스템은 χ와 η에 대해 두 개의 일차 및 두 개의 이차 PDE로 감소하여 수치적 통합에 더 효율적인 형태로 전환된다.
• 기존의 FLUIDEX 소프트웨어를 통해 입증된 바와 같이, 원래의 장 방정식을 풀지 않고도 직접 수치 시뮬레이션을 수행할 수 있다.
• 추가 보존 법칙을 액션에 통합함으로써 향후 안정성 분석의 기초를 제공한다.