보존형 Allen‑Cahn 방정식을 위한 개선된 DUGKS

초록

본 논문은 보존형 Allen‑Cahn 방정식(CACE)을 풀기 위해 기존 DUGKS의 미세플럭스 계산이 1차 정확도에 머무르는 문제를 지적하고, 셀 인터페이스에서의 분포함수 값을 2차 다항식 재구성(parabolic reconstruction)으로 보정함으로써 2차 정확도를 확보한 새로운 DUGKS 스키마를 제안한다. 제안 방법은 원형 인터페이스의 대각선 이동, Zaleska 디스크 회전, 원형 인터페이스 변형 등 3가지 벤치마크 실험에서 기존 DUGKS 대비 인터페이스 위치와 형태를 더 정확히 재현함을 보인다.

상세 분석

본 연구는 다상 흐름 시뮬레이션에서 인터페이스를 묘사하는 확산‑계면 접근법 중 하나인 보존형 Allen‑Cahn 방정식(CACE)에 초점을 맞춘다. 기존 Lattice Boltzmann 방식은 격자와 CFL 조건에 제약을 받아 수치 확산이 크게 발생하는 반면, Discrete Unified Gas Kinetic Scheme(DUGKS)은 유한체적 형태와 시간·공간 격자 간의 자유로운 선택이 가능해 연속체 흐름에 유리하다. 그러나 기존 DUGKS는 셀 인터페이스에서의 미세플럭스를 선형 재구성(linear reconstruction)으로 계산함에 따라, 실제 목표인 2차 정확도 대신 1차 정확도에 머무른다. 이는 특히 CACE와 같이 힘(term) 항이 포함된 경우, 혹은 충돌 모델의 1차 모멘트가 보존성을 갖지 않을 때 수치적 오차가 누적돼 인터페이스 위치가 왜곡되는 원인이 된다.

논문은 이러한 문제를 정량적으로 분석하고, 인터페이스를 지나는 특성선 상에서 분포함수 값을 구할 때 2차 다항식 재구성(parabolic reconstruction)을 적용하면, 특성선 끝점에서의 분포함수 값을 정확히 2차 정확도로 복원할 수 있음을 증명한다. 구체적으로, 기존 DUGKS는 셀 중심값을 단순히 선형 보간해 인터페이스 면에 전달했지만, 제안된 방법은 이웃 셀들의 값을 이용해 2차 다항식 계수를 구하고, 이를 통해 인터페이스 면에서의 평균값을 계산한다. 이 과정에서 충돌 항과 외부 힘 항을 트라페조이달(rule)과 중점법(mid‑point)으로 처리해 시간적 정확도 역시 2차로 유지한다.

수치 실험에서는 (1) 원형 인터페이스를 대각선으로 평행 이동시키는 경우, (2) Zaleska 디스크를 회전시키는 경우, (3) 원형 인터페이스를 변형시키는 경우를 선택하였다. 각 실험에서 제안된 DUGKS는 기존 DUGKS 대비 L2 오차가 현저히 감소하고, 인터페이스의 기하학적 특성(곡률, 면적 보존 등)이 더 정확히 유지된다. 특히 회전 및 변형 테스트에서 보존형 Allen‑Cahn 방정식의 질량 보존 특성이 유지되면서도 수치 확산이 크게 억제되는 것을 확인하였다.

이 논문의 주요 기여는 다음과 같다. 첫째, DUGKS에서 미세플럭스 계산이 1차 정확도에 머무르는 근본 원인을 명확히 규명하였다. 둘째, 파라볼릭 재구성을 도입해 인터페이스 면에서의 분포함수 값을 2차 정확도로 복원함으로써 전체 스키마의 정확도를 향상시켰다. 셋째, 보존형 Allen‑Cahn 방정식에 특화된 균형 분포함수와 힘 항을 설계해, 기존 LBE 기반 모델이 갖는 격자·CFL 제약을 탈피하면서도 질량 보존과 인터페이스 형상 유지라는 물리적 요구조건을 만족시켰다. 마지막으로, 제안된 방법이 실제 복잡한 다상 흐름 시뮬레이션에 적용 가능함을 벤치마크 실험을 통해 실증하였다. 향후 연구에서는 3차 재구성, 고차 속도 격자(D3Q27) 적용, 그리고 고마흐수(Ma) 영역에서의 안정성 검증 등을 통해 더욱 일반화된 DUGKS 프레임워크를 구축할 여지가 있다.