高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

基于N-S方程的翼型气动优化设计

本文简要回顾了气动优化设计的最新发展，并采用连续伴随方法对粘性条件下的翼型气动外形进行了优化设计。设计过程中，采用B样条曲线的控制顶点作为设计变量，通过求解伴随方程得到目标函数关系设计变量的梯度，最后用最速下降法和BFGS变尺度法对二维翼型进行了外形优化，验证了本方法直接、快速的设计能力。     

基于分区的隐式求解二维不可压NS方程的并行实现

本文通过对目前常用并行算法的分析,基于区域分解的思想,给出了适用于隐式求解原始变量二维不可压NS方程压力修正方法的并行算法,并在网络异构编程PVM环境下进行了数值验证,讨论了计算结果,结果显示对任务级的并行,并行效率也同样随粒度的增加而提高.另外此方法还可推广到三维复杂流场的计算问题.     

基于数据链的空战对抗建模及增援决策分析

研究了空战双方力量的变化问题.将飞机的空战能力指数转化为平均战斗力水平,进而以蓝彻斯特方程为基础,考虑数据链对空战的影响和作战中有增援力量的加入,建立了对抗双方战机数量变化的微分方程模型,给出了其数值计算方法,指出了离散计算的时间步长所应满足的条件,并根据增援模式的不同,将模型细化为不同形式.仿真实验表明,存在有效增援时间范围,一次性增援和匀减速增援具有时间优势,并且通过构造评价函数有效衡量了这两者具有较高的增援效率,为空中作战应用和决策提供了参考依据.      

随机粗糙面的双站散射模型分析

以先进积分方程模型(AIEM)为工作模型,从理论上系统分析了随机粗糙面的双站散射特征,探讨了不同粗糙程度下散射系数在方位平面上的谷点位置,并利用数值仿真模拟、实测值对理论预测的散射系数进行了说明与验证.结果显示,散射系数在方位向上的谷点位置与粗糙程度和极化特征密切相关.利用数值仿真模拟,可从物理的角度更系统深入地分析统计粗糙表面的散射特性,从而有利于设计以后的更有效的双站探测模式.      

基于CUDA的超声二维声场EFIT仿真

随着图形处理器（GPU）的快速发展,基于计算设备统一构架（CUDA）可以方便地将并行计算技术应用于超声声场数值仿真计算,极大地提升计算效率。阐述了弹性动力学有限积分算法（EFIT）的原理,在采用CPU实现带吸收边界的钢材料二维点源激励声场仿真的基础上,基于GPU实现了仿真模型的并行计算,介绍了GPU程序的设计流程和参数优化方法,包括纹理内存使用、吸收边界优化和数据传输优化。对比了相同条件下CPU和GPU仿真计算的耗时和平均计算效率,定量分析了GPU对于EFIT模型效率的提升。比对结果表明,EFIT具有良好的并行计算条件,采用并行计算方法能够有效提升模型计算速度,对于复杂声场仿真应用具有广阔的应用前景。      

求解奇异最优问题的控制变量转移法

形如integral from t_0 to t_fx~T(Qxdt)二次型性能指标的最优化问题,是一类较为常见的奇异控制问题。本文就这类奇异控制提出一种新的求解方法,利用正交变换进行控制变量转移,将系统原控制向量置于重新构造的哈密顿函数H之外,选择状态变量中的一部分作为系统新的控制向量,且该新的控制向量与H不再呈线性关系,于是可以运用一般的求解非奇异最优控制问题的极大值原理来求解奇异最优控制和奇异最优轨线。本文提出的求解方法较为简便,所得结果与其他方法所得结果完全一致。     

求解代数Riccati方程的一种新算法

本文讨论了Hamiltonian矩阵在辛相似变换下的标准形,由此给出代数Riccati方程存在非负定解的一个充分条件,提出了求解代数Riccati方程的一种新的算法。该方法节省运算量,且精度较高,尤其适合于解阶数不太高、系数矩阵为满阵的Riccati方程。最后给出了一个数值例子,并将该方法与其他方法作了比较.