COMSOL培训复习课程.ppt

求解器参数，高级 瞬态或参数求解中清理内存 单独设定各因变量的精度，提高收敛性，如： u 1e-4 sigma 1e6 集成刚度矩阵过程中的网格单元数量单位 得到收敛和一致的初始条件的技巧 对瞬态问题，设置一个很短时间内的缓冲 使用参数求解器来加强高度非线性问题的收敛性 使用瞬态求解器得到稳态解 如果是多物理场问题使用分离式求解器 使用高度非线性选项 4、COMSOL典型算例分析 多孔介质中的有效扩散率 PID浓度控制 复制网格 通过复制网格可以在不同的边界上产生一致的网格 自动检查目标面上源网格的取向 拉伸和旋转2D网格 在2D中绘制几何并剖分网格 把网格拉伸或旋转形成3D结构 网格导入 NASTRAN 数据文件 最通用的网格文件格式 一阶和两阶单元 从网格开始建模 跳过创造几何结构和网格剖分步骤 网格统计 检查各单元类型的单元数目 检查自由度数目 检查网格质量 可单独检查各子域、边界或边上的统计信息 网格可视化 显示基于坐标的单元（比如，横截面） 使用彩色图案显示网格质量 网格框架 轻松地在不同网格间转换 手动调整网格，应用于求解器设定 自适应网格 稳态计算中，根据结果调整网格 模型库Heat Transfer ModuleProcess and Manufacturingcontinous casting 移动网格案例：电化学抛光 问题描述 电化学抛光：利用金属电化学阳极溶解原理进行修磨抛光 简化的2D模型由两个电极和复合电解液构成。正极有一个凸点，表示表面缺陷。模拟了一段时间后凸点和周围的电极材料的损耗 模型使用了传导介质DC 和瞬态移动网格 (ALE) Model Library Path: COMSOL_Multiphysics/Electromagnetics/electrochemical_polishing 几何模型 电场边界条件 30V v=0 移动网格边界 Vx=0 Vx=0 dx=0, dy=0 Vn=-K*Jn_dc K是比例系数，Jn_dc是法向电流密度 计算结果 自适应网格的设定 在求解器参数设定对话框中设定 选中 线性求解器 直接求解器 UMFPACK, SPOOLES, TAUCS, PARDISO等 易于使用，鲁棒性，占用内存大 适于处理小规模问题，高度非线性和多物理场问题 迭代求解器 GMRES, FGMRES, Conjugate Gradient, BiCGSTAB等 占用内存少，更多的选择，调整比较困难 应用于特定的物理场，如，EM，CFD等 需要预处理器，网格框架，平滑器等 直接线性求解器 直接线性求解器通过一步“求逆”得到结果u=K-1F： 把Ku=F分解成LUu=F，所以L和U是容易求逆并且具有鲁棒性 u=U-1L-1F 等同于Gaussian消去法 优点：鲁棒性强 缺点：内存开销大 直接线性求解器 UMFPACK 对一般的非对称矩阵是鲁棒和高效的 要么计算成功，要么运算内存不足 SPOOLES 利用对称矩阵 内存使用比UMFPACK有效，但计算速度较慢 PARDISO 利用对称矩阵 和UMFPACK类似，但使用内存比SPOOLES少 共享内存式并行处理 在矩阵分解过程中不需要选主元从而节省内存，这导致不精确的因子 由于支持并行的折中处理，不是100%的鲁棒性 TAUCS 非常适合于对称，正定矩阵 迭代求解器 对于规模较大的问题（单元数多，自由度大），直接求解器计算会出现内存不足 矩阵分解是很耗内存的 L和U比K具有更多的非0元素 如何避免分解？ 迭代求解器： 不形成L和U 精细的迭代策略 对每一步测试是否r=Ku-F~0（即是否左侧等于右侧） r 称为残差(residual) 迭代求解器和预处理器 为了在合理的计算时间内达到收敛，迭代求解器需要一个好的初始估算值 利用预处理器 预处理器M是K的近似值，预处理后的系统变为 M-1Ku=M-1F ? Au=B A=M-1K, B=M-1F 预处理后的系统收敛较快，Au=B 比 Ku=F 更容易（快）求解 通常，迭代方法根据前面的残差（r=Ax-b）对u进行较小的改变 迭代线性求解器 GMRES 在前面所有搜索方向上最小化残差，直到重新开始 如何调整重新求解前迭代步数（默认为50） 更节省内存 --- 减小 得到较好的鲁棒性 --- 增加 FGMRES GMRES的一个灵活的变种 能有效地处理更多类的预处理器 比GMRES开销2倍多的内存 Conjugate Gradient 对称正定问题 在计算时比GMRES更快、内存使用效率更高 预处理器 不完全 LU （Incomplete LU） 最具有鲁棒性 内存要求大 代数多网格 （Algebraic Multigr