비압축성 유체와 얇은 변형체를 위한 위상·기하 정확 결합

초록

본 논문은 라그랑지안 입자 기반 Voronoi 다이어그램을 클리핑하고 ‘스티칭’ 과정을 통해 얇은 구조물 주변의 유체 도메인 연결성을 정확히 보존하는 새로운 격자 생성 기법을 제안한다. 이 기법은 유체‑구조 경계에서 압력 라그랑주 승수를 이용해 속도 경계조건을 강제하고, 구조에 압력력을 직접 전달함으로써 양방향 강한 결합을 구현한다. 결과적으로 얇은 막이나 협소한 통로를 포함한 복잡한 씬에서도 인위적인 누수나 차단 없이 정확한 유동을 재현한다.

상세 분석

이 연구는 기존 유체‑구조 결합 방법이 얇은 코디멘션(공변 차원) 구조물과의 상호작용에서 겪는 ‘누수(leakage)’와 ‘인위적 씰링(sealing)’ 문제를 근본적으로 해결한다는 점에서 혁신적이다. 핵심 아이디어는 라그랑지안 유체 입자들을 Voronoi 사이트로 사용해 자연스럽게 비정형 셀을 생성하고, 고체 경계면을 클리핑 평면으로 삽입해 유체 셀을 고체와 정확히 맞닿게 만든다. 클리핑 과정에서 발생하는 ‘고아 셀(orphaned cells)’은 인접 입자에게 재할당하는 스티칭 알고리즘을 적용해 유체 도메인의 경로 연결성을 보존한다. 이때 셀 재배치는 위상 보존을 전제로 하며, 입자 샘플링이 부족한 협소 통로에서도 유효한 흐름 경로를 복원한다는 점이 특히 주목할 만하다.

생성된 비정형 Voronoi 메시는 유체의 압력 투영(pressure projection) 단계에서 직접 사용된다. 압력 포아송 방정식의 이산화는 유체‑유체와 유체‑고체 인터페이스 모두를 포함하는 면에 대해 수행되며, 고체 면에 적용되는 압력 라그랑주 승수는 속도 경계조건을 강제하고 동시에 고체에 물리적 압력력을 전달한다. 따라서 ‘속도-압력’ 쌍을 동일 메시에 두어 양방향 결합을 구현함으로써, 기존의 Immersed Boundary나 Cut‑Cell 방식이 겪는 인터페이스 정밀도 저하와 자유도 폭발 문제를 회피한다.

또한 저자들은 동적 위상 변화(예: 심장 판막의 개폐) 상황에서 압력 게이지(pressure gauge)의 일관성을 유지하기 위한 전략을 제시한다. 유체 도메인이 여러 연결 컴포넌트로 분리·재결합될 때 각 컴포넌트의 압력 영(pressure nullspace)를 조정해 전역적인 압력 기준을 보존한다. 이는 수치적 안정성과 물리적 정확성을 동시에 확보하는 중요한 설계 선택이다.

실험에서는 얇은 막을 통과하는 흐름, 다층 얇은 시트가 형성하는 협소 통로, 변형 가능한 심장 판막 등 다양한 시나리오를 테스트했다. 특히 얇은 판막이 정상·병변 상태에서 보이는 역류 현상을 정확히 재현했으며, 고체가 매우 얇아 격자 셀보다 작을 경우에도 유체가 차단되지 않고 정상 흐름을 유지한다는 점을 시각적으로 입증했다. 전체적으로 이 방법은 입자 수에 비례하는 고정된 자유도를 유지하면서도, 복잡한 기하와 위상 변화를 정확히 포착한다는 점에서 기존 방법에 비해 확장성 및 효율성이 크게 향상되었다.