汽液相变对流场结构影响的数值分析论文.pdfVIP

汽液相变对流场结构影响的数值分析 王美利1， 罗喜胜2， 杨基明1 (1、中国科学技术大学，力学和机械工程系，合肥230027) (2、埃茵霍温理工大学，应用物理系，荷兰) 摘要：本文主要以数值方法研究非定常流场中的汽液相变对流场结构的影响。在含有可凝结 气体的跨音速或超声速流场中，由于可凝结气体发生凝结并向周围环境释放热量，从而对流 场局部进行加热而导致流场结构发生变化。在数值计算中，本文采用基于有限体积方法的2一D 非定常欧拉方程，并结合水蒸气结核与凝结生长模型对流场进行模拟。结果表明由于水蒸气 的凝结释放热量对流场局部进行不均匀加热，使得流场参数、结构等都发生明显的变化。 1．背景 相变是一个非常普遍的、基本的物理过程，其中包括凝结、蒸发等现象。在我们的日常 生活中凝结和蒸发是一种很普遍的现象，如地球上的水蒸发到大气中，当水蒸汽的相对饱和 度达到过饱和或大气中有足够多的灰尘粒子时，就会发生凝结形成云等现象。但在没有灰尘 粒子存在的情况下，凝结核的形成是通过水蒸汽本身自发形成的，我们通常称这种过程为同 质结核过程。Volmer和Weber首先给出了一个经典的凝结核的形成公式[1】。在这种情况下， 相对饱和度就成为发生凝结的重要参数。一般情况下，相对饱和度S是压力和温度的函数， 且定义为S=P，／p,仞，Pv为水蒸汽的分压，风为给定温度的水蒸汽的饱和压力。当环境温 度下降到某一温度时，饱和度S增大并达到过饱和即SI，这时水蒸汽会发生凝结现象，并 向周围环境释放热量，使得环境局部温度升高。所以把这种物理过程与空气动力学过程相结 合，在很多工程与技术上的应用是非常重要的，如喷管流、涡轮机、飞机在大气层中飞行和 实验流体力学等。 超音速喷管流中凝结现象的研究开始于30年代。当时学者们在用湿空气来启动超音速喷 管时，发现喷管中存在不可解释的扰动。后来由Prandtl[2]首先在圆弧形喷管实验中显示了 一种称“X．shock”的激波并正确的解释这种激波是水蒸汽凝结的结果。之后各种更深入的喷 管中的水蒸气凝结现象的研究得到了发展【3．6】。在国内，对流场中的凝结现象研究得比较少， 杨顺华等在超燃研究中研究了水蒸汽凝结在燃烧加热风洞中的应用【7】。 在大气层中往往含有很大比例的水蒸气，当飞机在大气层中做跨音速或超音速飞行时由 于水蒸汽发生凝结并向周围气体释放潜热，使机翼周围的温度、压力发生扰动，因而影响飞 机的升力、阻力等空气动力学性能。在这方面工作上，Schnerr用无粘的2．D欧拉方程并结合 数值求解大气层中湿空气绕机翼的跨音速问题，并给出升力和阻力 Hill的经典结核理论[8]， 的变化，结果表明，在与干空气情况下相比，由于潜热的影响，机翼的升力下降60％，压力 阻力提高200％[9]。 在含有凝结现象的可压缩流问题上虽然已经研究了很久，但直到现在模型理论仍然不够 成熟，现有的结核模型理论在某些条件下与实验数据还有不相符合的现象，所以结核和凝结 过程与流场的相互作用还没有得到很好的解决。在本文中，我们主要用数值方法研究湿氮气 在变截面管中的凝结现象以及凝结释放的潜热对流场的影响。为此，我们做如下假设：1、两 种气体的混合在流场中各自遵守理想气体状态方程；2、在凝结生长过程中，凝结的水滴直径 很小，它所占的体积可以忽略不计；3、小水滴随气体一起流动，并且不发生速度滑移。在计 算过程中，我们采用2．D无粘的非定常欧拉方程，并结合经典的结合理论[8】和水滴增长模型 [10】实现流场的数值模拟【1l】。 330 汽液相变对流场结构影响的数值分析 2．数值方法 在数值计算中，我们采用2-D非定常欧拉方程和理想气体的状态方程，并结合结核理论 和水滴生长模型，所以整个控制方程可表示为： aU押aG． 百+面+面2s ‘1’ 其中： 0 p pu pv 0