燃气非平衡流再生冷却流动传热数值模拟

为准确预测液体火箭发动机推力室身部再生冷却换热状况,采用数值模拟方法,对燃气、推力室壁和超临界气氢进行三维流动和换热耦合计算。采用6组分9步反应动力学模型计算燃气的非平衡化学反应,采用DO模型计算燃气辐射换热,考虑超临界气氢物性随温度和压力的变化。获得了室壁温度场、燃气及冷却剂流场。结果表明,Redlich-Kwong方程、Peng-Robinson方程、Lucas法、TRAPP法能分别准确计算超临界氢的密度、定压比热容、粘度、导热系数,采用燃气非平衡流计算所得流场值更符合实际情况。     

转子叶片径向受限的“冲击-气膜出流”冷却结构流动与换热

为揭示转子叶片径向受限的"冲击-气膜出流"冷却结构流动换热规律,以某型双层壁叶片肋化分割形成的冷却单元为研究对象,通过数值模拟的方式,对冲击雷诺数Rej,旋转数Ro,无因次温比(Tw-Tf)/Tw等参数变化下流场和换热特性变化规律展开研究。结果表明:在哥氏力和离心力作用下,受限空间内存在射流偏转、径向二次流动以及二次冲击等现象;流动的径向受限可抑制射流偏转,强化冲击换热;相同的旋转数Ro下,逆转向冲击(叶背区)换热努赛尔数Nu比顺转向冲击(叶盆区)高8%。在研究的参数范围内,数值模拟和试验结果说明径向受限周向出流结构能有效的抑制旋转对换热的削弱。     

冷却系统阻力对汽车总阻力影响的数值预测模型研究

采用CFD方法数值模拟了小型汽车发动机冷却系统对全车阻力的影响特性.汽车的内部及外部(发动机舱及舱内)均设置了数值模拟网格.分别计算了散热器及风扇对阻力的影响,采用本校实验结果(散热器压力损失)及制造商提供的数据(散热器及风扇的质量流率)确定初始及边界条件,为基于RNS的CFD提供输入.对汽车轮胎运动模式及行驶速度对散热系统阻力及汽车总阻的影响特性也进行了计算分析.预测结果表明,冷却系统阻力占汽车总阻的6.9％,与实测值吻合.     

Effect of flexibility on unsteady aerodynamics forces of a purely plunging airfoil

Introducing flexibility into the design of a vertically flapping wing is an effective way to enhance its aerodynamic performance. As less previous studies on the aerodynamics of vertically flapping flexible wings focused on the lift generated in a wide range of angle of attack·a 2D numerical simulation of a purely plunging flexible airfoil is employed using a loose fluid–structure interaction method. The aerodynamics of a fully flexible airfoil are firstly studied with the flexibility and angle of attack. To verify whether an airfoil could get aerodynamic benefit from the change in structure, partially flexible airfoil with rigid leading edge and flexible trailing edge were further considered. Results show that flexibility could always reduce airfoil drag while lift and lift efficiency both peak at moderate flexibility. When freestream velocity is constant, lift is maximized at a high angle of attack about 40° while this optimal angle of attack reduces to 15° in drag-balanced status. The airfoil drag reduction, lift augmentation as well as efficiency enhancement mainly attribute to the passive pitching other than the camber deformation. Partially deformed airfoil with the longest length of moderate flexible trailing edge can achieve the highest lift. This study may provide some guidance in the wing design of Micro Air Vehicle (MAV).     

Oxygen concentration variation in ullage influenced by dissolved oxygen evolution

To determine the oxygen concentration variation in ullage that results from dissolved oxygen evolution in an inert aircraft fuel tank, the CFD method with a mass transfer source is applied in the present study. An experimental system is also designed to evaluate the accuracy of the CFD simulations. The dissolved oxygen evolution is simulated under different conditions of fuel load and initial oxygen concentration in ullage of an inert fuel tank with stimulations of heating and pressure decrease. The increase in the oxygen concentration in ullage ranges from 0.82% to 5.92% upon stimulation of heating and from 0.735% to 12.36% upon stimulation of a pressure decrease for an inert ullage in the simulations. The heating accelerates the release of the dissolved oxygen from the fuel by increasing the mass transfer rate in the mass transfer source and decreasing the pressure, thereby accelerating the dissolved oxygen evolution by increasing the concentration difference between the gas and the fuel. The time constant that represents the oxygen evolution rate is independent of the initial oxygen concentration in ullage of an inert tank but depends closely on the fuel load, temperature and pressure. The time constant can be fitted using a polynomial equation relating the fuel load to temperature in the heating stimulation with an accuracy of 4.77%. Upon stimulation of a pressure decrease, the time constant can be expressed in terms of the fuel load and the pressure, with an accuracy of 5.02%.     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

层板冷却结构流阻特性的实验与数值模拟

在吸气式风洞上对两种不同内部流动结构的单层多孔层板进行了流阻实验，采用沿程阻力关系式得出了其流阻特性。应用软件FLUENT，采用隐式有限体积法求解Navier—Stokes雷诺时均湍流方程，对层板冷却结构的内部流动进行了数值模拟，湍流模型采用Realizable κ-ε双方程模型，近壁处湍流利用壁面函数法处理，采用SIM—PLE算法求解速度与压力的耦合。计算表明，层板内部流场结构十分复杂，射流冲击后在扰流柱前反卷形成驻涡。层板结构内部存在流速很低的区域。计算得出的流阻特性与实验结果进行了比较，两者符合较好。     

叶排间轴向间距对时序效应影响的数值研究

时序效应具有改进多级叶轮机效率的潜力。文中在已有认识基础上，对叶排间轴向间距对时序效应影响进行深入的二维非定常数值研究。算例采用1．5级静／动／静布局的涡轮叶栅。结果发现，叶排间轴向间距对时序效应自身影响很大，会使涡轮效率的变化幅度以及效率极值出现的位置发生明显变化。另外，处于两静叶排正中的动叶位置的前后移动或缩短静叶排间的距离均能带来时序效应影响下效率的整体提升。     

Theoretical properties of offset bipolar electric field solitary structures in space plasmas

The bipolar electric field solitary (EFS) structures have been frequently observed in the near Earth plasma regions, such as auroral zone, magnetopause, cusp regions, and magneto-tail. Sometimes these structures are observed as offset bipolar structures. In this paper, the properties of the offset bipolar EFS structures parallel to the magnetic field are studied with an ion fluid model in a cylindrical symmetry by considering electrostatic condition. The model results show that the offset bipolar EFS structures can develop from both ion-acoustic waves and ion cyclotron waves, and propagate along the magnetic field line in the space plasmas if plasma satisfies some conditions. The offset bipolar EFS structures can have both polarities. It will be first negative pulse and then positive pulse if the initial electric field E₀ < 0 or reverse in order if E₀ > 0. The amplitude of the offset bipolar EFS structures first decreases and then increases with the wave propagation velocity. Some results from our model are consistent with the observations.