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 公式应该没问题，它这个处理好像就是在径向有一个微弱的速度衰减。

[梦中飞翔的阿涛]

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