求解偏微分方程开源有限元软件deal.II学习–Step 12

Posted on 2016-09-17   |   In computational material science   |   暂无评论

引子

本例主要讲解MeshWorker框架和间断Galerkin方法，即DG。

• 用间断Galerkin法离散线性对流方程
• 使用MeshWorker::loop()来组装系统矩阵

本例主要关心的就是间断Galerkin法的循环，这相当复杂，因为必须分辨边界、常规内部边和有悬点的内部边。MeshWorker框架能够对所有的单元和边进行标准循环，它将分辨过程隐藏在了内部。
使用MeshWorker需要手动做两件事：一是针对特定问题写内部积分器，二是从该命名空间中选择类，然后将它们组合起来来解决问题。

要求解的问题是线性对流方程：

边界条件是：

这是入流边界，定义是：

这个方程是之前的Step9的守恒版本。
具体取值为：

，代表一个环形逆时针流场，在上，而在上

。
这里使用迎风间断Galerkin方法。
从这里就完全看不懂了。留坑待填。

#deal.II