关于OpenFOAM-9中limitTemperature的用法

[miying]

大家好，在使用OF-9中limitTemperature来限制温度，查找了tutorial中的案例，添加了fvConstrain字典文件，其主要使用代码如下：

limitT
{
    type            limitTemperature;
    selectionMode   all;
    min             200;
    max             4000;
}

但是计算过程会出现如下警告（不影响计算）：
--> FOAM Warning :
From function void Foam::fvConstraints::checkApplied() const
in file cfdTools/general/fvConstraints/fvConstraints.C at line 124
Constraint limitTemperature defined for field ea but never used

而且最终这个字典文件也没有起作用，温度还是出现了低于200K的值，想请问各位老师有遇到过类似问题的嘛，是如何解决的？

[李东岳]

你的ea是什么？

[miying]

@李东岳 感谢李老师回复。其实我也没太搞懂ea是什么，可能是湍流熵，看log文件与之对应的还有eau，按照求解器的命名规律应该是层流熵。但是有个问题是我明明是给T加的约束，但是却表现在了ea上，我也给压力p施加了类似边界，log文件对应显示表现到了p上。log文件具体代码如下：

Post-updating compressibleSystem:
DICPBiCGStab:  Solving for Ux, Initial residual = 0.000396322, Final residual = 1.53337e-14, No Iterations 1
DICPBiCGStab:  Solving for Uy, Initial residual = 0.999946, Final residual = 2.81432e-11, No Iterations 1
DICPBiCGStab:  Solving for ea, Initial residual = 7.24313e-05, Final residual = 1.3332e-14, No Iterations 1
DICPBiCGStab:  Solving for eau, Initial residual = 7.24313e-05, Final residual = 1.3332e-14, No Iterations 1
Creating fvConstraints from "system/fvConstraints"

Selecting finite volume constraint type limitPressure
    Name: limitp
    min 10129.8
    max 1.01298e+06

Selecting finite volume constraint type limitTemperature
    Name: limitT
    - selecting all cells
    - selected 75978 cell(s) with volume 0.0408519
DILUPBiCGStab:  Solving for omega, Initial residual = 2.97194e-06, Final residual = 1.22878e-16, No Iterations 1
DILUPBiCGStab:  Solving for k, Initial residual = 1, Final residual = 4.06328e-11, No Iterations 1

但是这个好像是求解器的问题，因为在/opt/openfoam9/tutorials/compressible/rhoSimpleFoam/aerofoilNACA0012算例也使用了类似limit限制，但是其发挥作用是在开始计算之前，之后只要不超出边界，该限制就不会记录在log中，而blastXiFoam却在每一步都运行且记录。可能是源代码的顺序问题，我跟开发者沟通一下。

[李东岳]

类似这种确实可以跟blastFoam那面沟通一下，没准是bug

[miying]

@李东岳 嗯嗯，现在得到的回复是，blastXiFoam 现在还不支持使用fvConstraint和fvModel.

[李东岳]

@miying