[논문 리뷰] Numerical simulation of solutions and moments of the smoluchowski coagulation equation

Dustin D. Keck, David M. Bortz|arXiv (Cornell University)|2013. 12. 27.

Advanced Mathematical Modeling in Engineering인용 수 4

한 줄 요약

이 논문은 유한요소법(FEM)과 유한체적법(FLFM)이 스몰وخ프스키 응집방정식을 수치적으로 해결하는 데 대해 정확도와 계산 비용을 비교한다. 솔루션의 모멘트 근사값과 계산 비용에 중점을 두며, 느린 응집 시스템에서는 FEM이 훨씬 낮은 계산 비용으로 뛰어난 정확도를 보이며, 빠른 응집에서는 제0차 모멘트에서 FLFM이 약간 더 정확하고, 제1차 모멘트에서는 일관되게 더 정확하다. 두 방법 모두 제2차 수렴성을 보인다.

ABSTRACT

Researchers have employed variations of the Smoluchowski coagulation equation to model a wide variety of both organic and inorganic phenomena and with relatively few known analytical solutions, numerical solutions play an important role in studying this equation. In this article, we consider numerical approximations, focusing on how different discretization schemes impact the accuracy of approximate solution moments. Pursuing the eventual goal of comparing simulated solutions to experimental data, we must carefully choose the numerical method most appropriate to the type of data we attain. Within this context, we compare and contrast the accuracy and computational cost of a finite element approach and a finite volume-based scheme. Our study provides theoretical and numerical evidence that the finite element approach achieves much more accuracy when the system aggregates slowly, and it does so with much less computation cost. Conversely, the finite volume method is slightly more accurate approximating the zeroth moment when the system aggregates quickly and is much more accurate approximating the first moment in general. Lastly, our study also provides numerical evidence that the finite element method (conventionally considered first order) actually belongs to a class of discontinuous Galerkin methods that exhibit superconvergence, or second order in our case.

연구 동기 및 목표

• 스몰وخ프스키 응집방정식을 해결하는 데 있어 유한요소법(FEM)과 유한체적법(FLFM)의 정확도와 계산 비용을 비교하는 것.
• 이산화 방법이 다양한 응집 역학 조건 하에서 솔루션 모멘트, 특히 제0차 및 제1차 모멘트 근사에 미치는 영향을 평가하는 것.
• 절단 파라미터 $x_{\text{max}}$와 격자 간격 $\Delta x$가 두 방법의 수치 정확도에 미치는 영향을 조사하는 것.
• 예를 들어 쿨터 카운터(제0차 모멘트) 또는 플로우 사이토메트리(제1차 모멘트) 데이터와 같은 실험 데이터 유형에 따라 방법 선택을 안내하는 것.
• FEM의 수렴 속도 및 초수렴 현상에 대한 이론적 및 수치적 증거를 제공하는 것.

제안 방법

• 연구는 Banks와 Kappel의 이전 작업을 기반으로 하되, Ackleh와 Fitzpatrick에 의해 확장된 유한요소법(FEM)을 사용하여 응집방정식의 공간 이산화를 수행한다.
• 비교를 위한 기준으로, $L^1$에서 제2차 수렴성을 보이는 것으로 알려진 Filbet와 Laurençot의 유량 방법(FLFM)이라는 유한체적 유형의 방법을 사용한다.
• 수치적 해는 상수 핵심 $K_A(x,y) \equiv 1$과 곱셈 핵심 $K_A(x,y) = xy$의 두 가지 응집 핵심에 대해 계산된다.
• FEM과 FLFM의 정확도는 미세 격자 기준 솔루션과 알려진 모멘트 값과의 비교를 통해 평가된다.
• 감도 및 수렴 행동 평가를 위해 $\Delta x$와 $x_{\text{max}}$를 변화시킨 격자 정밀도 연구가 수행된다.
• 모멘트 근사는 $M_i(f; x_1,x_2) = \int_{x_1}^{x_2} x^i f(t,x) \, dx$를 통해 계산되며, 여기서 $f$는 크기 분포를, $g = x f$는 부피 분포를 나타낸다.

실험 결과

연구 질문

• RQ1 스몰وخ프스키 응집방정식의 전체 솔루션 근사를 고려할 때 FEM과 FLFM의 정확도는 어떻게 비교되는가?
• RQ2 느린 및 빠른 응집 역학 조건 하에서 제0차 및 제1차 모멘트 근사에 대해 어느 방법이 더 정확한가?
• RQ3 실험 데이터가 작은 부피 범위로 제한된 경우, 절단 파라미터 $x_{\text{max}}$는 각 방법의 정확도에 어떤 영향을 미치는가?
• RQ4 수치적 증거에 따르면 초수렴을 보이지만, 기존에 제1차 수렴으로 간주되는 FEM은 실제로 초수렴을 보이는가?
• RQ5 동일한 정확도 수준을 확보할 때, 어느 방법이 더 높은 계산 효율성을 제공하는가?

주요 결과

• FEM은 $L^1[\mathbf{X}]$에서 제2차 수렴성을 보이며, 이는 FEM이 기존에 제1차 수렴으로 간주되지만 초수렴을 보일 수 있음을 시사하는 가설을 지지한다.
• 느린 응집 시스템($K_A(x,y) \equiv 1$)에서는 FEM이 FLFM와 유사한 솔루션 정확도를 달성하지만, 계산 비용은 크게 낮다.
• 느린 응집 조건에서는 FEM이 제0차 모멘트를 FLFM보다 더 정확하게 근사하며, 제1차 모멘트에서는 FLFM가 약간 더 정확하다.
• 빠른 응집($K_A(x,y) = xy$) 조건에서는 FLFM가 모든 경우에서 제1차 모멘트 근사에 더 정확하고, 제0차 모멘트에서도 약간 더 정확하다.
• 특히 곱셈 핵심 하에서 FEM은 FLFM보다 $x_{\text{max}}$ 선택에 덜 민감하여, 실험 데이터의 상한 탐지 한계가 제한된 경우 더 강건하다.
• 동일한 정확도를 확보할 때, FEM은 800개 격점에서 FLFM가 100개 격점에서 수행하는 것보다도 상당한 계산 비용 절감을 제공한다.