我有一个表示三维矩形空间的稀疏矩阵。在一些边界上,我知道值是什么(它是一个常量)。其他边界可以是反射边界、差分边界等。
我是不是应该把这个问题设置得好像所有的边界都是微分的,然后回过头来,把解向量b中的节点设为常量?
谢谢!
发布于 2010-06-25 23:07:45
在有限元方法中,您可以区别对待Dirchelet (值约束)和Neumann (导数约束)。通常,在不考虑边界条件的情况下组装矩阵,然后应用边界条件,然后进行LU分解来求解。
您可以通过修改组合矩阵和RHS矢量来应用边界条件。我必须知道更多细节才能确切地告诉你你需要做什么。
https://stackoverflow.com/questions/3119043
复制相似问题