帶有剪切應(yīng)力輸運(yùn)性質(zhì)的k-ε兩方程湍流模型構(gòu)造與應(yīng)用
發(fā)布時(shí)間:2021-02-12 10:22
由于傳統(tǒng)的壁面衰減函數(shù)峰值并不考慮雷諾數(shù)的變化,因此僅能在特定的雷諾數(shù)下還原邊界層的近壁衰減效應(yīng)。為了改善這個(gè)缺陷,一個(gè)結(jié)合不同湍流雷諾數(shù)的因子被引入壁面衰減函數(shù),從而使其峰值對(duì)雷諾數(shù)的變化能夠合理感知。另外,基于DNS數(shù)據(jù),將Bradshaw假設(shè)中的常數(shù)進(jìn)行了重新標(biāo)定,引入平滑的湍渦黏性系數(shù)切換函數(shù),最終形成了一個(gè)帶有剪切應(yīng)力輸運(yùn)性質(zhì)的k-ε兩方程湍流模型。通過(guò)平板、翼型、二維鼓包以及翼身組合體進(jìn)行了校驗(yàn),結(jié)果顯示新模型的精度較高,具有一定的工程實(shí)用價(jià)值。
【文章來(lái)源】:空氣動(dòng)力學(xué)學(xué)報(bào). 2019,37(03)北大核心
【文章頁(yè)數(shù)】:7 頁(yè)
【部分圖文】:
圖1基于DNS和標(biāo)準(zhǔn)k-ε模型求得的渦黏性系數(shù)μtFig.1Eddy-viscositycoefficientgotbyDNS
NS數(shù)據(jù),將三種常見(jiàn)湍流模型(Abidk-ε[17]、Yang-Shihk-ε[18]和Launder-Sharmak-ε[19])中的WDF函數(shù)通過(guò)式(4)求解出來(lái)。如圖2,隨著雷諾數(shù)的變化,三種模型的WDF函數(shù)并未隨著一同改變,而是被機(jī)械地限制在了值為1附近。通過(guò)表1可以看出,三種模型的WDF函數(shù)在自變量趨于無(wú)窮時(shí),峰值均為1,然而從DNS數(shù)據(jù)反求出的WDF峰值是隨著雷諾數(shù)不同而變化的。圖2基于DNS數(shù)據(jù)求解的不同模型的WDF函數(shù)[16]Fig.2WDFfunctionofdifferentmodelsbasedonaprior-DNSdata[16]第3期王偉等:帶有剪切應(yīng)力輸運(yùn)性質(zhì)的k-ε兩方程湍流模型構(gòu)造與應(yīng)用124
Rek/100-0.45)(1.2Rek/Reτ-0.408)·[1-exp(-Rek/42)]fhybrid=1-exp(-Rek/10)fμ=min(1-tanh(Rek/100-0.45)fμ1+max(0,fμ2),1.1)烅烄烆(5)式(5)的構(gòu)造理念在于———fμ1是將現(xiàn)有的WDF進(jìn)行聯(lián)合,以此達(dá)到近壁面附近更精確的衰減作用,通過(guò)fμ2控制峰值的增減。圖3為不同雷諾數(shù)槽道中Rek/Reτ的表達(dá)。如圖4所示,新構(gòu)造的WDF函數(shù)(a)Changefactorwiththechangeofy+(b)NewWDF’sstructurefunctionfμ2圖3新構(gòu)造的WDF帶入DNS數(shù)據(jù)中的變化Fig.3ChangesofnewlyconstructedWDFintoDNSdata圖4新WDF函數(shù)在不同雷諾數(shù)下的變化Fig.4VariationofnewWDFfunctionatdifferentReynoldsnumbers的峰值與雷諾數(shù)變化具有一種自適應(yīng)關(guān)聯(lián),其中fnc代表非均衡湍流識(shí)別函數(shù)。1.3改進(jìn)的湍渦黏性系數(shù)計(jì)算方式一般情況下,在k-ε兩方程模型中,湍渦黏性系數(shù)的確定方式為:μt=Cμfμk2/ε(6)式(6)是從均勻各向同性平衡湍流中推導(dǎo)而來(lái),因此對(duì)于帶有分離和逆壓梯度的流場(chǎng)適應(yīng)性不強(qiáng)。Menter[6]的研究表明,在逆壓梯度的流動(dòng)中,湍
【參考文獻(xiàn)】:
期刊論文
[1]A one-equation turbulence model for recirculating flows[J]. Yang Zhang,Jun Qiang Bai,Jing Lei Xu,Yi Li. Science China(Physics,Mechanics & Astronomy). 2016(06)
[2]提高k-ω SST模型對(duì)翼型失速特性的模擬能力[J]. 文曉慶,柳陽(yáng)威,方樂(lè),陸利蓬. 北京航空航天大學(xué)學(xué)報(bào). 2013(08)
[3]New perspective in statistical modeling of wall-bounded turbulence[J]. Fazle Hussain. Acta Mechanica Sinica. 2010(06)
本文編號(hào):3030701
【文章來(lái)源】:空氣動(dòng)力學(xué)學(xué)報(bào). 2019,37(03)北大核心
【文章頁(yè)數(shù)】:7 頁(yè)
【部分圖文】:
圖1基于DNS和標(biāo)準(zhǔn)k-ε模型求得的渦黏性系數(shù)μtFig.1Eddy-viscositycoefficientgotbyDNS
NS數(shù)據(jù),將三種常見(jiàn)湍流模型(Abidk-ε[17]、Yang-Shihk-ε[18]和Launder-Sharmak-ε[19])中的WDF函數(shù)通過(guò)式(4)求解出來(lái)。如圖2,隨著雷諾數(shù)的變化,三種模型的WDF函數(shù)并未隨著一同改變,而是被機(jī)械地限制在了值為1附近。通過(guò)表1可以看出,三種模型的WDF函數(shù)在自變量趨于無(wú)窮時(shí),峰值均為1,然而從DNS數(shù)據(jù)反求出的WDF峰值是隨著雷諾數(shù)不同而變化的。圖2基于DNS數(shù)據(jù)求解的不同模型的WDF函數(shù)[16]Fig.2WDFfunctionofdifferentmodelsbasedonaprior-DNSdata[16]第3期王偉等:帶有剪切應(yīng)力輸運(yùn)性質(zhì)的k-ε兩方程湍流模型構(gòu)造與應(yīng)用124
Rek/100-0.45)(1.2Rek/Reτ-0.408)·[1-exp(-Rek/42)]fhybrid=1-exp(-Rek/10)fμ=min(1-tanh(Rek/100-0.45)fμ1+max(0,fμ2),1.1)烅烄烆(5)式(5)的構(gòu)造理念在于———fμ1是將現(xiàn)有的WDF進(jìn)行聯(lián)合,以此達(dá)到近壁面附近更精確的衰減作用,通過(guò)fμ2控制峰值的增減。圖3為不同雷諾數(shù)槽道中Rek/Reτ的表達(dá)。如圖4所示,新構(gòu)造的WDF函數(shù)(a)Changefactorwiththechangeofy+(b)NewWDF’sstructurefunctionfμ2圖3新構(gòu)造的WDF帶入DNS數(shù)據(jù)中的變化Fig.3ChangesofnewlyconstructedWDFintoDNSdata圖4新WDF函數(shù)在不同雷諾數(shù)下的變化Fig.4VariationofnewWDFfunctionatdifferentReynoldsnumbers的峰值與雷諾數(shù)變化具有一種自適應(yīng)關(guān)聯(lián),其中fnc代表非均衡湍流識(shí)別函數(shù)。1.3改進(jìn)的湍渦黏性系數(shù)計(jì)算方式一般情況下,在k-ε兩方程模型中,湍渦黏性系數(shù)的確定方式為:μt=Cμfμk2/ε(6)式(6)是從均勻各向同性平衡湍流中推導(dǎo)而來(lái),因此對(duì)于帶有分離和逆壓梯度的流場(chǎng)適應(yīng)性不強(qiáng)。Menter[6]的研究表明,在逆壓梯度的流動(dòng)中,湍
【參考文獻(xiàn)】:
期刊論文
[1]A one-equation turbulence model for recirculating flows[J]. Yang Zhang,Jun Qiang Bai,Jing Lei Xu,Yi Li. Science China(Physics,Mechanics & Astronomy). 2016(06)
[2]提高k-ω SST模型對(duì)翼型失速特性的模擬能力[J]. 文曉慶,柳陽(yáng)威,方樂(lè),陸利蓬. 北京航空航天大學(xué)學(xué)報(bào). 2013(08)
[3]New perspective in statistical modeling of wall-bounded turbulence[J]. Fazle Hussain. Acta Mechanica Sinica. 2010(06)
本文編號(hào):3030701
本文鏈接:http://sikaile.net/kejilunwen/lxlw/3030701.html
最近更新
教材專著