轨道转移计算中的逐次消元迭代法

本文针对导弹飞行轨道转移的一个算例,叙述应用消元法思想构造数值求解非线性方程组的迭代算法的一般原理,具体算法步骤和一些计算结果。给出算法成功的理论依据,最后给出实现算法时的一些变形。     

用隐式多重网格法计算三维粘性流动

求解高雷诺湍流流动时,边界层法向网格间距较流向和展向相比非常小,因此边界层中存在高度伸缩的网络,这将大大降低多重网格的求解效率。通过谐调的处理多重网格过程的各个细节,既提高了解的精度,又克服了网格展弦比的影响,使得计算效率提高了6～7倍。     

翼型外形高气动效率/低可探测性的优化

 采用ｖａｎＬｅｅｒ矢通量分裂格式求解Ｅｕｌｅｒ方程的方法计算了绕翼型的气动特性;采用矢通量分裂方法计算了绕翼型的时域电磁散射场特性及雷达散射截面积（ＲＣＳ）;采用一种简单而有效的数值优化方法对流场解和电磁场解进行了翼型外形高气动效率／低可探测性的优化计算。算例结果表明,本方法提供了一种对翼型既可作气动优化设计亦可进行多学科综合设计的有效工具。     

基于非交错有限体积格式和 k-ω 湍流模型模拟叶栅湍流分离流动

发展了一种求解任意坐标系下二维Ｎ－Ｓ方程的非交错有限体积法和一种高阶ＴＶＤ对流格式。采用一种ｋ－ω湍流模型和几种不同的ｋ－ω模型对叶栅湍流分离流动进行了详细的分析。数值结果与得到的实验结果比较表明,此方法不仅精度高,而且具有较好的数值稳定性。     

考虑测量约束的非合作目标近距离快速导引控制方法

研究了考虑测量约束的非合作目标近距离快速导引控制方法。首先,基于C-W方程建立了解析形式的径切联合近距离快速导引控制模型,并采用微分修正方法将解析结果进行修正;然后推导了任意初始条件下,径切联合控制双脉冲与视场角的解析关系;最后基于相对偏心率矢量建立了径切联合控制的控后效果评估模型。仿真结果表明:提出的控制方法和控后效果评估模型有效。     

基于B-P本构模型的涡轮盘应力分析

针对Bodner-Partom统一粘塑性本构方程,采用Euler前差显式积分方法,编制了用户子程序,把该本构方程结合进了大型通用有限元软件ANSYS,进行了一维和二维的算例验证,并针对某涡轮盘进行了三维应力分析。该程序可考虑材料变形的率相关、循环加载等情况,适用于航空发动机高温部件非线性结构分析。      

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

热空腔-喷气流的组合辐射

讨论了热空腔(喷气发动机的加力筒)和喷气流的组合辐射问题, 推演了组合辐射的辐射传输微分方程, 导出了有边界辐射条件下的解。讨论了由多个温度不同的壁面组成且腔内充满高温燃气的热空腔辐射传输方程解的具体形式, 导出了它向后半球的光谱辐射强度的表达式;同时还推导了热空腔-喷气流组合辐射的光谱辐射强度的表达式。最后给出算例, 算例结果表明, 热空腔-喷气流组合辐射的空间分布具有与实测分布较一致的梨形包络, 其光谱分布也与实测的光谱分布的一致性较好。      

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。