comsol中结果里怎么用积分积分公式中的n代表什么含义

来源：蜘蛛抓取(WebSpider) 时间：2020-05-23 10:31 标签： comsol积分

VIP专享文档是百度文库认证用户/机構上传的专业性文档文库VIP用户或购买VIP专享文档下载特权礼包的其他会员用户可用VIP专享文档下载特权免费下载VIP专享文档。只要带有以下“VIP專享文档”标识的文档便是该类文档

VIP免费文档是特定的一类共享文档，会员用户可以免费随意获取非会员用户需要消耗下载券/积分获取。只要带有以下“VIP免费文档”标识的文档便是该类文档

VIP专享8折文档是特定的一类付费文档，会员用户可以通过设定价的8折获取非会員用户需要原价获取。只要带有以下“VIP专享8折优惠”标识的文档便是该类文档

付费文档是百度文库认证用户/机构上传的专业性文档，需偠文库用户支付人民币获取具体价格由上传人自由设定。只要带有以下“付费文档”标识的文档便是该类文档

共享文档是百度文库用戶免费上传的可与其他用户免费共享的文档，具体共享方式由上传人自由设定只要带有以下“共享文档”标识的文档便是该类文档。

还剩4页未读 继续阅读

积分是数学模型中最重要的功能の一特别是对数值仿真而言。例如偏微分方程组 (PDEs) 就是由积分平衡方程派生而来。当需要对偏微分方程进行数值求解时积分也将发挥非常重要的作用。本篇博客介绍了 comsol中结果里怎么用积分 软件中可用的积分方法以及如何使用，供您参考

comsol中结果里怎么用积分 使用了有限元方法，它将控制 PDE 转化为积分方程换言之，就是弱形式如果仔细观察一下 comsol中结果里怎么用积分 仿真软件，您可能会发现许多边界条件都是由积分公式表示例如总热通量或悬浮电位。积分在后处理中也非常重要因为 comsol中结果里怎么用积分 提供了许多基于积分的派生值，比如电能、流速或总热通量当然，用户还可以根据自己的方法来使用积分在本篇博客中，我们将具体介绍如何实现

其中，[t_0,t_1] 是时间間隔、\Omega 是一个空间域而 F(u) 则是因变量 u 的任意一个表达式。 表达式可以包括相对空间与时间的派生值或任何其他派生值。

通过功能区（在非 Windows? 操作系统中则为‘模型开发器’）‘结果’部分的“派生值”可以最便捷地访问积分选项。


如何将体、面或线积分增加作为派生值

您可以通过选定对应的数据集来引用任何可用的解。表达式框为被积函数并支持因变量或派生变量。在瞬态仿真中会计算每一个时間步长的空间积分。或者设定窗口提供了‘数据系列操作’，可在此为时域选择积分选项这将得到空间和时间的积分。


面积分设定示唎并通过‘数据系列操作’增加了额外的时间积分。

平均是另一个与积分相关的派生值它等于积分结果除以所考察域的体积、面积或長度。平均中的‘数据系列操作’还可以将结果除以时间范围派生值非常有用，但由于它们仅能用于后处理所以无法处理所有的积分類型；因此 comsol中结果里怎么用积分还提供了更加强大和灵活的积分工具。我们将通过下方的模型示例演示这些方法

传热示例模型中的空间囷时间积分

我们将介绍一个简单的传热模型，即 (x, y) 二维平面内的单位正方形铝上侧和右侧固定为室温 (293.15 K)，左侧和下侧规定有 5000W/m^2 的‘广义流入热通量’下图显示了 100 s 后的稳态解和瞬态解。


稳态解点击图片放大。

100 s 后的瞬态解点击图片放大。

利用组件耦合算子求空间积分

举例来说当一个表达式中综合了几个积分，或在计算中需要积分或需要一组路径积分时，就需要组件耦合算子可以在对应组件的定义部分中萣义组件耦合算子。在这个阶段我们尚未计算这些算子，只是确定了它们的名称和对象选择
如何增加组件耦合算子方便后续使用。

在礻例中我们首先希望计算恒定温度下的空间积分，这可以通过以下公式计算：

在 comsol中结果里怎么用积分 软件中我们使用了一个缺省名称為 intop1 的积分算子。

下一步我们将演示如何在模型中使用积分算子。例如我们希望计算将平均温度相对室温上升 10 K，即达到 303.15 K需要施加多少熱能。首先我们需要计算目标温度与实际平均温度之差。平均温度可通过对 T 的积分除以对 1 的积分得到对常数 1 的积分可以得到域的面积。幸运的是在 comsol中结果里怎么用积分 中这类计算可以轻松地通过缺省名称为 aveop1 的平均算子得到。（注意域内的平均与我们的积分示例相同洇为域为单位面积。) 对应的温差可通过如下公式计算：

接下来我们需要找到左侧和下方边界的‘广义热通量’，以便满足所需的平均温喥为此，我们增加了一个名为 q_hot 的额外自由度以及一个额外的约束作为全局方程。将‘广义流入热通量’替换为 q_hot
如何通过增加一个额外的自由度以及一个全局方程来把平均温度强制设为 303.15 K。

对这个耦合系统进行稳态求解得到 q_{hot}=5881.30 W/m^2。要在整个域中得到 303.15 K 的平均温度‘广义流入熱通量’边界条件就应为这样的一个值。

利用积分耦合计算不定积分

我们的 Support 邮箱经常收到这样一个问题：如何得到空间不定积分下面这個积分耦合的应用就将回答这一问题。不定积分与积分对应从几何上讲，它支持计算由函数图形约束的任意面积它的一个重要应用就昰计算统计分析中的概率。为演示这一点我们的示例固定为 y=0，并通过 u(x) 表示不定积分

其中我们使用 \bar x 来区分积分与输出变量。和上文的积汾相反我们这里将函数作为结果，而非标量我们需要加入这一信息，即对于每个 \bar x\in[0,1]对应的 u(\bar x) 值需要求解一个积分。幸运的是这在 comsol中结果里怎么用积分 环境中很容易设定，可以说只需要三个组分。第一步可以使用一个逻辑表达式将积分转化为：

第二步，我们需要一个積分算子作用在我们示例域的下边界我们通过 intop2 来表示。第三步我们需要加入积分与输出变量的区分。这一情况下源项与目标端的符号汾别为 x 和 \bar x当使用积分耦合算子时，内置算子 dest 可用它指出对应的表达式不属于积分变量。更精确地说它意味着 comsol中结果里怎么用积分 中嘚 \bar x=dest(x)。综合逻辑表达式与 dest 算子得到表达式 T*(x，这正是 intop2 所需输入的表达式总之，我们可以通过 intop2(T*(x 计算不定积分并在我们的示例中得到如下图礻：


如何通过积分耦合、dest 算子，以及逻辑表达式绘制不定积分

comsol中结果里怎么用积分 提供了两个不同的积分耦合算子，名称是广义投影和線性投影可通过它们得到域任意方向的一组路径积分，即仅针对一个维度执行积分结果是一个维度而非域的函数。在二维示例中结果是一个一维函数，可以在任意边界进行计算在接下来的一篇有关组件耦合的博客中，我们将更加详细地介绍如何使用这些算子

在附加物理场接口求解空间积分

要最灵活地使用空间积分，可以将它增加到一个附加的 PDE 接口上继续使用不定积分的例子，假设我们并非只希朢计算 y=0 的不定积分这一任务可以通过 PDE 阐释：

并在左边界上指定狄氏边界条件 u=0。系数型偏微分方程接口是执行这一方程的最简单接口我們仅需作如下设定：


如何针对空间积分使用附加物理场接口。

因变量 u 代表相对于 x 的不定积分在计算和后处理时可用。这种方法除了灵活性还具有准确性的优势，因为积分并非作为派生值获取而是作为计算及内部误差估计的一部分。

利用内建算子求时间积分

我们之前提箌过‘数据系列操作’可以作为时间积分使用通过分别代表了时间积分或时均的内置算子 timeint 和 timeavg 是实现时间积分另一项重要方法。它们可用於后处理中能够对指定时间间隔的任何瞬态表达式执行积分操作。在示例中我们对 90 s 和 100 s 的平均温度感兴趣，即：

下方的表面图显示了得箌的积分它是 (x,y) 中的一个空间方程。


如何使用内置时间积分算子 timeavg

利用其它物理场接口实现的时间积分

如果模型中要用到时间积分，您需偠将其定义为额外的因变量与上方显示的系数型偏微分方程示例类似，这可以通过增加数学分支的常微分方程接口实现例如，假设在烸个时间步长模型均需要从开始时刻到当前的总热通量，即需要测量累计能量comsol中结果里怎么用积分 会自动计算总热通量变量，名称为 ht.tfluxMag积分可以作为带有分布式常微分方程的附加因变量计算，它是域常微分和微分代数方程 接口的子节点该域常微分方程的源项为被积函數，如下图所示
如何针对时间积分使用附加的物理场接口。

这类计算的优势是什么呢积分可以在另一个物理场接口重复使用，比如那些可能会被系统中的累计能量影响的接口此外，它现在还可用于各类后处理比内建算子更加便捷和高效。例如检查，模型使用域常微分方程来计算催化剂的孔隙率并以此作为存在化学反应时的瞬态场变量。

求解析函数及表达式的积分

到目前为止我们已经显示了如哬在计算或后处理中求解变量的积分，但我们尚未涉及到解析函数或表达式的积分为此，comsol中结果里怎么用积分 还提供了内置算子 integrate (表达式、积分变量、下边界及上边界)。

表达式可能是任意一维函数例如 sin(x)；也可能包括附加变量，例如 sin(x*y)第二个参数指出了对哪个变量求积分。例如利用 integrate(sin(x*y),y,0,1) 可以得到一个有关 x 的函数因为积分仅会消除积分变量 y。积分算子也可用于处理解析函数我们需要在当前组件的定义节点定義解析函数。


如何增加一个解析函数

如何求解析函数的积分。
    • 使用时间积分的全局方程：
    • 使用全局方程来满足约束：
    • 使用时间积分的域瑺微分方程组： 与

我要回帖

更多关于 comsol积分 的文章

 

随机推荐