twoPhaseEulerFoam中的IATE模型系数



  • 各位前辈,在使用双流体IATE模型时发现,内部读入的表征碰撞聚集的IATE系数Crc、Cti好像并未参与计算,更改这两个系数,计算结果没有差异,计算代码如下,最后返回tR,但是赋值是给了R,这个是不是有问题,还请各位老师指点一下,谢谢

    Foam::tmp<Foam::volScalarField>//计算Cti,破碎项的
    Foam::diameterModels::IATEsources::turbulentBreakUp::R() const
    {
    
        tmp<volScalarField> tR
        (
            new volScalarField
            (
                IOobject
                (
                    "R",
                    iate_.phase().U().time().timeName(),
                    iate_.phase().mesh()
                ),
                iate_.phase().U().mesh(),
                dimensionedScalar("R", dimless/dimTime, 0)
            )
        );
    
        volScalarField R = tR();
    
        scalar Cti = Cti_.value();
        scalar WeCr = WeCr_.value();
        volScalarField Ut(this->Ut());
        volScalarField We(this->We());
        const volScalarField& d(iate_.d()());
    
        forAll(R, celli)
        {
            if (We[celli] > WeCr)
            {
                R[celli] =
                    (1.0/3.0)
                   *Cti/d[celli]
                   *Ut[celli]
                   *sqrt(1 - WeCr/We[celli])
                   *exp(-WeCr/We[celli]);
            }
        }
        return tR;
    }
    
    Foam::tmp<Foam::volScalarField>//计算Crc,碰撞聚集项的
    Foam::diameterModels::IATEsources::randomCoalescence::R() const
    {
    
        tmp<volScalarField> tR
        (
            new volScalarField
            (
                IOobject
                (
                    "R",
                    iate_.phase().U().time().timeName(),
                    iate_.phase().mesh()
                ),
                iate_.phase().U().mesh(),
                dimensionedScalar("R", dimless/dimTime, 0)
            )
        );
    
        volScalarField R = tR();
    
        scalar Crc = Crc_.value();
        scalar C = C_.value();
        scalar alphaMax = alphaMax_.value();
        volScalarField Ut(this->Ut());
        const volScalarField& alpha = phase();
        const volScalarField& kappai = iate_.kappai();
        scalar cbrtAlphaMax = cbrt(alphaMax);
     
    
        forAll(R, celli)
        {
            if (alpha[celli] < alphaMax - SMALL)
            {
    
                scalar cbrtAlphaMaxMAlpha = cbrtAlphaMax - cbrt(alpha[celli]);
    
                R[celli] =
                    (-12)*phi()*kappai[celli]*alpha[celli]
                   *Crc
                   *Ut[celli]
                   *(1 - exp(-C*cbrt(alpha[celli]*alphaMax)/cbrtAlphaMaxMAlpha))
                   /(cbrtAlphaMax*cbrtAlphaMaxMAlpha);
    
            }
        }
    
        return  tR;
    }
    
    
    Foam::tmp<Foam::volScalarField>//计算Cwe,这个有作用
    Foam::diameterModels::IATEsources::wakeEntrainmentCoalescence::R() const
    {
        return (-12)*phi()*Cwe_*cbrt(CD())*iate_.a()*Ur();
    }
    


  • 更改这两个系数,计算结果没有差异

    因为不一定处于破碎区。里面存在临界微博数的判断。



  • 东岳老师我通过添加输出得到了,无论对于Crc还是Cti这两项都是进行了计算,在Crc、Cti那个函数里面输出一个给定单元格R[1927]值,这时候的输出都是有值的,但是在下面函数里面输出时发现,Crc和Cti计算得到的R值均是0,是不是我上面提到的Crc和Cti的函数返回是有问题的,还是有什么其他东西,另外即使临界韦伯数判断,那对于Crc这个临界体积分数判断,应该大部分的单元格还是满足它的判定条件的,希望各位前辈指点一二,不胜感激!

      forAll(sources_, j)
        {
    Info << "\nStarting R\n" << endl;
    Info << R[1927] << endl;
            R -= sources_[j].R();
    Info << R[1927] << endl;
        }
    


  • 各位前辈,经过验证发现,对于IATE模型中表示气泡聚合和破碎的Crc和Cti并未参与计算,因为这两个函数的返回值是tR,而进行赋值计算的是R,经过连续输出发现,赋值计算的R无论是否为0,返回值tR始终为0,因此对于Crc、Cti无论改变为多大的值对于气泡的聚合破碎毫无作用。
    在将返回值改为R以后,由于给定算例的连续相速度始终为0,而Crc和Cti计算R表达式中用到了连续相速度,所以还会出现返回值一直为0的情况,若改变算例,连续相速度给一定值后,就会发现此时返回值是tR时,返回值还一直为0,但返回值是R时,发现此时返回值就发生了变化,粒径分布也就发生了变化。因此在想是不是这个程序返回值tR是写错了还是我有哪个地方搞错了,请各位老师指点一下。另外还遇到一个问题就是在linux系统下,将返回值改正为R以后会报错如下:(以上所说的测试是基于Windows版本的绿化版本),希望各位前辈给点建议,这个错误应该怎么处理,不胜感激!

    --> FOAM FATAL ERROR: 
    different patches for fvPatchField<Type>s
    
        From function void Foam::fvPatchField<Type>::check(const Foam::fvPatchField<Type>&) const [with Type = double]
        in file /opt/openfoam4/src/finiteVolume/lnInclude/fvPatchField.C at line 209.
    
    FOAM aborting
    
    #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    #1  Foam::error::abort() at ??:?
    #2  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::operator-=(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
    #3  Foam::diameterModels::IATE::correct() at ??:?
    #4  Foam::twoPhaseSystem::correct() at ??:?
    #5  ? at ??:?
    #6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    #7  ? at ??:?
    已放弃 (核心已转储)
    


  • @upc_nghtwoPhaseEulerFoam中的IATE模型系数 中说:

    由于给定算例的连续相速度始终为0

    这个不应该啊?为什么是0



  • @李东岳 这个地方我表述错了,但是在算例中即使把Crc、Cti调整到100000000这个量级对计算结果还是毫无影响,即使存在临界数的判断,Crc这个依靠相分数判断的碰撞聚集项肯定会发生变化。还有就是对R赋值,返回却是tR,这个地方对吗,我现在是这个地方不理解,谢谢东岳老师的解答。



  • 把这一行

    volScalarField R = tR();
    

    改成

    volScalarField& R = tR.ref();
    

    再试试



  • 怎么样 有进展么



  • 东岳老师不好意思,我现在试试,刚看到您的回复贴



  • @李东岳 东岳老师我刚刚试了一下,还是不行,会报同样的错误



  • 我测试了一下,可以编译



  • @李东岳 对的,可以编译,但是运行算例还是会报错,和我以前一样,编译通过,运行出错



  • 4楼那个错误么,我测试了下运行也没有错误,你是不是动别的东西了



  • 那老师方便把您改动好的代码发我一下吗,我对着检查再看看,谢谢老师



  • 把这一行

    volScalarField R = tR();
    

    改成

    volScalarField& R = tR.ref();
    

    就这个 没别的了



  • 那最后一行还是return tR吗?



  • :142:



  • 这种情况,我一会再测试一下Crc和Cti是否参与了计算,谢谢老师百忙之中抽空解答我的疑惑,谢谢:xinxin3:


Log in to reply