發(fā)布時(shí)間：2019-04-02 18:33

【摘要】：Navier-Stokes方程(簡稱N-S方程)具有高度的非線性特性,幾乎無法獲得精確解,因此采用數(shù)值方法來求解N-S方程則成為了研究的熱點(diǎn)。目前主要的數(shù)值方法有:有限元法、有限體積法和有限差分法等,但這些方法因存在網(wǎng)格易在迭代計(jì)算過程中出現(xiàn)網(wǎng)格畸變、計(jì)算不穩(wěn)定、數(shù)值振蕩等問題,從而限制了其計(jì)算精度以及效率的提高,因此為N-S方程的求解提供一種高效穩(wěn)定的數(shù)值計(jì)算方法是十分必要的。無網(wǎng)格Galerkin(Element-free Galerkin,EFG)法是一種新型且僅需要節(jié)點(diǎn)信息的計(jì)算方法,它可構(gòu)造出高階的場函數(shù),能擺脫不連續(xù)性對(duì)問題的約束,具有更高的光滑性等優(yōu)勢(shì),因此本文采用分離式的EFG法來求解N-S方程,并對(duì)其工程應(yīng)用進(jìn)行研究,其主要內(nèi)容如下:(1)通過采用罰因子法來處理N-S方程中的壓力項(xiàng),得到了二維N-S方程的EFG法的離散矩陣形式,通過算例驗(yàn)證了所提算法的可行性,探討了高斯點(diǎn)個(gè)數(shù)、縮放系數(shù)以及流體雷諾數(shù)對(duì)計(jì)算結(jié)果的影響規(guī)律,并與有限元結(jié)果進(jìn)行了對(duì)比分析,所得結(jié)果顯示EFG法具有較好的收斂速度和計(jì)算精度。(2)探討了EFG法中不同本質(zhì)邊界條件處理方式即拉格朗日乘子法(LMM)、罰函數(shù)法(PM)、耦合有限元法(FEM—EFG)對(duì)N-S方程計(jì)算結(jié)果的影響,并與有限元結(jié)果進(jìn)行比較可知,LMM法的計(jì)算結(jié)果與有限元解的相對(duì)誤差最小,而PM法結(jié)果的相對(duì)誤差最大,FEM—EFG法的結(jié)果與LMM法的結(jié)果十分接近,同時(shí)探討了PM法施加邊界條件時(shí)罰因子的選取對(duì)計(jì)算精度的影響,所得結(jié)果顯示調(diào)整罰因子的值可有效提高模擬計(jì)算精度。(3)針對(duì)翼型和彈丸流場中對(duì)流項(xiàng)(慣性項(xiàng))會(huì)引起數(shù)值不穩(wěn)定問題,將EFG法僅需要節(jié)點(diǎn)信息的優(yōu)勢(shì)與迎風(fēng)穩(wěn)定化算法相結(jié)合,構(gòu)造了一種計(jì)算點(diǎn)的偏心支持域,并將其應(yīng)用到NACA0012翼型和彈丸流場問題的模擬計(jì)算中,有效地避免了數(shù)值振蕩現(xiàn)象,并提高了計(jì)算精度。本文采用EFG法和采用偏心支持域的EFG法對(duì)二維N-S方程進(jìn)行了求解和分析,得到的結(jié)果表明EFG法具有計(jì)算穩(wěn)定、精度高、收斂快等特點(diǎn),所得結(jié)論對(duì)流體研究及其工程應(yīng)用具有一定的理論參考和應(yīng)用價(jià)值。
[Abstract]:The Navier-Stokes equation (abbreviated as Nus equation) has a high degree of non-linearity and can hardly obtain the exact solution. Therefore, the numerical method for solving the Nus equation has become a hot topic in the field of research. At present, the main numerical methods are: finite element method, finite volume method and finite difference method, but these methods have many problems such as grid distortion, computational instability, numerical oscillation and so on. Therefore, it is necessary to provide an efficient and stable numerical calculation method for solving the Navier-Stokes equation, which limits the accuracy and efficiency of the calculation. Meshless Galerkin (Element-free Galerkin,EFG) method is a new computing method which only needs node information. It can construct higher-order field functions, get rid of the constraints of discontinuity to the problem, and has the advantages of higher smoothness and so on. The main contents of this paper are as follows: (1) the penalty factor method is used to deal with the pressure term in the Navier-Stokes equation, and its engineering application is studied in this paper. The main contents are as follows: (1) the penalty factor method is used to deal with the pressure term in the Nus equation. The discrete matrix form of the EFG method for the two-dimensional Navier-Stokes equation is obtained. The feasibility of the proposed algorithm is verified by an example. The effects of the number of Gao Si points, the scaling coefficient and the Reynolds number of the fluid on the calculated results are discussed. Compared with the finite element results, the results show that the EFG method has good convergence rate and calculation precision. (2) the different methods of dealing with the essential boundary conditions in the EFG method, that is, Lagrangian multiplier method (LMM), are discussed. The effect of penalty function method (PM), coupled finite element method (FEM-EFG) on the calculation results of the Nus equation is compared with that of the finite element method. It is shown that the relative error between the LMM method and the finite element solution is the smallest. The relative error of PM method is the largest, and the result of FEM-EFG method is very close to that of LMM method. At the same time, the influence of penalty factor on the accuracy of calculation is discussed when PM method is applied with boundary conditions. The results show that adjusting the penalty factor can effectively improve the simulation accuracy. (3) the convective term (inertia term) in the airfoil and projectile flow field will cause the numerical instability problem. Combining the advantage of EFG method which only needs node information and upwind stabilization algorithm, an eccentricity support domain for calculating points is constructed and applied to the simulation of NACA0012 airfoil and projectile flow field, which effectively avoids the phenomenon of numerical oscillation. The calculation precision is improved. In this paper, the EFG method and the EFG method with eccentricity support domain are used to solve and analyze the two-dimensional Navier-Stokes equation. The results show that the EFG method has the characteristics of stable calculation, high precision and fast convergence. The conclusion has certain theoretical reference and application value for fluid research and its engineering application.
【學(xué)位授予單位】：湘潭大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2016
【分類號(hào)】：O241.82

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 劉紹球 ,李顯霖;有限元法淺說[J];系統(tǒng)工程與電子技術(shù);1981年02期

2 王月f;《水資源中有限元法》——第四屆國際水資源應(yīng)用有限元法會(huì)議會(huì)議錄介紹[J];地質(zhì)科技情報(bào);1983年04期

3 鮑四元;;本科生《有限元法》課程教學(xué)的現(xiàn)狀與若干改革研究[J];課程教育研究;2014年03期

4 石艷華;侯錄;席俊明;;有限元法在口腔醫(yī)學(xué)中的應(yīng)用[J];今日科苑;2007年08期

5 李昱君;;從語境的視角看有限元法的早期發(fā)展[J];自然辯證法通訊;2008年03期

6 魯建霞;茍惠芳;;有限元法的基本思想與發(fā)展過程[J];機(jī)械管理開發(fā);2009年02期

7 葛藤;;《有限元法》本科教學(xué)改革探討[J];科技信息;2010年13期

8 張泉樂;薛荔;;非線性問題的有限元法[J];武漢理工大學(xué)學(xué)報(bào);2010年15期

9 王小榮;林順洪;丁劍平;陳小亮;;有限元法課程的教學(xué)改革探討[J];重慶科技學(xué)院學(xué)報(bào)(社會(huì)科學(xué)版);2011年06期

10 徐建全;陳銘年;林大同;;《有限元法》本科課程綜合改革與實(shí)踐[J];機(jī)電技術(shù);2011年06期

相關(guān)會(huì)議論文 前10條

1 唐旭海;鄭超;張建海;;基于虛結(jié)點(diǎn)的多邊形有限元法[A];中國力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)'2009論文摘要集[C];2009年

2 趙俊峰;周慎杰;王炳雷;王錫平;;應(yīng)變梯度微梁的有限元法[A];力學(xué)與工程應(yīng)用（第十三卷）[C];2010年

3 Liu Bin;Xu Ran;Zhang Yihui;Fang Daining;;原子有限元法在多物理場研究中的新進(jìn)展(英文)[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

4 丁勝勇;邵國建;;基于多邊形有限元法的混凝土細(xì)觀模擬[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

5 王耀;黑寶平;楊在林;;用于非均勻介質(zhì)瞬態(tài)分析的動(dòng)態(tài)非均勻有限元法[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

6 金峰;方修君;;擴(kuò)展有限元法及與其他數(shù)值方法的聯(lián)系[A];第17屆全國結(jié)構(gòu)工程學(xué)術(shù)會(huì)議論文集（第Ⅰ冊(cè)）[C];2008年

7 鄧嬌;王宏志;傅向榮;鄧小環(huán);宋孟燕;葛麗娜;田歌;;單自由度時(shí)空有限元法[A];北京力學(xué)會(huì)第19屆學(xué)術(shù)年會(huì)論文集[C];2013年

8 徐春暉;李明瑞;;可離散有限元法中應(yīng)變度量的選擇[A];第五屆全國計(jì)算爆炸力學(xué)會(huì)議論文摘要[C];2012年

9 張謐;徐春輝;;由于連續(xù)體破壞問題研究的可離散有限元法[A];北京力學(xué)會(huì)第19屆學(xué)術(shù)年會(huì)論文集[C];2013年

10 張俊峰;李興國;婁國偉;;有限元法在單片電路耦合設(shè)計(jì)中的應(yīng)用[A];2007年全國微波毫米波會(huì)議論文集（下冊(cè)）[C];2007年

相關(guān)重要報(bào)紙文章 前1條

1 實(shí)習(xí)記者 南秀淵;有限元法研究的崢嶸歲月[N];新清華;2014年

相關(guān)博士學(xué)位論文 前10條

1 賈紅剛;復(fù)合材料斷裂分析的擴(kuò)展有限元法研究[D];西北工業(yè)大學(xué);2015年

2 張雙文;時(shí)域有限元法及其截?cái)噙吔鐥l件的研究[D];西南交通大學(xué);2008年

3 付振東;非均質(zhì)飽和多孔介質(zhì)準(zhǔn)靜態(tài)行為分析的耦合多尺度及耦合升尺度有限元法[D];大連理工大學(xué);2010年

4 林志良;比例邊界有限元法及快速多極子邊界元法的研究與應(yīng)用[D];上海交通大學(xué);2010年

5 卜令方;增強(qiáng)有限元法及其在巖土工程非連續(xù)變形分析中的應(yīng)用[D];浙江大學(xué);2014年

6 林峰;瞬變電磁場中時(shí)空有限元法的研究[D];沈陽工業(yè)大學(xué);2005年

7 曹鳳帥;比例邊界有限元法在勢(shì)流理論中的應(yīng)用[D];大連理工大學(xué);2009年

8 鄭恩希;非多項(xiàng)式最小二乘有限元法在幾種散射問題中的應(yīng)用[D];吉林大學(xué);2012年

9 王輝;基于徑向基函數(shù)的無網(wǎng)格數(shù)值方法及雜交Trefftz有限元法[D];天津大學(xué);2007年

10 李霄琳;非均質(zhì)材料的光滑多尺度有限元法研究[D];吉林大學(xué);2015年

相關(guān)碩士學(xué)位論文 前10條

1 李昱君;語境視角下的有限元法發(fā)展史[D];山西大學(xué);2008年

2 蘇偉江;基于有限元法的漁船碰撞仿真研究[D];大連海事大學(xué);2015年

3 王眾;三維隨機(jī)水質(zhì)的Galerkin有限元法及應(yīng)用[D];長春工業(yè)大學(xué);2016年

4 俞晨洋;基于有限元法的天線優(yōu)化設(shè)計(jì)及運(yùn)動(dòng)目標(biāo)散射特性研究[D];華東師范大學(xué);2016年

5 羅丹;計(jì)及非線性時(shí)同步發(fā)電機(jī)穩(wěn)態(tài)參數(shù)隨運(yùn)行條件變化的規(guī)律[D];華北電力大學(xué)(北京);2016年

6 王佳圓;基于比例邊界有限元法的彈塑性動(dòng)力響應(yīng)分析[D];浙江大學(xué);2016年

7 易若愚;二維N-S方程的EFG數(shù)值模擬及其應(yīng)用研究[D];湘潭大學(xué);2016年

8 李永利;邊光滑有限元法在汽車電磁兼容計(jì)算中的研究及應(yīng)用[D];湖南大學(xué);2016年

9 樊現(xiàn)行;光滑有限元法理論及算法研究[D];山東大學(xué);2013年

10 萬德濤;二維光滑擴(kuò)展有限元法研究[D];湖南大學(xué);2014年

，

本文編號(hào)：2452804