[논문 리뷰] A General Approach for Producing Hamiltonian Numerical Schemes for Fluid Equations

Colin J. Cotter|ArXiv.org|2005. 01. 26.

Numerical methods for differential equations참고 문헌 34인용 수 6

한 줄 요약

이 논문은 라그랑주 입자 운동과 오일러 메esh 기반 해밀토니안 평가를 조합하여 유체 방정식을 위한 준연속 해밀토니안 스킴을 구축하기 위한 일반적인 입자-메쉬 방법을 제시한다. 이 방법은 심플렉틱 시간 적분을 위한 캐논리컬 구조를 보장하며, 장시간에 걸쳐 O(Δt) 이내로 해밀토니안을 유지한다. 2D EP-Diff 방정식에 대한 실험을 통해 뛰어난 보존 성질을 입증하였다.

ABSTRACT

Given a fluid equation with reduced Lagrangian $l$ which is a functional of velocity $\MM{u}$ and advected density $D$ given in Eulerian coordinates, we give a general method for semidiscretising the equations to give a canonical Hamiltonian system; this system may then be integrated in time using a symplectic integrator. The method is Lagrangian, with the variables being a set of Lagrangian particle positions and their associated momenta. The canonical equations obtained yield a discrete form of Euler-Poincaré equations for $l$ when projected onto the grid, with a new form of discrete calculus to represent the gradient and divergence operators. Practical symplectic time integrators are suggested for a large family of equations which include the shallow-water equations, the EP-Diff equations and the 3D compressible Euler equations, and we illustrate the technique by showing results from a numerical experiment for the EP-Diff equations.

연구 동기 및 목표

• 유체 방정식의 캐논리컬 해밀토니안 준연속 이산화를 체계적으로 구축할 수 있는 접근법을 개발하는 것.
• 오일러 좌표계에서 비캐논리컬 리-포아송 구조를 이산화하는 데 어려움을 해결하기 위해 라그랑주 입자를 사용하여 캐논리컬 형태를 복원하는 것.
• 입자 기반 운동과 메쉬 기반 해밀토니안 평가를 조합하여 유체 시뮬레이션에서 해밀토니안의 장기적 보존을 가능하게 하는 것.
• 연속형 수류, EP-Diff, 압축성 에일러 방정식을 포함한 다양한 유체 모델에 적용 가능한 프레임워크를 제공하는 것.
• 2D EP-Diff 방정식에 대한 수치 실험을 통해 계산의 타당성과 안정성을 입증하는 것.

제안 방법

• 심플렉틱 시간 적분을 위한 캐논리컬 구조를 보장하기 위해 라그랑주 입자의 위치와 운동량을 캐논리컬 변수로 사용하는 것.
• 입자 기반 운동량과 밀도를 오일러 메쉬로 보간하여 해밀토니안을 격점의 합으로 평가하는 것.
• 이산 미분법 프레임워크를 사용해 이산 오일러-포아송 방정식의 기울기 및 발산 연산자를 표현하는 것.
• 감소된 라그랑주언의 레지플레스 변환을 통해 이산 해밀토니안 시스템을 유도하고, 입자-메쉬 보간을 통해 격자에 투영하는 것.
• 심플렉틱 시간 적분기—예를 들어 심플렉틱 오일러 방법—을 준연속 시스템에 적용하여 장시간에 걸쳐 해밀토니안을 유지하는 것.
• 이산 방정식이 이산 오일러-포아송 방정식의 형태를 유지하도록 보장하여 변분적 구조를 유지하는 것.

실험 결과

연구 질문

• RQ1 오일러 좌표계에서 감소된 라그랑주언이 있는 유체 방정식에 대해 캐논리컬 해밀토니안 준연속 이산화를 생성할 수 있는 일반적 방법을 개발할 수 있는가?
• RQ2 오일러 좌표계에서의 비캐논리컬 구조가 심플렉틱 시간 적분에 적합한 캐논리컬 형태로 어떻게 변환될 수 있는가?
• RQ3 입자-메쉬 보간은 해밀토니안 구조를 유지하면서도 해밀토니안의 격자 기반 평가를 가능하게 하는 데 어떤 역할을 하는가?
• RQ4 심플렉틱 시간 적분기는 준연속 시스템에서 장시간에 걸쳐 해밀토니안을 얼마나 잘 유지하는가?
• RQ5 이 프레임워크는 얕은 수면-α 및 3D 압축성 에일러 방정식을 포함한 다양한 유체 모델에 적용될 수 있는가?

주요 결과

• 장시간에 걸쳐 이산 해밀토니안은 초기 값으로부터 O(Δt) 이내에 유지되어 심플렉틱 시간 적분의 효과성을 확인한다.
• 2D EP-Diff 방정식에 대한 수치 실험은 폭주나 비물리적 진동이 없는 잘 조율된 해를 보여준다.
• 이 방법은 H¹ 공간에서의 약한 해를 허용하여 δ-함수와 같은 특이한 운동량 분포를 처리할 수 있다.
• 입자-메쉬 형식은 속도 장에 대해 유한요소를 사용할 수 있도록 지원하여 견고한 공간 이산화를 가능하게 한다.
• 이 방법은 라그랑주 변수가 얽히지 않은 채로 오일러 형태의 복잡한 미분 연산자를 포함할 수 있도록 한다.
• 이 프레임워크는 스칼라, 벡터, 텐서를 포함한 다양한 운반 물질을 가진 준직접곱 시스템으로 확장 가능하다.