- 締切済み
拡散方程式
∂u/∂t=D∇(の2乗)u を解くと u=N/4πDtEXP[-x(の2乗)+y(の2乗)/4Dt]となるようなんですが、なかなか解けません。解き方をどなたか教えて頂けないでしょうか？よろしくお願いします。
- toratora11
- お礼率20% (3/15)
- 物理学
- 回答数1
- ありがとう数3
- みんなの回答 （1）
- 専門家の回答
みんなの回答
- KENZOU
- ベストアンサー率54% (241/444)
拡散方程式は放物型の偏微分方程式と呼ばれていますね。これはGreen関数を使って解くやり方が一般的と思います。ただ、それをここで述べるのは非常に大変なので、解き方がかなり詳しく書かれている次の本をご紹介します。ちょっと複素関数論の知識が必要となるかも知れませんが、テキストにその辺の技術的なことも書かれていますので、一度図書館や本屋等でご覧になられてはいかがでしょうか。 小泉義晴著「現代工学のための量子物理学とグリーン関数」増補改訂版 現代工学社（P112～117）￥2500円 尚、コンピュータで数値計算されるのであれば下記ＵＲＬにプログラム（Fortran言語）が載っています。
関連するQ&A
- 二次元拡散方程式の一般解が求まりません
すみません、拡散方程式で解けない問題がありまして、どなたかご教授ください。 u(x,y,t)の位置(x,y)と時間(t)のみに依存する関数があり、 拡散方程式 du/dt=D*(d^2u/dx^2+d^2u/dy^2) （dは本来は偏微分のパーシャルdです。Dは定数） 一辺の長さが1.0の正方形を考えています。(0<x<1 , 0<y<1) 境界条件は、u(0,y,t)=0.0 , u(x,0,t)=0.0 ,u(1.0,y,t)=0.0 , u(x,1.0,t)=0.0 です。 初期条件は u(x,y,t)=10.0 です。 すみませんができれば解のみではなく方針までお答えいただけると幸いです。よろしくお願いします。
- 締切済み
- 数学・算数
- 二次元拡散方程式の一般解が求まりません
二次元拡散方程式の一般解が求まりません すみません、拡散方程式で解けない問題がありまして、どなたかご教授ください。 u(x,y,t)の位置(x,y)と時間(t)のみに依存する関数があり、 拡散方程式 ∂u/∂t=D*(∂^2u/∂x^2+∂^2u/dy^2) （Dは定数） (0<x<a , 0<y<b) 境界条件は、u(0,y,t)=0.0 , u(x,0,t)=0.0 ,u(a,y,t)=0.0 ,u(x,b,t)=0.0 です。 初期条件は u(x,y,0)=f(x,y) です。 変数分離 u(x,y,t)=X(x)Y(y)T(t) 代入後uで両辺を割る T´/(D*T)＝X´´/X+Y´´/Y 後はD*X´´/X＝α、D*Y´´/Y＝β （α、β、kは定数）ここで，k=-(α+β)とおく。 の３つの微分方程式を解いて初期条件、境界条件を用いて定数を決定します。 X(x)=Acos√αx+Bsin√αx Y(y)=Ccos√βy+Dsin√βy とおいて、境界条件を代入し X(0)=X(a)=0 Y(0)=Y(b)=0 X(a)=Bsin√αa=0 α=(nπ/a)^2 (n=1,2,・・・) Y(b)=Dsin√βb=0 β=(nπ/b)^2 (n=1,2,・・・) 境界条件u(0,y,t)=0.0 , u(x,0,t)=0.0 ,u(a,y,t)=0.0 ,u(x,b,t)=0.0がときのものは 一般解を求められました。 次に， 境界条件u(0,y,t)=0.0 , u(x,0,t)=0.0 ,u(a,y,t)=1.0 ,u(x,b,t)=0.0のときの一般解を求めたいのですが、上手く出来ません。 X(x)=Acos√αx+Bsin√αx Y(y)=Ccos√βy+Dsin√βy とおいて、境界条件を代入し X(0)=0 X(a)=1 Y(0)=Y(b)=0 X(a)=Bsin√αa=1 Y(b)=Dsin√βb=0 β=(nπ/b)^2 (n=1,2,・・・) X(a)=Bsin√αa=1をどう解けばいいのか分かりません。 ご教授お願いします。
- ベストアンサー
- 数学・算数
- 拡散方程式と波動方程式の初期値・境界値問題です。
よろしくお願いします。 １、正方形領域 0≦x≦π、0≦y≦πにおいてu(x,y,t)に関する以下の拡散方程式の初期値・境界値問題を考える。 u_t＝u_xx＋u_yy …(1) u(x,0,t)＝u(x,π,t)＝u(0,y,t)＝u(π,y,t)＝0 …(2) u(x,y,0)＝x(π－x)y(π－y) …(3) (1)m,nが自然数のとき、f_(m,n)(t)sin(mx)sin(ny)が(1)を満足する特解であるとする。（f_(m,n)：添え字） f_(m,n)(t)を求めよ。ただしf_(m,n)(0)＝1とする。 (2)(1)、(2)、(3)を満足する解をu(x,y,t)＝Σ(∞、m,n＝1)a_(m,n)f_(m,n)(t)sin(mx)sin(my)とする。 a_(m,n)を求めよ。 ２、同じ方法を用いて、正方形領域 0≦x≦π、0≦y≦πにおける以下の波動方程式の初期値・境界値問題の解u(x,y,t)を求めよ。 u_tt＝u_xx＋u_yy …(4) u(x,0,t)＝u(x,π,t)＝u(0,y,t)＝u(π,y,t)＝0 …(5) u(x,y,0)＝x(π－x)y(π－y)、u_t(x,y,0)＝0 …(6)
- 締切済み
- 数学・算数
- 拡散方程式の一般解が求まりません
すみません、拡散方程式で解けない問題がありまして、どなたかご教授ください。 u(y,t)の位置(y)と時間(t)のみに依存する関数があり、 拡散方程式 du/dt=D*(d^2u/dy^2) （dは本来は偏微分のパーシャルdです。Dは定数） 境界条件は、 u(±h,t)=Ucosωt （h,ωは定数） となっています。これだけの条件では解けないのでしょうか？？すみませんができれば解のみではなく方針までお答えいただけると幸いです。よろしくお願いします。
- 締切済み
- 数学・算数
- 微分方程式論について
微分方程式論について (1) 関数族{f_n|n=2,3,...} f_n:[0,1]→R f_n=(n^2)x (0≦x≦1/n) f_n=-(n^2)x+2n (1/n≦x≦2/n) f_n=0 (2/n≦x≦1) このとき、{f_n}はいかなる関数にも一様収束しないことをε論法で示せ (2) f:R×R≧0→R f(x,y)=√yと初期値問題 dy/dx=f(x,y) y(0)=0・・・(＊)について (1)f(x,y)はyに関して局所リプシッツ連続ではないことをε論法で示せ (2)定数c≧0に対して、関数 y_c:R→R≧0 y_c(x)=0 (x≦c) y_c(x)=1/4(x-c)^2 (x＞c) は(＊)の解であることを示せ (3) Gronwallの不等式{f,g,uは区間I上の連続関数でg(x)≧0とする} u(x)≦f(x)+∫_(x_0→x)g(t)u(t)dt (x＞x_0) ⇒ u(x)≦f(x)+∫_(x_0→x)g(t)u(t)exp{∫_(t→x)g(s)ds}dt (x＞x_0) を示せ。また、f(x)が非減少関数 ⇒ u(x)≦f(x)exp{∫_(x_0→x)g(t)dt} (x＞x_0) を示せ (4) 初期値問題 dy/dx=y y(0)=a についてPicardの逐次近似法により、解を具体的に構成せよ ただし、y(x)の定義域については考えなくてよい これらの解法を教えてください！ わからなくて困っています…
- 締切済み
- 数学・算数
- 微分方程式について
微分方程式について。 yやdy/dxの形ならば解けるのですが ちょっと変わった形になると解けずに困っております。 回答お願いします。 1 未知関数x(t),y(t)に関する微分方程式 x´(t)=y(t), y´(t)=-x(t)を 初期条件x(0)=a, y(0)=bの下で解け。 2 x=x(t)を変数tのC^∞級関数とする。 このとき、 d^2x/dt^2 +(dx/dt)^2 -4=0 を解け。 3 tの関数x(t)が次の微分方程式を満たすとする x´+x^2+a(t)x+b(t)=0 ただしx´=dx/dtである。 ・x(t)=u´(t)/u(t)のとき、関数u(t)の満たす微分方程式を求めよ。 ・微分方程式 x´=x(1-x)の一般解を求めよ。 長いですが回答お願いします
- ベストアンサー
- 数学・算数
- 支配方程式について
以下の問題を解くことができません。解き方がわかる方、教えてくださると幸いです。 (相図については大まかなイメージを教えてくださるだけでもありがたいです。) ※問題文に不備、不明瞭な点があったら、ご指摘ください。 個体数をy(t)とすると、支配方程式は次のように与えられる。 dy/dt=(a-∫[0~t]y(t)dt)y ここで、累積個体数をx(t)とおくと、 x(t)= ∫[0~t]y(t)dt すなわち、 y(t)=dx/dt したがって、もとの方程式は、 (d/dt)*(dx/dt)=(a-x)*dx/dt=a*(dx/dt)-x(dx/dt) となるが、 (d/dt)*{(x^2)/2}=x(dx/dt) に注意すれば方程式ではtで積分できる。 dx/dt=ax-(x^2)/2+C ※Cは積分定数 初期個数値をy0とおけばx(0)=0なので、上式のC=y0と知れる。よって、 dx/dt=y0+ax-(x^2)/2 すなわち、 y=y0+ax-(x^2)/2 下記の4問が問題です。 (1) dy/dx=0を与えると、そのときのyの値を求めよ。 (2) y=0となるときのxの値を求めよ。ただし、x>0のものを選べ。 (3) 以上の情報を頼りに、横軸にx、縦軸にyをとって相図を描け。 (4)相図を解釈せよ。
- ベストアンサー
- 数学・算数
- 微分方程式、初期値で微分すると？
u=u(x,y), v=v(x,y)は(x,y)∈R＾２のC＾1級関数とし、 （ξ、η）∈R＾２を初期値とする初期値問題 dx/dt=u(x,y), x(0)=ξ dy/dt=v(x,y), y(0)=η の解を x=x(t;ξ,η), y=y(t;ξ,η) とする。 ∂x/∂ξ｜[t=0]＝１、 ∂/∂t(∂x/∂ξ)｜[t=0]=∂u/∂x(ξ,η) を示せ。 って問題なんですけど、初期値（ξ）で偏微分？何のことか全然わかりません。xにt=0代入してからξで微分すれば１にはなりますけどそんなことしたらだめですよね？ 助けてください。
- ベストアンサー
- 数学・算数
