[decibelle]

Hi,

I work on a subsonic case for external aerodynamic. when I run my case with rhoPimpleFoam after few iteration my simulation crash and I can see this error:

[14] smoothSolver:  Solving for e, Initial residual = 9.39539e-06, Final residual = 9.39539e-06, No Iterations 0

[14]
[14] --> FOAM FATAL ERROR:
[14] Maximum number of iterations exceeded
[14]
[14]     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::hConstThermo<Foam::perfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
[14]     in file /Software/OpenFOAM/OpenFOAM-v1606+/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.
[14]
FOAM parallel run aborting
[14]
[14] #0  Foam::error::printStack(Foam::Ostream&)--------------------------------------------------------------------------
An MPI process has executed an operation involving a call to the
"fork()" system call to create a child process.  Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your MPI job may hang, crash, or produce silent
data corruption.  The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.

The process that invoked fork was:

  Local host:          cfdnode08 (PID 39476)
  MPI_COMM_WORLD rank: 14

If you are *absolutely sure* that your application will successfully
and correctly survive a call to fork(), you may disable this warning
by setting the mpi_warn_on_fork MCA parameter to 0.
--------------------------------------------------------------------------
 in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[14] #1  Foam::error::abort() in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[14] #2  Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy>::TEs(double, double, double) const in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so"
[14] #3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so"
[14] #4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfluidThermophysicalModels.so"
[14] #5  ? in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
[14] #6  __libc_start_main in "/lib64/libc.so.6"
[14] #7  ? in "/Software/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 14 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.
--------------------------------------------------------------------------

I checked my boundary condition and my fvSolution but I don't understand what's wrong.

Someone can help me?

I joined at this post the fvsolution file and my boundary condition

Thanks in advance

fvSolution.docx

U.docx

alphat.docx

epsilon.docx

k.docx

[decibelle]

this is the other file of my boundary condition

nut.docx

p.docx

T.docx

[guanjiang.chen]

Hi decibelle,

In the fvsolution, diagonal seems to be enough for the rho. In your pressure condition, zeroGradient inlet didn't not cooperate well with inletoutlet in my experiences. I would like recommend outletinlet or freestreampressure for the outlet.

Regards,
Guanjiang