CFD中文网

把原来流场的小眼睛点开

在openfoam7中，涉及燃烧的求解器把原先默认输出的Qdot设置成了默认不输出，并且增加了一个后处理函数Qdot，在有些情况下希望计算后后处理输出Qdot的值。

但是使用reactingFoam -postProcess -func Qdot后发现出现了Qdot的场，但是全是零

假如修改代码，在求解器中建立Qdot场并注册，并且在每个时间步更新，就像openfoam6那样，使用reactingFoam -postProcess -func "writeObjects(Qdot)"后也是发现出现了Qdot的场，但是全是零

怀疑是Qdot计算时是integrateReactionRate，在constant/combustionProperties中设置integrateReactionRate off;也没用

但是假如在controlDict文件中加入#includeFunc Qdot进行运行时后处理就没问题，可以输出正常的Qdot

所以有什么办法在计算后后处理得到Qdot的吗，原来openfoam6使用writeObjects是可以的

可能是windows和linux文件换行符的问题

@bestucan 有道理，但是无论是哪方面问题，都没有解决办法
因为平台是购买的机时，无法更换硬件
计算域十分简单就是个方块，也不太好优化分块方案了

超算平台：
两个节点，每个节点20核，共40核，centos系统，Intel编译，Intel-mpi，磁盘没有特别的设置
工作站：
共48核，使用40核，ubutnu，gcc编译，openmpi

计算相同的算例，一维爆轰管问题，计算效率如下

1. cpu时间和真实时间差了这么多，是因为通信还是磁盘读写的原因
2. 大家跑算例的时候也是相差这个量级吗
3. 有没有什么优化的方法

regions
(
    cylinderToCell
    {
        p1 (0.102 0 -1);
        p2 (0.102 0 1);
        radius 0.1;
        fieldValues
        (
            volScalarFieldValue T 300
            volScalarFieldValue p 1e5
        );
    }
);

以前用的

fixedValue可以设定定常非均匀场，参考一下codestream的资料

两个计算域从5000网格映射到300000网格，网格不相容，重叠的只是一小部分，选的mapNearest方法，运行特别慢，十多分钟才跑完。
感觉不应该这么慢吧，哪里出问题了？

@刘雄国 对，我是说没有功能上的作用

应该没啥用，单纯为了保证只引用一次此文件，类似program once

@Foamer24 抱歉，是divSchemes

@Foamer24 应该不是这个default吧，后面主要物理量都覆盖了

请问插值格式用的什么

@东岳 老师，对于fluent中的压力出口边界条件所说的extrapolated from the flow in the interior所说的extrapolated是这种外推边界条件吗，他是怎么实现的呢？还是说他只是zerogrident。

fluent中的压力出口边界条件，资料上写着

压力出口边界条件在流场出口边界上定义静压，而静压的值仅在流场为亚音速时使用。如果出口边界上流场达到超音速，那么边界上的压力将从流场内部通过插值得到。
其他流场变量均从流场内部插值获得。

这里说的插值是怎样的插值，导数为零？还是二阶导数为零？

@东岳 确实是，还有一种思路就是使用再向内一层单元的梯度作为边界梯度，可是好像openfoam的为一般非结构网格设计的网格结构得不到这个信息

在openfoam中有没有和fluent压力出口对应的边界条件，压力温度速度分别对应什么

在openfoam中，有没有一种边界条件，可以实现由流场内部向边界外推，保证边界二阶导数为零。
或者如何自己实现，关键是不知道如何获取物理量在边界体心上的梯度（当边界值还未给出）。