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

[ldymig]

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

在编译好论文中的求解器后进行计算，按照论文中的信息进行设置（部分信息），却发现流场怎么都算不收敛（根据我的经验一般计算发散都是流场而非颗粒场造成的）

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

[ldymig]

不知道为啥上传图片出错，我把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

[李东岳]

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

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

[ldymig]

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