COMSOL做不定积分

用COMSOL做定积分应该是比较容易的事情，今天介绍下如何做不定积分。举个简单的例子，本文用会到intop和dest算符。

通俗解释：函数g(r)是关于向量r在空间上的函数，欲求在不同r下g(r)的体积分情况。

为了方便验算，我们对一个球体进行体积积分。球体积和半径的关系是：

模型建立：

1. 建立3D稳态General Form PDE模型，ps.其实用那个模块无所谓。

2. 因为我们要计算球体模型，所以采用球坐标更方便，故再建立一个 Spherical System(Sys2)


3. 画个球出来，为了方便后处理建立一条实体边（半径）

4. PDE模块填入参数，为了方便演示，我们算一个固定值。将Source Term改为0，建立各个曲面的一个Dirichlet Boundary Condition，令其为1，计算域初始值为1.这样不管在没算，u的值一直为1.

5. Definitions下建立一个Intergration 1，默认标示符为 intop1，选择全域。

6. 定义变量 intSp = intop1(u*(sys2.r<=dest(sys2.r)))，其实和intop1(1*(sys2.r<=dest(sys2.r)))一样。

7.求解，后处理1D Plot

选择画的那条半径边，分别画出 intSp 和 4/3*pi*sys2.r^3

蓝线为计算值，绿色点图为理论值。

讲解已经非常详细，故不提供模型，本模型不提供建模指导。

本文源自：天乐树的微博