Color-gradient lattice Boltzmann modeling of wetting boundary condition on curved solid boundaries

이 논문은 고스트 노드(ghost nodes) 상의 질서 매개변수(order parameters)를 업데이트함으로써 컬러 그래디언트 격자 볼츠만법(color-gradient lattice Boltzmann method) 내 곡면 고체 표면에 대한 젖음 경계 조건을 도입하며, 이는 큰 밀도 및 점도 차이를 효과적으로 처리하고 가짜 전류(spurious currents)를 최소화하면서 정적 및 동적 접촉선 거동을 모두 정확하게 재현하도록 GPU 하드웨어에서 검증된 방식이다.

핵심

당신이 물방울이 곡면(예: 나뭇잎에 떨어지는 빗방울이나 유리잔을 타고 미끄러지는 기포)에 부딪힐 때 어떻게 행동하는지를 시뮬레이션하려고 한다고 상상해 보십시오. 이를 컴퓨터로 구현하기 위해, 과학자들은 **격자 볼츠만 방법(Lattice Boltzmann Method)**이라는 기법을 사용합니다. 이 방법을 아주 작은 타일들로 이루어진 거대한 보이지 않는 격자가 컴퓨터 화면을 덮고 있는 것이라고 생각하십시오. 각 타일은 약간의 "유체" 정보를 담고 있으며, 컴퓨터는 이 타일들을 단계별로 업데이트하여 유체가 어떻게 움직이는지 관찰합니다.

까다로운 부분은 경계 조건(boundary condition), 구체적으로 유체가 고체 벽과 맞닿을 때 어떻게 행동하느냐 하는 것입니다. 현실 세계에서 물은 벽에서 갑자기 멈추는 것이 아니라, 표면이 젖어 있는지(깨끗한 유리처럼) 또는 건조한지(왁스 칠을 한 자동차처럼)에 따라 특정한 각도(접촉각)를 형성합니다.

문제점: 기계 속의 "유령"

컴퓨터 시뮬레이션에서 고체 벽은 매끄러운 선이 아닙니다. 왜냐하면 정사각형 격자 타일로 만들어졌기 때문에 울퉁불퉁하기 때문입니다. 수학적 계산을 성립시키기 위해, 컴퓨터는 유체가 존재하지 않는 고체 벽 내부에서 유체가 무엇을 하고 있는지 알아내야 합니다. 이 가상의 지점들을 **"유령 노드(ghost nodes)"**라고 부릅니다.

이 유령 노드들에게 무엇을 할지 알려주는 기존 방식들에는 몇 가지 결함이 있었습니다:

• 때때로 유체가 아무런 힘 없이 스스로 움직이는 것처럼 보이는 가짜 "유령 전류(spurious velocities)"를 만들어냈습니다.
• 곡면을 처리하는 데 어려움을 겪었으며, 종종 평평한 벽만을 위해 설계된 것처럼 작동했습니다.
• 중립적인 각도(물이 퍼지지도 않고 구슬처럼 맺히지도 않는 상태)를 다루기 위해 별도의 복잡한 수학이 필요할 때가 있었습니다.

해결책: 유령들을 위한 새로운 규칙

저자들은 이 유령 노드들을 위한 더 단순하고 새로운 규칙을 도입했습니다.

비유: 유체가 어떤 "기분"(0은 기체, 1은 액체를 나타내는 색상으로 표현됨)을 가지고 있다고 상상해 보십시오. 현실 세계에서 이 기분은 표면을 가로지를 때 기체에서 액체로 부드럽게 변합니다.

• 기존 방식: 벽 뒤에 서 있는 사람의 기분을 짐작하기 위해 무작위로 추측을 외치는 것과 같았습니다.
• 새로운 방식: 저자들은 만약 벽 바로 바깥쪽(유체 쪽)에 있는 사람의 "기분"을 안다면, 그 부드러운 기분의 곡선을 벽을 통과하여 유령 노드까지 수학적으로 연장할 수 있다는 사실을 깨달았습니다. 그들은 단순히 이렇게 묻습니다: "만약 유체가 이곳에서 45도 각도를 형성하고자 한다면, 그 각도를 만들기 위해 유령 노드의 기분은 무엇이어야 하는가?"

이 새로운 규칙은 마치 매끄러운 다리와 같습니다. 이는 유체 방울의 자연스러운 형태를 고체 벽 바로 위에서부터 벽 내부로 약간까지 연장하여, 물방울이 벽과 이루는 각도가 과학자가 요청한 것과 정확히 일치하도록 보장합니다.

테스트 내용

이 새로운 규칙이 작동함을 증명하기 위해, 저자들은 매우 강력한 컴퓨터 칩(NVIDIA A100 GPU)을 사용하여 여러 시뮬레이션을 실행했습니다:

1. 정적인 물방울: 평평한 판과 곡면인 원기둥 위에 물방울을 놓았습니다. 그리고 물방울이 요청한 정확한 각도로 자리 잡는지 확인했습니다.
   • 결과: 새로운 규칙은 특히 각도가 매우 날카롭거나(구슬처럼 맺히는 경우) 매우 평평할 때(넓게 퍼지는 경우) 기존의 가장 좋은 방법보다 더 정확했습니다.
2. 떠 있는 입자: 기름과 물의 경계면에 떠 있는 원기둥을 시뮬레이션했습니다.
   • 결과: 이들의 방법은 이전보다 수면의 위치를 더 정확하게 계산했습니다.
3. 떨어지는 물방울: 물방울이 떨어져 원기둥에 부딪히며 튀고 퍼지는 모습을 시뮬레이션했습니다.
   • 결과: 물방울은 현실적으로 움직였으며, 새로운 규칙은 유체 내부에 이상한 가짜 움직임을 일으키지 않았습니다.

핵심 요약

• 정확성: 새로운 방법은 곡면을 훨씬 더 잘 처리하며, 벽이 평평하든 둥글든 유체의 각도를 올바르게 유지합니다.
• 안정성: "가짜 노이즈(spurious currents)"를 거의 생성하지 않으므로, 시뮬레이션 속 유체가 훨씬 더 자연스럽게 보입니다.
• 단순성: 접촉각이 정확히 90도(중립)일 때 발생하는 복잡한 수학 문제를 피할 수 있어, 이전 방식들이 겪었던 골칫거리를 해결했습니다.
• 속도: 현대적인 컴퓨터 칩(GPU)과 특정 프로그래밍 스타일을 사용하여 시뮬레이션을 매우 빠르게 실행했습니다. 그들은 약간 덜 정밀한 숫자 형식(단정도 정밀도, single precision)을 사용하는 것이 대부분의 테스트 결과에 영향을 주지 않으면서도 컴퓨터를 두 배 더 빠르게 실행한다는 것을 발견했습니다.

요약하자면, 저자들은 컴퓨터 시뮬레이션이 액체와 고체가 만나는 경계면을 어떻게 처리해야 하는지에 대한 더 나은 "규칙서"를 만들었습니다. 이를 통해 디지털 물방울이 곡면에서도 실제 물방울처럼 보이고 행동하도록 만들었습니다.

기술 요약

기술 요약: 곡면 고체 경계에서의 습윤 경계 조건을 위한 컬러-그레이디언트 격자 볼츠만 모델링

문제 정의
습윤 현상을 정확하게 모델링하는 것은 다공성 매질 흐름, 미세 유체 역학, 콜로이드 박막 건조와 같이 기술적으로 중요한 공정을 시뮬레이션하는 데 매우 중요합니다. 확산 계면(diffuse interface)을 사용하는 격자 볼츠만 방법(LBM) 프레임워크에서 곡면 고체 경계에 습윤 경계 조건을 구현하는 것은 상당한 어려움을 수반합니다. 기존 방식들은 인공적인 질량 소스 또는 싱크(sink)의 도입, 계면 근처에서의 가짜(parasitic) 속도 생성, 평면 기하 구조로의 제한, 그리고 노슬립(no-slip) 조건과 습윤 조건을 동시에 강제할 때 발생하는 불일치 등의 한계를 갖는 경우가 많습니다. 또한, 많은 방법론이 중립 접촉각($\theta = 90^\circ$)에 대해 특별한 처리를 요구하거나, 고체 경계 근처에서 오차를 유발하는 선형 보간법에 의존합니다.

방법론
저자들은 2차원 컬러-그레이디언트 LBM 프레임워크 내에서 곡면 고체 경계를 위해 특별히 설계된 습윤 경계 조건을 소개합니다. 이 방법론의 핵심은 고체 상 내부에 위치한 "고스트 노드(ghost nodes)"에서 질서 매개변수(order parameter, 컬러 필드 $\phi$)를 업데이트하는 것입니다.

• 고스트 노드 업데이트 규칙: 일정한 값을 규정하거나 선형 보간을 사용하는 대신, 이 방식은 평형 컬러 프로파일을 고체 상 내부로 확장합니다. 고스트 노드($x_G$)에서의 컬러 필드 값은 고체 법선 방향을 따라 $\Delta x$만큼 떨어진 유체 영역 내의 대응점($x_F$)에서의 컬러 필드를 기반으로 계산됩니다.
• 수학적 정식화: 컬러 분포가 평형 상태 근에 있다고 가정할 때, 고체 법선 방향의 질서 매개변수 미분은 평형 프로파일과 규정된 접촉각 $\theta$로부터 유도됩니다. 이 관계를 적분하면 다음과 같은 비선형 업데이트 규칙을 얻습니다:
  $$\phi_G = \frac{\phi_F}{\phi_F + (1 - \phi_F) \exp(\varepsilon)}$$
  여기서 $\varepsilon = -4\Delta x \cos \theta / W$ 이며, $W$는 계면 폭입니다.
• 주요 특징: 이 정식화는 고스트 노드의 컬러 필드가 0과 1 사이의 범위 내에 머물도록 보장하여 인공적인 질량 전달을 방지합니다. 결정적으로, 이 규칙은 이전의 방식들과 달리 $\theta = 9^\circ$ (즉, $\varepsilon = 0$)일 때 특별한 처리를 요구하지 않습니다. 이 방법은 커널 퓨전(kernel fusion)과 자동 미분을 활용하여 고성능을 구현한 NVIDIA A100 GPU 기반의 JAX 프레임워크를 사용하여 구현되었습니다.

주요 기여 및 결과
제안된 스킴은 여러 테스트 케이스에 대해 해석적 해(analytical solutions) 및 Fakhari 등[20]의 벤치마크 스킴을 대상으로 검증되었습니다:

1. 평판 위의 정적 액적: 제안된 스킴은 $30^\circ$에서 $150^\circ$ 사이의 범위에서 측정된 접촉각에 대해 Fakhari 등의 방식보다 일관되게 낮은 오차를 나타냈습니다. 또한 벤치마크 스킴에 비해 더 작은 가짜 전류(spurious currents, $5.74 \times 10^{-7}$ ~ $7 \times 10^{-7}$ 격자 단위)를 생성했습니다.
2. 실린더 위의 액적: 곡면 실린더 위에 놓인 액적에 대한 실험에서, 제안된 방법은 특히 극단적인 접촉각($30^\circ$ 및 $150^\circ$)에서 액적 정점 높이($h_{max}$)를 예측하는 데 있어 개선된 정확도를 보여주었습니다(벤치마크 대비 오차가 약 $0.76\%$에서 약 $0.35\%$로 감소).
3. 액체-기체 계면에 고정된 입자: 계면에 갇힌 원통형 배열의 주기적 시뮬레이션에서, 제안된 스킴은 극단적인 각도에서 수직 계면 변위 오차를 벤치마크의 약 $14\%$에서 약 $9\%$로 줄였습니다.
4. 동적 액적 충돌: 본 스킴은 높은 레이놀즈 수와 본드 수에서 실린더에 대한 액적 충돌을 성공적으로 시뮬레이션하였으며, 친수성 표면에서의 액적 퍼짐과 소수성 표면에서의 분리와 같은 복잡한 동적 거동을 포착하였습니다. 결과는 박막 두께에 대한 해석적 스케일링 법칙과 잘 일치하였습니다.
5. 계산 성능: GPU 구현 결과, 단정밀도(FP32) 계산은 배정밀도(FP64)보다 최소 2배 빨랐으나 정적 및 중간 정도의 동적 사례에 대해 정확도를 유지했습니다. 그러나 고에너지 시나리오인 고 레이놀즈/웨버 수 액적 충돌 테스트에서는 FP32가 불안정해졌으며, 이 경우 FP64가 필요함을 확인했습니다.

의의 및 주장
본 논문은 제안된 습윤 경계 조건이 곡면 고체 경계에서의 비혼합성 2상 유동 시뮬레이션을 위한 견고하고 정확한 대안임을 주장합니다. 고스트 노드 업데이트 규칙을 통해 평형 프로파일을 고체 상으로 확장함으로써, 이 방법은 다음과 같은 이점을 제공합니다:

• 밀도 및 점도 차이가 큰 모델을 처리할 수 있는 능력을 유지합니다.
• 상대적으로 작은 가짜 전류를 생성합니다(대중적인 의사 퍼텐셜 모델보다 몇 자릿수 낮은 수준).
• 다른 방법들에서 발견되는 중립 접촉각 관련 수치적 특이성을 피합니다.
• 곡면 기하 구조에서의 정적 및 동적 접촉선을 모두 포함하여 기존 스킴(특히 Fakhari 등)보다 우수한 정확도를 입증합니다.

저자들은 현재 구현이 단일 GPU 실행 및 2차원에 국한되어 있으나, 기초가 되는 물리 법칙은 위상장(phase-field) 및 자유 에너지 모델에 직접 적용 가능하며 3차원으로 확장될 수 있다고 언급했습니다. 본 연구는 특히 접촉선 모세관력 계산 및 이력 현상(hysteresis) 스킴의 발전을 결과적으로 결합할 때, 액체-기체 계면의 고체 입자 시뮬레이션에 관한 토대를 제공합니다.