Comsol_Multiphysics_求解器简介sjtu.pdfVIP

MEMS 交流网 Comsol Multiphysics求解器 ——摘记自《COMSOL Multiphysics基本操作指南和常见问题解答》 Comsol Multiphysics 提供了很多不同的求解器供用户选择，根据求解器类型 可分为如下几种 1、稳态Stationary：求解稳态问题，即偏微分方程只含有空间上的偏微分。 2 、瞬态Time-dependent ：求解瞬态问题，即偏微分方程中还含有时间上的 偏微分。 3、特征值Eigenvalue ：求解特征值问题 4 、参数化Parametric ：按照一系列参数求解稳态问题。 5、稳态分离式Stationary segregated：采用分离耦合方式按照一系列参数求 解稳态多物理场问题。 7、瞬态分离式Time dependent segregated ：采用分离式耦合方式求解瞬态多 物理场问题。 从求解原理来看，可以氛围直接式线性系统求解器和迭代式线性系统求解 器。典型的直接求解器包括：UMFPACK 、SPOOLS、PARDISO 以及 Cholesky （TAUCS ）等，迭代式求解器包括：GMRES、FGMRES 、共轭梯度法CG、BICGStab 以及集合多重网格法等，其中迭代式求解器还需要一些预条件求解器、后平滑器 和粗化求解器等的支持，例如，Vanka 、SOR、SSOR、SORU 一直前面提到过的 各种求解器均可以作为这些预处理器使用。 一、求解器简介 Comsol Multiphysics 是采用有限元方法求解偏微分方程的大型仿真软件，有 限元发发是通过离散的手段，将秋节余离散成一系列相互关联的小单元，最终形 成一个离散的线性代数方程组，然后通过求解器解这个方程组，得到单元节点上 的解。至于节点外的任意位置，则通过插值的方法来得到结果，即： (cu) f Ku F 式中：c——偏微分方程的系数项 f——源项，转换成线性方程组后： K——刚度系数矩阵； F——载荷向量 u——解向量，其熟练也称为自由度（DOF ） 对于这种方程的求解，最理想的状态当然就是直接通过求逆矩阵的方程得 到向量解：u K 1F 然而，这在真实情况下是可遇不可求的，因此也产生大量的求解方法或求 解器。例如上面提到的在Comsol Multiphysics 中使用的直接式求解器和迭代式求 解器。 另外，当系数项中出现了因变量的函数时：[c(u )u ] f Ku F 上面这个偏微分方程具有较强的非线性，增加了求解的难度，对于这类非 线性问题，一种有效的求解办法是通过迭代的方法，采用线性求解器反复迭代， u 对于一个近似线 最终得到正确的结果，即，首先以初始值 代入方程的系数项， 0 MEMS 交流网 性方程求解u K (u )1 F 再以得到的结果代入系数项，进行下一步近似线性方程 1 0 的求解u K (u )1 F,un K (un )1 F 直到满足精确要求 ，即un un1  ，其中 2 1 1 n 为迭代次数。 二、直接式求解器 直接式求解器通过一步“求逆”方法进行求解，其中最常用的手段是进行 LU 分解，将刚度矩阵分解为上三角（U ）和下三角（L ）矩阵，即K=LU ，然后 1 1 分别对这两个矩阵求逆，然后求解，即u U L F 。 目前，Comsol Multiphysics 中提供的直接求