跳到主要内容

九州大学 工学府 機械系専攻 2018年度 熱工学 問2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

円管内を流体が発達した層流で流れ、0≤x0\le x の加熱区間には内半径 r1r_1、外半径 r2r_2、一定熱伝導率 kk の環状ヒータがある。ヒータは一様な体積発熱 qvq_v をもち、外側は完全に断熱されている。軸方向温度勾配は小さく、ヒータの熱伝導は半径方向のみの定常一次元問題とみなせる。(55 点)

環状ヒータと円管内の流れ(公式原図)

ヒータの熱伝導方程式は

1rddr(krdTdr)+qv=0(r1≤r≤r2).(1)\frac1r\frac{d}{dr}\left(kr\frac{dT}{dr}\right)+q_v=0 \qquad(r_1\le r\le r_2). \tag{1}

管内壁温度 T1T_1 と流体の混合平均温度 TmT_m は xx 方向に変化するが、加熱開始点から十分離れた発達域の温度分布は

T−T1Tm−T1=611[3−4(rr1)2+(rr1)4](0≤r≤r1)(2)\frac{T-T_1}{T_m-T_1} =\frac6{11}\left[3-4\left(\frac r{r_1}\right)^2 +\left(\frac r{r_1}\right)^4\right] \qquad(0\le r\le r_1) \tag{2}

である。流体の定圧比熱、熱伝導率、質量流量を cp,kl,m˙c_p,k_l,\dot m とする。Nusselt 数 Nu=hd/kl\mathrm{Nu}=hd/k_l は一定で、d=2r1d=2r_1、熱伝達率 hh は T1−TmT_1-T_m で定義する。

(a) r1r_1 での熱流束 q1q_1 を r1,r2,qvr_1,r_2,q_v で表せ。

(b) (1) を解き、T2−T1T_2-T_1(T2=T(r2)T_2=T(r_2))を r1,r2,qv,kr_1,r_2,q_v,k で表せ。

(c) 入口温度を TinT_{\mathrm{in}} とし、加熱区間の Tm(x)T_m(x) を求めよ。

(d) (2) から流体側の壁面温度勾配を計算し、Nusselt 数を求めよ。

题目描述​

内外半径 r1,r2r_1,r_2 的环状加热器包围圆管,均匀体积发热 qvq_v,外壁绝热、导热率 kk 恒定,忽略加热器轴向导热。管内为充分发展层流,流体导热率 klk_l、比热 cpc_p、质量流率 m˙\dot m;径向温度分布由式 (2) 给出。求内壁热流密度、加热器两壁温差、混合平均温度随 xx 的变化,以及以直径 2r12r_1 为特征长度的 Nusselt 数。

Kai​

(a)​

単位長さ当たりの発熱は πqv(r22−r12)\pi q_v(r_2^2-r_1^2) で、すべて内壁から流体へ伝わる。流体へ向かう向きを正とすると

q1=qv(r22−r12)2r1.\boxed{q_1=\frac{q_v(r_2^2-r_1^2)}{2r_1}.}

rr の正方向を熱流束の正方向に選ぶ場合は、この値に負号を付ける。

(b)​

外壁の断熱条件 T′(r2)=0T'(r_2)=0 を使って (1) を積分すると

T′(r)=qv2k(r22r−r).T'(r)=\frac{q_v}{2k}\left(\frac{r_2^2}{r}-r\right).

再び r1r_1 から r2r_2 まで積分して

T2−T1=qvr222kln⁡r2r1−qv4k(r22−r12).\boxed{T_2-T_1= \frac{q_vr_2^2}{2k}\ln\frac{r_2}{r_1} -\frac{q_v}{4k}(r_2^2-r_1^2).}

(c)​

長さ dxdx の流体の熱収支は

m˙cpdTmdx=πqv(r22−r12).\dot m c_p\frac{dT_m}{dx}=\pi q_v(r_2^2-r_1^2).

Tm(0)=TinT_m(0)=T_{\mathrm{in}} より

Tm(x)=Tin+πqv(r22−r12)m˙cpx.\boxed{T_m(x)=T_{\mathrm{in}}+ \frac{\pi q_v(r_2^2-r_1^2)}{\dot m c_p}x.}

(d)​

(2) を rr で微分して壁面で評価すると

∂T∂r∣r=r1=2411r1(T1−Tm).\left.\frac{\partial T}{\partial r}\right|_{r=r_1} =\frac{24}{11r_1}(T_1-T_m).

流体に向かう熱流束は q1=kl(∂T/∂r)r1=h(T1−Tm)q_1=k_l(\partial T/\partial r)_{r_1}=h(T_1-T_m) なので

h=24kl11r1,Nu=h(2r1)kl=4811≃4.36.\boxed{h=\frac{24k_l}{11r_1},\qquad \mathrm{Nu}=\frac{h(2r_1)}{k_l}=\frac{48}{11}\simeq4.36.}