跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

突破有限元分析的极限

有限元法FEM分析变得日益复杂,同时有限元分析模型的大小和细节设计要求也在不断增加。由数百万个单元和数百万的自由度组成的有限元网格的模型已经变得司空见惯,然而模型的尺寸也在     

基于支持向量回归机和参数优化的空中交通管制员需求量预测

空中交通管制员作为落实空管安全工作的主体,是整个管制系统的重要环节.本文根据历年民航局相关数据,首先采用基于网格搜索和交叉验证的支持向量回归机(Support vector regression,SVR)预测算法,构建了飞行保障架次与管制员需求量之间的映射模型;然后通过模型预测了 7个地区空管局的管制员需求量;最后根据...     

基于光学几何测量的整体叶盘几何失谐建模与辨识

通过光学几何测量技术获取精确的叶片型面差异化信息（即几何失谐）建立整体叶盘的高保真动力学模型的方法，并进一步开展整体叶盘几何失谐辨识的研究。采用先进的三维结构蓝光扫描系统测量构建精确的叶片几何型面点云模型，然后采用网格变形技术，将谐调叶片有限元模型的表面节点自动投射至实测的点云表面，以回避传统逆向工程的实体模型重建环节，从而实现整体叶盘高保真动力学模型的快速构建。该模型可直接用于量化识别叶片几何失谐对其固有频率和振型的影响，其中各叶片“一弯”频率失谐量在2.1%以内，同时可以精确比对各叶片间的模态置信因子，因此可大幅提高整体叶盘建模和动力学仿真分析的准确性。     

非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用

常规工程应用中，非定常数值模拟（如多体分离）的计算量十分巨大，如果为了达到更高的计算精度，加密网格或者采用高精度方法将会使得计算量进一步增大，导致非定常数值模拟在CFD工程应用中成为十分耗时和昂贵的工作，因此，提高非定常数值模拟的可扩展性和计算效率十分必要。为充分发挥既有分布内存又有共享内存的多核处理器的性能和效率优势，对作者团队开发的非结构网格二阶精度有限体积CFD软件（HyperFLOW）进行了混合并行改造，在计算节点间采用MPI消息传递机制，在节点内采用OpenMP共享内存的MPI+OpenMP混合并行策略。首先分别实现了两种粒度（粗粒度和细粒度）的混合并行，并基于国产in-house集群采用CRM标模（约4 000万网格单元）定常湍流算例对两种混合并行模式进行了测试和比较。结果表明，粗粒度在进程数和分区数较少的小规模并行时具有效率优势，16线程时效率较高；而细粒度混合并行在大规模并行计算时具有优势，8线程时效率较高。其次，验证了混合并行在非定常计算情况下的可扩展性，采用机翼外挂物投放标模算例，分别生成3.6亿和28.8亿非结构重叠网格，采用对等的（P2P）网格读入模式和优化的重叠网格隐式装配策略，网格读入和重叠网格装配耗时仅需数十秒；采用3.6亿网格，完成了非定常状态效率测试及非定常分离过程的湍流流场计算，在in-house集群上12 288核并行效率达到90%（以768核为基准），在天河2号上12 288核并行效率达到70%（以384核为基准），数值模拟结果与试验结果符合良好。最后，在in-house集群上采用28.8亿非结构重叠网格进行了4.9万核的并行效率测试，结果显示，4.9万核并行效率达到55.3%（以4 096核为基准）。     

非结构网格体心梯度求解方法的精度分析

对于空间二阶精度的非结构网格求解器,无黏通量和黏性通量的计算都需要流场变量在单元体心处的梯度值。体心梯度求解的精度很大程度上决定了求解器整体的精度。本文推导了几种常用的梯度求解方法(包括格林高斯法,节点格林高斯法和最小二乘法)求解体心梯度的理论精度,并从数值方面对上述理论精度进行了验证。在数值测试过程中,通过以当前单元体心为基准进行坐标局部缩放的做法,保证了在非结构网格上做精度测试时,网格拓扑能够严格保持不变。理论推导和数值测试一致发现:使用保线性权插值到节点的节点格林高斯法、使用共面单元的最小二乘法和使用共点单元的最小二乘法,不管网格拓扑关系如何,都能保证梯度求解为一阶精度。而格林-高斯法和使用距离倒数权插值到节点的节点格林-高斯法在一般的网格中只有零阶精度,只有在特定的网格和权函数关系下,才能有一阶精度。最后用NACA 0012翼型亚声速绕流模拟进一步验证了梯度求解方法的收敛精度。     

机翼－机身－吊架－短舱的跨声速流计算

本文通过使用多块网格的嵌合体技术，计算了机翼－机身－吊架－短舱的跨声速流场。将机翼－机身－吊架－短舱流场分为机翼－机身和吊架－短舱两个子区域，根据嵌套网格结构建立机翼－机身－吊架－短舱的组合网格，分别对这两个子区域根据相应的边界条件独立求解，然后采用嵌合体技术传递两个子区域的干扰信息，通过两个流场的耦合迭代得到机翼－机身－吊架－短舱流场的计算结果，数值模拟结果与实验结果进行了比较，两者吻合较好。     

应重视对量仪测力的检定

接触式量仪都具有一定测力，测力太小会影响示值的稳定性，测力太大又将产生压变形，尤其是对于软材质薄的金属件及点接触形式的测量，变形更为明显。例如：用杠杆千分尺以三针法测量螺纹中径，便是如此。在相同压力下，三针与螺旋面的接触应力显然比三针与杠杆子分尺测量面间接触应力大，因而其接触变形也就大。压陷量的存在，必然使三针在牙廓间的位置下降，使实测M值小于未考虑测力时的M值，由此引起螺纹中径的测量误差。最近，有个单位检验部门，在用三针法比较测量M6螺纹塞规中径时，发现2把杠杆千分尺检定结果相差3μm，一把尺检定…     

连续波雷达抗网格式干扰

连续波(相控阵)雷达是一种正日益引起人们重视的抗干扰能力优良的雷达,而网格式干扰又是一种对目标指示雷达威胁相当大的转发式干扰。对连续波相控阵雷达抗网格式干扰的能力进行了分析,并提出了对其威胁严重的其它转发式干扰样式和可能采取的抗干扰技术措施。