翼-身组合体绕流的Euler方程数值模拟

用代数方法生成翼-身组合体H-O型网格,并用有限体积法研制出翼-身组合体绕流三维Euler方程计算程序。该方法的特点是改进了机身与机翼表面网格点分布,机翼后缘有后掠时也能保证后缘与网格线一致。程序除能按常规提供横流截面展向压强分布外,还能提供弦向压强分布。对NASA TND-712的翼-身组合体模型的计算结果与实验符合很好。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法．在背景网格建立和三维初始推进面生成的前提下，利用推进阵面法生成了翼身组合体的三维空间网格，并应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

基于非结构网格的翼身组合体绕流数值模拟

应用投影变换的方法生成三维翼身组合体的表面网格，提出了新的三维空间多层背景网格生成方法，在背景网格建立和三维初始推进面生成的前提下．利用推进阵面法生成了翼身组合体的三维空间网格．应用网格光顺技术和变换方法，改进了初始生成网格的质量。应用格心格式的有限体积法对翼身组合体绕流进行了Euler方程的数值模拟，得到了较为满意的计算结果。     

无限插值方法生成三维网格

 本文研究了生成三维网格的无限插值理论,应用它构筑了绕翼身组合体的三维C-O型网格,这是一种代数方法,适用于不同拓扑结构的网格生成,具有程序编写简单、计算时间短的优点。在这种方法生成的网格上进行了二维与三维的数值实验,结果表明,它是一种有效的网格生成方法。     

翼-身组合体跨音速守恒型全位势流AF2高效算法的应用

 本文用Joukowski交换将圆截面机身的每一横截面逐一转绘为垂直切缝,因而将翼-身组合体转绘为某单独机翼。我们用作者采用过的数值保角转绘法对转换的单独机翼生成O型贴体坐标网格,然后用逆转绘构造翼-身组合体绕流的O型贴体坐标网格。跨音速计算采用守恒型金位势方程,精确边界条件和AF2高效有限差分迭代算法。本文指出,只需少量修改,就可将单独翼分析程序扩展为可同时适用于翼-身组合体。本文考虑的组合体其机身为有限或无限长,圆形截面,机翼前后机身轴线允许弯曲;机翼有任意平面形状,但翼尖弦长不能太短。     

两种导弹型的翼-身组合体的涡系干扰

通过染色法和氢气泡法流态显示技术观察了迎角从零至25°的两种翼-身组合体的涡系干扰现象,一种是十字翼-身组合体,另一种是×字翼-身组合体。并且对两种组合体涡系干扰的差异进行了分析比较和讨论。为了比较,对单独外露弹翼也进行了实验观测。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

翼—身—尾组合体绕流的Euler方程计算

本文采用重迭网格技术和Ｅｕｌｅｒ方程计算翼－身－尾组合体绕流。对翼－身与尾－身部分采用各自的Ｈ－Ｏ型网格。Ｅｕｌｅｒ方程求解采用Ｊａｍｅｓｏｎ有有限体积法，即中心差分挖和显式Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身－尾模型为例，计算的空气动力特性与实验符合较好。     

翼－身－尾组合体绕流的Euler方程计算

本文采用重迭网格技术和ＥＪｅｒ方程计算翼－身－尾组合体绕流。对算身与尾。身部分采用各自的Ｈ。Ｏ型网格，Ｅ＊ｅ／方程求解采用１ａｍｅｓｏｎ的有限体积法，即中心差分近似和显式Ｒｕｎｇｅ·Ｋｕｔｔａ时间推进。采用前后区交替迭代使前后两区通过重迭区交换信息。本文用ＮＡＣＡＴＮ４０４１翼－身，尾模型为例，计算的空气动力特性与实验符合较好。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

涡破裂诱导的垂尾抖振数值模拟

 采用数值模拟方法研究了大后掠三角翼前缘涡破裂诱导的垂尾抖振问题,分析了大迎角条件下的垂尾抖振特性。采用Navier-Stokes方程求解非定常气动力、耦合结构动力学方程,建立了气动弹性方程,在时域内采用松耦合方式推进以得到垂尾结构响应。研究结果表明:涡破裂流的脉动频带覆盖了垂尾扭转模态的固有频率,诱发了垂尾抖振现象;与传统的颤振频域响应特性不同,垂尾抖振响应的各阶位移与加速度响应主频均位于各阶结构模态固有频率附近。此外,弯曲与扭转响应存在耦合效应,且耦合作用的频率与提取的垂尾表面气动载荷脉动频率一致。垂尾的位移响应由一阶弯曲模态主导,振幅不大;加速度响应主要由扭转模态产生,量级较大,使结构持续遭受严重的附加惯性载荷作用。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

无尾布局嵌入式舵面的大迎角纵向操纵能力研究

 针对某前掠翼翼身融合无尾布局由鸭面与尾舵组成的纵向基本控制舵面在大迎角状态操纵效率降低的问题,采用数值模拟方法研究一种机身下表面嵌入式新概念纵向操纵舵面实施大迎角纵向操纵补充的可行性。提出了嵌入式舵面的设计思想,研究了嵌入式舵面高度、偏度及其与尾舵组合时的相对位置等参数影响,提出了嵌入式舵面的设计原则、流动机理以及提供低头力矩增量的作用原理。研究结果表明:嵌入式舵面是无尾布局飞机大迎角纵向操纵的高效补充措施,单独使用,最大可提供约平衡10°迎角的低头操纵力矩,并对升阻特性影响很小;与尾舵组合使用,在研究迎角范围内(迎角α≤32°),可提供约6°迎角的低头平衡力矩增量,且对升阻性能产生有利影响。本文工作可为其他翼身融合无尾布局的气动舵面设计提供借鉴。     

基于共轭方程法的跨声速机翼阻力优化设计

机翼的气动特性从根本上决定了飞机巡航状态下的气动特性,经过良好设计的机翼会明显地提升全机的气动特性。连续共轭方程方法由于其求解梯度的计算量与设计变量数目无关,可以对复杂构型进行有效的优化设计。本文采用连续共轭方程方法对粘性条件下的机翼气动外形进行了优化设计。采用有限体积法进行流动控制方程和共轭方程的数值空间离散,多步Runge-Kutta方法进行时间推进,对于某三维机翼进行了有效的气动优化设计。     

运动突风作用下机翼-机身-尾翼亚音速非定常气动力数值计算

 本文研究了亚音速机翼-机身-尾翼复杂平面形状在运动突风作用下非定常气动特性。采用有限基本解法,在不同入射方向和倾斜角下,对不同复杂平面形状的非定常气动力进行了计算,计算结果与实验数据也较接近。     

三维不可压缩层流、转捩、湍流边界层的数值计算

本文提出了一种计算三维不可压缩边界层的数值解法,其中包括对层流边界层的微分法,对湍流边界层的积分法以及利用Orr-Sommerfeld方程计算转捩的数值解,当扰动振幅达到e~(13.5)时发生转捩。 按照上述方法,本文计算了ONERA M6后掠角为30°的机翼,计算结果分别和Cebeci和Cousteix的结果作了对比,二者吻合较好。     

基于动态重叠网格方法的尾翼对螺旋桨滑流的影响

螺旋桨飞机产生的滑流会对其扫掠过的部件产生显著干扰，研究尾翼部件对滑流的影响有助于将滑流与尾翼的干扰进行解耦分解。采用动态重叠网格方法模拟螺旋桨定轴转动，通过求解三维非定常雷诺平均Navier-Stokes （URANS）方程，数值模拟了某螺旋桨飞机带尾翼构型的有/无滑流状态，通过试验结果对计算方法的正确性进行了验证。在此基础上，分别开展了有/无尾翼构型的滑流计算，结果表明：扣除尾翼气动力后，有/无尾翼的升阻力变化规律基本一致，俯仰力矩由于机身后体修形不同呈线性平移关系；对比有/无尾翼空间切面的速度分布云图、不同空间位置和拉力系数下的下洗角和侧洗角变化曲线，发现尾翼对滑流的影响仅局限在其周围，不同拉力系数下尾翼的干扰规律也基本类似。通过研究认为，在飞机初期设计和选型阶段，螺旋桨滑流与尾翼的相互干扰，可简化为滑流单向对尾翼产生影响，尾翼对滑流的影响可以忽略。     

微型扑翼飞行器扑翼/尾翼气动干扰的数值研究

以二维刚性约束条件下的微型扑翼飞行器模型为研究对象,在动网格技术基础上,应用非定常数值分析手段对比分析了单翼/纵列翼布局的气动性能,深入研究了纵列翼缩减频率、扑翼—尾翼无量纲水平间距、来流攻角对其气动性能的影响.结果表明:①纵列翼尾翼对扑翼产生正效应干扰,相对于单翼布局,扑翼—尾翼无量纲水平间距为0.5倍翼型弦长时的纵列翼布局的推力系数和推进效率分别增加28.7%和5.7%;②缩减频率是影响推力的关键参数,随着缩减频率的增加,脱落涡的强度增加,推力系数增大.对于单翼、纵列翼两种布局模式,当缩减频率在1.0附近时推进效率达到最优;③对于纵列翼布局,在扑翼—尾翼无量纲水平间距为1.1倍翼型弦长时推进效率达到峰值;④在0°~20°来流攻角变化范围内,随着来流攻角的增加,升力系数增加,推力系数减小,当来流攻角大于9°时,两种布局的推力均为负值.      

翼身融合无人机外形优化设计研究

翼身融合布局无人机具有较好的升阻特性和隐身特性,拥有广泛的应用前景。为了提高无人机的续航能力,兼顾翼身融合无人机的气动特性和结构重量要求,选择无人机升阻比与全机面积作为优化目标,应用多目标优化方法研究翼身融合无人机的外形设计,提出一种针对翼身融合无人机的外形参数化优化设计方法,并进行实例验证。结果表明:外形优化可以提高无人机升阻比、减轻结构重量,从而获得合理的翼身融合无人机设计方案。