High-order implicit discontinuous Galerkin schemes for unsteady compressible Navier–Stokes equations

Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin(DG) discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta(IRK) scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.     

大型锻件窗框加工变形及工艺方案研究

窗框零件作为民用飞机上的重要承力构件,具有结构复杂、曲率大、难加工部位多、加工变形控制难等特点,而且采用锻件加工受毛坯初始残余应力影响较大。通过对典型窗框零件工艺特点以及加工关键技术进行研究,从毛坯残余应力分布、加工变形仿真预测、加工工艺过程优化等多方面进行研究,解决该类零件加工过程中变形处理的瓶颈问题,提出工程化解决方案,并形成窗框类零件加工典型工艺方案。     

MUSCL格式在跨声速转子流场计算中应用

以工程应用为背景,采用Spalart-Allmaras湍流模型和LU-SGS隐式解法,开发了基于MUSCL格式的CFD程序,并将程序应用到压气机跨声速转子内部流场的数值模拟中,研究了5种差分格式限制器的计算效果.与Rotor 67和Rotor 37转子试验数据的对比表明:当程序采用VanLeer限制器,且入口湍流粘性取值约为3.6倍的分子粘性时,收敛效果最好,能较准确地模拟出流场的总性能参数和激波结构.     

Godunov型显式大时间步长格式研究进展

综述了Godunov型显式大时间步长格式的研究进展。首先介绍了显式大时间步长格式的概念、分类和优势。然后重点阐述了Godunov型显式大时间步长格式的构造方法、高阶精度推广方法、多维问题推广方法和收敛特性、分辨率及计算效率等性能，展示了其在典型问题中的应用和验证。最后给出了Godunov型显式大时间步长格式研究进一步可能的发展方向。     

边界层内WENO格式特征重构精度的理论分析

WENO格式因其高精度以及良好的激波捕捉能力而被广泛地用于可压缩湍流的数值模拟过程中。为了在有强间断的流场中仍然能保证数值稳定，通常将流动物理量投影到特征空间后再进行WENO重构。然而这种特征重构方法在边界层内的精度问题仍需要进一步分析与研究。通过可压缩边界层自相似解构造平板边界层内法向的流动参数，并在该方向上进行WENO重构，得到求解半点数值通量的WENO非线性权重。通过分析实际得到的权重与理想权重的偏差大小发现在边界层流动中特征重构的精度明显低于分量重构。为了更深入地分析特征重构精度降低的原因，通过对曲线坐标系下的流动方程进行理论分析，推导得到特征重构时的特征变量形式。根据特征变量表达式以及自相似解中权重误差分布发现引起特征重构精度降低的主要原因是特征投影过程在半点处出现了额外的极值点，同时这种因特征投影过程而产生极值的特性将不会随着左右特征矩阵的不同选取而发生改变。进一步地，基于理论分析的普适性，对于任何光滑流场，只要使用特征重构的投影过程，都会在半点处产生极值从而使得WENO系列格式产生精度降低的现象。     

利用极值分布概型预测海峡地区风环境期望风速

根据《公路桥梁抗风设计规范》,对于特别重要的桥梁要求对桥址处的设计基本风速进行重新观测或参考实测风速样本进行推算得到,选取合理的基本设计风速.本文选取厦门一漳州跨海大桥附近的厦门市专业气象台和龙海市气象局360个月的极值10分钟平均时距的风速和风向标准值,采用3种极值分布概型推算出不同重现期内的设计风速.得出了厦门海峡桥址处100年一遇的期望风速为35.8 m/s,小于规范给出的厦门地区100年一遇基本设计风速39.7m/s.     

二维进气道不启动流场非定常特性的混合LES/RANS模拟

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)方程模拟方法结合5阶WENO(weighted essentially non-oscillatory)格式对马赫数为3的来流中、内收缩比为1.5的不启动状态下的二维进气道进行了计算,再现了不启动进气道中的非定常流场.计算结果表明:所采用的模拟方法对入口处的平均绝热壁温、摩擦速度和雷诺应力的计算精度较好,进气道不启动流场中激波波系和分离区存在大时空尺度的低频运动,其占主导的特征频率和典型的激波/湍流边界层干扰问题中激波和分离区的低频频率接近,且进气道出现了间歇性的启动状态.      

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

空间机械臂建模方法与控制策略研究

针对刚柔耦合空间机械臂动力学建模中对柔性体采用的传统描述方法（有限元法、模态综合法以及集中参数法等）并不足以精确描述柔性大变形的问题,采用绝对节点坐标法描述柔性体,采用自然坐标法描述刚性体,建立了末端带集中质量的双连杆柔性机械臂的动力学模型并且研究了机械臂的空间定位问题。结合广义α法以及工程上常用的Scaling技术,开发了计算程序,实现了动力学方程的高效精确数值求解。针对机械臂的空间定位以及柔性变形问题,提出了一种运动规划方案,采用PD控制策略,实现了机械臂的运动跟踪控制;仿真结果表明：提出的运动规划方案能有效地减弱机械臂的柔性变形。     

DEVELOPMENT OF HYBRID RANS / LES CODE AND ITS VALIDATION

A computational code based on the hybrid RANSLES approach is developed. The hybrid approach com- bines the delayed detached-eddy simulation（DDES） with an improved RANS-LES hybrid model aiming at wall modeling in LES（WMLES）. In the code, the convective flux is solved using the fourth-order skew-symmetric scheme so as to diminish the negative effect of numerical dissipation. The Spalart-Allmaras（S-A） model is applied as a suhgrid scale（SGS） model. To validate the developed code, homogeneous isotropic turbulence and turbulent channel flow are simulated and the results are compared with experimental data and DNS results. The results of the isotropic turbulence show that the fourth-order skew-symmetric scheme is adequate enough and the model works well coupling with the convective scheme. The results of the turbulent channel flow agree well with the DNS data, the predicted velocity profiles at Reynolds number from 178 to 2 700 match well with the Reichardt＇s law, and the organized vortical structures are well captured.