Generalized LBM

ELBM: 熵格子Alpha求解C++程序

Posted by on 05/5/2012 in Generalized LBM, LBM理论 | 0 comments

这分享我使用的求解ELBM中非线性Alpha方程的核心代码,希望这对感兴趣熵格子方法的朋友有用!欢迎大家使用! //----------------------------------------------------------- #include <math.h> //----------------------------------------------------------- //External variables used in the computation. //----------------------------------------------------------- #define TOL            5e-16 #define   MAXIT        200 #define   TOLL          1e-7 // double     log_[9]={1./8,0.5,0.5,0.5,0.5,2,2,2,2}; double     F_eval_r(double); void         FD_eval_r(double,double *,double *); double   ...

Read More