[논문 리뷰] A Locally Corrected Multiblob Method with Hydrodynamically Matched Grids for the Stokes Mobility Problem

Anna Broms, Mattias Sandberg외 1명|arXiv (Cornell University)|2022. 07. 22.

Lattice Boltzmann Simulation Studies인용 수 4

한 줄 요약

이 논문은 강체 입자에 대한 3차원 스토크스 이동도 문제를 해결할 때 정확도를 향상시키기 위해 유체역학적으로 일치하는 격자를 갖춘 局소 보정 다블록 방법을 제안한다. 이상적인 입자의 자기상호작용을 일치시키기 위해 블록 격자 구성을 최적화하고, 스토크시안 다이내믹스에서 영감을 얻은 국소 쌍보정을 적용하여 자기상호작용 오차와 근접 영역 오차를 모두 감소시킴으로써, 구형 및 축대칭 막대입자에 대해 조밀한 해상도에서도 고정밀 시뮬레이션을 가능하게 한다.

ABSTRACT

Inexpensive numerical methods are key to enable simulations of systems of a large number of particles of different shapes in Stokes flow. Several approximate methods have been introduced for this purpose. We study the accuracy of the multiblob method for solving the Stokes mobility problem in free space, where the 3D geometry of a particle surface is discretised with spherical blobs and the pair-wise interaction between blobs is described by the RPY-tensor. The paper aims to investigate and improve on the magnitude of the error in the solution velocities of the Stokes mobility problem using a combination of two different techniques: an optimally chosen grid of blobs and a pair-correction inspired by Stokesian dynamics. Optimisation strategies to determine a grid with a certain number of blobs are presented with the aim of matching the hydrodynamic response of a single accurately described ideal particle, alone in the fluid. Small errors in this self-interaction are essential as they determine the basic error level in a system of well-separated particles. With a good match, reasonable accuracy can be obtained even with coarse blob-resolutions of the particle surfaces. The error in the self-interaction is however sensitive to the exact choice of grid parameters and simply hand-picking a suitable blob geometry can lead to errors several orders of magnitude larger in size. The pair-correction is local and cheap to apply, and reduces on the error for more closely interacting particles. Two different types of geometries are considered: spheres and axisymmetric rods with smooth caps. The error in solutions to mobility problems is quantified for particles of varying inter-particle distances for systems containing a few particles, comparing to an accurate solution based on a second kind BIE-formulation where the quadrature error is controlled by employing quadrature by expansion (QBX).

연구 동기 및 목표

• 다블록 방법에서 지배적인 자기상호작용 오차를 줄이기 위해 블록 격자 구성 방식을 최적화하여 이상적인 입자의 유체역학적 반응을 일치시키는 것.
• 스토크시안 다이내믹스에서 영감을 얻은 저비용의 국소 쌍보정 기법을 사용하여 중간 거리 및 가까운 거리에 있는 입자들에 대한 정확도를 향상시키는 것.
• 자기상호작용에서 기초 오차를 최소화함으로써 굵은 블록 이산화로도 강체 입자에 대한 정확한 시뮬레이션을 가능하게 하는 것.
• 이질적인 분산계 및 복잡한 입자 형상에 적용 가능한 다블록 방법의 제어 가능한 정확도 프레임워크를 구축하는 것.
• 스토케스틱 다이내믹스에 적합하고 양의 정부호를 유지하는 보정된 이동도 행렬을 확보하여 향후 대규모 시뮬레이션을 지원하는 것.

제안 방법

• 이deal 입자의 이동도를 일치시키기 위해 블록 격자 기하구조(예: rt-격자, r- 및 t-격자)를 최적화하여 자기상호작용 오차를 최소화하는 것.
• 이동과 회전에 대해 별도로 격자를 최적화한 후 해를 조합하여 자기상호작용 오차를 감소시키는 병합 해법 전략을 사용하는 것.
• 스토크시안 다이내믹스에서 영감을 얻은 국소적이고 저비용의 쌍보정을 적용하여 입자 간의 근접 영역 유체역학적 상호작용을 보정하는 것.
• 오차 정량화를 위해 고정밀 기준으로 사용하기 위해 적분방정식의 제2종 경계적분방정식(BIE) 설정과 전개를 통한 적분법(QBX)을 사용하는 것.
• 입자 수가 적은 시스템(구형 및 막대형)을 대상으로 다양한 입자 간 거리와 구성 조건에서 오차를 평가하기 위한 수치 실험을 수행하는 것.
• 보정된 방법이 기준 해와 대비하여 원거리 및 근접 영역 모두에서 속도 정확도 향상을 정량화하기 위해 검증하는 것.

실험 결과

연구 질문

• RQ1 스토크스 이동도 문제에 대한 다블록 방법에서 자기상호작용 오차를 최소화하기 위해 블록 격자 구성 방식을 어떻게 최적화할 수 있는가?
• RQ2 스토크시안 다이내믹스에서 영감을 얻은 쌍보정 기법은 중간 거리에서부터 가까운 거리에 있는 입자들에 대해 얼마나 속도 오차를 감소시키는가?
• RQ3 이동과 회전에 대해 별도 격자를 사용하는 병합 해법 전략은 기존의 다블록 격자 대비 자기상호작용 오차를 상당히 줄일 수 있는가?
• RQ4 자기상호작용 오차와 다블록 방법에서의 쌍보정 효과 사이의 관계는 어떠한가?
• RQ5 특히 매끄러운 끝부분을 가진 축대칭 막대입자에 대해서는 입자 형상이 다블록 방법의 정확도에 어떻게 영향을 미치는가?

주요 결과

• 특히 막대입자에 적합한 rt-격자 구성 최적화를 통해 자기상호작용 오차를 매우 낮추어, 원거리 시뮬레이션에서 더 이상의 오차 감소가 어려운 오차 바닥값 수준으로 낮출 수 있었다.
• 병합 해법 전략은 이동과 회전에 대해 별도로 격자를 최적화함으로써 자기상호작용 오차를 상당히 감소시키며, 특히 큰 입자 간격에서 효과가 뚜렷했다.
• 쌍보정은 근접 영역 상호작용 오차가 자기상호작용 오차 바닥값을 초과할 경우에만 효과적이며, 따라서 자기상호작용이 이미 잘 일치된 상태일 때 가장 유익하다.
• 粗한 이산화 조건에서 병합 해법에 쌍보정을 적용하면 정확도 향상이 뚜렷하게 나타나며, 이는 낮은 자기상호작용 오차 덕분에 보정 효과가 명확하게 드러나기 때문이다.
• 모든 수치 실험에서 보정된 이동도 행렬이 양의 정부호임을 명시적인 고유값 계산을 통해 확인하여, 이는 스트로크스틱 다이내믹스에 응용 가능한 것을 뒷받침한다.
• 격자가 유체역학적으로 일치하고 쌍보정이 적절히 적용된다면, 블록 해상도가 낮더라도 스토크스 이동도 문제에 대해 정확한 해를 도출할 수 있다.