如何在controlDict中创建一个场，不需要修改求解器

[学流体的小明]

话说这个latex代码怎么转义不了
使用LES计算气泡槽道流，因为是多相流，某一点可能会有气泡经过，解析的雷诺应力$u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$会受到很大影响，所以我想对 $u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$ 进行加权处理。
具体方法是：生成张量场，计算每一个cell每一个瞬时的$\alpha u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$，$(1-\alpha )u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$，$\rho u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$，然后使用fieldAverage进行系综平均。
我想在controlDict中生成场$\rho u{\mathrm{\prime }}_{i}u{\mathrm{\prime }}_j$，然后每个时刻计算一下这个场，怎样通过functions{ }操作呢？因为是要每一个时间步都计算以获得正确的值，所以不太想按一定时间间隔输出数据之后再进行后处理。
当然也可以直接在solver的createFields.H中创建，我就是想找一个更通用的不用修改求解器的方法，直接很简便地在controlDict - functions{ }中实现。
毕竟functions{ }可以用writeObjects输出

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Version:  v2012
    \\  /    A nd           | Website:  www.openfoam.com
     \\/     M anipulation  |
-------------------------------------------------------------------------------
Description
    Writes out specified objects, e.g. fields, stored on the case database.

\*---------------------------------------------------------------------------*/

type            writeObjects;
libs            (utilityFunctionObjects);

objects         (<object names>);

writeControl    writeTime;

// ************************************************************************* //

[18190691340]

大佬解决了吗？最近也想这样搞

[学流体的小明]

编程解决的，直接改求解器

[18190691340]

@学流体的小明 在 如何在controlDict中创建一个场，不需要修改求解器 中说：

编程解决的，直接改求解器
也只能这样了

[coolhhh]

OpenFOAM-v2206/tutorials/incompressible/pimpleFoam/LES/planeChannel/setups.orig/common/system/controlDict，有一个计算Cf场的函数，可以参考下：

    Cf
    {
        type                coded;
        libs                (utilityFunctionObjects);
        name                Cf;
        timeStart           $/timeStart;
        writeControl        writeTime;
 
        codeExecute
        #{
            auto* CfPtr =
                mesh().getObjectPtr<volScalarField>("Cf");
 
            if (!CfPtr)
            {
                Info<< "Create skin-friction coefficient field" << nl;
                CfPtr = new volScalarField
                (
                    IOobject
                    (
                        "Cf",
                        mesh().time().timeName(),
                        mesh(),
                        IOobject::NO_READ,
                        IOobject::AUTO_WRITE
                    ),
                    mesh(),
                    dimless
                );
 
                regIOobject::store(CfPtr);
            }
 
            auto& Cf = *CfPtr;
 
            Info<< "Computing skin-friction coefficient field\n" << endl;
 
            const auto& tau =
                mesh().lookupObject<volVectorField>("wallShearStress");
 
            const dimensionedScalar Ubulk(dimVelocity, 17.55);
 
            Cf = mag(tau.component(0))/(0.5*sqr(Ubulk));
        #};
    }