FEniCSx入门——解Poisson方程

一、 Poisson方程

 可转换为线性弱形式

 其中，

u: trial function 是近似解

v: test function 是Poisson方程两边同乘的函数，用于转换为弱形式。

 

在Fenics中解PDE时，必须进行下面几个步骤：

1、将PDE问题转换成离散变分问题：即寻找u∈V，使得

 2、选择V和V_hat空间，换句话说就是确定网格和有限元单元的类型。

 3、定义Input data，如f，u0。

 

二、Code in Fenics

这里定义 f 和 u0 为：

 

下面是示例

 

 

 

 

参考：来自fenicsbook，代码实例有点老旧