복잡한 도메인 중간 레이놀즈 수 흐름을 위한 격자볼츠만 경계조건 선택

초록

본 연구는 혈관 내 중간 레이놀즈 수와 비정상 흐름을 모델링하기 위해 격자볼츠만법(LBM)에서 사용되는 경계조건을 비교한다. 단순 반사(bounce‑back)와 Guo‑Zheng‑Shi(GZS), Bouzidi‑Firdaouss‑Lallemand(BFL), Junk‑Yang(JY) 세 가지 고급 방법을 구현하고, 임의 방향의 원통 파이프와 곡선 파이프에서 Poiseuille, Womersley, Dean 흐름을 시뮬레이션하였다. 결과는 BFL과 GZS가 2차 공간 수렴을 보이며, BFL이 비정상 흐름에서 더 높은 정확도와 낮은 연산 비용을 제공함을 보여준다. JY는 고레일놀즈 수에서 불안정해 사용을 권장하지 않는다. 충돌 연산자와 속도 집합의 선택은 정확도에 큰 영향을 주지 않았다.

상세 분석

이 논문은 혈관과 같은 복잡한 생체 구조에서 격자볼츠만법(LBM)을 적용할 때 가장 큰 난관 중 하나인 경계조건 선택 문제를 체계적으로 조사한다. 기존에 가장 널리 쓰이는 단순 반사(bounce‑back, BB) 방식은 구현이 간단하고 계산 비용이 낮지만, 격자와 경계면이 비정렬될 경우 1차 정확도로 수렴한다는 한계가 있다. 이를 보완하기 위해 Guo‑Zheng‑Shi(GZS), Bouzidi‑Firdaouss‑Lallemand(BFL), Junk‑Yang(JY) 세 가지 고급 경계조건을 선택하였다. GZS는 힘과 질량 보존을 동시에 만족하도록 설계된 반사‑보정 방식이며, BFL은 경계와 격자 사이의 거리 정보를 이용해 선형 보간을 수행함으로써 2차 정확도를 달성한다. JY는 충돌 연산 단계에서 경계조건을 직접 적용하는 새로운 접근법이지만, 높은 레이놀즈 수에서 수치적 불안정성을 보인다.

연구팀은 오픈소스 혈류 시뮬레이션 플랫폼 HemeLB에 위 세 방법을 구현하고, 다양한 물리적 파라미터(레일놀즈 수 1–300, 워머슬리 수 4–12, 딘 수 100–200)를 갖는 파이프 흐름을 테스트하였다. 특히 파이프 축이 격자와 임의 각도로 배치된 경우를 선택해 경계조건의 일반성을 검증하였다. 정밀한 해석을 위해 Poiseuille 정류 흐름, Womersley 비정상 흐름, 그리고 곡선 파이프의 Dean 흐름에 대한 해석적 해와 비교하였다.

수치 결과는 BFL과 GZS가 모두 2차 공간 수렴을 보이며, 평균 절대 오차가 0.5 % 이하로 감소함을 확인했다. 반면 BB는 1차 수렴에 머물러 격자 해상도가 증가해도 오차 감소율이 제한적이었다. 비정상 흐름(Womersley)에서는 BFL이 GZS보다 약 30 % 낮은 오차를 기록했으며, 계산 시간도 약 15 % 빠른 것으로 나타났다. JY는 레일놀즈 수가 100을 초과할 때 급격히 발산했으며, 안정적인 시뮬레이션을 위해서는 매우 작은 시간 단계와 고해상도가 필요했지만, 전체적인 효율성은 크게 떨어졌다.

또한 충돌 연산자 선택(Lattice Bhatnagar‑Gross‑Krook, LBGK vs. Multiple Relaxation Time, MRT)과 속도 집합(D3Q15, D3Q19, D3Q27)의 변화가 결과에 미치는 영향을 조사했는데, 모든 경우에서 오차와 수렴 속도에 유의미한 차이가 없었다. 이는 복잡한 경계조건이 정확도에 미치는 영향이 충돌 연산자와 속도 집합보다 훨씬 크다는 것을 시사한다.

결론적으로, 복잡한 혈관 기하와 중간 레이놀즈 수, 비정상 흐름을 다루는 LBM 시뮬레이션에서는 BFL이 가장 실용적인 선택이다. BFL은 구현 난이도가 GZS와 비슷하면서도 비정상 흐름에서 더 높은 정확도와 낮은 연산 비용을 제공한다. JY는 특정 조건(극저 레이놀즈 수, 정류 흐름)에서는 사용할 수 있으나, 일반적인 혈류 시뮬레이션에서는 안정성 문제가 크게 작용한다. BB는 간단하지만 정확도가 요구되는 연구에서는 부적합하다.

이 연구는 LBM 기반 혈류 시뮬레이션에서 경계조건 선택이 결과 정확도와 계산 효율성에 미치는 영향을 정량적으로 제시함으로써, 향후 대규모 복잡 혈관 네트워크 모델링에 중요한 지침을 제공한다.