田麗亭 高永坤 王佳麗 閔春華 楊歷
(河北工業(yè)大學(xué)能源與環(huán)境工程學(xué)院 天津 300401)
三角小翼V形陣列換熱特性實驗研究
田麗亭 高永坤 王佳麗 閔春華 楊歷
(河北工業(yè)大學(xué)能源與環(huán)境工程學(xué)院 天津 300401)
在Re=4800~10200范圍內(nèi),實驗研究了矩形通道內(nèi)圓管側(cè)后方兩排三角小翼V形陣列布置的強(qiáng)化換熱特性,比較了三角小翼V形陣列布置和單排大三角小翼布置的換熱差異,研究了攻角、翼高和排列間距對換熱的影響。結(jié)果表明:三角小翼V形陣列的Nu數(shù)比單排大三角小翼平均增加4.9%,且隨著陣列中翼片前后排間距的增加,換熱呈增強(qiáng)趨勢。對于三角小翼V形陣列和單排大三角小翼,隨著小翼攻角的增加,通道換熱先增加后減小,最佳攻角均為α=60°,隨著小翼翼高的增加,通道換熱單調(diào)增加。實驗數(shù)據(jù)擬合得到了V形三角小翼陣列和單排大三角小翼的Nu數(shù)無量綱關(guān)聯(lián)式,為管翅式換熱器中三角小翼的布置提供了依據(jù)和參考。
管翅式換熱器;強(qiáng)化傳熱;縱向渦;實驗研究
管翅式換熱器在工業(yè)生產(chǎn)和日常生活中廣泛應(yīng)用,如汽車工業(yè)、空調(diào)制冷、石油化工、航空航天等諸多領(lǐng)域。通常管外側(cè)流過的工質(zhì)為氣體,在換熱器總熱阻中氣體側(cè)熱阻占據(jù)70%~90%。在低耗、環(huán)保、能源高效利用的要求下,強(qiáng)化氣體側(cè)的換熱顯得更為重要??v向渦強(qiáng)化傳熱技術(shù)是一種非常有效的無源強(qiáng)化傳熱技術(shù)[1],其流動換熱性能得到許多學(xué)者的研究。Tian L T等[2]和Wu J M等[3]對三角小翼和矩形小翼產(chǎn)生的縱向渦對矩形通道流動換熱的影響進(jìn)行了數(shù)值和實驗研究,表明三角小翼的綜合換熱性能優(yōu)于矩形小翼,且小翼攻角越大,換熱性能越好。文獻(xiàn)[4-6]對矩形通道內(nèi)其他形式的縱向渦發(fā)生器(梯形翼、斜截圓柱體、斜截橢圓柱體、八邊形翼等)的強(qiáng)化換熱特性進(jìn)行了實驗研究。由于縱向渦擾動使換熱增強(qiáng)的同時伴隨的阻力增加較小的優(yōu)點(diǎn),縱向渦強(qiáng)化傳熱技術(shù)也越來越多應(yīng)用于換熱器翅片上,O′Brien J E等[7]實驗研究發(fā)現(xiàn),在換熱管側(cè)后方的平直翅片上布置三角小翼,在低雷諾數(shù)下,翅片換熱增強(qiáng)2倍,在高雷諾數(shù)下,換熱增加也接近50%。Torii K等[8]提出一種把三角小翼對稱側(cè)置于換熱管旁邊的布置方式,這種方式可有效推遲換熱管表面邊界層分離,并使換熱管尾流覆蓋面積減少。Joardar A等[9-10]對七排換熱管的平直翅片圓管換熱器進(jìn)行了實驗和數(shù)值研究,三角小翼側(cè)置換熱管的布置方式使換熱器綜合換熱性能提高了32%。縱向渦發(fā)生器除應(yīng)用于平直翅片圓管換熱器中,文獻(xiàn)[11-14]還研究了三角小翼在圓管波紋翅片、百葉窗翅片、扁平管波紋翅片以及螺旋板換熱器中的應(yīng)用。上述文獻(xiàn)中縱向渦發(fā)生器都是強(qiáng)化了矩形通道內(nèi)或管外翅片間通道的換熱,Habchi C等[15]和車翠翠等[16]對圓管內(nèi)安置梯形翼渦發(fā)生器進(jìn)行了數(shù)值研究和PIV實驗,結(jié)果表明翼片對流體的擾動大大促進(jìn)了圓管內(nèi)對流換熱的強(qiáng)化。
三角小翼縱向渦發(fā)生器布置在管翅式換熱器中時,常見的布置方式是放置于換熱管側(cè)后方。大部分研究都是對管后一對三角小翼的強(qiáng)化換熱效果展開研究,較少研究多排小翼的強(qiáng)化換熱性能。在自然界大雁呈隊列飛行的啟發(fā)下,本文主要對換熱管側(cè)后方布置的三角小翼陣列進(jìn)行了實驗研究,對比了三角小翼V形陣列和單排大三角小翼的換熱差異,并研究了小翼攻角、翼高及排列間距對通道換熱的影響。
本實驗裝置為開式吸風(fēng)系統(tǒng),如圖1所示。主要由入口段、整流柵、前后穩(wěn)定段、實驗段、收縮段、流量測量段、離心式變頻風(fēng)機(jī)等幾部分組成。入口段為漸縮喇叭形入口,金屬網(wǎng)格整流柵使流經(jīng)的空氣得到梳理與均勻。實驗中空氣流量的調(diào)節(jié)使用變頻器控制風(fēng)機(jī)的轉(zhuǎn)速來實現(xiàn),流量測量由測量段上的渦街流量計測得。
圖1 實驗裝置示意圖Fig.1 Schematic sketch of the experimental setup
實驗段為一矩形通道,尺寸為800 mm×240 mm ×40 mm(長×寬×高),通道底面進(jìn)行電加熱。電加熱系統(tǒng)主要由220 V交流電源、變壓器、加熱膜、電壓表及電流表組成。實驗過程中調(diào)節(jié)變壓器,改變加熱膜兩端電壓來實現(xiàn)通道加熱功率的調(diào)整。加熱膜是一張800 mm×240 mm(長×寬)的聚酰亞胺加熱膜,加熱膜內(nèi)部導(dǎo)電發(fā)熱體是寬度為2 mm的金屬箔條,以間距2 mm蛇形盤旋排列,外部以聚酰亞胺薄膜為絕緣體。聚酰亞胺加熱膜厚度只有0.3 mm,熱容量小升溫迅速,加熱功率分布均勻。為了使實驗段通道底面獲得均勻的熱流,在加熱膜和通道底面間夾入了一塊厚度為2 mm的鋁板。加熱膜與鋁板、鋁板與通道底面間均涂有AS?608A型導(dǎo)熱硅脂以減小各層間接觸熱阻。為降低熱損失,整個實驗段最外側(cè)包裹橡塑保溫棉。
實驗段空氣的進(jìn)出口溫度及通道底面壁溫均采用標(biāo)定后的直徑0.2 mm的銅?康銅T型熱電偶測量。實驗段進(jìn)出口截面分別均勻布置了6根和8根熱電偶。由于圓管和三角小翼在實驗段呈對稱布置,測量壁溫的熱電偶只布置在通道中心線的一側(cè),通道底面外側(cè)共埋設(shè)了43根熱電偶。
實驗時圓管和三角小翼用耐熱膠垂直固定于通道底面上。圓管直徑為60 mm,圓心位于通道中心線上,距實驗段進(jìn)口70 mm。與來流成一定攻角的三角小翼陣列,對稱布置于圓管下游側(cè)后方。本文研究的三角小翼陣列呈V形排列,陣列由兩排三角小翼組成,每個三角小翼的翼長l=30 mm,翼高h(yuǎn)分別取20 mm、15 mm、12 mm,小翼與來流方向的攻角α分別取15°、30°、45°、60°、75°,前后排小翼間距s分別取0 mm、2 mm、4 mm。為了比較兩排三角小翼V形陣列和單排三角小翼對通道換熱的影響,本文也進(jìn)行了單排三角小翼的實驗,為了保證兩種情況下三角小翼的迎流面積相同,單排三角小翼的翼長取為2l=60 mm,其他參數(shù)與三角小翼陣列一致。實驗中小翼和圓管的相對位置如圖2所示。
圖2 小翼和圓管的相對位置Fig.2 Relative position between winglets and tube
2.1 數(shù)據(jù)處理
雷諾數(shù):
式中:uin為通道內(nèi)空氣進(jìn)口速度,m/s;de為通道進(jìn)口的當(dāng)量直徑,m;ν為空氣的運(yùn)動粘度,m2/s。
實驗中,加熱膜的加熱功率由加熱系統(tǒng)的電流和電壓計算:
通道內(nèi)空氣換熱量:
式中:cp為空氣的定壓比熱容,J/(kg·K);qm為通道內(nèi)空氣的質(zhì)量流量,kg/s;Tf,in和Tf,out分別為空氣的進(jìn)出口平均溫度,K。
理論上空氣換熱量Q2應(yīng)等于通道底面電加熱功率Q1,但是由于散熱損失和測量誤差等多種實驗條件的影響,二者并不相等。因此,取兩者的平均值Q作為計算換熱量:
熱平衡偏差:
通道底面的熱流密度:
通道底面的平均換熱系數(shù):
其中,對數(shù)平均溫差:
平均努塞爾數(shù):
式(6)~式(9)中:A為通道底面加熱面積,m2;Tw為通道壁面的平均溫度,K;λ為空氣的導(dǎo)熱系數(shù),W/(m·K)。
2.2 誤差分析
實驗中,通道內(nèi)空氣的體積流量由精度等級為1.0的渦街流量計測得,銅?康銅熱電偶的測溫誤差為±0.2℃,加熱膜的加熱功率由精度等級為0.25的電壓表和電流表測得,所有實驗工況的熱平衡偏差要求在±5%以內(nèi)。實驗中存在的系統(tǒng)誤差與測量誤差會影響實驗結(jié)果的準(zhǔn)確性,本文采用二次冪法對實驗結(jié)果進(jìn)行不確定度分析。本實驗參數(shù)范圍內(nèi),雷諾數(shù)、換熱系數(shù)和努塞爾數(shù)的不確定度分別為3%、10%、11%。
3.1 實驗系統(tǒng)校核
在進(jìn)行三角小翼V形陣列換熱實驗之前,首先對光通道和不帶三角小翼的圓管通道進(jìn)行了換熱實驗,其平均Nu數(shù)隨Re數(shù)的變化關(guān)系如圖3所示。為了對實驗系統(tǒng)進(jìn)行校核,引入光通道平均Nu數(shù)的經(jīng)驗關(guān)聯(lián)式:Dittus?Boelter式和Gnielinski式[17]。
圖3 光通道和圓管通道的實驗結(jié)果Fig.3 Experimental results of the smooth channel and the channel with tube
從圖3中可以得出:光通道平均Nu數(shù)的實驗結(jié)果與Dittus?Boelter式吻合較好,最大偏差3.6%,平均偏差1.9%,與Gnielinski式相比,實驗Nu數(shù)偏大,最大偏差16.4%,平均偏差10.5%。通過比較可知,光通道的實驗結(jié)果與經(jīng)驗關(guān)聯(lián)式的吻合情況總體良好,能夠滿足工程要求,從而驗證了實驗系統(tǒng)的可靠性和測量的準(zhǔn)確性。在實驗的Re數(shù)范圍內(nèi),光通道和圓管通道相比,圓管通道的平均Nu數(shù)比光通道增加了15.9%~19.5%。
3.2 三角小翼V形陣列對換熱的影響
圖4是三角小翼翼高h(yuǎn)=15 mm,攻角α=60°時,單排大三角小翼及三角小翼V形陣列布置通道平均Nu數(shù)隨Re數(shù)的變化規(guī)律。圖中以圓管通道平均Nu數(shù)為比較基準(zhǔn),圓管通道無三角小翼布置。從圖中可以看出,三角小翼的添加明顯增強(qiáng)了通道內(nèi)的換熱。流體掠過三角小翼產(chǎn)生的縱向渦加劇了下游流體的摻混,同時使得更多主流流體進(jìn)入到圓管下游的尾流區(qū),壓縮了換熱惡化的尾流區(qū)范圍。在實驗的Re數(shù)范圍內(nèi),和圓管通道相比,單排大三角小翼布置使通道內(nèi)的換熱平均增大16.7%~20.8%。在相同的Re數(shù)下,三角小翼V形陣列的平均Nu數(shù)都要高于單排大三角小翼,平均增加約4.9%。這是由于三角小翼陣列由兩排小翼前后呈“V”形組成,兩排小翼產(chǎn)生的渦旋相互疊加合成一個大的縱向渦,此縱向渦形成的二次流速度與單排大三角小翼形成的二次流相當(dāng),但流體流經(jīng)前后排小翼間隙時會形成一小股橫向沖刷,該小股橫向流可改善前排小翼背壁面處的局部流動死區(qū),排擠更多主流流體沖入圓管尾流區(qū),同時該股橫向流壓迫縱向渦的渦心更貼近下壁面,加劇了主流低溫流體和近壁面高溫流體的混合,使通道底面的流動邊界層和熱邊界層進(jìn)一步減薄,從而使通道換熱得到提升。對于三角小翼V形陣列來說,不同布置間距的換熱效果相差不大,隨著間距的增加,通道平均Nu數(shù)有增加的趨勢。間距s=4 mm的通道平均Nu數(shù)比間距s=0 mm大約高1.5%。
3.3 攻角對換熱的影響
圖5是三角小翼翼高h(yuǎn)=20 mm時,小翼不同攻角布置時對通道換熱的影響。從圖中看出,對于單排大三角小翼和三角小翼V形陣列,通道平均Nu數(shù)均隨著攻角α的增加而增加,當(dāng)α=60°時通道平均Nu數(shù)最大,傳熱性能最好。而當(dāng)攻角大于60°后,通道換熱反而有所減弱,α=75°的通道平均Nu數(shù)低于60°,與45°相差不多。隨著攻角的增加,三角小翼下游產(chǎn)生的縱向渦強(qiáng)度增加,渦旋變強(qiáng),渦旋波及到的通道高溫底面范圍加大,同時渦旋也吸卷更多流體進(jìn)入圓管尾流區(qū)域,減小了尾流區(qū)影響,加強(qiáng)了換熱。但當(dāng)攻角逐漸增大時,三角小翼下游產(chǎn)生的渦旋結(jié)構(gòu)中橫向渦的比重增加,縱向渦旋強(qiáng)度隨之減弱,同時,隨著主流向下游發(fā)展過程中,橫向渦比縱向渦衰減較快,橫向渦涉及到的強(qiáng)化換熱的范圍要小于縱向渦,這是換熱性能隨著攻角增加先上升后降低的原因??v向渦的發(fā)生和發(fā)展與發(fā)生器的攻角有很大的關(guān)系,存在一個縱向渦強(qiáng)度的最佳攻角范圍,實驗表明對于通道換熱來說,無論是單排大三角小翼,還是三角小翼V形陣列布置,最佳攻角均為α=60°。
圖4 三角小翼V形陣列對Nu數(shù)的影響Fig.4 Effect of the V?shaped array of delta winglets on Nu
3.4 翼高對換熱的影響
圖6是三角小翼攻角α=60°時,小翼翼高對通道換熱的影響。從圖中看出,對于單排大三角小翼和三角小翼陣列布置,通道的平均Nu數(shù)均隨著翼片高度的增加而增加。在三角小翼長度和攻角相同情況下,翼高越大,小翼對來流的阻擋越大,小翼前后的壓差隨著增大,下游產(chǎn)生的縱向渦強(qiáng)度越大,流動擾動越強(qiáng),同時擾流的影響范圍也就越大,更好的破壞邊界層的發(fā)展,使得換熱性能增加。
3.5 平均努塞爾數(shù)關(guān)聯(lián)式
根據(jù)單排大三角小翼共42組實驗點(diǎn)、三角小翼陣列共126組實驗點(diǎn),歸納實驗結(jié)果,通過多元線性回歸得到了單排大三角小翼和三角小翼陣列的傳熱特性準(zhǔn)則關(guān)聯(lián)式:
單排大三角小翼:
V形三角小翼陣列:
圖7為三角小翼通道的Nu數(shù)實驗結(jié)果與擬合關(guān)聯(lián)式結(jié)果的對比情況。由圖可知,Nu數(shù)的擬合關(guān)聯(lián)式和實驗結(jié)果吻合良好。對于單排大三角小翼和
三角小翼陣列的Nu數(shù)關(guān)聯(lián)式,經(jīng)統(tǒng)計分析可知,分別有97.6%和95.2%的關(guān)聯(lián)式結(jié)果分布在實驗結(jié)果±5%的范圍以內(nèi)。
圖5 攻角對Nu數(shù)的影響Fig.5 Effect of attack angle on Nu
圖6 翼高對Nu數(shù)的影響Fig.6 Effect of winglet height on Nu
在Re=4800~10200范圍內(nèi),對圓管側(cè)后方三角小翼V形陣列及單排大三角小翼布置的強(qiáng)化換熱特性進(jìn)行了對比實驗,并研究了小翼布置參數(shù)(攻角、翼高及排列間距)對通道換熱的影響。主要結(jié)論如下:
1)相同的迎流面積下,三角小翼V形陣列布置的通道換熱要優(yōu)于單排大三角小翼布置,增加約4.9%左右。同時,三角小翼V形陣列布置時,隨著前后排間距的增加,通道換熱有增加的趨勢。
2)單排大三角小翼和三角小翼V形陣列布置中,通道的換熱均隨著翼片攻角的增加先增加后減小,當(dāng)翼片攻角α=60°時,通道換熱最強(qiáng)。
3)單排大三角小翼和三角小翼V形陣列布置中,翼高越大,流動擾動越強(qiáng),通道換熱性能越好。
圖7 Nu數(shù)擬合結(jié)果與實驗結(jié)果的對比Fig.7 Comparison of experiment data and correlations for Nu
4)通過擬合實驗數(shù)據(jù),得到了包含翼高、攻角和翼間距在內(nèi)的單排大三角小翼和V形三角小翼陣列的傳熱特性準(zhǔn)則關(guān)聯(lián)式,并對其準(zhǔn)確性進(jìn)行了檢驗。
[1] Ahmed H E,Mohammed H A,Yusoff M Z.An overview on heat transfer augmentation using vortex generators and nanofluids:Approaches and applications[J].Renewable and Sustainable Energy Reviews,2012,16(8):5951?5993.
[2] Tian Liting,He Yaling,Lei Yonggang,et al.Numerical study of fluid flow and heat transfer in a flat?plate channel with longitudinal vortex generators by applying field synergy principle analysis[J].International Communication in Heat and Mass Transfer,2009,36(2):111?120.
[3] Wu J M,Tao W Q.Effect of longitudinal vortex generator on heat transfer in rectangular channels[J].Applied Ther?mal Engineering,2012,37:67?72.
[4] 齊承英,周國兵,曹惠玲,等.新型渦流發(fā)生器強(qiáng)化換熱實驗研究[J].工程熱物理學(xué)報,2002,23(Suppl.):173?176.(Qi Chengying,Zhou Guobing,Cao Huiling,et al.Experimental investigation of heat transfer enhancement of a new vortex generator[J].Journal of Engineering Ther?mophysics,2002,23(Suppl.):173?176.)
[5] Zhou Guobing,Ye Qiuling.Experimental investigations of thermal and flow characteristics of curved trapezoidal wing?let type vortex generators[J].Applied Thermal Engineer?ing,2012,37:241?248.
[6] 閔春華,孔祥飛,董江峰,等.矩形通道內(nèi)八邊形翼縱向渦發(fā)生器強(qiáng)化傳熱的試驗研究[J].動力工程學(xué)報,2010,30(8):607?611.(Min Chunhua,Kong Xiangfei,Dong Jiangfeng,et al.Experimental study on heat transfer in rectangular channel with longitudinal octagonal?shaped wing?type vortex generator[J].Journal of Chinese Society of Power Engineering,2010,30(8):607?611.)
[7] O’Brien J E,Sohal M S.Heat transfer enhancement for finned?tube heat exchangers with winglets[J].Journal of Heat Transfer ASME,2005,127:171?178.
[8] Torii K,Kwak K M,Nishino K.Heat transfer enhance?ment accompanying pressure?loss reduction with winglet?type vortex generators for fin?tube heat exchangers[J].In?ternational Journal of Heat and Mass Transfer,2002,45 (18):3795?3801.
[9] Joardar A,Jacobi A M.A numerical study of flow and heat transfer enhancement using an array of delta?winglet vortex generators in a fin?and?tube heat exchanger[J].Journal of Heat Transfer ASME,2007,129:1156?1167.
[10]Joardar A,Jacobi A M.Heat transfer enhancement by win?glet?type vortex generator arrays in compact plain?fin?and?tube heat exchangers[J].International Journal of Refrige?ration,2008,31:87?97.
[11]田麗亭,何雅玲,楚攀,等.不同排列方式下三角翼波紋翅片管換熱器的換熱性能比較[J].動力工程,2009,29(1):78?83.(Tian Liting,He Yaling,Chu Pan,et al.Heat transfer performance comparison of wavy finned tube heat exchanger with delta winglets under different ar?rays[J].Journal of Power Engineering,2009,29(1):78?83.)
[12]Huisseune H,T Joen C,Jaeger P D,et al.Performance enhancement of a louvered fin heat exchanger by using del?ta winglet vortex generators[J].International Journal of Heat and Mass Transfer,2013,56(1/2):475?487.
[13]Du X Z,F(xiàn)eng L L,Yang Y P,et al.Experimental study on heat transfer enhancement of wavy finned flat tube with longitudinal vortex generators[J].Applied Thermal Engi?neering,2013,50(1):55?62.
[14]張亞龍,劉升學(xué),卿德藩,等.曲面通道內(nèi)組合渦發(fā)生器的強(qiáng)化傳熱及結(jié)構(gòu)優(yōu)化研究[J].制冷學(xué)報,2014,35(3):96?101.(Zhang Yalong,Liu Shengxue,Qing De?fan,et al.Study on enhanced heat transfer and structural optimization of combination vortex generator in warped sur?face channel[J].Journal of Refrigeration,2014,35(3):96?101.)
[15]Habchi C,Russeil S,Bougeard D,et al.Enhancing heat transfer in vortex generator?type multifunctional heat ex?changers[J].Applied Thermal Engineering,2012,38:14?25.
[16]車翠翠,田茂誠.圓管內(nèi)置梯形翼片的流場特性PIV實驗[J].化工學(xué)報,2013,64(11):3976?3984.(Che Cui?cui,Tian Maocheng.PIV experiment on flow disturbance characteristics of embedded trapezoid winglets in tube[J]. CIESC Journal,2013,64(11):3976?3984.)
[17]楊世銘,陶文銓.傳熱學(xué)[M].3版.北京:高等教育出版社,1998:162?171.
田麗亭,女,博士,河北工業(yè)大學(xué)能源與環(huán)境工程學(xué)院,(022)60435795,E?mail:lttian@hebut.edu.cn。研究方向:強(qiáng)化傳熱與節(jié)能技術(shù)。
About the author
Tian Liting,female,Ph.D.,School of Energy and Environmental Engineering,Hebei University of Technology,+86 22?60435795,E?mail:lttian@hebut.edu.cn.Research fields:heat transfer en?hancement and energy saving technology.
Experimental Study on Heat Transfer Characteristics of the Delta Winglets with V?shaped Array
Tian Liting Gao Yongkun Wang Jiali Min Chunhua Yang Li
(School of Energy and Environmental Engineering,Hebei University of Technology,Tianjin,300401,China)
Water In the range of Re=4800~10200,the experiment was made to investigate the heat transfer characteristics of two rows of the delta winglets arrayed in V shape at the back side of the round tube in the rectangular channel.The heat transfer enhancement was com?pared between the delta winglets with V?shaped array and the single?row large delta winglets,and the effects of attack angle of delta wing?let,winglet height and spacing of winglet array were examined.Results show that Nu of the delta winglets with V?shaped array is 4.9% higher than that of the single?row large delta winglets,and the heat transfer in the channel has an increasing trend with the increase of the spacing in winglet array.For the delta winglets with V?shaped array and the single?row large delta winglets,with increasing the attack an?gle of the delta winglet,the heat transfer enhancement increases firstly and then decreases,and the optimal attack angle is α=60°,while,the heat transfer enhancement in both layouts monotonically increases with the increase of the winglet height.The correlations of Nu for the delta winglets with V?shaped array and the single?row large delta winglets are respectively derived by the regression analysis of ex?perimental data.These correlations can be referred for the arrangement of the delta winglets in the fin?and?tube heat exchanger.
fin?and?tube heat exchanger;heat transfer enhancement;longitudinal vortex;experimental study
TK124
A
0253-4339(2015)04-0029-07
10.3969/j.issn.0253-4339.2015.04.029
國家自然科學(xué)基金(51106041)資助項目。(The project was supported by the National Natural Science Foundation of China(No. 51106041).)
2014年11月23日