The flow structure near the windward side of V-shaped wings with an attached shock wave on the leading edges

The results of a numerical calculation of a symmetric flow of supersonic gas with the Mach number M=3 past the windward side of V-shaped wings with an opening angle =40° and apex angles =30, 45, and 90° are given. The possibility of the ascent of one or two Ferri points from the break point of the transverse contour of the wing is discovered and explained. It is shown that conical flow near wings of finite length need not exist in flow regimes corresponding to angles of attack at which a Ferri point ascends, while at angles of attack smaller and larger than a certain interval, conical flow will exist. The investigation is conducted by means of a numerical method of stabilization with an artificial viscosity. The longitudinal coordinate, relative to which the steady system of equations is hyperbolic, played the part of the time variable, usual for methods of stabilization. The numerical method constructed using the scheme of [1] is described in [2] and was successfully applied to the calculation of different regimes of supersonic flow past conical wings with supersonic leading edges [2–6]. In. the present investigation the calculation algorithm of [2] is modified and makes it possible to realize motion with respect to the parameter a, this being particularly important for the stabilization of the solution in the calculation of flow regimes for which regions with a total velocity Mach number close to unity arise in the flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–131, January–February, 1986.     

Effect of Heat Supply on the Stability of the Supersonic Swept Atiachment-Line Boundary Layer

The stability of a boundary layer with volume heat supply on the attachment line of a swept wing is investigated within the framework of the linear theory at supersonic inviscid-free-stream Mach numbers. The results of numerical calculations of the flow stability and neutral curves are presented for the flow on the leading edge of a swept wing with a swept angle χ=60° at various free-stream Mach numbers. The effect of volume heat supply on the characteristics of boundary layer stability on the attachment line is studied at a surface temperature equal to the temperature of the external inviscid flow. It is shown that in the case of a supersonic external inviscid flow volume heat supply may result in an increase in the critical Reynolds number and stabilization of disturbances corresponding to large wave numbers. For certain energy supply parameters the situation is reversed, the unstable disturbances corresponding to the main flow-instability zone are stabilized but another zone of flow-instability with small wave numbers and a significantly lower critical Reynolds number appears.     

Numerical simulation for nonstationary Mach reflection of a shock wave: A kinetic-model approach

A numerical simulation was performed for the process of formation of single Mach reflection on a wedge by solving a BGK type kinetic equation for the reduced distribution function with a finite difference scheme. The calculations were carried out for a shock Mach number 2.75 and wedge angle 25° in a monatomic gas, which corresponds to the conditions of single Mach reflection in the classical von Neumann theory. The calculations were performed for both diffuse and specular reflection of molecules at the wall surface. It is concluded that the diffuse reflection of molecules at the wall surface or the existence of the viscous or thermal layer is an essential factor for a nonstationary process at the initial stage of Mach reflection. Furthermore, the numerical results for diffuse reflection are found to simulate the experimental results very well, such as a transient process from regular reflection to Mach reflection along with shock propagation.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1990.     

Velocity field near a wing—Fuselage combination at an angle of attack in a supersonic stream

To investigate interference between the wing and fuselage at supersonic flight velocities, one can, besides numerical methods based on the exact equations of motion, make effective use of the theory of small perturbations [1]. This is the direction adopted, in particular, in [2–4], in which the problem is solved in the framework of linear theory. In [5], the results obtained in the first approximations are corrected by taking into account the following term in the expansion of the potential function in a series in a small parameter. The present paper considers the velocity field near an arbitrarily profiled wing with supersonic edges and the features due to the presence of the fuselage. A general expression is found for the singular term of the asymptotic expansion of the solution of the linear equation in the neighborhood of the Mach cone with apex at the point of intersection of the leading edge of the wing with the surface of the fuselage. A uniformly exact solution for the linear differential equation for the additional velocity potential is constructed. The position and intensity of the shock wave on the upper surface of the wing are determined. Analytic dependences and quantitiative estimates are obtained for the local downwashes below the wing in the region of the flow where the linear theory leads to the largest errors. The obtained results are important for the correct determination of the aerodynamic characteristics of aircraft in the three-dimensional velocity field produced by the wing-fuselage combination.Translated from Izvestlya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 136–148, November–December, 1980.I am grateful to M. F. Pritulo for discussing the results of the work.     

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.     

Numerical analysis on flutter of Busemann-type supersonic biplane airfoil

The Busemann-type supersonic biplane can effectively reduce the wave drag through shock interference effect between airfoils. However, considering the elastic property of the wing structure, the vibration of the wings can cause the shock oscillation between the biplane, which may result in relative aeroelastic problems of the wing. In this research, fluid–structure interaction characteristics of the Busemann-type supersonic biplane at its design condition have been studied. A theoretical two-dimensional structure model has been established to consider the main elastic characteristics of the wing structure. Coupled with unsteady Navier–Stokes equations, the fluid–structure dynamic system of the supersonic biplane is studied through the two-way computational fluid dynamics/computational structural dynamics (CFD/CSD) coupling method. The biplane system has been simulated at its design Mach number with different nondimensional velocities. Different initial disturbance has been applied to excite the system and the effects of the position of the mass center on the system’s aeroelastic stability is also discussed. The results reveal that the stability of the airfoil in supersonic biplane system is decreased compared with that of the airfoil isolated in supersonic flow and such stability reduction effect should be given due attention in practical design.     

Shock wave solution of the Boltzmann kinetic equation in a 13-moment approximation

We consider the classic problem of a one-dimensional steady shock-wave solution of the Boltzmann kinetic equation utilizing a new type of 13-moment approximation proposed by Oguchi (1997). The model, unlike previous ones, expresses the collision term in an explicit function of the molecular velocity. This enables us to examine directly the nature of the singularity of the distribution function to this particular problem caused by the vanishing molecular velocity. We can thus obtain moment integrals directly because of its explicit expression. The principal value is utilized for the moment integral to cope with the singularity, and we can have five relations for five unknown functions to be determined with respect to the coordinate x. These relations can be reduced to a first-order differential equation that is solved to provide the familiar smooth monotonic transition from the upstream supersonic state to the subsonic downstream state. Computed values of shock thickness for various shock Mach numbers agree well with existing results obtained by different methods to the certain Mach number beyond which no solution exists.Received: 17 May 2002, Accepted: 1 May 2003, Published online: 15 August 2003PACS: 51.10. + y     

Wing of simple planform in a sharp-edged gust

The entry of a wing into the zone of a sharp-edged gust is considered in the linear formulation. The case is studied when the wing velocity is supersonic and its edges satisfy the supersonic flow condition. The gust intensity is considered to be variable, and its edge may move into the undisturbed medium. Equations in finite form are obtained for the forces and moments for a rectangular wing of infinite span, and also for triangular wings with positive and negative sweep, for the case when the gust intensity varies linearly. Sudden envelopment of the wing and penetration of the wing into a gust whose edge is fixed relative to the undisturbed medium are considered.     

Peak heating for reattachment of separated flow on a spiked blunt-body

A forward facing spike attached to a hemisphere-cylinder reduces the aerodynamic drag and the heat flux at supersonic and hypersonic Mach numbers. A numerical simulation is carried out to examine the effects of freestream Mach number on the flowfield and the heat transfer over the spiked blunt-body. Axisymmetric compressible Navier–Stokes equations are solved using a finite volume discretization in conjunction with a multistage Runge–Kutta time stepping method. Lengths of the separated region on the spike are influenced by the freestream Mach number. The computed results show that the peak heat flux on the nose of the blunt-body is also influenced by the freestream Mach number. The peak pressure and the wall heat flux on the blunt-body increases with increasing freestream Mach number. The computed results are reasonable in agreement with experimental data from the literature. Received on 1 July 1999     

Velocity field generated by wing vibrations propagating along the surface

The velocity field generated by wing vibrations propagating along an elastic wing surface with finite velocity is studied.The gasdynamic problem is reduced to a mixed boundary-value problem with a moving boundary for the three-dimensional wave equation. The solution is obtained in closed form when the wing travels at supersonic velocity following an arbitrary law, the vibration propagation front is an arbitrary curve displacing along the wing surface, and the wing edges are supersonic.     

Evolution of Vortex Structures on the Leeward Side of a Delta Wing

Results of an experimental study of a supersonic flow around the leeward side of a delta wing are presented. The experiments are performed on three delta wings with leading–edge sweep angles = 68°, 73°, and 78° for Mach numbers M =2—4 and angles of attack = 0—22°. Data on the structure and position of internal shock waves are obtained; the size and location of primary and secondary vortices are found. New regimes of the flow around a delta wing are identified. The chart of flow regimes around delta wings is refined and extended.     

Infrared thermography investigations in transitional supersonic boundary layers

The study of boundary-layer transition in supersonic flows is conducted employing infrared thermography (IRT). Several models of swept wings are tested in a blow-down facility at Mach number 2.4. The effects of wing sweep and other parameters (angle of attack, leading-edge contour, presence/absence of surface roughness) are successfully observed. The transition front is clearly identified, demonstrating the utility of IRT for this type of study. The technique is particularly indicated for flows that are sensitive to surface alterations (roughness), such as transitional boundary layers, because it does not require interaction with the model or the flow under investigation. The additional advantage of no need for special apparatus, except for the infrared camera, makes IRT well suited for both wind-tunnel and in-flight testing. Practical problems and limitations encountered when dealing with IRT in high-speed flows are also discussed.     

适用于几何非线性气动弹性分析的结构动力学降阶模型方法研究

几何非线性是壁板颤振和大展弦比机翼气动弹性等问题的一个主要特征,在进行数值仿真分析时往往需要采用商业非线性有限元求解器,存在计算量大和耦合迭代策略不易控制等问题。本文发展了一种适用于几何非线性的结构动力学降阶模型(CSD-ROM),利用广义坐标的非线性多项式表征非线性内力,采用参数识别方法获取多项式系数,并通过增加额外的线性模态来改善模型预测精度。基于此方法,分别针对壁板颤振、切尖三角翼的CFD/CSD-ROM非线性颤振问题开展了时域响应分析。计算结果表明,通过CSD-ROM计算出的壁板颤振速度为590 m/s,颤振频率为174 Hz,与有限元结果误差分别为0.8%和1.7%。马赫数0.879时切尖三角翼的颤振动压预测结果为2.25 psi,与非线性有限元相比的误差为3.8%。本文采用的非线性和线性模态基底组合方法,在保证计算精度的基础上可有效降低训练样本数量,一定程度上可替代非线性有限元开展气动弹性分析。     

Submerged shocks in conical gas flows in the presence of a mach shock-wave configuration

The results of a numerical study of a new type of singularities in the Mach shock-wave structure realized in supersonic nonsymmetric conical flows over V-wings with a bow shock attached to the leading edges are presented. Within the framework of the ideal gas model we study the changes in the shock system on transition, with increase in the sweep angle, from the region of nonsymmetric Mach interaction of the shocks attached to the leading edges of the wing to the region of special flow patterns, where on the windward cantilever surface a rarefaction flow is realized rather than a flow with an internal shock. It is shown, in particular, that in the region with special wing flow patterns a Mach system of shocks with a submerged shock proceeding from the branch point above the windward cantilever may exist.     

Stationary flow past the lower surface of a piecewise planar delta wing with supersonic leading edges

The considered wing has any finite number of inflections in its plane with lines of inflection intersecting at the point of inflection of the leading edge. In the present paper, this generalizes the author's earlier work [1] on flow past the undersurface of a flat wing at unite angle of attack with finite angle of slip and supersonic leading edges. In [1], calculations were not given. The special case of flow without slip in the same situation was considered later in [2], However, this paper contains errors, indicated at the end of the present paper. The calculations given in [2] are not correct. In the quoted papers, the gas flow is assumed to be a perturbation of a homogeneous flow behind a plane oblique shock wave. Such flows are treated systematically in [3]. Here and in [1], we use and generalize the representation of the linearized conservation laws across the shock front as the conditions of a boundary-value problem for an analytic function of a complex variable as obtained in [4, 5]. Calculations are given of the pressure distribution over the span for a number of different flow regimes and the pressure coefficients in the middle of the wing are compared with a numerical solution presented partly in [6].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 80–90, September–October, 1979.I am very grateful to V. I. Lapygin for making available a large number of variants of his numerical solution, and to L. E. Pekurovskii for assistance in the calculations.     

Modelling of dynamical crack propagation using time-domain boundary integral equations

We study dynamic antiplane cracks in the time domain by the boundary integral equation method (BIEM) based on the integral equation for displacement discontinuity (or crack opening displacement, COD) as a function of stress on the crack. This displacement discontinuity formulation presents the advantage, with respect to methods developed by Das and others in seismology, that it has to be solved only inside the crack. This BIEM is, however, difficult to implement numerically because of the hypersingularity of the kernel of the integral equation. Hence it is rewritten into a weakly singular form using a regularization technique proposed by Bonnet. The first step, following a method due to Sladek and Sladek, consists in converting the hypersingular integral equation for the displacement discontinuity into an integral equation for the displacement discontinuity and its tangential derivatives (dislocation density distribution); the latter involves a Cauchy type singular kernel. The second step is based on the observation that the hypersingularity is related to the static component of the kernel; the static singularity is then isolated and can be expressed in terms of weakly singular integrals using a result due to Bonnet. Although numerical applications discussed in this paper are all for the antiplane problem, the technique can be applied as well to in-plane crack dynamics.

The BIEM is implemented numerically using continuous linear space-time base functions to model the COD on the crack. In the present scheme the COD gradient interpolation is discontinuous at the element nodes while the integral equations are collocated at the element midpoints. This leads to an overdetermined discrete problem which is solved by standard least-squares methods. We use the dynamic BIEM to study a set of problems that appear in earthquake source dynamics, including the spontaneous dynamic crack propagation for a very simple rupture criterion. The numerical results compare favorably with the few exact solutions that are available. Then we demonstrate that difficulties experienced with finite difference simulations of spontaneous crack dynamics can be removed with the use of BIEM. The results are improved by the use of singular crack tip elements.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程，结合广义Betti-Rayleigh互易等式与时域下的边界积分方程，推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解，将经过主部分析的积分核函数分解为静态和动态部分，其中动态积分核不具有奇异性。在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换，采用配置点法计算超奇异积分，获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序，得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律，数值结果稳定且收敛速度快。     

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann–Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integration method can be developed and adopted to attack complex flows with different Mach numbers. HPF parallel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarily with massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuille- channel flow and pressure-driven gas flows in two-dimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of micro-scale gas flows occuring in the Micro-Electro-Mechanical System (MEMS). The project supported by the National Natural Science Foundation of China (90205009 and 10321002), and the National Parallel Computing Center in Beijing. The English text was polished by Yunming Chen.     

Nondimensional frequency scaling of aerodynamically-tensioned membranes

Membrane wings have applications that involve low Reynolds number flyers such as micro air vehicles. The time-averaged and time-dependent deformations of the membrane affect the aerodynamic characteristics of the wing, primarily in the region beyond the maximum aerodynamic efficiency of the wing. This paper investigates an appropriate nondimensional vibration frequency scaling of a spanwise tensioned membrane with free (unattached) leading and trailing edges at low Reynolds numbers relative to nondimensional aeroelastic parameters. Silicone rubber membranes with varying spanwise pre-tension, aerodynamic tension (due to wing angle-of-attack and flow dynamic pressure), modulus of elasticity, span, and thickness are studied. Experimental results are compared to a proposed scaling that simplifies the aerodynamic loading as a uniform pressure distribution acting on the membrane. Data is further compared and discussed relative to previous published results of membrane wings with finite wing spans (three-dimensional flow) and fixed (rigid) leading edges.     

跨音速后掠翼抖振边界计算

本文提出一种确定跨音速后掠翼抖振边界的数值计算方法,现有的确定跨音速翼型抖振边界的F.Thomas 准则被推广到包括具有大后掠角的后掠翼,计算是对侧滑翼进行,其中用积分法对三维可压缩湍流边界层的计算是根据本文作者听发展的方法,对于跨音速压强分布是利用A.Eberle的解全速位方程的有限元素法给出,按本文方法计算出F-86A 飞机的抖振边界与相同雷诺数下飞行试验所得结果符合得很好。