Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been developed to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

Numerical Analyses of Discrete Gust Response for an Aircraft

Based on Navier-Stokes equations and structural and flight dynamic equations of motion, dynamic responses in vertical discrete gust flow perturbation are investigated for a supersonic transport model. A tightly coupled method was developed by subiterations between aerodynamic equations and dynamic equations of motion. First, under the assumption of rigid-body and single freedom of motion in the vertical plunging, the results of a direct-coupling method are compared with the results of quasi-steady model method. Then, gust responses for the one-minus-cosine gust profile arc analyzed with two freedoms of motion in plunging and pitching for the airplane configurations with and without the consideration of structural deformation.

跨音速翼型反设计的一种大范围收敛方法

求解跨音速翼型的反设计问题时,传统的梯度型方法一般均为局部收敛.为增大求解的收敛范围,依据同伦方法的思想,通过构造不动点同伦,将原问题的求解转化为其同伦函数的求解,并依据拟Sigmoid函数调整同伦参数以提高计算效率,进而构造出一种具有较高计算效率的大范围收敛反设计方法.数值算例以RAE2822翼型的表面压力分布为拟合目标,分别采用B样条方法,PARSEC方法及正交形函数方法等3种不同的参数化方法,并分别以NACA0012,OAF139及VR15翼型为初始翼型进行迭代计算.计算结果证明,该方法适用于多种参数化方法,且具有较好的计算效率,从多个不同的初始翼型出发,经较少次数迭代后,均能与目标翼型很好地拟合,是一种高效的大范围收敛方法.

Drag on an ellipsoid particle in free-molecular plasma flow

Based on the analysis of molecular gas dynamics, the drag and moment acting on an ellipsoid particle of revolution X-2/a(2) + Y-2/a(2) + Z(2)/c(2) = 1, as an example of nonspherical particles, are studied under the condition of free-molecular plasma flow with thin plasma sheaths. A nonzero moment which causes nonspherical particle self-oscillation and self-rotation around its own axis in the plasma flow-similar to the pitching moment in aerodynamics-is discovered for the first time. When the ratio of axis length c/a is unity, the moment is zero and the drag formula are reduced to the well-known results of spherical particles. The effects of the particle-plasma relative velocity, the plasma temperature, and the particle materials on the drag and moment are also investigated.

三角翼受迫俯仰滚转耦合运动的气动特性研究

基于刚性动网格的技术,选用B-L湍流模型,利用有限控制体积法对N-S方程进行数值离散,对76°大后掠三角翼的受迫俯仰滚转耦合运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析.计算结果表明:俯仰滚转耦合运动时,三角翼上表面的涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大.

Research progress on high-enthalpy and hypersonic flows

The research progress on high-enthalpy and hypersorlic flows having been achieved in the Institute of Mechanics, Chinese Academy of Sciences, is reported in this paper. The paper consists of three main parts: The first part is on the techniques to develop advanced hypersonic test facilities, in which the detonation-driven shock-reflected tunnel and the detonation-driven shock-expanded tube are introduced. The shock tunnel can be used for generating hypersonic flows of a Mach number ranging from 10 to 20, and the expansion tube is applicable to simulate the flows with a speed of 7 similar to 10km/s. The second part is dedicated to the shock tunnel nozzle flow diagnosis to examine properties of the hypersonic flows thus created. The third part is on experiments and numerical simulations. The experiments include measuring the aerodynamic pitching moment and heat transfer in hypersonic flows, and the numerical work reports nozzle flow simulations and flow non-equilibrium effects on the possible experiments that may be carried out on the above-mentioned hypersonic test facilities.

Experimental Demonstration Of A New Concept Of Drag Reduction And Thermal Protection For Hypersonic Vehicles

A new idea of drag reduction and thermal protection for hypersonic vehicles is proposed based on the combination of a physical spike and lateral jets for shock-reconstruction. The spike recasts the bow shock in front of a blunt body into a conical shock, and the lateral jets work to protect the spike tip from overheating and to push the conical shock away from the blunt body when a pitching angle exists during flight. Experiments are conducted in a hypersonic wind tunnel at a nominal Mach number of 6. It is demonstrated that the shock/shock interaction on the blunt body is avoided due to injection and the peak pressure at the reattachment point is reduced by 70% under a 4A degrees attack angle.

Euler Solution Using Adaptive Cartesian Grid With A Gridless Boundary Treatment

A quadtree-based adaptive Cartesian grid generator and flow solver were developed. The grid adaptation based on pressure or density gradient was performed and a gridless method based on the least-square fashion was used to treat the wall surface boundary condition, which is generally difficult to be handled for the common Cartesian grid. First, to validate the technique of grid adaptation, the benchmarks over a forward-facing step and double Mach reflection were computed. Second, the flows over the NACA 0012 airfoil and a two-element airfoil were calculated to validate the developed gridless method. The computational results indicate the developed method is reasonable for complex flows.

Numerical investigation on a transition prediction model

A new transition prediction model is introduced, which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition. Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil. Performance of this transition model for incompressible flows is studied, with numerical results consistent to experimental data. The requirement of grid resolution for this transition model is also studied.

压电纤维复合材料铺层用于翼面设计的驱动特性与刚度影响

压电纤维复合材料驱动器在形状控制、振动控制、颤振抑制与抖振控制等方面有广泛的应用前景. 首先简单介绍了压电应变驱动的比拟载荷方法, 并采用该方法讨论了压电陶瓷片状驱动器与压电纤维复合材料驱动器在驱动特性上的主要差异. 在此基础上, 对压电纤维复合材料在不同铺设方式、铺设角度与铺设层数下的驱动特性进行了分析, 在刚度影响方面展示了不同铺设角度下模型刚轴的移动. 分析结果表明: 对称铺设反向电场可以同时获得弯曲与扭转变形, 而反对称铺设同向电场主要获得扭转变形; 两种铺设方式下45°铺设角均获得最大弦向转角, 而0°铺设角将获得最大挠度; 多铺层可以增加驱动载荷, 但总体变形效果还取决于结构系统的刚度比例; 对称铺设方式下铺设角对结构刚轴移动的影响非常明显, 在气动弹性控制中应着重关注.

三维蛇形机器人巡视者Ⅱ的开发

蛇形机器人具有比传统移动机器人更强的运动能力,为实现机器人的三维运动而开发的蛇形机器人巡视者II是一个具有强驱动力和高机动性的三维蛇形机器人。它由具有3自由度的模块化单元组成,该单元具有俯仰、偏航和回转3自由度,单元的俯仰和偏航运动是通过由3个伞齿轮组成的差速机构耦合驱动。该单元的圆柱形外壳周围安装有一系列的被动轮来增加机器人运动灵活性。对蛇形机器人模块单元的自由度作了详细分析,并基于机器人的运动机动性设计三维蛇形机器人的单元结构。给出蛇形机器人的控制系统结构,并对耦合变量进行解耦。将蛇形机器人的蜿蜒运动和扭转运动两种基本步伐应用到巡视者II上,试验结果证明了该三维蛇形机器人具有很强的机动性和运动能力。