全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

航空交叉流式湿热交换器的热性能分析

 本文分析研究两种流体互不混合的湿空气交叉流式热交换器的热性能。按热质平衡原理,首先建立了一套描述其热湿性能的微分方程组,而后应用有限差分技术把它们化为差分格式进行计算,最后与实验结果进行了比较,获得较为满意的结果。     

三维控制翼周围粘性分离流动的数值模拟

本文首先分析了中心型隐式格式和迎风型隐式格式的优缺点,在此基础上发展了一种自由参数可确定的反扩散隐式格式。用此格式计算了绕三维控制翼的层流和湍流分离流动。出发方程是三维简化N-S方程,用Baldwin-Lomax湍流模型。结果表明,这种隐式格式效率高,能满意地刻划流场。     

直升机非线性运动方程及其数值分析

基于牛顿法建立了单旋翼带尾桨直升机的非线性运动方程，提出了一种求解直升机非线性运动方程的数值方法。建立方程时，用欧拉角描述直升机在空中的姿态，旋翼采用有挥舞铰且在铰上带有弹性约束的模型，桨盘处的入流采用线性模型，通过对一算例直升机的操纵响应进行数值仿真，将非线性与小扰动线性化的直升机运动方程进行了对比分析。结果表明：从操纵输入到随后的2s左右，非线性和小扰动线性化的直升机运动方程的计算结果非常接近，但在此之后，非线性运动方程的计算结果更能反映直升机的运动规律。     

求解常微分方程组的小波技术

研究发展一种新的求解常微分方程的数值方法—快速小波配置法 (FWCM) ,该方法与传统的时间推进或频率区间方法完全不同。快速小波配置法是将任意函数展开为小波基函数 ,用快速离散小波转换技术 (DWT) ,有效地构造常微分方程的近似解。计算过程中 ,在小波展开层次和自变量区间两个不同方面采用了多层自适应和多区间自适应处理技术 ,提高了计算过程的稳定性和收敛性 ,并且具有均匀的误差分布。特别是在常微分方程的长时间解方面 ,与Runge-Kutta方法比较 ,具有稳定的长时间性态。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。     

均匀来流中旋转振荡圆柱绕流的数值研究

通过求解二维非定常不可压缩 Ｎ-Ｓ方程,研究了均匀来流中旋转振荡圆柱的流动特性,系统地改变流动的振荡频率、振幅,分析了圆柱的受力特性及尾迹结构,发现了与横向、流向振荡的圆柱相似,旋转振荡也是流动控制的一个有效方法,有其自身的流动规律及特点。     

喷管分离流动的数值模拟

针对高空喷管地面试验时产生流动分离的现象,用Ｂｅａｍ －Ｗａｒｍ ｉｎｇ 近似因式分解法求解薄层Ｎ－Ｓ方程,对分离流场进行数值模拟。预示了分离点的位置以及喷管壁面压强分布等,可以为高空发动机地面试车提供参考, 并为进一步研究分离的影响因素打下基础     

绕大攻角翼型低速非定常流动的SUPG有限元数值模拟

采用既有良好稳定性又有较高精度的ＳＵＰＧ有限元方法,通过数值求解涡量一流函数形式的不可压缩流ＮＳ方程组,模拟大攻角翼型在高Ｒｅ数下几种典型情况的非定常绕流流场,同时验证了ＳＵＰＧ有限元方法在求解对流占优的二维非定常对流的－扩散方程时的有铲性。