renxiaosa00

> 亲爱的杨博士， 我是一名多孔介质多相流LBM学习者。你的工作非常出色，我们可以花更多的时间来研究物理问题本身。 当我运行文件夹“2phase”中的“lbm_solver_3d_2phase.py”时，发生了一些奇怪的事情。10000步后，在出口处出现侵入相，但此时侵入相并未突破。这是由于使用周期性边界条件时体力较小造成的吗？ > >  您好，当运行文件夹“2phase”中的“lbm_solver_3d_2phase. py”后，参数都是默认的吗？结果的展示是在paraview里面看的，但是只有8个不同形态，请问参数这些都是默认的吗还是有什么别的参数设置或者是我查看的方式有问题啊？我这边没有形成这样的流动形态。谢谢您！

> solid_np，phase_np = init_geo('./img_ftb131.txt' ，'./phase_ftb131.dat'）这是定义几何形状和初始水/油分布的线（phase_ftb131.dat，1=水，-1=油） 您好，请问您是如何将数据转换为需要的solid_np，phase_np = init_geo('./img_ftb131.txt' ，'./phase_ftb131.dat'）这是定义几何形状和初始水/油分布的线（phase_ftb131.dat，1=水，-1=油）这2个文件的呢？img_ftb131是先按照切片数为第一层循环，然后按照行为第二层循环，再按照列为第三层循环进行读取CT的数据嘛？同样的，phase_ftb131相分布数据，是同样的顺序，然后将第一个切片的像素为0的位置设置为1，其余位置是-1，其余切片不用设置1和-1嘛？

img_ftb131和phase_ftb131的matlab生成如下： clc;clear;close all; % 将二值化0和128的bmp转换为2D的imgftb300 % 将每一个CT大小的前6行赋值为1形成1D的phase_ftb300 ImgPath = "自己的路径" matFiles = dir(fullfile(ImgPath,'*.bmp')); imgftb300 = []; phase_ftb300 = []; phase_ftb300_i= ones(300, 300);% 数字自己决定大小 phase_ftb300_i(1:6,:) = -1;%原始测试数据是前6行 lens = size(matFiles); for...

> 杨博士您好， 我计算了水相和气相的饱和度，但是当我达到200步时，突然遇到错误，并且输出vtk文件也不正确。请问这个问题是什么原因造成的呢？ > >      What is the purpose of calculating the water lattices and saturations? This is to calculate these parameters whether it converges...

> 看起来模拟配置中有些不太可行......你能可视化初始空气和气体的分布吗？它们的粘度是多少？您用来驱动流体的体力是多少？ > [……](#) > On Sat, 2 Sept 2023 at 09:42, Vivia93 ***@***.***> wrote: Hello Dr. Yang, I calculated the saturation of the water and gas phases, but when...

我是计算二氧化碳进入水的两相流模拟过程，数据大小也是150*150*150

好的，谢谢您！ （1）“体力1e-3”指的是fx fy 嘛？ （2）周期边界加固定颜色指的是边界条件的bcxr嘛？比如说以下图片的设置这块嘛 

您好，我按照您的建议设置参数如下：  得到的phase结果迭代0时如下：  中间过程结果如下：   最后的结果是：  假设入口是二氧化碳蓝色，水是红色，但流动的时候，理论上应该是蓝色二氧化碳慢慢驱出红色水，但是这个怎么出现的是灰色的相呢？请问这个是我的结果计算错了嘛还是出现什么别的设置的问题了？最后谢谢老师，期待您的指导！

  我设置了threshold，但是好像选择的opacity不太对，只要这块不是1，图像就没了，似乎是没找对geometry部分opacity，请问这个是在哪里找到的呢？我搜的可视化方法基本上很少有这个介绍。

好的，谢谢您的指导，我属实是第一次接触多相流，真是太感谢您的指导了！就是我的出口总是在刚开始迭代的时候就产生相的变化，按理说刚开始二氧化碳进入孔隙的时候左边入口才开始慢慢移动，右边出口不该出现二氧化碳的，但是无论我怎么尝试，刚开始进入二氧化碳时，出口总是有二氧化碳存在，之前也听您解释过这个原因，但是改参数总也不知道哪些值往哪些方向改，就有些不知所措了