基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

内埋式弹舱舱门气动载荷计算分析研究

采用中心有限体积法进行空间离散和四步Runge-Kutta显式推进方法进行时间推进数值求解了非定常欧拉方程;采用Delaunay方法生成内埋式弹舱舱盖打开时的非结构网格.数值模拟了内埋式轰炸机弹舱单开舱门在分别打开8°和45°时的复杂流场.通过对上述复杂流场的分析,发现作用于舱门的气动载荷具有使其趋向于关闭的特性,且有随打开角度的增大而减小的趋势.     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

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

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

基于迹向控制量的编队星群相位控制

为了得到控制量与控制时机对相位改变的直观公式,基于Hill方程以及迹向控制的约束条件,定性分析共面绕飞编队中控制量、控制时机与编队卫星相位之间的关系,分析结果表明:当需要改变的相位角小于90°时,可实现最省燃料的控制;当需要改变的相位角大于90°时,可进行趋于最小控制量的控制;当在相对运动椭圆上下点实施控制时,相位无变化.     

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

Automatic differentiation based discrete adjoint method for aerodynamic design optimization on unstructured meshes

A systematic methodology for formulating, implementing, solving and verifying discrete adjoint of the compressible Reynolds-averaged Navier-Stokes (RANS) equations for aerodynamic design optimization on unstructured meshes is proposed. First, a general adjoint formulation is con-structed for the entire optimization problem, including parameterization, mesh deformation, flow solution and computation of the objective function, which is followed by detailed formulations of matrix-vector products arising in the adjoint model. According to this formulation, procedural components of implementing the required matrix-vector products are generated by means of auto-matic differentiation (AD) in a structured and modular manner. Furthermore, a duality-preserving iterative algorithm is employed to solve flow adjoint equations arising in the adjoint model, ensur-ing identical convergence rates for the tangent and the adjoint models. A three-step strategy is adopted to verify the adjoint computation. The proposed method has several remarkable features:the use of AD techniques avoids tedious and error-prone manual derivation and programming;duality is strictly preserved so that consistent and highly accurate discrete sensitivities can be obtained; and comparable efficiency to hand-coded implementation can be achieved. Upon the cur-rent discrete adjoint method, a gradient-based optimization framework has been developed and applied to a drag reduction problem.     

高压输电网设计中的二阶非线性方程的次谐波分枝

本文用Melnikov方法研究高压输电网设计中提出的非线性二阶受迫振动方程的次谐波分枝,分析了系统可能出现次谐解的条件和参数区域,以此为工程应用提供参考。     

旋转对管内流体流动的影响

本文描述了旋转管流的基本运动方程。在笛卡尔坐标系和柱坐标系中,针对两种基本的几何流动,推导出旋转情况下管内的流动方程,定义了影响旋转管流的流动和换热的两个无量纲准则数J和Ra,分析了二次流产生的原因,从而说明旋转对流动的影响以及旋转管流的运动特征。