[論文レビュー] A fourth-order maximum principle preserving operator splitting scheme for three-dimensional fractional Allen-Cahn equations

Dongdong He, Kejia Pan+1 more|arXiv (Cornell University)|Apr 19, 2018

Solidification and crystal growth phenomena参考文献 1被引用数 4

ひとこと要約

本稿では、三次元空間分数マクスウェル＝アレン＝カーン方程式に対する4次精度最大原理保存型オペレータ分割スキームを提案する。ストラング分割、分数拡散項におけるクランク＝ニコルソン交替方向陽解法（ADI）、4次精度有限差分、リチャードソン補外を組み合わせることで、時間刻みの制約下においても無条件安定性を達成し、離散的最大原理を保存する。2次元および3次元の数値実験により、4次精度および最大原理の満たし方が確認された。

ABSTRACT

In this paper, by using Strang's second-order splitting method, the numerical procedure for the three-dimensional (3D) space fractional Allen-Cahn equation can be divided into three steps. The first and third steps involve an ordinary differential equation, which can be solved analytically. The intermediate step involves a 3D linear fractional diffusion equation, which is solved by the Crank-Nicolson alternating directional implicit (ADI) method. The ADI technique can convert the multidimensional problem into a series of one-dimensional problems, which greatly reduces the computational cost. A fourth-order difference scheme is adopted for discretization of the space fractional derivatives. Finally, Richardson extrapolation is exploited to increase the temporal accuracy. The proposed method is shown to be unconditionally stable by Fourier analysis. Another contribution of this paper is to show that the numerical solutions satisfy the discrete maximum principle under reasonable time step constraint. For fabricated smooth solutions, numerical results show that the proposed method is unconditionally stable and fourth-order accurate in both time and space variables. In addition, the discrete maximum principle is also numerically verified.

研究の動機と目的

• 三次元空間分数マクスウェル＝アレン＝カーン方程式に対して、離散的最大原理を保存する高次精度数値スキームの開発。
• 時間・空間方向に4次精度を達成するとともに、無条件安定性を維持すること。
• 妥当な時間刻み制約下で、数値解が[−1, 1]の範囲に有界に保たれることを保証すること。
• 非線形ソルバーを回避し、エネルギー安定性のトレードオフを避ける線形で効率的かつ無条件安定な手法を提供すること。
• 離散的最大原理の保存およびスキームの高次収束性を数値的に検証すること。

提案手法

• ストラングの2次精度オペレータ分割を用いて、分数マクスウェル＝アレン＝カーン方程式を3つの部分問題に分解する。そのうち2つは解析的常微分方程式を含み、残り1つは3次元線形分数マクスウェル＝アレン＝カーン拡散方程式を含む。
• 3次元分数拡散項は、クランク＝ニコルソン交替方向陽解法（ADI）を用いて解かれる。これにより多次元問題が1次元問題の逐次解法に還元される。
• 空間分数微分を高次精度で離散化するために4次精度有限差分法を用いる。
• 時間積分の精度を4次に高めるためにリチャードソン補外を適用する。
• フーリエ解析を用いて、滑らかな解に対して全体のスキームが無条件安定であることを証明する。
• 離散的最大原理が数値的に保存されるように、時間刻みの制約を導出する。

実験結果

リサーチクエスチョン

• RQ1 3次元空間分数マクスウェル＝アレン＝カーン方程式に対して、離散的最大原理を保存する高次（4次）精度数値スキームを構築可能か？
• RQ2 分数マクスウェル＝アレン＝カーン方程式の存在下で、時間・空間方向に両方とも高次精度を達成しつつ、無条件安定性を維持することは可能か？
• RQ3 提案された線形オペレータ分割スキームは、妥当な時間刻み制約下で解の境界[−1, 1]を保つことができるか？
• RQ4 分数マクスウェル＝アレン＝カーン方程式の分数マクスウェル＝アレン＝カーン指数αが、相分離シミュレーションにおける界面厚さおよび粗大化ダイナミクスにどのように影響するか？
• RQ5 リチャードソン補外は、安定性や最大原理保存性を損なわずに、分割スキームの時間精度を有効に向上させることができるか？

主な発見

• 提案されたスキームは、フーリエ解析により、滑らかな解に対して無条件安定であることが証明された。
• 数値収束テストにより、人工的に生成された滑らかな解を用いて、時間および空間変数において4次精度が確認された。
• 離散的最大原理は数値的に検証された：時間刻みが導出された制約（例：α=1.7の場合Δt ≤ 0.4945）を満たす限り、解は[−1, 1]の範囲に有界に保たれる。
• 時間刻みが安定性制約を超える場合（例：Δt = 2）、解の最大値が1を超えることが観測され、離散的最大原理が破られる。
• 分数マクスウェル＝アレン＝カーン指数αを小さくすると、界面が薄くなり、相分離の粗大化が遅くなる。特にα=1.2では、最も不均質で細分化された相構造が得られる。
• 本スキームは、界面の進化および粗大化を高精度で捉え、有界な解を保つことに成功した。