10楼: Originally posted by nagami at 2013-06-15 08:05:43
用Newton-Leibniz公式，然后与gauss公式类似的处理方法，就能出来包含函数值与其一阶导数的边界项。散度定理是其高维类似物，同样的

11楼: Originally posted by yuahua at 2013-06-17 10:26:57
Newton-Leibniz公式不是用来求解积分的吗，我下面这样推倒了一下，没有出来一阶导数的边界项啊。

无标题.png
...

12楼: Originally posted by nagami at 2013-06-17 12:10:54
我自己写的，参考下，不保证正确性，离散完就是你的事情了。基就是一次Lagrange多项式。

1.jpg

2.jpg

3.jpg
...

13楼: Originally posted by yuahua at 2013-06-17 17:07:07
LZ，数学功底很好啊。还有个问题
1： 2.1.1公式中的n为法向量，如果是一维（x方向）问题的话，是不是直接把n换成x就行了呢？
2： 我们前面说的那个边界条件，就是我在上面7楼说的那个13.2.4  怎么和13.2.3联系起 ...

14楼: Originally posted by nagami at 2013-06-17 17:16:16
1:是的没错，区间，对于b点，n就是x；对于a点，n就是-x；在N-L公式中，已经隐晦的处理这个问题，f（b）-f（a），多出来的“-”就是n的方向；在这个公式里面我换好了；
2：你把Nuemann条件显示写出，直接带入就行。 ...

14楼: Originally posted by nagami at 2013-06-17 17:16:16
1:是的没错，区间，对于b点，n就是x；对于a点，n就是-x；在N-L公式中，已经隐晦的处理这个问题，f（b）-f（a），多出来的“-”就是n的方向；在这个公式里面我换好了；
2：你把Nuemann条件显示写出，直接带入就行。 ...

16楼: Originally posted by yuahua at 2013-06-17 22:12:21
更清楚的边界条件

2.jpg
...

12楼: Originally posted by nagami at 2013-06-17 12:10:54
我自己写的，参考下，不保证正确性，离散完就是你的事情了。基就是一次Lagrange多项式。

1.jpg

2.jpg

3.jpg
...

18楼: Originally posted by yuahua at 2013-06-17 22:31:09
对于这个瞬态PDE，您的Galerkin  method中，an是关于时间的函数，Φ是关于x的函数，  是不是这个意思吗？

1.jpg
...

19楼: Originally posted by nagami at 2013-06-17 22:35:05
对，这个你有经验，相当于把对偶形式往基上投影，每投影一次，产生n个系数值，关于an的；投影n次，得到系数矩阵...