本文選題：流體動(dòng)力學(xué) + 多相流��； 參考：《中國(guó)工程物理研究院》2016年博士論文

【摘要】：多相反應(yīng)流的數(shù)值模擬在許多應(yīng)用研究領(lǐng)域發(fā)揮著重要的作用,例如在武器彈藥設(shè)計(jì)與毀傷評(píng)估,重大爆炸災(zāi)害的防災(zāi)減災(zāi),高速推進(jìn)器技術(shù)以及天體物理中的超新星爆發(fā)等國(guó)防,民用和科學(xué)研究中都是不可或缺的研究手段.本文主要研究多相反應(yīng)流的數(shù)值模擬方法.多相反應(yīng)流指含有化學(xué)反應(yīng)的多相流,不同相對(duì)應(yīng)的狀態(tài)參數(shù)不同.本文首先從數(shù)值格式,邊界處理,網(wǎng)格技術(shù)和邊界處理等多個(gè)方面,綜述了微分方程數(shù)值模擬與爆轟史上的重要進(jìn)展和發(fā)展變化,特別是多相反應(yīng)流的歷史發(fā)展.對(duì)爆轟的經(jīng)典理論,如Chapman-Jouguet(CJ)模型和Zeldovich-von Neumann-Doring(ZN D)模型作了比較詳細(xì)的論述,介紹了爆轟過(guò)程從瞬間反應(yīng)到區(qū)分前導(dǎo)激波和反應(yīng)區(qū),再到二維多波相互作用而產(chǎn)生的胞格現(xiàn)象的整個(gè)過(guò)程.對(duì)于理想流體,本文利用Godunov格式,結(jié)合Harten-Lax-Leer-Contact (HLLC)類(lèi)型的Riemann求解器,構(gòu)造了相應(yīng)的數(shù)值方法,并應(yīng)用到一維和二維理想爆轟問(wèn)題的數(shù)值模擬,包括穩(wěn)態(tài)的與非穩(wěn)態(tài)爆轟問(wèn)題.數(shù)值模擬結(jié)果表明,這種算法魯棒性好,能夠有效捕捉爆轟波的結(jié)構(gòu)特征.同時(shí)研究了高精度的移動(dòng)網(wǎng)格數(shù)值方法在爆轟流體中的應(yīng)用,比較了他們與廣義黎曼問(wèn)題(GRP)方法在數(shù)值精度方面的差異.對(duì)于非理想流體,本文主要考慮了三種不同類(lèi)型的狀態(tài)方程：剛性狀態(tài)方程,Cochran-Chan(CC)狀態(tài)方程和Jones-Wilkins-Lee(JWL)狀態(tài)方程.針對(duì)這些非理想反應(yīng)流,由于歐拉方法數(shù)值模擬中,會(huì)出現(xiàn)不同相的混合單元,本文發(fā)展了一種物理量重構(gòu)法.應(yīng)用單元格內(nèi)不同物質(zhì)的物理量之間的關(guān)系,如混合密度與各自密度之間的關(guān)系,以及混合內(nèi)能與各物質(zhì)內(nèi)能之間的關(guān)系,根據(jù)熱力學(xué)平衡條件,建立了多個(gè)未知變量的方程組.通過(guò)求解這個(gè)方程組,重構(gòu)出混合單元內(nèi)的各相物理量,結(jié)合前述的HLLC求解器,計(jì)算出單元界面的數(shù)值通量.在求解方程組時(shí),由于相應(yīng)的代數(shù)方程的多解,給數(shù)值方法帶來(lái)一定的困難,本文提出一種“移動(dòng)跟蹤法”,能夠快速得到具有物理意義的解.本文將上述算法應(yīng)用到多相流的數(shù)值模擬中,給出了大量的一維和二維算例.數(shù)值結(jié)果表明,這種算法既能清晰地捕捉一維和二維的爆轟波的結(jié)構(gòu),又能比較準(zhǔn)確地捕獲多波的相互作用,得到的胞格邊界清楚,排列有序,三波點(diǎn)的特征明顯.這些都驗(yàn)證了該算法的有效性和可靠性.作為論文的結(jié)尾,總結(jié)討論了這些算法的特點(diǎn),并展望了下一步的工作可能遇到的困難和問(wèn)題.
[Abstract]:Numerical simulation of multiphase reaction flows plays an important role in many applied research areas, such as weapons and ammunition design and damage assessment, disaster prevention and mitigation of major explosive disasters,High-speed thruster technology and supernova explosion in astrophysics are indispensable for national defense, civil and scientific research.In this paper, the numerical simulation method of multiphase reaction flow is studied.Multiphase reaction flow refers to a multiphase flow containing chemical reactions with different state parameters.In this paper, the important progress and changes in numerical simulation of differential equations and detonation history, especially the history of multiphase reaction flow, are reviewed in terms of numerical scheme, boundary treatment, grid technology and boundary treatment.The classical theories of detonation, such as Chapman-Jouguett CJ model and Zeldovich-von Neumann-Doring(ZN Dmodel, are discussed in detail, and the whole process of detonation process from instantaneous reaction to distinguishing leading shock and reaction region to two-dimensional multiwave interaction is introduced.For the ideal fluid, the corresponding numerical method is constructed by using the Godunov scheme and the Riemann solver of Harten-Lax-Leer-Contact / HLLC type, and applied to the numerical simulation of one-dimensional and two-dimensional ideal detonation problems, including steady and unsteady detonation problems.The numerical simulation results show that the proposed algorithm is robust and can capture the structural characteristics of detonation waves effectively.At the same time, the application of high precision moving grid numerical method in detonation fluid is studied, and the difference between them and the generalized Riemannian problem (GRP) method in numerical accuracy is compared.For nonideal fluids, three different types of state equations are considered in this paper: the rigid state equation (Cochran-Chann CCS) and the Jones-Wilkins-Leewe (JWL) equation of state.For these non-ideal reaction flows, due to the existence of mixed elements with different phases in the numerical simulation of Euler's method, a physical quantity reconstruction method is developed in this paper.Based on the thermodynamic equilibrium conditions, the relationship between the physical quantities of different substances in a cell, such as the relation between the mixing density and the respective density, and the relationship between the mixed internal energy and the internal energy of each substance, is applied.The equations of several unknown variables are established.By solving the equations, the physical quantities of each phase in the mixed element are reconstructed, and the numerical flux of the interface of the unit is calculated by combining the HLLC solver mentioned above.In solving the equations, it is difficult to solve the numerical method because of the multiple solutions of the corresponding algebraic equations. In this paper, a "moving tracking method" is proposed, which can quickly obtain the solution with physical significance.In this paper, the algorithm is applied to the numerical simulation of multiphase flow, and a large number of one-dimensional and two-dimensional numerical examples are given.The numerical results show that this algorithm can capture the structure of one-dimensional and two-dimensional detonation waves clearly and accurately capture the interaction of multiple waves. The results show that the lattice boundary is clear, the arrangement is ordered, and the characteristics of the three wave points are obvious.These results demonstrate the validity and reliability of the algorithm.At the end of the paper, the characteristics of these algorithms are summarized and discussed, and the difficulties and problems that may be encountered in the next work are prospected.
【學(xué)位授予單位】：中國(guó)工程物理研究院
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2016
【分類(lèi)號(hào)】：O241.82

【相似文獻(xiàn)】

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

1 崔虎雄;于興娜;成天濤;徐永福;房文;陳建民;;SO_2與鐵氧化物多相反應(yīng)對(duì)粒子光學(xué)特性的影響[J];過(guò)程工程學(xué)報(bào);2009年S2期

2 ;[J];;年期

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

1 翁春生;金志明;袁亞雄;;高含能顆粒床多相反應(yīng)流動(dòng)的研究[A];面向21世紀(jì)的科技進(jìn)步與社會(huì)經(jīng)濟(jì)發(fā)展（下冊(cè)）[C];1999年

2 崔虎雄;于興娜;成天濤;徐永福;房文;陳建民;吳蓓莉;;SO_2與鐵氧化物多相反應(yīng)對(duì)粒子光學(xué)特性影響的實(shí)驗(yàn)?zāi)M研究[A];中國(guó)顆粒學(xué)會(huì)第六屆學(xué)術(shù)年會(huì)暨海峽兩岸顆粒技術(shù)研討會(huì)論文集（下）[C];2008年

3 張佳晨;陳海涵;陳建民;;堿性分子催化氧化鐵與羰基硫多相反應(yīng)研究[A];第五屆全國(guó)環(huán)境化學(xué)大會(huì)摘要集[C];2009年

4 魏進(jìn)超;;燒結(jié)過(guò)程中二惡英的形成及抑制技術(shù)[A];第5屆中國(guó)金屬學(xué)會(huì)青年學(xué)術(shù)年會(huì)論文集[C];2010年

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

1 黃苗榮 毛恭忠;揮灑環(huán)保新篇章[N];中國(guó)有色金屬報(bào);2003年

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

1 曾現(xiàn)洋;多相反應(yīng)流的數(shù)值模擬方法[D];中國(guó)工程物理研究院;2016年

2 周媛;結(jié)構(gòu)化多相反應(yīng)器單通道及床層尺度流動(dòng)與混合特性研究[D];北京化工大學(xué);2011年

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

1 楊東輝;數(shù)學(xué)建模在多相反應(yīng)動(dòng)力學(xué)上的應(yīng)用[D];河北師范大學(xué);2005年

2 陳朝暉;多相反應(yīng)色譜過(guò)程的數(shù)學(xué)模擬[D];浙江大學(xué);2006年

3 陳亞嵐;結(jié)構(gòu)化多相反應(yīng)器內(nèi)液相流動(dòng)及軸向擴(kuò)散性能的實(shí)驗(yàn)研究[D];北京化工大學(xué);2011年

4 杜建平;反應(yīng)物粒度對(duì)多相化學(xué)反應(yīng)影響的研究[D];太原理工大學(xué);2005年

，

本文編號(hào)：1769316