Temperature Field of a Falling Film with Constant Heat Flux Boundary Conditions and Its Influence on the Infrared Detection Range
-
摘要: 深入分析了利用恒热流边界条件下下降液膜(降膜)技术对军事目标垂直进行红外抑制的方法。建立了层流降膜的流动与传热模型,采用能量守恒法,忽略液膜入口段效应,得到了降膜表面的温度分布,简化了求解过程。通过与不同研究者结果做对比,验证了能量守恒法对于评估降膜在充分发展段换热过程的有效性。计算了探测器对采用液膜红外抑制技术后目标的识别距离,结果表明:增大液膜流动的流量、合理控制液膜入口温度可增强对特定目标的红外抑制效果。本人的研究方法和结论对军事目标的红外抑制领域具有中重要的参考意义。Abstract: A falling liquid film (FLM) method for suppressing the vertical infrared of military targets under constant heat flux boundary conditions is analyzed in detail. A flow and heat transfer model of a laminar falling film is established. The temperature distribution on the falling film surface is obtained using the energy conservation method and omitting the liquid film inlet effect, which simplifies the solution process. By comparing the results of different studies, the effectiveness of the energy conservation method for evaluating the heat transfer process of a falling film in a fully developed section is verified. The recognition distance of the detector to the target is calculated using liquid film infrared suppression technology. The results show that increasing the flow rate and sufficiently controlling the inlet temperature of the liquid film enhances the infrared suppression effect on a specific target. The research methods and conclusions of this study have important applications in providing a point of reference for infrared suppression of military targets.
-
0. 引言
垂直或倾斜壁面上的液体降膜流动具有很高的换热效率,在工业和军事领域都有重要的应用。德国学者Nusselt[1]最先建立了液膜流动速度与换热的解析解。在恒壁面热流密度条件和恒壁温条件下,蒋章焰等[2]、Gimbutis等[3-4]、Shmerler和Mudawwa[5-6]、Saouli等[7]、彭友顺等[8-9]对工质为水的过冷液膜进行了研究。其中,在Saouli等[7]的研究中,忽略液膜入口段效应,采用分离变量法,得到了恒热流下倾斜板层流降膜流动的温度分布。在彭友顺等[8-9]的研究中,采用积分法,得到了恒壁温边界条件下竖壁层流降膜流动的温度分布;进一步采用数值模拟法,得到了恒热流和恒壁温边界条件下竖壁层流降膜流动的温度分布的数值解[8-9]。将液体降膜引入到海上大型军事目标红外特征控制领域有着重要的意义。加拿大Davis公司进行了试验研究,并取得了显著的效果[10],但这些求解液膜温度场的方法非常繁琐,而且大多数液膜红外抑制的研究都是针对实验的定性研究,没有定量分析液膜红外抑制效果,不利于实际应用。
本文针对过冷降膜温度场及其对红外特征的抑制效果进行了进一步理论分析,采用能量守恒法,简化了求解液膜温度场的方法,得到了恒热流下层流降膜自由表面的温度分布,定量地研究了液膜红外抑制对探测距离的影响。
1. 控制方程
海上竖壁结构受到日晒、内部热源等影响会呈现出与环境显著的红外特征差异。通过布设液体降膜散热可有效降低其红外特征。
分析液膜在竖壁上流动与传热关系(如图 1所示),假设液膜在竖壁上进入充分发展流动,对于一个流体微元体,其重力ρg(Δ-y)dx与其黏性力μdu/dxdy应该相等,可得:
$$ \rho g\left( {\mathit{\Delta} - y} \right){\text{d}}x = \mu \frac{{{\text{d}}u}}{{{\text{d}}y}}{\text{d}}x $$ (1) 边界条件为:
$$ y=0\text{,}\,\,u=0 $$ (2) 根据液膜流动的Re=4Γ/μ,可得液膜内其他变量的分布情况[4],(3)~(6)式分别是液膜垂向速度、质量流量、厚度和液膜内的平均速度:
$$ u\left( y \right) = \frac{{\rho g{\mathit{\Delta} ^2}}}{\mu }\left[ {\frac{y}{\mathit{\Delta} } - \frac{1}{2}{{\left( {\frac{y}{\mathit{\Delta} }} \right)}^2}} \right] $$ (3) $$ {q_{\rm{m}}} = \frac{{{\rho ^2}g{\mathit{\Delta}^3}}}{{3\mu }} $$ (4) $$ \mathit{\Delta} = {(\frac{{3{\upsilon ^2}}}{{4g}}{\rm{Re}})^{\frac{1}{3}}} $$ (5) $$ {u_{\rm{m}}} = \frac{{\int_0^\mathit{\Delta} {\rho u{\text{d}}y} }}{{\rho \mathit{\Delta}}} = \frac{{g{\mathit{\Delta}^2}}}{{3\upsilon }} $$ (6) 式(3)~(6)中:Γ为单位宽度上液膜的质量流率;μ为流体动力粘度;υ为流体运动粘度;g为重力加速度;ρ为流体密度;Δ为液膜厚度。
忽略液膜的入口段效应,认为液膜处于充分发展状态。对液膜建立建立能量方程有:
$$ u\left( y \right)\frac{{\partial T}}{{\partial x}} = a\frac{{{\partial ^2}T}}{{\partial {y^2}}} $$ (7) 式中:$ T $为温度;$ a $为热扩散率。
根据热平衡分析有:
$$ \begin{gathered} {q_{\rm{w}}}{\text{d}}x + \int_0^\mathit{\Delta} {{\rho _{\rm{w}}}{c_{\rm{w}}}uT{\text{d}}y} = \int_0^\mathit{\Delta} {{\rho _w}{c_{\rm{w}}}uT{\text{d}}y} \hfill \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\; + \frac{{\text{d}}}{{{\text{d}}x}}\left( {\int_0^\mathit{\Delta} {{\rho _{\rm{w}}}{c_{\rm{w}}}uT{\text{d}}y} } \right){\text{d}}x \hfill \\ \end{gathered} $$ (8) 由于为常物性流动,平均温度为:
$$ {T_{\rm{f}}} = \frac{1}{{{u_{\rm{m}}}\mathit{\Delta}}}\int_0^\mathit{\Delta} {uT{\rm{d}}y} $$ (9) 化简后得到:
$$ \frac{{{\text{d}}{T_{\rm{f}}}}}{{{\text{d}}x}} = \frac{{3\upsilon {q_{\rm{w}}}}}{{\rho {c_{\rm{p}}}g{\mathit{\Delta}^3}}} $$ (10) 在恒热流条件下dT/dx=dTf/dx,定义无量纲坐标Y=y/Δ,把式(3)代入能量方程式(7)变为:
$$ \frac{{\rho g{\mathit{\Delta}^2}}}{\mu }\left( {Y - \frac{1}{2}{Y^2}} \right)\frac{{3\upsilon {q_{\rm{w}}}}}{{\rho {c_{\rm{p}}}g{\mathit{\Delta}^3}}} = a\frac{{{\partial ^2}T}}{{{\mathit{\Delta}^2}\partial {Y^2}}} $$ (11) 边界条件为:
$$ x=0\text{,}\,\,T=T_{0} $$ (12) $$ y=0\text{,}-\lambda \frac{\partial T}{\partial y}={q}_{\rm{w}} $$ (13) $$ y=\mathit{\Delta}\text{,}\text{ }\frac{\partial T}{\partial y}=0 $$ (14) 根据式(13)和式(14)可以得到:
$$ T = \frac{{{q_{\rm{w}}}\mathit{\Delta}}}{\lambda }\left( { - \frac{1}{8}{Y^4} + \frac{1}{2}{Y^3} - Y{\text{ + }}C} \right) $$ (15) 式中:C为积分产生的任意常数,并由边界条件获得。壁面温度Tw是跟x有关的,假设Tw已知,即:
$$ Y=0\text{,}T=T_{\rm{w}} $$ (16) 把式(16)代入式(15)得到:
$$ T = {T_{\rm{w}}}(x) + \frac{{{q_{\rm{w}}}\mathit{\Delta}}}{\lambda }\left( { - \frac{1}{8}{Y^4} + \frac{1}{2}{Y^3} - Y} \right) $$ (17) 把速度分布式(3)和温度分布式(17)代入平均温度式(9),可得平均温度与壁面温度的关系为:
$$ {T_{\rm{f}}}(x) = {T_{\rm{w}}}(x) - \frac{{17}}{{35}}\frac{{{q_{\rm{w}}}\mathit{\Delta}}}{\lambda } $$ (18) 根据式(10)和初始条件式(12),可以得到:
$$ {T_{\rm{f}}}(x) = {T_0} + \frac{{3\upsilon {q_{\rm{w}}}}}{{\rho {c_{\rm{p}}}g{\mathit{\Delta}^3}}}x $$ (19) 把式(18)和式(19)代入式(17),得到最后的温度分布为:
$$ \begin{array}{l} T = {T_0} + \frac{{17}}{{35}}\frac{{{q_{\rm{w}}}\mathit{\Delta}}}{\lambda } + \frac{{3\upsilon {q_{\rm{w}}}}}{{\rho {c_{\rm{p}}}g{\mathit{\Delta}^3}}}x \hfill \\ \;\; + \frac{{{q_w}\mathit{\Delta}}}{\lambda }\left( { - \frac{1}{8}{{\left( {\frac{y}{\mathit{\Delta}}} \right)}^4} + \frac{1}{2}{{\left( {\frac{y}{\mathit{\Delta}}} \right)}^3} - \frac{y}{\mathit{\Delta}}} \right) \hfill \\ \end{array} $$ (20) 恒热流条件下对流换热为:
$$ h=q_{\rm{w}}/(T_{\rm{w}}-T_{\rm{f}}) $$ (21) 根据平均温度与壁面温度的关系式(18),可以得到恒热流条件下液膜对流换热系数:
$$ h = \frac{{35}}{{17}}\frac{\lambda }{\mathit{\Delta}} $$ (22) 定义努塞尔数为:
$$ {\text{Nu = }}\frac{{h\mathit{\Delta}}}{\lambda } $$ (23) 把式(22)代入到式(23)得到:
$$ {\text{Nu}} = \frac{{35}}{{17}} $$ (24) 由式(24)可知,在忽略入口段效应及恒热流边界条件的前提下,液膜在壁面上的换热努塞尔数$ {\text{Nu}} $为常数。
从图 2和图 3可以看出在恒热流条件下在换热发展段对流换热和努塞尔数都为一条直线。从图 2中数据可以看出液膜的对流换热系数一般都在3000 W⋅m-2 K-1左右,所以液膜这种高传热系数能有效降低目标表面温度。
2. 液膜对目标表面温度场的影响
海上目标的竖壁结构受内部恒定热源的影响,可认为在壁面上具有恒热流边界条件。假设液膜工质为水,初始条件为液膜入口温度T0=20℃,竖壁热流强度qw=50 kW⋅m-2,质量流率Γ=0.25 kg⋅m-1⋅s-1,本文求得的守恒解,与Saouli等[7]的解和彭友顺等[8]的积分解进行对比,如图 4所示。
从图 4可以看出,本文所求守恒解的变化趋势与积分解和Saouli解一致。在液膜入口段x<0.1 m,守恒解与积分解存在较大差异,但与Saouli解吻合较好;在液膜换热充分发展段x>0.5 m,守恒解与积分解和Saouli解吻合较好,并且3个解的液膜温度均呈线性变化。其原因是由于守恒解模型和Saouli解模型都是忽略了液膜入口段的效应,认为液膜是处于换热充分发展段,但是积分解模型是考虑了液膜的入口段效应。因为入口段效应通常只占液膜流动过程一小部分,所以忽略液膜入口段效应,认为液膜处在换热充分发展段是合理的。所求结果与其他学者的结果验证了守恒法的有效性。
3. 液膜对红外探测作用距离的影响
在液膜作用下,目标壁面温度可以有效抑制。从红外探测的角度,可采用系统噪声等效温差(NETD)来评价液膜对竖壁红外探测作用距离的影响[10-11]:
$$ R = \left\{ {\begin{array}{*{20}{c}} {{{\left[ {\frac{{\mathit{\Delta}{J_{\mathit{\Delta}\lambda }} \cdot {\tau _{\rm{a}}}(R)}}{{{\rm{NETD}} \cdot \omega \cdot \mathit{\Delta}N \cdot {\rm{SNR}} \cdot {K_n}}}} \right]}^{0.5}}}&{(\frac{{{A_{\rm{t}}}}}{{{R^2}}} \geqslant \omega )} \\ {{{\left[ {\frac{{\mathit{\Delta}{J_{\mathit{\Delta}\lambda }} \cdot {\tau _{\rm{a}}}(R)}}{{{\rm{NETD}} \cdot \omega \cdot \mathit{\Delta}N \cdot {\rm{SNR}}}}} \right]}^{0.5}}}&{(\frac{{{A_{\rm{t}}}}}{{{R^2}}} < \omega )} \end{array}} \right. $$ (25) 式中:NETD为系统噪声等效温差;R为作用距离;τa(R)为大气透过率;$ \mathit{\Delta}{J_{\mathit{\Delta}\lambda }} $为目标与背景的波段辐射强度差;At为目标面积;Kn=ωt/ω;ω为瞬时视场立体角;ωt为目标张角;ΔN为辐射亮度差;SNR为系统信噪比。
目标与背景的波段辐射强度差为:
$$ \mathit{\Delta}{J_{\mathit{\Delta}\lambda }} = \frac{{\mathit{\Delta}{W_{\mathit{\Delta}\lambda }} \cdot {A_{\rm{t}}}}}{\pi } = \frac{{{C_{\mathit{\Delta}\lambda }} \cdot {W_{b\mathit{\Delta}\lambda }} \cdot {A_{\rm{t}}}}}{\pi } $$ (26) 由于舰船的背景比较复杂,此处选取20℃海水为目标的背景,计算了恒热流条件下舰船表面铺设液膜后得到不同入口温度、流量下目标的识别距离。
从图 5可以看出,随着入口温度的降低,液膜表面平均温度降低,探测器的作用距离先减小后增大。说明并不是液膜表面温度越低,红外抑制效果越好。红外抑制效果与背景温度密切相关,当目标温度与背景温度差值越小,抑制效果才越好。所以在目标表面铺设液膜的时候,要控制液膜表面温度与目标温度的差值,防止目标成为过冷目标。从图 6可以看出,增大液膜流量,液膜表面平均温度和探测器对目标的识别距离都迅速减小。
一般情况下,在面对敌方性能较为优良成像探测器、利于红外侦察与探测的天气状况(大气透射率较高)以及敌方探测设备对目标具有较低探测概率(探测概率仅为10%)下,我方目标在距离探测器10 km处仍然不被发现,就认定目标达到红外对抗的要求[12]。
所以在本节条件下,入口温度控制在11℃~16℃能满足红外抑制的要求。由于背景温度为20℃,液膜入口温度为20℃,所以加大流量可使液膜表面温度无限接近背景温度。综上所述,最佳方式是同时控制液膜的入口温度和流量,采用降低入口温度和加大流量能更好实现目标的红外抑制。
4. 结论
通过计算恒热流边界条件下降膜温度场及其作用距离,可以得出以下主要结论:
① 恒热流边界条件降膜流动过程中,忽略了液膜的入口段效应,认为液膜处于换热充分发展段,简化了方程的求解过程。采用守恒法所求结果与积分解、Saouli解吻合较好,验证了守恒法的有效性。
② 采用守恒法,得到了液膜自由表面的温度分布以及对流换热系数和努塞尔数。一般情况下对流换热系数都在3000 W⋅m-2⋅K-1左右,说明液膜这种高传热系数能有效降低竖壁目标表面温度。
③ 计算了探测器对采用液膜红外抑制技术后目标的识别距离。在军事目标的竖壁表面合理布膜,可以大大降低目标的可探测性,采用同时降低入口温度和加大流量能更好实现目标的红外抑制。
-
-
[1] Nusselt W. Die Oberflächenkondensation des Wasserdampfes[J]. Z. Ver. Deutch. Ing., 1916, 60: 541-546.
[2] 蒋章焰, 马同泽, 赵嘉琪, 等. 垂直管外降落液膜的流动和传热特性[J]. 工程热物理学报, 1988, 9(1): 70-74. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB198801018.htm JIANG Zhangyan, MA Tongze, ZHAO Jiaqi, et al. Hydrodynamics and heat transfer in thin water films falling down the outside of a vertical tube[J]. Journal of Engineering Thermophysics, 1988, 9(1): 70-74. https://www.cnki.com.cn/Article/CJFDTOTAL-GCRB198801018.htm
[3] Gimbutis G J, Gimbutyte S S, Einkunas S S. Heat transfer in a falling liquid film with large curvature[J]. Heat Transfer Research, 1993, 25(2): 216-219.
[4] Gimbutis G J, Gimbutyte S S, Einkunas S S. Heat transfer of falling liquid film under different heat flux distributions in the film[C]//Proc. of the 1st Baltic Heat Transfer Conference, Goteborg, Sweden, 1991: http://www.beoworks.com/recent-advances-in-heat-transfer-proceedings-of-the-first-baltic.pdf。
[5] Shmerler J A, Mudawwar L. Local evaporative heat transfer coefficient in turbulent free-falling liquid films[J]. Int. J. Heat Mass Transfer, 1988, 31(1): 731-742.
[6] Shmerler J A, Mudawwar L. Local heat transfer coefficient in wavy free-falling turbulent liquid films undergoing uniform sensible heating[J]. Int. J. Heat Mass Transfer, 1988, 33(1): 67-77.
[7] Saouli S. Second law analysis of laminar falling liquid film along an inclined heated plate[J]. International Communications In Heat And Mass Transfer, 2004, 31(6): 879-886. DOI: 10.1016/S0735-1933(04)00074-0
[8] 彭友顺, 杨立, 杜永成. 过冷降膜的温度场及红外隐身应用研究[J]. 红外与激光工程, 2012, 41(10): 2572-2577. DOI: 10.3969/j.issn.1007-2276.2012.10.003 PENG Youshun, YANG Li, DU Yongcheng, et al. Application of subcooling films in temperature field and infrared stealth technology[J]. Infrared and Laser Engineering, 2012, 41(10): 2572-2577. DOI: 10.3969/j.issn.1007-2276.2012.10.003
[9] 彭友顺, 张晓怀, 杨立, 等. 竖壁液膜温度分布数值模拟和红外抑制效果[J]. 哈尔滨工程大学学报, 2013, 34(2): 131-136. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201302002.htm PENG Youshun, ZHANG Xiaohuai, YANG Li, et al. Application of liquid film on ship surface in temperature distribution and infrared suppression effective[J]. Journal of Harbin Engineering University, 2013, 34(2): 131-136. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBG201302002.htm
[10] Thompson J, Vaitekunas D A, Brooking G. Signature management-the pursuit of stealth lowering warship signatures[C/OL]//Electromagnetic and Infrared: Proceedings of SMI Conference-Pursuit of Stealth, 2000: https://www.davis-eng.com/docs/papers/lowering_warship_signatures.11pdf.
[11] 陈翾. 基于海天背景的舰船红外成像仿真与隐身技术研究[D]. 武汉: 海军工程大学, 2010. CHEN Xuan. Research on Infrared Imaging Simulation and Stealth Technology of Ship Based on Sea-sky Background[D]. Wuhan: Naval University of Engineering, 2010.
[12] 路远, 凌永顺, 吴汉平, 等. 地面目标的红外被动测距研究[J]. 红外与毫米波学报, 2004, 23(1): 77-80. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH200401016.htm LU Yuan, LING Yongshun, WU Hanping, et. al. Study on passive distance measurement of ground objects by infrared radiation[J]. J. Infrared Millim. Waves, 2004, 23(1): 77-80. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYH200401016.htm
-
期刊类型引用(6)
1. 黄洁,杨晶媛,徐海华,普凌. 云南锗产业发展现状及对策建议. 中国科技产业. 2024(05): 53-55 . 百度学术
2. 陈应红,彭伟校. 含氟锗精矿氯化蒸馏固氟工艺研究. 广州化工. 2024(18): 154-156 . 百度学术
3. 张臣智,李文琪,李静川,滕祥青,代兵,徐梁格,朱嘉琦. 红外罩失效及修复研究进展. 表面工程与再制造. 2024(05): 11-31 . 百度学术
4. 向兴宇. 锗在太阳能电池中的应用. 现代工业经济和信息化. 2023(01): 170-172 . 百度学术
5. 方迪,苏春风,孙梦荷. 碱熔-电感耦合等离子体原子发射光谱(ICP-AES)法测定锗精矿中锗的含量. 中国无机分析化学. 2023(08): 869-874 . 百度学术
6. 许丰,郭建友,於国兵,徐平,张建. 基于虚拟刻度原理对放射性体源效率刻度方法的发展. 原子核物理评论. 2022(03): 352-358 . 百度学术
其他类型引用(5)