“欧拉-拉格朗日”喷管气固两相流模拟计算发散



  • 如题,我当前的研究方向是喷管中的气固两相流动(是xx研究所的喷管就不放图了.jpg),采用的是“欧拉-拉格朗日”的研究方法,采用的工具是openfoam+liggghts+CFDEM。目前正在参考论文https://www.researchgate.net/publication/317369528_CFDEM_Modelling_of_Particle_Heating_and_Acceleration_in_Cold_Gas_Spraying进行复现工作(论文中的求解器需要单独的编译后才能使用)

    在编译好论文中的求解器后进行计算,按照论文中的信息进行设置(部分信息),却发现流场怎么都算不收敛(根据我的经验一般计算发散都是流场而非颗粒场造成的)
    ![0_1591019287928_1591019232(1).jpg](正在上传 100%)

    我看到最近东岳老师也在做流固耦合相关的工作,但是算法并不用的这一套,不知道东岳老师有木有兴趣一起搭个伙呀!



  • 不知道为啥上传图片出错,我把Terminal终端中的错误信息直接粘出来

    [0] --> FOAM FATAL ERROR: 
    [0] Maximum number of iterations exceeded
    Foam::error::printStack(Foam::Ostream&)[0] 
    [0]     From function Foam::scalar Foam::species::thermo<Thermo, Type>::T(Foam::scalar, Foam::scalar, Foam::scalar, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar, Foam::scalar) const, Foam::scalar (Foam::species::thermo<Thermo, Type>::*)(Foam::scalar) const) const [with Thermo = Foam::eConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
    [0]     in file /home/smithlee/OpenFOAM/OpenFOAM-4.x/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.
    [0] 
    FOAM parallel run aborting
    [0] 
    [0] #0  Foam::error::printStack(Foam::Ostream&) at ??:?
    [2] #1  Foam::error::abort() at ??:?
    [0] #1  Foam::error::abort() at ??:?
    [2] #2  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
    [0] #2  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
    [2] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
    [0] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::eConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
    [0] #4   at ??:?
    [2] #4  mainmain in "/home/smithlee/CFDEM/CFDEMcoupling/platforms/linux64GccDPInt32Opt/bin/cfdemSolverRhoPimple"
    [0] #5  __libc_start_main in "/home/smithlee/CFDEM/CFDEMcoupling/platforms/linux64GccDPInt32Opt/bin/cfdemSolverRhoPimple"
    [2] #5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
    [0] #6   in "/lib/x86_64-linux-gnu/libc.so.6"
    [2] #6  ?? in "/home/smithlee/CFDEM/CFDEMcoupling/platforms/linux64GccDPInt32Opt/bin/cfdemSolverRhoPimple"
    --------------------------------------------------------------------------
    MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD
    with errorcode 1.
    
    NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
    You may or may not see output from other processes, depending on
    exactly when Open MPI kills them.
    --------------------------------------------------------------------------
     in "/home/smithlee/CFDEM/CFDEMcoupling/platforms/linux64GccDPInt32Opt/bin/cfdemSolverRhoPimple"
    [ubuntu:26692] 1 more process has sent help message help-mpi-api.txt / mpi-abort
    [ubuntu:26692] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
    


  • 这种可压缩的问题很常见。需要仔细检查下边界条件。

    是的,我下半年就要做可压缩喷管的模拟 :baobao: 如果进度正常的话,如果没有临时加活的话



  • @东岳 多谢东岳老师指点:tishizi: 我再改一改边界条件算一下,搞定了再发帖,用这套软件的道友太难找了:zoule:


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2