A 回答 (5件)
- 最新から表示
- 回答順に表示
No.5
- 回答日時:
断面積A(x).、熱伝導率κ(x)、温度T(x)、発熱密度H(x)とすると、解くべき微分方程式は
d/dx [A(x)κ(x) dT(x)/dx]=-A(x)H(x)
各部材の中ではA,κ,Hは定数なので、
i番目の部材の解は2次の係数が-H_i/2κ_iであるような2次多項式になります。
1次の係数と定数項(格部材ごとに定数を決めるので全部で6個)は境界条件から決める必要があり、
①アルミの端でj(x)=0である
②部材の接合部でT(x)が連続である(接合部2箇所)
③部材の接合部でj(x)が連続である(接合部2箇所)
④ステンレスの先端でT(x)=T_0
から決める事になります。
ここで位置xでの熱流をj(x)=A(x)κ(x)dT(x)/d(x)としています。
未知数も式の場所も6つなので、計算を頑張れば答えは出てくるはずです。
1次元熱伝導方程式による近似(太さ方向の温度勾配を無視する近似)が適切かはお書きの内容からは判断できないので、計算の目的に応じてご自身で判断してください。(それで良いという判断の元での質問には見えますが)
No.4
- 回答日時:
いくつかの部材が直列しているという設定ですね。
もし、どの部材も長さに比べて太さがうんと細いのであれば、温度分布は(ご質問の通り)1次元で近似的に表せて、それなら手計算でやれるだろう。すなわちL[0] = 0
L[k-1]<L[k] (k=1,2,...,N)
とする。N個の一様な部材が直列してできている物体で、
L[n-1]≦x<L[n]
の範囲では材質も断面も一様である。
n番目(n=1,2,...,N)の部材の単位長さあたりの発熱量をq[n]、単位長さあたりの比熱をR[n]、熱伝導率をD[n]とする(「単位長さあたりの」というのは、単に材質に依存するだけでなく、「断面積に比例する」という意味であることに注意)。
まずは全体の温度分布が平衡状態になったときにどういう解がありうるかを考えてみる。
xにおける単位長さあたりの熱量をQ(x)とすると、L[n-1]≦x<L[n]のとき、(xによる微分を ' で表すことにして)
Q''(x) = - D[n]q[n]
より
Q(x) = (-D[n]q[n]/2)x² + C[n,1]x + C[n,2]
ただしC[n,1], C[n,2]は定数。
また、xにおける温度をT(x)とすると、
T(x) = Q(x)/R[n]
だから、Tは「区分的に連続な2次関数」で表され、未知数は
C[1,1], C[2,1], ..., C[N,1]
C[1,2], C[2,2], ..., C[N,2]
の2N個である。
ただし、x<0へは熱が移動しないので
Q'(0) = 0
より
C[1,2] = 0
また x=L[N]では温度が一定値T₀なので
T(L[N]) = T₀
という境界条件と、部材の接続部で温度が一致するので
((-D[n]q[n]/2)(L[n-1])² + C[n,1]L[n-1] + C[n,2])/R[n] = ((-D[n]q[n+1]/2)(L[n])² + C[n,1]L[n] + C[n,2])/R[n]
であるという制約がつく。これらの制約条件によって未知数はN-1個に減る。
言い換えれば、N-1次元の自由度が残されているわけで、この自由度の範囲内で無限個の平衡解がある。
というわけで、平衡状態を考えるだけでは話が確定しない。例えば「全部の部材が温度T₀である状態から出発する」というような初期条件を与えて、平衡に至るまでの過渡的変化を考えることになるだろう(けれども、それはご質問には含まれていない)。もちろん、方程式に∂Q/∂tの項を付け加えるだけのことであり、初期条件が簡単なら手計算でも扱えるに違いない。
No.3
- 回答日時:
#1です。
訂正「特に→解くに」。ところで、「周辺は断熱条件」と書いていますがステンレスの端面以外は外部との熱の受授が無いという意味ですか?それなら熱流にとって太さ(高さと幅)は無関係です。三種類の材料の太さがあまり違わない(式で示す)という条件下なら一次元の熱伝導方程式で近似できます。太さが違うと接触面近傍は1次元では表現できませんまた、ステンレスの端面からしか熱が出ないとすると、大体の条件ではどんどん温度が上昇し定常状態にはならないでしょう。条件も考え直したらいかがですか。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 急募! 熱力学の問題です 物質の定熱容量Cpは(∂H/∂T)P 1. 偏微分の計算規則のみ使い、定熱 1 2021/12/03 19:02
- 物理学 自由電子モデルの質問 1 2021/10/21 09:28
- 工学 電磁気学の質問です。 電流による電磁束密度ベクトルポテンシャルの計算ですが、ベクトルポテンシャルのポ 1 2022/04/19 17:23
- 物理学 ニュートンの冷却法則と熱伝導方程式について 3 2023/03/05 19:51
- 工学 制御工学についてです。 1巡伝達関数Lが L=k/(s+1)(s+2)(s+3) である。kをゲイン 2 2023/01/31 09:28
- Excel(エクセル) キーボード入力が変になった?! 2 2021/12/28 09:59
- 物理学 マクスウェル方程式についての超基本的な質問です。 4 2021/10/24 23:32
- 物理学 直交座標系で表す熱伝導方程式と円筒座標系で表す熱伝導方程式の使い分けについて 4 2023/06/30 12:15
- 物理学 熱伝導方程式について 1 2023/03/22 20:55
- 物理学 物理の問題 3 2022/12/21 22:56
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
因数分解って何に役立つの?
-
未知数の数と必要な方程式の数...
-
何年生で習う範囲ですか?
-
エクセルでxを求めたいのですが!
-
2次関数 y=ax2+bx+cのxを求め...
-
3点を通る平面の方程式を行列...
-
xの5乗=1 の答えを教えてく...
-
高2数学の質問です。 円の方程...
-
2次関数と2次方程式の違い
-
円柱と円の方程式
-
円の方程式?円の関数じゃないの?
-
与えられた2数が和と積のとき...
-
数学の初歩的な質問です。 x^2 ...
-
実数係数4次方程式の判別式
-
z^3=1を満たす複素数を答えよ、...
-
数学の3大分野、代数・幾何・解析
-
カシオの関数電卓
-
2x3行列の逆行列の公式
-
数学IIの問題です。 kを定数と...
-
aを定数とするとき、次の方程式...
おすすめ情報