COMSOL3.5教程案例——使用PDE模式建立数学模型

《COMSOL3.5教程案例——使用PDE模式建立数学模型》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、w页码，1/8(W)基于公式建模的PDE模式：案例—基于公式建模案例—基于公式建模本节的模型求解在一个单元圆盘上的Poisson方程，首先设定常数源项，然后在圆心设定一个点源。两种情况下均已知解析解，可用来判定数值结果的精确性。这些模型作为在COMSOLMultiphysics中使用PDE的基于公式建模的介绍性案例，在COMSOLMultiphysics模型库的“基于公式建模”的章节中可以找到更多的相关案例。单元圆盘上的Poisson方程一个众所周知的经典PDE是Poisson方程在单元圆盘Ω上，f=1，其解析解为因此可以将COMSOLMultiphysics的数值

2、解与解析解进行对比。结果下面的图显示了误差(由COMSOLMultiphysics得到的数值解和解析解之间的差)：mk:@MSITStore:C:UsersAdministratorDesktopCOMSOL_LibDoc_Multiphysics.C...2011/5/26w页码，2/8(W)有限元解在求解域的中心很接近解析解，这是很自然的，因为解析解和单元形函数都是二阶多项式。然而，靠近边界的误差就比较大，因为边界上的单元形状与几何相似，局部和整体坐标不是线性。因此，尽管在局部坐标中是二阶多项式，形函数在x和y上并不是严格的二阶多项式。模型库路径：COM

3、SOL_Multiphysics/Benchmarks/poisson_unit_disk使用图形化用户界面建模下面说明如何使用图形化用户界面得到上面的图。模型导航视窗1在模型导航视窗，从空间维度中选择2D。2在应用模式列表中，打开COMSOLMultiphysics>PDE模式，然后打开古典偏微分方程。3选择泊松方程式。4在单元列表中选择拉格朗日–二次(缺省的单元类型)。5点击确定。选项和设定1在选项菜单选择轴/格点设定。2输入下表中的轴限制。3点击格点标签。4清除自动选择框，输入下面的格点间格和特别点：轴格点mk:@MSITStore:C:UsersAdmi

4、nistratorDesktopCOMSOL_LibDoc_Multiphysics.C...2011/5/26w页码，3/8(W)x最小-2x间格0.5x最大2特别xy最小-1.5y间格0.5y最大1.5特别y几何建模1在绘图工具条上点击椭圆形/圆形(以圆心)。2点击右键，画一个圆心在(0,0)，半径为1的圆盘。物理设定边界条件模型中的缺省边界条件是u=0，因此不需要做改变。求解域设定缺省f为1，因此也不需要做改变。网格1点击主工具条上的网格模式按钮来初始化和显示网格。2点击主工具条上的细化网格按钮。求解点击主工具条上的求解按钮。后处理和图形化要以3D表面的形

5、式观察结果，点击绘图工具条上的3D表面图。将FEM的结果与解析解对比：1点击绘图参数工具条按钮。2点击表面标签。3在表面数据的表达式编辑框中键入u-(1-x^2-y^2)/4。mk:@MSITStore:C:UsersAdministratorDesktopCOMSOL_LibDoc_Multiphysics.C...2011/5/26w页码，4/8(W)4点击确定。施加点源对于在圆心处有一个点源的单位圆上的Poisson方程，其正式表达式为：其中δ是在圆心的Diracδ分布。解析解为，在圆心处存在一个奇异点。最简单的描述点源的方法是使用一个额外的弱项，为此

6、，用户需要了解一些FEM理论和PDE的弱解的基础知识。要得到广义Poisson方程的弱解形式(其中f是任意的源项)，乘上一个试函数u，并在求解域内Ωtest积分：分布积分，引入Lagrange乘子λ处理边界上的约束，得到：mk:@MSITStore:C:UsersAdministratorDesktopCOMSOL_LibDoc_Multiphysics.C...2011/5/26w页码，5/8(W)这些方程中的各项是域Ω内或边界∂Ω上的积分。可以向方程中引入一个额外的弱项。要处理在Ω或∂Ω上的积分，COMSOLMultiphysics同样要处理来自于边和奇

7、异点的贡献。软件不是积分来自点的贡献，而是直接将它们累加到方程中。只要源项f是一个空间、解变量和时间的函数，COMSOLMultiphysics就可以离散给定的积分方程。不幸的是，无法以COMSOLMultiphysics能积分的函数形式表述Diracδ分布。根据Diracδ分布的定义，下式成立：因此，设定f为0，并在原点上添加一个弱项u来修正弱形式。下面的案例模型采用了局部精细化网格test来处理位于原点的奇异点的解析度。模型库路径：COMSOL_Multiphysics/Benchmarks/point_source使用图形化用户界面建模模型导航视窗1在模型