![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?e8efa67)
No.1ベストアンサー
- 回答日時:
変数分離型の解から出発します。
c=R(r)T(t)とおき
∂c/∂t=D/r^2*∂/∂r(r^2*∂c/∂r)
に代入する。以下拡散係数Dは定数とする。
RdT/dt=(DT/r^2)*d/dr(r^2*dR/dr)
(1/DT)dT/dt=(1/R)(1/r^2)d/dr(r^2*dR/dr)
変数分離ができたのでこの値をk(t,rによらない定数)とおくと
(1/DT)dT/dt=k (1)
(1/R)(1/r^2)d/dr(r^2*dR/dr)=k (2)
(1)より
dT/dt=kDT
解は
T=ce^(kDt) (3)
(2)より
d^2R/dr^2+(2/r)dR/dr-kR=0
k=1/p^2,r=puとおくと
d^2R/du^2+(2/u)dR/du-R=0
この解は変形球ベッセル関数(n=0)(url参照)
参考URL:http://ja.wikipedia.org/wiki/%E3%83%99%E3%83%83% …
No.2
- 回答日時:
#1です。
いくつかurlを見てみると私の用いたk<0として解いている例が多いようです。
拡散して希釈していくことを前提にしているのでしょうか。
ソースがあって一定の放出をしている場合は拡散によって閉空間の濃度は大きくなる場合もあるように思います。
ともあれ、球座標系におけるノイマン型境界地問題に触れたペーパーがありました。
参考URL:http://www.gfd-dennou.org/arch/spmodel/1d-radial …
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
Debug.Printで表示される内容を...
-
電流の時間微分、電圧の時間微分
-
物理の計算で m×dv/dt×v=d/dt{...
-
質問です。 磁束密度Bの一様な...
-
d^2r/dt^2の意味
-
ソレノイド(東大レベル)
-
v^2-v0^2=2ax 今日この式を習っ...
-
運動方程式の微分積分の計算
-
物理の問題です! ∫[mr{d^2θ...
-
物理の問題
-
微分積分のdの意味
-
d/dx=dt/dx * d/dt =d/dt * dt/...
-
摩擦を考慮したサイクロイド曲線
-
EXCEL上の数字を自動で振り分け...
-
ナビエストークスについてです。
-
蒸発速度
-
dx/dt=√(1-x^2)の一般解の求め...
-
直流磁界Hと直流電流密度J(A/m^...
-
単振動
-
物理 気球の問題
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
雨滴の運動質量が変化する落体...
-
d^2r/dt^2の意味
-
電流の時間微分、電圧の時間微分
-
質量流量の記号「・ の読み方を...
-
力学の雨滴の落下の問題です
-
ポテンシャルエネルギーから力...
-
物理の計算で m×dv/dt×v=d/dt{...
-
機械力学の問題です!!!
-
Debug.Printで表示される内容を...
-
最後のdv/dtは何でしょうか。
-
運動方程式を求めてください
-
微分積分のdの意味
-
運動方程式の微分積分の計算
-
EXCEL上の数字を自動で振り分け...
-
v^2-v0^2=2ax 今日この式を習っ...
-
d/dx=dt/dx * d/dt =d/dt * dt/...
-
微分記号“d”について
-
dx/dt=√(1-x^2)の一般解の求め...
-
加速度 a=dv/dt = (d^2 x) /dt^2
-
ラプラス方程式
おすすめ情報