齊保良楊秀芝汪晴晴

(山東建筑大學 信息與電氣工程學院，山東 濟南250101)

0 引言

太陽能熱水器是日常生活中必不可少的加熱裝置，出水溫度的穩(wěn)定性控制是滿足洗浴舒適性的關鍵，而混水閥是決定其性能高低的關鍵部件。目前，常用的混水閥有手動調(diào)節(jié)與熱傳感器混水閥兩大類。其中，手動調(diào)節(jié)混水閥結構簡單，但在壓差較大時會造成出水溫度調(diào)節(jié)困難，若進入混水閥的冷水較熱水壓力大，會經(jīng)常出現(xiàn)冷水通過混水閥流入熱水器的串水現(xiàn)象以及出水溫度不能及時調(diào)節(jié)到所需溫度等情況，并且調(diào)節(jié)過程長還會造成水力資源的浪費。而熱傳感器恒溫閥本質(zhì)上是一種機械式恒溫閥，其內(nèi)部結構復雜，連接和密封的部分較多而容易造成漏水，且使用壽命低[1-2]。

BIGLIARDI等[3]和王金海等[4]采用計算流體動力學軟件Fluent對壓力調(diào)節(jié)閥進行模擬仿真，研究了對閥門流量特性影響的因素；WU等[5]采用Fluent模擬仿真了汽車供油系統(tǒng)中的壓力調(diào)節(jié)閥，得到了壓力調(diào)節(jié)閥的流量隨壓力變化的曲線。隨著計算流體動力學CFD(Computational Fluid Dynamics)技術的發(fā)展，其優(yōu)越性越加凸顯，已廣泛應用于調(diào)節(jié)閥內(nèi)部流場的研究，為調(diào)節(jié)閥結構設計優(yōu)化帶來了很大的便捷，但有關太陽能熱水器混水閥在水壓變化情況下瞬態(tài)流場問題的研究相對較少，因此對太陽能熱水器混水閥瞬態(tài)流場的研究具有較大意義和價值。

文章在研究傳統(tǒng)太陽能熱水器混水閥的基礎上，設計改進了太陽能熱水器混水閥機械結構，以期能夠解決混水閥出水溫度調(diào)節(jié)困難的問題，實現(xiàn)水溫的精確控制；采用Fluent平臺對改進混水閥壓力平衡特性進行數(shù)值模擬研究，并從流體動力學方面驗證混水閥的效果。

1 太陽能熱水器恒溫混水閥的結構設計及其工作原理

1.1 太陽能熱水器恒溫混水閥的結構

文章設計的太陽能熱水器恒溫混水閥是在傳統(tǒng)混水閥的基礎上增加了水壓平衡單元，其主要由橡膠塞、連桿、波紋膜片等構成。恒溫混水閥剖面圖如圖1所示，波紋膜片將恒溫混水閥冷熱水腔室隔開，其一個邊界與閥芯連接，其余邊界則與冷熱水腔交界處的閥體壁面連接，波紋膜片與閥芯、閥體接觸的部分均使用密封膠墊密封。波紋膜片與帶橡膠塞的連桿連接，構成了水壓平衡單元。

設計水壓平衡單元結構時，熱水側橡膠塞位于熱水進水管與熱水腔連接口左側，且橡膠塞右側面與熱水進水管的管壁相切；冷水側橡膠塞與熱水側對稱設計。當冷水或熱水突然中斷時，為了確保橡膠塞完全阻塞水管，橡膠塞厚度需要稍大于進水管的直徑，同時考慮到橡膠塞隨波紋膜片運動的靈敏度需求，將橡膠塞厚度設置為進水管直徑的1.2倍。

1.2 太陽能熱水器恒溫混水閥的工作原理

熱水進水管水壓不變，當冷水進水管水壓變大時，冷水進入冷水腔的流量增加，此時波紋膜片受壓變形，通過連桿帶動橡膠塞向熱水腔側運動，從而使冷水進水口流量減少，保證流量與冷水壓變化之前冷水進入冷水腔的流量一致，從而平衡冷熱水水壓；同樣地，當冷水進水管水壓變小時，連桿隨波紋膜片受壓變形帶動橡膠塞向冷水腔側運動，使熱水流量減少，從而平衡冷熱水壓力。若冷水進水管壓力為零時，在熱水端壓力下，波紋膜片會帶動橡膠塞瞬間移動到極限位置，阻塞熱水進水管，避免燙傷；若熱水突然中斷，在冷水端壓力下，波紋膜片會帶動橡膠塞瞬間阻塞冷水進水管，切斷冷水，避免冷激。

圖1 太陽能熱水器恒溫混水閥整體外觀結構及其剖面圖

當冷水或者熱水壓力變化時，為了能及時改變冷、熱水進水口的流量，設計冷熱水側的橡膠塞側面與進水管的管壁為“相切”結構，其工作原理為:當冷水壓力增大的瞬間，波紋膜片將帶動連桿將橡膠塞向熱水側移動，冷水側的橡膠塞與進水管的管壁“相切”能更及時更大限度地縮小冷水進水口直徑，減少冷水的流量；由于熱水側橡膠塞在熱水進水管口的左側，向熱水側的移動不會改變熱水進水口大小與流量。因此，此結構能夠自動調(diào)節(jié)冷水端壓力變大時冷水的流量，從而維持出水溫度的穩(wěn)定性。

2 基于Fluent的恒溫混水閥壓力平衡特性分析方法

文章基于Fluent與前處理軟件ICEM CFD進行平衡分析，ICEM CFD為與Fluent配套的網(wǎng)格劃分軟件，允許用戶對網(wǎng)格進行修改(粗化/細化)，使模擬流動、湍流的能力與實際情況更接近[6]。

2.1 恒溫混水閥波紋膜片的設計與計算

2.1.1 波紋膜片的選擇

影響波紋膜片特性的因素主要為膜片波紋深度、膜片厚度、材料力學性能等[7-8]。在選擇波紋膜片時要依次確定波紋膜片型面和厚度。波紋膜片有三角形、梯形、正弦等波紋形狀，正弦紋波性能比梯形波紋性能好，而次于鋸齒紋波性能。但是在鋸齒波紋成形時，對波紋膜片拉伸比較大，存在應力集中的情況，并且波紋膜片在混水腔中承受交變力作用[9]，因此，為延長混水閥使用壽命，文章將波紋膜片型面制作成跨距大、高度小的正弦波紋。

波紋膜片剛度與厚度成正比，但是硬度過高會使膜片易受疲勞損傷，而厚度過小則機械強度過小，不能滿足恒溫混水閥耐壓強度要求[10-11]。具體波紋膜片厚度選擇可通過Fluent仿真平臺數(shù)值模擬波紋膜片在恒溫混水閥中的受力情況。

2.1.2 基于Fluent仿真平臺的波紋膜片撓度分析計算

(1)波紋膜片邊界條件的設定與參考點的布置

恒溫混水閥水腔內(nèi)水壓力隨冷水進水管壓力不斷變化，導致水腔內(nèi)部水流情況較為復雜，在Fluent軟件進行模型求解的過程中，選擇壓強作為波紋膜片的邊界條件，其余邊界條件見表1。對冷水水壓波動所引起的波紋膜片表面壓力分布情況進行數(shù)值模擬，波紋膜片上距離波紋膜片中心10 mm處的點所受的壓差最大，最能反映波紋膜片的耐壓情況，選此為參考點進行數(shù)據(jù)觀察。

表1 邊界條件表

(2)波紋膜片撓度的分析計算

波紋膜片邊界與恒溫混水閥冷熱水腔交界處內(nèi)壁連接，波紋膜片實際承受壓力是冷熱水腔內(nèi)壓力差。自來水壓力為0.1～0.6 MPa，通常取0.3 MPa。太陽能熱水器熱水壓力一般為0.05～0.6 MPa，＞0.6 MPa水壓的用戶一般需要安裝減壓閥[12]，因此太陽能熱水器混水閥冷熱水兩側的壓差為0～0.6 MPa。在Fluent軟件對波紋膜片進行模型求解時，考慮到安全裕度因素，選擇0.05～1 MPa作為模擬波紋膜片所承受的冷熱水壓力差，對波紋膜片周邊增加固定約束，在其兩側增加軸對稱約束。

根據(jù)設置的邊界條件進行迭代計算，并將求解得到的數(shù)據(jù)導入到函數(shù)繪圖軟件Origin軟件進行數(shù)據(jù)處理，描繪出參考點處(距離波紋膜片中心10 mm)不同厚度的波紋膜片壓強—撓度曲線，結果如圖2所示。

圖2 參考點處不同厚度波紋膜片的壓強—撓度曲線圖

由圖2可知，在相同壓力點下，隨著波紋膜片厚度增加，撓度逐漸降低。在滿足機械強度要求的前提下，紋波膜片厚度越薄，隨著冷熱水壓差的變化反應越靈敏，更有利于拉動連桿、帶動橡膠塞隨波紋膜片的形變而移動。厚度為0.2～0.5 mm的波紋膜片均能夠承受0.05～1 MPa的壓力，可以滿足恒溫混水閥內(nèi)冷熱水壓差為0～0.6 MPa時的機械強度要求，但是在滿足機械強度要求的前提下，波紋膜片厚度應選擇較小值為宜。因此，文章選擇厚度為0.2 mm的波紋膜片進行數(shù)值模擬計算。

2.2 恒溫混水閥計算模型的建立

標準k-?雙方程[13]為工程實際中最常用的湍流模型，其在湍流動能方程的基礎上引入一個關于湍流動能耗散率的方程。湍流動能k及耗散率?的微分方程分別由式(1)和(2)表示為

式中ρ為流體密度，kg/m3；k為湍流動能，k＝ui＋為湍流黏度為流體動力黏度，Pa·s；σk為湍流動能對應的普朗特系數(shù)，在Fluent軟件中默認值為1.0；Gk為由平均速度梯度引起的湍流動能，J；ε為湍流耗散率，ε＝ui＋為湍流動耗散率所對應的普朗特系數(shù)，在Fluent軟件中默認值為1.3；C1ε、C2ε為經(jīng)驗常數(shù)，在Fluent中默認值為C1ε＝1.44、C2ε＝1.92；Cμ為經(jīng)驗常數(shù)，Cμ＝0.09。

雷諾數(shù)Re為表征流體流動情況的無量綱數(shù)，其計算公式由式(3)[14]表示為

式中v為流場的特征速度，m/s；L為流場的特征長度，m。

參考曼德束MSJC混水閥的技術參數(shù)(出水口流量4 m3/h、通徑DN25等)，變量取值為ρ＝1 000 kg/m3、μ＝0.001 Pa·s、v＝2.27 m/s。由式(3)計算出雷諾數(shù)為227 000。已知雷諾數(shù)＞10 000數(shù)量級為湍流，在Fluent軟件進行模型求解的過程中，應該選擇處理湍流的標準k-?雙方程為計算模型。

3 基于Fluent的恒溫混水閥模擬仿真結果分析

Fluent軟件對恒溫混水閥的模擬仿真與對波紋膜片的處理相同，在此不再贅述。實驗中通過轉動混水閥手柄調(diào)節(jié)冷熱水混合比來調(diào)節(jié)出水溫度，恒溫混水閥的仿真結果主要通過靜壓力云圖與速度矢量圖進行描述，如圖3所示。圖中的不同顏色分別代表了壓強與流速的大小[15]，其值與光譜波長成正比，按紅橙黃綠青藍等顏色順序而逐漸減小。

(1)冷熱水進水壓力相同時恒溫混水閥仿真結果分析

由圖3(a)可知，冷熱水兩側的光譜波長呈中心軸對稱分布，表明冷熱水兩側的壓強對應相同。由圖3(b)可知，冷熱水兩側的光譜波長亦呈中心軸對稱分布，表明冷熱水側的水流速度對應相同。

圖3 冷熱水進水壓力相同時恒溫混水閥仿真結果圖

(2)冷水進水壓力＞熱水進水壓力時恒溫混水閥的仿真結果分析

由圖4(a)可知，冷水側的光譜波長大于熱水側相對位置的光譜波長，表明冷水側的水壓力大于熱水側對應位置的水壓力；波紋膜片向熱水側彎曲變形，混合出水口處冷熱水對應位置的光譜波長相同，表明在混合出水口處冷熱水對應位置的壓力相同。由圖4(b)可知，波紋膜片發(fā)生形變，混合出水口處冷熱水對應位置的光譜波長相同，表明在混合出水口處冷熱水對應位置的水流流速相同。

圖4 冷水進水壓力＞熱水進水壓力時恒溫混水閥的仿真結果圖

(3)冷水進水壓力＜熱水進水壓力時恒溫混水閥的仿真結果分析

由圖5可知，冷水側的光譜波長小于熱水側相對位置的光譜波長，波紋膜片向冷水側彎曲變形，但在混合出水口處冷熱水對應位置的光譜波長相同，表明在混合出水口處冷熱水對應位置的壓力和水流流速都相同。

圖5 冷水進水壓力＜熱水進水壓力時恒溫混水閥的仿真結果圖

(4)與傳統(tǒng)混水閥的仿真結果對比

當冷水水壓變化時，混水閥混合出水口處平均水壓隨冷水水壓的變化而變化。由圖6可知，當冷水進水口壓力變化時，改進混水閥出水口冷熱水兩側截面平均壓強發(fā)生突變，由于水壓平衡單元的調(diào)節(jié)作用，使混合出水口截面處冷熱水兩側平均壓強相同。而沒有水壓平衡單元調(diào)節(jié)作用的傳統(tǒng)混水閥，與改進混水閥對應截面區(qū)域側的平均壓強不同，且與改進混水閥相對應的截面區(qū)域冷熱水兩側壓強波動規(guī)律不同，使出水溫度隨冷水進水口水壓的變化而波動。

圖6 混水閥出水口處冷熱水兩側截面平均壓強曲線圖

根據(jù)表1設置的邊界條件進行求解，當熱水和冷水進水溫度設置在78.85、26.85℃時，迭代大約160次后，混水閥混合出水口截面平均溫度檢測曲線如圖7所示。由此圖可知，當冷水進水口壓力變大時，混水閥出水口截面平均水溫突然下降，但經(jīng)過水壓平衡單元的調(diào)節(jié)，混水閥出水溫度快速達到理想的穩(wěn)定狀態(tài)，其值約穩(wěn)定為38.5℃。

圖7 改進混水閥出水口截面平均溫度檢測曲線圖

綜上所述，當冷熱水進水壓力相同時，混水閥冷熱水兩側對應位置的壓強和水流速度相同。通過轉動混水閥手柄調(diào)節(jié)冷熱水混合比，調(diào)節(jié)出水溫度。當冷熱水進水壓力不同時，通過水壓平衡單元平衡冷熱水腔內(nèi)的水壓可以解決混水閥出水溫度隨外界水壓變化而變化的問題。

4 結論

文章設計了一種新型結構的太陽能熱水器恒溫混水閥，借助Fluent仿真平臺對改進恒溫混水閥進行了流場的數(shù)值模擬仿真，通過靜壓力云圖、速度矢量圖和混水閥出水口截面平均溫度檢測曲線分析，可以得到如下結論:

(1)當冷熱水進水壓力不相同時，水壓平衡單元可以平衡冷熱水水腔內(nèi)的水壓，保證冷水流量與冷水壓力變化前一致，即保證冷水流量不變。

(2)在熱、冷水入口溫度分別為76.85、26.85℃的情況下，當冷水壓力變大時，改進混水閥能使出水溫度快速調(diào)節(jié)到理想的穩(wěn)定狀態(tài)，該穩(wěn)定值約為38.5℃。