第七章 变量耦合：反问题、线积分、积分方程.docVIP

第七章 变量耦合：反问题、线积分、积分方程 和积分微分方程 W. B. J. ZIMMERMAN Department of Chemical and Process Engineering, University of Newcastle Street E-mail: w.zimmerman@shef.ac.uk 本章进一步深入探索耦合变量，考虑如何利用耦合变量求解反演方程和各种类型的积分方程。本章选择了四种重要的应用实例——使用激光探测位置，污染云层的传播，电容层析成像反演问题，纤维媒介中的非局部传热以及颗粒处理中的聚集平衡方程。我们将用到之前在COMSOL Multiphysics中没有用到的新特征——通过MATLAB与优化工具箱耦合，扩展网格，使用瞬态求解器作为稳态非线性模型的迭代工具，以及在耦合模式下选择激活/非激活多物理场模式的能力。后者对于单一耦合问题非常有用（如第3章中介绍的球形颗粒表面催化剂周围的水动力学问题）。为了处理本章中的积分方程，设置辅助域上的虚拟因变量。耦合变量在该域进行各种运算，但是因变量本身并不需要。因此，可以更好的调节对积分方程的有限元近似。 1. 引言 我们已经对边界和域积分较为熟悉了——可以在COMSOL Multiphysics的后处理菜单中选择。边界积分对计算各种表面量非常有用：例如静电体表面电荷，流体中的曳力等。域积分常用于平均和在域中定义的较高程度自由度的结合。这些特征是可信的，并给出了通过积分形式表达的有限元方法，伽辽金方法（见第2章），COMSOL Multiphysics包含了有效和精确的积分算法。因此，如果读者对任意积分阶数的数值积分或者常微分方程组（ODES）感兴趣，可以查阅MATLAB的内置函数（见第1章），这里不再赘述。在本章中将讨论积分方程及其理论的更多复杂应用，均采用COMSOL Multiphysics求解。线积分、积分方程以及积分－微分方程是应用数学研究的对象，也是COMSOL Multiphysics之所以具备变量耦合能力的基础。这也是为什么本章要取名为“多物理场耦合问题的伸展”。尽管题目选定了，但这只是说明性的，我们仍然以解决化工过程中的问题为目标，对COMSOL Multiphysics的特征进行深入了解。这里最重要的是使用激光雷达定位以及气体污染云层的扩散，聚集平衡方程（通常为微分方程的初值问题），电容层析的反演问题，以及纤维媒介中的非局部传热问题的计算。 扩展多物理场 在第一个方法中，我必须承认自己曾对我们要研究的扩展的多物理场问题保持怀疑的态度。最终，标量耦合变量吸引了我，并提供了撰写第4章的素材。它也促使我的课题组对微流网络这样一个新的领域产生了兴趣。目前，COMSOL Multiphysics提供了其他两种耦合变量的概念－拉伸和投影耦合变量。关于使用他们的例子均与后处理有关，例如，用于在域的截面函数性的表达求解结果以分析特定特征。 尽管如此，在一些较少出现的情况下，耦合变量（拉伸和投影）会应用到模型中与独立域变量同时求解。多个域、多尺度以及多过程模型对于工程数学和数学物理并不常见。通常情况下，模型是局部的－包含于一个特定的微分或偏微分方程组中，而方程组的边界条件和初始条件均较好的给出了定义。这些是局部连续性模型。从历史观点来说，已经从分析技术使用的角度分析了这方面的发展，在一些领域上已经涉及闭式模型的种类。即使是能够采用连续方法处理的模型，仍然需要根据离散方法进行近似。连续性颗粒水动力学[1]和离散方法[2]越来越受到关注，但是较老的方法，如分子动力学模拟[3]，蒙特卡罗方法[4]，微观水动力学[5]，分区自动化[6]和气体/等离子体的精确数值解[7]为连续系统/离散系统建立了桥梁。另一个技术策略基于满足偏微分方程约束的优化理论－罚不可行自由度。混合积分非线性规划[8]，遗传算法[9]以及遗传规划[10]都对混合离散/连续系统较为适用。COMSOL Multiphysics倾向于以PDE约束下的连续性模型进行方程描述。因此，从概念上讲，物理场伸展并不是处理较粗糙时的弥补性思考。它可以允许从基础上处理离散系统，而且不仅仅是偏微分方程约束的问题，而且也可以处理积分约束问题。实际上，耦合变量允许非局部和离散模拟。 在第3节（Scalar），第4节（projection）以及第5节（extrusion）中，我们重新回顾了耦合变量来研究如何使用COMSOL Multiphysics处理反演问题，线积分以及积分/积分－微分方程问题。 标量耦合变量 毫无疑问，标量耦合变量从概念上最容易理解。在第4章中，标量耦合变量用来连接非均相反应器中的再循环流体－反应器的输出，通过合理改变比例，重新进入原料流中。建立了零维单元和第二个几何体用于模拟缓冲容器，该容器获