【有限体積法】1次元定常移流拡散方程式を風上差分で解く C++コード付き

今回は1次元定常移流拡散方程式を有限体積法(風上差分)で解いていきます。C++コード付きです。1次元定常移流拡散方程式とは のような二階の常微分方程式のことをいいます。ここで、 は溶質の濃度、 は風などによる移流速度、 は拡散の効果を表す拡散係数と解釈することができます。ある溶質が移流や拡散のもとでどのよう…