1 引言
   隨著計(jì)算機(jī)技術(shù)、數(shù)值計(jì)算技術(shù)的發(fā)展，現(xiàn)代模擬仿真技術(shù)計(jì)算流體力學(xué)（computational　fluid　dynamics，CFD）也隨之而生。它是對(duì)純理論和純實(shí)驗(yàn)方法很好的促進(jìn)和補(bǔ)充。CFD作為一門(mén)新興學(xué)科，它力求通過(guò)數(shù)值實(shí)驗(yàn)替代實(shí)物實(shí)驗(yàn)，采用虛擬流場(chǎng)來(lái)模擬真實(shí)流場(chǎng)內(nèi)部的流體流動(dòng)情況，從而使得實(shí)驗(yàn)研究更加方便，研究場(chǎng)景更加豐富可編程。
FLUENT軟件提供了多種基于非結(jié)構(gòu)化網(wǎng)格的復(fù)雜物理模型，并針對(duì)不同物理問(wèn)題的流動(dòng)特點(diǎn)創(chuàng)建出不同的數(shù)值解法。用戶(hù)可根據(jù)實(shí)際需求自由選擇，以便在計(jì)算速度、穩(wěn)定性和精度等方面達(dá)到最佳，提高設(shè)計(jì)效率。
  關(guān)于渦街流量計(jì)的發(fā)生體數(shù)值模擬研究，主要集中在渦街發(fā)生體形狀和尺寸上。Yamasaki指出發(fā)生體的形狀與幾何參數(shù)和渦街流量計(jì)的流量特性（儀表系數(shù)、線(xiàn)性度、重復(fù)性、測(cè)量范圍）與阻力特性存在相當(dāng)大的關(guān)聯(lián)關(guān)系。S.C.Luo等人研究旋渦發(fā)生體尾緣形狀以及迎流角度對(duì)渦街性能的影響，在風(fēng)洞和水槽實(shí)驗(yàn)中，得出在全長(zhǎng)相等的情況下，旋渦強(qiáng)度隨尾緣夾角的增大而減小。彭杰綱等人在50mm口徑管道氣流量實(shí)驗(yàn)中，通過(guò)對(duì)不同尾緣夾角角度的旋渦發(fā)生體進(jìn)行實(shí)驗(yàn)研究，得出旋渦發(fā)生體尾緣的夾角為41.8°時(shí)具有很好的線(xiàn)性度。賈云飛等人通過(guò)對(duì)二維渦街流場(chǎng)中的壓力場(chǎng)進(jìn)行數(shù)值仿真研究，得出T形發(fā)生體產(chǎn)生的旋渦信號(hào)的強(qiáng)度要優(yōu)于三角柱發(fā)生體。
    渦街流量計(jì)利用流體振動(dòng)原理進(jìn)行流量測(cè)量。選取了應(yīng)力式渦街流量計(jì)進(jìn)行研究。它通過(guò)壓電檢測(cè)元件獲取電壓頻率，再根據(jù)流體流量與渦街頻率成正比得出被測(cè)流量。在過(guò)去的渦街流量計(jì)研究中，一直將研究重點(diǎn)放在真實(shí)流場(chǎng)實(shí)驗(yàn)中，但這需要重復(fù)更換口徑、調(diào)節(jié)流量，大大降低了工作效率。為解決此問(wèn)題，采用三維渦街流場(chǎng)數(shù)值分析的方法對(duì)內(nèi)部流場(chǎng)的變化進(jìn)行研究。
    通過(guò)FLUENT軟件對(duì)三維渦街流場(chǎng)進(jìn)行數(shù)值仿真，并將不同流速下的升、阻力系數(shù)進(jìn)行比較，驗(yàn)證數(shù)值仿真可行性。并通過(guò)改變管截面與截流面之間的夾角，在低、中、高速流速下，進(jìn)行取壓，最終得出隨著夾角的不同，信號(hào)強(qiáng)度不同。夾角在1°～7°范圍，對(duì)信號(hào)強(qiáng)度的衰減影響不大，超過(guò)7°以后對(duì)信號(hào)強(qiáng)度影響變大，并隨著流速的增加，趨勢(shì)越來(lái)越強(qiáng)。
2 升、阻力系數(shù)
    旋渦脫落時(shí)，流體施加給柱體一個(gè)垂直于主流的周期性交變作用力，稱(chēng)為升力。由于柱體兩側(cè)交替的釋放旋渦時(shí)，剛釋放完渦流的一側(cè)柱面，擾流改善，側(cè)面總壓力降低；將要釋放渦流的另一側(cè)柱面，擾流較差，側(cè)面總壓力較大，從而形成一個(gè)作用在三角柱上、方向總是指向剛釋放完渦流那一側(cè)的作用力，所以升力的交變頻率和旋渦的脫落頻率一致，升力的變化規(guī)律和旋渦的變化規(guī)律一致，因而通過(guò)監(jiān)視柱面上的升力變化規(guī)律，可以反映旋渦脫落規(guī)律。
    阻力系數(shù)反映的是柱體迎流方向上的作用力變化情況，每當(dāng)柱體兩側(cè)不管哪一邊的釋放旋渦一次，迎流方向上的作用力都會(huì)隨壓力變化有規(guī)律地變化一次，因此，升力系數(shù)變化的一個(gè)周期內(nèi)，阻力系數(shù)變化為兩個(gè)周期。
3 三維渦街流場(chǎng)模擬的可行性分析
3.1 幾何建模與網(wǎng)格劃分
圖1是在A(yíng)NSYS　Workbench中建立的三維渦街流量計(jì)幾何模型。其中管道口徑50mm，管道長(zhǎng)1000mm，旋渦發(fā)生體截流面寬度14mm，管截面與截流面夾角為α。

對(duì)幾何模型進(jìn)行非結(jié)構(gòu)網(wǎng)格劃分，作為數(shù)值模擬的載體，如圖2所示。

3.2 仿真參數(shù)設(shè)置
在FLUENT中，三維渦街流場(chǎng)參數(shù)設(shè)置如下：
1）流體：空氣（air）；
2）湍流模型：Renormalization－group（RNG）k－ε模型；
3）邊界條件
①流速入口邊界：根據(jù)需要設(shè)置不同流速、湍流動(dòng)能和耗散率；
②壓力出口邊界：零壓；
4）求解器：基于壓力的三維雙精度瞬態(tài)求解器；
5）數(shù)值計(jì)算過(guò)程：SIMPLE算法。
3.3 升、阻力變化頻率的計(jì)算結(jié)果及分析
圖3所示速度等值。三維渦街流場(chǎng)在夾角為0°，入口流速為5m/s的情況下的速度等值線(xiàn)圖。

    通過(guò)仿真模擬，圖4給出流速u(mài)＝5m/s時(shí)，作用在三角柱上的升力系數(shù)和阻力系數(shù)變化曲線(xiàn)。由圖5升力系數(shù)的FFT曲線(xiàn)可以看出其頻率為FL＝87.92Hz。從圖6阻力系數(shù)的FFT曲線(xiàn)可以看出其頻率為FD＝176.43Hz，約為升力系數(shù)變化頻率的2倍。

    為了驗(yàn)證將FLUENT用于渦街流量計(jì)的三維流場(chǎng)仿真的可行性，對(duì)不同流速下的升、阻力頻率進(jìn)行比較，如表1所示�？梢钥闯鲎枇ο禂�(shù)變化頻率是升力系數(shù)變化頻率的2倍，說(shuō)明用FLUENT進(jìn)行渦街流量計(jì)的三維仿真是可行的。

4 仿真結(jié)果
    基于上述通過(guò)升、阻力變化頻率的關(guān)系驗(yàn)證出利用FLUENT對(duì)三維渦街流場(chǎng)進(jìn)行仿真是可行的。本節(jié)應(yīng)用FLUENT對(duì)截流夾角、流速和信號(hào)強(qiáng)度之間的關(guān)系進(jìn)行了仿真研究。分別取7m/s、40m/s和70m/s的流速α的角度在0°～10°范圍內(nèi)取值（發(fā)生體的安裝偏差一般不會(huì)超過(guò)10°），進(jìn)行數(shù)值仿真。記錄信號(hào)強(qiáng)度，如表2所示。

    將表2的數(shù)據(jù)繪制成圖7，將圖7中流速為7m/s的數(shù)據(jù)放大如圖8所示。觀(guān)察圖7、8，可以直觀(guān)的反應(yīng)出夾角、流速與信號(hào)強(qiáng)度的關(guān)系變化。通過(guò)對(duì)比這3張圖可以看出，信號(hào)強(qiáng)度隨著夾角、流速的不同而不同。并從圖中得出結(jié)論：
1）渦街的信號(hào)強(qiáng)度與流速成正比，隨著流速的增加，旋渦脫落頻率信號(hào)強(qiáng)度會(huì)顯著增加。
2）在流速相同的情況下，隨著夾角的增大，信號(hào)強(qiáng)度逐漸減小，并隨著夾角的增大，信號(hào)強(qiáng)度的衰減程度也逐漸增大。夾角在1°～7°范圍，對(duì)信號(hào)強(qiáng)度的衰減影響不大，可忽略，超過(guò)7°以后對(duì)信號(hào)強(qiáng)度影響變大，不可忽略。
3）在夾角相同的情況下，隨著流速的增大，信號(hào)強(qiáng)度衰減趨勢(shì)越來(lái)越明顯。

5 結(jié)論
    流場(chǎng)仿真在渦街流量計(jì)的設(shè)計(jì)和完善中正變得越來(lái)越重要，它通過(guò)理論支持指導(dǎo)仿真的可實(shí)施性，并將仿真結(jié)論用于實(shí)驗(yàn)中，提高效率。通過(guò)模擬三維渦街流場(chǎng)三角柱繞流現(xiàn)象，將升、阻力頻率進(jìn)行對(duì)比，驗(yàn)證了可將FLUENT用于三維渦街流場(chǎng)的仿真中。并從不同流速和不同截流夾角兩方面分別考慮，對(duì)比分析了三維渦街信號(hào)的信號(hào)強(qiáng)度，得出夾角在1°～7°范圍，對(duì)信號(hào)強(qiáng)度的影響不大，超過(guò)了7°以后影響變大。從而為以后的實(shí)驗(yàn)做出理論指導(dǎo)。進(jìn)一步的研究可以通過(guò)對(duì)不同形狀的旋渦發(fā)生體取不同截流夾角和不同流速進(jìn)行仿真對(duì)比研究。

以上內(nèi)容來(lái)源于網(wǎng)絡(luò)，如有侵權(quán)請(qǐng)聯(lián)系即刪除!