A Small Question #16
Comments
|
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 when you try to visualize use vtr?
If you want to export some specific data, you can go to the last part of
the code, you will find numpy format data:
pointData={ "Solid": np.ascontiguousarray(solid.to_numpy()),
"rho": np.ascontiguousarray(rho.to_numpy()[0:nx,
0:ny,0:nz]),
"phase": np.ascontiguousarray(psi.to_numpy()[0:
nx,0:ny,0:nz]),
"velocity": (np.ascontiguousarray(v.to_numpy()[0
:nx,0:ny,0:nz,0]), np.ascontiguousarray(v.to_numpy()[0:nx,0:ny,0:nz,1]),np.
ascontiguousarray(v.to_numpy()[0:nx,0:ny,0:nz,2]))
}
…On Fri, 7 Jul 2023 at 09:27, ViviviVirgil ***@***.***> wrote:
Dear author, hello. I am very interested in your developed work, but I
encountered some issues during the trial. 1. After running the
"lbm_solver_3d_2phase.py" file in the "2 phase" module, I generated a vtr
file, but I am unable to open and read it in ParaView. Could you please
advise me on how to view the simulation results? 2. Additionally, I would
like to output specific frames of the model for further research. Do you
have any suggestions?
—
Reply to this email directly, view it on GitHub
<#16>, or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQFR3YC6PJEQUG73IGTXO7CAJANCNFSM6AAAAAA2BPXDAM>
.
You are receiving this because you are subscribed to this thread.Message
ID: ***@***.***>
|
|
Thank you for your response. I only made modifications to the iteration count in "lbm_solver_3d_2phase.py" in the "2phase" directory, setting it to 8001 iterations, and also changed the number of iterations for exporting intermediate results to 1000. Here is my environment configuration. However, when visualizing the results in ParaView, no images are displayed. |
|
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 the geometry is
not displayed
…On Sat, 8 Jul 2023 at 10:48, ViviviVirgil ***@***.***> wrote:
Thank you for your response. I only made modifications to the iteration
count in "lbm_solver_3d_2phase.py" in the "2phase" directory, setting it to
8001 iterations, and also changed the number of iterations for exporting
intermediate results to 1000. Here is my environment configuration.
However, when visualizing the results in ParaView, no images are displayed.
[image: 654e50680354015463ae4f0ab67c694]
<https://user-images.githubusercontent.com/33839347/251965603-9af9d52d-7bfe-4e57-b8a2-9b28e3ee11f0.png>
[image: 83603465e12a0879cd9ce580999d4c5]
<https://user-images.githubusercontent.com/33839347/251965628-72b5edb1-c826-4d88-9d71-76f36f8bdeae.png>
[image: de5a06994ef8cccd5e0eb4c35c47e20]
<https://user-images.githubusercontent.com/33839347/251965635-16ed113a-bdb8-4187-8512-c31fa3e68282.png>
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQF3XKUOVM5SDNVPGEDXPEUIPANCNFSM6AAAAAA2BPXDAM>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
Okay......，thanks!I have iterated the program 80000 times according to the instructions above, and imported the resulting 8 vtr files into Paraview. However, I did not achieve the desired fluid flow effect. I would like to know if this is a configuration or Paraview software operation issue, or if I need to further study the LBM principles and modify the program. |
|
|
你可以多探索paraview里的功能，时间比较长了，我记得大概是设置的问题
…---- 回复的原邮件 ----
| 发件人 | ***@***.***> |
| 日期 | 2023年10月07日 13:52 |
| 收件人 | ***@***.***> |
| 抄送至 | ***@***.***>***@***.***> |
| 主题 | Re: [yjhp1016/taichi_LBM3D] A Small Question (Issue #16) |
好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。
您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题
—
Reply to this email directly, view it on GitHub, or unsubscribe.
You are receiving this because you authored the thread.Message ID: ***@***.***>
|
这个是不是迭代步数的问题呀，目前只有8个结果文件，流动起来也只流动1/6处，是不是设置多一些迭代步数，得到更多的结果，才会形成文章中这样从头流动到尾部的结果呀？ |
|
请问你说的想要的效果是视觉上的效果还是流体分布的位置和example图上不一样？
如果是视觉效果，那是paraview使用上面的问题，vtr 文件只提供可视化的数据，怎么渲染取决于如何操作paraview
（例子图片钟可视化的是固体部分（使用了半透明的设置）和non-wetting phase 的分布，不是速度或者压力）
如果是这些流体分布位置不一样，那是模拟初始化数据的设置问题，可以通过变化压力差和接触角还有表面张力大小调整，如果你的模拟毛细数很小，比如<1e-5（意味着毛细力主导，超过粘性作用很多）
那么模拟会非常的慢，80000步还远远不够， 如果想更快的看到变化，可以加大毛细数，让它在1e-4 到 1e-3这个数量级，模拟会快很多
…On Sat, 7 Oct 2023 at 06:52, renxiaosa00 ***@***.***> wrote:
好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。
[image: 图像]
<https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png> [image:
图像]
<https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png>
您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
您好，非常感谢您的回信，我是刚开始接触LBM，Python倒是会用，关于您的回复，我有一些问题： |
|
你好，对于你的问题:
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:
请问你说的想要的效果是视觉上的效果还是流体分布的位置和example图上不一样？ 如果是视觉效果，那是paraview使用上面的问题，vtr
文件只提供可视化的数据，怎么渲染取决于如何操作paraview （例子图片钟可视化的是固体部分（使用了半透明的设置）和non-wetting
phase 的分布，不是速度或者压力）
如果是这些流体分布位置不一样，那是模拟初始化数据的设置问题，可以通过变化压力差和接触角还有表面张力大小调整，如果你的模拟毛细数很小，比如<1e-5（意味着毛细力主导，超过粘性作用很多）
那么模拟会非常的慢，80000步还远远不够， 如果想更快的看到变化，可以加大毛细数，让它在1e-4 到 1e-3这个数量级，模拟会快很多
… <#m_3962391688291200852_>
On Sat, 7 Oct 2023 at 06:52, renxiaosa00 *@*.*> wrote:
好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。
[image: 图像]
https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png
<https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png>
[image: 图像]
https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png
<https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png>
您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题 — Reply to this email directly, view it on
GitHub <#16 (comment)
<#16 (comment)>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE
<https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE>
. You are receiving this because you commented.Message ID: @.*>
您好，非常感谢您的回信，我是刚开始接触LBM，Python倒是会用，关于您的回复，我有一些问题：
（1）视觉上的效果目前是红蓝颜色的展示，就是按照原来的默认参数运行的结果，流体从入口到出口流动过程非常少，只有整个影像的1/6
（2）流体分布的位置和example图上不一样，我将原来的80000改成800000后，得到81个结果，观察流体的流动，从入口到迭代步骤停止时，比之前的8个结果图要流动的更长远一些。
（3）初始参数的设置上，压力差、接触角、表面张力大小这些是实验获取的还是依据什么参数决定的呢？模拟毛细数这个参数是怎么根据数据进行统计的呢？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQCABINAY3CYMLX6MVLX6IBG3AVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRHA4TAMZVGU>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
您好，请问要实现文章中这样的两相流数值情况，相关的参数是什么样子的呢？因为我看代码中有的参数是文章中已经提到的参数，别的参数都是什么样的呢？如果方便的话，我的微信是rxs137,欢迎您的指导，非常感谢您！ |
|
如果你多试几种参数就可以找到了，数值模拟和实验也类似，根据实验结果调整再试。比如你举得不容易推进去，就降低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
<https://www.sciencedirect.com/science/article/abs/pii/S0309170808000511>
可以参考这篇文章就知道这里面参数的意义了。另外i毛细数是一个无量纲数，你可以Google一下毛细管数就怎么计算了，比较简单就知道了:-)
…… <#m_-7977376210493536912_>
On Sun, 8 Oct 2023 at 02:50, renxiaosa00 *@*.*> wrote:
请问你说的想要的效果是视觉上的效果还是流体分布的位置和example图上不一样？ 如果是视觉效果，那是paraview使用上面的问题，vtr
文件只提供可视化的数据，怎么渲染取决于如何操作paraview （例子图片钟可视化的是固体部分（使用了半透明的设置）和non-wetting
phase 的分布，不是速度或者压力）
如果是这些流体分布位置不一样，那是模拟初始化数据的设置问题，可以通过变化压力差和接触角还有表面张力大小调整，如果你的模拟毛细数很小，比如<1e-5（意味着毛细力主导，超过粘性作用很多）
那么模拟会非常的慢，80000步还远远不够， 如果想更快的看到变化，可以加大毛细数，让它在1e-4 到 1e-3这个数量级，模拟会快很多 …
<#m_3962391688291200852_> On Sat, 7 Oct 2023 at 06:52, renxiaosa00 @.>
wrote:
好的……，谢谢！我按照上面的说明迭代了程序80000次，并将生成的8个vtr文件导入到Paraview中。然而我并没有达到想要的流体流动效果。我想知道这是配置问题还是Paraview软件操作问题，还是需要进一步研究LBM原理并修改程序。
[image: 图像]
https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png
<https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png>
https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png
<https://user-images.githubusercontent.com/33839347/252137951-91c38163-dc1a-4d2f-89de-4614c641ef64.png>
[image: 图像]
https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png
<https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png>
https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png
<https://user-images.githubusercontent.com/33839347/252138068-d5f98fa7-0b6e-484f-92d5-6ccbfddc21a0.png>
您好，请问您是怎么解决这个问题的呢？我也遇到了相同的问题 — Reply to this email directly, view it on
GitHub <#16 <#16> (comment)
<#16 (comment)
<#16 (comment)>>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE
<https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE>
https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE
<https://github.com/notifications/unsubscribe-auth/AJEDKQB5H4OMHLEODJK5UOLX6DUYDAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRGYYTCMZWHE>
. You are receiving this because you commented.Message ID: @.>
您好，非常感谢您的回信，我是刚开始接触LBM，Python倒是会用，关于您的回复，我有一些问题：
（1）视觉上的效果目前是红蓝颜色的展示，就是按照原来的默认参数运行的结果，流体从入口到出口流动过程非常少，只有整个影像的1/6
（2）流体分布的位置和example图上不一样，我将原来的80000改成800000后，得到81个结果，观察流体的流动，从入口到迭代步骤停止时，比之前的8个结果图要流动的更长远一些。
（3）初始参数的设置上，压力差、接触角、表面张力大小这些是实验获取的还是依据什么参数决定的呢？模拟毛细数这个参数是怎么根据数据进行统计的呢？ —
Reply to this email directly, view it on GitHub <#16 (comment)
<#16 (comment)>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQCABINAY3CYMLX6MVLX6IBG3AVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRHA4TAMZVGU
<https://github.com/notifications/unsubscribe-auth/AJEDKQCABINAY3CYMLX6MVLX6IBG3AVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJRHA4TAMZVGU>
. You are receiving this because you commented.Message ID: @.***>
您好，请问要实现文章中这样的两相流数值情况，相关的参数是什么样子的呢？因为我看代码中有的参数是文章中已经提到的参数，别的参数都是什么样的呢？如果方便的话，我的微信是rxs137,欢迎您的指导，非常感谢您！
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQD6K4PA4X26HCQHLKLX6PC6FAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONJSG4YDENZQGQ>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
您好，这个数值稳定是怎么看的呢？是看程序运行过程中打印在控制台的参数列表嘛？ |
|
您好，（1）速度对应的是velocity，相对应的是phase，但是压力对应的是哪个参数哦？ |
|
计算全部流体element的，如果三维就是整个geometry注意排除固体部分，另外你是要算二维切片中的流动？那你得把它弄成伪三维，就是把这个一样的切片叠五层做成一个比如说150
X 150 X 5的结构计算
renxiaosa00 ***@***.***> 于 2023年10月18日周三 上午3:32写道：
… 您好，（1）速度对应的是velocity，相对应的是phase，但是压力对应的是哪个参数哦？
（2）速度最大值还有压力的最大最小这三个值，是指的模拟的某一刻，统计各自的三维立方体内150*150*
150大小的岩石内所有像素还是统计的各自所在当前的二维界面150*150当前切面呀？比如说，模拟到30w步，是计算三个值的每一个变量的三维立方体的值，还是计算每一个值所在的当前迭代步之下的二维切面的值呢？
我看单相流可以通过初始化一个lb3d.get_max_v(),两相流的话，应该就是直接统计numpy了吧？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQHKXEJMN7XZGR24KJDX745UJAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXGUYTOMZVGQ>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
我的数据是三维岩石，150乘以150乘以150的，统计速度最大值还有压力的最大最小这三个值的时候是，整个三维数据的流体一起统计，而不是说流体流到某一时刻到达一个切面，此切面150乘以150所在的二维图片上的流体速度，对吧？ |
|
压力是rho 密度换算来的，LBM 基于低雷诺数的假设，没有解poisson 方程 使用EOS 来计算压力
…On Wed, 18 Oct 2023 at 08:32, renxiaosa00 ***@***.***> wrote:
我的数据是三维岩石，150乘以150乘以150的，统计速度最大值还有压力的最大最小这三个值的时候是，整个三维数据的流体一起统计，而不是说流体流到某一时刻到达一个切面，此切面150乘以150所在的二维图片上的流体速度，对吧？
还有就是压力好像没看到计算值呀？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQDWGU6OBJJEE7FHLXTX76AY7AVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXHAZTSNBTHA>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
rho结果，无论粘度是否一样，也不管capA大小如何，我测试的3组实验中均没有什么变化，是不是我paraview设置的问题呀？还是本身就是变化不大呀？ |
|
如果你用的默认的color map
看到的是黄色的，说明中间有那一步simulation其实已经发散了，出现NaN了，所以你看到的都一样，应该是初始设置有些问题
…On Wed, 18 Oct 2023 at 09:20, renxiaosa00 ***@***.***> wrote:
rho结果，无论粘度是否一样，也不管capA大小如何，我测试的3组实验中均没有什么变化，是不是我paraview设置的问题呀？还是本身就是变化不大呀？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQDDNCX3AUSSI5EDDG3X76GNTAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXHEZTGNZRGI>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
我测试的3组实验基本上从0步到第80个结果这块，一直没变化，一直是下面的这种状态 |
|
因为你把固体流体一起显示了，所以差别看不出来了，因为流体部分的rho 变化大概在在0.99-1.01 ，但是固体是0.0， 所以如果都显示出来看不出区别
你可以先把流体部分抠出来，只显示这部分的rho， rescale 一下color bar 应该能看到变化
另外LBM的压力可能不是很准，没有接poisson 方程
…On Wed, 18 Oct 2023 at 09:28, renxiaosa00 ***@***.***> wrote:
我测试的3组实验基本上从0步到第80个结果这块，一直没变化，一直是下面的这种状态
[image: 0020e1790bf348df43bbd7e5b664aa2]
<https://user-images.githubusercontent.com/36992106/276171711-dfc3adb1-20ea-4f67-b999-1453290c4b41.jpg>
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQCLTCM2YDLYIMC56Y3X76HMPAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXHE2DQOJYHE>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
|
您好，在img_ftb131.txt中岩心数据都是由255和0组成的，请问应该怎样才能把流体部分抠出来？ |
|
你的问题是如何只显示孔隙部分？
RiverNight666 ***@***.***> 于 2024年3月21日周四 上午3:10写道：
… 您好，在img_ftb131.txt中岩心数据都是由255和0组成的，请问应该怎样才能把流体部分抠出来？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQH7JFIQ23ZUGE3HP2TYZJFSNAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDAMJRGEYDSOJVGY>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
老师，您好，我在做气水两相的流动模拟，我根据压力9MPa和温度15℃在NIST Chemistry WebBook上查到甲烷的密度是90.65Kg/m^3,水的密度是999.1Kg/m^3,根据相似性准则，两相流体的单位换算rho(physcial)(CH4)/rho(physcial)(water)=rho(LBM)(CH4)/rho(LBM)(water)=0.09，如果水在格子下的密度是1，那么甲烷的密度是0.09，我看见老师您博士论文中使用的葵烷和盐水的模拟，LBM单位制下密度都是1， |
|
密度比大的需要特殊的LBM, 由于以前读博时候做的是岩石里面的渗流，雷诺数非常小，所以inertia 可以忽略，也就是密度比就不重要，可以近似用1替代
xiwei gao ***@***.***> 于 2024年8月6日周二 上午7:18写道：
… 因为你把固体流体一起显示了，所以差别看不出来了，因为流体部分的rho 变化大概在在0.99-1.01 ，但是固体是0.0，
所以如果都显示出来看不出区别 你可以先把流体部分抠出来，只显示这部分的rho， rescale 一下color bar 应该能看到变化
另外LBM的压力可能不是很准，没有接poisson 方程
… <#m_4044327581349811504_>
On Wed, 18 Oct 2023 at 09:28, renxiaosa00 *@*.*> wrote:
我测试的3组实验基本上从0步到第80个结果这块，一直没变化，一直是下面的这种状态 [image:
0020e1790bf348df43bbd7e5b664aa2]
https://user-images.githubusercontent.com/36992106/276171711-dfc3adb1-20ea-4f67-b999-1453290c4b41.jpg
<https://user-images.githubusercontent.com/36992106/276171711-dfc3adb1-20ea-4f67-b999-1453290c4b41.jpg>
— Reply to this email directly, view it on GitHub <#16 (comment)
<#16 (comment)>>,
or unsubscribe
https://github.com/notifications/unsubscribe-auth/AJEDKQCLTCM2YDLYIMC56Y3X76HMPAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXHE2DQOJYHE
<https://github.com/notifications/unsubscribe-auth/AJEDKQCLTCM2YDLYIMC56Y3X76HMPAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTONRXHE2DQOJYHE>
. You are receiving this because you commented.Message ID: @.*>
老师，您好，我在做气水两相的流动模拟，我根据压力9MPa和温度15℃在NIST Chemistry
WebBook上查到甲烷的密度是90.65Kg/m^3,水的密度是999.1Kg/m^3,根据相似性准则，两相流体的单位换算rho(physcial)(CH4)/rho(physcial)(water)=rho(LBM)(CH4)/rho(LBM)(water)=0.09，如果水在格子下的密度是1，那么甲烷的密度是0.09，我看见老师您博士论文中使用的葵烷和盐水的模拟，LBM单位制下密度都是1，
1.jpg (view on web)
<https://github.com/user-attachments/assets/53d944bc-8cd2-4ebc-bddd-8f2a39a3631c>
想请教老师对于密度比大的我们应该怎么去处理？
—
Reply to this email directly, view it on GitHub
<#16 (comment)>,
or unsubscribe
<https://github.com/notifications/unsubscribe-auth/AJEDKQB4N2GZLAI5EEIJWETZQBTCZAVCNFSM6AAAAAA2BPXDAOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDENZQGQ3TAMRUGY>
.
You are receiving this because you commented.Message ID:
***@***.***>
|
Dear author, hello. I am very interested in your developed work, but I encountered some issues during the trial. 1. After running the "lbm_solver_3d_2phase.py" file in the "2 phase" module, I generated a vtr file, but I am unable to open and read it in ParaView. Could you please advise me on how to view the simulation results? 2. Additionally, I would like to output specific frames of the model for further research. Do you have any suggestions?
The text was updated successfully, but these errors were encountered: