圓臺型螺旋能量樁換熱特性的模擬分析與驗證

劉亞姣，黃光勤※，盧 軍，楊小鳳，莊春龍

?

圓臺型螺旋能量樁換熱特性的模擬分析與驗證

劉亞姣1，黃光勤1※，盧 軍2，楊小鳳1，莊春龍1

（1. 陸軍勤務學院軍事設施系，重慶 401311； 2. 重慶大學城市建設與環境工程學院，重慶 400045）

為了減少地源熱泵的熱干擾、提高其換熱效率，該文提出了一種新型“圓臺型螺旋能量樁（truncated cone helix energy pile，CoHEP）”。同時為了更加準確地模擬新型圓臺型螺旋能量樁的換熱特性，綜合考慮初始土壤溫度的不均勻性和土壤上表面動態邊界條件，建立了三維數值模型，分析了不同錐角下圓臺型螺旋能量樁的換熱特性。結果表明：圓臺型螺旋能量樁沿水流方向可以分成4個換熱階段：入口階段→熱短路階段→小溫差階段→出口階段。圓臺型螺旋能量樁的底部熱干擾現象較為嚴重，且錐角越大，底部熱干擾越明顯。相反，由于頂部螺旋半徑較大，熱干擾效應較弱。圓臺型螺旋能量樁的單位管長換熱量隨著錐角的增加而線性增大，當系統運行時間為12 h時，錐角從0增加到10°再增加到20°，單位管長換熱量增長率分別為2.54%和3.53%。新型圓臺型螺旋能量樁單位管長換熱量大于傳統圓柱型螺旋能量樁，20°錐角條件下的單位管長換熱量比傳統圓柱型螺旋能量樁高了約6.16%。通過對新型圓臺型螺旋能量樁的換熱特性進行探析，相對于傳統圓柱型螺旋能量樁，提升了換熱效率，同時為后續的推廣應用和工程設計奠定了理論基礎。

熱泵；換熱；數值模擬；螺旋能量樁

0 引 言

螺旋能量樁是一種應用于地源熱泵系統的新型的地下熱交換器[1]。它通常埋設在建筑物的混凝土樁基礎上，與建筑結構相融合。與傳統的U型和W型能量樁相比[2-4]，螺旋能量樁在同一深度下具有較大的傳熱面積[5-7]。對于傳統的螺旋能量樁，換熱管纏繞在圓柱壁面上，在本文中被稱為“圓柱型螺旋能量樁（CyHEP，cylinder helix energy pile）”。

目前關于圓柱型螺旋能量樁的研究主要集中在進行一定程度簡化的解析解模型上。Man等[8]針對傳統圓柱型螺旋能量樁提出了無限長和有限長的圓柱熱源模型，并推導得出其溫度場表達式，但是該模型不能分析管壁溫度的變化規律。Zhang等[9]以圓柱熱源模型為基礎，提出了考慮不同螺距的線圈熱源模型，并推導計算出了線圈熱源模型解析解形式的表達式。該模型可以計算分析管壁的溫度變化。Wang等[10]考慮了能量樁內外巖土熱物性差異的影響，并在此基礎上提出了圓柱熱源模型和線圈

熱源模型的解析解表達式。Li等[11]綜合了圓柱熱源模型和線圈熱源模型，提出了考慮圓柱型螺旋能量樁復雜結構的三維解析解模型，使解析解模型更加貼合實際。通過分析可以發現，盡管解析解模型的求解計算較為便捷，但是由于各類假設條件的存在，導致與實際工況的差異較大。比如解析解模型并沒有考慮螺旋管內部水的流動，只是假定了螺旋管的單位管長發熱量[12-13]；同時也認為土壤的初始溫度是均勻一致的，這與實際情況不相符；另外，幾乎所有解析解模型都沒有考慮土壤上表面動態環境的影響。

除了解析解模型外，一些學者還研究了圓柱型螺旋能量樁的數值解模型。Bezyan等[14]建立了圓柱型螺旋能量樁的三維數值傳熱模型，并利用商業數值模擬軟件FLUENT 6.3.26計算分析了不同螺距下的換熱性能。但是該模型假設土壤初始溫度是均勻的，并且土壤上表面邊界溫度為恒定值。Jalaluddin等[15]利用計算流體力學軟件ANSYS FLUENT14.5對圓柱型螺旋能量樁進行三維數值模擬，并分析了其換熱和沿程壓降規律。Xiang等[16]針對圓柱型螺旋能量樁，通過對三維模型進行降維簡化，建立了新的數值解模型，該模型包括了管內流體域的一維瞬態子模型和管外固體域的一維瞬態子模型。通過分析發現，圓柱型螺旋能量樁的三維數值解模型較為復雜，研究人員總是做出一定的簡化處理，導致與實際情況有一定的偏差。例如，大多數數值解模型沒有考慮土壤初始溫度的不均勻性和土壤上表面的動態環境參數[17-20]；還有一些研究人員把復雜的三維模型簡化為一維或二維模型。

雖然與U型和W型能量樁相比，螺旋能量樁在同 一深度下具有較大的傳熱面積。但是仍然有相關研究發 現[21-23]，由于相鄰螺旋線圈之間的距離較小，傳統圓柱型螺旋能量樁的熱干擾作用較為劇烈，極大地限制了換熱效率的提升。但是目前研究人員主要是研究傳統圓柱型螺旋能量樁的數學模型，較少關注其熱干擾效應和沿程換熱特性。因此，為了減少熱干擾、提高換熱效率，本文提出了一種新型圓臺型螺旋能量樁，如圖1所示。圓臺型螺旋能量樁作為地源熱泵系統中最核心的地下換熱部件，可依托建筑基礎，埋設于建筑樁基內部，可較大程度節約鉆孔費用，減少系統初投資。但通過對實際工程調研分析發現，由于建筑樁基數量有限，現有實際工程多采用樁基埋管輔以室外鉆孔埋管才能滿足負荷需求，即在室外空地鉆孔埋設“圓臺型螺旋能量樁”，本文所研究圓臺型螺旋能量樁就是基于此背景工況開展的，需要考慮上表面動態環境參數的影響。

針對該新型圓臺型螺旋能量樁，課題組在前期進行了初步的研究[24]。主要是通過建立三維解析解模型來對其熱干擾和換熱特性進行研究。但是存在如下幾點假設：土壤初始溫度均勻；土壤上表面邊界條件為定值；螺旋管設定一個恒定的單位長度發熱量。顯然這些假設會導致與實際情況存在偏差。更重要的是，該解析解模型不能用于分析螺旋管的沿程換熱特性。

在前期工作的基礎上，為了更加準確地模擬新型圓臺型螺旋能量樁的換熱特性，本文綜合考慮初始土壤溫度的不均勻性和土壤上表面動態邊界條件，建立了三維數值模型。并進行了相關試驗，驗證模型的準確性。

1 圓臺型螺旋能量樁的數值解模型

該新型“圓臺型螺旋能量樁（CoHEP，truncated cone helix energy pile）”的幾何特點在前期成果中已做了相關介紹[24]。前期研究發現，圓臺型螺旋能量樁采用上下顛倒埋設的2種方式，其單位管長換熱量相差不大。同時考慮到挖孔鉆井施工的可行性及安全性，采用上大下小的埋設方式，一方面便于施工開挖，另一方面可以使井壁具有一定的傾斜角度，降低了開挖過程中井壁坍塌的風險，所以本文提出的圓臺型螺旋能量樁均采用上大下小的埋設方式，具體結構如圖1所示，螺旋管纏繞在圓錐壁面上，且螺旋半徑隨著埋深的增加而線性減小，螺旋線遵循以下方程：

由于圓臺型螺旋能量樁的幾何結構復雜，實際尺寸的三維模擬計算工程量龐大，對硬件要求較高。為了提高模擬計算的效率，本文基于相似原理建立了圓臺型螺旋能量樁的縮比尺寸模型。這將有助于實現圓臺型螺旋能量樁的高效長期模擬。同時，還模擬計算了螺旋管沿程水溫變化，進而分析其沿程換熱特性。

注：hb為地表面距能量樁底面的距離，m；ht為地表面距能量樁頂面的距離，m；rt為能量樁頂面螺旋半徑，m；rb為能量樁底面螺旋半徑，m；b為深度方向的螺距，m；d為軸向相鄰管段間的距離，m；θ為錐角，(°)。

1.1 縮比模型

本文采用相似性原理對原型尺寸模型進行縮小。由相似性原理可知，為了使原型和縮比模型的熱性能相似，2個模型的相似性指標必須等于1。因此在設計縮比模型時，在原型基礎上采用長度尺度系數對其進行類似換算，土壤和流體的熱物理參數與原型相同。根據傳熱相似理論，由于縮比模型是按照同一長度尺度系數對原型進行等比例縮小，且原型和縮比模型的土壤和流體的熱物性參數均相同，所以經過計算可知原型和縮比模型中的Nu、Fo和Re數是相等的。余延順等針對地源熱泵系統原型和縮比模型之間的關系進行探究，得到二者之間的時間關系是[25]：

因此，通過對原型尺寸進行縮小，與原型模擬的長時間計算周期相比，可以在相對較短的時間內獲得計算結果，大大減少對計算機配置的要求。通過以上分析，根據相關標準規范及工程實際確定原型尺寸[26]，進而確定縮比模型相應的主要參數，如表1所示。

在本文中，模型均由軟件Solidworks構建，并導入軟件Workbench中劃分網格。由于螺旋管周圍的溫度變化較大，需要更高的計算精度，管周圍的網格進行了部分加密處理。同時，隨著半徑的增大，土壤溫度變化越來越小，所以網格大小沿徑向逐漸增大。同時進行網格無關性驗證，針對縮比模型選取4種不同的網格尺度大小分別進行網格劃分，得到4種不同的網格數量（253萬、350萬、397萬和434萬）如圖2所示，從圖中可以看出，當網格數量從253萬逐步增大到434萬時，分別模擬運行12 h后出口溫度值的變化不大，說明選取的4種網格數量對計算結果的影響很小，可認為253萬網格數量已經達到網格無關，因此后續模擬選取253萬的網格數量進行模擬計算。

表1 原型和縮比模型的主要參數

圖2 網格無關性驗證

模擬的主要參數如表2所示，各參數均是依據縮比試驗平臺測得，與驗證試驗一致。其中導熱系數通過“球體法”測得，比熱容通過“混合量測法”測得，螺旋管的相關參數由生產廠商提供[27]。

表2 數值模擬的主要參數

1.2 控制方程

對于地埋管換熱器，管內流體與管壁之間屬于強迫對流換熱，流體處于湍流狀態。本文中假設土壤為均勻各向同性的固體，同時不考慮地下水滲流的影響。

由于土壤被視為均質各向同性的固體，且不考慮地下水滲流的影響，所以土壤中只存在導熱作用。土壤中的熱控制方程如式（4）所示：

式中，ρ為土壤密度，kg/m3；c為土壤比熱容，J/(kg·K)；λ為土壤導熱系數，W/(m·K)；T為土壤溫度，℃。

1.3 初始及邊界條件

由于螺旋能量樁埋深較淺，通常埋設于地表以下10～20 m處，因此地埋管的傳熱受土壤上表面環境條件的影響較大。特別地，由于新型圓臺型螺旋能量樁的頂部螺旋直徑比傳統的圓柱型螺旋能量樁大，靠近土壤上表面的換熱管面積更大，因此，外界動態邊界條件對新型圓臺型螺旋能量樁的換熱影響不容忽視。在開展的數值模擬驗證試驗中，記錄了試驗運行12 h (720 min)期間土壤上表面處的動態環境溫度，并通過數據擬合得到土壤上表面動態環境溫度的函數表達式為：

初始土壤溫度受外界氣溫和太陽輻射的共同作用，因此地溫與深度和時間存在相應的關系，嚴重影響著地埋管的換熱性能，是理論計算中最重要的參數之一。特別是對于螺旋能量樁，由于其埋深較淺，初始土壤溫度分布不均勻，考慮其深度方向的溫度不均勻性將使模型更加貼合實際情況。在開展的數值模擬的驗證試驗中，記錄了初始時刻在圓柱形土壤箱（直徑1.8 m×高1.5 m）深度方向上的土壤溫度分布數據，并通過擬合得到了初始時刻土壤溫度的函數表達式如下：

式中T為初始時刻（2017年9月21日9:08）土壤溫度，℃；為地下深度，m。

根據上述擬合函數方程，利用C語言編寫自定義函數程序（UDF，user defined function）對土壤上表面動態環境溫度和初始土壤溫度進行定義，并編譯納入數值計算模型。對于模型的其它邊界條件，將圓柱形土壤模型的頂面、底面和四周定義為wall，螺旋管入口定義為Velocity-Flow-Inlet，出口定義為Pressure-Outlet，螺旋管和土壤定義為Solid，螺旋管內部定義為Fluid。

1.4 FLUENT軟件設置

本文使用商用CFD模擬軟件ANSYS_FLUENT進行數值模擬計算，由于螺旋管長20 m，遠遠大于螺旋管內徑0.005 m，屬于細長區域傳熱模型，選擇3D雙精度求解器進行計算。根據驗證試驗確定管內流體、螺旋管和土壤的熱物性參數（見表2）。依據實際工程情況將埋設有圓臺型螺旋能量樁的圓柱型土壤箱（見圖3a）的頂面定義為傳熱邊界，底面和四周定義為絕熱邊界條件以避免外界的影響。求解采用SIMPLEC算法和二階迎風格式進行計算。

2 數值模擬的試驗驗證

為了驗證第1節所建立數值模擬模型的準確性開展了驗證試驗，試驗系統如圖3a所示，包括埋有圓臺型螺旋能量樁的圓柱形土壤箱（直徑1.8 m×高1.5 m）、流量計（型號：LZB-6，量程：6～60 L/h，精度：2.5）、水泵（型號：MP-6R）、恒溫水箱、熱電偶（型號：K型）和數據采集儀（型號：SH-X）等。系統中的水在恒溫水箱中被持續加熱，然后通過水泵輸送到圓臺型螺旋能量樁，在能量樁中熱量由熱水散入到低溫土壤中，數據采集儀記錄出口水溫并監測其變化情況。測試從2017年9月21日9:08開始，運行12 h。有關參數設定如表1和表2所示。此外，圓柱形土壤箱的外測和外部管道均包裹著保溫材料，以減少外界環境的干擾[28-30]。

為了驗證數值模型的準確性，將通過數值模擬計算得到的出口水溫與數據采集儀記錄的實際出口水溫進行比較，如圖3b所示。從圖中可以看到，出口水溫隨時間逐漸升高并最終趨于穩定。模擬值與實測值吻合較好，兩者之間的最大差異出現在160 min時刻，經計算差異為1.5%，在數值模擬驗證的允許誤差范圍5%以內[31]。由以上分析可知，模擬結果與試驗測試結果吻合較好，數值模擬計算的準確性得到驗證。

3 結果與分析

3.1 換熱能力分析

圖3 數值模擬的試驗驗證

圖4所示為不同錐角條件下圓臺型螺旋能量樁換熱能力隨時間的變化曲線。圖4a顯示了3種錐角條件下模擬出口水溫隨運行時間的變化情況，可以發現出口水溫隨時間逐漸增大并最終達到穩定狀態。相應的，圖4b顯示了3種錐角條件下單位管長換熱量隨運行時間的變化情況，結果表明，單位管長換熱量隨著時間的推移而逐漸減小，最終達到穩定狀態。分析造成這一趨勢的主要原因為：系統運行初期，管內流體與土壤的溫差最大，換熱較為充分，因此在系統運行初期單位管長換熱量較大。然而隨著系統的持續運行，能量樁的熱影響范圍逐漸擴大，土壤的熱阻成為傳熱過程的主要矛盾，導致傳熱效率降低。經過一段時間的換熱后，土壤中的熱量持續積聚，管內流體與土壤的溫差持續減小，導致單位管長換熱量持續降低。之后當土壤溫度場趨于穩定時，換熱性能也逐漸穩定。

表3 不同錐角圓臺型螺旋能量樁尺寸參數

圖4 圓臺型螺旋能量樁換熱能力隨時間的變化曲線

如圖5所示，單位管長換熱量隨著錐角的增大而線性增大。當系統運行時間分別為4、8、12 h時，錐角為0°時，單位管長換熱量在不同時刻分別為4.46、3.34、2.76 W/m。錐角為10°時，單位管長換熱量在不同時刻分別為4.52、3.42、2.83 W/m。在20°錐角情況下，單位管長換熱量在不同時刻分別為4.74、3.52、2.93 W/m。隨著錐角的增大，圓臺型螺旋能量樁單位管長換熱量也就相應增大。當系統運行時間為12 h時，錐角從0°增加到10°再增加到20°，單位管長換熱量增長率分別為2.54%和3.53%。分析可知，隨著圓臺型螺旋能量樁錐角的增大，其單位管長換熱量也隨之增加，且呈線性增長趨勢?？梢娫谝欢ǚ秶鷥仍龃髨A臺型螺旋能量樁的錐角可增加單位管長換熱量，提高其換熱能力。

圖5 單位管長換熱量與錐角之間的關系

3.2 沿程換熱特性分析

圖6所示為不同錐角（0°、10°、20°）條件下單位管長換熱量沿管長方向的變化情況。根據圓臺型螺旋能量樁的幾何特性，結合圖6所示的沿程換熱特性，可以將螺旋管沿水流方向分為4個不同且連續的換熱階段：

1）入口階段S1（約0～2.5 m）：該階段位于螺旋能量樁的頂部，與上層的覆土區域之間接觸，當高溫流體經過這一階段時，熱量除了沿經向傳遞外，還有一部分沿軸向傳入無熱干擾存在的覆土區域。同時，由于圓臺型能量樁頂部的螺旋半徑最大，半徑方向的熱干擾最小，熱堆積最小，所以該階段換熱能力最大。

2）熱短路階段S2（約2.5～10 m）：隨著流體流過入口階段，由于土壤的熱容量有限，熱量積聚在土壤中無法快速散去，螺旋管之間沿軸向存在熱短路現象，導致熱量主要沿經向傳遞，所以該階段的換熱性能有所下降。

3）小溫差階段S3（約10～17.5 m）：當熱短路階段結束后，管道內的流體溫度持續的減小，導致流體與周圍土壤間的換熱溫差減小，換熱效率降低。同時，由于圓臺型螺旋能量樁下部螺旋半徑較小，由于熱量的積聚，在經向和軸向均存在嚴重的熱干擾效應。因此在該階段換熱能力進一步下降。

4）出口階段S4（約17.5～20 m）：和入口階段相類似，該階段的螺旋管與能量樁下部區域的土壤相接觸，軸向熱干擾相對較小，有一部分熱量沿軸向傳入能量樁下方土壤。因此在該階段換熱能力稍有提升。

從圖6可以看出，在入口階段和熱短路階段，0°錐角下螺旋管的單位管長換熱量最小，20°錐角下螺旋管的單位管長換熱量最大。這是因為20°錐角條件下圓臺型螺旋能量樁的頂部螺旋半徑最大，徑向熱干擾較小，換熱能力較大。相反，在小溫差階段和出口階段，0°錐角下的單位管長換熱量最大，20°錐角下的單位管長換熱量最小。原因在于圓臺型螺旋能量樁在20°錐角條件下的底部螺旋半徑最小，經向的熱干擾較大，導致傳熱效率較低。

注：S1—入口階段；S2—熱短路階段；S3—小溫差階段；S4—出口階段。

綜上所述，新型圓臺型螺旋能量樁的頂部螺旋半徑大于傳統圓柱型螺旋能量樁，導致徑向熱干擾小，換熱能力強。雖然圓臺型螺旋能量樁底部的熱干擾現象嚴重，換熱能力低，但入口階段和熱短路階段的流體溫度（換熱溫差）高于后2個階段，因此前個階段的換熱能力遠遠大于后個階段的換熱能力，承擔了整個螺旋能量樁的主要換熱量。所以整體而言，錐角越大，換熱性能越優越。而且，在深度方向的螺距（）相同的條件下，圓臺型螺旋能量樁軸向相鄰管段之間的距離（）明顯大于傳統圓柱型螺旋能量樁，且錐角越大，越大，可以有效地減少軸向熱干擾，提高換熱效率。因此，新型圓臺型螺旋能量樁的單位管長換熱量大于傳統圓柱型螺旋能量樁，且20°錐角條件下單位管長換熱量最大，比0°條件下（傳統圓柱型螺旋能量樁）高了6.16％。

3.3 熱響應分析

圖7所示為不同錐角條件下的管壁溫度分布云圖。從圖中可以看出，與0°錐角相比，20°錐角條件下的高溫流體幾乎全部分布在螺旋半徑較大的能量樁上部，由于上部區域的熱干擾效應相對較小，同時高溫流體可以直接與上方覆土區域進行熱量交換，提升了能量樁整體的換熱性能。相比之下，0°錐角下的高溫流體流過更多的螺旋圈數，幾乎到達了能量樁的中部，導致換熱效果較差。

圖8所示為不同錐角條件下的縱向對稱面的溫度分布云圖。從圖中可以看出，圓臺型螺旋能量樁的底部熱干擾現象較為嚴重，且錐角越大，底部熱干擾越明顯。這是因為圓臺型螺旋能量樁底部的螺旋半徑較小，由于熱量的積聚，在經向和軸向都存在熱短路現象。相反，由于頂部螺旋半徑較大，熱干擾效應較弱。另外從圖中可以明顯看出，隨著錐角的增大，軸向相鄰管段間的土壤高溫區域明顯減小，這是因為在相同深度方向螺距（）的條件下，錐角越大，軸向相鄰管段間的距離（）越大，可以有效地減少軸向熱干擾。

圖7 管壁溫度分布

圖8 縱向對稱面的溫度分布

4 結 論

為了減少熱干擾、提高換熱效率，提出了一種新型“圓臺型螺旋能量樁（truncated cone helix energy pile，CoHEP）”。同時為了更加準確地模擬研究新型圓臺型螺旋能量樁的換熱特性，本文綜合考慮初始土壤溫度的不均勻性和土壤上表面動態邊界條件，利用數值模擬的方法分析研究了不同錐角下圓臺型螺旋能量樁的換熱特性。通過分析得到如下結論：

1）圓臺型螺旋能量樁沿水流方向可以分為四個換熱階段：入口階段→熱短路階段→小溫差階段→出口階段。入口階段的換熱能力最強，在熱短路階段和小溫差階段換熱能力持續下降，出口階段有所回升。

2）圓臺型螺旋能量樁的底部熱干擾現象較為嚴重，且錐角越大，底部熱干擾越明顯。相反，由于頂部螺旋半徑較大，熱干擾效應較弱。此外，錐角越大，在相同深度方向螺距（）的條件下，軸向相鄰管段間的距離（）越大，軸向熱干擾作用越小。

3）在相同的管長（L）、深度方向螺距（）和能量樁高度（h）的情況下，增大錐角將增大能量樁頂部螺旋半徑（r），導致能量樁上部熱干擾作用減小，同時使更多的高溫流體位于能量樁上部，和覆土區域直接接觸，換熱能力增大。此外，隨著錐角的增大，軸向相鄰管段間的距離增大，軸向熱干擾作用減小。圓臺型螺旋能量樁的單位管長換熱量隨著錐角的增大而線性增加，當系統運行時間為12 h時，錐角從0°增加到10°再增加到20°，單位管長換熱量增長率分別為2.54%和3.53%。

4）新型圓臺型螺旋能量樁上部的熱干擾遠小于傳統圓柱型螺旋能量樁，且有更多的高溫流體位于能量樁上部。此外，在深度方向螺距（）相同的條件下，新型圓臺型螺旋能量樁軸向相鄰管段間的距離（）明顯大于傳統圓柱型螺旋能量樁，可有效降低軸向熱干擾。因此，新型圓臺型螺旋能量樁單位管長換熱量大于傳統圓柱型螺旋能量樁。且20°錐角條件下的單位管長換熱量比傳統圓柱型螺旋能量樁高了6.16％。

[1] Hamada Yasuhiro, Saitoh Hisashi, Nakamura Makoto, et al. Field performance of an energy pile system for space heating[J]. Energy and Buildings, 2007, 39(5): 517－524.

[2] Cecinato Francesco, Loveridge Fleur A. Influences on the thermal efficiency of energy piles[J]. Energy, 2015, 82: 1021－1033.

[3] Gao Jun, Zhang Xu, Liu Jun, et al. Thermal performance and ground temperature of vertical pile-foundation heat exchangers: A case study[J]. Applied Thermal Engineering, 2008, 28 (17-18): 2295－2304.

[4] Park Hyunku, Lee Seung-Rae, Yoon Seok, et al. Evaluation of thermal response and performance of PHC energy pile: Field experiments and numerical simulation[J]. Applied Energy, 2013, 103: 12－24.

[5] Zarrella Angelo, Capozza Antonio, De Carli Michele. Analysis of short helical and double U-tube borehole heat exchangers: A simulation-based comparison[J]. Applied Energy, 2013, 112: 358－370.

[6] Zarrella Angelo, De Carli Michele. Heat transfer analysis of short helical borehole heat exchangers[J]. Applied Energy, 2013, 102: 1477－1491.

[7] Zarrella Angelo, De Carli Michele, Galgaro Antonio. Thermal performance of two types of energy foundation pile: Helical pipe and triple U-tube[J]. Applied Thermal Engineering, 2013, 61(2): 301－310.

[8] Man Yi, Yang Hongxing, Diao Nairen, et al. A new model and analytical solutions for borehole and pile ground heat exchangers[J]. International Journal of Heat and Mass Transfer, 2010, 53(13-14): 2593－2601.

[9] Zhang Wenke, Yang Hongxing, Lu Lin, et al. Investigation on heat transfer around buried coils of pile foundation heat exchangers for ground-coupled heat pump applications[J]. International Journal of Heat and Mass Transfer, 2012, 55 (21): 6023－6031.

[10] Wang Deqi, Lu Lin, Cui Ping. A novel composite-medium solution for pile geothermal heat exchangers with spiral coils[J]. International Journal of Heat and Mass Transfer, 2016, 93: 760－769.

[11] Li Min, Lai Alvin C. K. Heat-source solutions to heat conduction in anisotropic media with application to pile and borehole ground heat exchangers[J]. Applied Energy, 2012, 96: 451－458.

[12] Zhang Wenke, Yang Hongxing, Lu Lin, et al. The research on ring-coil heat transfer models of pile foundation ground heat exchangers in the case of groundwater seepage[J]. Energy and Buildings, 2014, 71: 115－128.

[13] Zhang Wenke, Yang Hongxing, Lu Lin, et al. The analysis on solid cylindrical heat source model of foundation pile ground heat exchangers with groundwater flow[J]. Energy, 2013, 55: 417－425.

[14] Bezyan Behrad, Porkhial Soheil, Mehrizi Abbasali Aboui. 3-D simulation of heat transfer rate in geothermal pile- foundation heat exchangers with spiral pipe configuration[J]. Applied Thermal Engineering, 2015, 87: 655－668.

[15] Jalaluddin, Akio Miyara. Thermal performance and pressure drop of spiral-tube ground heat exchangers for ground-source heat pump[J]. Applied Thermal Engineering, 2015, 90: 630－637.

[16] Xiang Yong, Su Hua, Gou Wenshi, et al. A new practical numerical model for the energy pile with spiral coils[J]. International Journal of Heat and Mass Transfer, 2015, 91: 777－784.

[17] Ghasemi-Fare Omid, Basu Prasenjit. A practical heat transfer model for geothermal piles[J]. Energy and Buildings, 2013, 66: 470－479.

[18] Li Huai, Nagano Katsunori, Lai Yuanxiang. A new model and solutions for a spiral heat exchanger and its experimental validation[J]. International Journal of Heat and Mass Transfer, 2012, 55(15-16): 4404－4414.

[19] Loveridge Fleur, Powrie William. Temperature response functions (G-functions) for single pile heat exchangers[J]. Energy, 2013, 57: 554－564.

[20] Loveridge Fleur, Powrie William. G-Functions for multiple interacting pile heat exchangers[J]. Energy, 2014, 64: 747－757.

[21] Park Sangwoo, Lee Dongseop, Choi Hyun-Jun, et al. Relative constructability and thermal performance of cast-in-place concrete energy pile: Coil-type GHEX (ground heat exchanger) [J]. Energy, 2015, 81: 56－66.

[22] Park Sangwoo, Sung Chihun, Jung Kyoungsik, et al. Constructability and heat exchange efficiency of large diameter cast-in-place energy piles with various configurations of heat exchange pipe[J]. Applied Thermal Engineering, 2015, 90: 1061－1071.

[23] Yang Weibo, Lu Pengfei, Chen Yongping. Laboratory investigations of the thermal performance of an energy pile with spiral coil ground heat exchanger[J]. Energy and Buildings, 2016, 128: 491－502.

[24] Huang Guangqin, Yang Xiaofeng, Liu Yajiao, et al. A novel truncated cone helix energy pile: Modelling and investigations of thermal performance[J]. Energy and Buildings, 2018, 158: 1241－1256.

[25] 余延順，張少凡，馬娟，等. 土壤耦合熱泵系統模型試驗臺設計[J]. 南京理工大學學報(自然科學版)，2010，34(5)：613－617. Yu Yanshun, Zhang Shaofan, Ma Juan, et al. Design of scale-model experiment platform for ground-coupled heat pump system[J]. Journal of Nanjing University of Science and Technology (Natural Science), 2010, 34(5): 613－617. (in Chinese with English abstract)

[26] CECS 344: 2013，地源熱泵系統地埋管換熱器施工技術規程[S].

[27] 孔磊. 螺旋型地埋管換熱器換熱性能的研究[D]. 揚州：揚州大學，2015. Kong Lei. Study on Heat Transfer Performance of Helix Buried Heat Exchanger[D]. Yangzhou: Yangzhou University, 2015.

[28] Lim Kyoungbin, Lee Sanghoon, Lee Changhee. An experimental study on the thermal performance of ground heat exchanger[J]. Experimental Thermal and Fluid Science, 2007, 31(8): 985－990.

[29] Luo Jin, Rohn Joachim, Xiang Wei, et al. Experimental investigation of a borehole field by enhanced geothermal response test and numerical analysis of performance of the borehole heat exchangers[J]. Energy, 2015, 84: 473－484.

[30] Vashistha Chaitanya, Patil Anil Kumar, Kumar Manoj. Experimental investigation of heat transfer and pressure drop in a circular tube with multiple inserts[J]. Applied Thermal Engineering, 2016, 96: 117－129.

[31] William L. Oberkampf, Timothy G. Trucano. Verification and validation in computational fluid dynamics[J]. Progress in Aerospace Sciences, 2002, 38(3): 209—272.

[32] CJJ-101, 埋地聚乙烯給水管道工程技術規程[S].

Simulation analysis and verification of heat transfer characteristics of truncated cone helix energy pile

Liu Yajiao1, Huang Guangqin1※, Lu Jun2, Yang Xiaofeng1, Zhuang Chunlong1

(1.401311,; 2.400045,)

In order to reduce the thermal interference and to improve the heat transfer efficiency, a novel “truncated cone helix energy pile (CoHEP)” was proposed in this study. And in order to simulate the thermal performance of the novel CoHEP more accurately, a three-dimensional numerical model was established with consideration of the inhomogeneous initial soil temperature and the dynamic boundary condition above the soil surface. The heat transfer characteristic of the novel CoHEP of different cone angles was studied. In addition, in order to verify the accuracy of the numerical model, a verification test was conducted. The results showed that the novel CoHEP can be divided into four heat transfer stages along the flow direction: entrance stage → thermal short circuit stage → small temperature difference stage → exit stage. Among them, the heat transfer capacity of the entrance stage was the strongest, and the heat transfer capacity continued to decrease during the thermal short circuit stage and the small temperature difference stage, eventually the heat transfer capacity rebounded during the exit stage. The thermal interference at the bottom of the CoHEP was more serious, and the larger the cone angle, the more serious the thermal interference at the bottom. In contrast, the thermal interference effect was weaker at the top of the CoHEP due to the larger helix radius at the top. The heat flux per unit pipe length of the CoHEP increased linearly with the increase of cone angle. That was because with the same pipe length (L), pitch in the depth direction () and energy pile height (h), increasing the cone angle would increase the top radius (r) of the CoHEP, leading to small thermal interference in the upper part of the CoHEP. At the same time, more high-temperature fluid was located in the upper part of the CoHEP which directly contacted with the covered soil area. Thus the heat transfer capacity increased. When the system operating time was 12 h, the cone angle increased from 0°to 10°to 20°, the increasing rate of the heat flux was 2.54% and 3.53%, respectively. The thermal interference in the upper part of the novel CoHEP was much smaller than that of the traditional CyHEP, and more high-temperature fluid was located in the upper part of the energy pile which was good for heat transfer. In addition, the distance between the adjacent pipes in the axial direction () of the novel CoHEP was significantly larger than that of the traditional CyHEP under the same pitch in the depth direction (), which can effectively reduce the axial thermal interference. Thus the heat flux per unit pipe length of the novel CoHEP was greater than the traditional CyHEP. And when the cone angle was 20°, the heat flux per unit pipe length of the novel CoHEP was 6.16% higher than that of the traditional CyHEP.

heat pump system; heat transfer; numerical simulation; helix energy pile

劉亞姣，黃光勤，盧 軍，楊小鳳，莊春龍.圓臺型螺旋能量樁換熱特性的模擬分析與驗證 [J]. 農業工程學報，2018，34(7)：227－234. doi：10.11975/j.issn.1002-6819.2018.07.029 http://www.tcsae.org

Liu Yajiao, Huang Guangqin, Lu Jun, Yang Xiaofeng, Zhuang Chunlong. Simulation analysis and verification of heat transfer characteristics of truncated cone helix energy pile[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(7): 227－234. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2018.07.029 http://www.tcsae.org

2017-12-05

2018-02-26

國家自然科學基金項目（51706243）；國家自然科學基金項目（51478058）；重慶市基礎科學與前沿技術研究項目（cstc2016jcyjA0496）

劉亞姣，男，河南舞陽人，博士，主要從事建筑節能與淺層地熱能利用方面的研究。Email：lyj199164@126.com

黃光勤，男，重慶人，博士、講師，主要從事建筑節能與淺層地熱能利用方面的研究。Email：hgq880818@163.com

10.11975/j.issn.1002-6819.2018.07.029

TK124

A

1002-6819(2018)-07-0227-08