非定常N－S方程的线化及其差分求解

在微幅简谐振荡的假设下，对计算平面内非定常Ｎ－Ｓ方程进行了线比，并在ＭａｃＣｏｒｍａｒｋ于１９８０年提出的半隐格式的基础上，发展了一种适用于求解计算平面内线化Ｎ－Ｓ方程的差分格式．用该差分格式求解计算平面内的线化Ｎ－Ｓ方程时，可以避免求解五对角方程；且时间步长可以取得较大，计算效率高．     

超音通流风扇叶栅流场数值模拟

用矢通量分裂法求解Ｎ－Ｓ方程和Ｅｕｌｅｒ方程，对超音通流风扇叶栅流场进行了数值模拟，由计算结果给出的物理图象表明，叶栅通道内没有强激波，亦未出现分离回流，全叶栅流场均为超音速，说明计算的两种叶型符合设计要求，其性能有待实验检验。     

空气龄的数值解法

提出了用ＣＦＤ法求解空气龄的具体方法，其步骤为先求解Ｎ－Ｓ方程，后求解空气龄。为求解Ｎ－Ｓ方程，采用了改进型ＨＳＭＡＣ法，并给出了压力和速度同时修正的方法。为求解空气龄，推出了空气龄输送方程式和给出了空气龄和速度成分的边界条件。计算机模拟验证了求解方法的有效性并由空气龄求出了通风效率。     

高空大气成份的数值分析

本文给出一种用有限差分法求解高空大气成分数密度的方法。基本方程是二阶非线性偏微分方程组,采用隐式差分格式,为了节省存贮量和避免求解大型非线性方程组的复杂性,对于每一积分时间步,采用分组线化Gauss—Seidel迭代法,以其迭代过程的收敛值来近似非线性和隐式差分式。在本文给定的边界条件下,将边界条件的差分式联立于内点的差分方程组中,构造了一种三对角线矩阵的追赶算法。以本算法编制的计算程序对于五种大气成分,在几种不同状态条件下进行试算,和现有的观测资料比较,结果是令人满意的。     

Euler方程数值模拟绕悬停旋翼桨叶的流动

本研究应用格心格式有限体积龙格－库塔时间推进法求解三维欧拉方程,模拟绕旋翼桨叶的流动,民流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角,选择Ｏ－Ｈ型式生成绕桨叶贴体的弦向和展向网格,通过算为进一步寻求有效的Ｅｕｌｅｒ或ＮＳ方程数值解积累了经验。     

共轭梯度－快速傅里叶变换解齿状结构散射场

利用共轭梯度－快速傅里叶变换法（ＣＧ－ＦＦＴ）求解了谐振区齿状散射目标的雷达散射截面（ＲＣＳ）．通过共轭梯度－快速傅里叶变换求解电场积分方程（Ｅ－ＦＩＥ），给出齿状散射体的雷达散射截面．适当选取基函数和检验函数，将电场积分方程离散化，从而获致准确度较高的ＣＧ－ＦＦＴ实用程序，节省内存，提高了计算速度，缩短了计算时间．     

共轭梯度—快速傅时叶变换解齿状结构散射场

利用共轭梯度－快速傅里叶变换法（ＣＧ－ＦＦＴ）求解了谐振区齿状散射目标的雷达散射截面（ＲＣＳ）。通过共轭梯度－快速傅里叶变换求解电场积分方程（Ｅ－ＦＩＥ），给出齿状散射体的雷达散射截面，适当选取基函数的检验函数，将电场积分方程离散化，从而获致准确度较高的ＣＧ－ＦＦＴ实用程序，节省内存，提高了计算速度，缩短了计算时间。     

气动热CFD计算的格式效应研究

通过对N-S方程的求解计算了二维圆柱绕流及三维绕钝双锥高超声速流动过程中的气动热问题.对于求解方程,在空间离散的格式上,选用了3种上风格式——Roe的FDS,Van Leer的FVS和AUSM+格式以及一种中心差分格式,将各种格式所得的热流结果与实验结果进行了比较.对于AUSM+格式,还应用了minmod limiter 和 双 minmod limiter 两种限制器,并对结果进行了对比.      

Lax-Friderichs格式在磁流体模拟中的改进和应用

磁流体数值模拟是空间物理研究的重要手段.采用具有TVD（Total Variation Diminishing）特性的Lax-Friderichs差分格式求解了GLM-MHD（Generalized Lagrange Multiplier-Magnetohydrodynamics）方程.为降低格式的数值耗散,引入耗散修正系数对算法的通量计算过程进行改进.二维Rotor算例和磁云-电流片相互作用算例的模拟结果表明,GLM-MHD方法可以有效控制磁场散度误差,相对于泊松校正法可以节省一半以上的计算时间.在不破坏格式稳定性基础上,耗散修正系数有效降低了算法的数值耗散.      

吴方法在平面并联机构位置正解中的应用

采用吴方法对平面并联机构位置正解问题进行了研究.吴方法是一种求解非线性方程组的数学机械化方法,采用这种方法任何非线性方程组都可以在有限步内得到解决.在给出了吴方法基本原理的基础上,对本问题进行了求解,并将原始方程组转化成为一个三角化的方程组.其中单变量方程的次数为6次,说明平面并联机构可以有6个不同的位姿.最后用数值实例进行了验证,给出了计算结果.吴方法在这一问题中的应用,为求解其它机构学难题提供了新途径.     

求解Navier—Stokes方程的指数型差分法

我们以Burger's方程为例,推导了指数型差分格式,并证明了这格式是无条件稳定和高精度的。此外,我们利用这格式来求解一维单个和多个Navier-Stokes方程,取得了令人满意的结果。把这格式进一步推广,得到求解二维Navier-Stokes方程组的格式,并从理论上证明了这一推广格式的一系列优点,还以此来求解实践问题,得到了与实验和理论非常符合的结果。     

刚性方程的一个算法

在文献[3]中构造了一个与二阶Gear方法等价的差分格式,并采用简单迭代求解。本文详细论述了差分格式在实际计算中应采取的若干措施及程序框图,最后给出了数值例子。     

差分格式的优化组合

用偏微分方程组解实际的三维紊流流动问题时,如何选择各方程中对流项和扩散项的差分格式将对计算的稳定性、准确性与经济性影响很大.对一个实际的流动问题,计算机模拟验证了迎风差分、中心差分、混和格式、PLDS(Power Law Differencing Scheme)、QUICK(二次迎风插值)和OPTIMAL(Oscillation Preventing Technical algorithm based on the second-order Interpolation Method for Advection Linked convection)的各种组合对计算的稳定性、准确性与经济性的影响.对中心差分应用在雷诺时均方程中出现的数值解的振荡问题提出了在局部网格使用迎风差分的具体对策.对QUICK应用在κ,ε输运方程上出现的数值解的发散问题提出了QUICK和迎风差分格式交替使用的具体方法.      

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

极大似然法测定p—S—N曲线的置信度

对极大似然法测定的Ｐ－Ｓ－Ｎ赋予置信度，提出了极大似然法珠Ｐ－Ｓ－Ｎ曲线的置信限曲线公式，通过该公式确定的Ｐ－Ｓ－Ｎ能使估算的疲劳寿命具有置信度一给出了单点法测定出构件Ｐ－Ｓ－Ｎ的数据自理。     

绕大钝头短体外形跨音速大迎角湍流流动计算

采用数值模拟求解N-S方程方法计算与分析了绕大钝头短体外形的跨音速大迎角时的湍流流动.数值方法中选用Van Leer 分裂格式离散无粘通量项,并构造了一种Limiter函数以保证TVD(Total Variation Diminishing)性质,采用中心差分格式来离散粘性通量项.计算中采用了B-L湍流模型.算例定性地表明绕流的三维分离具有准周期性的非定常特性.     

刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

三维抛物化N-S方程全隐格式矢通量分裂求解

采用全隐格式的矢通量分裂算法,在通用曲线坐标系内求解三维无因次抛物化强守恒型N-S(Navier-Stokes)方程.所建立的计算程序,可计算较为复杂边界的三维流场,如求解亚、跨、超音速共存的流场,分区模拟内外流共存流动及其混合剪切尾流.给出了一轴对称亚跨超音速喷管内流场的计算结果,包括计算网格,流矢分布,压力、温度和马赫数等值线分布.也给出了一喷管内外流及其混合尾流的计算结果,包括计算网格,流矢分布,无因次压力、密度和温度等值线分布.     

Euler方程的分裂型通量分裂双时间步隐式方法

传统隐式方法有格式复杂、计算量大等缺点,在Euler方程的差分离散过程中,利用算子分裂思想,结合通量分裂法、双时间步法等隐式离散方法,构造了一种更简单的分裂型隐式计算方法.通过对典型空气动力学问题的计算,检验了该方法的有效性和可靠性,并对其性能做了具体讨论.该方法具有稳定性好、时间步长约束小等隐式格式的普遍优点,同时具有格式简单、程序易实现等优点;避免了传统隐式方法单步推进时的方程组常规求解及矩阵求逆过程,计算量小;比LU-SGS方法收敛速度快.      

固体火箭发动机燃烧室三维流动数值计算

采用ｋε双方程湍流模型,以ＳＩＭＰＬＥ 计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流ＮＳ 方程．在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构．计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡．