- ベストアンサー
球座標系の拡散方程式と成長モデルについて
- 球座標系の拡散方程式とは、球状の気泡の成長を表すモデルです。
- 式(1)は、半径がRのときの濃度cの半径方向の変化を表しています。
- 式(2)は、球座標の拡散方程式の径方向成分であり、積分法で解析的に求めることが可能です。
- みんなの回答 （3）
- 専門家の回答
質問者が選んだベストアンサー
こんにちは。 やはり両辺を体積分したのではないですか？ > d/dt((4pi*r^3*P)/(3RT))=4piR^2D*∂c/∂r |_r=R ・・・(1) この右辺は、まさしく(2)の右辺を体積分したものになっています。 ∫r^2dr∫dΩ (右辺) を実行すると、∫dΩ = 4π なので、すぐに導かれます。∫drの上限は R ですから、D r^2*∂c/∂r に R を代入したものになりますね。 これはガウスの積分定理を適用したのと同じことです。 左辺のほうは、∫dV ∂c/∂t = ∂/∂t (∫dV c(r,t) ) のようにして体積分を実行できます。 ところで、(2) で、(∂/∂r)* のように「*」で掛け算になっているのはおかしいですよ。微分演算子は数ではないので、掛け算ではありません。
その他の回答 (2)
- aquarius_hiro
- ベストアンサー率53% (194/360)
こんにちは。 そうですね… とりあえず、わかることだけ。 (2) の右辺は、c(r,t) について、div(grad c) を極座標で書いたものであり、 一方、(1) は、grad c のr方向の成分の r=R での値なので、 おそらく、(2) か、そこから出てくる式を、球状の体積を V として、体積分したのではないですか？ ガウスの積分定理 ∫dV div( grad c ) = ∫dS [grad c]_r を使います。右辺の [grad c]_r は、ベクトル grad c の r方向の成分を表します。それで、球対称な問題では、球面上で [grad c]_r は一定なので、右辺は S [[grad c]_r]_{r=R} になり、(1) に比例します。 > ・(2)式の右辺は微分してD(∂^2c/∂r^2 + (2/r)*(∂c/∂r))と変形できるとおもいますが、積分するのだからこんなことしなくてもいいのでしょうか。 体積分すれば ∫r^2 dr が出てきますので、右辺はそのままで積分できますよね。ちょうど1/r^2の因子と、r^2 dr の r^2 が打消しあいますので、微分をばらさないほうが良いと思います。 また、この体積分で dr が実行できますが、それは積分定理を使ったのと同じことです。 > ・定常状態として(2)式の左辺を0として、適当に境界・初期値条件を入れて解いたということでしょうか。それとも非定常状態としてそのまま左辺も右辺も積分するのでしょうか（それができるのかどうかもわかりませんが）。 境界条件を使っていることは間違いないですが、定常状態かどうかは、情報が少なすぎてわかりません。
補足
ありがとうございます。１．の方の補足に(1)式を記述しました。上記のアドバイスをフォローするのにちょっと時間がかかりますので、何かその他の追加があればよろしくお願いします。
- rabbit_cat
- ベストアンサー率40% (829/2062)
＞∂c/∂r |_r=R ・・・(1) ＞という項が出てきます。 どう出てくるのですか？ 少なくとも一つ目の式というもの自体を示してくれないと答えようがないです。
お礼
追加です。(1),(2)式のほかに気泡界面における力のバランス式がでてきます。 dr/dt = r/4y(P-2G/r)・・・(3) y:粘性係数,G:気泡の表面エネルギー
補足
すみませんでした。 d/dt((4pi*r^3*P)/(3RT))=4piR^2D*∂c/∂r |_r=R ・・・(1) となります。左辺は気泡中のガス粒子数（気体の状態方程式より）の時間変化です。
お礼
ありがとうございました。おかげさまで導出ができました。 ＞ところで、(2) で、(∂/∂r)* のように「*」で掛け算になっているのはおかしいですよ。微分演算子は数ではないので、掛け算ではありません 知りませんでした。以後気をつけます。