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

LBM: 格子单位与物理单位间的转换

Posted by on 04/8/2012 in LBM FAQ | 45 comments

1. 对于无量纲系统 考虑无量纲系统:假定格子数为 ,时间步数为 , 离散空间和时间尺度 和 分别给出如下: 此时,和的关系给出如下: 同时,我们有 那么,格子速度和黏性决定如下: . 由于无量纲速度为1, 所以我们有 对于,真实物理系统如法炮制,只需根据物理,重新定义,然后得到. 2. 真实物理系统 网友问题 物理参数为:=101325Pa =101225Pa, 压强差为100Pa Lx=25mm Ly=1mm,求出特征速度=13.9m/s, Re=920 空气 粘度为1.5X(10^-5),划分为250X10 问题:如何确定格子参数? <1>....

Read More

Luo L.S.与 Karlin I.V.之间的激烈争论

Posted by on 12/21/2011 in LBM理论 | 2 comments

在Lyon参加ICMMES-2011期间,我录制了Luo L.S.与 Karlin...

Read More

LBM: 相变临界参数

Posted by on 12/19/2011 in BGK-LBM, LBM FAQ | 4 comments

这里提到相变,主要指一阶相变。在日常生活中,我们会常常遇到的,比方水变为水汽。这类相变最重要的性质是系统的对称性发生骤变,系统在两种状态之间转化。此类问题以密度作为序参数,涉及非连续性行为。针对LBM,Shan-Chen在1993年提出了多相及多组分模型。对于此类问题,使用LBM模拟多相流的朋友大多有个疑惑,就是如何确定状态方程中的模型参数,在LBM: Palabos-OpenLB群中已经做过解释,在此,我对临界点再做解释,以备初学者了解。 Shan-Chen状态方程如下: 式中: . 临界点的确定: ,...

Read More

LBM常见问题汇总

Posted by on 11/29/2011 in LBM FAQ | 1 comment

我们会根据各位同行的需要,结合自身的理解和经验,整理出一系列的LBM常见问题及其解答。 有兴趣的同行,也可以直接回复本帖,提交各自关心的问题。

Read More

LBM与湍流

Posted by on 11/25/2011 in BGK-LBM, MRT-LBM | 2 comments

...

Read More