CFD中文网

东岳老师，想参加！hujinhua@mail.bnu.edu.cn

大家好，

当对每个网格的volScalarField p 场 求梯度时，

1. 在OpenFOAM中 fvm::grad(p )，是按照实际的xyz求的，还是网格倾斜的x'y'z'求的呢？
   如下图所示, 当网格为（1）倾斜的非结构化网格（平行四边形）,（2）倾斜的结构化网格

2. 对于（1）这样倾斜的非结构化网格（平行四边形）来说，
   （a）x'y'z'轴是如图这么画的吗 (但这样导致∠x' O y' 不等于 90°)？
   （b）那它的面法向梯度，是朝哪个方向呢？

算例为土柱冻结算例，起始土柱内部温度为6.7℃，土柱顶部施以温度为-6℃的循环液流(the circulating fluid)，使土柱从上到下逐步冻结，下图为示意图，

温度边界条件：
变热通量上边界条件/a variable heat flux upper boundary condition：
qh=hc(TCoolant-TTop), hc=28 W m-2 K-1
其中：TCoolant是循环液流的温度值(-6 ℃)，而TTop是土柱顶部温度，我理解的是每一次时间步长计算后都会更新，所以顶部通量是一个动态更新的过程。
1、之前有看过关于fixed heat flux boundary的帖子，采用fixedGradient边界条件代替，但没找到关于 variably heat flux boundary condition
2、 采用swak4Foam？
希望能给我一些建议，共同探讨下。非常感谢！

@general 在 OpenFoam定义新边界条件问题---通量边界 中说：

您好，您最后这个边界条件成功了吗？是变通量边界条件吗？谢谢！

请问您绝热壁面的类型，设置的是哪个啊

@东岳 第一次发帖求助，就能得到东岳老师的回复，真的非常感谢！
（1）是根据laplacianFOAM稍微改了下的求解器，但是如果是探针probe的话是可以运算出结果的。那cellSource不运行function的话，我需要怎么改下求解器吗
（2）还想请教下，算例边界条件设置的问题，左图是压力，右图是温度。
左图压力上下边界条件均为zero flux，右侧问题上下右边界条件为zero conductive flux。对于zero flux和zero conductive flux，边界条件是不是应设置为symmetryPlane，还是zeroGradient。我计算出的结果差别很大，不知道问题出在哪里，所以想确定下是不是我边界条件设置的不对。

第一次发帖求助，就能得到东岳老师的回复，真的非常感谢！
（1）是根据laplacianFOAM稍微改了下的求解器，但是如果是探针probe的话是可以运算出结果的。那cellSource不运行function的话，我需要怎么改下求解器吗
（2）还想请教下，算例边界条件设置的问题，左图是压力，右图是温度。
左图压力上下边界条件均为zero flux，右侧问题上下右边界条件为zero conductive flux。对于zero flux和zero conductive flux，边界条件是不是应设置为symmetryPlane，还是zeroGradient。我计算出的结果差别很大，不知道问题出在哪里，所以想确定下是不是我边界条件设置的不对。

@东岳 谢谢东岳老师能回复我，开心。下面是你说的时间步log吧。我设的时间步长是240s。

Create time

Create mesh for time = 0

SIMPLE: convergence criteria
    field P	 tolerance 1e-05
    field U	 tolerance 1e-05
    field nuTilda	 tolerance 1e-05
Reading field p
Reading transportProperties
Calculating 
Time = 240
Calculating the Saturation
Calculating the pressure head
DICPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0
Calculating the velocity
Calculating lamT
DILUPBiCG:  Solving for T, Initial residual = 0.550522, Final residual = 8.33185e-07, No Iterations 48
ExecutionTime = 34.74 s  ClockTime = 41 s

Time = 480


Calculating the Saturation


Calculating the pressure head

DICPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0

Calculating the velocity


Calculating lamT

DILUPBiCG:  Solving for T, Initial residual = 0.0427664, Final residual = 8.02042e-07, No Iterations 45
ExecutionTime = 37.47 s  ClockTime = 44 s

Time = 720


Calculating the Saturation


Calculating the pressure head

DICPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0

Calculating the velocity


Calculating lamT

DILUPBiCG:  Solving for T, Initial residual = 0.0180619, Final residual = 8.65532e-07, No Iterations 41
ExecutionTime = 39.39 s  ClockTime = 46 s

Time = 960


Calculating the Saturation


Calculating the pressure head

DICPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0

Calculating the velocity


Calculating lamT

DILUPBiCG:  Solving for T, Initial residual = 0.0128672, Final residual = 9.31796e-07, No Iterations 39
ExecutionTime = 40.97 s  ClockTime = 47 s

Time = 1200


Calculating the Saturation


Calculating the pressure head

DICPCG:  Solving for p, Initial residual = 0, Final residual = 0, No Iterations 0

Calculating the velocity


Calculating lamT

DILUPBiCG:  Solving for T, Initial residual = 0.0098683, Final residual = 8.24145e-07, No Iterations 38
ExecutionTime = 42.45 s  ClockTime = 49 s

@东岳 啊谢谢回复!我计算后，算例文件夹下没有postProcess文件，按理说应该有的吧，不知道您有什么建议。

我想监测温度场的最小值随时间的变化情况。能否直接有工具输入最小值的文件？
我用的OpenFOAM5.x，尝试采用cellSource来输出，将functions加入至controlDict文件中。但运行计算后并没有输出结果，百思不得其解。
您的任何建议，都万分感谢！

functions
{
    cellSource1
    {
        type              cellSource;
        function          ObjectLibs ("libfieldFunctionObjects.so");
        enabled           true;
        outputControl     outputTime;
        log               true;
        valueOutput       true;
        source            all;
        operation         min;
        fields
        (
            T
        );
    }
}

@myler 用这个代码是将最大值输出至cellSource1的文件夹吗？我也把这个代码添加至了controlDict里面，但是并没有看到最大值的输出结果啊