CFD中文网

@香柏树 你这个问题解决了吗

已经找到解答：https://cfd-china.com/topic/1553/并行-获取某个量达到最大值的这个网格中的其它值?_=1698117342935

在controlDict中我用fieldMinMax求最值，如果单核计算得到的最大温度是正确的，但是并行得到的是某个块上的，。如果想从多个块中筛选出最大值，应该怎么做，我网上看说用reduce，Pstream::，哪位老师有碰到过这种问题？

@氕氘-川 已解决，利用cantera在把cti转为xml

请教一下大家，cantera可以把chemkin机理文件转为cti，yaml文件，那.xml的燃烧机理文件怎么获得,或者怎么转啊？
提前感谢大佬的回复。

@李东岳老师，这个是求某个标量场内最值的吧？ 我是想先求出最大温度处的位置，然后根据这个位置在求出该位置处的速度

openfoam怎么获得最大温度位置处的速度或者密度？

openfoam或者paraview能不能监测两个等值线交点上值随时间变化啊？就比如下面这个图是两个contour线，先找到交点位置，然后在看这个交点位置的其他值。

@李东岳 效果非常好，非常感谢李老师

请教一下大家，openfoam比如要算这种： 在openfoam内部计算可以实现吗？
如果用python，是不是还是要依靠paraview脚本啊？

@氕氘-川 批量保存的脚本老报错： ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
vtkCSVWriter (000001F6DAFC08C0): Unable to open file: /-0.1150.csv
ERROR: In C:\bbd\df0abce0\build\paraview\src\paraview\ParaViewCore\VTKExtensions\Default\vtkCSVWriter.cxx, line 82
vtkCSVWriter (000001F6DAFCC8C0): Unable to open file: /-0.0650.csv
...

这是我的脚本：
import numpy as np
from paraview import numpy_support as ns
import paraview as pv
import paraview.simple as ps

afoam = FindSource('a.foam')

for y in np.arange(-0.115,0.115,0.05):

# Make a cut of the wing to have a local 2D profile
Slice1 = Slice( Input=afoam )
Slice1.SliceType = "Plane"
Slice1.Triangulatetheslice = 0
Slice1.SliceOffsetValues = [0.0]
Slice1.SliceType.Origin = [0.0, y, 0.0]
Slice1.SliceType.Normal = [0.0, 1.0, 0.0]
descriptiveStatistics1 = DescriptiveStatistics(Slice1)
descriptiveStatistics1.VariablesofInterest = ['T']
descriptiveStatistics1.AttributeMode = 0
Render()
writer = CreateWriter(str(y)+".csv",OutputPort(descriptiveStatistics1, 1))
writer.FieldAssociation = "Cells"
writer.UpdatePipeline()
del writer
#SaveData(str(y)+".csv", proxy=OutputPort(descriptiveStatistics1, 1), FieldAssociation='Cells')
descriptiveStatistics1.UpdatePipeline()

大家好，我切了很多slice求每个切片上平均值，现在求出来后，因为只需要mean这个值，我怎样可以批量导出这个数值？

这是我保存的方式：
Render()
writer = CreateWriter(str(y)+".csv"descriptiveStatistics1)
writer.FieldAssociation = "Points"
writer.UpdatePipeline()
del writer
但是这样保存会说writer.FieldAssociation = "Points"里面的points不对，然后我去了之后就会报错。

@xpqiu 太感谢您了 您教的好细， 非常感谢 ，对我帮助非常大 ，谢谢您

问一下大家，我下面这个图，如何把流向（x方向）的速度叠加起来成一条线的数据，再做平均啊？有没有处理它的脚本或者paraview处理方法啊？

@trojason 你看一下这个https://zhuanlan.zhihu.com/p/32441569

@李东岳 两种方法，一种是用胖节点可以一次性画好，另一种是先画差网格（千万级别），之后多次refine就可以。我的总共用时3个小时，总共网格是2.6亿

@bestucan 好好好 谢谢您

@xpqiu 感谢您的回复。
我分的processors太多 转的时候是不是也需要很久的时间啊？

@李东岳 别别别，不是大佬 ，李老师这说的我都无地自容了，我用的超算。现在可以了，每次只导入一个标量就不奔溃了