虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。     

跨声速开槽壁超临界翼型洞壁干扰的N—S方程模拟

本文采用二维可压缩非定常Ｎ－Ｓ方程模拟了带有上下开槽壁的跨声速二维管风洞中超临界翼型绕流的流场。网格生成中内层采用解双曲型偏微分方程，外层采用代数生成法。开槽壁采用只有压力梯度没有顶板运动的Ｃｏｕｅｔｔｅ流动的解析解近似。推导出槽间流动与开闭比的三次方关系。上下壁边界层可由无滑移条件直接模拟，也可用平板边界层近似。本文对驻室压力作了几种近似处理，结果表明，所给Ｎ－Ｓ方程的边界条件与风洞实验多接近一     

基于超椭圆方程和序列响应面法的回转壳开孔形状优化

含有多种功能开孔的轻质回转壳结构是航空航天结构系统内的常用支撑构件,其开孔形状直接影响结构的静动态性能。以回转壳结构为对象,基于超椭圆方程和坐标映射变换推导了回转壳开孔边界的参数化表达,开展了开孔形状动力学优化研究。为提高结构优化计算的精度、效率和收敛性,提出了准等弧长方法和基于均匀设计的序列响应面近似建模方法(SRSM),以分别实现空间超椭圆曲线的精确逼近、减少结构有限元分析成本和加快迭代收敛。以非支配排序遗传算法Ⅱ(NSGA-Ⅱ)作为响应面模型求解算法,结合有限元分析构建了回转壳开孔形状优化设计流程,开展了最大化结构一二阶频率带隙的典型回转壳结构开孔形状优化设计。结果表明,基于超椭圆方程和序列响应面法的开孔优化方法获得了有效改进结构动态特性的回转壳开孔形状,对开展计算耗时工程结构形状优化设计具有一定应用价值。     

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

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

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

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

3种铝合金材料动态性能及其温度相关性对比研究

 对3种铝合金2024-T351、7050-T7451和LY12-CZ进行了温度在77～573 K的静、动态压缩（应变率10^-3～6 000/s）和拉伸（应变率10^-3～3 000/s）试验,分别得到了3种铝合金材料的应力 应变关系和失效应变。对比研究结果表明,随温度升高3种材料的塑性流动应力降低,应变率敏感性增加。最后基于Johnson-Cook模型,拟合了用以预测铝合金材料塑性流动应力的模型参数,其预测结果与试验结果吻合较好。     

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

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

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

错频转子系统动力特性研究

根据某级压气机叶／盘转子的实际结构，建立了叶／盘带短轴的有限元模型。采用NASTRAN结构分析软件，得到了模化协调转子所需的数据，构建了多自由度叶／盘转子系统动力学模型，并在此基础上通过改变多自由度模型中叶片刚度实现错频。采用基尔法求解错频转子系统的运动方程。获得了此多自由度系统的动态响应，并对不同错频方案进行了对比分析。     

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。