Calculation of the flow of an ideal fluid past a wing of finite thickness

The problem of irrotational flow past a wing of finite thickness and finite span can be reduced by Green's formula to the solution of a system of Fredholm equations of the second kind on the surface of the wing [1]. The wake vortex sheet is represented by a free vortex surface. Besides panel methods (see, for example, [2]) there are also methods of approximate solution of this problem based on a preliminary discretization of the solution along the span of the wing in which the two-dimensional integral equations are reduced to a system of one-dimensional integral equations [1], for which numerical methods of solution have already been developed [3–6]. At the same time, a discretization is also realized for the wake vortex sheet along the span of the wing. In the present paper, this idea of numerical solution of the problem of irrotational flow past a wing of finite span is realized on the basis of an approximation of the unknown functions which is piecewise linear along the span. The wake vortex sheet is represented by vortex filaments [7] in the nonlinear problem. In the linear problem, the sheet is represented both by vortex filaments and by a vortex surface. Examples are given of an aerodynamic calculation for sweptback wings of finite thickness with a constriction, and the results of the calculation are also compared with experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 124–131, October–December, 1981.     

Calculation of unsteady supersonic flow past a two-dimensional cascade of plates ween vorticity inhomogeneities of the flow act on it

In the framework of the linear theory of small perturbations the problem of unsteady subsonic flow past a two-dimensional cascade of plates has been considered in a number of papers. Thus, the unsteady aerodynamic characteristics of a cascade of vibrating plates were calculated in [1] by the method of integral equations, while the same method was used in [2, 3] to calculate the sound fields that are excited when sound waves Coming from outside or vorticity inhomogeneities of the oncoming flow act on the cascade. The problem of a two-dimensional cascade of vibrating plates in a supersonic flow was solved in [4, 5]. In [4] the solution was constructed on the basis of the well-known solution of the problem of vibrations of a single plate, while in [5] a variant of the method of integral equations was used which differed slightly from the usual formulation of this method [1–3]. The approach proposed in [5] is used below to calculate the unsteady flow past a two-dimensional cascade of plates in the case when vorticity inhomogeneities of a supersonic oncoming flow act on it. Equations are obtained for the strength of the unsteady pressure jumps arising in such a flow and the vortex wakes shed from the trailing edges of the plates. Examples of the calculations illustrating the accuracy of the method and its possibilities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp, 152–160, May–June, 1986.     

Plane incompressible fluid flow past an array of arbitrary profiles vibrating with arbitrary phase shift

We consider the problem of the vibration of an array of arbitrary profiles with arbitrary phase shift. Account is taken of the influence of the vortex wakes. The vibration amplitude is assumed to be small. The problem reduces to a system of two integral Fredholm equations of the second kind, which are solved on a digital computer. An example calculation is made for an array of arbitrary form.A large number of studies have considered unsteady flow past an array of profiles. Most authors either solve the problem for thin and slightly curved profiles or they consider the flow past arrays of thin curvilinear profiles [1].In [2] a study is made of the flow past an array of profiles of arbitrary form oscillating with arbitrary phase shift in the quasi-stationary formulation. The results are reduced to numerical values. Other approaches to the solution of the problem of unsteady flow past an array of profiles of finite thickness are presented in [3–5] (the absence of numerical calculations in [3, 4] makes it impossible to evaluate the effectiveness of these methods, while in [5] the calculation is made for a symmetric profile in the quasi-stationary formulation).     

Calculation of flow around a lifting surface in the case of large deformations

Many studies have been made of the nonstationary flow of an ideal incompressible fluid around a lifting surface. The present state of the numerical methods of solution of this problem is reviewed in [1]. The present paper studies three-dimensional nonstationary flow around a lifting surface which undergoes deformation and behind which a wake vortex surface is formed. The lifting and wake vortex surfaces are represented in parametric form. The metrics of these surfaces are used, and the introduced vortex function is approximated by bicubic splines. For the convenient application of the theory developed here to the flapping flight of insects, for which it is sometimes difficult to distinguish the lateral and trailing edges of the wings, the following terminology is introduced. The part of the edge of the lifting surface from which the wake vortex surface is shed is called the trailing edge. The remaining part is called the leading edge. On the leading edge, the velocity has a singularity. Test calculations have demonstrated the effectiveness of the method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 72–79, July–August, 1980.     

Mathematical Modeling of Time-Dependent Aerodynamic Loads on an Airfoil with Dynamic Deflection of a Control on the Trailing Edge

The linear problem of the time-dependent inviscid flow past a thin symmetric airfoil with a control on its trailing edge deflected in accordance with an arbitrary law is considered. The aerodynamic loads on the airfoil are calculated. The intensity of the vortex wake shed from the airfoil is determined by numerically solving a Volterra integral equation of the first kind. Questions of the mathematical modeling of the time-dependent aerodynamic loads in a form convenient for the joint solution of the problems of aerodynamics and flight dynamics are also considered. The results of the modeling are compared with the numerical solutions obtained.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, 2005, pp. 157–169.Original Russian Text Copyright © 2005 by Khrabrov.     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.     

Asymptotic solution to the problem of viscous flow in the neighborhood of the axis of a vortex sheet

An asymptotic solution is constructed to the problem of the flow of a viscous incompressible fluid in the neighborhood of the axis of a vortex sheet generated by flow separation from sharp edges of a delta wing of small aspect ratio at large values of the Reynolds number and small angles of attack.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 57–65, January–February, 1984.     

Flow of a vortical stream in the neighborhood of the critical point

A study was made of the axisymmetric flow of a viscous incompressible fluid in the neighborhood of the critical point of an obstacle when steady-state vortices oriented in the direction of the angular coordinate are introduced into the oncoming flow. A solution is presented of the equation for the transfer of a vortex in the case of an external flow containing a single largesize vortex in the low-frequency part of the spectrum. Using a finite integral Hankel transform, the problem is reduced to the solution of a system of ordinary differential equations. It is shown that a sufficiently large-size vortex can have a considerable effect on the structure of viscous flow near an obstacle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 85–89, November–December, 1970.     

Vortex flows of a viscous heat-conducting gas in a slightly divergent tube with volume energy supply

The results of a numerical investigation of viscous vortex flow in a slightly divergent tube with thermal energy supplied to the flow are presented. The initial stage of vortex flow development is considered for two different longitudinal velocity distributions simulating the velocity profiles in jet-like and wake-like vortex flows in the vicinity of the vortex axis. The first type of flow can be considered as a model for the near-axis region of the vortex formed in the flow around a delta wing at incidence. The second type can serve as a model for the near-axis region of the trailing vortex downstream of a high-aspect-ratio wing. The development of the two flows is studied for a constant area tube, a slightly divergent tube, and in the case of thermal energy supply from a volume energy source at a constant wall temperature.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 90–97, September–October, 1996.     

Hypersonic flow past a wing at large angles of attack with detached shock wave

The thin shock layer method [1–3] has been used to solve the problem of hypersonic flow past the windward surface of a delta wing at large angles of attack, when the shock wave is detached from the leading edge (but attached to the apex of the wing) and the velocity of the gas in the shock layer is of the same order as the speed of sound. A classification of the regimes of flow past a delta wing at large angles of attack has been made. A general solution has been obtained for the problem of three-dimensional hypersonic flow past the wing allowing for nonequilibrium physicochemical processes of thermal radiation of the gas at high temperatures.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 149–157, May–June, 1985.     

The damped natural oscillations of a gas flowing past a cascade of flat plates

We solve the problem of the natural oscillations of a gas flowing past a cascade of flat plates under the Joukowsky-Chaplygin condition that the velocity at the trailing edge of the profiles is finite. In this case part of the energy of the oscillating gas is consumed in the formation of a trailing vortex. The corresponding eigenvalues of the problem are complex and so the natural oscillations of the gas are damped. The computational results are compared with the results of experimental investigation of acoustic resonance in flow past a cascade of flat plates obtained by Parker [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 84–88, September–October, 1970.     

A double vortex sheet model of the viscous flow near the trailing edge of a lifting aerofoil

Classical definitions of boundary layer mass and momentum flux deficiency thicknesses can lead to gross errors when applied to measurements near a trailing edge where the flow curvature in the free stream is appreciable. This paper presents a double vortex sheet model as a development from the single vortex sheet model of Helmholtz and others. Two bound vortex sheets define a potential function which can describe a flow with the same mass and momentum flux deficiencies as the viscous regions. The bound nature of these sheets allows the modelling of the integral properties of these regions while retaining the advantages of a potential flow. The application to the flow near the trailing edge of a lifting aerofoil is given     

On the nonlinear theory of the wing in a plane unsteady flow

The main aspects of the nonlinear theory of the wing in a plane unsteady fluid flow are generalized on the basis of the author’s previous results. An initial-boundary problem for complex velocity is formulated. A system of differential equations with conditions at points of vortex wake shedding is presented, which allows a large class of problems to be solved correctly. The Cauchy problem is solved by using a standard discretization procedure. The boundary-value problem is reduced at each time step to singular integral equations of the first and second kind. The accuracy of solving these equations by the method of discrete vortices and by the method of panels is compared. Specific features of pressure calculations in the case of a separated flow around the airfoil contour are discussed     

Contribution to the theory of thin-profile trailing edge separation

The conditions of nonsymmetric trailing edge flow with separation are investigated. Solutions of the equations for the interaction zone in the neighborhood of the trailing edge of a thin profile at an angle of attack of the order O(Re^–1/16) in the separated flow regime are constructed numerically. It is shown that for this zone a solution exists up to a certain angle of attack. In all the regimes the value of the friction on the upper surface at the very end of the trailing edge remains a positive quantity. The solution of the equations in the separated flow regimes is found to be nonunique. The flow over the leading edge is assumed to be unseparated, and the separation at the trailing edge, if present, is assumed to be localized in the interior of the boundary layer. The flow over a Kutta profile at zero angle of attack is taken as an example. In this case the satisfaction of the Chaplygin-Joukowsky condition at the trailing edge ensures smooth flow over both the trailing and leading edges.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 55–59, July–August, 1989.     

Some Specific Features of Integral Equations with the Cauchy Kernel on a Closed Contour in Hydrodynamic Problems

Singular integral equations of the first and second kind with the Cauchy kernel on a limiting narrow closed contour are theoretically considered. The initial equations are found to become different on the limiting contour. This singularity of integral equations with the Cauchy kernel does not allow boundary-value problems of the flow around thin airfoils to be solved correctly; therefore, a system consisting of integral equations of the first and second kind is proposed for solving such problems. The results of the present study are tested against an exact solution of the problem of the flow past a flat plate.     

An inviscid model for vortex shedding from a deforming body

An inviscid vortex sheet model is developed in order to study the unsteady separated flow past a two-dimensional deforming body which moves with a prescribed motion in an otherwise quiescent fluid. Following Jones (J Fluid Mech 496, 405–441, 2003) the flow is assumed to comprise of a bound vortex sheet attached to the body and two separate vortex sheets originating at the edges. The complex conjugate velocity potential is expressed explicitly in terms of the bound vortex sheet strength and the edge circulations through a boundary integral representation. It is shown that Kelvin’s circulation theorem, along with the conditions of continuity of the normal velocity across the body and the boundedness of the velocity field, yields a coupled system of equations for the unknown bound vortex sheet strength and the edge circulations. A general numerical treatment is developed for the singular principal value integrals arising in the solution procedure. The model is validated against the results of Jones (J Fluid Mech 496, 405–441, 2003) for computations involving a rigid flat plate and is subsequently applied to the flapping foil experiments of Heathcote et al. (AIAA J, 42, 2196–2204, 2004) in order to predict the thrust coefficient. The utility of the model in simulating aquatic locomotion is also demonstrated, with vortex shedding suppressed at the leading edge of the swimming body.      

Direct and inverse problems of flow over a wing of finite span in the linear formulation

A correspondence between the solutions of the direct and the inverse problem for wing theory is established for a wing of finite span in the framework of linear theory on the basis of solution of a wave equation in Volterra form for supersonic flow and solution of the Laplace equation in the form of Green's formula for subsonic flow. For the direct problem in the case of supersonic flow an expression is derived for finding the load on the wing with maximal allowance for the wing geometry. In the inverse problem for supersonic and subsonic flows, expressions are derived for finding the wing geometry from given values of the load on the wing and the variation of the load along the span of the wing. The solution of the inverse problem is presented in the form of integrals that converge for interior points of the wing surface in the sense of the Cauchy principal value, the wing surface being represented as a vortex surface of mutually orthogonal vortex lines. The conditions of finiteness of the velocities on the edges are discussed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 114–125, September–October, 1979.     

Time marching integral equation method for unsteady transonic flows

In this paper, we have proposed a time marching intregral equation method which does not have the limitation of the time linearized integral equation method in that the latter method can not satisfactorily simulate the shock-wave motions. Firstly, a model problem—one dimensional initial and boundary value wave problem is treated to clarify the basic idea of the new method. Then the method is implemented for 2-D and 3-D unsteady transonic flow problems. The introduction of the concept of a quasi-velocity-potential simplifies the time marching integral equations and the treatment of trailing vortex sheet condition. The numerical calculations show that the method is reasonable and reliable.     

Application of a vortex lattice numerical model in the calculation of inviscid incompressible flow around delta wings

The flow over a flat plate delta wing at incidence and in sideslip is studied using vortex lattice models based on streamwise penelling. For the attached flow problem the effect of sideslip is simulated by modifying the standard vortex lattice model for zero sideslip by aligning the trailing vortices aft of the wing along the resultant flow direction. For the separated flow problem a non-linear vortex lattice model is developed for both zero and non-zero sideslip angles in which the shape and position of the leading edge separation vortices are calculated by an iterative procedure starting from an assumed initial shape. The theoretical values are compared with available theoretical and experimental results.     

基于雨燕翅膀的仿生三角翼气动特性计算研究

针对低雷诺数微型飞行器的气动布局,设计出类似雨燕翅膀的一组具有不同前缘钝度的中等后掠(Λ=50?)仿生三角翼.为了定量对比研究三角翼后缘收缩产生的气动效应,设计了一组具有同等后掠的普通三角翼.为了深入研究仿生三角翼布局的前缘涡演化特性以及总体气动特性,采用数值模拟方法详细地探索了低雷诺数(Re=1.58×104)流动条...