[논문 리뷰] Accurate evaluation of integrals in slender-body formulations for fibers in viscous flow

Anna‐Karin Tornberg|arXiv (Cornell University)|2020. 12. 23.

Micro and Nano Robotics참고 문헌 9인용 수 4

한 줄 요약

이 논문은 점탄성 섬유에서의 비국소 섬유체 이론에 나타나는 유한부 및 거의 특이 적분을 고정밀도로 평가하기 위한 수치적 방법을 제시한다. 제품 적분과 사전 계산된 수정된 구분법 가중치를 조합함으로써, 추가 비용 없이 특이 적분에 대해 스펙트럼 정밀도를 달성하고, 관련 기법을 통해 섬유 표면 근처에서 높은 정밀도로 거의 특이 스토크스렛 적분을 처리한다.

ABSTRACT

A non-local slender body approximation for slender flexible fibers in Stokes flow can be derived, yielding an integral equation along the center lines of the fibers that involves a slenderness parameter. The formulation contains a so-called finite part singular integral, and can in the case of several fibers or evaluation of the flow field require the evaluation of nearly singular integrals. We introduce a numerical technique to accurately and efficiently evaluate the finite part integral. This technique can be applied combined with any panel based quadrature rule and will add no additional cost except for a small precomputation of modified quadrature weights. We also show how a related technique that was recently introduced can be applied for the evaluation of the nearly singular integrals.

연구 동기 및 목표

• 비국소 섬유체 이론에서 섬유의 스토크스 흐름에 대해 유한부 적분을 평가하기 위한 강력하고 효율적인 수치적 방법을 개발하는 것.
• 섬유가 서로 가까이 있을 때 또는 평가 점이 섬유 중심선 근처일 경우 특이하거나 거의 특이한 적분을 정확히 계산하는 데 도전하는 것.
• 특이성을 분석적으로 다루는 구분 기법을 제공하여, 추가 계산 비용을 초과하지 않는 한도에서 높은 정밀도를 보장하는 것.
• 특히 섬유 표면 근처에서 발생하는 유체 속도 평가에 나타나는 거의 특이 적분에 대해 이 방법을 확장하는 것.
• 나선형과 원형 고리와 같은 벤치마크 기하 구조에서의 검증을 통해 기존의 구분법에 비해 뛰어난 수렴성과 오차 제어 성능을 입증하는 것.

제안 방법

• 유한부 적분을, 매끄러운 다항식 인수와 특이 핵심으로 분리하여 각 패널에 대해 해석적으로 통합할 수 있도록 재구성한다.
• 매끄러운 인수의 다항식 근사가 사용되며, 계수는 사전 계산 단계에서 한 번만 풀어내는 전치 바르멘드 시스템을 통해 결정되어 수정된 구분법 가중치를 생성한다.
• 수정된 구분법 가중치는 모든 패널에 재사용되며, 평가 시 추가 비용 없이 효율적이고 정확한 특이 적분 평가를 가능하게 한다.
• 거의 특이 스토크스렛 적분의 경우, 유사한 제품 적분 접근법이 적용되며, 근을 찾는 과정에서 유도된 복소수 쌍을 사용해 인수를 분리한다.
• 해결된 적분의 해석적 평가에는 재귀 공식이 사용되며, 정밀도를 유지하기 위해 철저한 수치적 구현이 이루어진다.
• 이 방법은 가우스-레지앙드어 구분법 규칙을 사용하여 구현되었으며, 고정밀 기준 솔루션과 함께 3차원 섬유 기하 구조에서 검증되었다.

실험 결과

연구 질문

• RQ1 적분의 피적분함수가 대각선에서 특이성이 있을 경우 섬유체 이론의 유한부 적분을 어떻게 고정밀도로 평가할 수 있는가?
• RQ2 특이성을 분석적으로 다루면서도 효율성을 유지하고 추가 계산 비용을 피할 수 있는 구분 방법을 설계할 수 있는가?
• RQ3 평가 점이 섬유 중심선에 가까울 경우, 제안된 방법은 거의 특이 적분에 대해 어떻게 성능을 발휘하는가?
• RQ4 다양한 섬유 기하 구조와 밀도 함수에 대해 이 방법의 수렴 특성은 어떻게 되는가?
• RQ5 이 방법은 3차원 스토크스 흐름에서 자유단을 가진 다수의 상호작용 섬유에 적용될 수 있는가?

주요 결과

• 제안된 방법은 사전 계산된 수정된 구분법 가중치를 통해 특이성을 분석적으로 다룸으로써, 유한부 적분에 대해 스펙트럼 정밀도를 달성한다.
• 나선형 섬유 근처에서의 유체 속도 평가에서 최대 오차는 동일한 패널 수를 가진 표준 가우스-레지앙드어 구분법에 비해 수개의 주기로 감소한다.
• $ M = 8 $ 패널일 경우, 특수 구분법은 섬유 경계까지 낮은 오차를 유지하지만, 표준 구분법은 특이성 근처에서 오차를 줄일 수 없다.
• 이 방법은 빠른 수렴을 보이며, 최대 오차는 패널 수 $ M $ 증가에 따라 급격히 감소한다. 특히 매끄러운 힘 밀도 함수에 대해 뚜렷한 경향을 보인다.
• 거의 특이 적분을 위한 기법은 평가 점이 섬유에 매우 가까이 있을 경우에도 효과적으로 작동하며, 거리가 $ 2.2 \times 10^{-3} $ 이하일지라도 오차는 여전히 작게 유지된다.
• 사전 계산 단계는 $ n = 16 $ 인 $ n $ 개의 전치 바르멘드 시스템을 $ n \times n $ 크기로 풀고, $ n^2 $ 개의 수정된 가중치를 저장함으로써 수행되며, 이는 모든 패널에 대해 재사용되며 런타임 비용 없이 가능하다.