|
[Sponsors] |
about calculations of Courant number in openfoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
August 2, 2013, 06:14 |
about calculations of Courant number in openfoam
|
#1 |
Senior Member
Join Date: Jan 2013
Posts: 372
Rep Power: 14 |
Hi All,
I am wondering why the factor 0.5 is needed in the expression of Courant number in Openfoam? scalar CoNum = 0.5*gMax(sumPhi/mesh.V().field())*runTime.deltaTValue(); scalar meanCoNum = 0.5*(gSum(sumPhi)/gSum(mesh.V().field()))*runTime.deltaTValue(); Removing the factors seems correct when the above expressions are compared with the definition of CFL number. Thank you very much if someone can give some comments. |
|
August 2, 2013, 08:11 |
|
#2 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
For a cubic cell, you will add up the magnitude of the velocities on the six faces (because of surfaceSum). So by the factor 0.5, you make sure you use the average value of the velocity of opposite face when calculating (u_x/dx+u_y/dy+u_z/dz)*dt.
|
|
August 2, 2013, 11:49 |
|
#3 |
Senior Member
Join Date: Jan 2013
Posts: 372
Rep Power: 14 |
Thank you very much! Here the velocity is the normal velocity for each face that encloses a complete cell. Is my understanding correct? If yes, how can we directly add the velocity to get the average values? If the cell is tetrehral, how is this averaging performed? Thank you!
|
|
August 2, 2013, 11:52 |
|
#4 |
Senior Member
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21 |
Let me remind you that phi=U_f . A_f, so normal velocity times area. This is why you divide by volume. How is the Courant number defined for tet-cells?
|
|
August 2, 2013, 11:57 |
|
#5 |
Senior Member
Join Date: Jan 2013
Posts: 372
Rep Power: 14 |
In openfoam, it seems that the Courant number is not particular for tet cells. So in that case, how does the factor 0.5 come from? Thanks.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[mesh manipulation] Mesh Refinement | Luiz Eduardo Bittencourt Sampaio (Sampaio) | OpenFOAM Meshing & Mesh Conversion | 42 | January 8, 2017 12:55 |
decomposePar pointfield | flying | OpenFOAM Running, Solving & CFD | 28 | December 30, 2013 15:05 |
dynamic Mesh is faster than MRF???? | sharonyue | OpenFOAM Running, Solving & CFD | 14 | August 26, 2013 07:47 |
Fluent 6.4 courant number | Aris Nikolopoulos | FLUENT | 0 | May 6, 2008 08:52 |
LES near wall model & courant number | kasim | CFX | 5 | March 16, 2008 18:23 |