演化方程的有限元分析.doc

第二章 发展方程的有限元分析 W. B. J. ZIMMERMAN，B. N. HEWAKANDAMBY Department of Chemical and Process Engineering, University of Newcastle Street E-mail: w.zimmerman@shef.ac.uk 科学研究和工程应用中的偏微分方程（PDE）多源自复杂的平衡方程。常见的偏微分方程主要来自质量守恒、动量守恒、组分守恒和能量守恒定律。由于这些守恒定律是整个域上的积分方程，所以在连续性假设下，偏微分方程很容易用有限元方法近似描述。本章介绍了COMSOL Multiphysics中典型的三种不同类型“时间－空间”系统偏微分方程——椭圆方程，抛物线方程和双曲线方程。本章还对有限元方法进行了总体介绍，结合应用实例讲解有限元方法精确计算的特性，更深层次的内容将在后续章节中引出。 1. 简介 在科学研究和工程应用中常会遇到满足守恒定律约束的偏微分方程，通常以积分形式出现。所有由质量守恒、动量守恒、组分守恒和能量守恒控制的传递现象都会产生连续逼近的偏微分方程。相信化工人员对传热、传质和动量传递现象不会感到陌生。 与前一章COMSOL Multiphysics化工实例中介绍的零维、一维空间系统相比，化工课程中通常不会出现超过二维或三维的偏微分方程计算。从文献[1]中找到一个非常珍贵的例子。实际上，很多常见的化工模型和公式都是由实际过程中更高维数的动力学过程简化而来的。流体动力学中的阻力系数，传热传质系数，多相催化的Thiele模型，精馏塔设计中的McCabe-Thiele图等许多描述高维数系统传递现象或非稳态动力学过程的技术，都是半经验性的方法，也许可以用偏微分方程来描述这些过程，但是由于基本物理、化学现象的复杂性，这些方程通常很难求解。所以对于初步的设计计算，这些快速计算的简化方法很受欢迎，但是对于细节设计、设计翻新、过程分析和优化过程，只有简化方法是不够的。在基础科学研究过程中，这些方法仍然不断地从化学工程师传给生物学家、材料学家等，逐渐应用到各个领域的工作中。但是计算流体动力学（CFD）的出现彻底改变了这些方法在传导模型中的地位。虽然CFD在传导现象可视化和量化方面具有独特优势，但是唯象方法对于描述分布系统模型仍然有非常重要的作用。 COMSOL Multiphysics不是一个“商业CFD软件”，但是也可以做一些CFD计算。它包含一些通用的CFD软件包，在某些模型支持上具有其独特的优点。对于CFD方法，大多数过程工程师希望有对湍流和燃烧模型的支持。但是COMSOL Multiphysics不同，擅长多物理场计算。除了CFD的传统传递现象，COMSOL Multiphysics还包含了电动力学、磁动力学、结构力学等实例模式，可以对这些现象进行模拟计算。实际上COMSOL Multiphysics的最大优点在于——首先，用户自定义编程容易，可以轻易建立用户自定义的模型，通过改变变量系数、边界条件、初始条件，还可以同时耦合多个物理场（甚至不在同一个域上的物理场）；其次，基于MATLAB（或COMSOL Script），可以实现对复杂模型模拟的所有编程功能，可以把COMSOL Multiphysics看作是一个方便的高阶有限元编程和分析软件。上一章中，我们介绍了一些用户自定义编程分析的功能，本章将介绍COMSOL Multiphysics在更高维数偏微分方程系统有限元建模上的核心优势。多物理场、扩展多物理场和无PDE约束的处理将留给后续章节。 2． 偏微分方程 偏微分方程可以根据阶数、边界条件类型和线性度（线性、非线性和准线性）进行分类。令人惊讶的是，大多数科学研究和工程应用中碰到的偏微分方程都是二阶的，即最高导数项是二阶偏导数。这是偶然的么？这一巧合往往是在各自领域分别解释的。但是，最近Frieden[2]证明所有已知物理定律都可以由最小Fisher信息原理推导出，而该原律通常将一个二阶项作为最高项引入相应的物理定律中——包括从Schrodinger波函数方程到经典电动力学等各种物理定律。所以二维和三维的二阶“空间－时间”系统求解和分类在科学研究和工程中有着广泛的应用，这非常重要！出于这个原因，且有限元方法本质上非常适合处理二阶系统，所以有限元是有着广泛应用的技术。 本章我们将集中介绍二维和三维空间的二阶系统。这里有三个典型的例子，在课本中经常见到，它们是： 下面用COMSOL Multiphysics对这几种偏微分方程进行建模，同时在每个例子中我们将引入一个新的或特殊的COMSOL Multiphysics功能。这些椭圆、抛物线和双曲线方程可以归结为一种通用的线性二阶偏微分方程，具有两个独立变量和一个非独立变量：