LBM-IBM-DEM进阶之路4----Two-relaxation time (TRT) LBM 双松弛时间LBM（之1）

CFD升级打怪的龙猫

已于 2024-07-03 21:53:11 修改

阅读量1k

点赞数 17

文章标签：算法学习

于 2024-07-03 21:42:13 首次发布

本文链接：https://blog.csdn.net/2401_82853022/article/details/140152653

版权

本节我们来看双松弛时间LBM，TRT LBM的执行方法。

双松弛时间（TRT）模型是对传统的单松弛时间（SRT）和多松弛时间（MRT）模型的一种改进。TRT的核心思想在于：

松弛时间的分离：TRT模型将粒子分布函数分解为对称和反对称分量，并分别使用不同的松弛时间（ $\tau^+$ 和 $\tau^-$ ）进行松弛。这种方法允许更精细地控制扩散和对流过程的数值稳定性和准确性。
更高的数值稳定性：通过分别调节对称和反对称分量的松弛时间，TRT模型可以在保持高数值精度的同时，提高系统的数值稳定性。

1. 背景及初始化

TRT LBM中溶质的热力学状态由一个Q维粒子分布函数 $f_q(\mathbf{r}, t)$ 定义，其中q=0,…,Q−1。该函数在每个格点( $\mathbf{r}$ )和每个离散时间( $t$ )上定义。格点通过一组与格子轴和对角线对齐的离散速度 $\mathbf{c}_q$ 相互连接。对所有离散速度求和粒子分布得到溶质浓度 $C$ ：

$C = \sum_{q=0}^{Q-1} f_q$

在每个节点上，粒子分布函数可以分解为对称和反对称分量，

$f_q = f_q^+ + f_q^-$

其中， $f_q^+ = \frac{f_q + f_{\bar{q}}}{2}$ ， $f_q^- = \frac{f_q - f_{\bar{q}}}{2}$

TRT LBM的物理解释:

对称分量反映扩散：对称分量 $f_q^+$ 反映了粒子分布的平均行为，可以与经典的扩散现象相对应。
反对称分量反映对流：反对称分量 $f_q^-$ 反映了粒子分布的偏移行为，可以与经典的对流现象相对应。

整数 $\bar{q}$ 是指向 $\mathbf{c}_q$ 反方向的速度索引 $(\mathbf{c}_{\bar{q}} = -\mathbf{c}_q)$ 。

对于零速度的其余粒子 $(\mathbf{c}_0 = 0)$ ， $f_0^+ = f_0 = C - \sum_{q=1}^{Q-1} f_q$ ， $f_0^- = 0$

为了简化对所有速度 $\mathbf{c}_q$ 的求和，假设速度按以下顺序排列：

$\bar{q} = q + \frac{(Q-1)}{2}$ ， $1 \leq q \leq \frac{(Q-1)}{2}$ ，这里举例：对于D2Q5, 3 = 1 + (5-1)/2; 4 = 1 + (5-1)/2;

$f_0^+ = C - 2 \sum_{q=1}^{(Q-1)/2} f_q^+$

2. 碰撞迁移方程

一个节点上的粒子通过非零离散速度 $\mathbf{c}_q$ $(\mathbf{c}_q \neq 0 \ldots, Q-1)$ 移动到邻近节点，或者对于 $\mathbf{c}_0$ 保持在节点上。一旦粒子到达邻近节点，结果分布通过TRT碰撞算子放松到平衡状态：

$f_q(\mathbf{r} + \mathbf{c}_q, t + 1) = f_q(\mathbf{r}, t) - \frac{1}{\tau^+}(f_q^+ - e_q^+) - \frac{1}{\tau^-}(f_q^- - e_q^-)$

其中 $e_q^+$ 和 $e_q^-$ 是平衡粒子分布的对称和反对称分量 $(e_q = e_q^+ + e_q^-)$ ； $\tau^+$ 和 $\tau^-$ 分别是对称和反对称松弛参数； $f_q(\mathbf{r}, t)$ 表示在时间 $t$ 时，位置 $\mathbf{r}$ 处的粒子分布函数； $\frac{1}{\tau^+}(f_q^+ - e_q^+)$ 是对称分量的碰撞项， $\frac{1}{\tau^-}(f_q^- - e_q^-)$ 是反对称分量的碰撞项。

最低0.47元/天解锁文章