人造卫星运动方程数值解中的几个主要问题

本文详细阐述了人造卫星运动方程数值解中的数值稳定化和步长均匀化问题,给出利用能量关系或采用辛算法控制沿迹误差的快速增长和自变量变换解决大偏心率轨道积分过程中的变步长问题,并用计算实例验证了方法的有效性。     

导弹运动方程的逆解析解

针对在采用反馈线性化理论设计导弹飞行控制系统时控制器算法对导弹逆模型的要求 ,推导了导弹运动方程的逆解析解 ,包括导弹质心运动的运动学方程及其动力学方程的逆解析解、姿态运动的运动学方程及其动力学方程的逆解析解 ,详细给出了逆模型的求解过程 ,为一类飞行控制系统运用全局线性化设计奠定了基础。最后 ,结合工程应用进行了仿真计算 ,仿真结果表明逆解析解的推导是正确的     

TVD格式计算一维Euler方程间断解存在的几个理论问题

本文讨论了TVD格式按本文第二节A段意义推广到拟线性方程组的几个理论问题。指出:1.对于Riemann间断分解,Euler方程的解和Riemann不变量都没有TVD性质,因此不能将TVD格式自然推广;2.TVD格式是提高分辨率的一种途径,它不能保证得到物理解,能得到物理解的TVD格式一定是耗散守恒格式,因此将TVD格式推广时,必须保证格式是耗散的,否则仍会得到非物理解。并给出了一阶迎风格式,TVD2,UNO2得到的若干非物理解算例;3.非物理解现象只是发生在具有剧烈膨胀的单波区的计算上,而某些TVD格式若会得到非物理解,则一定发生在特征值变号的单波区的计算上。     

二维横向喷流的N-S方程数值解

让一股射流横向注入超音速流产生所需要的干扰效应是一种飞行器飞行控制的常用手段,如航天飞机和机动导弹上的反作用控制系统(RCS).所以,对射流与主流相互干扰的某些流动特性应该有个基本的了解.本文通过应用MacCormack显格式和Baldwin-Lomax修正的代数湍流模型求解二维RANS方程对带横向射流的绕后台阶的M_∞为2.19和6.00的超音速外干扰流场进行了数值计算.与无喷射流的情形相比,回流涡旋区由2个增加到4个,且底部压力有大的回升,喷口宽度越大,回升效应越显著.另外,数值计算表明,在压力梯度变化比较平缓的区域附加一四阶人工粘性能有效地抑制数值振荡,加快收敛.     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

人卫数值法定轨中运动方程右函数间断的处理

人卫精密定轨中受摄星历（或称精密星历,即状态转移）,可由分析解或数值解相应的定轨方法亦有分析法定轨与数值法定轨之称。对于后者,在一般情况下,现有的常微分方程数值法（或称积分器）已能满足精度要求,除长弧定轨中遇到的困难外,还有一类问题值得注意,即地影“间断”和轨道机动飞行导致的间断。对这种间断的处理关系到如何在保证星历精度的前提下提高计算效率的问题。本文针对这一问题,给出了相应的处理措施和计算方法,     

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

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

二维分离流中的NS方程数值计算方法

本文采用显式多步龙格－库塔方法求解了二维可压非常完全ＮＳ方程，为了使计算适用于一般的风格，采用有限体积法进行几何离散，研究定常问题计算中的一些加速收敛措施，如：当地时间步长、残值光顺技术等。讨论了模拟分离流问题中的技术问题，计算结果表明本文的计算结果与国外文献是吻合的。     

管道中超声速旋涡运动的数值模拟

本文利用ＥＮＮ格式通过求解ＮＳ方程，数值模拟了超声速旋涡在管道中沿其轴向的演变规律。结果表明，当流动处于顺压区，在涡轴附近，截面流线为不稳定的螺旋点形态；当流动处于逆压区，在涡轴附近，截面流线为稳定的螺旋点形态；当从顺压区过渡到逆压区，涡心附近出现不稳定极限环；反之，出现稳定极限环。数值模拟结果同张涵信的拓扑分析理论十分相符。     

绕流问题的区域收缩法及其数值解

本文提出了一种适于绕流问题的数值方法,它考虑了流场中不同区域内不同的物理现象,给出了外部边界条件的良好逼近。使用这种方法可以节省大量的计算机存贮单元和计算时间,也可以推广到三维绕流问题的数值计算。     

固液相变问题的边界元法数值解

相变问题的数值计算方法有多种,本文用边界元法解应用工程中常见的相变导热问题。     

跨音速小扰动方程的多重解

本文利用Engquist-Osher格式,离散跨音速小扰动方程(TSD)和边界条件,选择不同初场计算NACA0012翼型二维定常位势流场,得到了与试验值较符合的结果,同时用Murman-Cole非守恒和守恒两格式计算并进行了比较。着重用数值实验进一步探讨多重解问题,发现并总结出了几点很有参考价值的规律和现象。     

KdV方程的数值解及数值模拟

文章在已构造的一种分片光滑的Lagrange型插值多项式空间中利用Galerkin法对非线性KdV方程进行求解。分析了数值格式的线性稳定性并进行了数值模拟,结果表明格式精确有效。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

解不可压N—S方程的PCIMG法

本文提出了一种解不可压N-S方程的预测、修正、迭代、多重网格法。这种算法程序简单,结果精确及收敛快速。     

飞机穿越下冲暴流安全起飞运动方程的研究

推导了存在风切变时在气流坐标系（相对风轴系）中飞机质心空间三维运动方程和在体轴系中飞机的旋转运动方程，指出了按航迹坐标系，（绝对风轴系）建立移动运动方程时某些不便之外。为了进行比较，给出了分别应用上述两组运动方程计算的波音－７２７飞机起飞行轨迹的不同结果。