求解大型对称特征值问题的块Chebyshev-Lanczos方法

本文提出了计算大型对称矩阵若干个最大或最小特征对的块Chebyshev迭代法,讨论了块Chebyshev迭代法对块Lanczos方法的应用,给出了块Chebyshev-Lanczos方法。计算实践表明块Chebyshev-Laaczos方法比块Lanczos方法和Chebyshev-Lanczos方法都优越。     

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

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

黑体辐射反问题的正则化Lanczos方法

黑体辐射反问题就是利用测量的黑体辐射能量谱确定黑体的区域温度分布.本文将该问题离散化为欠定线性不适定问题,提出求解欠定线性不适定问题的正则化Lanczos方法.该方法基于Lanczos双对角化过程,用一系列小型线性方程组逼近原不适定问题,应用截断SVD正则化方法使迭代稳定化,并用L-曲线确定相应的正则化参数.数值结果说明了新方法的有效性,并说明了在观测数据出现误差的情况下新方法也能有效地重构区域温度分布.     

求解大型对称线性方程组的循环收缩Lanczos算法

向Krylov子空间中加入一些模接近于零的特征值对应的特征向量能够加快收敛速度，事实上，对于这些模接近于零的特征值对应的特征向量，可以用Krylov子空间方法得到，并且在新的Krylov子空间形成的过程中，近似特征向量的近似度会不断提高，特别在标准Krylov子空间方法中，如果因为这些特征向量而减缓了收敛速度，则随着这些特征向量的近似度的提高，用增广Krylov子空间方法解线性方程组的收敛速度会明显加快。Lanczos算法是求解大型对称不定线性方程组的有效方法之一。但在计算过程中由于Lanczos向量失去正交性减慢了收敛速度。本文根据增广Krylov子空间方法提出循环收缩Lanczos算法，新算法充分利用Lanczos过程所得到的谱信息，确定预处理，从而加速Lanczos算法的收敛速度。     

求解大规模矩阵问题的Krylov子空间方法

求解大规模矩阵问题包括线性方程组和特征值问题等是计算数学和科学工程计算中的重大课题，最近几年，其研究工作取得了许多重大进展。文中给出大型线性方程组和特征值问题Krylov子空间方法若干进展的一个概述，其中包括作者对这些问题的研究成果。涉及的专题包括求解大型线性方程组的共轭梯度法、SYMMLQ算法、MINRES算法、GMRES算法、Lanczos双正交化算法、QMR算法以及这些算法的块格式；求解大对称特征值问题的Lanczos算法和块Lqnczos算法；求解大型非对称特征问题的Lanczos算法、Arnodi算法以及这些算法的推广。讨论求解大规模矩阵问题的加速技术和预处理技术。了一些有待进一步研究的问题。     

求解大型非对称线性方程组的UNSYMMLQ方法

本文给出了求解大型非对称线性方程组的Lanczos方法的一个判据,提出了求解非对称方程组Ax=b的UNSYMMLQ方法,它是Paige和Saunders求解对称线性方程组的SYMMLQ方法的推广。文中描述并讨论了一些数值试验。     

一种新的尾迹-地面干扰修正方法

在使用涡尾迹方法对处于地面效应状态下的旋翼进行气动计算时,会出现由于部分涡丝运动至地面下方而导致的尾迹迭代不易收敛、尾迹结构计算不够准确等问题。针对这一问题,基于空气的低速不可压理论,提出了一种新的尾迹-地面干扰修正办法。该方法与之前的方法相比,能将地面附近的涡丝修正到更准确的位置。算例表明,引入该方法的计算模型,获得的尾迹位置总体精度提高,尾迹迭代的收敛速度也有所增加。     

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

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

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

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

加速SIMPLE流场算法收敛的探讨

为加速SIMPLE算法收敛,采用了双线追赶迭代技术加速线性计算部分的收敛,并用外推技术从稀网格解外推出密网格解的逼近值作为后者的初场,减少了迭代次数。这两种方法结合可明显提高计算效率。文中的两个算例的机时节约率分别为48％和24.5％。     

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

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

亚声速飞机外挂物气动载荷和分离轨迹工程估算方法

本文简要介绍了一种计算亚声速飞机外挂物气动载荷和分离轨迹的工程方法。本方法分别采用源汇模型和涡格模型模拟母机的体积效应和法向力效应,采用迭代方法计及母机-外挂物之间的二阶干扰。在气动中心高速所研制的战术导弹气动特性工程计算方法的基础上,以迎角沿弹身轴线和翼片变化的流动条件代替该方法中的均匀来流条件,而建立了非均匀流场中外挂物气动载荷的计算方法。最后采用四阶的 Adams 数值积分方法求解六自由度运动方程而得到外挂物的分离轨迹。与国内外其它计算方法相比,本方法具有适用范围广、迅速、方便、实用等特点。本方法对一系列算例进行了计算,其结果与风洞实验数据具有令人满意的一致性。     

设计灵敏度分析的迭代模态法

在结构动态有限元数学模型修正、结构动力学修改以及结构动态设计中，灵敏度分析是十分重要的一环，往往也是主要计算工作量所在。如何在保证精确度前提下减少灵敏度分析工作量和计算时间，具有重要意义．本文提出的迭代模态法，可以大大提高特征向量对设计参数导数（模态灵敏度）的计算效率。本文还从理论上证明经典模态法和修改模态法只是迭代模态法的特殊情况。文中还给出了ＦＯＲＫ，３Ｄ－ＦＲＡＭＥ两种结构为实例的计算机仿真，并与已有两种模态方法进行比较。理论分析与计算机仿真表明，本方法不仅能充分保证特征向量灵敏度分析的精度，而且大大提高了计算速度．     

一种高度并行的Schwarz型混乱松弛法及其收敛性证明

文[1—2]把混乱松弛思想引入到Schwarz交替法中,构造了一种Schwarz型混乱松弛法。但这个方法在进行第n+1步迭代时,在拟边界上必须要用到第n步迭代的值,从而影响了算法的并行性,得不到相应的同步或异步MIMD并行算法、为此,本文给出一种高度并行的Schwarz型混乱松弛法,这个方法包括了Schwarz交替法及其相应的同步和异步MIMD并行算法。对于二阶线性与非线性微分方程Dirichlet问题,本文应用微分方程极值原理证明了该方法的收敛性。     

梁的三维空间大转动的有效处理方法

有限转动角的计算和迭加方法与线位移不同。为了能利用修正的Lagrangian法及带有动坐标的迭代法来分析三维空间大转动梁的受力和变形，研究了三维空间大转动梁的共轴转动机理，引入了伪矢量的概念．严格地应用数学和弹性力学理论导出了相关的公式，给出了迭代过程。解决了分析空间大转动梁的带有动坐标的迭代法所涉及到的有关有限转动的两个关键问题：(1)单元畸变中节点转角位移的计算；(2)迭代过程中新的一轮转角位移的确定。文中所给出的公式和迭代过程在计算机上易于操作，数值结果也证明了公式的正确性。