粘性項だけ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ステップの費用の大半を食う。
[論文レビュー] 30 K でブシネスクが嘘をつくとき — Demou(2018) NOB自然対流と圧力分割
·5 min read
水のΔT=30 Kでは補正|1-ρ0/ρ|の最大が0.0117で、圧力Poissonは定係数のまま残ります。