后处理问题请教

[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 https://github.com/yjhp1016/taichi_LBM3D/issues/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 https://github.com/yjhp1016/taichi_LBM3D/issues/33#issuecomment-2383257896, 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 https://github.com/yjhp1016/taichi_LBM3D/issues/33#issuecomment-2383593133, or unsubscribe https://github.com/notifications/unsubscribe-auth/AJEDKQFJYSIJUTJEZZC42N3ZZFYKZAVCNFSM6AAAAABPBV6Y4GVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDGOBTGU4TGMJTGM . You are receiving this because you commented.Message ID: @.***>