@李东岳

FoamFile
{
    version     2.0;
    format      ascii;
    class       scalarField;
    location    "constant/../postProcessing";
    object      totalLiquid;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

1(0.2)

麻烦老师再帮我看一下，为什么每次都只会输出一个时刻的值，不会随时间增加输出的

@hongjiewang 在 OpenFOAM小代码 中说：

@李东岳

  totalLiquid
  {        
        libs            (utilityFunctionObjects);
        type            coded;
        name            totalLiquid;
        enabled         true;
        writeControl    timeStep;
        writeInterval   1;
        
        codeOptions
        #{
            -I$(LIB_SRC)/meshTools/lnInclude
        #};

        codeExecute
        #{

            const volScalarField& b =
                mesh().lookupObject<volScalarField>("alpha.liquid");

            IOList<scalar> liquidFraction
            (
                IOobject
                (
                    "liquidFraction",
                    mesh().time().constant(),
                    "../postProcessing",
                    mesh(),
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                1
            );
            liquidFraction[0] = b.weightedAverage(mesh().V()).value();
        #};
    }

东岳老师，我这样添加到controlDict里，运行后postProcessing里并没有出现liquidFraction的值。

需要修改两个部分，
1.在下方添加 .write（）
2.将codeExecute改为codeWrite
之后就可以得到想要的相含量结果。

@李东岳

  totalLiquid
  {        
        libs            (utilityFunctionObjects);
        type            coded;
        name            totalLiquid;
        enabled         true;
        writeControl    timeStep;
        writeInterval   1;
        
        codeOptions
        #{
            -I$(LIB_SRC)/meshTools/lnInclude
        #};

        codeExecute
        #{

            const volScalarField& b =
                mesh().lookupObject<volScalarField>("alpha.liquid");

            IOList<scalar> liquidFraction
            (
                IOobject
                (
                    "liquidFraction",
                    mesh().time().constant(),
                    "../postProcessing",
                    mesh(),
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                1
            );
            liquidFraction[0] = b.weightedAverage(mesh().V()).value();
        #};
    }

东岳老师，我这样添加到controlDict里，运行后postProcessing里并没有出现liquidFraction的值。

@李东岳 在 OpenFOAM小代码 中说：

IOField<scalar> utau
(
IOobject
(
"utau",
runTime.constant(),
"../postProcessing",
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
scalarField(totalFSize,0.0)
);

请问老师，如果我想要在postProcessing中输出一个标量，他只是一个数，并不是场量，我应该怎么定义他的类型。上面这种方法是不是只适应于场量的输出。我想输出的是下面c的数值。

const volScalarField& b = mesh().lookupObject<volScalarField>("alpha.liquid");
scalar c= b.weightedAverage(mesh().V()).value();

@Samuel-Tu 您好，您例子中输出的是一个场量，如果我只是想把一个数值输出到postProcessing中，应该定义一个什么样的类型（类似IOField）。

const volScalarField& b = mesh().lookupObject<volScalarField>("alpha.liquid");
scalar c= b.weightedAverage(mesh().V()).value();

请问，我如何把终端输出文件中的liquid fraction单独输出到一个log文件，并实时输出他的变化曲线。

同时请教一下，我想把计算过程中的liquid fraction单独输出到一个文件中，应该如何操作

模型：
在下壁面同时进行加热并控制只在界面处蒸发的情况下，在界面相变的时候，出现了一些更容易相变的点，使界面锯齿状，请问这种现象应该怎么解释。

@史浩 按照老师说的用alpha识别相界面，不管我设置的那个横截面是否把网格切分，比如我设置z轴0.0002m一下为液体区域，当z轴方向50个网格时，初始计算他的液相含量是0.2；当我z轴方向51个网格时，初始计算他的液相含量变成了0.196078，导致0文件中alpha非0即1，因此并没有识别出相界面。

非常感谢老师的回复

@史浩 在 界面相变 中说：

相含量的变化不一定只和源项有关，还可能和宏观流动、求解误差等相关

从运算出来的test1看，即使规定它不相变，网格内的相含量也是有变化的，这也导致我采用梯度是否为0判断相界面的不准确性，和老师所说的一致。目前我将test1==0修改为一个范围值，也不知道是否正确。

不用相分数，来判断相界面

整个计算空间是1✖50✖1mm（网格划分为50✖250✖50），初始时利用setFields将1✖50✖0.2mm区域设置为液体区，即0时刻，每个网格的液体含量非0即1，所以可能没办法利用相分数判断0时刻的相界面。