【数理モデル】線型反応移流拡散方程式とマルサスモデル(Malthusモデル)
前回1次元Fisher-KPP方程式(フィッシャーKPP方程式)を有限要素法(Galerkin法)で解く話を書いたのですが、そのときに線型反応移流拡散方程式とマルサスモデル(Malthusモデル)の関係について気がつきました。せっかくなので文章として残しておきます。
まず、1次元の線型反応移流拡散方程式は
のような偏微分方程式です。ここで、 は物質の濃度、 は流速、 は物質の広がり具合(拡散)を表す拡散係数、 が反応により減少する度合いを表す係数と解釈することができます。ここで です。後で説明しますが、偏微分方程式が一本しかないとき、 だと解が発散してしまいます。システムの場合(偏微分方程式が何本かある場合)は必ずしも正でなくてもよかったと思います。初期条件は で与えます。 は時間方向と空間方向ともに変化するので二変数関数 となっています。ある物質が移流、拡散、反応のもとどのような分布になるかを表しています。
さて、Fisher-KPP方程式は拡散方程式とLogisticモデルを組み合わせたものと考えることができるのでした。
これと同様にして、線型反応移流拡散方程式は移流拡散方程式とMalthusモデルの組み合わせたものと考えることができます。まず反応項を右辺へ移項します。
これは移流拡散方程式
とMalthusモデル
を組み合わせたものです。すなわち移流拡散によって物質の空間的な広がりを、Malthusモデルによって(その場所における、減少するような)物資の反応をあらわしているのです。もし だとすると、 となり、指数関数的に増大してしまい解が発散してしまいます。なので条件として がついているのです。移流拡散方程式とMalthusモデルの組み合わせと考えるとすっきり理解できます。
参考
- 作者: 二宮広和,三村昌泰,竹内康博,森田善久
- 出版社/メーカー: 共立出版
- 発売日: 2014/07/09
- メディア: 単行本
- この商品を含むブログ (2件) を見る
- 作者: 柳田英二
- 出版社/メーカー: 東京大学出版会
- 発売日: 2015/06/26
- メディア: 単行本
- この商品を含むブログを見る