Matlab で物理的な流れの問題 (パイプ内の境界層の流れ) を解決しようとしています。そのためには、次のような ODE 系を解かなければなりません。
function dy = ODEFUN(t,y)
dy = dp + c;
end
問題は、dp も不明であり、次のステップでいくつかの物理的条件を満たすように設定する必要があることです。dp の明示的な方程式はありませんが、次の時間ステップ t+1 で正しいかどうかを確認することしかできません。正しくない場合は、前のタイム ステップに戻って dp を調整する必要があります。
ODE15s のような標準の ODE ソルバーはこれを行うことができないか、少なくとも私はそれらを伝える方法を見つけていないようです。このアルゴリズムを実装する可能性はありますか?
助けてくれてありがとう！