CFD中文网

落

@capillaryFix 是的，我求解了Eqn.30 -> Eqn.34 ->Eqn.33，但是Eqn.33的结果用在哪里了？

落

@capillaryFix ，您好。我仔细研究了您的文章，有两点疑惑，向您请教：

最后一步是关于capillary flux的过滤是基于原始的$F_{st,f}$，过滤结束的$\phi^{filt}_{cf}$是放在这里吗？

solve
        (
            UEqn
         ==
            fvc::reconstruct
            (
                (
                   - ghf*fvc::snGrad(rho)
                    - fvc::snGrad(p_rgh)
                ) * mesh.magSf()+phi_filt_cf
            ),
            mesh.solver(U.select(finalIter))
        );

和

surfaceScalarField phig
    (
        (
          - ghf*fvc::snGrad(rho)
        )*rAUf*mesh.magSf()+ phi_filt_cf*rAUf
    );

这步$F^{f}{st,f} = F{st,f} − F^{filt}_{st,f}$不知道怎么植入进方程，是不是关于capillary flux的过滤应该基于这步过滤的结果

谢谢、谢谢

落

@知行是的，好像并没有这个包。我是在paraview绿色版里找到python的相关包，直接复制过去的

落

所以，这个是和版本相关的，我2022R2 遇到相同的error

还有好多warning

落

@落花风放两张根据@田畔的风大佬程序出的图，真的超棒。

落

@田畔的风感谢，是我看的不仔细了。还有两点疑问：

     data = paraview.servermanager.Fetch(casefoam)
     data = dsa.WrapDataObject(data)

这个data变量并没有被用到，注释掉也没有影响(我是这样)，data起什么作用呢？

Path("slice").mkdir(parents=True, exist_ok=True)

建这个空文件夹是？

落

@田畔的风膜大佬！研究了半晚上。几点小疑问，向您请教：

reader = OpenFOAMReader(FileName="PEE.foam", CaseType=0)

caseType=0 是控制什么的？这句导致我的case读不到任何一个时间步。删除之后可以读到。

logData = np.where(
    mesh.point_data['driftDensity'] < 1e-12,
    -12.0,
    np.log10(mesh.point_data['driftDensity']),
)

这句是控制什么的呢？我print(mesh.dict),发现没有‘driftdensity’
我用的是六面体网格，怎么划分成三角形，你的代码里有相关的吗？（我刚理解到第二个问题，后面还没看，肝不动了）

多谢分享，多谢多谢

落

@chszkc 感谢分享，在尝试使用这个包的过程中遇到点小麻烦，向您请教。
首先，我的计算域是这样的：

我尝试画出温度等值线图，代码如下

#%%
from fluidfoam import readmesh
sol = './T293/massFlowRate0.412'
x, y, z = readmesh(sol)
from fluidfoam import readscalar

timename = 'latestTime'
T = readscalar(sol, timename, 'T')
# %%
import numpy as np
import matplotlib.pyplot as plt

# Define plot parameters
fig, ax = plt.subplots(figsize=(7, 3), dpi=100)
plt.rcParams.update({'font.size': 10})
plt.xlabel('x')
plt.ylabel('y')
# Plots the contour of sediment concentration
levels = np.arange(30, 300, 10)

plt.tricontourf(x, y, T, cmap=plt.cm.Reds, levels=levels)

结果却是

请问您有类似经历吗
先谢谢您

CFD中文网

落花风

帖子