CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Courant number absurdly high - HELP

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 17, 2014, 23:28
Default Courant number absurdly high - HELP
  #1
Member
 
JuNNioR
Join Date: Feb 2014
Location: Brazil
Posts: 44
Rep Power: 12
jrsilvio_ver is on a distinguished road
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);
}
}
jrsilvio_ver is offline   Reply With Quote

Old   February 18, 2014, 17:15
Default
  #2
Member
 
JuNNioR
Join Date: Feb 2014
Location: Brazil
Posts: 44
Rep Power: 12
jrsilvio_ver is on a distinguished road
No suggestion?
jrsilvio_ver is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


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


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