求解大型非对称特征问题的块Arnoldi方法

本文提出了计算大型非对称矩阵若干个模最大或模最小特征值以及相应特征向量的块Arnoldi方法。研究了块Arnoldi方法的收敛率,推广了Saad关于Arnoldi方法收敛率的一些结果,给出了块Arnoldi方法收敛率的一些估计。提出并讨论了由该方法所产生的数值结果。     

一维ADI-FDTD方法的数值色散分析

给出了交替变向隐式时域有限差分法(ADI-FDTD)的一维差分格式，利用增长矩阵推导了数值色散关系式。将其与过去研究中给出的色散关系式进行了比较。利用牛顿迭代法对其进行了全面的分析，发现了时间步长的选取与网格、色散误差的大小有一定的联系。     

非线性方程迭代算法收敛性分析

从几何上分析了迭代格式xn+ 1=φ(xn)所产生的序列收敛于方程根的收敛条件 ,导出了不依赖函数可导性判断迭代序列收敛的收敛定理 ,给出了产生收敛的迭代序列的技巧。     

二维不可压分离流计算的约束迭代法

本文提出约束迭代法以改进不可压分离流计算中的迭代收敛性。在每个时间步的迭代中,对涡量传输方程的边界涡量,预先进行近似估计。因此它相当于一种隐式的边值算法。 两个翼型的算例表明,约束迭代法具有高的收敛速度。由于无需寻找最优松弛参数,因而节省人力,便于使用。     

A Linear-correction Least-squares Approach for Geolocation Using FDOA Measurements Only

  A linear-correction least-squares(LCLS) estimation procedure is proposed for geolocation using frequency difference of arrival (FDOA) measurements only. We first analyze the measurements of FDOA, and further derive the Cram閞-Rao lower bound (CRLB) of geolocation using FDOA measurements. For the localization model is a nonlinear least squares(LS) estimator with a nonlinear constrained, a linearizing method is used to convert the model to a linear least squares estimator with a nonlinear constrained. The Gauss-Newton iteration method is developed to conquer the source localization problem. From the analysis of solving Lagrange multiplier, the algorithm is a generalization of linear-correction least squares estimation procedure under the condition of geolocation using FDOA measurements only. The algorithm is compared with common least squares estimation. Comparisons of their estimation accuracy and the CRLB are made, and the proposed method attains the CRLB. Simulation results are included to corroborate the theoretical development.     

含误差的齿面接触分析初始点确定算法

研究含误差的螺旋锥齿轮齿面接触分析中初值的选择与确定问题,基于几何等量关系,计算出离初始啮合点P较近的一点P',计算P'的过程不涉及两个齿面在安装坐标系下的旋转等操作,相比直接计算P所求解的方程简单,然后用迭代法计算出初始啮合点P.计算实例表明新方法对于迭代初值的敏感性下降,对迭代初值要求降低,提高了误差齿面接触分析(ETCA)的计算稳定性与效率,该方法使得ETCA的算法更加的实用、计算中初值的选取更加简单.      

空域数据重排的后多普勒自适应处理方法

传统的后多普勒自适应处理方法,如因子法(FA)和扩展因子法(EFA)虽然能大大降低自适应处理时的运算量和独立同分布样本的需求量,但由于实际中均匀训练样本数目的限制,当天线阵元数进一步增大时,FA和EFA抑制杂波和检测动目标的能力会显著恶化。针对这一问题,提出了一种空域数据重排的后多普勒自适应处理方法。该方法将多普勒滤波后的空域数据重排为一行列数相近的矩阵,空域滤波器权系数也表示成可分离的形式,从而得到一双二次代价函数,利用循环迭代的思想求解权系数。实验表明该方法具有快速收敛,所需训练样本少的优点,尤其在大阵列、小样本条件下该方法抑制杂波的性能明显优于FA和EFA。     

弹塑性接触问题边界元缩聚法分析

给出了利用隐含对称条件的二维弹塑性边界元分析公式,并应用于弹塑性接触问题分析。根据接触条件列出了弹塑性接触问题迭代求解的耦合方程。利用接触问题局部非线性的特点,将缩聚法用于迭代求解。基于拟线性、双重迭代法的思想,解决了两重非线性相耦合的问题。最后给出了一个算例。     

基于雅可比矩阵精确计算的GMRES隐式方法在间断Galerkin有限元中的应用

为改善高阶间断Galerkin有限元方法(DG)时间推进效率,在三维非结构网格下针对该方法建立了并行广义最小残差(Generalized Minimal Residual,GMRES)隐式时间迭代方法,GMRES方法基于科学计算工具包PETSc中的Krylov子空间求解器实现。为进一步提高GMRES的计算效率,发展了方程组右端项残值雅可比精确计算方法,针对无黏通量Roe格式和黏性通量BR2(Bassi Rebay 2)黏性计算方法,分别解析给出其对守恒变量多项式自由度的雅可比矩阵。基于建立的方法首先采用NACA0012翼型研究了GMRES的重启次数及收敛参数对方法收敛性影响,然后采用无黏及黏性算例对比研究了基于雅可比矩阵不同计算方法的GMRES计算效率,同时对比研究了雅可比矩阵完全近似求解下GMRES和LU-SGS(Lower Upper-Symmetric Gauss-Seidel)的计算效率。结果表明,建立的基于右端项残值雅可比矩阵精确求解的GMRES方法能够大幅提高不同精度DG方法的CFL(CourantFriedrichs-Lewy)数,相比前面提到的其它方法具有更高的计算效率,其收敛速度实现量级以上的提高。     

空心阴极的发射热特性研究

空心阴极的发射体在电子发射时会处于动态热平衡状态,这种发射热特性对发射体的利用率和工作寿命有重要意义。为研究空心阴极发射热特性及相关热设计方法,建立一种流-热耦合数值模型,将等离子体流场计算与组件热分布计算进行耦合,并采取一种逆向迭代方法来收敛计算,数值模型在验证试验下的计算误差低于11.3%。在此基础上,对发射体温差的机理以及阴极顶孔径对发射体温差的影响规律进行数值分析,主要结论为:羽状模式下的核心电离区较发散,导致发射体温差较点状模式高出50 K左右;适当增加阴极顶孔径可对发射体温差有降低作用。