[论文解读] Benchmark Computation of Morphological Complexity in the Functionalized Cahn-Hilliard Gradient Flow

Andrew Christlieb, Keith Promislow+4 more|arXiv (Cornell University)|Jun 8, 2020

Solidification and crystal growth phenomena被引用 4

一句话总结

本文针对两亲性体系相分离过程中的形态复杂性，引入功能化Cahn-Hilliard（FCH）梯度流的基准计算，采用二阶时间离散化与傅里叶伪谱空间方法。比较了基于BDF2的四种格式——IMEX、SAV、ETD与全隐式PSD，发现SAV与IMEX在计算效率方面表现最佳，而PSD虽计算量更大，但实现了最高精度。

ABSTRACT

Reductions of the self-consistent mean field theory model of amphiphilic molecules in solvent can lead to a singular family of functionalized Cahn-Hilliard energies. We modify these energies, mollifying the singularities to stabilize the computation of the gradient flows and develop a series of benchmark problems that emulate the "morphological complexity" observed in experiments. These benchmarks investigate the delicate balance between the rate of absorption of amphiphilic material onto an interface and a least energy mechanism to disperse the arriving mass. The result is a trichotomy of responses in which two-dimensional interfaces either lengthen by a regularized motion against curvature, undergo pearling bifurcations, or split directly into networks of interfaces. We evaluate a number of schemes that use second order BDF2-type time stepping coupled with Fourier pseudo-spectral spatial discretization. The BDF2-type schemes are either based on a fully implicit time discretization with a PSD nonlinear solver, or upon IMEX, SAV, ETD approaches. All schemes use a fixed local truncation error target with adaptive time-stepping to achieve the error target. Each scheme requires proper "preconditioning" to achieve robust performance that can enhance efficiency by several orders of magnitude.

研究动机与目标

• 开发具有物理解释的基准问题，以复现两亲性体系相分离过程中观察到的形态复杂性。
• 研究两亲性物质在界面处到达速率与能量最小化重分布机制之间的微妙平衡。
• 在自适应时间步长与固定局部截断误差条件下，评估并比较多种二阶时间积分格式在刚性FCH梯度流中的性能。
• 通过测量达到目标全局误差所需快速傅里叶变换（FFT）调用次数，实现计算效率的公平比较。
• 证明能量衰减必须与精度相结合，因为在高度非线性系统中，仅能量稳定性可能掩盖错误结果。

提出的方法

• 提出一种光滑化、非奇异的功能化Cahn-Hilliard（FCH）能量形式，以稳定梯度流的数值计算。
• 采用二阶后向差分公式（BDF2）时间离散化，结合傅里叶伪谱空间离散化，以实现高精度与谱收敛性。
• 实现四种数值格式：(1) 全隐式预处理最速下降法（PSD），(2) 线性隐式IMEX，(3) 标量辅助变量（SAV），(4) 指数时间差分（ETD）用于对比。
• 采用自适应时间步长策略，固定局部截断误差目标，确保所有格式的一致精度。
• 在所有格式中引入预处理技术，以增强对刚性非线性系统的鲁棒性与效率。
• 采用SAV方法确保无条件修正能量稳定性，并在格式下可证明能量衰减特性。

实验结果

研究问题

• RQ1 在自适应时间步长下，不同二阶时间积分格式在解析FCH梯度流形态复杂性方面表现如何？
• RQ2 在使用不同数值格式求解刚性FCH系统时，精度与计算成本（以FFT调用次数衡量）之间的权衡如何？
• RQ3 仅能量稳定性是否足以确保准确的模拟结果？在复杂相场系统中，精度是否必须与能量衰减相辅相成？
• RQ4 两亲性物质的输送速率如何影响界面形态的选择——如曲率驱动运动、珠串化（pearling）或网络形成？
• RQ5 预处理在FCH模拟中对非线性和线性隐式格式的性能提升程度如何？

主要发现

• 全隐式BDF2-PSD格式在固定局部截断误差下实现了最小的全局离散化误差，表明其在所测试格式中精度最高。
• IMEX与SAV格式计算效率最高，达到相同全局误差水平时，所需快速傅里叶变换（FFT）调用次数显著少于PSD。
• SAV格式性能与IMEX极为接近，仅因辅助变量系统导致FFT调用次数增加约1.4倍。
• ETD格式作为对比引入，但在计算效率方面仍逊于IMEX与SAV方法。
• 预处理对实现稳健性能至关重要，在所有格式中均实现数个数量级的效率提升。
• 基准问题成功复现了实验中观察到的三类响应行为——曲率驱动运动、珠串化分岔及曲线分裂成网络结构。