我用多网格法解了二阶偏微分方程后,发现其数值解都收敛到了负数,而其实际的物理意义应该是非负数的。因此要对其边界点作处理,要用到无通量情况下的边界条件。有那位大侠帮我,没齿难忘!!!提供有关信息也给分!!!!

解决方案 »

  1.   

    用MATLAB把!
    然后把它嵌近来!
      

  2.   

    我们这个学期学了偏微分方程数值解法
    具体的讲要看是什么方程
    你用的应该是有限元吧
    我们后面讲到了椭圆型和抛物型方程的有限元方法
    但我没用过
    我编了很多差分方法的程序,关键是收敛性和稳定性的判断
    最好选用稳定性好的方法,各种隐格式虽然计算复杂些,但稳定性都是比较好的
    如果求解区域是矩形,那用差分就比较方便了
    若是很复杂的曲边区域,最好还是用有限元如果用软件,其实matlab的功能很弱的,方法也少得很,不专业,很多时候由于方程和格式性能,不是增加或细分节点就能解决的,因为本来格式可能就不收敛到准确解,最好还是自己选用好的方法。不过,做有限元的有很多通用的软件,像ANSYS之类的,或者更专业的商业软件,工程方面很多人都用这些东西,虽然还是不如自己编程放心,但省事多了。