[논문 리뷰] Hydrodynamic simulations on a moving Voronoi mesh

Volker Springel|arXiv (Cornell University)|2011. 09. 10.

Fluid Dynamics Simulations and InteractionsEngineering참고 문헌 2인용 수 21

한 줄 요약

이 논문은 이동 가능한 메쉬 생성 점들의 베이저노이 타일링을 사용하여 고정밀, 갈릴레이 불변성, 왜곡 없는 라그랑주 유체역학을 달성하는 새로운 이동 베이저노이 메쉬 방법을 소개한다. 이 방법은 비정규 메쉬 위에서 정확한 리만 해법을 사용하는 2차 정확도의 유한체적 Godunov 스킴을 조합하여, 고정 유럽형 또는 SPH 방법보다 더 낮은 해상도에서 더 높은 정확도와 수렴성을 보이며, 충격파, 접촉 불연속성, 대규모 유동이 있는 문제에서 특히 뛰어나다.

ABSTRACT

At the heart of any method for computational fluid dynamics lies the question of how the simulated fluid should be discretized. Traditionally, a fixed Eulerian mesh is often employed for this purpose, which in modern schemes may also be adaptively refined during a calculation. Particle-based methods on the other hand discretize the mass instead of the volume, yielding an approximately Lagrangian approach. It is also possible to achieve Lagrangian behavior in mesh-based methods if the mesh is allowed to move with the flow. However, such approaches have often been fraught with substantial problems related to the development of irregularity in the mesh topology. Here we describe a novel scheme that eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite volume approach, based on a second-order Godunov scheme with an exact Riemann solver. A particularly powerful feature of the approach is that the mesh-generating points can in principle be moved arbitrarily. If they are given the velocity of the local flow, a highly accurate Lagrangian formulation of continuum hydrodynamics is obtained that is free of mesh distortion problems, while it is at the same time fully Galilean-invariant, unlike ordinary Eulerian codes. We describe the formulation and implementation of our new Voronoi-based hydrodynamics, and we discuss a number of illustrative test problems that highlight its performance in practical applications.

연구 동기 및 목표

• 기존의 수치 유체역학 방법의 근본적 한계, 즉 SPH에서의 충격파 해상도 부족과 고정형 유크레안 코드에서의 갈릴레이 불변성 결여 문제를 해결하기 위해.
• 메쉬 뒤틀림이나 위상적 잡음 없이 라그랑주 행동을 달성하는 메쉬 기반 방법을 개발하기 위해.
• SPH의 기하학적 유연성과 메쉬 기반 방법의 정확성 및 보존 성질을 결합하기 위해.
• 강한 밀도 대비와 대규모 유동을 포함한 천체 유동의 적응형 해상도와 정확한 시뮬레이션을 가능하게 하기 위해.

제안 방법

• 메쉬는 이동 가능한 점들의 집합에 대한 베이저노이 타일링으로 구성되며, 유동에 따라 지속적으로 변형되어 메쉬 품질을 유지한다.
• 유체역학은 비선형 보존 법칙에 대한 정확한 리만 해법을 사용하는 2차 정확도의 유한체적 Godunov 스킴으로 해결된다.
• 메쉬 생성 점들이 국소 유체 속도에 따라 이동하여 거의 라그랑주 형식을 달성하고, 이동 오차를 최소화한다.
• 베이저노이 타일링이 자연스럽게 메쉬 뒤틀림이나 뱃줄 모양 셀을 방지하므로 위상적 안정성이 보장된다.
• 메쉬 운동이 유체의 운동을 따라가므로 메쉬 운동이 유체의 운동에 따라 유도되어 모든 속도에 의존하는 수치적 잡음 요소가 제거되어, 방법은 완전히 갈릴레이 불변성이 보장된다.
• 비정규 메쉬 프레임워크 덕분에 적응형 해상도와 곡면 경계 조건을 자연스럽게 지원한다.

실험 결과

연구 질문

• RQ1 베이저노이 타일링 기반의 이동 메쉬 방법이 다차원 유체역학에서 갈릴레이 불변성을 달성하고 메쉬 뒤틀림을 피할 수 있는가?
• RQ2 표준 유체역학 테스트 문제에서 이동 베이저노이 메쉬의 정확도가 고정형 유크레안 및 SPH 코드보다 어떻게 비교되는가?
• RQ3 대규모 유동 또는 접촉 불연속성이 있는 시뮬레이션에서 이동 메쉬가 이동 오차를 얼마나 줄이고 수렴성을 향상시키는가?
• RQ4 강한 밀도 대비와 복잡한 유동 형태(예: 켈빈-헬름홀츠 불안정성 또는 클라우드 파괴)를 포함한 문제에서 이 방법이 높은 정확도를 유지할 수 있는가?

주요 결과

• 고정형 메쉬 유크레안 코드보다 낮은 해상도에서 클라우드 파괴 테스트에서 수렴 결과를 달성하여, 더 뛰어난 수치적 효율성을 나타낸다.
• 이동 메쉬 방법은 이동 오차 감소와 갈릴레이 불변성 덕분에 고정형 메쉬 코드보다 항상 더 높은 잔류 클라우드 질량 분율을 예측한다.
• SPH보다 충격파 및 접촉 불연속성의 해상도가 향상되었으며, 클라우드 파괴 시뮬레이션에서 질량 손실을 과소평가하는 SPH의 문제를 해결한다.
• 코드는 인위적 점성 없이 공간 및 시간에서 2차 정확도를 보이며, SPH에서 흔히 발생하는 비현실적인 표면 장력 효과를 피한다.
• 3D에서 코드는 해상도 요소당 SPH보다 약 2배 느리고, 고정형 메쉬 유크레안 코드보다 3~4배 느리지만, 이 비용은 낮은 해상도에서의 높은 정확도로 상쇄된다.
• 기하학적 탄력성과 위상적 안정성 덕분에, 곡면 경계와 강한 대규모 유동을 포함한 복잡한 흐름에서도 방법은 강력하다.