|
[Sponsors] |
Setting up: Cube, 2 holes, high pressure inside |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
October 30, 2009, 12:27 |
Setting up: Cube, 2 holes, high pressure inside
|
#1 |
Member
Dennis Rown
Join Date: Oct 2009
Posts: 51
Rep Power: 17 |
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. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Pressure BC for combustion chamber | Giuki | FLUENT | 1 | July 19, 2011 12: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 |