用保角转绘坐标解二维和轴对称可压流Euler方程

本文推导出用保角转绘正交曲线坐标表示的二维与轴对称Euler方程,方程中的未知量以曲线坐标方向的速度分量表示,并有不含未知量导数的源项。方程组的形式与直角坐标系中的接近,因而简化了计算。本文将该方程组用于计算二维任意翼型和旋成体轴对称绕流,采用分辨率较高的Eberle特征平均显式有限体积差分算法求解,在不同区域中用不同的显式推进步数加速收敛。本方法已分别对翼型和旋成体轴对称绕流编制了计算程序,可用于自由流为亚声速、跨声速和超声速的绕流计算。大量计算结果表明,本方法精度好,激波分辨率高,编程序简单。     

不可压Navier—Stokes方程组的SUPG有限元数值解

本文从定常不可压NaVier-Stokes方程组出发,构造了SUPG加权剩余公式。为保证数值解的精度,本文对速度取八节点插值,保留了摄动项中的二阶导数项。从用本文方法所做的算例来看,计算结果是令人满意的。     

热弹性复合材料层合板的响应分析

利用混合变分原理,考虑了热平衡方程与导热方程中变量的对偶关系,通过增加非齐次正则方程的维数,成功地将热弹性体机一热耦合问题的非齐次Hamilton正则方程转化为能独立求解的齐次Hamilton正则方程．所以大大简化了问题的求解。并通过数值实例研究了在温度栽荷和动力栽荷作用下及考虑阻尼的四边简支复合材料层合板的响应问题。     

Legendre-Gauss拟谱法求解最优控制问题

 提出一种新的求解基于常微分方程(ODE)和微分代数方程(DAE)的最优控制问题的数值方法。本方法基于直接配置法,通过Legendre-Gauss拟谱法同时离散化状态变量和控制变量把最优控制问题转化为一个非线性规划问题。与传统的直接转换法相比,本方法具有精度高、计算量小、结构简单的特点,而且可以求解最优控制“多相”问题。数值结果表明,本方法是一种通用的精度较高的最优控制直接数值求解法,可用于求解ODE/DAE最优控制问题。     

一种GPS整周模糊度的解算方法

初始整周模糊度的求解是利用GPS载波相位进行测量的关键问题,文中提出了一种GPS整周模糊度快速解算方法。该方法：首先,对双差观测方程中指向卫星的方向阵进行QR分解;然后,通过比较估计出的基线长度和真实的基线长度建立起整周模糊度初始搜索范围;最后,综合利用快速剔除不合理模糊度解的检验方法以得到最终的模糊度解。文中对所提出的算法进行了理论推导,并进行了基于实测数据的实验与分析。实验结果表明,该算法计算量小,快速准确,适合整周模糊度的快速动态求解。     

航天器用低温多层隔热性能计算方程

航天器常用有效发射率表征多层的隔热性能,其经验值范围为0.02～0.04,经验值仅与多层单元数有关,不随温度变化.有效发射率经验值的适用条件是多层的热面温度约-10℃～50℃,多层的冷面不照太阳.使用条件偏离得越多,有效发射率的经验值导致的计算偏差就越大.针对该问题,提出采用辐射项加导热项的混合传热模型,取代传统的纯辐...     

基于简单WENO-间断Galerkin的Euler方程自适应计算

为了得到Euler方程的高精度、高分辨率数值解,介绍了间断Galerkin方法、三角形单元上简单WENO限制器的基本原理以及基于自适应网格加密的激波捕捉方法。将简单WENO限制器-间断Galerkin方法应用到曲边四边形单元上,通过单元边界上高斯积分点的坐标来搜索相邻单元从而得到相邻单元的单元编号,实现了基于“问题单元”的局部网格加密自适应计算。对若干典型问题进行编程计算,结果表明,简单WENO限制器可以应用到曲边四边形单元上,且可适用于局部网格加密时具有“悬挂节点”的非结构网格上的激波捕捉。      

高超声速飞行器刚体／弹性体耦合动力学建模

高超声速飞行器广泛采用升力体、乘波体等气动布局和轻质复合材料、薄壁结构等,导致结构振动与刚体运动频率非常接近,给飞行器制导控制系统设计带来了巨大挑战.针对该类飞行器的特点,考虑结构的横向位移,将机身前后体简化为于质心处固联的2根悬臂梁,并从统一的能量观点出发,基于拉格朗日方程与虚功原理,在纵向平面推导出适合高超声速飞行器的刚体/弹性体耦合动力学模型.通过对比耦合模型与传统刚体模型的极点分布情况,发现结构振动与刚体短周期模态紧密耦合,离心力的引入影响了高度与长周期模态,对高超声速飞行器航迹运动的作用不可忽视.最后分析了飞行速度与结构阻尼变化对耦合模型动态性能的影响.结果证明飞行速度对刚体运动模态影响显著,而结构阻尼的变化主要改变弹性模态.     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

高精度的对流差分格式

对于广义对流差分格式,本文应用两种方法(截断误差分析法及事后估计误差法)导出了使其误差最小的最优α值。和以前的方法不同,它们是从所研究的对流-扩散方程本身出发导得的,而无需事先知道方程的解析解。大量算例结果表明,截断误差分析法有很高的精度。事后估计误差法的误差,除个别点外,也都比一般格式的误差要小或者小得多。但是,对于复杂的问题,前者的应用可能会有困难,而后者却很易实现。