matlab 周期性边界,Matlab PDE工具箱周期性边界条件

我能够在2016a中实现周期性边界条件，以便在pde工具箱中使用。这花了不少时间，因为没有记录太多的函数和数据结构。请注意我的域是一个三维立方体，但这里的方法也适用于2D。 (我已经测试过抛物线和双曲问题)。似乎应该通过改变PDE工具箱的“H”矩阵来实现周期边界条件，但是我还没有能够得到这种方法来处理椭圆，抛物线或双曲问题。

下面是用于实现抛物线和双曲线问题周期性边界条件(使用PDE工具箱线的方法)的方法：

使用点的Delaunay三角网创建一个网格(确保接触域的边界的节点具有在相反面上“匹配”的节点)。请注意，不能使用pde工具箱为3D立方体生成的网格节点，因为它们没有适当的周期性结构。

将周期性边界条件视为时间相关的狄利克雷边界条件。你可能会遇到一些麻烦让matlab在时间依赖的边界条件下处理(我不得不在这里忽略细节)。然而，在这里的是，我在applyBoundaryCondition()函数所使用的功能如下：

函数[bcMatrix]

％用于楠= tdependentdiri(地区，州)，以 “欺骗” MATLAB成处理的边界条件为时间依赖的狄利克雷

如果(isnan(state.time))

bcMatrix = NaN;

别的

bcMatrix = state.u;

端

端

更新FAS数据(内generateMesh从genmeshinternal()())。不幸的是，我们无法看到genmeshinternal())的代码。然而，fas是包含六个条目的单元格，其中每个条目都是2xN矩阵。每个条目的第二行。第一行指定一个面在域边界上的四面体。第二行指定哪个节点(四面体)位于域的内部。编写一个算法来更新网格的fas。实际上，您需要更新以下行(在generateMesh()内)中的每个变量(使用网格数据)：

[nodes，fas，tet，Hmax，Hmin] = genmeshinternal(self。几何形状，HMAX，HMIN，geomOrder);

请注意，节点和tet由delaunay三角剖分照顾。

节点和边界值(随时间变化)的“周期性映射”应在callValueFuncOnFace()内完成。基本上你需要在它的周期性合作伙伴的位置(由pbc_indx这里表示)来改变边界上当前解决值的值：

state.u = self.uN(：，pbc_indx); ％self.uN保存当前时间的解决方案

我把这个正确的行后面写着： appRegion = self.applicationRegion(xyzPts(：，i)，faceNormals(：，i));

可能有一个更简单的方法，但至少这个最终的作品。