A volume penalization method for solving conjugate scalar transport with interfacial jump conditions

기술 요약: 계면 점프 조건이 있는 공액 스칼라 수송을 위한 체적 페널티 방법

문제 정의
계면 점프 조건(예: 미세 채널 열전달, 계면에서의 화학 반응, 상변화 시스템 등)을 포함하는 공액 스칼라 수송(예: 열 또는 질량 전달)은 열 및 화학 공정에서 널ily하게 나타난다. 이러한 시나리오에서는 계면에서의 스칼라 값 자체(예: 온도 또는 농도 점프) 및/또는 그 플럭스(flux)의 불연속성이 존재할 수 있다. 체 바디 적합 격자(Body-Fitted Mesh, BFM) 시뮬레이션은 별도의 상 영역을 명시적으로 결합함으로써 이러한 조건을 처리할 수 있지만, 복잡한 기하학적 구조를 위한 격자 생성이 까다롭다는 단점이 있다. 임머시드 경계법(Immersed Boundary Methods, IBMs), 특히 체적 페널티 방법(Volume Penalization Method, VPM)은 복잡한 격자 생성 없이 임의의 기하학적 구조를 카테시안 격자에 삽입할 수 있는 방법을 제공한다. 그러나 기존의 VPM 방식은 노이만 경계 조건(Neumann boundary conditions)과 일반적인 계면 점프 조건을 처리하는 데 어려움이 있다. 전통적인 접근 방식은 고체 영역 내에 비물리적인 스칼라 장을 생성하거나, 추가적인 복잡한 결합 절차 없이는 플럭스 점프나 스칼라 불연속성을 정확하게 강제하지 못하는 경우가 많다.

방법론
저자들은 일반적인 계면 점프 조건을 가진 공액 스칼라 수송을 해결하기 위해 VPM 프레임워크 내에서 새로운 처리를 제안한다. 방법론은 다음과 같이 세 가지 주요 단계로 진행된다:

1. 레벨 셋 표현: 임의의 유체-고체 계면을 레벨 셋 함수($\phi_0$)로 표현하고, 이를 두께 $2\delta$의 얇은 계면 영역을 통해 0(유체)에서 1(고체)로 전이되는 상 지표($\phi$)로 변환한다.
2. 노이만 경계 조건 처리: 노이만 조건을 가진 단일 상 이송-확산 문제의 경우, 저자들은 지배 방정식에 추가될 발산 형태(divergence form)의 소스 항 $\nabla \cdot (q_w \phi \mathbf{n})$을 유도한다. 고체 체적 전체에 걸쳐 페널티 항을 통해 미분값을 강제하는 기존 방식과 달리, 이 새로운 항은 계면 영역 내에서만 0이 아닌 값을 갖는다. 이는 규정된 스칼라 플럭스 $q_w$가 정확하게 적용되도록 보장하는 동시에, 고체 내부의 분자 확산율을 0으로 유지하여 고체 상으로의 비물리적인 확산을 방지한다.
3. 점프가 있는 공액 수송으로의 일반화: 스칼라 점프($c_f = \alpha c_s$)와 플럭스 점프를 모두 처리하기 위해, 저자들은 "등가 스칼라(equivalent scalar)" $h = \alpha c$를 도입한다. 이 변환은 스칼라 불연속성을 제거하여 두 상 모두에 대해 통일된 지배 방정식을 사용할 수 있게 한다. 플럭스 점프 조건은 스칼라 점프의 효과와 규정된 플럭스 점프를 모두 포함하는 등가 플럭스 $q^*_w$로 재구성된다. 등가 스칼라 $h$에 대한 통일된 지배 방정식은 다음과 같다:
   $$\frac{\partial h}{\partial t} + \mathbf{u} \cdot \nabla h = \nabla \cdot (D \nabla h) + \nabla \cdot (q^*_w \phi \mathbf{n})$$
   여기서 $D$는 상 의존적 확산율이다.

이 스케마는 오픈 소스 유한 체적 코드인 OpenFOAM (버전 8)에 구현되었으며, 공간 항에는 2차 선형 이산화(second-order linear discretization)를, 과도 항에는 오일러 암시적(Euler implicit) 방식을 사용하였다.

주요 기여

• 새로운 소스 항 공식화: 본 논문은 계면 플럭스 점프를 정확하게 강제하면서 고체 영역 내에서 비물리적인 해를 생성하지 않는 발산 형태의 소스 항을 도입한다.
• 통일된 지배 방정식: 등가 스칼라 및 플럭스를 정의함으로써, 유체 및 고체 상에 대한 지배 방정식을 통일하여 명시적인 계면 결합이나 별도의 도메인 솔루션 없이도 문제를 해결한다.
• 일반 점프 조건으로의 확장: 스칼라 값과 플럭스의 동시 점프를 처리할 수 있도록 VPM을 성공적으로 확장하였으며, 이는 기존의 유한 체적 기반 임머시드 경계 스킴에서 부족했던 기능이다.
• 오차 분석: 저자들은 이 스케마가 계면 곡률이 존재하는 경우에도 격자 크기에 대해 1차 정확도를 달in함을 보여주는 분석적 오차 분석을 제공한다. 이는 기존 VPM의 특성과 일치한다.

결과
제안된 방법은 세 가지 뚜로 된 테스트 케이스를 통해 검증되었다:

1. 1차원 확산: 노이만 경계 조건을 가진 과도 확산 문제가 해결되었다. 결과는 해석적 해와 잘 일치하였으며, 평균 상대 오차는 약 1.25%였다. 결정적으로, 기존 VPM 스킴과 달리 본 방법은 고체 내부에서 0의 스칼라 장을 유지하여 비물리적인 확산을 방지하였다. 오차는 확산율 변화에 관계없이 안정적이었던 반면, 기존 스킴은 확산율이 높아짐에 따라 오차가 증가하는 모습을 보였다.
2. 유체-고체 결합 확산: 스칼라 및 플럭스 점프가 모두 존재하는 사각형 도메인 내 실린더 시뮬레이션이 수행되었다. VPM 결과는 chtMulti_regionFoam을 사용한 체 바디 적합 격자(BFM) 시뮬레이션과 비교되었다. 확산율 비와 점프 크기가 변하는 4가지 테스트 케이스 전반에 걸쳐 VPM과 BFM 결과 사이의 상대적 편차는 3.0% 이내(평균 1.46%)였다.
3. 유체-고체 결합 이송-확산: 계면 점프가 있는 실린더 주변의 정상 상태 층류 흐름($Re=40$)이 시뮬레이션되었다. VPM 결과는 BFM 벤치마크와 다시 한번 강력한 일치를 보였으며, 상대적 편차는 1.59%였다. 이 방법은 계면 근처의 스칼라 분포와 플럭스 점프를 정확하게 재현하였다.

의의 및 주장
본 논문은 제안된 VPM이 복잡한 계면 기하학 및 점프 조건을 가진 공각 스칼라 수송 문제에 대해 체 바디 적합 격자를 대체할 수 있는 강력하고 효율적인 대안을 제공한다고 주장한다. 주요 의의는 단일 통일 지배 방정식을 사용하여 결합된 유체-고체 문제를 해결하면서 국부적인 점프 조건을 정확하게 강제할 수 있다는 점에 있다. 저자들은 이 접근 방식이 노이만 조건에 대한 전통적인 VPM 구현에서 흔히 발생하는 비물리적 아티팩트를 피할 수 있다고 언급한다.

저자들은 이 연구를 기초적인 단계로 겸손하게 위치시키며, 복잡한 기하학적 구조를 다룰 수 있는 이 방법의 능력이 거친 벽면(rough walls) 및 복잡한 열전달 표면에 대한 향후 응용에 적합함을 시사한다. 또한, 이 VPM을 복사, 난류 및 공액 열전달 문제에서의 형상 및 위상 최적화를 위한 어드조인트(adjoint) 방법과 결합할 수 있는 잠재력을 강조하였으나, 이러한 구체적인 응용은 현재 연구의 즉각적인 결과라기보다 향후 연구 방향으로 식별되었다.