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

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

20世纪后30年中国西北西部降水年代际变化机理分析

在全球气候变暖背景下1980年代中期中国西北西部气候发生了年代际突变,降水明显增加且主要是极端降水的贡献. 本文用欧洲中期预报中心再分析资料ERA-40计算了瞬变涡动、环流、水汽输送及其散度等的年代际变化并进行了分解分析. 结果表明, 1987年后大气环流年代际变化形成"东高西低"环流背景场,西风带向该区的水汽输送增加, 水汽辐合增强,其中主要是定常波水汽输送的贡献, 其结果是造成该区及周边空气水汽含量上升, 有利于降水或极端降水的增加. 夏季该区瞬变涡动增强且造成水汽辐合, 增加了降水天气出现的频率; 冬季瞬变涡动减弱, 天气过程减少, 说明降雪增加的原因应该是空气水汽含量增加, 导致降雪强度增加. 因此该区降水的增加即有中高纬度大气变暖变湿的贡献,也有瞬变涡动的动力作用. 未来西北西部降水的趋势主要取决于瞬变涡动与大气水汽含量变化的竞争．     

基于瞬变流和遗传算法的管道泄漏辨识

输水系统管道的泄漏影响瞬变水击波波形的畸变和衰减特性,对泄漏的定位和定量是泄漏辨识研究的重要环节.考虑水力瞬变中非恒定摩阻的影响,将泄漏参数反映在建立的管道瞬变流时域数学模型中,给出离散网格和求解方法,分析不同泄漏参数对瞬变水击波的影响,用遗传算法结合瞬变流时域模型的特征线法最小化管道测点压力实测和计算响应值的均方误差来辨识泄漏参数,其中实测信号首先经过预滤波和小波去噪处理.通过模型试验对瞬变流数学模型和泄漏反问题分析模型作了率定和验证,结果表明此法能有效的辨识出小泄漏参数,并且检测信号主动,阀门扰动可操作性强.     

现行电力系统热经济性状态方程

本文用系统工程观点,定义了热力发电厂热经济性系统状态方程。指出,状态方程的结构与系统的结构一一对应,系统的热经济性取决于系统的结构及系统所处热力学状态。该方程的出现为电厂热经济性的计算、分析与综合提供了新的基础和研究方法。     

Phonon spectrum and related thermodynamic properties of microcrack in bcc-Fe

The phonon spectrum and the related thermodynamic properties of microcracks in bcc-Fe are studied with the recursion method by using the Finnis--Sinclair (F--S) N-body potential. The initial configuration of the microcracks is established from an anisotropic linear elastic solution and relaxed to an equilibrium by molecular dynamics method. It is shown that the local vibrational density of states of the atoms near a crack tip is considerably different from the bulk phonon spectrum, which is closely associated with the local stress field around the crack tip; meanwhile, the local vibrational energies of atoms near the crack tip are higher than those of atoms in a perfect crystal. These results imply that the crack tip zone is in a complex stress state and closely related to the structure evolution of cracks. It is also found that the phonon excitation is a kind of local effect induced by microcracks. In addition, the microcrack system has a higher vibrational entropy, which reflects the character of phonon spectrum related to the stress field induced by cracks.     

Thermal post-bunkling analyses of functionally graded material rod

The non-linear governing differential equations of immovably simply supported functionally graded material (FGM) rod subjected to thermal loads were derived. The thermal post-buckling behaviors of FGM rod made of ZrO2 and Ti-6Al-4Vwere analyzed by shooting method. Firstly, the thermal post-buckling equilibrium paths of the FGM rod with different gradient index in the uniform temperature field were plotted, and compared with the behaviors of the homogeneous rods made of ZrO2 and Ti-6A1-4V materials, respectively. For given value of end rotation angles, the influence of gradient index on the thermal post-buckling behaviors of FGM rod was discussed. Secondly, the thermal post-buckling characteristics of the FGM rod were analyzed when the temperature difference parameter is changed while the bottom temperature parameter remains constant, and when the bottom temperature parameter is changed while the temperature difference parameter remains constant, and compared with the characteristics of the two homogeneous material rods.     

瞬变平面热源法测定常压下四氢呋喃水合物的导热系数

采用瞬变平面热源法,测定了常压、温度233.15～273.15K下四氢呋喃水合物的导热系数,其范围在0.45～0.54W·m-1·K-1之间。实验结果表明,常压下四氢呋喃水合物的导热系数呈现随温度升高而增大的趋势,与玻璃体的变化特性相同。常压下四氢呋喃水合物的导热系数与高压下的数值相比差别很小,但常压下测得的四氢呋喃水合物导热系数的玻璃体变化特性比高压下明显。     

The dynamical feature of transition of a Hamiltonian system to a dissipative system

The mechanism of generation and annihilation of attractors during transition from a Hamiltonian system to a dissipative system is studied numerically using the dissipative standard map. The transient process related to the formation of attracting basins of periodic attractors is studied by discussing the evolution of the KAM tori of the standard map. The result shows that as damping increases, attractors are mainly generated from elliptic orbits of the Hamiltonian system and annihilated by colliding with unstable periodic orbits originating from the corresponding hyperbolic orbits of the Hamiltonian system. The transient process also exhibits the general feature of bifurcation.     

新世纪的绿色环保能源

 我们伴随着日益严重的世界性能源短缺和环境污染走进了２１世纪。原油价格不断上涨的浪潮还在升温,而地球上的化石燃料如石油、煤炭等这些不可再生能源的储量在快速减少。那么,寻找无污染且可持续发展的绿色环保能源来取代化石燃料已是人类的共同课题之一。按照我国的２１世纪发展计划,即到２０５０年进入中等发达国家的行列,那么,我国那时的年人均ＧＤＰ要达到５０００至６０００美元（美国、日本２０世纪８０年代的水平）。设想那时的单位产值能耗可降至０．４１吨标准煤／千美元ＧＤＰ,可以推算出２０５０年一次能源总需求量约为３４．４亿吨标准煤,按１６亿人口计算,人均约为２．１５吨标准煤。     

考虑若干影响的一元蒸汽凝结流数值解

本文在一元蒸汽凝结流动量方程中引入了滑移系数Ｂ,考虑两相滑移的影响。并对Ｋｅｌｖｉｎ方程进行了实际气体修正。发展了表格插值的方法以高精度和快速计算水蒸汽热力性质。