renxiaosa00

我的颜色边界采用的默认的代码： bc_psi_x_left, psi_x_left = 1, -1.0 bc_psi_x_right, psi_x_right = 0, 1.0 bc_psi_y_left, psi_y_left = 0, 1.0 bc_psi_y_right, psi_y_right = 0, 1.0 bc_psi_z_left, psi_z_left = 0, 1.0 bc_psi_z_right, psi_z_right = 0,...

我看后面备注是 1 = constant value on the boundary, value = -1.0 phase1 or 1.0 = phase 2，那么是不是bc_psi_x_right, psi_x_right = 1, 2.0啊？毕竟刚开始的初始条件右边出口按理说得是phase 2-水，左边入口才是-1对应的二氧化碳吧？是不是变量我理解错了，比如bc_psi是设置颜色的，psi是设置别的什么参数的？其实我没太搞懂pis_x/y/z, psi_solid分别是什么意思。

> 感谢您的详细测试和说明。在我们的帮助下，我将 gcc 降级到版本 9，它解决了 ub22.04 上的问题。 How to set the stl file to apply multiphase?

作者您好，我想请教一下，这个是怎么输入单波段的CT数据的啊？因为我看数据都是默认初始化一个二维数组进行处理的啊，并没有外界的输入数据啊

好的，谢谢作者哦~

> 好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。   您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题

> 你可以多探索paraview里的功能，时间比较长了，我记得大概是设置的问题 > [……](#) > ---- 回复的原邮件 ---- | 发件人 | ***@***.***> | | 日期 | 2023年10月07日 13:52 | | 营养 | ***@***.***> | | 抄送至| ***@***.***>***@***.***> | | 主题...

> 请问你说的想要的效果是视觉上的效果还是流体分布的位置和example图上不一样？ 如果是视觉效果，那是paraview使用上面的问题，vtr 文件只提供可视化的数据，怎么渲染取决于如何操作paraview （例子图片钟可视化的是固体部分（使用了半透明的设置）和non-wetting phase 的分布，不是速度或者压力） 如果是这些流体分布位置不一样，那是模拟初始化数据的设置问题，可以通过变化压力差和接触角还有表面张力大小调整，如果你的模拟毛细数很小，比如 […](#) > On Sat, 7 Oct 2023 at 06:52, renxiaosa00 ***@***.***> wrote: 好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。 [image: 图像] [image: 图像] 您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题 — Reply to this...

> 好，对于你的问题： 1：现在默认的设置是毛细力主导的（毛细管主导）所以流动非常缓慢，毛细力远大于粘性力，所以速度尤其是润湿阶段的速度粒度，所以你不怎么看得，这个是正常的2：程序默认参数不是生成那个图的参数，那个图的毛细数想要更小，需要大的表面支撑出来加大的体力运行很长时间才可以得到，毛细数小时候其实对应的dt也细，另外真实实验钟毛数小实验本身物理时间也很长（界面移动也很慢），所以模拟需要的时间更长了。具体参数需要根据你的需求来定，如果你觉得移动的慢，可以加大体力，或者减小表面张力CapA，另外粘度也有一定影响 3：这个代码就是这篇文章的复现：https: [//www.sciencedirect.com/science/article/ abs/pii/S0309170808000511](https://www.sciencedirect.com/science/article/abs/pii/S0309170808000511) 可以参考这篇文章就知道这里面参数的意义了。另外i毛细数是一个无量纲数，你可以Google一下毛细管数就怎么计算了，比较简单就知道了:-) > [……](#) > On Sun, 8 Oct 2023 at 02:50, renxiaosa00 ***@***.***> wrote: 请问你说的想要的效果是视觉上的效果还是流体分布的位置和example图上不一样？ 如果是视觉效果，那是paraview使用上面的问题，vtr 文件只提供可视化的数据，怎么渲染取决于如何操作paraview （例子图片钟可视化的是固体部分（使用了半透明的设置）和non-wetting phase 的分布，不是速度或者压力） 如果是这些流体分布位置不一样，那是模拟初始化数据的设置问题，可以通过变化压力差和接触角还有表面张力大小调整，如果你的模拟毛细数很小，比如 wrote: 好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。 [image:...

> 如果你多试几种参数就可以找到了，数值模拟和实验也类似，根据实验结果调整再试。比如你举得不容易推进去，就降低CapA 表面张力，或者加大bodyforce, 或者加大接触角COS值， 如果数值不稳定了，就稍微增大粘度，减小bodyforce, 减小CapA, 都可以帮助提高数值稳定性 > […](#) > On Mon, 9 Oct 2023 at 10:56, renxiaosa00 ***@***.***> wrote: 好，对于你的问题： 1：现在默认的设置是毛细力主导的（毛细管主导）所以流动非常缓慢，毛细力远大于粘性力，所以速度尤其是润湿阶段的速度粒度，所以你不怎么看得，这个是正常的2：程序默认参数不是生成那个图的参数，那个图的毛细数想要更小，需要大的表面支撑出来加大的体力运行很长时间才可以得到，毛细数小时候其实对应的dt也细，另外真实实验钟毛数小实验本身物理时间也很长（界面移动也很慢），所以模拟需要的时间更长了。具体参数需要根据你的需求来定，如果你觉得移动的慢，可以加大体力，或者减小表面张力CapA，另外粘度也有一定影响 3：这个代码就是这篇文章的复现：https: //[www.sciencedirect.com/science/article/](http://www.sciencedirect.com/science/article/) abs/pii/S0309170808000511 可以参考这篇文章就知道这里面参数的意义了。另外i毛细数是一个无量纲数，你可以Google一下毛细管数就怎么计算了，比较简单就知道了:-) …… On Sun, 8...