codedFixedValue这样写对吗？

尚善若水

请教各位老师，在文献中看到，射流入口用了双曲正切函数分布，我用codedFixedValue这样写对吗？
我编译可以通过，但是 plot over line看了下分布，和直接给速度定值喷射好像没有差别。

文献定义

文献中的速度分布

    inletJet
    {
        type            codedFixedValue;
        value           uniform (0 0 0); 

        name            tanhInletJet;    

        code
        #{
            const fvPatch& boundaryPatch = patch();
            const vectorField& Cf = boundaryPatch.Cf();
            vectorField& field = *this;

            scalar Uj_y = 1250, r_0 = 0.001, theta_0 = 0.00009;   

            forAll(Cf, faceI)
            {
                const scalar x = Cf[faceI][0]; 
                const scalar z = Cf[faceI][2]; 
                const scalar r = sqrt(x*x+z*z); 

                field[faceI] = vector(0,Uj_y*0.5*(1-tanh(1/(4*theta_0)*((r/r_0)-(r_0/r)))),0); 
            }
        #};

        // codeInclude
        //#{
        //    #include "fvCFD.H"
        //    #include <cmath>
        //    #include <iostream>
        //#};

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

wangfei9088

@尚善若水代码没错，是公式（3.2）的原因。分析一下：
当 r=r0 时，tanh[..]=0, w/Wj=0.5;当r=0时，tanh[...]=-1,w/Wj=1。这个曲线是可以用excel或matlab画出的。
我觉得这个公式对应的应该是图中Elliptic (major axis)这条曲线，实质上也就是你说的和直接给速度定值喷射一样。要得到Circular曲线也可以，把公式（3.2）中的r除以2即可。

尚善若水

@wangfei9088 公式应该没问题，它这个处理好像就是在径向有一个微弱的速度衰减。

梦中飞翔的阿涛

大佬，我刚刚发了一个帖子，能不能帮我看一下，关于煤粉燃烧的。之前你帮我回答过相关问题，感谢了。