Setting up: Cube, 2 holes, high pressure inside
 

Hello!

I am trying to set up a very simple cube (6*6*6m), of which two sides are holes. Inside the cube I set a field of 1m³ with very high pressure via setFields to simulate a detonation. (10,000,000 to 10,000,000,000 Pa)

But it always crashes after about a millisecond.

I try to solve this with the rhopSonicFoam solver. Since this is my first 3D case, I am not sure if I set the right boundary conditions for the walls and the holes. First, in the blockMeshDict I declare the :

patch outlet
(
(1 2 6 5)
(0 4 7 3)
)
wall walls
(
(0 1 5 4)
(5 6 7 4)
(3 7 6 2)
(0 3 2 1)
)

The boundary conditions in the boundary file are:

outlet
{
type patch;
nFaces 1800;
startFace 78300;
}
walls
{
type wall;
nFaces 3600;
startFace 80100;
}

Now I am really unasure what kinds of bondary conditions I have to set up in the /0 folder. Here's what I assume:

p)

outlet
{
type fixedValue;
value uniform 100000;
}
walls
{
type zeroGradient;
}

T)

outlet
{
type fixedValue;
value uniform 300;
}
walls
{
type fixedValue;
value uniform 300;
}

U)

outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
walls
{
type fixedValue;
value uniform (0 0 0);
}

Is this setup correct or did I set something wrong?



The error is:
ExecutionTime = 181.34 s ClockTime = 182 s

Time = 0.000486941

Max Courant Number = 0.486019
deltaT = 1.75373e-19
DILUPBiCG: Solving for rho, Initial residual = 4.76267e-15, Final residual = 2.21084e-28, No Iterations 1
DILUPBiCG: Solving for rhoUx, Initial residual = 1.77789e-13, Final residual = 7.59548e-27, No Iterations 1
DILUPBiCG: Solving for rhoUy, Initial residual = 2.84341e-13, Final residual = 1.21384e-26, No Iterations 1
DILUPBiCG: Solving for rhoUz, Initial residual = 1.34266e-13, Final residual = 6.23753e-27, No Iterations 1
DILUPBiCG: Solving for rhoE, Initial residual = 2.98644e-09, Final residual = 1.4146e-22, No Iterations 1
DICPCG: Solving for p, Initial residual = 4.46359e-15, Final residual = 7.24712e-31, No Iterations 1
#0 Foam::error::printStack(Foam::Ostream&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 void Foam::divide<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhopSonicFoam"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::operator/<Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhopSonicFoam"
#6 main in "/home/simulation/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/rhopSonicFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 _start at /usr/src/packages/BUILD/glibc-2.9/csu/../sysdeps/x86_64/elf/start.S:116
Gleitkomma-Ausnahme


Greets