某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

A transient bulk flow model with circular whirl motion for rotordynamic coefficients of annular seals

Bulk flow model with perturbation simplification has been used to calculate rotordynamic coefficients in annular seals which have significant influences on the dynamic behavior of rotors in turbomachinery. In this work, a transient bulk flow model with arbitrary rotor motion is developed, and the boundary conditions and friction factor in the model are calibrated with steady Computational Fluid Dynamics (CFD) analysis. The numerical solution scheme is developed based on the finite element method to obtain the transient reaction force in the seal clearance. With a periodic circular rotor orbit, the transient forces at multiple whirling frequencies are used to evaluate the rotordynamic coefficients. The leakage flowrate of CFD analysis has good agreement with experimental results and the calibrated parameters in bulk flow model are dependent on operating conditions. Although CFD calibration improves the accuracy of the perturbed bulk flow model, the direct damping is overestimated and the cross-coupled damping is underestimated. Compared with the perturbed model, the predictions of the transient bulk flow model are more agreeable with the experiment.     

Interface flux reconstruction method based on optimized weight essentially non-oscillatory scheme

Aimed at the computational aeroacoustics multi-scale problem of complex configurations discretized with multi-size mesh, the flux reconstruction method based on modified Weight Essentially Non-Oscillatory (WENO) scheme is proposed at the interfaces of multi-block grids. With the idea of Dispersion-Relation-Preserving (DRP) scheme, different weight coefficients are obtained by optimization, so that it is in WENO schemes with various characteristics of dispersion and dissipation. On the basis, hybrid flux vector splitting method is utilized to intelligently judge the amplitude of the gap between grid interfaces. After the simulation and analysis of 1D convection equation with different initial conditions, modified WENO scheme is proved to be able to independently distinguish the gap amplitude and generate corresponding dissipation according to the grid resolution. Using the idea of flux reconstruction at grid interfaces, modified WENO scheme with increasing dissipation is applied at grid points, while DRP scheme with low dispersion and dissipation is applied at the inner part of grids. Moreover, Gauss impulse spread and periodic point sound source flow among three cylinders with multi-scale grids are carried out. The results show that the flux reconstruction method at grid interfaces is capable of dealing with Computational AeroAcoustics (CAA) multi-scale problems.     

Numerical simulation of broken pin effects on the flow field and cooling performance of a double-wall cooling configuration

In this paper, the flow characteristics of the double wall structure are presented and the effect of the broken pin size on the cooling performance and flow field of the double wall configuration is investigated. A periodic plate model with seven units is adopted, and there are an impingement hole and a film hole in each unit. Under five blowing ratios, six different sizes of the broken pin are compared, and the double wall configuration without broken pins is taken as the baseline. The results show that if the broken pins height is too small, the cooling effectiveness usually cannot be improved. With the presence of broken pins with a height of more than 0.4, the effectiveness is improved due to the enhancement of reattachment and recirculation of coolant. With the increase of the broken pin height, the cooling effectiveness increases. However, the increase of the diameter does not always improve the cooling performance, since the limiting effect of the wall jet. In this study, Case 6 with the largest broken pin always has the best cooling performance, but also the largest flow resistance. In Case 6 temperature is reduced by almost 15 K compared to the baseline, and more areas have relatively higher cooling effectiveness.     

A hybrid multidimensional Riemann solver to couple self-similar method with MULTV method for complex flows

Since proposed, the self-similarity variables based genuinely multidimensional Riemann solver is attracting more attentions due to its high resolution in multidimensional complex flows. However, it needs numerous logical operations in supersonic cases, which limit the method’s applicability in engineering problems greatly. In order to overcome this defect, a hybrid multidimensional Riemann solver, called HMTHS (Hybrid of MulTv and multidimensional HLL scheme based on Self-similar structures), is proposed. It simulates the strongly interacting zone by adopting the MHLLES (Multidimensional Harten-Lax-van Leer-Eifeldt scheme based on Self-similar structures) scheme at subsonic speeds, which is with a high resolution by considering the second moment in the similarity variables. Also, it adopts the MULTV (Multidimensional Toro and Vasquez) scheme, which is with a high resolution in capturing discontinuities, to simulate the flux at supersonic speeds. Systematic numerical experiments, including both one-dimensional cases and two-dimensional cases, are conducted. One-dimensional moving contact discontinuity case and sod shock tube case suggest that HMTHS can accurately capture one-dimensional expansion waves, shock waves, and linear contact discontinuities. Two-dimensional cases, such as the double Mach reflection case, the supersonic shock / boundary layer interaction case, the hypersonic flow over the cylinder case, and the hypersonic viscous flow over the double-ellipsoid case, indicate that the HMTHS scheme is with a high resolution in simulating multidimensional complex flows. Therefore, it is promising to be widely applied in both scholar and engineering areas.     

降落伞数值模拟中的计算流体力学研究进展简介

文章介绍了降落伞数值模拟的基本思想。着重对降落伞数值模拟中的计算流体力学研究进展进行了综述 ,还对近年来的降落伞数值模拟应用进行了总结 ,突出了近两年来的发展变化。最后提出了降落伞流场的数值模拟尚需改进的若干方面     

旋翼螺旋桨/机翼巡航状态气动干扰规律研究

为优化倾转旋翼机气动布局,研究旋翼螺旋桨与机翼间复杂的气动干扰,在考虑有限翼展三维效应、翼尖小翼情况下,通过改变机翼外翼段的长度,研究了旋翼螺旋桨与机翼间的气动干扰规律.使用计算流体力学的三维非定常计算方法获得了巡航状态下旋翼螺旋桨与机翼的流场特征及气动特性参数,研究了螺旋桨滑流对机翼的上洗、下洗作用,以及机翼对螺旋桨...     

基于强化学习的多无人机避碰计算制导方法

针对大量固定翼无人机在有限空域内的协同避碰问题,提出了一种基于多智能体深度强化学习的计算制导方法.首先,将避碰制导过程抽象为序列决策问题,通过马尔可夫博弈理论对其进行数学描述.然后提出了一种基于深度神经网络技术的自主避碰制导决策方法,该网络使用改进的Actor-Critic模型进行训练,设计了实现该方法的机器学习架构,...     

Effect of Coriolis and centrifugal forces on flow and heat transfer at high rotation number and high density ratio in non orthogonally internal cooling channel

Numerical predictions of three-dimensional flow and heat transfer are performed for a two-pass square channel with 45° staggered ribs in non-orthogonally mode-rotation using the second moment closure model. At Reynolds number of 25,000, the rotation numbers studied were 0,0.24, 0.35 and 1.00. The density ratios were 0.13, 0.23 and 0.50. The results show that at high buoyancy parameter and high rotation number with a low density ratio, the flow in the first passage is governed by the secondary flow induced by the rotation whereas the secondary flow induced by the skewed ribs was almost distorted. As a result the heat transfer rate is enhanced on both co-trailing and co-leading sides compared to low and medium rotation number. In contrast, for the second passage, the rotation slightly reduces the heat transfer rate on co-leading side at high rotation number with a low density ratio and degrades it significantly on both co-trailing and co-leading sides at high buoyancy parameter compared to the stationary, low and medium rotation numbers. The numerical results are in fair agreement with available experimental data in the bend region and the second passage, while in the first passage were overestimated at low and medium rotation numbers.     

压气机旋转失速的一个气动声学模型

 根据广义Lighthill方程建立了压气机旋转失速的气动声学模型。提出了一种应用声学方法获取压气机旋转失速特征参数(旋转失速传播速度、数目)的新途径并结合试验结果进行了具体的分析。此外,在简化假设下,用数值法分析了旋转失速特征参数对声音幅值变化的影响。