[논문 리뷰] Calculation of interface curvature with the level-set method

Karl Yngve Lervåg|arXiv (Cornell University)|2011. 01. 01.

Fluid Dynamics and Heat Transfer참고 문헌 7인용 수 4

한 줄 요약

이 논문은 계면 곡률과 법선 벡터를 계산하기 위한 강건하고 구현이 간편한 이산화 방법을 제안한다. 특히 다공성 변화가 발생하는 영역, 예를 들어 물방울의 융합과 같은 경우에 유용하다. 기하학적 인식 곡률 추정과 뾰족한 부분에서의 방향 차이를 적응적으로 적용함으로써, 기존의 유한차분 방법에서 발생하는 비물리적인 곡률 피크와 잘못된 압력장 문제를 해결하여, 위상 전이가 발생하는 이중상 유동의 정확한 시뮬레이션을 가능하게 한다.

ABSTRACT

The level-set method is a popular method for interface capturing. One of the advantages of the level-set method is that the curvature and the normal vector of the interface can be readily calculated from the level-set function. However, in cases where the level-set method is used to capture topological changes, the standard discretization techniques for the curvature and the normal vector do not work properly. This is because they are affected by the discontinuities of the signed-distance function half-way between two interfaces. This article addresses the calculation of normal vectors and curvatures with the level-set method for such cases. It presents a discretization scheme that is relatively easy to implement in to an existing code. The improved discretization scheme is compared with a standard discretization scheme, first for a case with no flow, then for a case where two drops collide in a shear flow. The results show that the improved discretization yields more robust calculations in areas where topological changes are imminent.

연구 동기 및 목표

• 위상 변화가 발생할 경우 기존의 유한차분 방법을 사용할 때 곡률과 법선 벡터 계산의 부정확성을 해결하기 위해.
• 근접한 계면 사이의 뾰족한 부분에서 부호거리 함수의 불연속성으로 인해 발생하는 곡률 피크 문제를 해결하기 위해.
• 이미 존재하는 레벨셋 코드베이스에 큰 아키텍처 변경 없이도 쉽게 통합할 수 있는 강건하고 정확한 곡률 이산화 방법을 개발하기 위해.
• 정적 및 동적 이중상 유동 사례, 특히 융합 근처에서 표준 방법과의 비교를 통해 향상된 방법의 유효성을 검증하기 위해.
• 향상된 방법이 곡률 오류로 인한 잘못된 압력장 문제를 방지함으로써 물리적인 융합을 가능하게 함을 보여주기 위해.

제안 방법

• 맥클린 및 루엔그럽(2011)의 방향 차이 방법을 사용하여 법선 벡터를 계산함으로써, 기존 기울기 방법보다 뾰족한 부분을 더 잘 다룰 수 있도록 함.
• 맥클린 등(2012)의 기하학적 인식 곡률 이산화 방법을 도입하여, 레벨셋 함수의 이阶 도함수에 의존하지 않고 국소 계면 기하학을 적합시켜 곡률을 추정함.
• 고정된 오일러 메esh 위에서 유한차분 프레임워크에 향상된 곡률 이산화 방법을 적용함으로써, 기존 레벨셋 구현에 최소한의 변경만을 요구함.
• 표면장력력을 연속 표면력(CSF) 방법으로 모델링한 2차원 비압축성 이중상 유동 해석기에서 방법을 구현함.
• 부호거리 레벨셋 함수를 사용하고 재초기화를 적용하여, 수치 오차를 줄이고 거리 함수 성질을 유지함.
• 고정 격자 유한차분 방법을 사용하여 라우아르-나비에-스토크스 방정식과 레벨셋 이송을 해결하며, 매 시간단계에서 곡률과 법선 벡터를 계산함.

실험 결과

연구 질문

• RQ1 향상된 곡률 이산화 방법이 근접한 두 계면 사이의 뾰족한 부분에서 비물리적인 곡률 피크를 제거하는가?
• RQ2 근접한 두 물방울 사이의 박막 영역에서 표준 방법이 유도하는 곡률 피크로 인한 압력장은 어떻게 영향을 미치는가?
• RQ3 표준 방법이 실패하는 전단류 흐름 시뮬레이션에서 향상된 방법이 물리적인 물방울 융합을 가능하게 하는가?
• RQ4 곡률 정확도와 안정성 측면에서 향상된 방법이 기존의 유한차분 방법과 비교해 어떻게 성능을 발휘하는가?
• RQ5 향상된 곡률 이산화 방법은 큰 수정 없이도 기존의 레벨셋 코드에 강건하고 실용적으로 통합 가능한가?

주요 결과

• 기존의 유한차분 곡률 이산화 방법은 근접한 두 물방울 사이의 뾰족한 부분에서 심각한 곡률 피크를 생성하여 비물리적인 압력 왜곡을 유도한다.
• 향상된 곡률 방법은 뾰족한 부분에서 곡률 피크를 성공적으로 제거하여 매끄럽고 물리적으로 일관된 곡률 분포를 제공한다.
• 표준 방법에서 곡률 피크로 인해 발생하는 잘못된 압력장은 융합을 억제하여 전단류 흐름 사례에서 물방울 융합을 방지한다.
• 향상된 이산화 방법을 사용할 경우, 박막 영역의 압력장은 중심에서 높고 가장자리에서 낮은 물리적으로 정확한 분포를 보이며, 융합을 촉진하는 외향 유동을 유도한다.
• 향상된 방법 시뮬레이션에서는 융합이 성공적으로 발생하지만, 표준 방법에서는 곡률 오류로 인한 비물리적 힘으로 인해 융합이 차단된다.
• 향상된 곡률 이산화 방법은 최소한의 코드 변경만으로도 기존의 레벨셋 구현에 실용적이고 강건한 개선을 가능하게 한다.