可壓縮多項(xiàng)流問題的數(shù)值研究及應(yīng)用
發(fā)布時(shí)間：2018-08-29 19:11
【摘要】：在可壓縮多相流的數(shù)值模擬研究中,帶自由界面的Cut-Cell方法已經(jīng)由Chang等人(2013)[1]所發(fā)展,我們使用基于此方法所發(fā)展起來的MuSiC~+程序,數(shù)值研究了包含氣體-氣體,氣體-液體以及氣體-固體介質(zhì)的Richtmyer-Meshkov(RM)界面不穩(wěn)定,驗(yàn)證了氣-液介質(zhì)小擾動(dòng)振幅下的界面增長率在前期符合Yang等人(1994)[2]提出的線性化模型,后期符合ZhangSohn(1996,1997)[3,4]提出的非線性模型:模擬了有著廣泛應(yīng)用前景的微射流問題,通過修改Peters等人2013[5]的模型,得到了一個(gè)適用于由激波誘導(dǎo)的微射流演化過程中最大速度與初始接觸角,激波強(qiáng)度的半經(jīng)驗(yàn)關(guān)系式;通過數(shù)值擬合獲得了一個(gè)適用于高M(jìn)ach數(shù)情形的射流最大速度的模型關(guān)系式�？紤]到在MuSiC~+程序的實(shí)現(xiàn)中,Chang等人通過Level Set函數(shù)值來演化物質(zhì)界面,并引入小網(wǎng)格單元的切分和融合,以此來刻畫自由界面,在界面演化的過程中通過對(duì)界面速度的傳播及插值來完成對(duì)界面附近網(wǎng)格點(diǎn)的速度的賦值,這樣的處理在界面發(fā)生嚴(yán)重變形尤其是在有界面拓?fù)浣Y(jié)構(gòu)變化發(fā)生的一些問題時(shí),對(duì)一個(gè)網(wǎng)格被界面多次切分的現(xiàn)象還需考慮。通過繼承MuSiC~+程序中的網(wǎng)格切分的思想,參考Hu等人(2006)提出的一種守恒界面方法[6]的一些界面處理過程,以及引入Ghost fluid的方法(Fedkiw 1999,2001,2002[7-9];Liu 2003,2005[10,11]),在界面處我們提出 了一種 real-ghost mixing 的切分網(wǎng)格狀態(tài)的處理方法,這種方法很好的處理了在含有物質(zhì)界面拓?fù)浣Y(jié)構(gòu)變化的一些多相流問題的演化模擬。新發(fā)展的數(shù)值方法(我們稱之為一種基于切割網(wǎng)格和ghost fluid的可壓縮多相流求解方法,簡(jiǎn)記為CCGF)通過對(duì)一維激波管問題(Air-SF6,Air-Helium,Water-Air)的計(jì)算,以及與解析解的比較,能夠看出是準(zhǔn)確的;通過對(duì)大量二維經(jīng)典問題(Air-SF6以及Air-Helium的Richtmyer-Meshkov的界面不穩(wěn)定問題,Air-Helium以及Air-R22的激波氣泡相互作用問題,水下氣泡在強(qiáng)激波作用下的破碎問題,水下爆炸問題,空氣中激波撞擊液柱及雙液柱問題)的模擬研究,以及比較相關(guān)實(shí)驗(yàn)結(jié)果和之前的數(shù)值研究結(jié)果,能夠看出目前的數(shù)值方法是可信的。使用MuSiC~+程序與CCGF程序,數(shù)值研究了氣體-氣體介質(zhì)的RM不穩(wěn)定性,給出了界面演化的分布圖,界面增長率的變化曲線圖,分析了引起界面增長率震蕩的原因,并和已有的理論結(jié)果(Yang等人線性化模型結(jié)果,ZhangSohn的非線性模型(ZS),Sadot等人非線性模型(SEA)[12],DimonteRamaprabhu的非線性模型(DR)[13])以及數(shù)值結(jié)果(Holmes等人[14,15],茅德康等人[16])做比較。通過考慮不同初始擾動(dòng)振幅下的界面增長,進(jìn)一步驗(yàn)證在小擾動(dòng)情形下三個(gè)非線性模型與數(shù)值結(jié)果是一致的;對(duì)于大擾動(dòng)情形,DR模型與數(shù)值結(jié)果有很好的一致性,ZS模型過低的估計(jì)了界面增長,SEA模型則高估界面的擾動(dòng)增長。為了研究更為復(fù)雜的流動(dòng)問題(多介質(zhì)液滴的相互碰撞和穿透問題),我們?cè)噲D同時(shí)使用基于MuSiC~+程序的數(shù)值方法和基于CCGF的數(shù)值方法來解決這類問題,考慮在多界面情況下兩種方法共同使用分別處理不同界面的流動(dòng),甚至也考慮在兩相流的不同演化階段采用不同的方法來模擬研究,目前我們已經(jīng)能夠?qū)崿F(xiàn)在兩相流問題中兩種不同方法的相互轉(zhuǎn)化使用,通過對(duì)水下激波作用下氣泡破碎問題的模擬,能夠看出兩種方法的相互轉(zhuǎn)化來模擬一個(gè)問題是成功的。兩種方法共存模擬多界面問題是未來的工作。
[Abstract]:In the numerical simulation of compressible multiphase flow, the Cut-Cell method with free interface has been developed by Chang et al. (2013)[1]. Using the MuSiC~+ program developed based on this method, the Richtmyer-Meshkov (RM) interfacial instability including gas-gas, gas-liquid and gas-solid media has been numerically studied. The interfacial growth rate under small disturbance amplitude of gas-liquid medium conforms to the linearized model proposed by Yang et al. (1994) [2] in the early stage and to the nonlinear model proposed by Zhang Sohn (1996, 1997) [3, 4] in the later stage. The problem of micro-jet with wide application prospects is simulated. By modifying the model of Peters et al. [5] in 2013, a model suitable for shock induced flow is obtained. A semi-empirical relationship between the maximum velocity and the initial contact angle and the shock intensity in the evolution process of a guided micro-jet is obtained by numerical fitting. Considering the realization of MuSiC~+ program, Chang et al. used Level Set function to evolve the material interface, and introduced a new model to describe the maximum velocity of a guided micro-jet. In the process of interface evolution, the velocity of grid points near the interface is assigned by the propagation and interpolation of the velocity of the interface. In this way, the interface is severely deformed, especially when there are some problems of interface topological structure changes. By inheriting the idea of grid segmentation in MuSiC~+ program, referring to some interface processing procedures of Hu et al. (2006) and introducing Ghost fluid (Fedkiw 1999, 2001, 2002 [7-9]; Liu 2003, 2005 [10, 11]), we propose a new method at the interface. The real-ghost mixing method, which deals well with the evolutionary simulation of some multiphase flow problems involving changes in the topology of the material interface, is presented. Comparing with the analytical solutions, it can be seen that the calculation of the one-dimensional shock tube problem (Air-SF6, Air-Helium, Water-Air) is accurate; through a large number of two-dimensional classical problems (Air-SF6 and Air-Helium's Richtmyer-Meshkov interface instability problem, Air-Helium and Air-R22's shock bubble interaction problem, underwater bubble in strong excitation The present numerical method is believable by comparing the experimental results with the previous numerical results. The RM instability of gas-gas medium is numerically studied by using the MuSiC~+ and CGF codes. In this paper, the distribution of interface evolution and the change curve of interface growth rate are given, and the reasons causing the oscillation of interface growth rate are analyzed. The results are compared with the existing theoretical results (Yang et al. linear model, Zhang Sohn's nonlinear model (ZS), Sadot et al.'s nonlinear model (SEA) [12], Dimonte Ramaprab's nonlinear model (DR) [13]) and numbers. Numerical results (Holmes et al. [14,15], Mao Dekang et al. [16]) are compared. Considering the interface growth under different initial disturbance amplitudes, it is further verified that the three nonlinear models are consistent with the numerical results in the case of small disturbance; for the case of large disturbance, the DR model is in good agreement with the numerical results, and the ZS model underestimates the interface. In order to study more complex flow problems (collision and penetration of droplets in multi-media), we attempt to solve these problems simultaneously by using numerical methods based on MuSiC~+ program and CGF. In dealing with the flow at different interfaces, different methods are even considered to simulate the two-phase flow at different stages of evolution. At present, we have been able to realize the mutual conversion of two different methods in the two-phase flow problem. To simulate a problem is successful. The two method of coexistence and Simulation of multiple interface problems is the future work.
【學(xué)位授予單位】：中國科學(xué)技術(shù)大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2016
【分類號(hào)】：O359
本文編號(hào)：2212138
[Abstract]:In the numerical simulation of compressible multiphase flow, the Cut-Cell method with free interface has been developed by Chang et al. (2013)[1]. Using the MuSiC~+ program developed based on this method, the Richtmyer-Meshkov (RM) interfacial instability including gas-gas, gas-liquid and gas-solid media has been numerically studied. The interfacial growth rate under small disturbance amplitude of gas-liquid medium conforms to the linearized model proposed by Yang et al. (1994) [2] in the early stage and to the nonlinear model proposed by Zhang Sohn (1996, 1997) [3, 4] in the later stage. The problem of micro-jet with wide application prospects is simulated. By modifying the model of Peters et al. [5] in 2013, a model suitable for shock induced flow is obtained. A semi-empirical relationship between the maximum velocity and the initial contact angle and the shock intensity in the evolution process of a guided micro-jet is obtained by numerical fitting. Considering the realization of MuSiC~+ program, Chang et al. used Level Set function to evolve the material interface, and introduced a new model to describe the maximum velocity of a guided micro-jet. In the process of interface evolution, the velocity of grid points near the interface is assigned by the propagation and interpolation of the velocity of the interface. In this way, the interface is severely deformed, especially when there are some problems of interface topological structure changes. By inheriting the idea of grid segmentation in MuSiC~+ program, referring to some interface processing procedures of Hu et al. (2006) and introducing Ghost fluid (Fedkiw 1999, 2001, 2002 [7-9]; Liu 2003, 2005 [10, 11]), we propose a new method at the interface. The real-ghost mixing method, which deals well with the evolutionary simulation of some multiphase flow problems involving changes in the topology of the material interface, is presented. Comparing with the analytical solutions, it can be seen that the calculation of the one-dimensional shock tube problem (Air-SF6, Air-Helium, Water-Air) is accurate; through a large number of two-dimensional classical problems (Air-SF6 and Air-Helium's Richtmyer-Meshkov interface instability problem, Air-Helium and Air-R22's shock bubble interaction problem, underwater bubble in strong excitation The present numerical method is believable by comparing the experimental results with the previous numerical results. The RM instability of gas-gas medium is numerically studied by using the MuSiC~+ and CGF codes. In this paper, the distribution of interface evolution and the change curve of interface growth rate are given, and the reasons causing the oscillation of interface growth rate are analyzed. The results are compared with the existing theoretical results (Yang et al. linear model, Zhang Sohn's nonlinear model (ZS), Sadot et al.'s nonlinear model (SEA) [12], Dimonte Ramaprab's nonlinear model (DR) [13]) and numbers. Numerical results (Holmes et al. [14,15], Mao Dekang et al. [16]) are compared. Considering the interface growth under different initial disturbance amplitudes, it is further verified that the three nonlinear models are consistent with the numerical results in the case of small disturbance; for the case of large disturbance, the DR model is in good agreement with the numerical results, and the ZS model underestimates the interface. In order to study more complex flow problems (collision and penetration of droplets in multi-media), we attempt to solve these problems simultaneously by using numerical methods based on MuSiC~+ program and CGF. In dealing with the flow at different interfaces, different methods are even considered to simulate the two-phase flow at different stages of evolution. At present, we have been able to realize the mutual conversion of two different methods in the two-phase flow problem. To simulate a problem is successful. The two method of coexistence and Simulation of multiple interface problems is the future work.
【學(xué)位授予單位】：中國科學(xué)技術(shù)大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2016
【分類號(hào)】：O359
【參考文獻(xiàn)】
相關(guān)期刊論文 前10條
1 劉軍;馮其京;周海兵;;柱面內(nèi)爆驅(qū)動(dòng)金屬界面不穩(wěn)定性的數(shù)值模擬研究[J];物理學(xué)報(bào);2014年15期
2 劉金宏;鄒立勇;曹仁義;廖深飛;王彥平;;繞射激波和反射激波作用下N_2/SF_6界面R-M不穩(wěn)定性實(shí)驗(yàn)研究[J];力學(xué)學(xué)報(bào);2014年03期
3 羅喜勝;王顯圣;陳模軍;翟志剛;;可控肥皂膜氣柱界面與激波相互作用的實(shí)驗(yàn)研究（英文）[J];實(shí)驗(yàn)流體力學(xué);2014年02期
4 施紅輝;肖毅;杜凱;吳宇;;用垂直激波管研究多層流體界面上的Richtmyer-Meshkov不穩(wěn)定性[J];中國科學(xué)技術(shù)大學(xué)學(xué)報(bào);2013年09期
5 劉金宏;譚多望;柏勁松;黃文斌;鄒立勇;張旭;;激波管實(shí)驗(yàn)研究非均勻流場(chǎng)R-M不穩(wěn)定性[J];實(shí)驗(yàn)力學(xué);2012年02期
6 施紅輝;杜凱;王超;章利特;賈會(huì)霞;董若凌;;不同密度梯度的多層流體界面上的Richtmyer-Meshkov不穩(wěn)定性研究[J];實(shí)驗(yàn)流體力學(xué);2011年05期
7 M.A.ULLAH;高文斌;茅德康;;Numerical simulations of Richtmyer-Meshkov instabilities using conservative front-tracking method[J];Applied Mathematics and Mechanics(English Edition);2011年01期
8 何長江;周海兵;杭義洪;;爆轟驅(qū)動(dòng)金屬鋁界面不穩(wěn)定性的數(shù)值分析[J];中國科學(xué)(G輯:物理學(xué) 力學(xué) 天文學(xué));2009年09期
9 ;Numerical investigation of Richtmyer-Meshkov instability driven by cylindrical shocks[J];Acta Mechanica Sinica;2006年01期
10 ;Two-phase Flow Patterns in a Square Mini-channel[J];Journal of Thermal Science;2004年02期
，本文編號(hào)：2212138
本文鏈接：http://sikaile.net/shoufeilunwen/jckxbs/2212138.html
最近更新
教材專著
