不完全LU分解预条件BICGSTAB算法实现感应测井二维FDFD快速正演模拟

在柱坐标系下推导了二维感应测井差分格式,采用频率域有限差分方法求解感应测井正演问题。针对差分近似得到的线性方程组系数矩阵是大型稀疏复系数病态矩阵求解困难等问题,采用不完全LU分解预条件的稳定双共轭梯度(BICGSTAB)算法求解该线性方程组。研究结果表明,本算法具有速度快、精度高和稳定性好等优点,能有效提高感应测井正演模拟的效率和精度。

不完全LU分解预处理的BICGSTAB算法在大地电磁二维正演模拟中的应用

基于双二次插值的有限单元法求解大地电磁二维正演问题,以不均匀网格剖分为基础,推导出大地电磁响应的计算公式。针对有限单元法最后形成一个线性方程组,系数矩阵是大型稀疏的带状对称正定复系数矩阵,并且其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难等问题,采用不完全LU分解(即上三角与下三角分解)处理的稳定双共轭梯度算法(BICGSTAB算法)求解该线性方程组,通过对层状介质和二维模型电磁响应进行计算,获得二维大地电磁的视电阻率曲线和阻抗相位曲线。研究结果表明,BICGSTAB算法具有速度快、精度高和稳定性好等优点。

基于预处理BICGSTAB法的电力系统潮流并行计算方法

为实现大规模电力系统潮流的准确、快速求解,以非精确牛顿法为基础,提出一种基于CPU-GPU异构平台的电力系统潮流并行计算方法。修正方程组的求解是牛拉法潮流计算中最为耗时的部分,提升修正方程组的求解效率可有效提升潮流计算效率。为此,根据雅可比矩阵的不对称不定性,采用稳定双正交共轭梯度(bi-conjugate gradient stabilized,BICGSTAB)法进行修正方程组的求解。进一步,为改善BICGSTAB法的收敛性,根据雅可比矩阵的稀疏性和类对角占优性,提出一种改进PPAT(Preconditioner with sparsity Pattern of AT,PPAT)预处理器和改进Jacobi预处理器相结合的两阶段预处理方法,并对雅可比矩阵进行预处理,提升BICGSTAB法的收敛性能。然后,将上述潮流算法移植到CPU-GPU异构平台,实现电力系统潮流的并行求解。最后,通过不同测试系统算例对所提方法进行验证、分析。结果表明,所提潮流并行计算方法可实现电力系统潮流的准确、快速求解。

预处理的BiCGSTAB算法及其在涡流求解中的应用(英文)

针对涡流场有限元分析形成的大型稀疏复线性方程组,提出了预处理稳定双共轭梯度法(PBiCGSTAB).利用二维正交链表结构实现系数矩阵的全稀疏存储,采用无填充的不完全LU分解对方程组进行预处理以降低系数矩阵的条件数.将稳定双共轭梯度法从实数领域扩展到复数领域,并利用它求解预优过的复线性方程组.基于有限元法和PBiCGSTAB算法编制直线电机性能分析软件,并对一电机堵转的情况进行仿真,将计算结果和ANSYS计算出的结果进行了比较,证明了该软件的正确性.并且通过求解器比较发现,在相同精度条件下PBiCGSTAB算法只需要BiCG算法三分之一的迭代步数,证明了该算法的快速性.

多自由度块行压缩存储技术及大型稀疏方程组的求解

在地球物理电磁勘探领域有限元数值模拟中,最后都会得到一个大型稀疏的复系数线性方程组,受计算机内存空间的限制,必须根据有限元刚度矩阵的稀疏性对其进行压缩存储。由于电磁场有限元计算的自由度大都在三个以上,因而提出了适合多自由度的块按行压缩稀疏存储方案,并通过存储格式的转换,把块按行压缩方式转换成流行的,大型稀疏矩阵的行压缩存储格式,以便于求解。用求解大型稀疏方程组的Krylov子空间方法中的稳定双共轭梯度(Bicg-stab)方法,收敛速度快,精度高,而且稳定性好,结合ilu预处理技术,可以大大提高求解大型稀疏方程组的效率。

多重网格算法在地球物理电磁三维正演中的应用

对使用多重网格方法解决时谐电磁问题进行了研究。在地球物理勘探中针对深部的目标体,只需研究相对低频的电磁波特性,大地中的光速电磁波可忽略不计,主要的研究对象是感应电流的扩散效应。对控制方程使用有限积分法进行离散化,然后使用多重网格法迭代求解。在均匀网格条件下,多重网格法显示了出色的收敛性,在非均匀网格中表现不佳,尽管多收敛速度可以通过使用双共轭梯度稳定法得到补偿,但处理大型非均匀网格的效率仍然不尽如人意。

基于双共轭梯度稳定法的数值保角变换计算法

根据复变函数中的保角变换,研究了基于模拟电荷法的单连通区域的数值保角变换法.针对约束方程的求解问题采用双共轭梯度稳定法,求得模拟电荷,构造了近似保角变换函数.用数值实验说明了本文算法的有效性,并将该保角变换法应用于图像变换中,将图像通过保角变换并显示出变换后的图像.