[논문 리뷰] High efficiently numerical simulation of the TDGL equation with reticular free energy in hydrogel

Jun Han, Hui Zhang외 1명|arXiv (Cornell University)|2017. 05. 29.

Material Dynamics and Properties참고 문헌 13인용 수 4

한 줄 요약

이 논문은 망간 자유 에너지가 포함된 시간에 의존하는 기니즈부르크-랜다우(TDGL) 방정식을 사용하여 고분자 마이크로구리 복합체(MMC) 수중고체의 상분리 현상을 시뮬레이션하기 위해 두 가지 $L^2$-안정성 있고 2차 정확도를 갖는 반음성 유한차분 스킴을 제안한다. 적응형 시간스텝 전략이 도입되었으며, 이는 자유 에너지가 급격히 감소하는 초기 단계에는 작은 스텝을, 변화가 미미한 후기 단계에는 큰 스텝을 사용함으로써 최대 80%의 CPU 시간 절감을 이끌어내면서도 높은 정확도를 유지한다.

ABSTRACT

In this paper, we focus on the numerical simulation of phase separation about macromolecule microsphere composite (MMC) hydrogel. The model equation is based on Time-Dependent Ginzburg-Landau (TDGL) equation with reticular free energy. We have put forward two $L^2$ stable schemes to simulate simplified TDGL equation. In numerical experiments, we observe that simulating the whole process of phase separation requires a considerably long time. We also notice that the total free energy changes significantly in initial time and varies slightly in the following time. Based on these properties, we introduce an adaptive strategy based on one of stable scheme mentioned. It is found that the introduction of the time adaptivity cannot only resolve the dynamical changes of the solution accurately but also can significantly save CPU time for the long time simulation.

연구 동기 및 목표

• 망간 자유 에너지를 포함한 TDGL 방정식에 의해 지배되는 MMC 수중고리의 상분리 현상을 시뮬레이션하기 위한 강력한 수치적 스킴 개발.
• $L^2$-노름 안정성 스킴을 통해 수치적 안정성과 에너지 감소 성질 확보.
• 장시간 시뮬레이션의 높은 계산 비용 문제를 해결하기 위해 적응형 시간스텝 전략 도입.
• 수치 실험을 통해 적응형 스킴의 효율성과 정확도 검증.

제안 방법

• MMC 수중고리의 망간 자유 에너지를 기반으로 변분 원리를 사용하여 TDGL 방정식 유도 및 단순화.
• 두 가지 2차 정확도를 갖는 반음성 유한차분 스킴 제안: 하나는 고차원 선형 시스템의 켤레 기울기 해법을 포함하고, 다른 하나는 비선형 시스템에 대한 뉴턴 반복을 사용.
• 총 자유 에너지의 변화율에 따라 시간스텝을 조정하는 적응형 시간스텝 전략 구현 — 자유 에너지가 급격히 감소할 땐 작은 스텝, 변화가 느릴 땐 큰 스텝.
• 스킴의 에너지 감소 성질을 활용하여 후기 시뮬레이션 단계에서 큰 시간스텝을 안전하게 사용할 수 있음을 정당화.
• 고정된 시간스텝과 적응형 시간스텝 간의 CPU 시간, $L^2$-오차, 자유 에너지 변화 곡선을 비교하기 위한 수치 실험 수행.
• 세밀한 고정 시간스텝($\Delta t = 0.001$)과의 비교를 통해 총 자유 에너지 곡선과 $L^2$-오차를 기반으로 적응형 전략의 타당성 검증.

실험 결과

연구 질문

• RQ1 망간 자유 에너지를 포함한 단순화된 TDGL 방정식에 대해 $L^2$-안정성과 2차 정확도를 갖는 반음성 스킴을 구성할 수 있는가?
• RQ2 MMC 수중고리의 상분리 장시간 시뮬레이션에서 정확도를 훼손하지 않고 계산 비용을 어떻게 줄일 수 있는가?
• RQ3 자유 에너지 동역학에 기반한 적응형 시간스텝 전략이 '거칠게-매끄럽게' 변화하는 패턴을 보이는 시스템에서 시뮬레이션 효율성을 향상시키는가?
• RQ4 적응형 대비 고정 시간스텝을 사용할 경우 CPU 시간 절감과 정확도 향상의 정량적 성과는 어떠한가?

주요 결과

• 적응형 시간스텝 전략은 $t = 500$ 시점에서 $\Delta t = 0.001$ 대비 약 80%의 CPU 시간 절감을 이끌어내었으며, 각각 2809.4초 대비 14341.4초를 기록했다.
• 적응형 전략을 사용할 경우 $t = 200$ 시점에서의 $L^2$-오차는 0.000325로, 고정 스텝($\Delta t = 0.05$: 0.0376, $\Delta t = 0.01$: 0.0083)보다 유의미하게 낮았다.
• 적응형 전략 하에서의 총 자유 에너지 변화 곡선은 $\Delta t = 0.001$의 결과와 매우 유사하여, 후기 단계에서 큰 시간스텝을 사용함에도 불구하고 높은 정확도를 유지함을 확인했다.
• 스킴은 에너지 감소 및 $L^2$-안정성 성질을 유지하여, 후기 시뮬레이션 단계에서 큰 시간스텝을 안전하게 사용할 수 있도록 했다.
• 초기 단계에서 상분리가 급격히 발생하고 자유 에너지가 크게 감소하므로, 초기에는 작은 시간스텝이 타당하다.
• 적응형 방법은 후기 단계에서 $\Delta t = 0.001$보다 최대 100배 큰 시간스텝을 사용함으로써 높은 효율성을 달성했으며, 오차 증가율은 극히 미미했다.