弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

跨声速翼身组合体欧拉方程计算

本文利用Ｊａｍｅｓｏｎ等提出的有限体积法，迭代求解三维定常欧拉方程，探讨了加快收敛的途径，并用这种计算了大展弦比翼身组合体跨声速定常无粘绕流，计算结果与有关实验及文献进行了比较，获得了满意结论。     

挠性多体航天器姿态动力学建模与分析

针对带有多个挠性附件以及天线转动机构的卫星等多体航天器．利用伪坐标Lagrange方程建立了系统的姿态动力学模型。通过仿真计算．分析了星体与天线运动的耦合特性．以及挠性附件干扰对卫星姿态的影响．验证了数学模型的正确性。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

非结构网格上N—S方程解法研究及在多段翼型绕流中的应用

在非结构网格上求解N－S方程时，靠近物面的网格伸展比相当大，用有限体积法进行显式求解时不易收敛并且雷诺数有一定的限制。本文采用的人工耗散项中含有网络的方向尺度，因而在不同区域选取不同：在粘性区内该项比物理粘性项小得多，在粘性区外与Jameson的人工耗散项相同。生成非结构网格采用的是层推进技术与传统的阵面推进法相结合的方法，这样能保证在粘中性区内网格有一定的方向规则，而在粘性区外网格完全非结构化。     

同心圆柱和同心球壳中格林函数

用分离变量、傅里叶变换和留数定理导出了同轴线内外导体之间的格林函数；用分离变量和正交函数傅里叶级数展开导出二个不同半径的金属同心球壳中格林函数，获得的格林函数解析表达式能够方便地用于求解同轴线中和同心球壳中具有分布源的非齐次波动方程。     

非线性规划的人工释能及单变量曲径寻优算法

给出了求解一般不等式约束的非线性规划问题的一个常微分程的解法。其一维搜索的路径是约束央面上的一条最短线,其方程是由变分法建立的一组常微分方程的初值问题所确定的,。在初始点位于可行域内部时,采用人工释能法来求得下一个改进的可行点。数值例子表明核算法具有较了的计策效果。     

翼型周线布线性涡的有限基本解法

在翼型周界线上布线性变强度分布涡,在一定边界条件下建立求解涡强度的线性方程组。经过试探,确定了比较合适的控制点位置。为满足库塔条件而令后缘处涡强度为零,但导致了方程组的矛盾性。采用最小二乘法求解矛盾方程组,可获得误差最小的唯一解。本方法简便,易于应用。计算结果与试验值或精确解相比,有良好的准确性,可用于一般工程计算。     

任意形状壳体非线性非K－L普适理论

从理性力学的观点导出任意形状壳体非线性非Ｋ－Ｌ理论的普适方程，由此即得Ｋ－Ｌ理论的普适方程，以及各种近似方程，并可将现有文献上的方程作为其近似特例而包含在其中．     

用泰勒级数求解非线性代数和微分方程组

提出一种求解非线性代数方程和非线性常微分方程的新方法。所论方程组被转换成增量形式，而解被表示成泰勒级数。在代数方程（多项式形式）的情况下，此算法归结为求解一系列线性方程组，而且在每一增量步中系数矩阵是不变的；在常微分方程组的情况下，此算法归结为一组递归算式，并不要求解方程组。此方法固有的自动走步功能，可保证得到收敛的解，从而大大节省计算时间。