- 締切済み
拡散方程式の一般解が求まりません
piro2dogの回答
- piro2dog
- ベストアンサー率33% (1/3)
変数分離法を適用してみてはいかがでしょうか。 ※変数分離法を適用するための条件についてはご自身で調査して下さい。 変数分離法とは、u(y,t)が yのみに依存する関数Y(y)と、 tのみに依存する関数T(t)との積で表せると仮定して解く方法です。 この場合、拡散方程式に、u(y,t)=Y(y)T(t)を代入すると、 Y(y)dT/dt=D(T(t)d^2Y/dy^2) ※dは常微分のd となります。 これを変形すると T'/T=D*(Y''/Y) ※U'はUのtでの微分、Y''はYのyでの二階微分 となり、左辺がtのみの関数、右辺がyのみの関数となり、 結局、左辺も右辺も定数であると分かります。 なので、T'/T=D*(Y''/Y)=C(定数)とおくと、 T'=CT Y''=(C/D)Y と単純な二つの常微分方程式が得られます。 あとはこれを解くだけです。
noname#95312
関連する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をどう解けばいいのか分かりません。 ご教授お願いします。
- ベストアンサー
- 数学・算数
- 方程式の解についてです
U*(ξ)=∫(-∞～ξ)U(τ)dτ (ξ∈R） これを微分して dU*(ξ)/dξ=U(ξ) d^2 U*(ξ) / dξ^2=dU(ξ)/dξ これらを以下の方程式☆に代入します d^2 U(ξ) /dξ^2 + ξ/2•dU(ξ)/dξ－αU(ξ)=0 --(☆) （α:任意定数) そうするとdU(ξ)/dξ+ξ/2•U-(α+1/2)U*(ξ) ----(★) となりますよね。 ★をξで微分すると d^2 U(ξ) /dξ^2 + ξ/2•dU(ξ)/dξ+1/2U(ξ)-(α+1/2)U(ξ) これは☆より＝０となります。 ここで、 U(-∞)=(dU/dξ)(-∞)=0ならば★=0 ということなのですが、これはどういうことでしょうか？ U(-∞)=(dU/dξ)(-∞)=0という条件は何故必要なのでしょうか？ どなたか解説よろしくお願い致します。
- 締切済み
- 数学・算数
- 微分方程式の解法を教えてください！
常微分方程式の解法はどんなものがあり、どのような場合に適用すれば解けるでしょうか。 解法を覚えても、それが適用される場合についての判断ができません。教えてください！ 以下の場合だとどのように解けばよいでしょうか。 （１）d^2x/dt^2＋ω^2x=0の一般解の求め方。（ωは定数） （２）dx/dt=-c^2y、 dy/dt=c^2x の一般解の求め方。（cは定数） （３）dx/dt=u、 du/dt=-kx-cu＋f(t) (k,cは定数) のとき (1)f(t)=0のとき、t=0でx=x0のもとでの解を求め る。 (2)f(t)=cosωtのときの解。
- 締切済み
- 数学・算数
- 偏微分方程式の問題です。準線形方程式 u・(∂u/
偏微分方程式の問題です。準線形方程式 u・(∂u/∂x) + (∂u/∂y) = 1 の解で、初期曲線がx0(s)=s、y0(s)=2s、u0(s)=-s (0≦s≦1)で与えられるものを求めよ。 自分で解いてみたら途中までこんな感じになりました。 >>dx/dt = u、dy/dt = 1、du/dt = 1を解くと x=ut+s、y=t+2s、u=t-s おそらくこの連立方程式を解くのかと思うのですがうまく行きません。どうすれば良いのでしょうか。回答お願いいたします！
- ベストアンサー
- 数学・算数
- 微分方程式の一般解
単刀直入に質問いたします。 du/dt=(1-u)u 簡単なようですが上の式の一般解が分かりません。 変数分離形だと思い解いてみたのですが、変数を分離して両辺積分してまとめると、 log|u/(1-u)|+C=t となったのですが、むりやりu=の形に直すとあまりいい形になりません。 この問題の後に、 (2)u(0)=aのもとで上記方程式を解け。但しaを実定数とする。 (3)(2)で求めた解がt≧0全体で定義されるためのaの満たすべき条件をもとめよ。 という問題が続くので、一般解をできるだけきれいな形で出しておきたいと思い質問いたしました。 ちなみに自分で無理やり解いた解は、 u=(±exp(t)+C)/(±exp(t)+C+1) (Cは積分定数) です。 間違っている、これでよい等どんなことでも良いので助言をお願いします。
- ベストアンサー
- 数学・算数
- 連立微分方程式と特殊解について
dx/dt=-3x-y, dy/dt=4x+2yの特殊解が定数A、B,mを用いて、x=Aexp(mt), y=Bexp(mt)と表されるとして、微分方程式の一般解を求める方法を教えてください。
- ベストアンサー
- 数学・算数
- 偏微分方程式の問題です。
aを定数とするとき、次の偏微分方程式を解け。 du/dt+a・du/dt＝0 ただし、初期条件を以下とする。 u(x,0)＝bx^-x （b：定数） （ヒント） u(x,t)=g(x)h(t)と変数分離できることを仮定してよい。 解答・解説お願いします。 考え方の提示のみではなく、答えまでよろしくお願いいたします。
- ベストアンサー
- 数学・算数
- 反応拡散方程式の定常解について
反応拡散方程式 dU/dt = f(U,V)+ΔU dV/dt = g(U,V)+ΔV (fとgは非線形な関数でΔは空間1次元のラプラシアン) には、条件が整えば空間的に非一様な定常解（時間変化しない解） が存在します。 このような空間非一様性を満たす定常解 0 = f(U,V)+ΔU 0 = g(U,V)+ΔV をノイマン境界条件のもとで数値的に求めようとして、 ラプラシアン(Δ)を離散化した次の連立方程式 0 = f( U<i>, V<i> ) + [U<i+1>-2U<i>+U<i-1>]/h^2 0 = g( U<i>, V<i> ) + [V<i+1>-2V<i>+V<i-1>]/h^2 (i=0～N) をニュートン法で解いてみました。 結果、 初期値(u(x)とv(x)の形)によってはまぐれで空間非一様な定常解らしき状態に収束することがあるのですが、ほとんどの場合発散してしまうか、一様な解（u(x)=0,v(x)=0）に収束してしまいます。1000回ほどランダムに初期値を変えて試行したところ空間非一様＜らしき＞解に収束したのはたった数回でした。また、離散格子点を多くとった場合（例えばN=200）には、まったく収束しなくなりました。直感では、初期値さえうまくとれば収束しそうな気がしたのですが、、 質問： 上の反応拡散方程式のような非線形連立偏微分方程式の定常解を、数値的に求める手段というのはあるのでしょうか？ アドバイスよろしくお願いします。
- ベストアンサー
- 数学・算数
- 連立微分方程式の解き方
d^2x/dt^2=f(t) d^2y/dt^2=g(t, y, dy/dt, z, dz/dt, d^2z/dt^2) d^2z/dt^2=h(t,d^2y/dt^2, z, dz/dt) というような連立微分方程式を解きたいのですが，どのような方法で解くことができるのでしょうか？ ルンゲ・クッタは適用できますか？ （d^2y/dt^2の式にd^2z/dt^2が変数として出てきたりしているので，わからなくなってしまいました．） 宜しくお願いします．
- ベストアンサー
- 数学・算数
お礼
変数分離法はこちらもやってみたのですが、 その場合 T=exp(Ct) と表現されてしまうため、境界条件とうまく適合されません。 Cを虚数としてRe(T)（Tの実部分）としてみたところそれ以降どう扱えばいいかわからず頓挫してしまいました。 貴重なご意見ありがとうございました。