从湍流模型的correct函数中，输出浮力项、流线弯曲修正项等

[lizhisongsjtu]

对SST湍流模型作修正，引入了浮力项N2和流线弯曲修正项Fr1。希望在输出k、omega的时候，同步输出N2和Fr1.请问要怎么写？

尝试了N2.write();，这个语句把每个dt的结果全都输出了。

请问要怎么改写？

    // Calculate the Brunt-Vaisala frequency
    volScalarField N2 = gField_&fvc::grad(rho1)/rho1;
    N2.write(); // this exports data every time setp.
    //curvature correction
    const scalar cr1_ = 1.0;
    const scalar cr2_ = 2.0;
    const scalar cr3_ = 1.0;
    volScalarField Fr1    
    (
        max
    	(
            min
            (
                (scalar(1.0) + cr1_)*2.0*rStar /
                (scalar(1)+rStar) *
                (scalar(1.0) - cr3_*atan(cr2_*rTilda))-
                cr1_
                ,
                scalar(1.25)
            )
            ,
			scalar(0.0)
        )
    );

    // Turbulent kinetic energy equation
    tmp<fvScalarMatrix> kEqn
    (
        fvm::ddt(alpha, rho1, k_)
      + fvm::div(rhoPhi, k_)
      - fvm::laplacian(alpha*rho1*DkEff(F1), k_)
     ==
        alpha()*rho1()*Pk(G) * Fr1() // curvature correction
      - fvm::SuSp((2.0/3.0)*alpha()*rho1()*divU, k_) 
      - fvm::Sp(alpha()*rho1()*epsilonByk(F1, F23), k_)
      - fvm::Sp(nut*alphaBS_*N2*rho1/max(k_,this->kMin_),k_) // Buoyancy production term

      + rho1*kSource()
      + fvOptions(alpha, rho1, k_)
    );

[李东岳]

volScalarField N2("N2", gField_&fvc::grad(rho1)/rho1);

然后在算例层面，通过输出缓存场的方法输出 http://dyfluid.com/code.html

[lizhisongsjtu]

可能是版本不同的问题，在OF2206中，添加缓存没有成功。

用另一种方法，算是曲线完成任务了。

在H文件中，添加 volScalarField NN;，在C文件初始化的列表中，加入

 NN
(
    IOobject
    (
        "NN",
        this->runTime_.timeName(),
        this->mesh_,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
     ),
    this->mesh_
),

然后在恰当的位置，完成赋值 NN = N2 即可。

不过有一个缺陷，需要再初始文件夹中设定 NN 的初始文件。

(先解决问题，麻烦就麻烦点吧)