CFD中文网

@知行 是的，好像并没有这个包。我是在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)

结果却是

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