1

これは簡単な質問かもしれませんが、勾配の外側にある空間的に変化する係数で PDE を解きたい場合、使用する正しい FiPy 構文は何ですか? これまでに見たすべての例は、勾配内の係数についてのみ説明しています。

例えば:

d/dt(シグマ) = (1/r) d/dr (r^0.5 d/dr(nu シグマ r^0.5))

(数値要因は無視しています)

Sigma(t,r) について解きたいと思います。d/dr の前で (1/r) を処理するにはどうすればよいですか?

この単純な方程式はマッサージできるので、勾配の外側にある空間的に変化する係数について心配する必要はありません (または単に係数を時間導関数項内に移動するだけです)。私が解決しようとしている実際の問題であり、そのトリックはもはや有効ではありません。たとえば、方程式が次のようになっている場合はどうすればよいですか。

d/dt (変数) = f(r) d^2/dr^2 (変数) + g(r) d/dr (変数)

どんな助けでも大歓迎です!

4

1 に答える 1