許鵬飛 唐 豪

?

空調(diào)服近體流動與換熱研究

許鵬飛 唐 豪

（南京航空航天大學能源與動力學院 南京 210016）

為了探討空調(diào)服系統(tǒng)的流動與換熱規(guī)律，基于人體外形創(chuàng)建了體表-服裝幾何模型。采用數(shù)值仿真的方法研究了橫向型管路系統(tǒng)的流動分配特征，探討了體表平均溫度、內(nèi)外溫差關(guān)于入口速度、入口溫度的變化關(guān)系。結(jié)果顯示：管路質(zhì)量流量分配存在突變現(xiàn)象；體表上半部分溫度低于下半部分溫度；體表平均溫度、內(nèi)外溫差與入口速度成二次函數(shù)關(guān)系，與入口溫度成線性關(guān)系；入口溫度每升高1℃，體表平均溫度升高0.34℃，而內(nèi)外溫差降低0.33℃。

空調(diào)服；流動換熱；幾何模型；數(shù)值仿真

0 引言

進入新時代，人民對美好生活的向往已成為全社會所為之奮斗的目標。在這樣的大背景下，軍警、消防、采礦等特種行業(yè)惡劣的工作環(huán)境越來越得到人們的關(guān)注。特種行業(yè)從業(yè)人員常常須忍受37℃以上的高溫，例如深井采礦，有的工作面上隅角溫度可達37～40℃[1]。身處如此環(huán)境中開展工作，作業(yè)人員的身體健康安全時刻面臨著威脅。環(huán)境溫度過高，人體產(chǎn)熱不能及時散發(fā)，只能積蓄在體內(nèi)，導致體溫上升，產(chǎn)生熱應(yīng)激，這樣一來作業(yè)人員注意力下降，體力逐漸消耗殆盡，不僅降低了工作效率，而且更容易引發(fā)生理異常，重者甚至威脅生 命[2,3]。因此，有必要研制出一套空調(diào)服系統(tǒng)，改善人體微環(huán)境，即使身處惡劣境況也能保障人體熱舒適性，從根本上杜絕了安全事故的發(fā)生。

近年來，國內(nèi)外相關(guān)學者對空調(diào)服進行了大量的研究，取得了碩果累累的成就。劉靜[4]等人提出了一種基于風扇陣列驅(qū)動的空調(diào)服，其原理是在衣服中布置若干微小型風扇，每個風扇體積介于1mm3到1cm3，風扇運轉(zhuǎn)，增強空氣與人體表面的對流換熱，達到降溫冷卻的目的。曾彥彰[5]等人又對該種空調(diào)服進行了理論計算和實驗測試，結(jié)果表明：在室內(nèi)環(huán)境下，相比于無風扇時36℃的體表溫度，身穿空調(diào)服可將體表溫度降至29℃；而在太陽輻射環(huán)境下，無風扇時體表溫度是40.3℃，身穿空調(diào)服可將體表溫度降至30℃。李伊潔[6]等人介紹了國內(nèi)外通用的熱舒適評價標準，闡述了人體熱舒適研究的發(fā)展趨勢。王濤[7]等人在實驗中嘗試向液冷服注入TiO2納米流體充當冷卻液，并將其散熱性能與冷卻液為水的情形作對比，發(fā)現(xiàn)納米流體冷卻服散熱性能與流體質(zhì)量流量、進口溫度和納米顆粒體積分數(shù)有關(guān)。Maurissa[8]等人針對現(xiàn)有冷卻服裝的固有缺陷，提出將熱電制冷器應(yīng)用于氣冷服，并對此進行了實驗測試。結(jié)果表明：熱電制冷裝置可在低功率狀態(tài)下順利運行，同時空氣溫度可被降低約10℃，完全可以達到預期目的。王軍[9]建立了人體微環(huán)境模型，通過數(shù)值模擬的方法對四種典型氣流組織模式進行了分析，并得到了各自的室內(nèi)人體微環(huán)境需求目標的實現(xiàn)特征。

本文基于ANSYS Fluent數(shù)值仿真，重點研究了空調(diào)服近體流動與換熱規(guī)律，推導了體表平均溫度、內(nèi)外溫差關(guān)于入口速度、入口溫度的函數(shù)關(guān)系，確立了設(shè)計和優(yōu)化空調(diào)服應(yīng)遵循的基本原則，同時也為個體冷卻裝備的研究提供思路和借鑒。

1 模型及數(shù)值方法

1.1 物理模型

圖1 人體幾何模型

本文只研究人體上身軀干，不涉頭部及四肢。采用機械設(shè)計軟件UG NX10.0建立人體幾何模型，如圖1所示。

圖1為成年男子上身模型，所涉及的尺寸數(shù)據(jù)如表1所示。

表1 人體模型尺寸數(shù)據(jù)

基于人體模型架構(gòu)，進而建立體表-服裝幾何模型，如圖2所示。模型中包括外基礎(chǔ)服裝、空氣夾層、體表和管路。

人體與外界環(huán)境是通過皮膚進行熱量交換，因此將皮膚，即人體表面設(shè)定為幾何模型內(nèi)邊界，而不考慮人體內(nèi)部的物理、化學變化。圖2中模型的最外層為外基礎(chǔ)服裝，為了模型的簡化，本文在此忽略其厚度。服裝層與體表之間留有30mm的空氣夾層，用于配置管路系統(tǒng)。本文選用橫向型管路系統(tǒng)，如圖3所示。

圖3 橫向型管路

管路由4根橫向回路及1根縱向主軸構(gòu)成，每根管路直徑均為22mm，管路表面共開有19個小孔，孔徑16mm。按照從左至右，自上而下的原則依次為每個出口編號，即從1號至19號。冷氣流通過各個小孔流出，與皮膚進行對流換熱，促進人體散熱，從而降低體溫。

1.2 數(shù)學模型

本文針對穩(wěn)態(tài)、常物性、不可壓縮流體，數(shù)學模型控制方程為[10]：

湍動能方程：

湍動能耗散方程：

1.3 邊界條件

體表-服裝幾何模型的外表面是忽略厚度的服裝層，它與外界環(huán)境發(fā)生對流換熱，在此定義環(huán)境溫度為37℃，服裝表面對流換熱系數(shù)為10W/(m2·K)；體表源源不斷地向外界散發(fā)熱量，定義其單位面積散熱量為50W/m2；橫向型管路冷氣入口定義為速度入口，即velocity-inlet；幾何模型中的頸口、袖口、底面出口均定義為壓力出口，即pressure-outlet。

1.4 數(shù)值方法

本文數(shù)值方法基于有限容積法（Finite Volume Method, FVM），采用ANSYS ICEM CFD對模型劃分非結(jié)構(gòu)化網(wǎng)格，網(wǎng)格總數(shù)為3314677。湍流模型采用標準模型。壓力和速度的耦合關(guān)系采用Simple算法。湍動能方程及湍流耗散率方程的對流項采用二階迎風格式進行離散，動量方程與能量方程的對流項采用Quick格式進行離散。

2 結(jié)果與討論

本文采用ANSYS Fluent 16.0作為仿真工具，分別模擬了空調(diào)服在不同入口速度、入口溫度條件下近體流動與換熱狀況，歸納總結(jié)了空調(diào)服使用過程中的規(guī)律特征。

2.1 體表溫度分布

如圖4所示為入口速度v=14m/s、入口溫度T=18℃時的體表溫度分布云圖。

（a）前身

（b）后背

圖4 體表溫度分布

Fig.4 Temperature distribution on body surface

從圖中可以發(fā)現(xiàn)，溫度最高區(qū)域出現(xiàn)在體表下半部分，位于腰部兩側(cè)；溫度較低區(qū)域出現(xiàn)在體表上半部分，位于前胸及后背脊椎位置。這種溫度分布特點與管路結(jié)構(gòu)設(shè)計有關(guān)。

結(jié)合圖5分析，即管路流線圖，氣流進入主軸管后，更多地匯聚于最上側(cè)回路，因此從上側(cè)回路出口流出的氣流占多數(shù)。

圖5 管路流線圖

體表溫度分布特點也可從管路各出口質(zhì)量流量分配得到印證。

圖6 管路質(zhì)量流量分配

2.2 入口速度對制冷性能的影響

圖7 體表平均溫度隨入口速度的變化

控制入口溫度T不變，改變?nèi)肟谒俣?i>v，得到T在不同定值條件下v對空調(diào)服制冷性能的影響。圖7展示了體表平均溫度隨入口速度的變化情況。

隨著入口速度的增大，體表溫度逐漸降低，但降低的幅度越來越小，即體表溫度下降越緩慢。以入口溫度T=14℃為例，當入口速度v從8m/s增至9m/s時，體表平均溫度T降低0.94℃，而當入口速度v從13m/s增至14m/s時，體表平均溫度T降低0.43℃，降幅僅為前者的45.74%。

圖8展示了內(nèi)外溫差隨入口速度的變化情況。

圖8 內(nèi)外溫差隨入口速度的變化

內(nèi)外溫差Δ是指衣內(nèi)空氣層平均溫度與外界環(huán)境溫度的差值。觀察圖8發(fā)現(xiàn)，隨著入口速度的增大，內(nèi)外溫差逐漸升高，但上升的幅度越來越小，即內(nèi)外溫差上升越緩慢，這一點與體表平均溫度的變化情況相反但類似。同樣以入口溫度T=14℃為例，當入口速度v從8m/s增至9m/s時，內(nèi)外溫差Δ升高0.55℃，而當入口速度v從13m/s增至14m/s時，內(nèi)外溫差Δ升高0.26℃，增幅僅為前者的47.27%。

根據(jù)變化曲線圖，利用多項式擬合的方法，分別得到體表平均溫度、內(nèi)外溫差隨入口速度的變化關(guān)系，即：

其中，b、ΔT均為截距，大小與入口溫度有關(guān)。

2.3 入口溫度對制冷性能的影響

控制入口速度v不變，改變?nèi)肟跍囟?i>T，得到v在不同定值條件下T對空調(diào)服制冷性能的影響。

圖9 體表平均溫度隨入口溫度的變化

圖10 內(nèi)外溫差隨入口溫度的變化

圖9、10所示分別為體表平均溫度、內(nèi)外溫差隨入口溫度的變化情況，很容易發(fā)現(xiàn)，兩者均隨入口溫度線性變化。入口溫度每升高1℃，體表平均溫度將升高0.34℃，而內(nèi)外溫差將下降0.33℃。

最后，通過線性擬合的方法，得到相應(yīng)的函數(shù)表達式，即：

T=0.3381T+b（8）

=-0.3297T+b（9）

其中，b、ΔT均為截距，大小與入口速度有關(guān)。

3 總結(jié)與展望

本文基于人體外形創(chuàng)建了體表-服裝幾何模型，通過數(shù)值仿真研究了空調(diào)服近體流動與換熱規(guī)律，所得結(jié)論主要有：

（1）體表溫度分布呈現(xiàn)“上冷下熱”的特點；

（2）管路系統(tǒng)各出口質(zhì)量流量分配存在突變現(xiàn)象；

（3）體表平均溫度、內(nèi)外溫差關(guān)于入口速度成二次函數(shù)關(guān)系，且入口速度越大，其對空調(diào)服制冷性能的影響越小；

（4）體表平均溫度、內(nèi)外溫差關(guān)于入口溫度成線性關(guān)系，入口溫度每升高1℃，體表平均溫度將升高0.34℃，而內(nèi)外溫差將下降0.33℃。

由于本文只討論了橫向型管路系統(tǒng)，且表面開口數(shù)量固定，因此，有關(guān)管路系統(tǒng)及開口數(shù)量對空調(diào)服制冷性能的影響有待進一步的研究。

[1] 袁亮.淮南礦區(qū)礦井降溫研究與實踐[J].采礦與安全工程學報,2007,(3):298-301.

[2] Malchaire J, Kampmann B, Havenith G, et al. Criteria for estimating acceptable exposure times in hot working environments: a review[J]. International Archives of Occupational & Environmental Health, 2000,73(4):215- 220.

[3] Nag P K, Nag A, Ashtekar S P. Thermal limits of men in moderate to heavy work in tropical farming[J]. Indus-trial Health, 2007,45(1):107.

[4] 中國科學院理化技術(shù)研究所.基于微/納米風扇陣列驅(qū)動的人體降溫空調(diào)服[P].北京:CN200520114514.6. 2007-03-28.

[5] 曾彥彰,鄧中山,劉靜,等.基于微型風扇陣列系統(tǒng)的人體降溫空調(diào)服[J].紡織學報,2007,28(6):100-105.

[6] 李伊潔,劉何清,劉天宇,等.國內(nèi)外通用室內(nèi)環(huán)境熱舒適評價標準的分析與比較[J].制冷與空調(diào),2017,31(1): 14-22.

[7] 王濤,王亮,林貴平,等.TiO2納米流體在液冷服上的應(yīng)用實驗研究[J].浙江大學學報(工學版),2016,(4):681-690.

[8] D’Angelo M, D’Angelo J, Almajali M, et al. Augmented cooling vest system subassembly:Design and analysis[J]. Energy Conversion & Management, 2014,79(79):140- 145.

[9] 王軍.室內(nèi)人體微環(huán)境的需求目標實現(xiàn)特征數(shù)值分析[J].制冷與空調(diào),2016,30(3):278-282.

[10] 陶文銓.數(shù)值傳熱學(第2版)[M].西安:西安交通大學出版社,2001:1-4.

Research on the Flow and Heat Transfer Around Body of Air-Conditioning Garment

Xu Pengfei Tang Hao

( College of Energy and Power Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 210016 )

In order to explore the regulation of the flow and heat transfer in air-conditioning garment system, a geometric model of body surface-garment was created based on the body shape. By numerical simulation, the flow distribution characteristic in horizontal pipeline system was studied and the relationships between average temperature on body surface, temperature difference between inside and outside and inlet velocity, inlet temperature were discussed. The results show that: there is mutation phenomenon during the flow distribution of pipeline; the temperature of upper half part on body surface is lower than the temperature of lower half part on body surface; the relationships between average temperature on body surface, temperature difference and inlet velocity are quadratic functions, as linear relationships with inlet temperature; when inlet temperature increases by 1℃, average temperature on body surface will increase by 0.34℃, temperature difference will decrease by 0.33℃.

air-conditioning garment; flow and heat transfer; geometric model; numerical simulation

O551

B

南京市321計劃重點項目（2013B03004）

許鵬飛（1992-），男，碩士研究生，研究方向為空調(diào)服及其應(yīng)用，E-mail：1660494071@qq.com

唐 豪（1960-），男，教授，博士生導師，主要從事近體調(diào)溫技術(shù)的研究，E-mail：hao.tang@nuaa.edu.cn

2018-01-16

1671-6612（2018）05-555-05