- 遊びで1Dで浮遊砂のみの河床変動計算を書いていたのですが，射流時に飛んじゃいます．
- 正確に言うと，初期に一定勾配を与えた場合は，E=Dとなるまで河床が低下し，その状態の平衡状態が形成されます．
- 一方，河床に凹凸を与えた場合，一定時間経過後，濃度に振動が始まり，それに併せて，河床高および水面に振動が始まります．下図の感じです．
- 色々探ってみると，濃度の保存則の移流項が原因な感じです．確かに濃度の勾配は大きいと思いますが，そんなに非線形性は強いと思わないですが...
- 可能な限り陰的に取り扱ってみたのですがダメでした．原因不明です．
- 濃度のような移流拡散方程式でそんな話って聞いたことあります?完全な陰解法にしないと駄目ですかね．
- 一応コードを置いておきますので興味があれば．VBAです．グチャグチャです．
LuaLaTeXのサンプル
作ったので置いておきます．
windows環境ですいません．
