怎么把一个标量写成一个标量场

[win]

OpenFOAM怎么可以把一个标量写成一个标量场啊？

[李东岳]

volScalarField tauc_
(
    IOobject
    (
        "tauc",
        mesh_.time().timeName(),
        mesh_,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh_,
    dimensionedScalar(dimTime, 1.0)
),

[win]

谢谢教授

[Exthan]

@李东岳 老师，请问如果一个场量要与一个标量比较大小，该如何比较？

if (this->nut() / this->transport_.nu() <  scalar(0.2))

[cresendo]

@exthan 你这样将一个场除以一个场，得到的是一个新的场，是没有办法和一个标量比较的

[Exthan]

@cresendo 是的，如何让场量里的对应值相除呢，就是取场量每个位置的值进行相除运算呢？

 forAll(Prt_,celli){
    if (this->nut()(celli)/ this->transport_.nu()(celli) <  scalar(0.2))
    ...

错误信息是个temp
.
应该如何修改 才能进行比较？

[Exthan]

@李东岳 老师，求指导

[Exthan]

@exthan

[cresendo]

@exthan 尝试一下

this->nut()()[facei]

[Exthan]

@cresendo 是的，这部分修改已经完成了。将湍流普朗特数从标量改为场量之后，湍流模型调用好像并没有用这个新的湍流普朗特数？请问还需要修改其他地方吗？

[cresendo]

@exthan 所以之前就是这边语法的问题，对吗？关于湍流普朗特数，本人不是太了解，可能没法给出有效的建议

[Exthan]

@cresendo 是的，之前调用值不正确？

[李东岳]

@exthan 在 怎么把一个标量写成一个标量场 中说：

将湍流普朗特数从标量改为场量之后，湍流模型调用好像并没有用这个新的湍流普朗特数？请问还需要修改其他地方吗？

怎么看出来的

[Exthan]

@李东岳 输出的最大最小Prt均为同一个值，并没有变化。存在情况1：新prt模型失效；2：流场刚好都在一个判断条件下。所以如何将Prt以一个文件形式输出？

在solver里面加入可以输出吗？

Prt = tranport.Prt();