这种东西通常要看模型的物理意义来决定外边界。

跟你采用的网格类型有关，如果是点中心的，直接就是第一个网格点吧，如果是块中心的，可以在内边界左边在多虚设一个网格，或者直接采用第一个网格近似，你可以看看张烈辉编的哪本油藏数值模拟的书，上面有具体的处理方法。其实这类模型我感觉你直接求解析解不就得了，为什么要数值解。

引用回帖:

4楼: Originally posted by 毅丝不挂 at 2017-09-11 09:49:01
跟你采用的网格类型有关，如果是点中心的，直接就是第一个网格点吧，如果是块中心的，可以在内边界左边在多虚设一个网格，或者直接采用第一个网格近似，你可以看看张烈辉编的哪本油藏数值模拟的书，上面有具体的处理 ...

非常感谢！这里其实我想问的主要是柱坐标下边界条件的处理。我查了很多，大部分的差分方法只是针对直角坐标系，极坐标系很少。而且关于极坐标系下的边界条件如果要是按数学物理方法去处理的话相比直角坐标也更麻烦，内边界要加一个有界条件，角度部分还要加一个周期条件。如果用拉普拉斯变换的话，对于这里的模型还好一点，反演很简单，如果要是稍微变一变，那反演起来就比较麻烦了，搞不好还得数值反演。
所以就想看看如果纯粹是用数值解的话该怎么弄，是不是更简单一些。谢谢
，