非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。     

多重网格在突扩燃烧室流场计算中的应用

多重网格法是一种具有快速收敛性特点的新计算技术,该法在求解偏微分方程时用一系列逐步加密或减疏的网格去离散求解区域,不同粗细网格可以消除不同波长的误差,从而加速了收敛。本文将此法应用于计算突扩燃烧室流场。SIMPLE算法和紊流k-ε模型被采用来预估旋流器后气流速度u,v,w和紊流参数k,ε分布以及回流区尺寸。计算结果表明,在达到同样计算精度条件下,与单层网格相比,此法可以加快收敛速度,缩短机时,减少计算费用,从而使数值计算在工程应用方面得到进一步推广。     

Banach空间中椭圆方程Dirichlet问题的迭代解

以半序理论为工具，在Ｂａｎａｃｈ空间中用上下解方法研究了椭圆方程Ｄｉｒｉｃｈｌｅｔ问题的迭代解，并证明了存在两个迭代序列分别一致收敛于该问题的最小解和最大解。     

轴对称进气道外罩跨音速非守恒型全位势流AF3算法

本文针对轴对称进气道外罩跨音速绕流问题的特点,采用了贴体坐标网格;并对这类网格拓扑构造了三种Baker的隐式近似因式分解AF3有限差分迭代算法。通过理论分析与数值实验,找到了既能计算亚临界情形,又能计算超临界情形的AF3格式;并就此格式讨论了适合于进气道绕流特性的AF3迭代格式中间变量的提法,具体提出了格式实施的相应方法,研制了可供分析任意轴对称进气道外罩跨音速特性的计算程序。由于采用了任意曲线坐标系上的轴对称非守恒型全位势方程适合于用不同方法生成的计算网格,因此,本文的工作还有可能进一步推广到计算更复杂的轴对称进气道跨音速绕流问题。计算表明,本方法收敛快,结果好。     

解不可压N-S方程的预测修正迭代法

本文发展了作者早先提出的一种算法,使计算效率提高了20～108％。在雷诺数Re为100、400,和1000,网格节点为15×15的Cavity流场中进行了数值试验。苦取速度收敛误差为10~(-3),则在IBM PC/XT微机中所需的计算迭代次数和时间分别为42次、77次、86次和10分、21分、23分。     

基于泊松方程的二元翼型网格生成技术

利用求解Posson 方程的方法来划分计算网格，使用Hilgenstock 方法来 控制网格间距和网格正交度。对于安装在风洞中的二元翼型，在进行Navier－Stokes 方程 计算时，需要在物面附近的网格间距很小，这导致网格划分时计算收敛速度很慢。 本文对 Hilgenstock 方法进行了改进，大大提高了网格划分的计算收敛速度。     

基于预处理和网格自适应的风力机翼型气动力计算

用预处理后的Navier-Stokes方程数值模拟低马赫数下风力机翼型的静态和动态绕流流动。为了提高计算的精度,运用了网格自适应技术:初始网格的流场计算收敛后,采用流场速度的紊乱度捕捉湍流区域;采用直接分裂初始父单元的H型网格加密方法对湍流区域进行局部加密。基于初始网格与自适应网格的计算结果与实验结果的对比表明,网格自适应方法提高了计算的精度。     

旋翼翼型低速动态失速研究

为了提高旋翼翼型动态失速模拟的精度,基于动态混合网格技术和ALE形式的RANS控制方程,构建了一套可用于低速流场中旋翼翼型动态失速分析的计算方法。采用守恒变量形式的低速预处理技术,解决了由于特征值差异过大引起的收敛困难问题;在物面采用层推进泊松方程光顺法生成结构网格,以获得较好贴体性和正交性;采用分离流中应用广泛的k-ωSST湍流模型捕捉深失速下流场的大分离特性。计算结果表明该计算方法可以有效地分析不同马赫数下的旋翼翼型动态失速,收敛精度有不小于两个数量级的提升。针对低速流场不同马赫数下深失速的流场特征的计算分析表明,马赫数对动态失速的迟滞特性具有明显的规律性影响。     

求解对称矩阵极端特征的Chebyshev—PBL方法

为了加速预处理块Lanczos方法的收敛法，本文采用组合Chebyshev迭代和预处理块Lanczos方法，提出了求解大型对称稀疏矩阵极端特征的一种新方法－Chebyshev－PBL方法。数值结果表明，新方法对计算大型对称稀疏矩阵的几个最大（或最小）特征值是有效的。     

特征向量导数计算的扰动法

评述了有关特征向量导数的计算方法。综合了它们的优缺点后,提出了两种计算特征向量导数的扰动迭代法。方法的概念简单、实施容易。算例表明,迭代的收敛率较高,适于工程应用。     

关于Newton迭代公式的几个改进

通过对Newton迭代公式进行改进,本文构造了三种新的迭代公式。迭代公式I是一种单步迭代公式,在单根附近具有二阶收敛速度,且无须求函数的导数值;迭代公式II也是一种单步迭代公式,在单根附近具有三阶收敛速度;迭代公式III是一种两步迭代公式,具有至少三阶收敛速度,虽然该公式形式比较复杂,但是具有计算时不需求函数的导数值的优点。此外,证明了三种新的迭代公式的收敛性。最后,通过数值实验验证了三种迭代公式的有效性。     

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

Riccati和Lyapunov矩阵代数方程解的一种新的简便的快速算法

本文讨论了Riccati和Lyapunov矩阵代数方程解的矩阵符号函数算法,证明了牛顿迭代公式对矩阵符号函数的收敛性。同时,指出了Riccati矩阵代数方程解P的相互等效的四个公式,并给出了计算这两种方程解的程序框图和算例。     

双流混合边界的正交贴体网格生成

通过数值实验表明，用泊松方程变换形式生成的贴体坐标网格，其边界点的β值与源项呈连续单调的函数关系。根据这种关系，可以用方程迭代求根方法中的增值寻根法和对分法，使网格边界正交化。由薄壁结构相隔的两种流体混合边界正交网格生成的实例表明，这种方法对于生成适应于红外抑制器全流场计算的网格是很有效的。     

圆柱壳板非线性屈曲的分叉解

采用商用有限元分析软件MSC/Nastran,对两直边铰支、其余两边自由的各向同性圆柱壳板中心受横向集中载荷下的非线性屈曲和后屈曲过程进行了分析,简单地采用了不对称网格技术就可获得该基准问题的分叉屈曲解,计算得到分叉屈曲解的极限载荷比基准解极限载荷小10%左右。对采用了多种不对称网格划分得到的模型进行了分析,以研究不同网格划分对最终结果的影响。结果表明,虽然不对称网格技术十分简单,但是也有一定的局限性,文中给出了能得到分叉屈曲解的网格划分方法,还分析了类似的复合材料圆柱壳板问题。研究结果表明,即使采用对称网格划分也能够得到该问题的分叉屈曲解,并阐述了其原因。     

动边界问题无网格算法及其动态点云

提出基于Delaunay图映射的点云运动处理技术.该技术使点云随物体运动再生成问题变成简单的映射关系,无需迭代,且能保持点云结构不变,故有效地避免了因点云结构变动而产生的附加计算量.运用该技术,结合非定常Euler方程双时间步无网格算法,成功地模拟了翼型的跨声速非定常流动.其中物理时间迭代采用二阶隐式格式,伪时间遮代采用四步龙格一库塔显式格式.文中给出了不同翼型的计算结果,并与实验结果进行了比较,取得了令人满意的结果.     

直升机载火控雷达抗干扰捷变波形优化算法

针对直升机载火控雷达面向抗干扰的捷变波形优化问题,基于奇异值分解(Singular value decomposition,SVD)和循环算法提出了适用于多脉冲和短码长的捷变相位编码波形设计方法。本文首先通过构造具有目标区域低旁瓣模糊函数的带约束四次型优化问题,然后根据SVD分解将四次型转化为二次型,给出了迭代收敛的优化算法;同时在循环算法基础上给了一种四次型的循环迭代收敛算法。最后仿真结果表明,两种算法在设计捷变波形优化性能上接近,但SVD分解算法收敛速度更快,循环计算算法具有更高的运算速度。     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。