为了提高内燃机的燃烧效率并降低污染物的排放，需要对燃烧过程所涉及的主要流动过程如燃油雾化、液滴蒸发、油气混合以及化学反应等进行优化设计。闪急沸腾喷雾作为一种能够有效改善燃油雾化质量的手段，在内燃机中得到了越来越多的应用。在该喷雾过程中， 过热液滴会因为内部气泡的不断生长而完全破碎并产生大量微小液滴，这种效应不仅能够有效改善喷雾流场的宏观特性：缩短液核区长度、增大喷雾扩张角度、缩短液滴贯穿距离，还能够增强液滴的微观特性，包括减少液滴粒径、提高液滴蒸发速率、促进燃油与空气的快速混合等。但目前对于闪急沸腾喷雾的应用还存在一些实际问题，比如在接近饱和状态下，过热管道内部流动将出现不稳定状况；在高过热状态下，多喷嘴喷雾结构将出现严重的坍缩现象，喷雾结构的穿透深度将大大增加甚至出现撞壁现象，因此有必要对闪急沸腾喷雾过程以及过热流动现象进行深入的研究。
       过热燃料的闪急沸腾喷雾是一个复杂的物理过程。 目前对过热喷雾的研究主要以实验方法为主，研究手段包括高速摄像机、纹影仪、米氏散射仪、 激光诱导荧光仪、激光多普勒测速、相位多普勒测速仪等。利用上述实验装置，人们对喷雾的穿透深度与宽度、液滴的粒径与速度分布、蒸气的浓度分布等喷雾流场特征量进行了观测分析，获得了大量实验数据，但是诸多现象背后的物理机制仍未清楚。且工程上运用过热喷雾现象进行产品设计存在着巨大的需求，而充分理解过热喷雾现象的物理机制对产品的优化设计十分关键， 因此有必要开展这方面的数值模拟研究工作。 然而现有的过热流动和闪蒸雾化模型充满唯象论色彩，缺乏真实物理背景的支持，有必要对其开展针对性的改进工作，以满足工程实际和研究需求。

       根据过热流动和闪急沸腾喷雾过程的主要流动特点，本文分别利用一维气泡生长模型、拉格朗日粒子方法和连续欧拉方法对过热流体中的气泡生长过程、外部闪蒸喷雾过程和过热管道内部流动过程进行了一系列的数值模拟研究工作。我们首先对过热流体中的气泡生长过程进行了详尽研究，运用直接数值离散方法和基于动态热边界层的数值方法对处于常压低过热度、常压中高过热度和低压中等过热度条件下的过热液体中的气泡生长状况进行了数值求解和模型验证工作。通过数值模拟计算获得了在相同过热条件下气泡在过热液体和过热液滴两种情形下的生长规律，并进行了对比研究。通过对过热液滴中的气泡生长动力学方程的小扰动分析，得到了过热液滴-气泡系统中的扰动发展方程，并对水、正戊烷、正己烷等三种液体在不同过热状态下的破碎情况进行了数值预测分析。对过热液滴的蒸发机制进行了研究，并初步探讨了液滴破碎与蒸发过程之间的关系。针对现有的破碎模型对过热液滴破碎过程描述不够清晰、 并且缺少足够的理论基础的现状，在前期研究的基础上发展了两种破碎模型来反映不同过热度下的过热液滴破碎特征：在低过热度下，基于现有的气动力破碎模型，将气泡的存在对液滴受力的影响考虑进去，发展出了过热条件下的液滴气动力破碎模型；在高过热度下，利用过热液滴-气泡系统中的扰动增长率发展模型来预测过热液滴的临界破碎状态，并利用能量守恒原理和能量转换系数来对子液滴特性参数进行计算，建立了一个考虑了液滴“微爆”效应的热力学破碎模型。针对复杂的过热流动现象，在开源平台上搭建了一个基于 VOF 方法的求解器，在其构造过程中考虑了可压缩效应、不可冷凝气体输运以及气液相变过程，并通过推导将多组分能量输运方程改造成了以温度为变量的输运方程形式。针对当前的过热相变模型大多基于均匀气泡数密度假设的现状，在 VOF 求解器中发展了基于气泡成核与输运过程的过热相变模型，在其构造过程中分别采用了异质成核模型、气泡数密度输运模型、过热气泡生长速率模型以及气液交界面密度估算公式来对应相关物理过程，同时研究了过热相变模型中的参数取值对数值模拟结果的影响规律， 并根据过热管道内部和外部的相变特征改进了原始的ELSA 模型，从而实现对过热管道内外流场的统一预测。 度下的过热液滴破碎特征：在低过热度下，基于现有的气动力破碎模型，将气泡的存在对液滴受力的影响考虑进去，发展出了过热条件下的液滴气动力破碎模型；在高过热度下，利用过热液滴-气泡系统中的扰动增长率发展模型来预测过热液滴的临界破碎状态，并利用能量守恒原理和能量转换系数来对子液滴特性参数进行计算，建立了一个考虑了液滴“微爆”效应的热力学破碎模型。针对复杂的过热流动现象，在开源平台上搭建了一个基于 VOF 方法的求解器，在其构造过程中考虑了可压缩效应、不可冷凝气体输运以及气液相变过程，并通过推导将多组分能量输运方程改造成了以温度为变量的输运方程形式。针对当前的过热相变模型大多基于均匀气泡数密度假设的现状，在 VOF 求解器中发展了基于气泡成核与输运过程的过热相变模型，在其构造过程中分别采用了异质成核模型、气泡数密度输运模型、过热气泡生长速率模型以及气液交界面密度估算公式来对应相关物理过程，同时研究了过热相变模型中的参数取值对数值模拟结果的影响规律， 并根据过热管道内部和外部的相变特征改进了原始的ELSA 模型，从而实现对过热管道内外流场的统一预测。

       通过上述研究工作， 我们获得了关于过热流动和闪急沸腾喷雾现象的丰富研究成果。首先通过数值计算发现了气泡在过热液体和过热液滴两种物理模型下的不同生长规律：在低过热度和较大的液滴直径情形下，气泡在过热液体和过热液滴中的生长规律基本相同；但在高过热度和较小的液滴直径的情形下，过热液滴中的气泡会出现二次加速生长现象，出现这种现象的主要原因是处于惯性生长阶段的气泡在生长到与液滴初始尺寸相当时，其所受惯性力会有明显的放大作用。过热液滴-气泡系统中的扰动发展方程为一元三次方程，方程的实数解表征了在气泡不同生长时刻系统内的扰动增长情况。通过对扰动率的时间积分，可以获得系统扰动量的变化情况，当其满足临界破碎准则时，过热液滴-气泡系统发生破碎，这对于建立合适的过热液滴破碎模型具有重要意义。 过热液滴的破碎模型和蒸发模型对于可靠的闪急沸腾喷雾模拟来说缺一不可，能量转换系数的取值大小会影响到喷雾流场的计算结果，对于本文算例来说 0.05-0.1之间的取值较为合适。 通过与实验结果的对比研究，验证了本文发展的过热液滴破碎模型能够较好的捕捉不同喷管构型和不同过热度下地闪急沸腾喷雾现象，并可以用于高过热度下的喷雾坍缩现象研究。 本文研究认为喷雾结构附近的卷吸结构随着液滴过热程度的提高而迅速增强， 并造成邻近油束乃至整个油束结构之间的相互干扰过程，是导致多喷嘴喷雾结构在高过热度下出现坍缩现象的主要原因。通过大量的数值计算，验证了本文发展的过热相变模型能够对不同过热状态下的管道出口流量和管道内部的压力分布、气泡含量分布做出很好的预测，并且能够较为合理地预测过热管道内部的成核、蒸发等不同的过热流动现象，从而为流场分析工作提供很好的支撑。在过热相变模型中， 异质成核因子的大小对不同过热状态下的管道内部流场计算有重要影响，其大小能在很大程度上决定过热流场的成核水平，并决定成核过程发生的临界过热度，本文建议的合理取值为 1e-6。通过数值计算发现在不同过热度下小长径比的过热管道分别表现出两种不同的内部流场特征： 在低过热度下，管道内的过热液体整体处于亚稳定状态，表现出典型的过冷空化流动特征；在流场过热度超过某一临界值后，过热管道内部会有成核过程的发生，流场表现出典型的外部沸腾特征。对于大长径比的管道，当其来流温度接近来流压力所对应的饱和温度时，过热管道内部将完全气化，进入到所谓的内部闪急沸腾状态。