[Calmly]

Hello FOAMers!

I am trying to simulate something like an explosion in an open field. It's a 2D simulation with the explosive on the ground. The explosive is simulated as a "balloon" with high pressure.

Although, this is the error I get:

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

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

Reading field U

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;
}

No finite volume options present

fluxScheme: Kurganov

Starting time loop

Mean and max Courant Numbers = 0.00026747 0.0990834
deltaT = 1.1999e-09
Time = 1.1999e-09

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Ux, Initial residual = 7.38269e-06, Final residual = 5.18693e-20, No Iterations 5
smoothSolver: Solving for Uy, Initial residual = 7.38239e-06, Final residual = 4.48882e-20, No Iterations 5
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for e, Initial residual = 1.46889e-05, Final residual = 1.36414e-19, No Iterations 5


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

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:erfectGas<Foam::specie> >; Type = Foam::sensibleInternalEnergy; Foam::scalar = double; Foam::species::thermo<Thermo, Type> = Foam::species::thermo<Foam::hConstThermo<Foam:er fectGas<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:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#4 ? at ??:?
#5 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6 ? at ??:?
Aborted

The B.C. are:

epsilon

dimensions [0 2 -3 0 0 0 0];

internalField uniform 1000000;

boundaryField
{
// sides
// {
// type fixedValue;
// value uniform 266000;
// }
sides
{
type inletOutlet;
inletValue uniform 1000000;
value uniform 1000000;
}
symmetry
{
type symmetryPlane;
}
ground
{
type epsilonWallFunction;
value uniform 1000000;
}
defaultFaces
{
type empty;
}
}

k

dimensions [0 2 -2 0 0 0 0];

internalField uniform 2400;

boundaryField
{
// inlet
// {
// type fixedValue;
// value uniform 1000;
// }
sides
{
type inletOutlet;
inletValue uniform 2400;
value uniform 2400;
}
symmetry
{
type symmetryPlane;
}
ground
{
type kqRWallFunction;
value uniform 2400;
}
defaultFaces
{
type empty;
}
}

p

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
sides
{
type zeroGradient;
}
symmetry
{
type symmetryPlane;
}
ground
{
type zeroGradient;
}

defaultFaces
{
type empty;
}
}

U

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
sides
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
symmetry
{
type symmetryPlane;
}
ground
{
type fixedValue;
value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
}

fvSchemes

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

fluxScheme Kurganov;

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default none;
div(tauMC) Gauss linear;
div(phi,epsilon) Gauss linearUpwind grad(epsilon);
div(phi,k) Gauss linearUpwind grad(k);

}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
// reconstruct(rho) Minmod;
// reconstruct(U) MinmodV;
// reconstruct(T) Minmod;
flux(rhoU) linear;
dotInterpolate(S,tauMC) linear;
interpolate(muEff) linear;
interpolate(rho) linear;
}

snGradSchemes
{
default corrected;
}

I would really appreciate if you could help me!

[dli]

[QUOTE=Calmly;664350]Hello FOAMers!

Hi, Calmly

Have you solved your problem?
I'm not sure where the exactly the problem are, but with this high compressible problem, I would choose the upwind or the limitedLinear in the grad schemes.

[Taataa]

The boundary conditions doesn't seem right to me. On the sides I would suggest to use fixedFluxPressure for p and inletOutlet or slip for U, assuming that sides are far enough from the region of interest.

pressureInletOutletVelocity is used when pressure is prescribed in the boundary.

Also, I suggest to use upwind for schemes at first to make sure model works. Once you get reasonable results you can use higher order schemes.