coalChemistryFoam改写

[Rachel0096]

我在对coalChemistryFoam改写，想把它和随机场方程进行结合。简单来说是对于连续相用一系列随机场进行计算，颗粒相不用（单独计算）。但这就遇到了一个问题，我结合着输运方程来说明：

    for(label i=0; i<Y.size(); i++)
    {
        {
            volScalarField Yi   = Ysf[nsf][i];
            volScalarField Ymean = Y[i];

            WienerTermYi = pow(2*gamma,0.5)*( fvc::grad(Yi) & dW_sf);
            fvScalarMatrix YiEqn
            (
                  fvm::ddt(rho, Yi)
                + mvConvection->fvmDiv(phi, Yi)
                - fvm::laplacian(turbulence->muEff()/Sct, Yi)
                ==
                  coalParcels.SYi(i, Ymean)
                + combustion->R(Ymean)
                + fvOptions(rho, Ymean)
		        + (rho*WienerTermYi/deltaT)
            );
            YiEqn.relax();
            fvOptions.constrain(YiEqn);
            YiEqn.solve(mesh.solver("Yi"));
            YiEqn.solve();
            fvOptions.correct(Yi);

            Yi.max(0.0);
	        Ysf[nsf][i] = Yi;
            Yt += Yi;
        }
    }

在这个方程里，Yi指的是随机场的组分项（比如名为N2__01），Ymean指的是平均场的组分项（也就是组分场，名为N2）。颗粒的求解不考虑随机场，来自颗粒的源项被输入到每一个随机场中。但这里就会有一个问题，在计算时，会出现类似incompatible fields for operation [N2] + [N2__01]的报错，我的理解是因为这两个名称不同所导致的。这个问题要怎么去解决？我想的是，或许可以把这三项在方程外部导入到一个临时变量，再输入到方程里，但是不知道怎么用代码去实现。求各位大佬指点。