單側(cè)加熱方形通道內(nèi)超臨界水傳熱研究

引 言

再生冷卻一般利用燃料的熱沉冷卻燃燒室壁面，然而隨著飛行速度的增加，高超聲速飛行器燃燒室壁面的散熱面臨極大挑戰(zhàn)。常規(guī)的解決方法是通過熱管理提高燃料熱沉的利用效率[1]，然而僅依靠提高燃料熱沉的利用效率無法滿足冷卻要求，需要增加冷卻劑以輔助燃料冷卻燃燒室。常用冷卻劑有正癸烷、航空煤油RP-3、水和其他高密度烴。與其他冷卻劑相比，水具有成本低、來源廣和換熱能力強的特點，作為優(yōu)質(zhì)冷卻劑廣泛應(yīng)用于各種換熱器。由于燃燒室壁溫高，常壓下的水容易出現(xiàn)相變，而超臨界水的密度近似液體的密度量級，同時又具有與氣體相近的擴散能力，具有非常好的流動性和傳輸性能，超臨界水十分適合作為飛行器燃燒室壁面冷卻劑。

超臨界水在核燃料堆冷卻換熱研究較為深入，Shang等[2]數(shù)值模擬研究了直徑對超臨界水在水平圓管中流動的影響，發(fā)現(xiàn)直徑影響傳熱是由于邊界層浮力的變化，小管徑更有利于傳熱。范辰浩等[3]實驗研究了小管徑內(nèi)超臨界水的傳熱惡化特性，發(fā)現(xiàn)在高熱通量時傳熱惡化的產(chǎn)生是由于流動層流化的發(fā)展，而管內(nèi)徑向的工質(zhì)溫度和物性的巨大差異促進了流動層流化的發(fā)展。Wang等[4]基于大渦模擬研究了圓形管道中超臨界水在無重力、豎直向上和豎直向下流動下的傳熱，發(fā)現(xiàn)與無重力流動相比，向下流動的傳熱略有改善，向上流動的傳熱出現(xiàn)惡化，提出由于浮力的影響而產(chǎn)生這些差異。

Bai等[5]數(shù)值模擬了豎直向上圓形管道中超臨界水非均勻加熱對流動換熱的影響，發(fā)現(xiàn)與均勻加熱方式相比，非均勻加熱方式表現(xiàn)出顯著的差異，截面溫度分布沿周向呈現(xiàn)出較大的不均勻性，傳熱強化僅出現(xiàn)在局部區(qū)域。Gao等[6]對超臨界水在非均勻加熱的水平圓管中的流動換熱進行了數(shù)值模擬，分析了超臨界水的非均勻傳熱特性及機理，并討論了側(cè)向二次流和浮升力效應(yīng)對傳熱作用開始的判據(jù)，發(fā)現(xiàn)二次流與熱物性之間的相互影響對傳熱具有重要影響。許多文獻對超臨界水傳熱強化和惡化產(chǎn)生的機理進行了分析[7-11]，發(fā)現(xiàn)流體物性變化對超臨界水傳熱具有重要影響。文獻中對超臨界水換熱關(guān)聯(lián)式的適用性也進行了研究[12-17]，發(fā)現(xiàn)超臨界水的關(guān)聯(lián)式需要根據(jù)不同工況的換熱特點進行修正。

現(xiàn)有研究主要關(guān)注超臨界水在圓形通道內(nèi)的流動傳熱，非圓結(jié)構(gòu)如類矩形、類三角形和方環(huán)形等結(jié)構(gòu)也存在相關(guān)研究[18-21]，而高超聲速飛行器燃燒室壁面的冷卻通道具有方形單側(cè)加熱特點，目前缺乏超臨界水在方形單側(cè)受熱通道內(nèi)流動傳熱的相關(guān)研究，亟需獲得其流動傳熱規(guī)律。

本文對超臨界水在方形冷卻通道的流動換熱進行數(shù)值模擬研究，通過分析流場和超臨界水熱物性變化，闡明傳熱惡化產(chǎn)生和恢復(fù)的機理，并驗證超臨界流體常用傳熱關(guān)聯(lián)式在單側(cè)加熱方形通道內(nèi)的適用性，最后，提出非對稱凹槽和雙通道強化結(jié)構(gòu)，防止通道傳熱惡化，提升通道綜合換熱性能并進行場協(xié)同原理分析。

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

1.1 計算模型

航空發(fā)動機冷卻通道緊貼高溫壁面，冷卻介質(zhì)在多個微小通道內(nèi)流動，考慮模型的周期性，研究的模型可以簡化為單側(cè)加熱的方通道，如圖1所示。方通道當(dāng)量直徑為D，壁厚a=2 mm，方形通道加熱長度為Lt，為得到較均勻的速度分布并防止回流，模擬過程中增加長度Ld=Lu=0.1 m的入口段和出口段。

圖1

圖1   方形通道示意圖

Fig.1   Schematic of square channel

為提升通道換熱性能，提出了多種強化通道。圖2（a）~（c）分別表示對稱凹槽、倒角凹槽和非對稱倒角凹槽的結(jié)構(gòu)。其中加熱長度Lt=0.5 m，凹槽半徑R=2 mm，槽深d=1 mm，凹槽間距為25 mm，凹槽數(shù)目為16?？紤]入口處換熱較強，無須強化，第一個凹槽布置在距離入口50 mm處，其余凹槽均勻分布于內(nèi)壁面靠近加熱面處。凹槽加倒角結(jié)構(gòu)如圖2（b）、（c）所示，其中圓角半徑r=1 mm，r1=2 mm，r2=0.6 mm。非對稱倒角結(jié)構(gòu)由曲線與6個點（P0，P1，P2，P3，P4，P5）確定。此外，也研究了非對稱凹槽最低點偏下游以及偏上游對通道流動換熱的影響。

圖2

圖2   凹槽結(jié)構(gòu)示意圖

Fig.2   Schematic of grooves

雙通道結(jié)構(gòu)的布置是以1 mm厚度的壁面平均分割流體域，考慮兩種布置形式，即上下布置（下層靠近加熱端）與并排布置。二者均為對稱結(jié)構(gòu)，且結(jié)構(gòu)參數(shù)一致。

1.2 邊界條件與求解方法

采用質(zhì)量流量入口、壓力出口與均勻熱流邊界條件，計算工況如表1所示。基于商業(yè)軟件Fluent求解，速度與壓力的耦合采用SIMPLE算法求解，控制方程中對流項的離散采用二階迎風(fēng)格式，擴散項采用QUICK格式。當(dāng)連續(xù)性方程、動量方程和湍流方程殘差值小于1.0×10-5，能量方程的殘差值小于1.0×10-8，且不隨計算步數(shù)變化時認(rèn)為結(jié)果收斂。超臨界水熱物性通過調(diào)用NIST數(shù)據(jù)庫獲得。

表1   計算工況參數(shù)

Table 1  Conditions of numerical simulation

1.3 模型驗證

在靠近固體和流體區(qū)域的內(nèi)壁附近采用精細(xì)網(wǎng)格，確保無量綱距離y+小于1，以滿足SST k-ω模型的要求。經(jīng)網(wǎng)格的獨立性考核，最終采用627萬的網(wǎng)格數(shù)量進行計算。

由于目前缺乏對方形通道中超臨界水傳熱的實驗研究，首先對超臨界水在圓形管道中豎直向上流動的數(shù)值模擬結(jié)果與Yamagata等[22]的實驗結(jié)果進行對比，實驗段由AISI 316型不銹鋼管光管制成，內(nèi)徑為7.5 mm，壁厚為2 mm。實驗工況為：壓力P=24.5 MPa，質(zhì)量流速G=1260 kg·m-2·s-1，熱通量q=233 kW·m-2。如圖3所示，計算結(jié)果與文獻中的數(shù)據(jù)吻合很好，最大誤差為0.59%，表明本計算模型能夠較好地處理超臨界水的流動傳熱現(xiàn)象。

圖3

圖3   實驗結(jié)果與數(shù)值模擬結(jié)果的比較

Fig.3   Comparison between available experimental and present numerical results

2 結(jié)果分析與討論

2.1 傳熱惡化機理分析

本節(jié)對超臨界水在豎直向上流動時出現(xiàn)傳熱惡化的工況（case 1）進行了分析。由于研究對象為超臨界水的豎直向上流動，而此時流動方向與重力方向相反，二次流僅僅依靠流體密度差產(chǎn)生，影響十分微小，因此不考慮二次流動的影響。

Koshizuka等[23]對傳熱惡化時的表面?zhèn)鳠嵯禂?shù)做出如下定義:

h/hDB<0.3$h/h_{DB}<0.3$(1)

其中hDB為Dittus-Boelter關(guān)聯(lián)式[24]計算的表面?zhèn)鳠嵯禂?shù)。

圖4表示了數(shù)值模擬與Dittus-Boelter關(guān)聯(lián)式計算的內(nèi)壁面局部表面?zhèn)鳠嵯禂?shù)的比值隨流體平均焓的變化，可以發(fā)現(xiàn)傳熱惡化的起始點A、傳熱惡化點B和傳熱惡化恢復(fù)點C，這是后續(xù)分析的三個基準(zhǔn)點。

圖4

圖4   豎直向上與無重力流動內(nèi)壁面h/hDB沿程變化

Fig.4   h/hDB on the inner wall with and without gravity

需要注意的是，當(dāng)流體的焓值過高時，會再次出現(xiàn)h/hDB < 0.3的現(xiàn)象，這是由于流體溫度靠近擬臨界溫度，比定壓熱容增加劇烈，Prandtl數(shù)過大，導(dǎo)致hDB偏高，因此前面ABC三點才是傳熱惡化產(chǎn)生的區(qū)域。

圖5、圖6展示了傳熱惡化區(qū)域A、B和C點的流場及物性分布，分別從邊界層厚度與近壁區(qū)湍動能兩個方面來解釋傳熱惡化產(chǎn)生和恢復(fù)的原因。為防止局部高溫出現(xiàn)，主要關(guān)注單側(cè)加熱下的受熱面，故圖5、圖6主要表示了加熱面單側(cè)的一半流場與物性分布。其中，A到B的過程被定義為傳熱惡化加劇的過程，由圖5（a）、（b）可以看出，從A到B點黏性底層的急劇升溫導(dǎo)致流體動力黏度增加，動力黏度增加導(dǎo)致邊界層中黏性底層增厚，如圖5（c）所示。由于黏性底層中傳熱方式以導(dǎo)熱為主，而在傳熱惡化區(qū)，與緩沖層與對數(shù)層相比，黏性底層內(nèi)熱導(dǎo)率整體處于較低的水平[25]，綜合影響下增加了黏性底層整體熱阻，降低了邊界層流體的流動換熱能力。而對于B到C過程，黏性底層內(nèi)流體動力黏度降低，邊界層厚度恢復(fù)到與A點相似的水平，從而使C點換熱能力提升，傳熱惡化恢復(fù)到正常傳熱狀態(tài)。

圖5

圖5   邊界層厚度的影響

Fig.5   Effect of boundary layer thickness

圖6

圖6   近壁區(qū)湍動能的影響

Fig.6   Effect of near-wall turbulent kinetic energy

圖6（a）可以看到，相比于A、C兩點，B點處黏性底層與緩沖層內(nèi)密度均處于較低水平，緩沖層內(nèi)尤為明顯，因此從A到B點的過程中，低密度的流體在浮力作用下加速上升，緩沖層內(nèi)出現(xiàn)明顯的熱加速現(xiàn)象，該現(xiàn)象使得B點處緩沖層與對數(shù)層內(nèi)的速度分布比較平緩，從而使得該處湍動能明顯降低，導(dǎo)致傳熱惡化發(fā)生，如圖6（b）、（c）所示。而從B到C點過程中，對數(shù)層流體密度進一步減小，流體速度增加，形成M形速度分布，增強了湍流的強度，導(dǎo)致湍動能增加，從而使傳熱惡化恢復(fù)到正常傳熱水平。

綜上所述，導(dǎo)致傳熱惡化的原因為以下兩點：（1）黏性底層內(nèi)黏度增大使得邊界層局部增厚； （2）黏性底層與緩沖層內(nèi)密度減小導(dǎo)致湍動能降低。上述耦合影響導(dǎo)致了傳熱惡化的出現(xiàn)。

2.2 常用關(guān)聯(lián)式的驗證

文獻中存在較多超臨界流體換熱的關(guān)聯(lián)式，主要是在傳統(tǒng)Dittus-Boelter關(guān)聯(lián)式的形式上引入壁面物性參數(shù)修正，以提高關(guān)聯(lián)式精度。選擇了6個常用于超臨界流體換熱的關(guān)聯(lián)式以驗證其對方通道單側(cè)加熱超臨界水的適用性。式（2）~式（7）分別來自于Gupta等[26]、Mokry等[27]、Jackson[28]、Bishop等[29]、Kim等[30]和Fan等[31]，關(guān)聯(lián)式中計算的物性是根據(jù)模擬得到的流體平均溫度與相應(yīng)固體壁面溫度，由NIST數(shù)據(jù)庫查出對應(yīng)的物性參數(shù)。

Nub=0.01Re0.89bPrbˉˉˉˉˉˉ?0.14(ρwρb)0.93(λwλb)0.22(μwμb)?1.13$N{u}_b=0.01R{e}_b^{0.89}{\overline{P{r}_b}}^{-0.14}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{0.93}{\left(\frac{{\lambda }_w}{{\lambda }_b}\right)}^{0.22}{\left(\frac{{\mu }_w}{{\mu }_b}\right)}^{-1.13}$(2)Nub=0.0061Re0.904bPrbˉˉˉˉˉˉ0.684(ρwρb)0.564$N{u}_b=0.0061R{e}_b^{0.904}{\overline{P{r}_b}}^{0.684}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{0.564}$(3)Nub=0.0183Re0.82bPrb0.5(ρwρb)0.3(cˉpcp,b)nn=?????????????0.4,TwTpc≤1或TbTpc≥1.20.4+0.2(TwTpc?1),Tb≤Tpc≤Tw0.4+0.2(TwTpc?1),1≤TbTpc≤1.2$\begin{array}{c}N{u}_b=0.0183R{e}_b^{0.82}P{r_b}^{0.5}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{0.3}{\left(\frac{{\bar{c}}_p}{c_{p,b}}\right)}^n\\ n=\left(\begin{array}{c}0.4,\frac{T_w}{T_{pc}}\le 1或\frac{T_b}{T_{pc}}\ge 1.2\\ 0.4+0.2\left(\frac{T_w}{T_{pc}}-1\right),T_b\le T_{pc}\le T_w\\ 0.4+0.2\left(\frac{T_w}{T_{pc}}-1\right),1\le \frac{T_b}{T_{pc}}\le 1.2\end{array}\right)\end{array}$(4)Nub=0.0069Re0.9bPrbˉˉˉˉˉˉ0.66(ρwρb)0.43(1+2.4Dx)$N{u}_b=0.0069R{e}_b^{0.9}{\overline{P{r}_b}}^{0.66}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{0.43}\left(1+2.4\frac{D}{x}\right)$(5)Nub=2.0514Re0.928bPr0.742b(ρwρb)1.305(μwμb)?0.669(cˉpcp,b)0.888(q+)0.792$N{u}_b=2.0514R{e}_b^{0.928}P{r}_b^{0.742}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{1.305}{\left(\frac{{\mu }_w}{{\mu }_b}\right)}^{-0.669}{\left(\frac{{\bar{c}}_p}{c_{p,b}}\right)}^{0.888}{\left(q^{+}\right)}^{0.792}$(6)Nub=0.0061Re0.904bPrbˉˉˉˉˉˉ0.684(ρwρb)0.564(μwμb)?0.184$N{u}_b=0.0061R{e}_b^{0.904}{\overline{P{r}_b}}^{0.684}{\left(\frac{{\rho }_w}{{\rho }_b}\right)}^{0.564}{\left(\frac{{\mu }_w}{{\mu }_b}\right)}^{-0.184}$(7)

圖7展示了各個關(guān)聯(lián)式的驗證結(jié)果，綜合來看，Gupta、Kim和Jackson關(guān)聯(lián)式與計算數(shù)據(jù)誤差最大，大部分結(jié)果的驗證誤差均大于30%。Bishop關(guān)聯(lián)式有部分點的驗證誤差大于30%，Mokry關(guān)聯(lián)式在小熱通量和大質(zhì)量流速時與計算數(shù)據(jù)吻合較好，驗證誤差在30%以內(nèi)，但其他工況與計算數(shù)據(jù)的符合較差，驗證誤差均大于30%。對于Fan關(guān)聯(lián)式，除高質(zhì)量流速（case 9）和大管徑（case 10）的工況的部分?jǐn)?shù)據(jù)外，絕大多數(shù)數(shù)據(jù)驗證結(jié)果誤差均在30%以內(nèi)，預(yù)測精度較高。根據(jù)2.1節(jié)的分析，傳熱惡化的發(fā)生主要由近壁區(qū)黏度和密度的變化導(dǎo)致的，而Fan關(guān)聯(lián)式引入了壁面密度和黏度的修正，因此該關(guān)聯(lián)式精度較高，推薦使用Fan關(guān)聯(lián)式來預(yù)測單側(cè)加熱方形通道內(nèi)超臨界水的換熱性能。

圖7

圖7   關(guān)聯(lián)式計算與數(shù)值模擬結(jié)果的比較

Fig.7   Comparison between the correlations and present numerical results

2.3 不同結(jié)構(gòu)強化換熱效果比較

為防止豎直向上流動工況出現(xiàn)傳熱惡化，本節(jié)對比了2.1節(jié)所述多種強化換熱結(jié)構(gòu)的強化換熱效果。圖8為G=500 kg·m-2·s-1時不同強化換熱結(jié)構(gòu)的加熱底面平均溫度隨加熱長度的變化?？梢钥闯?，對傳熱惡化工況，5種強化換熱結(jié)構(gòu)均可降低壁面溫度，其中，四種凹槽結(jié)構(gòu)和并排的雙通道均可以避免傳熱惡化，顯著降低壁溫。并排雙通道、非對稱倒角下游凹槽、倒角凹槽、非對稱倒角上游凹槽和對稱凹槽，在光滑通道傳熱惡化點處的壁溫下降幅度分別為18.55%、16.12%、14.76%、14.65%和13.29%，強化換熱效果為：并排雙通道>非對稱倒角下游凹槽＞倒角凹槽>非對稱倒角上游凹槽>對稱凹槽，上下雙通道強化換熱效果最差，不能避免傳熱惡化的產(chǎn)生。

圖8

圖8   加熱底面平均溫度隨加熱長度的變化

Fig.8   Average temperature on the heated surface against heating length

為綜合評價強化換熱結(jié)構(gòu)的換熱能力和阻力損失，引入等泵功條件下綜合換熱評價因子PEC。

PEC=Nu/Nu0(f/f0)1/3$PEC=\frac{Nu/N{u}_0}{(f/f_0{)}^{1/3}}$(8)

圖9、圖10分別表示了豎直向上流動不同強化換熱結(jié)構(gòu)Nu/Nu0與PEC隨質(zhì)量流速的變化。可以看出非對稱倒角下游凹槽綜合換熱性能最優(yōu)，非對稱倒角下游凹槽在質(zhì)量流速G=500~1000 kg·m-2·s-1時，PEC有微小的下降，但維持在一個較高的水平，且Nu/Nu0比值維持較高水平，說明在該流量區(qū)間，隨著流量增加非對稱倒角下游凹槽對流體的擾動增加，相較于光滑管，強化傳熱效果能維持較高水平。對稱凹槽與倒角凹槽結(jié)構(gòu)也具有一定的綜合強化換熱性能，并排雙通道僅僅在質(zhì)量流速G=500 kg·m-2·s-1時有一定的強化效果，其他情況下PEC<1，綜合換熱性能較差。除了上下通道，其他通道在G=500 kg·m-2·s-1時PEC能達到最大數(shù)值，是因為在該流量下，強化結(jié)構(gòu)對流體的擾動影響最大，Nu/Nu0達到最大值，此時強化傳熱效果最佳。上下雙通道的Nu/Nu0較小，強化換熱效果有限，其PEC在所有質(zhì)量流速下均小于1，綜合換熱性能不如光管。

圖9

圖9   不同強化換熱結(jié)構(gòu)的Nu/Nu0

Fig.9   Nu/Nu0 of different enhanced structures

圖10

圖10   不同強化換熱結(jié)構(gòu)的PEC

Fig.10   PEC of different enhanced structures

綜上所述，在q=1000 kW·m-2、G=200~1500 kg·m-2·s-1的工況下，雙通道結(jié)構(gòu)換熱效果有限，不適用于超臨界水的強化換熱；凹槽結(jié)構(gòu)均具有一定的強化換熱效果，其中非對稱倒角下游凹槽結(jié)構(gòu)強化換熱效果最佳。

2.4 強化換熱機理分析

本節(jié)對凹槽結(jié)構(gòu)強化換熱的機理進行詳細(xì)分析。圖11為G=500 kg·m-2·s-1時四種凹槽結(jié)構(gòu)附近切面的流線圖，可以看到凹槽區(qū)域產(chǎn)生明顯的旋渦，導(dǎo)致該區(qū)域的流動分離，從而破壞邊界層的發(fā)展，增加流體的擾動。流體在凹槽下游位置與加熱壁面再接觸，使得該處邊界層變薄，從而強化換熱。由于旋渦區(qū)域僅出現(xiàn)在凹槽結(jié)構(gòu)內(nèi)部，對主流流動影響較小，因此不會帶來較大的壓降。值得注意的是，凹槽內(nèi)部的旋渦雖然增強了流體的擾動，但由于其速度較小，因此換熱能力較弱，導(dǎo)致凹槽內(nèi)部區(qū)域換熱效果較差，強化的部位只發(fā)生在凹槽下游的再附著點附近。因此，在保證流動分離的前提下，應(yīng)盡可能減小旋渦強度以降低阻力。對稱凹槽結(jié)構(gòu)由于壁面結(jié)構(gòu)的突變，產(chǎn)生了較強的流動分離，導(dǎo)致在凹槽區(qū)域出現(xiàn)了一大一小兩個旋渦，如圖11（a）所示。較大的旋渦可以使主流產(chǎn)生流動分離從而強化換熱，但較小的旋渦對換熱并沒有積極作用，反而由于增強了流體的擾動而增大該處的流動阻力。與對稱凹槽相比，倒角凹槽壁面附近流線更為順滑，消除了凹槽區(qū)域較小的旋渦，使得流阻降低。同時，在凹槽下游部位，再附著區(qū)域由尖角變?yōu)槠交^渡的圓角，相當(dāng)于增加了強化區(qū)域，因此倒角凹槽比對稱凹槽具有更高的傳熱性能。非對稱倒角上游凹槽結(jié)構(gòu)上游陡峭而下游平緩，產(chǎn)生了一大一小兩個旋渦，由于有利于強化傳熱的旋渦產(chǎn)生于下游，而下游的區(qū)域較為狹窄，所以非對稱倒角上游強化傳熱效果一般。非對稱倒角下游凹槽結(jié)構(gòu)壁面下游陡峭而上游平緩，增強了流體的擾動，產(chǎn)生了一大三小的四個旋渦，大的旋渦對倒角凹槽結(jié)構(gòu)影響范圍更大，強化換熱效果更為明顯，小的旋渦則增加了流動阻力，綜合來看大的旋渦影響范圍更廣，因此非對稱倒角凹槽結(jié)構(gòu)強化換熱效果更好，推薦使用非對稱倒角凹槽結(jié)構(gòu)以改善傳熱惡化。

圖11

圖11   凹槽附近剖面流線圖

Fig.11   Streamlines around various grooves

為進一步揭示傳熱強化機理，引入Guo等[32]提出的場協(xié)同原理。一方面是溫度梯度場與流場的協(xié)同關(guān)系，即溫度梯度與速度的協(xié)同角β，協(xié)同角β越小，對流傳熱能力越好；另一方面是速度矢量與速度梯度矢量的協(xié)同關(guān)系，即速度與速度梯度的協(xié)同角α，協(xié)同角α越大，流動阻力越小。

圖12、圖13均為沿著流動的中心面剖面圖，分別表示G=500 kg·m-2·s-1時對稱凹槽、倒角凹槽和非對稱倒角凹槽在Z=0.25 m附近的協(xié)同角β與α的分布?？梢钥闯?，協(xié)同角β值較小的范圍：非對稱倒角下游凹槽>倒角凹槽>對稱凹槽>非對稱倒角上游凹槽，說明非對稱倒角下游凹槽結(jié)構(gòu)的傳熱能力最強。此外，對稱凹槽結(jié)構(gòu)中，β值較小的區(qū)域主要集中在凹槽內(nèi)部，但由于該區(qū)域流速較小，換熱能力較差，因此對主流換熱能力的強化效果有限。而倒角凹槽與非對稱凹槽結(jié)構(gòu)中β值較小的區(qū)域主要集中在凹槽與主流區(qū)域交界處，該處的換熱強化有利于凹槽區(qū)域和主流之間的熱量交換，且非對稱凹槽β值較小的區(qū)域向主流區(qū)沿伸幅度大于倒角凹槽，因此，非對稱倒角凹槽結(jié)構(gòu)更有利于傳熱的強化。

圖12

圖12   凹槽附近的協(xié)同角β云圖

Fig.12   Distribution of synergy angle β around various grooves

圖13

圖13   凹槽附近的協(xié)同角α云圖

Fig.13   Distribution of synergy angle α around various grooves

由圖13（a）可以看出，對稱凹槽上游和下游部位存在較大的低α區(qū)域，這是由于對稱凹槽結(jié)構(gòu)尖角的存在，流體在流動分離點和再附著點附近產(chǎn)生比較強烈的流動方向轉(zhuǎn)變[圖11（a）]，使得該處流動阻力嚴(yán)重增加。對凹槽結(jié)構(gòu)采用倒角后，可以使流動更平滑，避免較強的流動分離及流動方向轉(zhuǎn)變，從而顯著改善這兩個區(qū)域的協(xié)同性，使局部阻力顯著降低。采用倒角的結(jié)構(gòu)中低α區(qū)域影響區(qū)域：非對稱倒角下游凹槽>倒角凹槽>非對稱倒角上游凹槽，說明非對稱倒角下游凹槽結(jié)構(gòu)流動阻力最大，非對稱倒角下游凹槽結(jié)構(gòu)由于下游壁面變陡，增強了流體的擾動，強化換熱的同時也增加了局部阻力，非對稱倒角上游凹槽雖然降低了局部阻力但強化換熱效果一般，綜合PEC考慮，非對稱倒角下游凹槽結(jié)構(gòu)強化換熱效果最佳。

3 結(jié) 論

建立了超臨界水在單側(cè)加熱方通道內(nèi)流動傳熱計算模型，基于數(shù)值模擬研究，分析了超臨界流體常用換熱關(guān)聯(lián)式對超臨界水在單側(cè)加熱方管中的適用性，計算結(jié)果可進一步為發(fā)動機再生冷卻通道的設(shè)計提供參考。同時詳細(xì)分析了傳熱惡化產(chǎn)生的機理，比較不同強化換熱結(jié)構(gòu)（雙通道和凹槽）的優(yōu)劣及機理分析。主要結(jié)論如下。

（1）引起單側(cè)加熱方管內(nèi)超臨界水流動傳熱惡化主要有以下兩點原因：①黏性底層內(nèi)黏度增大使得邊界層局部增厚；②黏性底層與緩沖層內(nèi)密度減小導(dǎo)致湍動能降低。上述耦合影響導(dǎo)致傳熱惡化的出現(xiàn)。

（2）綜合比較了六種關(guān)聯(lián)式，發(fā)現(xiàn)Fan關(guān)聯(lián)式與計算結(jié)果符合程度最高。傳熱惡化的發(fā)生主要由近壁區(qū)黏度和密度的變化導(dǎo)致，而Fan關(guān)聯(lián)式引入了壁面密度和黏度的修正，因此推薦采用Fan關(guān)聯(lián)式預(yù)測單側(cè)加熱方形通道內(nèi)超臨界水的換熱性能。

（3）在傳熱惡化產(chǎn)生的區(qū)域，綜合對比了雙通道和凹槽結(jié)構(gòu)的強化傳熱效果，發(fā)現(xiàn)雙通道結(jié)構(gòu)強化傳熱效果有限，不適用于超臨界水的強化傳熱；凹槽結(jié)構(gòu)具有較好的強化傳熱效果，且引起的壓力損失小，其中非對稱倒角下游凹槽結(jié)構(gòu)強化傳熱效果最佳。

（4）場協(xié)同原理分析發(fā)現(xiàn)，對于非對稱倒角下游的凹槽結(jié)構(gòu)，其β值較小的區(qū)域占比大，且低α區(qū)域增加不明顯，說明其強化傳熱的同時，流動阻力增加不顯著，因此非對稱倒角下游凹槽結(jié)構(gòu)綜合傳熱性能最佳。

符 號 說 明