利用codeStream设置相分数初始场

chengan.wang

大家好，
我想利用codeStream设置相分数初始场，比如盐水-淡水按照线性分层变化，但是编译总出错，1-Solution.pdf ，请大家提供点其他参考资料吧！谢谢

internalField #codeStream
{
    codeInclude
    #{
        #include "fvCFD.H"
    #};

    codeOptions
    #{
        -I$(LIB_SRC)/finiteVolume/lnInclude \
        -I$(LIB_SRC)/meshTools/lnInclude
    #};

    codeLibs
    #{
        -lmeshTools \
        -lfiniteVolume
    #};

    code
    #{
        const IOdictionary& d = static_cast<const IOdictionary&>(dict);
        const fvMesh& mesh = refCast<const fvMesh>(d.db());
        scalarField alpha.air(mesh.nCells(), 1.0);
        scalarField alpha.water(mesh.nCells(), 0.0);        
        scalarField alpha.saltwater(mesh.nCells(), 0.0);        
 
        scalar D = 0.0252;
        scalar ymin = -4.48*D;
        scalar ymax = 5.82*D; 
        scalar alpha_min = 0;
        scalar alpha_max = 1;       

        forAll(alpha.water,i)
        {
            const scalar x = mesh.C()[i][0];
            const scalar y = mesh.C()[i][1];
            const scalar z = mesh.C()[i][2];
            if ( (y >= ymin) && (y <= ymax) ) 
               {
                   alpha.water[i] = alpha_min+(alpha_max-alpha_min)*(y-ymin)/(ymax-ymin));
                   alpha.saltwater[i] = 1.0-alpha.water[i];
                   alpha.air[i] = 0;
               }                
        }  

    #};
};

队长别开枪

@chengan-wang 变量名不能包含dot

chengan.wang

@队长别开枪您好，我的算例涉及到三种相分数，alpha.air，alpha.water，alpha.saltwater，我该怎么定义呢？

Tens

codeStream可以在求解器里设置的吗，怎么还涉及到编译。我一直以为只能在初始条件设置，想问下codeStream是写在哪个文件里

Tens

@tens 如果是在初始条件里写的：1.假设在alpha.water文件里写的，不确定是否可同时对saltwater和air初始化。2.writeEntry这句没写，不确定是否有影响

CFD中文网

利用codeStream设置相分数初始场