计算亚音速偏转操纵面机翼升力分布的ALGOL-60源程序

本文提供了计算亚音速偏转操纵面机翼的定常和非定常升力分布的算法语言程序。用本程序计算的结果与其他方法[2]、[3]相比较有相同的精确度。     

带有传动机构的翼段颤振半主动抑制

尝试了将磁流变阻尼器安装于机翼操纵面传动机构上的布局，并由此建立了非定常气动力下机翼-操纵面-传动机构的气动弹性动力学方程，利用on—off控制算法对系统进行颤振抑制，并且研究了控制延迟时间、控制电压以及阻尼器滞回宽度等参数对控制效果的影响。计算机仿真结果表明，控制延迟时间和控制电压对颤振临界速度有较明显的影响，而阻尼器的滞回宽度的影响不大，在本文的阻尼器安装布局和控制策略下，系统的颤振临界速度至少提高了17％。     

类Starship飞行器大迎角动态特性数值研究

再入飞行器的动态稳定性是其任务成败的关键因素,因此有必要对新型类Starship飞行器大迎角动态特性进行深入研究。为辨识动导数,采用基于分区弹簧近似法动网格技术的非定常N-S方程求解器开展强迫振动的数值仿真。有限体积离散采用Roe格式及SST湍流模型,时间方向采用全隐式GMRES方法进行迭代求解。在高超声速弹道外形（Hyper ballistic shape,HBS）标模验证的基础上,获得了类Starship飞行器典型大迎角状态下不同飞行马赫数、重心位置、操纵面偏转角、迎角和减缩频率的俯仰动导数变化规律：典型状态下,俯仰动导数超声速时为负,亚声速时为正;随不同因素的变化,Ma=0.3和Ma=5.0时的俯仰动态特性差异明显;高速下操纵面上偏会导致阻尼减小,且后翼影响更显著。计算结果表明了本文方法在宽速域、大空域复杂外形飞行器动态特性辨识中的应用价值。     

处理动边界问题的无网格/直角网格混合算法

把无网格/直角网格混合算法发展用于涉及动边界的非定常流动计算问题.采用无网格/直角网格混合算法,计算区域整体采用直角网格,只在物面附近嵌入局部无网格.对于动边界问题,本文让主体直角网格保持固定,用原混合算法中的局部无网格区跟随运动边界,使得动边界物面附近有不变的点云结构,避免了点云重构;无网格区与直角网格交界面随边界运动跟踪确定,由于主体基于直角网格,方便了交界面的跟踪.用发展的混合算法,结合非定常Euler方程双时间推进方法,先对翼型作俯仰振荡的非定常跨声速流动进行了数值模拟,计算结果与现有文献和实验数据进行了比较与验证;结合刚体运动方程,成功模拟了二维外挂物投放过程,算例展示出发展的混合算法适合处理大位移运动边界问题.     

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。     

直升机旋翼桨叶动态气动载荷计算方法

为研究桨叶上的气动力，用动态入流模型计算的诱导速度，挥舞变形运动带来的相对气流，以及由于桨叶扭转和操纵线系变形带来的桨距角变化综合计入翼型气动环境，然后用Ｌｅｉｓｈｍａｎ和Ｂｅｄｄｏｅｓ发展的非定常气动模型计算了翼型的气动力。同时，考虑了桨尖形状对气动力的影响。最后，采用状态空间法对方程进行了离散化处理，以适合于计算机计算。编制了相应的计算程序，用于计算桨叶的气动载荷及其变形，并用算例分析了本方法的适用性。     

动边界问题无网格算法及其动态点云

提出基于Delaunay图映射的点云运动处理技术.该技术使点云随物体运动再生成问题变成简单的映射关系,无需迭代,且能保持点云结构不变,故有效地避免了因点云结构变动而产生的附加计算量.运用该技术,结合非定常Euler方程双时间步无网格算法,成功地模拟了翼型的跨声速非定常流动.其中物理时间迭代采用二阶隐式格式,伪时间遮代采用四步龙格一库塔显式格式.文中给出了不同翼型的计算结果,并与实验结果进行了比较,取得了令人满意的结果.     

舰面流场对直升机着舰时悬停操纵的影响

用N-S方程对某型军舰舰面流场流态进行了初步的数值模拟,采用实验数据修正计算结果,并将计算得到的流场叠加到某无人直升机旋翼流场中,初步研究了由于机库等钝体存在而引起的“陡壁”效应对直升机舰面起降时操纵量的影响。最后的计算结果表明,飞行甲板上方由于机库存在引起的垂向气流及其分布对舰载直升机悬停操纵特性影响较大,而侧向分量影响较小。文中着重建立适于工程应用的甲板流场计算简化方法,以用来方便地确定甲板流场和无人直升机着舰时的操纵需求,也可用于舰载直升机起降包络线计算来指导舰载直升机飞行训练     

舰面流场对直升机上舰时悬停操纵的影响

用Ｎ－Ｓ方程对某型军舰舰面流场流态进行了初步的数值模拟，采用实验数据修正计算结果，并将计算得到的流场叠加到某无人直升机旋翼流场中，初步研究了由于机库等钝体存在而引起的“陡壁”效应对直升机舰面起降时操纵量影响。最后的计算结果表明，飞行甲板上方由于机库存在引起的垂向气流及其分布对舰载直升机悬停操纵特性影响较大，而侧向分量影响较小。文中着重建立适于工程应用的甲板流场计算简化方法，以用来方便地确定甲板流场     

前飞状态直升机旋翼/机身非定常气动干扰的分析

一个计算旋翼/机身/尾迹间非定常气动干扰的分析方法是建立在二阶升力线/全展自由涡模型和机身面元模型基础之上的。通过迭代机身在桨盘平面、尾迹定位点的诱导速度和旋翼/尾迹在机身表面的诱导速度，形成一个耦合的分析模型。在分析中计入了非定常项。作为算例，对两种孤立机身表面的平均压强系数分布以及旋翼机身组合体中机身表面的非定常压强系数分布进行了计算，其结果与实验值相吻合。     

应用Plemelj积分生成二维通道内的正交网格

首先选用的解析函数由速度模以及流向角组成，在已知通道边界上的流向角时，求得通道边界上的速度。由势函数和流函数的定义，求出壁面上的势函数及流函数。其次将解析函数用势函数表示，运用Ｐｌｅｍｅｌｊ边界积分关系的内点公式，求得通道内任意点的势函数和流函数。由此，从等势线和等流函数值线，可形成流体力学数值计算用的正交网格。文中给出两个算例，该方法无需迭代。     

目标电磁散射特性的混合法研究

从Ｓｔｒａｔｔｏｎ－Ｃｈｕ公式出发，利用阻抗边界条件，得出阻抗边界条件积分方程，将目标表面分成不同的区域，有些区域用高频解析方法求出其表面电流，有些区域利用矩量法解阻抗边界条件积分方程得表面电流。本文以橄榄体（包括介质涂敷情形）为例，计算其电磁散射特性，混合法计算结果与实测结果或其他方法的计算结果比较吻合。     

电磁力连续控制圆柱绕流态变化的研究

作用于电介质溶液流体边界层上的Lorentz电磁力可以改变流体边界层的结构,从而控制流体绕流的形态。通过理论分析和数值模拟确定了实验控制的关键参数,交替分布的电极和磁极包覆在圆柱体的表面置于电介质溶液中,简单调整电磁力的分布可以方便地控制圆柱尾迹的流动结构,实现了电磁力消涡和增涡的连续控制。     

关于风洞中用尖劈和粗糙元模拟大气边界层的讨论

尖劈和粗糙元广泛地应用于风洞试验中的大气边界层模拟,该技术成功模拟了不同地貌特征的平均风速和湍流度剖面。随着风工程研究的深入,了解尖劈和粗糙元模拟过程中的作用机理有助于准确地模拟各种大气边界层湍流功率谱和尺度特性。试验表明:尖劈利用其迎风平板的分离流产生湍流涡旋,迎风板的宽度决定了涡旋的大小和湍流脉动强度,同时迎风板阻塞比沿高度递减产生近似线性的风速剖面;粗糙元用于模拟实际地面的摩擦效应,调整平均风速和湍流度的剖面分布。遗憾的是,尖劈下宽上窄的结构特点决定了该技术模拟的湍流功率谱和积分尺度的高度变化律与实际大气边界层相反。基于对模拟机理的认识,异型尖劈上部形状有助于模拟大比例模型试验要求的湍流风场。     

不可压缩Navier－Stokes方程的基本解法

利用不同阶次的Ｈｅｌｍｈｏｔｚ算子方程的基本解，建立了时间离散不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的完全边界积分表示式。对俯仰谐振的ＮＡＣＡ００１２翼型所计算得的升力和阻力的滞后现象，与实验结果相一致     

多层介质涂覆腔体电磁散射特性的IBC/FEM分析

将Leontovich阻抗边界条件(IBC)与矢量有限元法(EB-FEM)相结合对内壁涂覆多层介质开口腔体的电磁散射特性进行分析。将IBC和腔体口面边界积分方程分别转化为第三类边界条件形式，从而导出边值问题的相应泛函。表面输入阻抗通过广义反射系数得到，计入了介质内部电磁波的多次反射，提高了IBC／FEM混合法的精度。数值结果表明这种IBC／FEM方法在分析涂覆多层介质开口腔体电磁散射问题中的便捷性与精确性。     

Potential Method in the Theory of Thermoelasticity with Microtemperatures for Microstretch Solids

The linear equilibrium theory of thermoelasticity with microtemperatures for microstretch solids is con- sidered. The basic internal and external boundary value problems （BVPs） are formulated and uniqueness theorems are given. The single-layer and double-layer thermoelastic potentials are constructed and their basic properties are established. The integral representation of general solutions is obtained. The existence of regular solutions of the BVPs is proved by means of the potential method （boundary integral method） and the theory of singular integral e- quations.     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

基于PIV技术的压力场重构算法实现与研究

介绍了有限容积法、直接积分法和Poisson方程法3种基于PIV瞬时速度场重构压力场的基本原理以及相应的计算方法，选取管流突扩流场和偏置方块绕流流场两个不可压缩流场的瞬时速度场数据，采用上述3种压力场重构算法，分别研究了图像噪声、速度场精度、插值算法以及边界条件的类型与精度对重构压力场的影响。最后针对管流突扩过程第20ms的流场，给出了3种重构算法下的压力场云图以及对应的CFD模拟结果。研究表明，有限容积法和直接积分法容易受到噪声的影响而产生剧烈震荡，但是可以在较大的速度场误差范围内保持较高的精度，通过采用双线性插值可以获得更高精度的重构压力场；Poisson方程法不易受到噪声的影响而产生震荡，同时在高精度PIV速度场下的优势较为突出，通过采用双三次差值可以获得更高精度的重构压力场；混合边界条件仅仅测定边界上有限个点的压力值，就获得了接近狄利克雷边界条件下重构压力场的精度，远高于诺依曼边界条件；边界条件的误差严重降低重构压力场的精度，其影响程度比速度场误差还要大。     

DIRECT NUMERICAL SIMULATIONS OF MULTIPHASE FLOWS

INTRODUCTIONMultiphase and multifluid flows are commonin many natural and technologically importantprocesses. Rain,spray combustion,spray paint-ing,and boiling heat transfer are just a few ex-amples.While it is the overall,integral charac-teristics of such flows that are of most interest,the global behavior is determined to a large de-gree by the evolution of the smallestscales in theflow.The combustion of sprays,for example,depends on the size and the number density ofthe drops.Generally…