求大神帮看下程序

huangtiefei

是这样的，我把何雅玲老师书上的方腔驱动流的程序该为Fortran语言，在计算的时候，误差在减小几十步之后就开始不对了，最后会出现错误，调试了很久，不知道错误在哪里，求大神指导，程序见附件。

qjx1278

我最近也刚接触LBM方法，将后面那个顶盖驱动方腔转化成fortran语言，发现误差也越来越大，最后导致发散，怎么回事呢

Entropy

Lattice Boltzmann Method Fundamentals and Engineering Applications with Computer Codes
看看这么本书，书后有很多Fortran90的例子，有BGK和MRT，也有一些解对流扩散方程的算例。整本书的结构也是非常不错的。可以在sina下到电子版的。