跳到主要内容

東京大学 工学系研究科 機械工学専攻 2015年度 機械工学第1部 問題1 II

Author​

祭音Myyura (Based on SN's answer refined with GPT 6 Astra)

Description​

厚さ δ\delta、面積 AA の薄膜ヒーターを一定温度 T0T_0 の熱浴に密着させる。接触熱通過率は hh、ヒーターの熱拡散率は aa、熱伝導率は λ\lambda。ヒーター内部の温度 TT は一様で、周囲への熱損失を無視する。

  1. 初め T=T0T=T_0 とし、時刻0から一定の発熱量率 QQ を与える。熱収支式を立て、T=T0+(Q/Ah)[1−exp⁡{−hat/(δλ)}]T=T_0+(Q/Ah)[1-\exp\{-hat/(\delta\lambda)\}] を示せ。
  2. ヒーター温度を Ts0+ΔTs0sin⁡ωtT_{s0}+\Delta T_{s0}\sin\omega t とし、その面に十分細長い棒を完全熱接触させる。棒は側面断熱、熱拡散率 asa_s、遠端温度 Ts0T_{s0} で、ヒーター温度への影響は無視する。接触面を x=0x=0 とし、熱伝導方程式と境界条件を示せ。周期応答 Ts=Ts0+menxsin⁡(ωt+px)T_s=T_{s0}+m e^{nx}\sin(\omega t+px) の m,n,pm,n,p を求めよ。
  3. 間隔 Δx\Delta x の2点で測った位相差 Δϕ\Delta\phi、または温度振幅比 rsr_s を用い、asa_s を求めよ。
  4. ヒーター温度を正弦波状に制御するための角周波数の条件を δ,h,a,λ\delta,h,a,\lambda で示し、理由を説明せよ。

题目描述​

薄膜加热器与恒温热浴接触,已知厚度、接触传热系数、热扩散率和导热系数。建立恒功率加热的集中热容方程并求温度响应。再把加热器温度周期性变化施加于侧面绝热的半无限长细棒,求温度波的衰减与相位,利用两点相位差或振幅比反求热扩散率,并讨论加热器跟随温度调制所需的频率条件。

Kai​

(1)​

体積熱容量は λ/a\lambda/a なので

λaAδdTdt=Q−hA(T−T0),T(0)=T0.\frac{\lambda}{a}A\delta\frac{\mathrm dT}{\mathrm dt}=Q-hA(T-T_0),\qquad T(0)=T_0.

時定数 τ=δλ/(ha)\tau=\delta\lambda/(ha) を用いると

T=T0+QAh(1−e−t/τ).\boxed{T=T_0+\frac{Q}{Ah}(1-e^{-t/\tau})}.

(2)​

棒の微小区間の熱収支とフーリエの法則より

∂Ts∂t=as∂2Ts∂x2,Ts(0,t)=Ts0+ΔTs0sin⁡ωt,Ts(∞,t)=Ts0.\boxed{\frac{\partial T_s}{\partial t}=a_s\frac{\partial^2T_s}{\partial x^2}},\quad T_s(0,t)=T_{s0}+\Delta T_{s0}\sin\omega t,\quad T_s(\infty,t)=T_{s0}.

与式を代入すると n2=p2n^2=p^2、2asnp=ω2a_snp=\omega。遠方で減衰する条件 n<0n<0 より、定常周期応答は

m=ΔTs0,n=p=−ω2as,\boxed{m=\Delta T_{s0},\qquad n=p=-\sqrt{\frac{\omega}{2a_s}}},
Ts=Ts0+ΔTs0e−kxsin⁡(ωt−kx),k=ω2as.\boxed{T_s=T_{s0}+\Delta T_{s0}e^{-kx}\sin(\omega t-kx),\qquad k=\sqrt{\frac{\omega}{2a_s}}}.

(3)​

上流側振幅を下流側振幅で割った比を rsr_s とすると

Δϕ=kΔx,ln⁡rs=kΔx.\Delta\phi=k\Delta x,\qquad \ln r_s=k\Delta x.

したがって

as=ω(Δx)22(Δϕ)2=ω(Δx)22(ln⁡rs)2.\boxed{a_s=\frac{\omega(\Delta x)^2}{2(\Delta\phi)^2} =\frac{\omega(\Delta x)^2}{2(\ln r_s)^2}}.

(4)​

発熱量の変調に温度を十分追従させる条件は、熱応答時間が変調時間より十分短いこと、すなわち

ωτ=ωδλha≪1.\boxed{\omega\tau=\frac{\omega\delta\lambda}{ha}\ll1}.

このとき T−T0≃Q(t)/(hA)T-T_0\simeq Q(t)/(hA) で、振幅の減衰と位相遅れが小さい。薄膜内部の温度一様性にも ωδ2/a≪1\omega\delta^2/a\ll1 を要する。

なお、線形な集中熱容量モデルでは、任意の ω\omega でも過渡項の消滅後の正弦波入力に対する応答は正弦波であり、その振幅比は 1/1+(ωτ)21/\sqrt{1+(\omega\tau)^2}、位相遅れは tan⁡−1(ωτ)\tan^{-1}(\omega\tau) である。

Reference​