伴随凝结的无粘可压缩流动研究
【摘要】:伴随凝结的可压缩流动是一个非常普遍的流动现象,在传统工业如汽轮机、风洞及新兴技术如超紫外光刻技术等中应用广泛。气体转变为液体的凝结过程,伴随着凝结潜热的释放,急剧改变了流场的热力学性质。流场的改变直接影响了设备的运行状态,如汽轮机运行效率降低、风洞偏离设计参数等。本文研究凝结的自激发自抑制的基本过程对流动的影响程度,主要包括两方面内容:一是基本流动如一维管流、二维喷管流、绕角流中的凝结,着重于讨论、评估凝结引起的波现象;二是风洞两个典型应用中的凝结,分别是燃烧加热风洞中水蒸气的凝结,常规高超风洞中氮气的凝结,着重于分析凝结对流场参数的影响。这两方面均有较强的基础研究背景和工业价值。
本文首先从理论分析凝结的气动效应入手,讨论凝结启动过程中的非定常波系。凝结过程中包含了潜热释放、气体消耗、液滴增加三个子过程。通常研究中仅仅考虑加热过程,这一近似极大地简化了研究的复杂度,但同时也带来一定误差。我们对这一假设重新评估,发现潜热释放和气体消耗两个过程共同决定了非定常波的表现。在评估中,分别使用一维加热模型和减质模型讨论了潜热释放和气体消耗两过程中产生的非定常波的波强。在分析了典型物质水、壬烷凝结时的热量添加比率与质量减少比率之间的关系后,本文发现气体消耗引起的非定常波不容忽视。进一步地,通过综合考虑加热与减质两个因素,给出了关于凝结引起的非定常波的波强的合理评估。结果发现,综合考虑热质影响与仅考虑热影响所得到的非定常波的波强相差较大,对于水蒸气差别约为15%~20%,对于壬烷则约为30%~60%。当凝结气体的速度低于来流流速时,非定常波的差别将更大。
从气动理论的分析可以总结出,水蒸气的凝结中加热过程占据主要地位。为验证上述分析,我们首先考虑了喷管中凝结的两个典型算例及激波管中凝结的启动过程,发现其中的非定常波系与加热分析所得到的三个波系一一吻合。在激波管中凝结的研究中,更进一步地,我们讨论了无限长激波管中的凝结现象,引入了启动阶段、振荡阶段、渐近阶段三个阶段的概念,并分析了不同初始参数情况下启动阶段和振荡阶段的表现。在渐近阶段,随着流场不断演化,可以发现一个压力、温度的平台。从平衡凝结概念出发,我们求得了同质凝结的理论渐近解,与数值结果相比十分吻合。受激波管中凝结振荡的启发,我们进一步讨论了二维绕角膨胀流(Prandtl-Meyer流动)中的凝结。其中有界区域内的凝结,与喷管中凝结有类似之处。随着初始过饱和度的不断提高,可以依次发现亚临界、稳定凝结激波、振荡激波等现象。而在半平面区域绕角流动中,凝结在空间中呈现较强的二维效应,沿径向振荡型分布。可以发现,激波管中凝结随时间的振荡与半平面区域Prandtl-Meyer流动中凝结随径向的振荡有一定的时空相似性。
在完成关于基本流动中凝结的讨论之后,我们开展了风洞中凝结现象的研究。风洞中的凝结现象直接影响了流场的品质及模型实验的精准度。然而由于种种原因,国外的相关报道并不详实。本文开展的工作中,详细讨论了燃烧加热风洞中水蒸气的凝结,初步分析了常规高超风洞中氮气的凝结。燃烧加热风洞方面,着重讨论大扩张比喷管流动中所伴随的非平衡凝结过程对风洞试验段流场参数的影响,揭示凝结效应对燃烧加热风洞运行影响的物理机理。研究发现,由于燃烧加热风洞的试验气体中含有大量水蒸气,在喷管膨胀过程中,凝结是不可避免的。由于其高总温特点,在数值模拟时应考虑高温气体效应。水蒸气的凝结,改变了风洞试验段参数:静温提高,静压提高,马赫数降低。参数研究表明:较低的总温会导致更剧烈的凝结,结果会出现凝结后喷管出口静温反而更高的情况。此外,相同总温条件下,随着初始水蒸气含量的增大,凝结强度先增大后减小。氮气凝结方面,则尝试发展了可适用于模拟氮气同质凝结的数值程序。通过与以往文献对比,我们发现同质凝结理论并不能精确预测氮气的凝结。这是由于实验过程中试验气体里不可避免地混有极少量其它气体如水蒸气、二氧化碳等,在膨胀过程中沸点高的气体会首先形成极小的液滴,这些异质的凝结核提供了氮气异质成核的基础。因此,在评估实际风洞中氮气凝结时,兼顾考虑异质成核才能给出更准确的结论。
通过本文的讨论,我们得到了一系列有意义的结果:对凝结过程的气动分析提出了新的见解,强调了质量消耗的重要性,为深入分析气体凝结指出了新的方向;详细讨论了激波管中的凝结现象,给出的理论渐近解很好地解释了其流动特点;分析了有界区域内Prandtl-Meyer流动中凝结与喷管中凝结的类似性,提出了半平面区域Prandtl-Meyer流动中凝结的空间振荡分布;详细讨论了燃烧加热风洞中水蒸气的凝结,为燃烧加热风洞的发展提供了有价值的参考;初步探讨了氮气的凝结,为进一步研究氮气凝结奠定了基础,指明了方向。
|
|
|
|
1 |
李开泰,黄艾香;可压缩流动流函数微分方程的可解性和有限元解的误差估计[J];西安交通大学学报;1980年04期 |
2 |
张谋进,黄灿明;用于粘性可压缩流动数值计算的SIMPLE方法[J];航空动力学报;1997年01期 |
3 |
施卫平,胡守信,阎广武;计算可压缩流体流动的LB多速度模型[J];吉林大学自然科学学报;1996年02期 |
4 |
何冰;封卫兵;张武;武频;白文;李立;;基于非结构网格的Gas-Kinetic方法[J];计算机辅助工程;2009年01期 |
5 |
张池平,冯国泰,王民智;可压缩矢量形式湍动能与湍应力方程[J];哈尔滨工业大学学报;1992年06期 |
6 |
芮洪兴;双重介质油藏单相可压缩流动的特征差分格式及理论分析[J];应用数学学报;1994年01期 |
7 |
黄典贵;;统一的对流扩散型可压缩流体力学方程与解法[J];工程热物理学报;2006年03期 |
8 |
应展烽;陈志华;徐晓磊;董刚;范宝春;;可压缩混合层流动的三维大涡模拟[J];南京理工大学学报(自然科学版);2008年03期 |
9 |
高瑞泽;阎超;;超音速底部流动的大涡模拟研究[J];航空计算技术;2009年03期 |
10 |
景思睿;胡庆康;;叶栅可压缩流动反问题的一种数值解法[J];西安交通大学学报;1990年01期 |
11 |
张士杰,袁新,叶大均;全三维粘性流场计算的快速收敛方法[J];清华大学学报(自然科学版);2000年10期 |
12 |
赵兴艳,苏莫明,苗永淼;高分辨率激波捕捉格式及其CFD应用[J];计算物理;2001年02期 |
13 |
张锡文,李亨,姚朝晖;多孔介质/纯流体耦合区域内可压缩气体的流动[J];化工学报;2003年09期 |
14 |
陈榴;戴韧;杨爱玲;陈康民;;一种可压缩与不可压缩气体流动的统一算法及其应用[J];力学季刊;2009年02期 |
15 |
刘卫平;苏本庆;席德科;杨新铁;;可压缩流动声干涉现象也具有迈克尔逊-莫雷效应[J];机械科学与技术;2007年09期 |
16 |
钟诚文;李凯;孙建红;卓从山;解建飞;;基于格子Boltzmann方法的可压缩翼型绕流模拟(英文)[J];Transactions of Nanjing University of Aeronautics & Astronautics;2009年03期 |
17 |
邹滋祥;任意回转面上可压缩流动透平叶栅壁面最佳流速分布的计算方法[J];工程热物理学报;1980年04期 |
18 |
王平,朱自强,吕晓斌;迎风格式在二维非结构网格中的应用[J];北京航空航天大学学报;2000年05期 |
19 |
赵兴艳,苏莫明,苗永淼;非线性双曲型Euler方程的数值分析[J];甘肃工业大学学报;2000年03期 |
20 |
祁大同,王尚锦;奇点法在叶轮机械任意回转面亚音速可压缩流动计算中的应用[J];应用力学学报;1989年02期 |
|