本發(fā)明涉及一種溫度場(chǎng)和熱流的重構(gòu)方法，具體涉及一種利用非接觸溫度測(cè)量對(duì)含內(nèi)部復(fù)雜邊界結(jié)構(gòu)體的溫度場(chǎng)和熱流同時(shí)重構(gòu)方法。
背景技術(shù)：
隨著航天航空技術(shù)的發(fā)展，飛行器的速度不斷提升。隨著飛行速度的提高，劇烈的氣動(dòng)加熱作用使得飛行器承受極高的熱流密度，研制高效的熱防護(hù)系統(tǒng)是高超聲速飛行器發(fā)展急需解決的問題，而準(zhǔn)確獲取熱流密度和結(jié)構(gòu)體三維溫度場(chǎng)是設(shè)計(jì)熱防護(hù)系統(tǒng)的前提和基礎(chǔ)。
當(dāng)受熱表面熱流較高，熱流計(jì)直接測(cè)量壁面的熱流密度具有較大的誤差。當(dāng)受熱結(jié)構(gòu)的結(jié)構(gòu)簡(jiǎn)單(如前緣、平板等)時(shí)，常規(guī)的導(dǎo)熱反問題計(jì)算方法可以解決這類問題，即可通過非接觸測(cè)量壁面的溫度，借助導(dǎo)熱反問題方法，反算出受熱壁面的熱流密度。但是，由于高超飛行器發(fā)動(dòng)機(jī)的結(jié)構(gòu)復(fù)雜，并且含有內(nèi)部冷卻通道，此時(shí)，常規(guī)的導(dǎo)熱反問題計(jì)算方法無法直接計(jì)算。為了實(shí)現(xiàn)結(jié)構(gòu)體三維溫度場(chǎng)和熱流密度的預(yù)測(cè)，需要探索新型的測(cè)試方法。
技術(shù)實(shí)現(xiàn)要素：
針對(duì)上述問題，本發(fā)明的目的是提供一種利用非接觸溫度測(cè)量對(duì)含內(nèi)部復(fù)雜邊界結(jié)構(gòu)體的溫度場(chǎng)和熱流同時(shí)重構(gòu)方法。
為實(shí)現(xiàn)上述目的，本發(fā)明采取以下技術(shù)方案：一種含內(nèi)部復(fù)雜邊界結(jié)構(gòu)體的溫度場(chǎng)和熱流同時(shí)重構(gòu)方法，其特征在于，包括以下步驟：
1)獲取結(jié)構(gòu)體的受熱表面溫度w1，作為導(dǎo)熱微分方程的求解邊界條件：
受熱表面溫度w1：
式中，λ為結(jié)構(gòu)體的導(dǎo)熱系數(shù)；t為結(jié)構(gòu)體的溫度場(chǎng)；t表示計(jì)算時(shí)間；q1表示結(jié)構(gòu)體的熱流密度；n為熱流密度的方向；
2)獲取結(jié)構(gòu)體的內(nèi)部局部溫度t1，t2，…，作為導(dǎo)熱微分方程的求解限定條件：
內(nèi)部局部溫度t1，t2，...：t(ω1,t)＝t1；t(ω2,t)＝t2，...(2)
式中，ωi表示結(jié)構(gòu)體的計(jì)算區(qū)域，i＝1,2,3,…；
3)根據(jù)流體與冷卻通道壁面的熱傳遞的熱流值建立冷卻通道壁面上的熱平衡方程，作為導(dǎo)熱微分方程的求解邊界條件：
內(nèi)部邊界s1，s2，...上的熱平衡方程：
式中，si為結(jié)構(gòu)體的內(nèi)部邊界，i＝1,2,3,…；hi為結(jié)構(gòu)體不同內(nèi)部邊界上的對(duì)流傳熱系數(shù)，i＝1,2,3,…；tf為結(jié)構(gòu)體內(nèi)部冷卻通道的流體溫度；
4)利用結(jié)構(gòu)體的受熱表面溫度w1，內(nèi)部邊界s1，s2，…上的熱平衡方程，以及內(nèi)部局部溫度t1，t2，…，并結(jié)合導(dǎo)熱微分方程，獲得完整的控制方程組，即式(1)-式(5)：
導(dǎo)熱微分方程：
初始條件：t(ω,t)＝t0(5)
式中，ρ為結(jié)構(gòu)體的密度；cp為結(jié)構(gòu)體的熱容；x和y表示結(jié)構(gòu)體的坐標(biāo)；t0為結(jié)構(gòu)體初始時(shí)刻的溫度；
5)求解上述控制方程組，最后同時(shí)確定結(jié)構(gòu)體的溫度場(chǎng)t和熱流密度q1。
在進(jìn)行上述步驟1)時(shí)，采用基于非接觸輻射測(cè)量方法獲取結(jié)構(gòu)體的受熱表面溫度。
在進(jìn)行上述步驟2)時(shí)，在結(jié)構(gòu)體內(nèi)部開若干個(gè)微孔，插入熱電偶或熱電阻，通過接觸式測(cè)得結(jié)構(gòu)體的內(nèi)部局部溫度。
在進(jìn)行上述步驟3)時(shí)，首先根據(jù)結(jié)構(gòu)體內(nèi)部冷卻通道的流體種類、溫度、壓力和流速等狀態(tài)參數(shù)，選擇合適的流體對(duì)流換熱計(jì)算準(zhǔn)則關(guān)聯(lián)式；然后根據(jù)冷卻通道壁面邊界上的熱平衡關(guān)系，建立結(jié)構(gòu)體不同內(nèi)部邊界s1，s2，…上的熱平衡方程。
在一個(gè)優(yōu)選的實(shí)施例中，流體對(duì)流換熱計(jì)算準(zhǔn)則關(guān)聯(lián)式采用dittus-boelter公式：nuf＝0.023re0.8pr0.4。
在進(jìn)行上述步驟5)時(shí)，通過根據(jù)共軛梯度法求解控制方程組。
本發(fā)明由于采取以上技術(shù)方案，其具有以下優(yōu)點(diǎn)：1、本發(fā)明采用基于非接觸輻射測(cè)量方法獲取結(jié)構(gòu)體高溫表面溫度場(chǎng)分布，避免了傳統(tǒng)接觸式測(cè)溫法在溫度場(chǎng)測(cè)量以及環(huán)境適應(yīng)性等方面的局限性。2、本發(fā)明中對(duì)結(jié)構(gòu)體的內(nèi)部開一些微孔，用于測(cè)量結(jié)構(gòu)體的內(nèi)部局部溫度，此溫度將用于導(dǎo)熱微分方程求解時(shí)的限定條件，可以有效地提升方程的求解精度。3、本發(fā)明根據(jù)結(jié)構(gòu)體內(nèi)部冷卻通道的流體種類、溫度、壓力、流速等狀態(tài)參數(shù)，利用流體對(duì)流換熱計(jì)算準(zhǔn)則關(guān)聯(lián)式，建立結(jié)構(gòu)體不同內(nèi)部邊界上的熱平衡方程，用于導(dǎo)熱微分方程的邊界條件，使得方程組滿足封閉求解條件。4、本發(fā)明利用前述方法，獲得結(jié)構(gòu)體的受熱表面溫度、內(nèi)部邊界上的熱平衡方程、壁面內(nèi)部局部溫度，通過求解導(dǎo)熱微分方程，即可同時(shí)獲得結(jié)構(gòu)體的溫度場(chǎng)和熱流密度。
附圖說明
圖1本發(fā)明提供的含內(nèi)部復(fù)雜邊界結(jié)構(gòu)體的溫度場(chǎng)和熱流同時(shí)重構(gòu)方法示意圖。
具體實(shí)施方式
下面結(jié)合附圖和實(shí)施例對(duì)本發(fā)明進(jìn)行詳細(xì)的描述。然而應(yīng)當(dāng)理解，附圖的提供僅為了更好地理解本發(fā)明，它們不應(yīng)該理解成對(duì)本發(fā)明的限制。
本發(fā)明提供的含內(nèi)部復(fù)雜邊界結(jié)構(gòu)體的溫度場(chǎng)和熱流同時(shí)重構(gòu)方法，其包括以下步驟：
1)如圖1所示，采用基于非接觸輻射測(cè)量方法獲取結(jié)構(gòu)體的受熱表面溫度w1，作為導(dǎo)熱微分方程的求解邊界條件：
受熱表面溫度w1：
式中，λ為結(jié)構(gòu)體的導(dǎo)熱系數(shù)；t為結(jié)構(gòu)體的溫度場(chǎng)；t表示計(jì)算時(shí)間；q1表示結(jié)構(gòu)體的熱流密度；n為熱流密度的方向。
2)在結(jié)構(gòu)體內(nèi)部開若干個(gè)微孔，插入熱電偶或熱電阻，通過接觸式測(cè)得結(jié)構(gòu)體的內(nèi)部局部溫度t1，t2，…，作為導(dǎo)熱微分方程的求解限定條件：
內(nèi)部局部溫度t1，t2，...：t(ω1,t)＝t1；t(ω2,t)＝t2，...(2)
式中，ωi表示結(jié)構(gòu)體的計(jì)算區(qū)域，i＝1,2,3,…。
3)結(jié)構(gòu)體內(nèi)部冷卻通道的壁面與冷卻通道內(nèi)的流體進(jìn)行熱傳遞，可以根據(jù)流體與冷卻通道壁面的熱傳遞的熱流值建立冷卻通道壁面上的熱平衡方程，作為導(dǎo)熱微分方程的求解邊界條件，具體過程為：首先根據(jù)結(jié)構(gòu)體內(nèi)部冷卻通道的流體種類、溫度、壓力和流速等狀態(tài)參數(shù)，選取合適的流體對(duì)流換熱計(jì)算準(zhǔn)則關(guān)聯(lián)式(比如dittus-boelter公式：nuf＝0.023re0.8pr0.4)；然后根據(jù)冷卻通道壁面邊界上的熱平衡關(guān)系，建立結(jié)構(gòu)體不同內(nèi)部邊界s1，s2，…上的熱平衡方程：
內(nèi)部邊界s1，s2，...上的熱平衡方程：
式中，si為結(jié)構(gòu)體的內(nèi)部邊界，i＝1,2,3,…；hi為結(jié)構(gòu)體不同內(nèi)部邊界上的對(duì)流傳熱系數(shù)，i＝1,2,3,…；tf為結(jié)構(gòu)體內(nèi)部冷卻通道的流體溫度。
4)利用結(jié)構(gòu)體的受熱表面溫度w1，內(nèi)部邊界s1，s2，…上的熱平衡方程，以及內(nèi)部局部溫度t1，t2，…，并結(jié)合導(dǎo)熱微分方程，就可以獲得完整的控制方程組，即式(1)-式(5)：
導(dǎo)熱微分方程：
初始條件：t(ω,t)＝t0(5)
式中，ρ為結(jié)構(gòu)體的密度；cp為結(jié)構(gòu)體的熱容；x和y表示結(jié)構(gòu)體的坐標(biāo)；t0為結(jié)構(gòu)體初始時(shí)刻的溫度。
5)通過根據(jù)共軛梯度法(cgm)或者其他算法求解上述完整的控制方程組，最后同時(shí)確定結(jié)構(gòu)體的溫度場(chǎng)t和熱流密度q1。需要說明的是，由于cgm方法為本領(lǐng)域技術(shù)人員的公知常識(shí)，因此求解過程不再贅述。
上述各實(shí)施例僅用于說明本發(fā)明，其中各部件的結(jié)構(gòu)、連接方式和制作工藝等都是可以有所變化的，凡是在本發(fā)明技術(shù)方案的基礎(chǔ)上進(jìn)行的等同變換和改進(jìn)，均不應(yīng)排除在本發(fā)明的保護(hù)范圍之外。