CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Setting up: Cube, 2 holes, high pressure inside

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2009, 12:27
Question Setting up: Cube, 2 holes, high pressure inside
  #1
Member
 
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 7
shangzung is on a distinguished road
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:rintStack(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:perator/<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

Last edited by shangzung; October 30, 2009 at 12:52.
shangzung is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pressure BC for combustion chamber Giuki FLUENT 1 July 19, 2011 11:35
Setting dynamic pressure as the inlet BC ruined FLUENT 1 March 26, 2008 06:30
pressure gradient term in low speed flow Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 11:52
Too high pressure force ? Tim Franke Main CFD Forum 1 March 15, 2000 00:36
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 01:43.