如何设置边界条件(Dirichlet和neumann)？

Question

我需要帮助理解PETSc例子的有限元法与SNES非线性求解器。是否有任何好的材料或来源来理解这些例子？(我已经阅读了用户手册，但没有我需要的信息)。

我将用PetscFE和DMPlex (用于非结构化网格)编写一个包含PetscFE和DMPlex的代码。该程序将处理具有非线性热源的传热方程.

我想选择具有两个dirichlet条件的边界条件

1. 矩形的上、下边界。
2. 两个neumann条件(矩形上的左右边界)。

的示例12，用户应该定义：

quadratic_u_2d作为精确解(为了提高解的收敛性) f_0 (测试函数项的被积(与负载向量‘f’相关)、f_bc_0 (测试函数项的边界被积(与通量‘g’相关))(我认为这可能是边界条件) f_1 (测试函数梯度项的被积(与u_xd相关)) f_bc_1 (测试函数梯度项的边界积分)(我认为它是零，我认为，g3 (测试函数梯度和基函数梯度项)(我认为这个项在任何情况下都是固定为1的)。

但是，没有办法单独设置每个边界条件。请给我一些更好的建议。

谢谢你读我的问题..。

Answer 1

将Petsc应用于流体流动问题的FVM方法，通过修改矩阵A及其相应的RHS值来设置边界条件。