Jianhui Yang

那应该是设置上有些问题让simulation crash了。。。 粘度不要太小，尽量使用bodyforce来驱动流体，压力边界使用上需要特别小心，数值上不是很稳定 On Wed, 6 Dec 2023 at 12:45, Vivia93 ***@***.***> wrote: > 杨教授，您好！我运行了单相流的程序，为什么我运行到24步以后，突然就不行了呢？运行过程中也没有报错。 > VTV.IZ6SRMLBR.FPN.J7Z5G.png (view on web) > > ***@***.*** (view on web) > > > —...

判断稳定这个事确实有一定的主观性，在我处理的样本里面，均质性比较好的砂岩，基本上跑一段时间你如果画出湿润相和非湿润相的速度随时间变化曲线就能看出来基本变化很小了，就是稳定了，但是也有些尺度跨度很大的碳酸盐，跑了很久还是能看到速度曲线一直波动很大，如果能看出来事周期性变化，那也可以认为事稳定了，如果不是，可能就还得继续让流动继续一段时间，或者可视化下油水分布，看看是不是分布合理，人工判断下。这个和样本的复杂度是高度相关的 On Wed, 26 Jun 2024 at 15:19, Ffffffffffly ***@***.***> wrote: > 杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。 > > 按照您提出的方法，首先开展排驱和渗吸模拟，获得一定饱和度下的初始相分布，再施加体积力和周期边界，待稳定后根据两相达西计算相渗。（博士论文P116：“After > the permeability of both phases converged, the relative permeability for a > desired...

JBN 是实际生成种经常使用的方法，好处就是省时省力，但是我们认为稳态更准确更能反应由于是油气运移初期的正确分布等等。确定就是测量非常费时费力。可根据需要选择 On Wed, 26 Jun 2024 at 15:49, Ffffffffffly ***@***.***> wrote: > > 判断稳定这个事确实有一定的主观性，在我处理的样本里面，均质性比较好的砂岩，基本上跑一段时间你如果画出湿润相和非湿润相的速度随时间变化曲线就能看出来基本变化很小了，就是稳定了，但是也有些尺度跨度很大的碳酸盐，跑了很久还是能看到速度曲线一直波动很大，如果能看出来事周期性变化，那也可以认为事稳定了，如果不是，可能就还得继续让流动继续一段时间，或者可视化下油水分布，看看是不是分布合理，人工判断下。这个和样本的复杂度是高度相关的 > … > On Wed, 26 Jun 2024 at 15:19, Ffffffffffly *@*.*> wrote: > 杨老师，您好。我最近在做一些LBM计算稳态相渗的工作，阅读了您的博士论文，有些困惑想要请教您。...

https://core.ac.uk/download/pdf/76988721.pdf On Wed, 3 Jul 2024 at 13:58, Sam Ultraman ***@***.***> wrote: > JBN > 是实际生成种经常使用的方法，好处就是省时省力，但是我们认为稳态更准确更能反应由于是油气运移初期的正确分布等等。确定就是测量非常费时费力。可根据需要选择 > … > On Wed, 26 Jun 2024 at 15:49, Ffffffffffly *@*.*> wrote: >...

Hi, thanks for your interest, vtr file is naturally supported by Paraview, maybe there is some configuration problems so the simulation was not properly run? What error did you get...

I see, did you click apply button? then you should have the small eye like this: [image: image.png] it means displayed, if there is no eye over there, it means...

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

你好，对于你的问题: 1: 现在默认的设置是毛细力主导的 （capillary dominate）所以流动非常缓慢，毛细力远大于粘性力，所以速度尤其是wetting phase的速度很小，所以不怎么看得出来，这个是正常的 2： 程序默认参数不是生成那个图的参数，那个图的毛细数好像更小，需要大的表面张力外加大的bodyforce 运行很长时间才可以得到，毛细数小的时候其实对应的dt 也很小，另外真实实验钟毛细数小的实验本身物理时间也很长（界面移动也很缓慢），所以模拟需要的时间更长了。具体参数需要根据你的需求来定，如果你觉得移动的慢，可以加大bodyforce, 或者减小表面张力CapA, 另外粘度也有一定影响 3： 这个代码就是这篇文章的复现： https://www.sciencedirect.com/science/article/abs/pii/S0309170808000511 可以参考这篇文章就知道这里面参数的意义了。另外i毛细数是一个无量纲数，你可以Google 一下capillary number 就知道怎么计算了，比较简单 :-) On Sun, 8 Oct 2023 at 02:50, renxiaosa00 ***@***.***> wrote:...

如果你多试几种参数就可以找到了，数值模拟和实验也类似，根据实验结果调整再试。比如你举得不容易推进去，就降低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/ > abs/pii/S0309170808000511 > > 可以参考这篇文章就知道这里面参数的意义了。另外i毛细数是一个无量纲数，你可以Google一下毛细管数就怎么计算了，比较简单就知道了:-) > ……...

数值稳定就是程序会正确运行，不会出现非物理的极大值或者极小值，你可以通过打印整个速度或者压力的极大极小值来看，比如如果计算发散了，很容易出现1e10这类特别大的速度，这就说明计算发散了，需要调整输入的参数 renxiaosa00 ***@***.***> 于 2023年10月12日周四 上午2:47写道： > 如果你多试几种参数就可以找到了，数值模拟和实验也类似，根据实验结果调整再试。比如你举得不容易推进去，就降低CapA > 表面张力，或者加大bodyforce, 或者加大接触角COS值， 如果数值不稳定了，就稍微增大粘度，减小bodyforce, 减小CapA, > 都可以帮助提高数值稳定性 > … > On Mon, 9 Oct 2023 at 10:56, renxiaosa00 *@*.*> wrote: 好，对于你的问题： >...