Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    关于二维后台阶模型的壁面摩擦阻力系数Cf的计算

    OpenFOAM
    2
    2
    1481
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • A
      ahujeffery last edited by 李东岳

      各位大神,刚刚接触OpenFOAM不久,使用pisoFOAM,采用k-e湍流模型计算二维后台阶的算例,雷诺数为37500,计算后与实验对比数据结果,台阶后侧y方向的速度剖面与实验较为符合。
      在计算下壁面摩擦阻力系数 Cf 时,使用tecplot对流场的速度进行微分计算,进而求得切应力$\tau$,除以$0.5\rho\mathbf{U}\mathbf{U}$后得到Cf,其与实验结果有较大差异。
      不知道是否是tecplot计算的问题?OpenFOAM对于这种切应力有没有什么计算的方法。
      同样的算例在fluent中求解后的Cf结果与实验对比尚可。
      望能够得到大家的解答!

      1 Reply Last reply Reply Quote
      • 李东岳
        李东岳 管理员 last edited by 李东岳

        OpenFOAM对于这种切应力有没有什么计算的方法。

        剪切应力的定义是固定的,用什么软件的方程都是一样的。你要从流场中获得$\tau$?印象中OpenFOAM里面没有计算tau的现成的工具。如果你不确定,你可以和Fluent比较其他量看看。

        如果你要计算tau,你可以写个后处理小工具。比如这个

        bool Foam::functionObjects::Lambda2::calc()
        {
            if (foundObject<volVectorField>(fieldName_))
            {
                const volVectorField& U = lookupObject<volVectorField>(fieldName_);
                const tmp<volTensorField> tgradU(fvc::grad(U));
                const volTensorField& gradU = tgradU();
                const volTensorField SSplusWW
                (
                    (symm(gradU) & symm(gradU))
                  + (skew(gradU) & skew(gradU))
                );
                return store
                (
                    resultName_,
                   -eigenValues(SSplusWW)().component(vector::Y)
                );
            }
            else
            {
                return false;
            }
        }
        
        

        https://coding.net/u/dyfluid/p/OpenFOAM-5.x/git/blob/master/src/functionObjects/field/Lambda2/Lambda2.C

        线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
        CFD高性能服务器 http://dyfluid.com/servers.html

        1 Reply Last reply Reply Quote
        • First post
          Last post

        CFD中文网 | 东岳流体 | 京ICP备15017992号-2
        论坛登录问题反馈可联系 li.dy@dyfluid.com