[andreas0209@hotmail.com]

Hi everybody

I am trying to simulate a particle flow along a tube (like the vertical channel).
The flow is compressible and exactly here is my problem.
If I keep my boundary condition like a incompressible flow (same inlet and outlet condition the simulation works good. As soon as I change the parameter to a compressible flow I got problems with the thermo physics.

--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double) const) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3  Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::calculate() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#4  Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#5  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
Aborted (core dumped)

For my case I actually don't need thermo pysical properties, why I wanted to switch them off. But exactly this is obviously not simple. The best is to recompile the solver without the thermo physical properties. I am quite desperate and I was wondering if somebody could help me with that?

Should be not complicated for someone who is more experienced like me. The simpleReactingParcelFoamIPM foder is attached.

Regards,
Andreas

[chpjz0391]

Originally Posted by andreas0209@hotmail.com

Hi everybody

I am trying to simulate a particle flow along a tube (like the vertical channel).
The flow is compressible and exactly here is my problem.
If I keep my boundary condition like a incompressible flow (same inlet and outlet condition the simulation works good. As soon as I change the parameter to a compressible flow I got problems with the thermo physics.

--> FOAM FATAL ERROR: 
Maximum number of iterations exceeded

    From function thermo<Thermo, Type>::T(scalar f, scalar T0, scalar (thermo<Thermo, Type>::*F)(const scalar) const, scalar (thermo<Thermo, Type>::*dFdT)(const scalar) const, scalar (thermo<Thermo, Type>::*limit)(const scalar) const) const
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.3.0/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 76.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::error::abort() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2  Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::T(double, double, double, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double, double) const, double (Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>::*)(double) const) const in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#3  Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::calculate() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#4  Foam::heRhoThermo<Foam::rhoReactionThermo, Foam::SpecieMixture<Foam::reactingMixture<Foam::polynomialTransport<Foam::species::thermo<Foam::hPolynomialThermo<Foam::icoPolynomial<Foam::specie, 8>, 8>, Foam::sensibleEnthalpy>, 8> > > >::correct() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#5  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7  
 in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/simpleReactingParcelFoam"
Aborted (core dumped)

For my case I actually don't need thermo pysical properties, why I wanted to switch them off. But exactly this is obviously not simple. The best is to recompile the solver without the thermo physical properties. I am quite desperate and I was wondering if somebody could help me with that?

Should be not complicated for someone who is more experienced like me. The simpleReactingParcelFoamIPM foder is attached.

Regards,
Andreas

have you solved ur problem, I met the same question as yours.

[andreas0209@hotmail.com]

Originally Posted by chpjz0391

have you solved ur problem, I met the same question as yours.

Unfortunately, I didn't. Sorry.

[chpjz0391]

you said when you kept your boundary condition like an incompressible flow, the simulation worked well, in my case , I just need to simulate the incompressible flow for steady-state, but after I changed my BC ,it can only work for 45 steps. I am quite new of OF, so I have no idea what's wrong with my case. Could you tell me how did you set your BC as the incompressible. Could you help me with that? thank you very much The attachment is my case.

simpleReactingParcelFoam Case.tar.gz

[chpjz0391]

Originally Posted by andreas0209@hotmail.com

Unfortunately, I didn't. Sorry.

This issue almost drives me crazy. I need to solve this problem before next Monday, could you give me some advice?

[Chrisi1984]

Hi,

I had a very quick look inside your case. I realizied that you are using p=0 as outlet condition. So I think you assumed this to be the relative pressure like it is common for incompressible cases. But very likely you need to define the pressure level for that solver with absolute pressue so e.g. for environment conditions p=100000. (Check also tutorial case here: https://github.com/OpenFOAM/OpenFOAM...calChannel/0/p

Kind regards and good luck
Chrisi

[chpjz0391]

Originally Posted by Chrisi1984

Hi,

I had a very quick look inside your case. I realizied that you are using p=0 as outlet condition. So I think you assumed this to be the relative pressure like it is common for incompressible cases. But very likely you need to define the pressure level for that solver with absolute pressue so e.g. for environment conditions p=100000. (Check also tutorial case here: https://github.com/OpenFOAM/OpenFOAM...calChannel/0/p

Kind regards and good luck
Chrisi

Thank you very much chrisi, I found that mistake, I will fix it and try again