|
[Sponsors] |
February 17, 2014, 23:28 |
Courant number absurdly high - HELP
|
#1 |
Member
JuNNioR
Join Date: Feb 2014
Location: Brazil
Posts: 44
Rep Power: 12 |
Dear, good night.
I implemented the Energy Equation in the icoFoam solver. However the solution is not converging, basically the Courant number is absurdly high. Could anyone help me find the error? The boundary conditions are reasonable? I greatly appreciate everyone's attention. ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- blockMeshDict FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 1; vertices ( (0 0 0) //Pt 0 (0.4 0 0) //Pt 1 (0.4 0.6 0) //Pt 2 (0 0.6 0) //Pt 3 (0 0 5) //Pt 4 (0.4 0 5) //Pt 5 (0.4 0.6 5) //Pt 6 (0 0.6 5) //Pt 7 (0.4 0.8 0) //Pt 8 (0 0.8 0) //Pt 9 (0.4 0.8 5) //Pt 10 (0 0.8 5) //Pt 11 (0.4 1 0) //Pt 12 (0 1 0) //Pt 13 (0.4 1 5) //Pt 14 (0 1 5) //Pt 15 (0.6 0 0) //Pt 16 (0.6 0.6 0) //Pt 17 (0.6 0 5) //Pt 18 (0.6 0.6 5) //Pt 19 (0.6 0.8 0) //Pt 20 (0.6 1 0) //Pt 21 (0.6 0.8 5) //Pt 22 (0.6 1 5) //Pt 23 (1 0 0) //Pt 24 (1 0.6 0) //Pt 25 (1 0 5) //Pt 26 (1 0.6 5) //Pt 27 (1 0.8 0) //Pt 28 (1 1 0) //Pt 29 (1 0.8 5) //Pt 30 (1 1 5) //Pt 31 ); blocks ( hex (0 1 2 3 4 5 6 7) (20 30 20) simpleGrading (1 1 1) // Bloco 1 hex (3 2 8 9 7 6 10 11) (20 10 20) simpleGrading (1 1 1) // Bloco 2 hex (9 8 12 13 11 10 14 15) (20 10 20) simpleGrading (1 1 1) // Bloco 3 hex (1 16 17 2 5 18 19 6) (10 30 20) simpleGrading (1 1 1) // Bloco 4 hex (8 20 21 12 10 22 23 14) (10 10 20) simpleGrading (1 1 1) // Bloco 5 hex (16 24 25 17 18 26 27 19) (20 30 20) simpleGrading (1 1 1) // Bloco 6 hex (17 25 28 20 19 27 30 22) (20 10 20) simpleGrading (1 1 1) // Bloco 7 hex (20 28 29 21 22 30 31 23) (20 10 20) simpleGrading (1 1 1) // Bloco 8 ); edges ( arc 6 19 (0.5 0.5585786 5) arc 19 22 (0.6414214 0.7 5) arc 22 10 (0.5 0.8414214 5) arc 10 6 (0.3585786 0.7 5) arc 17 2 (0.5 0.5585786 0) arc 20 17 (0.6414214 0.7 0) arc 8 20 (0.5 0.8414214 0) arc 2 8 (0.3585786 0.7 0) ); boundary ( saida { type patch; faces ( (0 3 2 1) (1 2 17 16) (16 17 25 24) (3 9 8 2) (17 20 28 25) (9 13 12 8) (8 12 21 20) (20 21 29 28) ); } entrada { type patch; faces ( (4 5 6 7) (5 18 19 6) (18 26 27 19) (7 6 10 11) (19 27 30 22) (11 10 14 15) (10 22 23 14) (22 30 31 23) ); } paredeadiabatica { type wall; faces ( (0 4 7 3) (3 7 11 9) (9 11 15 13) (26 24 25 27) (27 25 28 30) (30 28 29 31) (15 14 12 13) (14 23 21 12) (23 31 29 21) (4 0 1 5) (5 1 16 18) (18 16 24 26) ); } queimador { type wall; faces ( (10 8 20 22) (6 19 17 2) (19 22 20 17) (6 2 8 10) ); } ); mergePatchPairs ( ); ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- P FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 101325; boundaryField { paredeadiabatica { type zeroGradient; } queimador { type zeroGradient; } entrada { type zeroGradient; } saida { type fixedValue; value uniform 101325; } } ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- T FoamFile { version 2.0; format ascii; class volScalarField; object T; } //**********************************// dimensions [0 0 0 1 0 0 0]; internalField uniform 300; boundaryField { entrada { type fixedValue; value uniform 300; } saida { type zeroGradient; } paredeadiabatica { type zeroGradient; } queimador { type fixedValue; value uniform 1200; } } ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- U FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { entrada { type fixedValue; value uniform (0 0 -1); } saida { type fixedValue; value uniform (0 0 -0.2); } queimador { type fixedValue; value uniform (0 0 0); } paredeadiabatica { type fixedValue; value uniform (0 0 0); } } |
|
February 18, 2014, 17:15 |
|
#2 |
Member
JuNNioR
Join Date: Feb 2014
Location: Brazil
Posts: 44
Rep Power: 12 |
No suggestion?
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Sudden jump in Courant number | NJG | OpenFOAM Running, Solving & CFD | 7 | May 15, 2014 13:52 |
dynamic Mesh is faster than MRF???? | sharonyue | OpenFOAM Running, Solving & CFD | 14 | August 26, 2013 07:47 |
High Courant number | arunsmec | Main CFD Forum | 4 | March 28, 2013 06:04 |
Stable boundaries | marcoymarc | CFX | 33 | March 13, 2013 06:39 |
very high courant number | dsanza | OpenFOAM Running, Solving & CFD | 2 | August 17, 2011 02:29 |