[enginpower]

Hello all,
I'm trying to simulate butane gas flowing through the geometry in the following picture sonicFoam.jpg which reach supersonic flow in some throats connecting in tube and out tube.
All patches not indicated in the picture, are signed like walls.

I want to simulate real gas, and I calculated all parameters to get it using Peng-Robinson equation of state as you can check in the files: shared.zip

You can also find inside them all applied boundary conditions. Actually it is the first time I'm using sonicFoam solver, so I'm not so expert about that.
Running the simulation I got this error, but it's not clear to me where is the problem:

Maximum number of iterations exceeded functions are : 

15
(
abort
coded
patchProbes
probes
psiReactionThermoMoleFractions
removeRegisteredObject
residuals
rhoReactionThermoMoleFractions
setTimeStep
sets
surfaces
systemCall
timeActivatedFileUpdate
writeDictionary
writeObjects
)



    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::janafThermo<Foam::PengRobinsonGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy>]
    in file /home/ubuntu/OpenFOAM/OpenFOAM-4.1/src/thermophysicalModels/specie/lnInclude/thermoI.H at line 66.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::janafThermo<Foam::PengRobinsonGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#4  ? at ??:?
#5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6  ? at ??:?
Aborted (core dumped)

All other similar errors in the forum are different from that because it lists some functions but I've not understood where to set them. And what exeeds 15 iterations? Not clear this point.

Please give me some useful ways to get a good simulation of that, mainly on the boundary condition settings

Thanks

[enginpower]

I solved some errors.
Now it runs but I think it doesn't go over the first iteration. Here the log file is:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : sonicFoam
Date   : Apr 21 2017
Time   : 14:19:44
Host   : "engine-VirtualBox"
PID    : 10986
Case   : /home/engine/OpenFoam/2DdobleTube
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          janaf;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

Creating field kinetic energy K

No MRF models present

No finite volume options present


Starting time loop

Time = 5e-07

Courant Number mean: 0.576416 max: 108.333
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0892473, No Iterations 1001
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.31323, No Iterations 1001
DILUPBiCG:  Solving for e, Initial residual = 1, Final residual = 9.86494e-10, No Iterations 189

It's not clear to me if it performs the first iteration or not. But what does it mean

DILUPBiCG:  Solving for Ux, Initial residual = 1, Final residual = 0.0892473, No Iterations 1001
DILUPBiCG:  Solving for Uy, Initial residual = 1, Final residual = 0.31323, No Iterations 1001
DILUPBiCG:  Solving for e, Initial residual = 1, Final residual = 9.86494e-10, No Iterations 189

in the last lines? Does it runs 189 iterations?
What can I do to understand what does not work in the simulation? Finally I chose a kepsilon model as starting model, but the boundary conditions are almost as set in the attached files (see above).
Hope it can help who wants to suggest me something interesting to learn new things

Thanks anyway