1次熱伝導方程式を差分法で解くという問題をC言語で書こうと言う問題を友達からまる投げされたのですが、解き方が全く分かりませんん;;
というのも差分法などの本を見てもよくわからない記号ばかりで。。。。
1次熱伝導方程式においで(∂u/∂t = κ(∂^2u/∂x^2))
初期条件として
u(x,0)=30
K=0.08
?tは指定されていませんでした(解くに何でも可っぽいです)
境界条件として
ux(L,t)=60℃ Lは1
u(0,t)=700℃
でも文章をよんでると有限の場合
棒の両端の温度が0に保たれているとき
u(0,t)=u(L,t)=0
棒の両端が断熱されているとき
ux(0,t)=ux(L,t)=0
今回の境界条件の場合はどうといていったらいいのか迷子になっているのですがどうかご教授してもらえたらありがたいです。
No.1ベストアンサー
- 回答日時:
1)領域を確認します。
0≦x≦L
0≦t≦T(Tは時間ステップ×刻み数ですが100ステップぐらいを意識しておきます。)
2)図を描きます。
x方向を横軸に、t方向を縦軸にとります。
tを10に刻んで下から上に(時間軸)i=1,2,...11
xを10に刻んで左から右にj=1,2,.....11に番号を振ります。
t-x平面を格子で覆います。格子点(節点)はu(i,j)(i=1,11:j=1,11)で表します。
刻み幅をΔt、Δxで表します。
3)境界条件、初期条件を与えられている、x-t平面の境界を赤で塗り、境界値、初期値を記入します。
直ちに気付くことは(1,1)で境界条件と初期条件が整合していないことです。この不整合はすぐに累積して解は発散します。よってこのままでは解けません。
ここではこの不整合は解消されたと仮定して話を進めます。
4)1次熱伝導方程式
(∂u/∂t = κ(∂^2u/∂x^2)) (1)
を差分表示します。一回微分の時間軸は2点差分、2回微分のx軸は3点差分で表します。その結果
(u(i+1,j)-u(i,j))/Δt=κ(u(i,j-1)-2u(i,j)+u(i,j+1))/(Δx)^2 (2)
5)格子の真ん中付近で(2)で表される4点を太く書いて、じっくりと眺めます。
今やろうとしているのはコンピューターによる計算です。計算は式の右辺の各項が既知でそれらを加減乗除して計算結果を今まで未知だった左辺の値として入れてやります。つまり4点 u(i+1,j)、u(i,j)、u(i,j-1)、u(i,j+1)のどれが既知でどれが未知かということです。この4点のセットを動かしながら考えると結局
u(i+1,j)が未知、u(i,j)、u(i,j-1)、u(i,j+1)が既知になることが解ります。従って
u(i+1,j)=u(i,j)+(κΔt/(Δx)^2)(u(i,j-1)-2u(i,j)+u(i,j+1)) (3)
が計算式です。
t=0、x=0から右方向にX=L(j=11)まで動かし、次にt=Δt(i=2)の線上を左から右に動かしてu(i+1,j)の値を次々に求めていきます。
プログラムは2重のループの計算になるでしょう。
何回も試行錯誤するkとになると思います。
うまく行ったらメッシュを細かくして精度を上げてください。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 移流熱拡散方程式の解き方 フーリエ変換 1 2022/08/15 15:25
- 物理学 ニュートンの冷却法則と熱伝導方程式について 3 2023/03/05 19:51
- 工学 電磁気学の質問です。 電流による電磁束密度ベクトルポテンシャルの計算ですが、ベクトルポテンシャルのポ 1 2022/04/19 17:23
- 物理学 熱力学 エントロピー 断熱自由膨張 熱力学第2法則 クラウジウスの不等式 2 2022/07/14 12:58
- 物理学 熱伝導方程式について 1 2023/03/22 20:55
- 大学受験 助けてください。 大学入試まであと38日なのですが英語が2割です…。 第1志望は5割位なのですが、2 8 2023/01/03 19:37
- 数学 二次関数 2次方程式ax^2+bx+c=0を分解する際 a>0という条件を付ける と参考書に書いてあ 0 2023/02/08 06:45
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- 数学 数学微分方程式の問題です。次に書く問題を教えて欲しいです。上端を固定された長さlの棒の先に質量mの質 2 2022/04/29 21:27
- 高校 対数方程式につきまして 4 2022/05/05 07:55
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
物理の有効数字2桁
-
bcc(面心立方体)の四面体位置で...
-
保存力→ポテンシャル→力 の計...
-
高圧ガス製造保安責任者 甲種...
-
有限積分法について
-
pHについてです
-
1次熱伝導方程式を差分法で解...
-
算数の問題がわかりません。
-
高さHのビルの屋上から初速v₀で...
-
∫f(x)dx と∫dxf(x)の違いは何で...
-
1/(2×3.14) × √(0.5×9.8/0.4) ≒...
-
斜面を上る運動 摩擦なし
-
局所ヌセルト数から平均ヌセル...
-
写真の問題の解き方を調べたと...
-
Lagrangeの運動方程式→Lorentz...
-
FEMのモード法と直接法の違...
-
気体分子の確率分布
-
小学生五年生 分数の足し算引き...
-
この計算この答えであってますか?
-
電磁気学の面積分
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
bcc(面心立方体)の四面体位置で...
-
物理の有効数字2桁
-
二次元流れにおいて、流体のx、...
-
てこの計算問題 輪軸の応用問題
-
熱力学の問題で、(6)について、...
-
pHについてです
-
局所ヌセルト数から平均ヌセル...
-
保存力→ポテンシャル→力 の計...
-
整数で答えよ、のように書かれ...
-
高校の理科の答えの書き方
-
水素原子の波動関数の直交性に...
-
高さHのビルの屋上から初速v₀で...
-
大学物理の電位差計式抵抗
-
算数の問題がわかりません。
-
単振動の運動方程式の解
-
たたみ込み積分って名前つける...
-
変数分離が成功したからといっ...
-
斜面を上る運動 摩擦なし
-
FEMのモード法と直接法の違...
-
1リットルは1辺が1cmの立方体が...
おすすめ情報