LBM 模拟多孔介质流动边界处理

城市阳光

刚接触LBM，使用matlab编程，然后遇到边界问题，D2Q9  模型中，9个格子速度是在常规坐标系下建立的，表示颗粒和孔隙的矩阵，在matlab中，矩阵读取是从上往下，从左向右读取的，跟我们常规的坐标系不一样，是否需要把格点速度转换成matlab中坐标下的速度，问题不知道是否描述清楚了，希望大家给些指导。

seanwillian

不需要。你可以实际的处理比对一下。

城市阳光

你好，因为是多孔介质，存在很多内部边界（矩阵中，颗粒用1表示，孔隙用0表示）我需要判断是否是内部边界，即是否是颗粒，若是颗粒（1），则执行反弹边界，假设遇到颗粒1，坐标为 （i,j），它的上面3个粒子均是孔隙，反弹速度应该是2，5，6,（按照常规坐标系建立的速度），那么上面3个粒子的坐标是按照常规坐标系建立的设计（跟速度坐标系一致），还是按照matlab 矩阵坐标系建立，谢谢你的指点，方便的话，能否留个邮箱或者联系方式，想请教一下

seanwillian

（i,j）和实际的（x，y）是一一对应的。只需要将行操作看作是对x的处理，列操作是对y处理就可以了。总之，只要这个对应关系确定了。无论你用什么来写代码都是一样的。