后处理问题请教

[sph-student]

老师您好
对于一名预期明年必要且刚接触到LBM模拟的学生而言，您的代码简洁且实用性强，这里想请教一下您的后处理使用的是什么软件？
是否是paraview？

[yjhp1016]

对，这个库输出是vtu格式的，paraview原生支持的，可以直接读取可视化 sph-student ***@***.***> 于 2024年9月29日周日 下午1:50写道：

…

老师您好 对于一名预期明年必要且刚接触到LBM模拟的学生而言，您的代码简洁且实用性强，这里想请教一下您的后处理使用的是什么软件？ 是否是paraview？ — Reply to this email directly, view it on GitHub <#33>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AJEDKQFDHFN6MDPDMB5IA43ZY7ZSVAVCNFSM6AAAAABPBV6Y4GVHI2DSMVQWIX3LMV43ASLTON2WKOZSGU2TIOJRHEYTEOI> . You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>

[sph-student]

谢谢您的解惑，我拟开展礁灰岩（碳酸盐岩，孔隙率56%左右）的水气驱替模拟，目前已经将数字岩心模型转化为txt格式和dat格式，并尝试了几次计算，但是效果并不理想。这里想要请教老师您关于无量纲参数计算的问题。
物理参数：
水的动力粘度（20°C）：μ_water ≈ 1.002 × 10⁻³ Pa·s
空气的动力粘度（20°C）：μ_air ≈ 1.81 × 10⁻⁵ Pa·s
水的密度：ρ_water ≈ 998 kg/m³
空气的密度：ρ_air ≈ 1.204 kg/m³
模型参数：
网格大小：232 × 232 × 232
实际物理尺寸：0.6 × 0.6 × 0.6 cm³
单个网格尺寸（Δx）：0.6 cm / 232 ≈ 2.5862 × 10⁻⁵ m
---------由于是刚接触LBM算法，后续的计算是否是--------
时间步长估算：
假设马赫数 Ma = 0.1，声速 c_s ≈ 343 m/s（其中马赫数您建议取什么）
Δt = Δx / (√3 * 343) ≈ 4.3558 × 10⁻⁸ s
对于水：
τ_water = 3 * (1.002 × 10⁻³ / 998) / ((2.5862 × 10⁻⁵)² / 4.3558 × 10⁻⁸) + 0.5 ≈ 0.5224
niu_l = (0.5224 - 0.5) / 3 ≈ 0.00746
对于空气：
τ_air = 3 * (1.81 × 10⁻⁵ / 1.204) / ((2.5862 × 10⁻⁵)² / 4.3558 × 10⁻⁸) + 0.5 ≈ 0.5004
niu_g = (0.5004 - 0.5) / 3 ≈ 0.000135
通过模拟我想要得到类似img文件夹下ket_drain.gif计算结果
但是我的模拟结果并不是很好

[yjhp1016]

一般的LBM的密度比和粘度比需要在一附近，要不然容易引起数值不稳定的问题，在低雷诺数下，惯性作用可以忽略，所以密度比可以近似为一，在速度很慢的毛细作用主导的渗流中，粘性力远小于毛细里，可以近似用1左右的粘度比。所以取决于你的问题的雷诺数和毛系数

…

On Mon, 30 Sept 2024 at 14:51, sph-student ***@***.***> wrote: 谢谢您的解惑，我拟开展礁灰岩（碳酸盐岩，孔隙率56%左右）的水气驱替模拟，目前已经将数字岩心模型转化为txt格式和dat格式，并尝试了几次计算，但是效果并不理想。这里想要请教老师您关于无量纲参数计算的问题。 *物理参数：* 水的动力粘度（20°C）：μ_water ≈ 1.002 × 10⁻³ Pa·s 空气的动力粘度（20°C）：μ_air ≈ 1.81 × 10⁻⁵ Pa·s 水的密度：ρ_water ≈ 998 kg/m³ 空气的密度：ρ_air ≈ 1.204 kg/m³ *模型参数：* 网格大小：232 × 232 × 232 实际物理尺寸：0.6 × 0.6 × 0.6 cm³ 单个网格尺寸（Δx）：0.6 cm / 232 ≈ 2.5862 × 10⁻⁵ m ---------由于是刚接触LBM算法，后续的计算是否是-------- *时间步长估算*： 假设马赫数 Ma = 0.1，声速 c_s ≈ 343 m/s（其中马赫数您建议取什么） Δt = Δx / (√3 * 343) ≈ 4.3558 × 10⁻⁸ s 对于水： τ_water = 3 * (1.002 × 10⁻³ / 998) / ((2.5862 × 10⁻⁵)² / 4.3558 × 10⁻⁸) + 0.5 ≈ 0.5224 niu_l = (0.5224 - 0.5) / 3 ≈ 0.00746 对于空气： τ_air = 3 * (1.81 × 10⁻⁵ / 1.204) / ((2.5862 × 10⁻⁵)² / 4.3558 × 10⁻⁸) + 0.5 ≈ 0.5004 niu_g = (0.5004 - 0.5) / 3 ≈ 0.000135 通过模拟我想要得到类似img文件夹下ket_drain.gif计算结果 但是我的模拟结果并不是很好 — Reply to this email directly, view it on GitHub <#33 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AJEDKQGH4SRGKBLOABLVVCDZZFJPFAVCNFSM6AAAAABPBV6Y4GVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDGOBTGI2TOOBZGY> . You are receiving this because you commented.Message ID: ***@***.***>

[sph-student]

收到 谢谢老师解惑，根据您的意思，我是否需要主要考虑 fx capa 以及压力差等参数的修改

[sph-student]

感谢老师！ 您的代码 最后的模拟结果惊艳到我们小组所有人了 虽然我们还需要进一步验证模拟的准确性，但是实在是太好看了

[yjhp1016]

很高兴能帮到你们，请注意两项流对网格resolution要求比较高，最细的孔道最好有七八个格子以上。 sph-student ***@***.***> 于 2024年9月30日周一 下午4:58写道：

…

rho1.png (view on web) <https://github.com/user-attachments/assets/a5082fe6-512c-4781-9e97-a3dada58e575> 感谢老师！ 您的代码 最后的模拟结果惊艳到我们小组所有人了 虽然我们还需要进一步验证模拟的准确性，但是实在是太好看了 — Reply to this email directly, view it on GitHub <#33 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AJEDKQFJYSIJUTJEZZC42N3ZZFYKZAVCNFSM6AAAAABPBV6Y4GVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDGOBTGU4TGMJTGM> . You are receiving this because you commented.Message ID: ***@***.***>