粘性項だけ2次にしたら誤差が1.82倍に増えた — semi-implicit Gearスキームが対流項を外挿する理由
·7 min read
BDF2に対流項をN^nのまま組み合わせると1次に落ち、誤差はむしろ増えます。外挿2N^n−N^(n−1)の1行で2次が戻ることをコードで確かめます。
圧力を方程式から消した男 — Chorinの射影法とフラクショナルステップ
·6 min read
予測直後の発散はおよそ1e-2、FFT射影の後は1e-13まで落ちる。その圧力Poissonが1ステップの費用の大半を食う。