4 扩展多物理场.docVIP

第四章 扩展多物理场 W. B. J. ZIMMERMAN1，P. O. MCHEDLOV-PETROSSYAN2*, G. A. KHOMENKO2 1Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom 2Laboratoire d’Oceanografi, Cotiere du Littoral, ELICO, Universite du Littoral Coted’ Opale, MREN, 32, Avenue Foch, 62930, Wimereux, France *Permanent address: NSS Kharkov Institute of Physics and Technology, Ukraine E-mail: w.zimmerman@shef.ac.uk 扩展多物理场是COMSOL Multiphysics的一个特性，从概念上讲它非常复杂。扩展多物理场的概念将两个或者更多的计算域通过变量耦合联接在一起，这些变量可以是设定边界条件的变量，也可以是子域偏微分方程的系数。这些耦合变量可以通过对整个域内、边界值的子域（或边界）进行积分得到。扩展多物理场通常出现在物理现象的多尺度建模过程中——大尺度的模型通常耦合了较简单的变量或更低空间维数的小尺度模型。扩展多物理场广泛存在于过程工程中，但是各种现象多发生在不同域上，所以至少需要通过入口和出口条件进行联接，但是更多则是通过更加精细的过程积分联接起来。所以，模拟过程的所有的优化、设计、改进和控制都可以归结为减轻扩展多物理场的过程。用Simulink积分可以使得一些运行单元可以用分布偏微分方程模型处理，而其它的则用集中参量处理，但是对于非平凡层面的耦合，还是需要扩展多物理场建模。 1．简介 如果说上一章讲的多物理场对于读者是个新概念，那么扩展多物理场就更加陌生了。COMSOL Multiphysics化工模型库中充满了扩展多物理场的例子——尤其是化工模型库[1]中燃料电池模型的Maxwell-Stefan扩散和Monolith反应器这两个例子。本章通过标量变量耦合的简单联接方式来引入扩展多物理场的概念。现在COMSOL Mulitphysics通过常微分方程设置功能引入了辅助状态变量，但是从本质上讲，它不是给简单扩展多物理场引入了二级概念域。耦合变量可以轻易的联接相同的域或者不同的子域。第二章中的周期性边界条件（孤波例子）就是一个简单的扩展多物理场例子。 第七章对三种类型的变量耦合方式做了更详细的介绍，通过这些例子介绍了在COMSOL Multiphysics中是如何用有限元方法来处理卷积积分、积分方程、积分微分方程和逆问题的。第九章介绍了如何将扩展多物理场的概念应用到遗传基因网络建模中的，这里的“节点”是微尺度流动计算的连接点，通过“长”的微槽连接。 正如文献[1]所述，一开始的时候，我以为扩展多物理场是关于多尺度耦合模型的。这是多相流动/异相系统研究的前沿问题，因为消失的相可以看作是域中的一个点约束，但是向各个方向都有流动信息。通常希望通过变量来处理这种约束，即通过建模使小尺度现象依赖于已知的大尺度流动，反之亦然，这样来实现耦合。通常小尺度现象有其自身的复杂性，例如Grammatika和Zimmerman[2]的微尺度静力学模拟，同时考虑了体相动力学的求解，所以把小尺度动力学模拟的简单函数形式附加到大尺度现象中去。使用变量强制同步方法有一些缺点，但是都通过“过渡简化”集中了起来。幸运的是，这种模型可以通过同时在两个尺度上逼近的物理系统来实验验证。传统的湍流模型也高度依赖于多尺度建模过程。由于有专门的多尺度建模技术，使得扩展多物理场看起来也不是一个非常有用的功能。它处理复杂模型的优点需要通过高性能计算才能体现出来。 经过很长时间我才意识到化工过程中充满了扩展多物理场的例子。首先让我们从COMSOL的意义上给出扩展多物理场的概念：一个模型如果需要对场变量在两个或更多的逻辑不相交域上加以描述，那么它就是扩展多物理场模型。它们可能不是物理上不相交的域，因为为了保证问题的求解，物理场必须在每个域上耦合。COMSOL允许用户使用几个不同的几何结构/实例模型来组成扩展多物理场模型。 所以为什么说化工过程中处处充满了扩展多物理场现象呢？参看图1。 图1 运行单元线性排列的流程图 对于图1中的流程： “5bar和30℃条件下的环丙烷以1kmol/h的速度进入系统。在进入连续搅拌器之前，通过热交换器加热到500℃反应温度。反应器容积为2m3并保持反应温度为500℃不变。异构反应