- ベストアンサー
移流拡散方程式の微分が含まれていない解を教えて下さい
- 移流拡散方程式の微分が含まれていない解を教えて下さい。
- 土壌中の化学物質量が、時間とともに土壌深度に拡散、移流するのを考えています。エクセルで変化を観察したいと思っています。
- 移流拡散方程式を計算する方法や時間の設定について教えていただけないでしょうか。
- みんなの回答 （3）
- 専門家の回答
質問者が選んだベストアンサー
- ベストアンサー
素人のコメントです。 （１）地表で x = x0 とすると、x < x0 には土が存在しないにも関わらず、x = 0 から上下に対称的に拡散する場合の解 p(x,t)=(1/√4πDt)*EXP(-x^2/4Dt) （＊） を使ってよいのだろうか、という疑問をもちました。 （２）深さ方向の拡散だけでよいのだろうか。水平方向への拡散は無視できる状況なのだろうかとも思いました。 （３）もし（＊）式が適用できるのであれば、移流の効果については http://thesis.ceri.go.jp/center/doc/geppou/ceri/0005006080.pdf が参考になるかもしれません。その（１５）式を見ると、流速を u として p(x,t)=(1/√4πDt)*EXP{-(x-ut)^2/4Dt} となるように思えます。（ただし、いまの場合の流速の意味がよくわかりません。） （４） >tはどうやって設定すればよいのでしょうか。 （＊）式にある x^2/(Dt) は無次元の量だと思いますから、x^2/D と同じ単位を使われたらよいのではないでしょうか。
その他の回答 (2)
>この式も移流によってピークがずれますか？ erfc((x-vt)/2√Dt) の部分はそうだと思いますが、その後ろの項については私にはわかりません。
お礼
お返事ありがとうございました。 x-vtが入ってるので、ピークがすれるような気はしています。 もういちどエクセルで頑張ってみようと思います。 本当に助かりました。ありがとうございました。
>p(x,t)=(1/√4πDt)*EXP{-(x-ut)^2/4Dt}をエクセルで計算しましたが、やはりピークがでる深度が変化しなかったです。 計算がどこかで間違っていると思います。上の式は、ピークの位置が x = u t に従って移動することを示しています。
お礼
すみません、できました！xの間隔を変えたらできました。お騒がせしました。。 ちなみに p(x,t)=Po*1/2*((erfc((x-vt)/2√Dt))+exp(vx/D)*(erfc((x+vt)/2√Dt))) という解を http://www.ekouhou.net/%E5%9C%B0%E7%9B%A4%E6%94%B9%E8%89%AF%E5%B7%A5%E6%B3%95/disp-A,2008-267016.html で、見つけたのですが、（土壌中のイオン濃度の変化の式なのでおそらくこちらを使うべきなのかなと） この式も移流によってピークがずれますか？ エクセルでやってみたのですか、ピークが下方にずれていかないのです。。。 私の計算ミスでしょうか。
お礼
お返事ありがとうございました。 （１）についてですが、確かにその通りです。。対照的な拡散の式を使っていたからピークが移動しなかたことに気付きました。 p(x,t)=(1/√4πDt)*EXP{-(x-ut)^2/4Dt}をエクセルで計算しましたが、やはりピークがでる深度が変化しなかったです。 もう一度土壌中の移流拡散方程式を探してみます。 重要なことに気付かせていただき、本当にありがとうございました。