阿依先古力·熱合曼 居來提·買提肉孜 佐熱姆·拜科日 李乾坤 于洋

摘 要：對(duì)沙療室傳熱系統(tǒng)的傳熱效果和沙體層的溫度分布情況，按照現(xiàn)有的最佳加熱系統(tǒng)自然對(duì)流-多向加熱采用標(biāo)準(zhǔn)k-？著湍流模型建立大型和小型兩種規(guī)格尺寸的三維物理模型，并按照基礎(chǔ)實(shí)驗(yàn)數(shù)據(jù)與模擬結(jié)果，不同邊界條件基礎(chǔ)下對(duì)系統(tǒng)的溫度場(chǎng)、速度場(chǎng)進(jìn)行數(shù)值模擬和數(shù)據(jù)對(duì)比。經(jīng)過研究得出了自然對(duì)流-多向加熱對(duì)大型沙療室加熱系統(tǒng)的傳熱效率無明顯的高，沙體層垂直溫度分布較均勻，以及強(qiáng)迫對(duì)流對(duì)多向加熱的系統(tǒng)溫度變化沒有明顯影響。以此大型沙療室無法達(dá)到溫度分布更為均勻等原因，將沙療室的布置改成為由沙療床組成的治療場(chǎng)所，通過對(duì)比分析得知了沙療床對(duì)沙療的過程和傳熱效率帶有所提高。

關(guān)鍵詞：k-？著模型；傳熱效率；數(shù)值模擬；沙療床

中圖分類號(hào)：TP391.9 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2095-2945（2019）21-0001-05

Abstract： In view of the heat transfer effect of the heat transfer system of the sand therapy room and the temperature distribution of the sand layer， based on the existing optimal heating system natural convection-multidirectional heating， the standard k-？著 turbulence model is used to establish three-dimensional physical models of large and small sizes. According to the basic experimental data and simulation results， the temperature field and velocity field of the system are simulated and compared under different boundary conditions. Through the study， it is concluded that the natural convection-multidirectional heating has no obvious effect on the heat transfer efficiency of the large sand therapy chamber heating system， the vertical temperature distribution of the sand layer is uniform， and the forced convection has no obvious effect on the temperature change of the multidirectional heating system. Due to the reason that the large sand therapy room cannot achieve more uniform temperature distribution and other reasons， the layout of the sand therapy room is changed into a treatment place composed of sand therapy bed. Through comparative analysis， it is found that the sand therapy bed has improved the process and heat transfer efficiency of sand therapy.

Keywords： k-？著 model； heat transfer efficiency； numerical simulation； sand therapy bed

引言

沙療目前是國(guó)際健身治病潮流而且有顯著療效的新療法。埋沙治療中[1]，將患者的膝蓋部位或者全體埋入在沙體中，因此患病處受到沙粒中各類微量元素、礦物質(zhì)和沙粒之間的相互移動(dòng)刺激患處，由此產(chǎn)生的熱量傳送給皮膚、肌肉、股骨等組織，推動(dòng)患處的血液循環(huán)[2]，在世界傳統(tǒng)民族醫(yī)學(xué)中沙療的這種通過傳熱來治病的特性具有重要地位。具體來說是吐魯番維醫(yī)中通過埋沙來治關(guān)節(jié)病歷時(shí)長(zhǎng)久。吐魯番具備著世界上最好的沙療資源和沙療場(chǎng)所，吐魯番無論是在氣候溫度的優(yōu)勢(shì)方面或者是含有豐富的磁鐵礦物質(zhì)和微量元素的田沙方面，有其他地方所不能具備的特殊優(yōu)勢(shì)[3]，是我國(guó)著名而唯一的沙療圣地。

盡管我國(guó)吐魯番是著名的沙療圣地[4]，沙療因受到氣候、環(huán)境、時(shí)間等自然因素的影響，一直沒能實(shí)現(xiàn)隨時(shí)能夠進(jìn)行治療的需求。到目前為止有關(guān)研究人員對(duì)優(yōu)化加熱系統(tǒng)并且提高系統(tǒng)傳熱效果做相關(guān)性研究分析，通過數(shù)值模擬和實(shí)驗(yàn)研究得出了能讓加熱系統(tǒng)傳熱效率提升的自然對(duì)流-多向加熱方式。現(xiàn)如今為世界傳統(tǒng)民族醫(yī)學(xué)中的精髓之一埋沙治療又有了功能最為優(yōu)異的新治療結(jié)構(gòu)沙療床。沙療床在加熱效率、能耗、加熱速度和經(jīng)濟(jì)等方面更有優(yōu)勢(shì)。因此沙療床在未來會(huì)是埋沙治療中的應(yīng)用范圍最廣泛，效率最優(yōu)異的一種沙療結(jié)構(gòu)。

1 實(shí)驗(yàn)與數(shù)據(jù)分析

1.1 實(shí)驗(yàn)

目前對(duì)埋沙治療有關(guān)研究人員建立小型沙療平臺(tái)進(jìn)行研究得出了最佳的加熱系統(tǒng)-自然對(duì)流-多向加熱方式。按照現(xiàn)有的埋沙治療加熱系統(tǒng)和有關(guān)研究結(jié)果，在吐魯番大型室內(nèi)沙療室經(jīng)過對(duì)5cm～30cm深度的沙層進(jìn)行做實(shí)際實(shí)驗(yàn)和模擬實(shí)驗(yàn)采集了溫度數(shù)據(jù)，從實(shí)驗(yàn)數(shù)據(jù)發(fā)現(xiàn)加熱時(shí)間同樣的控制在7小時(shí)到9小時(shí)之內(nèi)實(shí)際實(shí)驗(yàn)數(shù)據(jù)和模擬實(shí)驗(yàn)數(shù)據(jù)基本吻合，一直均勻增長(zhǎng)，溫度差較小，從而得知大型室內(nèi)沙療室的加熱溫度情況基本達(dá)到了人體最需要的最合適的溫度42.5℃[5]。以下為在大型沙療室數(shù)據(jù)采集當(dāng)中的數(shù)據(jù)采集示意圖（如圖1）和實(shí)驗(yàn)數(shù)據(jù)（如表1）和（如表2）。

1.2 數(shù)據(jù)分析

通過在實(shí)際沙療室做模擬和實(shí)際實(shí)驗(yàn)測(cè)出了以上數(shù)據(jù)。從這些數(shù)據(jù)可以看出同一個(gè)時(shí)間段沙子深度同樣厚度時(shí)模擬和實(shí)際實(shí)驗(yàn)數(shù)據(jù)之間的差別不大，隨著加熱時(shí)間的增高深度不同的沙體層溫度也逐漸的增高，雖然模擬實(shí)驗(yàn)所得出的數(shù)據(jù)稍微高一些，但基本接近實(shí)際實(shí)驗(yàn)的數(shù)據(jù)，因此我們初步能斷定自然對(duì)流-多向加熱系統(tǒng)大型沙療室基本達(dá)到最佳的加熱效果，各項(xiàng)數(shù)據(jù)基本合理[2]。為更明確的分析實(shí)驗(yàn)數(shù)據(jù)的變化情況，下面建立了兩項(xiàng)實(shí)驗(yàn)數(shù)據(jù)的XL曲線圖（如圖2）。

從曲線圖可以明確的看到，系統(tǒng)加熱到7小時(shí)至9小時(shí)之間的溫度變化情況模擬實(shí)驗(yàn)和實(shí)際實(shí)驗(yàn)漸漸上升的狀態(tài)，沙體垂直溫度層的一直均勻上升狀態(tài)可以滿足埋沙治療所需要，對(duì)人體無損害沙體溫度。

2 數(shù)值計(jì)算

2.1 大型沙療室的數(shù)值計(jì)算

2.1.1 建立物理模型與生成網(wǎng)格

按照現(xiàn)有的小型沙療平臺(tái)的加熱系統(tǒng)，建立幾何尺寸為5m*8m*2.7m的大型沙療室物理模型，并對(duì)它進(jìn)行劃分網(wǎng)格、定義邊界條件用FLUENT軟件分析出沙療室加熱系統(tǒng)的溫度場(chǎng)和速度場(chǎng)變化。為確保和提高數(shù)值計(jì)算結(jié)果的精確度，對(duì)物理模型要定義合理的邊界條件。數(shù)值計(jì)算當(dāng)中采用能量方程、DO模型以及標(biāo)準(zhǔn)k-？著湍流模型[6-7]。大型沙療室的物理模型（如圖3）。對(duì)于劃分網(wǎng)格傳熱系統(tǒng)物理模型在現(xiàn)有實(shí)驗(yàn)平臺(tái)的基礎(chǔ)之下，為了避免矩陣式[8]計(jì)算所帶來的麻煩、在網(wǎng)格劃分技術(shù)方面采用了沒有矩陣式計(jì)算的方式-結(jié)構(gòu)化和非結(jié)構(gòu)化相互結(jié)合方式。因此對(duì)沙體層相應(yīng)的是六面體形狀為主的結(jié)構(gòu)化網(wǎng)格方式；對(duì)于空氣層也同樣是六面體結(jié)構(gòu)化網(wǎng)格方式，而輻射加熱燈周圍的近距離應(yīng)用非結(jié)構(gòu)化網(wǎng)格方式，并且消除“臺(tái)階效應(yīng)”[9]對(duì)物理模型制造的網(wǎng)格劃分精度方面的影響，這樣網(wǎng)格劃分的效果會(huì)更加的接近貼體坐標(biāo)法（BFC）[10]的標(biāo)準(zhǔn)。網(wǎng)格劃分（如圖4）。

2.1.2 定義邊界條件

網(wǎng)格劃分的下一步是對(duì)物理模型進(jìn)行制定邊界條件，將室內(nèi)沙療室沙體層、沙體層側(cè)邊以及頂部熱源一律可定義為wall，對(duì)于計(jì)算區(qū)域的屬性設(shè)置定義的時(shí)候，系統(tǒng)內(nèi)地面的沙體層可定義為Solid，屋頂邊上的空氣層可定義為Fluid。然后將物理模型導(dǎo)入FLUENT軟件，為了保證數(shù)值計(jì)算的順利完成，對(duì)網(wǎng)格劃分進(jìn)行檢查。為了提高數(shù)值計(jì)算的準(zhǔn)確度要注意這幾個(gè)部分：（1）實(shí)驗(yàn)平臺(tái)中沙子底層

隔熱膜有防火作用，所以數(shù)值計(jì)算中各面設(shè)置成絕熱體；（2）加熱后氣體溫度差導(dǎo)致的密度差會(huì)受重力的影響，因

此浮升力不能忽略；（3）將風(fēng)扇簡(jiǎn)化成面再計(jì)算能實(shí)現(xiàn)極

高精確度的結(jié)果；（4）系統(tǒng)內(nèi)空氣服從牛頓流體條件，因此各面上空氣也全服從不滑動(dòng)的條件；（5）加熱系統(tǒng)是全密

封狀態(tài)，可不用考慮實(shí)驗(yàn)平臺(tái)里小縫隙導(dǎo)致的對(duì)流問題；以下為提高數(shù)值計(jì)算精確度對(duì)系統(tǒng)設(shè)定的邊界條件流程圖（如圖5）。

輻射選項(xiàng)選用為DO模型，并且將各參數(shù)設(shè)置為默認(rèn)值。在沙療中研究的是因?yàn)樯丑w層溫度場(chǎng)跟隨加熱時(shí)間所發(fā)生的變化，之所以以非定常計(jì)算為主要計(jì)算方法。設(shè)置邊界條件時(shí)，室內(nèi)沙療室周圍四面墻和底層與外面環(huán)境不會(huì)互換熱量，可一律定義為絕熱面，屋頂16支加熱燈為輻射性熱源。本文當(dāng)中所有的數(shù)值計(jì)算方式一律選用SIMPLE壓力計(jì)算法[11]。

2.1.3 溫度場(chǎng)和速度場(chǎng)分析

按照有關(guān)研究人員通過重復(fù)實(shí)驗(yàn)研究出的加熱效率最佳的自然對(duì)流-多向加熱系統(tǒng)的現(xiàn)狀，對(duì)幾何尺寸為5m*8m*2.7m的大型沙療室物理模型達(dá)到實(shí)際沙療場(chǎng)所需的溫度為止進(jìn)行數(shù)值模擬。因?yàn)樵趯?shí)際沙療場(chǎng)所通過實(shí)際實(shí)驗(yàn)所得的數(shù)據(jù)當(dāng)中達(dá)到理想溫度所需的時(shí)間為7小時(shí)到9小時(shí)。因此，數(shù)值計(jì)算當(dāng)中同樣的把時(shí)間控制在7小時(shí)到9小時(shí)之內(nèi)，并且對(duì)系統(tǒng)溫度場(chǎng)變化云圖進(jìn)行分析。加熱系統(tǒng)溫度場(chǎng)云圖（如圖6）。

根據(jù)數(shù)值計(jì)算計(jì)算結(jié)果對(duì)空氣層并且沙體層的溫度場(chǎng)變化進(jìn)行分析，從溫度場(chǎng)云圖分布情況可看到：系統(tǒng)內(nèi)因?yàn)橛休椛浼訜釤簟⒓訜岚?、沙表層的電暖線，空氣層的溫度分布情況比較集中在中間。同樣的將加熱時(shí)間設(shè)定在7小時(shí)到9小時(shí)之間，數(shù)值計(jì)算到7小時(shí)計(jì)算結(jié)果無變化，顯示出了較好的加熱效果。但沙體層垂直溫度梯度情況比較明顯，系統(tǒng)溫度分布狀態(tài)是偏高偏低。從此可以斷定大型沙療室溫度場(chǎng)情況跟實(shí)際實(shí)驗(yàn)對(duì)比有一定的出入和欠缺。

同樣的理論基礎(chǔ)之下，對(duì)大型沙療室物理模型把時(shí)間控制在9小時(shí)之內(nèi)的條件下進(jìn)行數(shù)值計(jì)算并且對(duì)系統(tǒng)速度場(chǎng)變化進(jìn)行分析。以下為加熱系統(tǒng)速度場(chǎng)云圖（如圖7）：

根據(jù)數(shù)值計(jì)算計(jì)算結(jié)果對(duì)空氣層并且沙體層的速度場(chǎng)變化進(jìn)行分析，從速度場(chǎng)云圖分布情況可看到：因?yàn)檩椛浼訜釤艉蜕丑w表面層溫度稍微高一些，以此在空氣中存在明顯的浮升力?？諝庵械臒釟饬魃仙问捷^為不規(guī)則，再因?yàn)橛龅街車鷫Φ淖柚?，在上升?dāng)中且會(huì)出現(xiàn)較為大一些的渦流形狀。多向加熱時(shí)在系統(tǒng)內(nèi)溫度分布情況雖然比較均勻、但也存在一些溫度梯度狀況，因此輻射加熱燈近距離流體跟隨遇熱溫上升，出現(xiàn)的渦流形狀比較明顯，以此確定7小時(shí)的加熱時(shí)間對(duì)系統(tǒng)內(nèi)空氣流動(dòng)速度的影響稍微大一些。

2.2 沙療床的數(shù)值計(jì)算與特性

2.2.1 沙療床的數(shù)值計(jì)算

關(guān)于沙療到今日為止有了各種先進(jìn)、無欠缺并且實(shí)惠方便的新科技。首先是以前的室外沙療變成了人工加熱的室內(nèi)沙療室；其次加熱系統(tǒng)存在欠缺的室內(nèi)沙療室優(yōu)化成了具有最佳加熱系統(tǒng)的室內(nèi)沙療室；然后在此正要研究的是-加熱方式為底部加熱和四周加熱方式的、具有完美加熱系統(tǒng)的、設(shè)備很先進(jìn)、體積適當(dāng)?shù)男　o欠缺、實(shí)惠又便攜式的沙療床。它還有稱之為“簡(jiǎn)單的健康大使”。

因此對(duì)沙療床做進(jìn)一步的數(shù)值計(jì)算分析。以下為用UG7.0軟件中建立的幾何尺寸為2.2m*1.2m*0.6m的沙療床物理模型（如圖8）和與沙療室同樣邊界條件下用Gambit軟件進(jìn)行劃分網(wǎng)格結(jié)果圖（如圖9）。

同樣的理論基礎(chǔ)之下，對(duì)沙療床進(jìn)行數(shù)值模擬分析。以下圖為沙療床的數(shù)值模擬分析結(jié)果，溫度場(chǎng)云圖（如圖10）和速度場(chǎng)云圖（如圖11）所示。

從溫度場(chǎng)云圖分布情況可得知：空氣層的溫度分布情況非常的均勻，沙體層垂直方向無任何溫度梯度情況，這說明沙子的加熱速度平穩(wěn)，傳熱效率高。從速度場(chǎng)云圖分布情況可得知：氣流上升形式很規(guī)則，無出現(xiàn)明顯的渦流形狀，加熱過程中空氣流速很均勻，模擬計(jì)算結(jié)果很理想。

2.2.2 沙療床的特性

相比于沙療室，沙療床在各方面更有優(yōu)勢(shì)，比如：（1）外觀優(yōu)美精致，總體積小、通常標(biāo)準(zhǔn)的尺寸規(guī)格為2.2m*1.2m*0.6m、方便移動(dòng)；（2）衛(wèi)生，不會(huì)污染環(huán)境；（3）床體材料純屬木質(zhì)，對(duì)人體和系統(tǒng)沒有負(fù)面影響；（4）壽命長(zhǎng)；（5）加熱速度快；

（6）電壓220V、能耗1800W、功率適當(dāng)、節(jié)省電和時(shí)間；（7）床內(nèi)沙子的更換周期長(zhǎng)，短時(shí)間內(nèi)不會(huì)容易發(fā)生質(zhì)變；（8）加熱-控制系統(tǒng)很完整，具備智能溫度調(diào)節(jié)-控制器、排氣孔、感溫探頭、床體內(nèi)外固定結(jié)構(gòu)；（9）溫度調(diào)節(jié)方便，使用容易；（10）加熱系統(tǒng)用的材料是遠(yuǎn)紅外頻譜加熱板相當(dāng)?shù)陌踩?，用電、價(jià)格、加熱原理等方面更有優(yōu)勢(shì)、能耗低、傳熱均勻。

3 結(jié)束語

現(xiàn)如今這樣人類越來越注重自身健康問題的時(shí)候，沙療床尚未普遍的使用，在價(jià)格方面不是任何一個(gè)人能接受得了的價(jià)格。如果價(jià)格方面更有優(yōu)勢(shì)，用電情況更為節(jié)約型，治療效果最為好的治療結(jié)構(gòu)將來會(huì)是另一種需求。

綜上所述，根據(jù)對(duì)于沙療室的實(shí)際實(shí)驗(yàn)對(duì)比結(jié)果、沙療室與沙療床的數(shù)值計(jì)算結(jié)果，可得知沙療床傳熱效率高、能耗損失低、加熱系統(tǒng)安全，相比于沙療室沙療床是傳熱效率最佳，目前為最好的、最先進(jìn)的一個(gè)埋沙治療結(jié)構(gòu)。本文的研究?jī)?nèi)容在沙療床未來的應(yīng)用當(dāng)中定會(huì)有重要意義。

參考文獻(xiàn)：

[1]孔春元.沙療對(duì)關(guān)節(jié)炎動(dòng)物模型股骨強(qiáng)度及股動(dòng)脈血流影響的試驗(yàn)研究與有限元分析[D].新疆：新疆大學(xué)，2015.

[2]李艷娜.不同加熱方式對(duì)系統(tǒng)傳熱效率的影響[D].新疆：新疆大學(xué)，2017.

[3]Hsieh K J，Lien F S.Numerical modeling of buoyancy-driven turbulent flows in enclosures[J].International Journal of Heat and Fluid Flows，2004，25：659-670.

[4]居來提·買提肉孜，阿依先古力·熱合曼，等.不同加熱方式對(duì)傳熱效率的影響[J].機(jī)械設(shè)計(jì)與研究，2018（6）：36-40.

[5]張春廣.沙療室加熱過程的實(shí)驗(yàn)研究與數(shù)值模擬[D].新疆：新疆大學(xué)，2014.

[6]王福軍.計(jì)算流體動(dòng)力學(xué)分析：CFD軟件原理與分析[J].北京：清華大學(xué)出版社，2004.

[7]陶文銓.數(shù)值傳熱學(xué)（第2版）[M].西安：西安交通大學(xué)出版社，2001.

[8]趙穎杰.大型飛機(jī)座艙流場(chǎng)的數(shù)值計(jì)算[D].南京航空航天大學(xué)，2009.

[9]王雷，欽蘭云，佟明，等.快速成型制造臺(tái)階效應(yīng)及誤差評(píng)價(jià)方法[J].沈陽工業(yè)大學(xué)學(xué)報(bào)，2008，30（3）：82-85.

[10]陳軒，梁國(guó)柱.三維貼體曲線坐標(biāo)變換通用方法與程序開發(fā)[J].電子技術(shù)及信息科學(xué)，2013（11）：159-162.

[11]劉波，李忠媛，張濤.一種基于三角形非結(jié)構(gòu)化網(wǎng)格SIMPLE算法的程序設(shè)計(jì)[J].計(jì)算力學(xué)學(xué)報(bào)，2015（6）：111-115.