comsolmultiphysics和数值分析基础.docVIP

第一章 COMSOL MULTIPHYSICS 及数值分析基础 W. B. J. ZIMMERMAN，B. N. HEWAKANDAMBY Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom E-mail: w.zimmerman@shef.ac.uk 本章主要介绍 COMSOL Multiphysics 在零维和一维模型数值分析方面的几 个关键内容。这些内容包括求根、步进式数值积分、常微分方程数值积分和线性 系统分析。这几乎是所有的化工过程数学分析方法。下面通过 COMSOL Multiphysics 中的一些常见化工过程应用实例来介绍这些方法，包括：闪蒸、管 式反应器设计、扩散反应系统和固体中热传导。 1．简介 本 章 内 容 很 多 ， 可 以 分 为 几 个 不 同 的 目 标 。 首 先 介 绍 了 COMSOL Multiphysics 的主要工作特性；其次介绍了如何使用这些特性来分析一些简单的， 位于零维空间、一维空间或“空间－时间”系统中的化工问题。本章还希望通过 展示 COMSOL Multiphysics 和 MATLAB 工具在化工过程分析中的强大功能，激 发读者对使用 COMSOL Multiphysics 进行建模与仿真的兴趣。 由于 COMSOL Multiphysics 不是一个通用的问题求解工具，所以一些目标需 要迂回实现。作者在使用 FORTRAN、Mathematica 和 MATLAB 解决化工问题方 面有着丰富的教学经验，并用这些工具实现过这里所有的例子。而且，扩展化工 问题的数值分析也已经在 POLYMATH[1]中实现，这似乎只在化工委员会的 CACHE 项目中使用过。 本书前一版已经介绍过在零维空间中求解非线性代数方程和与时间有关的 常微分方程的内容。从概念上讲，零维域就是一个简单的有限元。通过研究某一 特定有限元中的变化对理解有限元方法非常有用。但是，COMSOL Multiphysics 通过独立对话框设置，使得零维几何方程和与时间相关的常微分方程求解变得非 常简单。所以本章将同时采用这两种方法求解这些例子。 2．方法 1：求根 典型的数值分析课程会讲解多种求根方法，但是从实际经验来看，只有两种 算法非常有用——二分法和牛顿法。我们这里没有列出所有方法，而是重点考虑 为什么求根是最有效的数值分析工具。在线性系统中求根非常简单，但是对于非 线性系统这就是一个挑战，而所有感兴趣的动力学问题几乎都是非线性系统。对 非线性系统的求根起源于对反函数的描述。为什么呢？因为对于大多数非线性函 数，“正向”y=f(u)很好表示，但是它的反函数 u=f-－1(y)可能不能显式表示、多值 （无意义）或根本不存在。如果反函数存在的话，求解反函数其实就是求根的过 程——求解满足 F(u)=0 的 u 等价于求解 F(u)=f(u)-y=0。因为大多数数值分析的 目标是在系统约束下计算求解，所以这也等价于对所有的约束取反。COMSOL Multiphysics 拥有求解非线性问题的核心函数——femnlin，本节主要介绍用它求 解零维非线性问题。 femnlin 函数使用牛顿方法求解，由于只有一个变量 u，牛顿法通过对一阶 倒数 F (u)迭代来求根。该方法首先估计函数的斜率范围，然后再逼近根。该斜 率可以通过理论分析（牛顿－拉夫逊方法）和数值（正割法）方法求得。如果能 用任何一种方法求得斜率，就可以用泰勒定律来逼近根。其基本思想就是使用目 前猜测值 u0 的泰勒展开式： f u ? f u ? u ?u f u ? (1) ( ) ( ) ( ) ( ) 0 0 0 该公式可以化简，忽略(u-u0)的高阶项，计算根如下： f (u ) u ? u ? 0 (2) 0 f (u ) 0 这个方法可以快速地扩展到多维求解空间，例如将 u 看作未知矢量，“被 f (u ) 除”看作“乘以 f 的雅克比矩阵的逆”。下一节介绍 COMSOL Multiphysics 0 中的求根过程。 2.1 求根：COMSOL Multiphysics 非线性求解器的应用实例 如上节所述，求根本身是一个“零维”活动，至少对于“空间-时间”系统 多维未知矢量 u 来说是这样的。COMSOL 多物理场没有零维模式，所以我们临 时采用一维模式。这在方面增加了我们不需要的冗余功能。但是由于问题规模较 小，COMSOL Multiphysics 编码效率高，且现代微处理器的运算速度快，这点就 不成为问题了。 启动 MATLAB 并在命令