基于全隱式算法數(shù)值模擬空腔流動振蕩特性研究.pdf

1、空腔流動是航空工程實際中常遇到的問題，如飛行器武器艙、起落架艙等等，具有很高的工程應(yīng)用價值，另外，由于空腔流動涉及到非定常流、流體動力不穩(wěn)定性、聲—渦相互作用等許多基本的前沿問題，因而也具有廣泛的理論研究意義。因此，不論理論方面還是工程應(yīng)用方面，空腔流動問題目前已成為國際上空氣動力學(xué)領(lǐng)域的一個研究熱點。首先，在多塊結(jié)構(gòu)網(wǎng)格上，本文研究和發(fā)展了Euler方程與Navier-Stokes方程的全隱式無分裂算法?；谌[式無分裂方法，本文數(shù)值

2、模擬研究了空腔流動自激勵引起的剪切層振蕩特性。同時，本文還探究了多塊結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)上的并行計算方法。其主要內(nèi)容包括以下幾個方面： 本文基于多塊結(jié)構(gòu)網(wǎng)格，研究和發(fā)展了Euler方程與Navier-Stokes方程的全隱式無分裂算法。對流項的離散運用Roe格式，粘性項的離散采用中心格式。在每一次隱式時間迭代中，運用GMRES*方法求解隱式離散引起的稀疏線性方組，預(yù)處理方法采用零階塊不完全分解方法。對于非定常計算，時間層離散采用虛

3、時間方法。為了減少內(nèi)存需求以及矩陣與向量之間的運算操作數(shù)，本文發(fā)展了Jacobian 矩陣一種的逼近方法，這可以大大減小工作量。 在PVM并行環(huán)境下，研究了全隱式無分裂方法的多塊分布式并行計算算法。為提高并行效率，本文給出了一種負載平衡方法，負載平衡前后的墻上時間、并行效率的比較驗證了本文方法的有效性。針對本文研究所用的大型并行機特性，基于最小二乘原理，給出了一種估算計算時間成本和通訊時間成本的方法，并行算法的并行性能統(tǒng)計表明

4、該估算方法的有效性及合理性。另外，本文還探討了非結(jié)構(gòu)網(wǎng)格的并行生成技術(shù)，改進了“波陣面”區(qū)域分裂算法以使子網(wǎng)格及其邊界更有益于網(wǎng)格并行生成；并將改進的“波陣面”區(qū)域分裂算法應(yīng)用于非結(jié)構(gòu)網(wǎng)格區(qū)域的劃分，探討了非結(jié)構(gòu)網(wǎng)格上求解Euler方程的并行算法。 基于全隱式無分裂方法，首先，本文數(shù)值研究了某開式空腔的超音速、跨音速、亞音速流動，詳細探討了入口邊界條件對空腔流自持性振蕩的影響，與層流速度型相比，入口條件為紊流速度型和特征邊界型得

5、到的結(jié)果能夠吻合。然后，本文研究了有物體存在的空腔流動的剪切層自持性振蕩特性，并與無外掛物存在時的情況進行類比，壓強與時間的變化曲線、空腔底部流向的靜壓系數(shù)分布以及空腔后壁面法向的靜壓系數(shù)分布表明，空腔內(nèi)的壓力脈動與外掛物所處位置是相關(guān)的。最后，本文討論了腔內(nèi)有不同長度平板存在的空腔流動的剪切層自持性振蕩，腔內(nèi)平板將剪切層與空腔隔開，減弱了剪切層振蕩與空腔壓強脈動的耦合作用，隨著腔內(nèi)平板長度的增加，剪切層振蕩的程度有所減?。换谔卣鏖L度