“欧拉-拉格朗日”喷管气固两相流模拟计算发散
-
如题,我当前的研究方向是喷管中的气固两相流动(是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