变截面超音速汽液两相流升压过程的研究

本文针对超音速汽液两相流在变截面通道中的升压过程进行了实验研究,得到了变截面混合腔中超音速汽液两相流的压力变化规律,通过实验结果的分析得出了变截面混合腔中的渐缩部分的压力分布与出口压力无关、变截面通道的渐扩部分为一个单相流体扩压管、压力突变发生在变截面混合腔喉部和随着出口压力升高而激波强度增大、变截面超音速汽液两相流升压技术具有定流量特征等结论。     

环周进汽型超声速汽液两相流升压装置性能计算分析

对环周进汽型的变截面通道内超声速汽液两相流升压装置进行了实验及理论研究,实验中进汽压力为0.15～0.4MPa,进水压力为0.2～0.6 MPa。实验结果表明在不同的汽水参数条件下,混合腔内压力与温度分布呈现出相似的规律。在同一工况下,激波前混合腔内各点的压力基本保持不变,随着凝结激波的产生,压力突然增大。激波过程中蒸汽几乎全部凝结,激波过后温度分布趋于平缓。并在实验结果的基础上分别建立了水喷嘴、蒸汽喷嘴、混合腔内两相区和扩散段的数学模型,其预测的装置出口压与实验值之间的误差小于15%。     

单纤光纤探针测量空泡份额的实验研究

1前言气（汽）液两相流广泛应用于工业过程中。空泡份额是气（汽）液两相流的重要参数，它与流场、压力、热流密度和流型等密切相关。由于目前的理论计算模型还有较大的局限性，因此实验测量是研究气（汽）液两相流空泡份额的最主要的研究手段。经过近几十年的研究，各国学者开发出了许多有价值的空泡份额测量方法，其中包括平均空泡份额测量方法和局部空泡份额测量方法，这些测量方法本身各有各的局限性和一定的针对性。本文研究水平管束间汽液两相流流动和沸腾传热特性，流道中加热管柬的存在对流体流动产生强烈的影响，使管束间汽液两相…     

汽液两相混合物的加速与激波的热力学分析

本文分析了单相流和两相流通过喷管实现超音速流动的过程,建立数学模型并求解了两相流动在流道内加速至超音速的过程和超音速两相流动产生激波前后的热力学参数,得出了两相流动激波前后在不同升压比下参数的变化,两相流动激波前后的最高升压比与波前马赫数的关系式,分析了激波前后的熵产和 损。     

气固两相介质音速研究

本文分析计算了不同气固两相介质的音速,与实验结果比较得到了较好的吻合,验证了普朗特的理论分析.结合汽液两相音速分析表明,求解两相介质音速时可以不考虑相变的影响.以两相音速为基础的激波结果与理想气体激波结果的趋势吻合,较基于单相音速的激波更合理,说明了两相介质音速结果的正确性.     

毛细抽吸回路中反向式蒸发器内汽、液两相流动的理论分析

本文在建立物理模型的基础上，采用了一个二维计算模型对CPL的单个蒸发器的管内两相流动状态进行了分析和数值计算，分别得出汽、液相流体的速度、压力的分布．并完成系统的过程分析．     

纳米流体脉动热管的流动与传热性能研究

建立了纳米流体脉动热管的分析模型,将相变传热项合理引入了汽塞能量微分方程;液塞动量方程中考虑了剪切力的影响;纳米流体的物性采用了当量处理方法.采用数值迭代方法进行了求解,得到了汽塞压力、温度、质量的变化波形,分析了波形的频率,进而解释了初始条件、重力等对脉动热管的流动与传热影响的机理.     

自发凝结流动数值模拟方法及其在Laval喷管中的应用

本文对存在自发凝结的湿蒸汽两相流动建立了完全欧拉坐标系统下的数理模型。采用考虑了真实流体性质的 LU-SGS-GE隐式时间推进算法和改良型高精度、高分辨率MuSCL TVD差分格式求解存在自发凝结的汽液两相流动控制方程组。文中水及水蒸汽性质数据全部取自IAPWS-IF97国际标准公式。对某Laval缩放喷管内的湿蒸汽自发凝结流动的数值模拟结果表明,本文所采用的数理模型及计算方法是有效和可靠的。     

低压、低含汽量蒸汽-水两相临界流

一、引言 两相临界流研究高温、高压系统破口事故时,汽(气)-液两相流体的最高排放速度,是核电安全中的重要课题。由于两相流动的复杂性,两相临界流理论还没有到成熟阶段。不同试验者的数据,尤其是在低含汽量(x<0.02)时相差较大。     

ρ-VOF:一种可清晰模拟自由界面的两相流方法

文章通过对EFM(effective field modeling)模型进行简化, 消除了原模型的非守恒性项和非双曲性特性项, 发展了一种基于密度的气液两相流模拟方法: ρ-VOF方法.利用体积分数信息对控制单元内的自由界面进行重构, 得到了控制单元内流体的空间分布, 并采用AUSM+-up格式获得考虑气液流体接触间断信息的对流通量.新方法可统一处理激波间断和接触间断的相互作用, 保持自由界面的尖锐性, 并且其计算量与自由界面的空间复杂度无关.最后, 数值模拟了液体激波管气液激波管和气体激波跨二维液滴传播等问题, 并与文献结果进行对比, 验证了本方法在气液两相流模拟中的准确性.      

气液两相流压力波色散特性实验研究

设计了可调频式压力扰动源的气液两相流压力波实验装置,实验研究了垂直上升管内气液两相流泡状流、弹状流压力波的色散规律。实验结果表明,对泡状流,在实验范围内,压力波的传播速度及其衰减跟扰动频率有关,随着扰动频率的增加,波速及其衰减都增加;工质的流速对压力波的色散特性没有影响。结合数值模拟结果,验证了泡状流压力波色散特性的临界频率现象,即高于临界频率,压力波色散特性消失,本文分析了相应的物理机制。对弹状流,压力波同样具有典型的色散特性,已有研究结果还不能预测其色散规律。     

Calculation of shock wave propagation in water containing reactive gas bubbles

The entry of a shock wave from air into water containing reactive gas (stoichiometric acetylene–oxygen mixture) bubbles uniformly distributed over the volume of the liquid has been numerically investigated using equations describing two-phase compressible viscous reactive flow. It has been demonstrated that a steady-state supersonic self-sustaining reaction front with rapid and complete fuel burnout in the leading shock wave can propagate in this bubbly medium. This reaction front can be treated as a detonation-like front or “bubble detonation.” The calculated and measured velocities of the bubble detonation wave have been compared at initial gas volume fraction of 2 to 6%. The observed and calculated data are in satisfactory qualitative and quantitative agreement. The structure of the bubble detonation wave has been numerically studied. In this wave, the gas volume fraction behind the leading front is approximately 3–4 times higher than in the pressure wave that propagates in water with air bubbles when the other initial conditions are the same. The bubble detonation wave can form after the penetration of the shock wave to a small depth (~300 mm) into the column of the bubbly medium. The model suggested here can be used to find optimum conditions for maximizing the efficiency of momentum transfer from the pressure wave to the bubbly medium in promising hydrojet pulse detonation engines.     

悬浮RDX炸药粉尘爆轰的数值模拟

用两相流模型对悬浮RDX炸药粉尘爆轰波进行了数值模拟。RDX炸药颗粒在爆轰波阵面后的高温高速气流中加速并升温,颗粒表面发生熔化。参考液滴在高速气流作用下剥离的效应,假设炸药熔化部分在高速气流的作用下发生剥离,破碎成极小的颗粒,瞬时发生分解反应,释放出能量支持爆轰波传播。数值模拟了在不同粒径和浓度的悬浮RDX炸药粉尘中爆轰波的发展与传播过程,得到了爆轰波流场中气-固两相的物理量分布,并确定了爆轰波参数。在较低的RDX粉尘浓度条件下,爆轰波阵面压力的峰值曲线出现振荡。当RDX粉尘浓度在80~150 g/m3时,数值模拟得到的爆轰波阵面压力峰值曲线的振荡是规则的;当RDX粉尘浓度为70 g/m3时,爆轰波阵面压力峰值曲线出现不规则振荡。     

超高压下冲击波速度直接测量技术

针对超高压下透明材料的高压离化机理,分析了透明材料中冲击波直接诊断技术的基本方法. 利用Drude-自由电子气模型,分析了不同冲击压力下冲击波阵面反射率的变化. 从理论上比较了不同探针光波长反射率的区别,发现探针光波长为660 nm时比探针光波长为532 nm时获得的冲击波阵面反射率要高. 对探测器"致盲"问题也进行了研究. 通过分析反射信号的时间顺序和强度大小,发现"致盲"效应是由X光对透明窗口离化引起的. 同时,发现方波驱动脉冲平台的前沿到达时刻和X光离化效应出现的时刻相同,冲击波信号到达时刻晚于X光离化时刻. 通过实验结果,得到蓝宝石中冲击波速度为35 km/s时,其波阵面的反射率约为40%. 通过理论分析和实验数据比对的方法,验证了蓝宝石中的减速曲线. 给出了加蓝宝石窗口后的测速公式. 经过和实验对比,确认了测速公式的正确性. 关键词： 冲击波 光学诊断 成像 干涉仪     

Peculiarities of evolution of shock waves generated by boiling coolant

Simulation of compression wave generation and evolution at the disk target was performed for the case of explosive-type boiling of coolant; the boiling is initiated by endwall rupture of a high-pressure pipeline. The calculations were performed for shock wave amplitude at different times and modes of pipe rupture. The simulated pressure of a target-reflected shock wave is different from the theoretical value for ideal gas; this discrepancy between simulation and theory becomes lower at higher distances of flow from the nozzle exit. Comparative simulation study was performed for flow of two-phase coolant with account for slip flow effect and for different sizes of droplets. Simulation gave the limiting droplet size when the single-velocity homogeneous flow model is valid, i.e., the slip flow effect is insignificant.     

极端条件两相界面流与反应流机理、模型与算法研究进展

清华大学喷雾燃烧与推进实验室长期专注于极高速、强可压和高瞬变等极端条件下的两相流和反应流前沿科学问题,并致力于应用基础研究成果解决航空航天动力与推进系统的关键技术难题。综述了实验室近些年在极端条件下两相流动和含化学反应流动物理机理、数理模型与数值算法等方面的研究进展。首先,介绍了实验室发展的耦合高瞬变相变过程的强可压缩气液两相界面流的数理模型和高精度数值方法,以及针对激波受气液（曲）界面约束情况下,描述非定常激波透射/反射（如波角、波强等物理量关系）的激波动力学分析方法。其次,基于上述模型、算法与分析方法,实验室研究了激波液滴相互作用、高速液滴撞击壁面等一系列问题,解析了上述高瞬变过程中复杂波系与界面的时空演化过程。以被激波或壁面冲击的液滴内流体空化初生为例,揭示了曲界面汇聚膨胀波诱导流体空化的机理,推导了预测空化初生位置的理论公式。最后,介绍了面向发动机燃烧室内的强可压缩两相喷雾反应流动,实验室开发了基于Euler-Lagrange框架的高性能数值仿真软件TURFsim,并成功用于真实复杂几何结构的航空发动机燃烧室和超声速燃烧室的数值模拟。以典型的超声速混合层流动数值模拟为例,总结了斜激波增强混合层混合特性的规律及其物理机理,获得了极限条件火核生成及火焰的传播模式与机理,详细分析了液雾弥散与蒸发、小激波和局部爆震波的时空演化特性,提出使用"第三Damk?hler数Da_Ⅲ"定量表征燃烧模式,应用该无量纲参数成功进行了局部准等容燃烧过程的辨识与演化分析。上述研究结果对于航空航天发动机燃烧室复杂物理过程的理解及工程设计具有重要价值。      

Nonequilibrium Zeldovich-von Neumann-Doring theory and reactive flow modeling of detonation

This paper discusses the Nonequilibrium Zeldovich-von Neumann-Doring (NEZND) theory of self-sustaining detonation waves and the Ignition and Growth reactive flow model of shock initiation and detonation wave propagation in solid explosives. The NEZND theory identified the nonequilibrium excitation processes that precede and follow the exothermic decomposition of a large high explosive molecule into several small reaction product molecules. The thermal energy deposited by the leading shock wave must be distributed to the vibrational modes of the explosive molecule before chemical reactions can occur. The induction time for the onset of the initial endothermic reactions can be calculated using high pressure-high temperature transition state theory. Since the chemical energy is released well behind the leading shock front of a detonation wave, a physical mechanism is required for this chemical energy to reinforce the leading shock front and maintain its overall constant velocity. This mechanism is the amplification of pressure wavelets in the reaction zone by the process of de-excitation of the initially highly vibrationally excited reaction product molecules. This process leads to the development of the three-dimensional structure of detonation waves observed for all explosives. For practical predictions of shock initiation and detonation in hydrodynamic codes, phenomenological reactive flow models have been developed. The Ignition and Growth reactive flow model of shock initiation and detonation in solid explosives has been very successful in describing the overall flow measured by embedded gauges and laser interferometry. This reactive flow model uses pressure and compression dependent reaction rates, because time-resolved experimental temperature data is not yet available. Since all chemical reaction rates are ultimately controlled by temperature, the next generation of reactive flow models will use temperature dependent reaction rates. Progress on a statistical hot spot ignition and growth reactive flow model with multistep Arrhenius chemical reaction pathways is discussed. The text was submitted by the authors in English.     

Influence of air pressure on the performance of plasma synthetic jet actuator

Plasma synthetic jet actuator(PSJA) has a wide application prospect in the high-speed flow control field for its high jet velocity.In this paper,the influence of the air pressure on the performance of a two-electrode PSJA is investigated by the schlieren method in a large range from 7 k Pa to 100 k Pa.The energy consumed by the PSJA is roughly the same for all the pressure levels.Traces of the precursor shock wave velocity and the jet front velocity vary a lot for different pressures.The precursor shock wave velocity first decreases gradually and then remains at 345 m/s as the air pressure increases.The peak jet front velocity always appears at the first appearance of a jet,and it decreases gradually with the increase of the air pressure.A maximum precursor shock wave velocity of 520 m/s and a maximum jet front velocity of 440 m/s are observed at the pressure of 7 k Pa.The averaged jet velocity in one period ranges from 44 m/s to 54 m/s for all air pressures,and it drops with the rising of the air pressure.High velocities of the precursor shock wave and the jet front indicate that this type of PSJA can still be used to influence the high-speed flow field at 7 k Pa.     

Influence of ambient pressure on the performance of an arc discharge plasma actuator

The arc discharge plasma actuator (ADPA) has wide application prospects in high‐speed flow control because of its local heating effect and strong disturbance. In this paper, the influence of ambient pressure, which ranges from 3 to 20 kPa, on the performance of a two‐electrode ADPA is investigated by a schlieren system. The duration of the arc heated region, as well as its area, is extracted by image processing. As the ambient pressure increases, different flow field evolutions occur. The duration of the ADPA heated region increases with the ambient pressure. The maximum duration reaches 1185.3 µs at 20 kPa. The velocity of the discharge‐induced blast shock wave first decreases gradually and then remains at 345 m/s for all air pressures. The blast shock wave has a higher velocity at lower pressures when it is freshly produced. A maximum blast shock wave velocity of 582 m/s is observed at the pressure of 7 kPa. The arc heated region is not sensitive to ambient pressure, but the deposited energy from the arc increases when the pressure increases.