- ベストアンサー
拡散方程式
こんにちは。1次元の拡散方程式で、時間tを求めたいのですが、解き方が分かりません。どなたか教えていただけませんでしょうか。ちなみに初期条件はt=0で0<x<∞、C=0、all for tのときx=0、C=C0(シーゼロ)らしいです。
- yatomiryou
- お礼率85% (6/7)
- 物理学
- 回答数1
- ありがとう数1
- みんなの回答 (1)
- 専門家の回答
質問者が選んだベストアンサー
大抵の偏微分方程式のテキストに載っていると思いますが、さしずめ このサイト↓を読んで勉強してください。 http://www2.kobe-u.ac.jp/~iwayama/teach/kisoIII/2005/chap5.pdf
関連するQ&A
- 拡散方程式を解いてください!!!
拡散方程式を解いてください!!! 一次元の半無限体への拡散問題です。 初期の状態は濃度がC(x=0, t=0)=C*(定数)とC(x, t=0)=0(x>0)で、プラス方向に拡散します。 境界条件はlim x→∞ C(x, t)=0です。 お願いします。
- 締切済み
- 科学
- 拡散方程式のプログラミング
拡散方程式を種々の初期条件、境界条件の下にプログラミング化(できればfortranがいいですが、その他でも大歓迎です) まず、前回の質問たくさんの方に返事していただきました。 本当にありがとうございます!!!前回の質問を整理しました。 厚さがq(μm)、底面が気体が透過しないように設定されたペイントの表面の酸素濃度C(t)を時間幅0<=t<=t0(最大時間t0まで)周期的に変えます。 ここで、拡散方程式 σC/σt=Dm*(σC^2)/(σZ^2) C:気体濃度、Z:厚さ方向、t:時間、Dm(μm2/s):拡散係数 というのがあり、これを中心差分近似したものが、 (C(i,t+1)-C(i,t))/Δt=(C(i+1,t)-2C(i,t)+C(i-1,t))/(ΔZ^2) となります。この拡散方程式に加えて、初期条件、境界条件を以下のように加えます。 ・初期条件 (初期状態は厚さ方向に一様の真空状態) →t=0;C=0(この式は正しいかわかりません↓↓) ・境界条件 ペイントの底面では厚さ方向にて濃度が同じ。 →Z=h;σC/σZ=0 この条件を離散化したものが、 →(C(1,t)-C(0,t))/ΔZとなります。 また、境界条件といたしまして、ペイント表面の濃度変化を正弦関数C(t)=Ksin2πft (K:定数、f:周波数、t:時間) として与えます。 ここで、補足をいたしますと、 ペイント表面の気体濃度を変化させると、その表面の気体はペイントの拡散係数に応じてペイント内部に拡散します(ペイント底面は気体は透過しない)。 ここで求めたい問題は、ペイント内のある時間での気体の厚さ方向での濃度分布、かつ、その厚さ方向の濃度を合計したペイント全体の気体濃度、さらに、一番求めたいのはペイント全体の気体濃度の時間変化(0<=t<=t0まで)プログラミング化するのですが。。。 様々な方のアドバイスを聞きたくてたくさん質問しています。。。申し訳ありません。
- 締切済み
- 物理学
- 拡散方程式(フィックの第2法則)
次の問題の解き方がわかりません。 どなたか手助けをいただければ大変ありがたいです。どうかよろしくお願いします。 「問 拡散方程式(フィックの第2法則)を記せ。また、時刻t=0において、水が入った直方体の容器の1つの面上に均一の濃度で存在していた色素分子が時間とともに内部に拡散していく過程について、面からの距離xの時刻tにおける色素の濃度c(x,t)を表す式を記せ。(指数部分のみでよい) また、時間Tが経過した後の、色素分子の平均拡散距離はどの程度になるか。」 どうか知恵をお貸しください。 拡散方程式は∂c/∂t=D(∂^2c/∂x^2)だと思うのですが・・・
- 締切済み
- 化学
- 拡散方程式の数値計算(有限差分法)による解き方
下記の拡散方程式を数値計算(有限差分法)で解き方について教えてください。 拡散方程式 D∂^2C(x,y,t)/∂x^2+u∂C(x,y,t)/∂x=∂C(x,y,t)/∂t 境界条件 D∂C(x,y,t)/∂x=(K/β-u)×C(x,y,t)-KC'(y,t) at x=0 初期条件 C(x,y,0)=Co at t=0 質量保存則 ∂C'(y,t)/∂y=4/vd(K/β×C(0,y,t)-KC'(y,t) ) ---------------------------------------------------------- また、可能であれば、有限差分法以外にも数値解法(有限要素法、有限体積法など)があると思いますが、拡散方程式を解く場合、どの解法が一般的に適しているのか教えてください! よろしくお願い申し上げます。
- 締切済み
- 数学・算数
- ノイマンの制限条件における拡散方程式の解について
以下の境界条件, ∂c(-d,t)/∂t=∂c(d,t)/∂t=0(-d≦x≦d) 初期条件 c(0,0)=1,c(x,0)=0(x≠0) における拡散方程式の解は c(x,t)=d/2+1/d*Σexp(-(Dπ^2*m^2*t)/d^2)cos(mπx/d) で宜しいでしょうか. お分かりの方計算過程も教えていただけたら幸いです.
- 締切済み
- 物理学
- 二次元拡散方程式の一般解が求まりません
すみません、拡散方程式で解けない問題がありまして、どなたかご教授ください。 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 です。 すみませんができれば解のみではなく方針までお答えいただけると幸いです。よろしくお願いします。
- 締切済み
- 数学・算数
お礼
connykellyさん こんにちはyatomiryouです。拡散方程式を教えていただきありがとうございました。是非参考にさせていただきます!