Near-field fluctuations and far-field noise of a three-element airfoil system by a discrete vortex method

A circulation-based discrete vortex method is used on a three-element airfoil system. Kutta conditions and Kelvin’s circulation theorem are additional conditions required for this method to determine the circulation distributions on each element and to determine vortex shedding. Discrete shed vortices are introduced near the four sharp edges to represent the sharp-edge vortex shedding caused by unsteady flow separation. The computational procedure warrants neutrally stable solutions of the self-sustained fluctuating flowfield that can provide broad-band spectral information for far-field noise predictions. The near-field vortex method directly calculates the parameters used in an asymptotic formula for far-field sound computation that attributes the noise sources to vortex interactions among the shed vortices and the surface circulations of the three-element airfoil system. The far-field noise characters are then analyzed and compared to the experimental data in the literature.     

翼涡干扰前缘开孔被动控制数值研究

开孔方法是一种简单的流动被动控制方法.为找到一种有效降低桨涡干扰效应的被动控制方法,以NACA 0012翼型作为研究对象,建立了4种前缘开孔的模型.在不同来流速度、涡的强度和干扰距离条件下,对4种前缘开孔模型和无孔的基准翼型进行了二维平行桨涡干扰（翼涡干扰）数值模拟,对比了升力系数的变化.结果表明:前缘开孔可以降低翼涡干扰效应,但对翼型升力系数有一定的影响;宽度为2.5%弦长的直孔能在翼型升力系数损失较小的情况下有效地降低翼涡干扰效应,且适用范围较广.     

Contribution to the propeller aerodynamic interference

Propeller wake can significantly change the flowfield at the downstream lifting surfaces and therefore influence its aerodynamic coefficients. The numerical model of the propeller presented here is using discrete vortices to form vortex sheet that is leaving each blade. Model is also applicable for combination of lifting surface and propeller using undeveloped propeller vortex sheet in determination of aerodynamic interference of the propeller on the downstream lifting surface, wing or tail for small angle of attack. This low computational cost numerical model is suitable for implementation in component build–up method used in preliminary estimation of aerodynamic coefficients for different propeller aircraft configurations.     

Simulation of Aerodynamic Performances of Flexible Flapping Wing Airfoils北大核心CSCD

与固定翼相比,在低速、小Reynolds数条件下,扑翼飞行具有显著的气动性能优势,受到越来越多的重视。然而,目前对扑翼翼型的研究以刚性翼型为主,对柔性翼型气动性能认识还不清楚。该文建立了柔性椭圆翼型的流固耦合仿真模型,分析了不同风速、迎角下柔性椭圆翼型的周围流场、变形以及气动性能。仿真结果表明,较刚性翼型,柔性翼型延缓了尾涡脱落时间,有效降低升力扰动振荡频率;柔性翼型显著抑制了尾流流场的扰动,降低升力扰动振荡幅值,合适的弹性模量翼型使得扰动振荡完全消除。研究结果可为软飞行器气动设计提供参考。     

Numerical simulation of laminar flow past a circular cylinder

The present paper focuses on the analysis of two- and three-dimensional flow past a circular cylinder in different laminar flow regimes. In this simulation, an implicit pressure-based finite volume method is used for time-accurate computation of incompressible flow using second order accurate convective flux discretisation schemes. The computation results are validated against measurement data for mean surface pressure, skin friction coefficients, the size and strength of the recirculating wake for the steady flow regime and also for the Strouhal frequency of vortex shedding and the mean and RMS amplitude of the fluctuating aerodynamic coefficients for the unsteady periodic flow regime. The complex three dimensional flow structure of the cylinder wake is also reasonably captured by the present prediction procedure.     

Numerical bifurcation analysis of static stall of airfoil and dynamic stall under unsteady perturbation

By the finite element method combined with Arbitrary-Lagrangian-Eulerian (ALE) frame and explicit Characteristic Based Split Scheme (CBS), the complex flows around stationary and sinusoidal pitching airfoil are studied numerically. In particular, the static and dynamic stalls are analyzed in detail, and the natures of the static stall of NACA0012 airfoil are given from viewpoint of bifurcations. Following the bifurcation in Map, the static stall is proved to be the result from saddle-node bifurcation which involves both the hysteresis and jumping phenomena, by introducing a Map and its Floquet multiplier, which is constructed in the numerical simulation of flow field and related to the lift of the airfoil. Further, because the saddle-node bifurcation is sensitive to imperfection or perturbation, the airfoil is then subjected to a perturbation which is a kind of sinusoidal pitching oscillation, and the flow structure and aerodynamic performance are studied numerically. The results show that the large-scale flow separation at the static stall on the airfoil surface can be removed or delayed feasibly, and the ensuing lift could be enhanced significantly and also the stalling incidence could be delayed effectively. As a conclusion, it can be drawn that the proper external excitation can be considered as a powerful control strategy for the stall. As an unsteady aerodynamic behavior of high angle of attack, the dynamic stall can be investigated from viewpoint of nonlinear dynamics, and there exists a rich variety of nonlinear phenomena, which are related to the lift enhancement and drag reduction.     

Numerical study of helicopter blade–vortex mechanism of interaction using the potential flow theory

The blade–vortex interaction (BVI) phenomenon plays a key role in the rotorcraft aerodynamics. Numerical investigations of BVI using classical CFD approaches are computationally expensive. In the present research we propose a numerical approach, based on the potential flow theory, for the numerical investigation of helicopter blade–vortex mechanism of interaction. This approach overcomes the computational expenses posed by the CFD techniques. The influence of vertical miss distance, angle of attack, airfoil camber, and vortex strength on the helicopter blade–vortex mechanism of interaction is subject of investigation. The study reveals that the magnitude of the aerodynamic coefficients decreases with the increase of vertical miss distance and angle of attack, and the decrease of vortex strength and core size.     

Computational fluid dynamics simulation of aerodynamic performance and flow separation by single element and slatted airfoils under rainfall conditions

This study investigated the effects of rainfall on flow separation and the aerodynamic performance of single element and slatted NACA 0012 airfoils by using a mathematical model developed with the commercial computational fluid dynamics solver ANSYS FLUENT 18.2. A two-way momentum coupled Eulerian–Lagrangian multiphase approach was used to simulate the formation of the water film layer on the airfoil's surface. According to the results, very low values of the lift-to-drag ratio at low angles of attack reflected severe degradation of the aerodynamic performance of the airfoil in the presence of water accumulated on its surface. The impact of rain droplets on the leading-edge slat surface led to less water accumulating on the main section of the airfoil. In particular, the maximum water film mass concentrated on the airfoil surface decreased from 15 g to 1 g compared with the single element airfoil. Hence, the thickness of the water film layer was not sufficiently large to significantly affect the aerodynamic coefficients of the slatted airfoil, especially the maximum lift coefficient, compared with the thicker water film layer on the single element airfoil. In addition, the use of slats clearly enhanced the aerodynamic coefficients and increased the stall angle from 13° to 22° in dry conditions, and from 16° to 24° in rainy conditions. Slats also significantly decreased the boundary layer thickness and delayed the separation at higher angles of attack.     

开敞式索膜结构的气动失稳临界风速

考虑几何非线性的影响,利用无限薄的旋涡薄层模拟气流在结构表面形成的扰动,由非稳态Bernoulli方程和环量定理将空气压力表示成旋涡密度的函数;然后由涡格法结合耦合边界条件求出旋涡密度,根据系统的稳定性判据得到结构发散失稳临界风速的解析表达式．通过三维开敞式膜结构的计算分析,发现空间膜结构的曲率是影响结构气动失稳临界风速的主要因素．     

Unsteady Flow and its Influence to Aerodynamic and Aeroacoustic of VAWT's

A study is reported of the influence of unsteady flow on the aerodynamics and aeroacoustics of vertical axis wind turbines by numerical simulation. The combination of aerodynamic predictions with a discrete vortex method and aeroacoustic predictions based on Ffowcs Williams-Hawkings equation is used to achieve this goal. The numerical results show that unsteady flow of the turbine has a significant influence on the turbine aerodynamics and can lead to a decrease in generated noise as compared to the conventional horizontal axis wind turbine at the similar aerodynamic performance. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical analysis of the dynamics of distributed vortex configurations

A numerical algorithm is proposed for analyzing the dynamics of distributed plane vortex configurations in an inviscid incompressible fluid. At every time step, the algorithm involves the computation of unsteady vortex flows, an analysis of the configuration structure with the help of heuristic criteria, the visualization of the distribution of marked particles and vorticity, the construction of streamlines of fluid particles, and the computation of the field of local Lyapunov exponents. The inviscid incompressible fluid dynamic equations are solved by applying a meshless vortex method. The algorithm is used to investigate the interaction of two and three identical distributed vortices with various initial positions in the flow region with and without the Coriolis force.     

Mathematical and numerical treatment of instabilities of non-axisymmetric confined vortices under the Dirichlet boundary conditions

This paper is concerned with an efficient analytical and numerical method for computing instability bounds of non-axisymmetric confined vortices. We propose a new approach based on a recently mathematical method of meshless collocation, when the shifted orthogonal Chebyshev base is adapted to satisfy the boundary conditions. Appropriately modified Tollmien–Schlichting modes used for the linear stability analysis offer the main advantage of dealing with a steady perturbation model of a confined vortex that leads to a variety of information describing the behavior of the full three-dimensional unsteady flow.     

An Improved Hybrid Adjoint Method in External Aerodynamics Using Variational Technique for the Boundary Integral Based Optimal Objective Function Gradient

An improved hybrid adjoint method to the viscous, compressible Reynold-Averaged Navier-Stokes Equation (RANS) is developed for the computation of objective function gradient and demonstrated for external aerodynamic design optimization. In this paper, the main idea is to extend the previous coupling of the discrete and continuous adjoint method by the grid-node coordinates variation technique for the computation of the variation in the gradients of flow variables. This approach in combination with the Jacobian matrices of flow fluxes refrained the objective function from field integrals and coordinate transformation matrix. Thus, it opens up the possibility of employing the hybrid adjoint method to evaluate the subsequent objective function gradient analogous to many shape parameters, comprises of only boundary integrals. This avoids the grid regeneration in the geometry for every surface perturbation in a structured and unstructured grid. Hence, this viable technique reduces the overall CPU cost. Moreover, the new hybrid adjoint method has been successfully applied to the computation of accurate sensitivity derivatives. Finally, for the investigation of the presented numerical method, simulations are carried out on NACA0012 airfoil in a transonic regime and its accuracy and effectiveness related to the new gradient equation have been verified with the Finite Difference Method (FDM). The analysis reveals that the presented methodology for the optimization provides the designer with an indispensable CPU-cost effective tool to reshape the complex geometry airfoil surfaces, useful relative to the state-of-the-art, in a less computing time.     

基于线性涡模型的风洞孔壁干扰特性分析

发展了一种适用于二元翼型试验洞壁干扰特性的评估和修正方法.基于Prandtl-Glauert速度势方程和布置在模型及洞壁表面的线性涡,采用迭代方法计算了风洞孔壁对翼型表面压力分布特性的影响,分析了不同孔壁透气特性参数的影响规律和量值,利用与国外参考结果及风洞试验结果的对比确定了该方法的准确性.结果表明,孔壁对翼型绕流的影响主要反映在上翼面吸力峰和最大厚度位置之间,使压力系数减小,积分后的升力系数降低,且随着孔壁透气特性参数的增大,洞壁干扰由实壁特性向开口特性发展,洞壁干扰、影响量急剧增大.与传统方法相比,该方法计算快速,结果可靠,同时具备试验前评估的能力,可用于亚临界范围内翼型表面压力的快速估算,以及翼型试验的洞壁干扰修正.     

径向基函数参数化翼型的气动力降阶模型优化

基于小扰动和弱非线性假设,提出了一种基于气动力降阶模型和径向基函数参数化的翼型优化方法.其主要方法是用径向基函数参数化翼型扰动;通过CFD辨识参数扰动对翼型气动力影响的降阶模型核函数;基于叠加法建立了参数变化对翼型气动力影响的降阶模型;最后基于该气动力降阶模型计算并优化翼型升阻特性.NACA0012翼型优化的结果表明基于气动力降阶模型的优化方法是可行的,可以极大地提高翼型优化速度.     

High performance computation for DNS/LES

The paper is focused on high-order compact schemes for direct numerical simulation (DNS) and large eddy simulation (LES) for flow separation, transition, tip vortex, and flow control. A discussion is given for several fundamental issues such as high quality grid generation, high-order schemes for curvilinear coordinates, the CFL condition for complex geometry, and high-order weighted compact schemes for shock capturing and shock–vortex interaction. The computation examples include DNS for K-type and H-type transition, DNS for flow separation and transition around an airfoil with attack angle, control of flow separation by using pulsed jets, and LES simulation for a tip vortex behind the juncture of a wing and flat plate. The computation also shows an almost linear growth in efficiency obtained by using multiple processors.     

Third-order accurate conservative method on unstructured meshes for gasdynamic simulations

A third-order accurate finite-volume method on unstructured meshes is proposed for solving viscous gasdynamic problems. The method is described as applied to the advection equation. The accuracy of the method is verified by computing the evolution of a vortex on meshes of various degrees of detail with variously shaped cells. Additionally, unsteady flows around a cylinder and a symmetric airfoil are computed. The numerical results are presented in the form of plots and tables.     

Discrete Adjoint based Sensitivity Analysis for Optimal Active Flow Control of High-Lift Configurations

Blowing and suction type of active flow control techniques can be used to delay the flow separation on the flap and to enhance the aerodynamic performance of high-lift configurations. Effective separation control and maximum enhancement in the mean lift coefficient are achieved by finding the optimal actuation parameters. The optimal set of actuation parameters can be obtained by combining the gradient based algorithms with discrete adjoints. In the present work, an unsteady discrete adjoint incompressible RANS solver is developed for the optimal active separation control. The adjoint solver is applied to the test case of active flow control on the flap of a 2D high-lift configuration. Sensitivity gradients are presented to demonstrate the accuracy of unsteady adjoint RANS solver. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of the unsteady aerodynamic interaction of two plane airfoil cascades

A method for the calculation of unsteady aerodynamic interaction of two plane airfoil cascades that are in relative motion in a subsonic flow of ideal gas is developed. This interaction provides a two-dimensional approximation of the flow in a stage of an axial turbomachine. The method is based on the reduction of the problem to the calculation of the unsteady flow in a single interblade passage of each of the cascades. The calculation uses generalized space-time periodicity relations corresponding to the unsteady process of interest. The calculation is based on the direct numerical integration of the non-stationary gas dynamics equations with the use of the finite difference Godunov-Kolgan-Rodionov scheme of the second approximation order with respect to time and space. The calculation procedure includes the determination of the acoustic fields that are generated by the stage in the incident flow and in the flow behind it. The results of the calculations that illustrate the accuracy of the numerical solution and the capabilities of the method are presented.