[논문 리뷰] Implementation of on-site velocity boundary conditions for D3Q19 lattice Boltzmann

Martin Hecht, Jens Harting|TU/e Research Portal (Eindhoven University of Technology)|2008. 11. 27.

Lattice Boltzmann Simulation Studies인용 수 10

한 줄 요약

이 논문은 3차원 D3Q19 격자 보울츠만 방법에 대해 국소적으로 적용 가능한 현장 속도 경계 조건을 제안하며, 격자 정렬과 무관하게 임의의 유량 방향을 허용한다. 이 방법은 2차 정확도를 보장하고 수치적 슬립을 제거하며, 속도를 0으로 설정함으로써 자연스럽게 완전고정 조건을 포함한다. 다양한 시험 케이스를 통한 검증을 통해 복잡한 기하구조에서 높은 정확도와 강건성을 입증하였다.

ABSTRACT

On-site boundary conditions are often desired for lattice Boltzmann simulations of fluid flow in complex geometries such as porous media or microfluidic devices. The possibility to specify the exact position of the boundary, independent of other simulation parameters, simplifies the analysis of the system. For practical applications it should allow to freely specify the direction of the flux, and it should be straight forward to implement in three dimensions. Furthermore, especially for parallelized solvers it is of great advantage if the boundary condition can be applied locally, involving only information available on the current lattice site. We meet this need by describing in detail how to transfer the approach suggested by Zou and He to a D3Q19 lattice. The boundary condition acts locally, is independent of the details of the relaxation process during collision and contains no artificial slip. In particular, the case of an on-site no-slip boundary condition is naturally included. We test the boundary condition in several setups and confirm that it is capable to accurately model the velocity field up to second order and does not contain any numerical slip.

연구 동기 및 목표

• 3차원 D3Q19 격자 보울츠만 방법에 대해 일반적이고 국소적으로 적용 가능한 속도 경계 조건을 개발하여, 임의의 유입 방향을 허용한다.
• 수치적 슬립을 제거하고 속도장 재구성에서 2차 정확도를 확보한다.
• BGK 감쇠 시간에 의존하지 않는 경계 조건을 제공하여 다공성 매체 및 마이크로플루이딕 장치의 시뮬레이션에서 정확한 벽 위치를 가능하게 한다.
• Zou와 He의 접근법을 D2Q9에서 D3Q19로 일반화하여, 정상이 아닌 유량에 대해 횡방향 동량 보정을 포함한다.
• 카르테시안 평면과 일치하지 않는 벽을 포함한 복잡한 기하구조, 예를 들어 마이크로믹서나 다공성 물질에서의 정확한 시뮬레이션을 가능하게 한다.

제안 방법

• D2Q9에 대한 Zou와 He(1997)의 접근법을 D3Q19 격자로 확장하여, 임의의 유입 방향을 갖는 유량 경계 조건을 유도한다.
• 경계 노드의 분포 함수를 알려진 밀도와 속도 성분을 사용하여 국소적 복구 단계를 적용한다.
• 비평형 부분에 대해 반사 규칙을 적용하고, 동역학적 모멘텀 균형을 유지하기 위해 $ N^{x}_{y}, N^{x}_{z}, N^{y}_{x}, N^{y}_{z} $ 항을 통해 횡방향 동량 보정을 추가한다.
• 유입/유출 경계에서의 밀도 및 분포 함수 표현식을 속도 성분과 이웃 분포 함수 기반으로 유도한다.
• 속도 성분과 국소 분포 값에 의존하는 보정 항을 포함하여 평형 분포와의 일致성을 확보한다.
• 유입 및 유출 경계 모두에 대칭적 형태를 사용하며, 각 면(예: x=0, x=max, y=0, y=max)에 대해 별도의 식을 설정하여 국소 적용 가능성을 보장한다.

실험 결과

연구 질문

• RQ1 Zou와 He 경계 조건 접근법은 3차원 D3Q19 격자 보울츠만 시뮬레이션에 대해 임의의 유량 방향으로 일반화될 수 있는가?
• RQ2 제안된 방법은 수치적 슬립 없이 속도장 재구성에서 2차 정확도를 달성하는가?
• RQ3 현재 격자 위치의 정보만을 사용하여 국소적으로 적용 가능한가? 이는 대규모 시뮬레이션에서 효율적 병렬 처리를 가능하게 하는가?
• RQ4 벽이 격자 축과 일치하지 않는 복잡한 기하구조에서 이 방법은 어떻게 성능을 발휘하는가?
• RQ5 속도를 0으로 설정함으로써 완전고정 조건이 자연스럽게 포함되며, BGK 감쇠 시간에 의존하지 않는가?

주요 결과

• 제안된 경계 조건은 수치 수렴 테스트를 통해 속도장 재구성에서 2차 정확도를 확보하였다.
• 적절한 횡방향 동량 보정 덕분에, 정상이 아닌 유량 방향에서도 수치적 슬립이 발생하지 않았다.
• 국소적 적용이 가능하여 현재 격자 위치의 정보만으로도 충분하며, 대규모 시뮬레이션에서 효율적 병렬 처리를 지원한다.
• 속도 벡터를 0으로 설정함으로써 완전고정 조건이 자연스럽게 포함되며, 벽 위치가 BGK 감쇠 시간에 영향을 받지 않는다.
• 채널의 임의의 정렬 방향을 가진 흐름 및 다공성 매체를 포함한 다양한 설정에서의 검증을 통해 다양한 기하구조에서 강건성과 정확도를 입증하였다.
• 경계에서의 분포 함수 표현식(예: 식 51–54, 65–68)은 격자 보울츠만 평형과 동역학 모멘텀 보존을 보장한다.