
[Sponsors] 
April 1, 2011, 07:04 
fvScheme and NonOrthogonal Correctors

#1 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Hi all,
I have to run a simulation with a mesh of prisms and hexas. This is my checkMesh: Create polyMesh for time = 0 Time = 0 Mesh stats points: 101178 internal points: 0 faces: 259884 internal faces: 118566 cells: 69762 boundary patches: 6 point zones: 0 face zones: 0 cell zones: 1 Overall number of cells of each type: hexahedra: 29640 prisms: 40122 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Face vertices OK. Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology obstacle 1560 3120 ok (nonclosed singly connected) bottom 49 100 ok (nonclosed singly connected) right 68 138 ok (nonclosed singly connected) top 49 100 ok (nonclosed singly connected) left 68 138 ok (nonclosed singly connected) defaultFaces 139524 101178 ok (nonclosed singly connected) Checking geometry... Overall domain bounding box (0 0 0) (0.00447 0.0062 5e05) Mesh (nonempty, nonwedge) directions (1 1 0) Mesh (nonempty) directions (1 1 0) All edges aligned with or perpendicular to nonempty directions. Boundary openness (1.63471e19 2.90979e18 4.90414e20) OK. Max cell openness = 2.19738e16 OK. Max aspect ratio = 13.3858 OK. Minumum face area = 8.1356e11. Maximum face area = 7.04863e09. Face area magnitudes OK. Min volume = 4.0678e15. Max volume = 3.16076e13. Total volume = 1.21253e09. Cell volumes OK. Mesh nonorthogonality Max: 31.2939 average: 4.74797 Nonorthogonality check OK. Face pyramids OK. Max skewness = 0.495827 OK. Mesh OK. What are the best setting for fvScheme and nonorthogonal corrector? Is the option "limited 0.5" a good choise? Something like this: Code:
gradSchemes { default cellLimited Gauss linear 1; } laplacianSchemes { default Gauss linear limited 0.5; } snGradSchemes { default limited 0.5; } And for the divergence what i have to choose? (Iam using interFoam) divSchemes { div(rho*phi,U) ?? div(phi,alpha) ?? div(phirb,alpha) ?? } and finally for the non orthogonal corrector in fvSolution, what is best choise? Thanks for any help andrea 

April 1, 2011, 08:19 

#2 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Your mesh has a good quality, so probably it would not be necessary to limit so much the gradients, neither to use a lot of non orthogonal correctors. I would suggest this setting as a first try:
GradSchemes: cellMDLimited Gauss linear 0.5 divSchemes(for "rho*phi,U"): Gauss linearUpwindV cellMDLimited Gauss linear 0.5 divSchemes(others): Gauss upwind LaplacianSchemes: the same as in your post SnGradSchemes: as above nNonOrthogonalCorrectors: 1 Regards V. 

April 1, 2011, 08:31 

#3 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Thanks for reply Vesselin,
what is the difference between cellMDLimited and cellLimited? However I will try your settings to see if it's good! Thanks again andrea 

April 1, 2011, 08:37 

#4  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Quote:
Regards V. 

April 5, 2011, 05:50 

#5 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Hi Vesselin,
i tried your suggestions in a simple testcase: the relaxation of a square drop. With an orthogonal mesh i can find the analytical solution without any problems. Then i tried to run the same case with a mesh made of prisms, but the results are very bad. The drop starts to deform in a nonphysical way (it should be circular) until it exits from the domain. Here my "checkMesh allTopology allGeometry". as you can see there is one failed check, while only "checkMesh" is ok (i checked it also in the previous post geometry and i got the same error, the mesh is build up with gmsh). maybe is that the problem? or have i to increase the limit on the schemes? Time = 0 Mesh stats points: 32734 internal points: 0 edges: 113611 internal edges: 15891 internal edges using one boundary point: 0 internal edges using two boundary points: 15891 faces: 113134 internal faces: 48146 cells: 32256 boundary patches: 2 point zones: 0 face zones: 0 cell zones: 1 Overall number of cells of each type: hexahedra: 0 prisms: 32256 wedges: 0 pyramids: 0 tet wedges: 0 tetrahedra: 0 polyhedra: 0 Checking topology... Boundary definition OK. Point usage OK. Upper triangular ordering OK. Topological cell zipup check OK. Face vertices OK. Faceface connectivity OK. <<Writing 476 cells with with single nonboundary face to set twoInternalFacesCells Number of regions: 1 (OK). Checking patch topology for multiply connected surfaces ... Patch Faces Points Surface topology Bounding box atm 476 952 ok (nonclosed singly connected) (0 0 0) (0.06 0.06 0.001) defaultFaces 64512 32734 ok (nonclosed singly connected) (0 0 0) (0.06 0.06 0.001) Checking geometry... Overall domain bounding box (0 0 0) (0.06 0.06 0.001) Mesh (nonempty, nonwedge) directions (1 1 0) Mesh (nonempty) directions (1 1 0) All edges aligned with or perpendicular to nonempty directions. Boundary openness (2.27697e19 9.96175e20 0) OK. Max cell openness = 2.1596e16 OK. Max aspect ratio = 10.4958 OK. Minumum face area = 5.28134e08. Maximum face area = 7.05771e07. Face area magnitudes OK. Min volume = 5.28134e11. Max volume = 2.09285e10. Total volume = 3.6e06. Cell volumes OK. Mesh nonorthogonality Max: 30.8774 average: 6.37461 Nonorthogonality check OK. Face pyramids OK. Max skewness = 0.521316 OK. Min/max edge length = 0.000311967 0.001 OK. All angles in faces OK. Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1 All face flatness OK. Cell determinant (wellposedness) : minimum: 0 average: 4.26433e35 ***Cells with small determinant found, number of cells: 32256 <<Writing 32256 underdetermined cells to set underdeterminedCells Failed 1 mesh checks. End Thanks andrea 

April 5, 2011, 08:54 

#6  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Quote:
Regards V. 

April 5, 2011, 09:26 

#7 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Hi Vesselin.
I have attached my test case. The structured grid has 22550 cells (all hexa, i found convergence to the analytical solution also with 3600 cells) and the unstructured one has 32256 cells (all prisms). The BC are the same in both case. I have also included the fvScheme and fvSolution for both cases. The dimensions of the domain are 6cmx6cm. The analytical solution is a drop of radius 2cm. as you can see from the movie the solution of the unstructured case is very bad. Any suggestion would be appreciated! andrea 

April 5, 2011, 09:26 

#8 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
ops, i forgot the .zip!
andrea 

April 5, 2011, 09:36 

#9 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
There is still no .zip attachement...


April 5, 2011, 09:40 

#10 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
sorry, was too large! Now is ok
andrea 

April 5, 2011, 10:12 

#11 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Well, first of all the structured mesh seems much finer, and this is not "fair": you'll have to refine a lot the unstructured mesh and give it a second try with the same fvSchemes and fvSolution settings (but with cellMDLimited option instead of cellLimited). Secondly, if the results are still far from the structured mesh response, you can try to use also on the (fine) unstructured the following divSchemes (leaving the linearUpwindV scheme on velocity):
div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; Regards V. 

April 5, 2011, 10:46 

#12 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Thanks again for your patience.
I found good results also on a "less" fine structured grid (3600 hexa, 60X60 as in the attached). The unstructured mesh is almost 9 times more defined. Note that the computational cost drastically increase in the unstructured case (few second for the structured mesh, 1h and a half for the unstructured case). So of course i can refine my unstructured mesh but is that the solution? I mean, Is the problem related to the "resolution" of the mesh, or to the selection of the right discretization schemes?..Or is there something else you have to take into account when you want to use an unstructured grid? from what I know, OF should work fine on unstructured mesh... However I will try to make the unstructured grid more fine and see what happens! Thanks Andrea 

April 5, 2011, 10:46 

#13 
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 14 
Hello Andrea,
Hi V, nice to hear you again. As far as I know, tet mesh will never be comparable in term of convergence and accuracy to a hexa mesh. This is due to the way interpolation schemes work, thus nothing will change using a very very fine tet mesh. The only way to achieve the same accurcay is to use the reconCentral scheme, which is provided on 1.6ext. Andrea, maybe this discussion is useful for you: pressure eq. "converges" after few time steps. It applies to a mainly tet mesh. V also contributed to it. cheers mad 

April 5, 2011, 12:02 

#14 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Hi mad,
really really interesting thread!!..thanks for the hint. I would like to run some test using your and V. suggestions, changing the fvScheme and fvSolution and see what happens. However i'm using interFoam and i do not know if the right schemes for simpleFoam are good also for interFoam, but I'll do a check! thanks again andrea 

April 6, 2011, 02:54 

#15  
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 14 
Quote:
mad 

April 6, 2011, 03:50 

#16 
Senior Member
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 182
Rep Power: 9 
Hi,
just a question from my side. Has anyone tried the recently published pointLinear scheme? Quote from https://github.com/OpenCFD/OpenFOAM...2918f2e1bb7904 "surfaceInterpolation: New pointLinear scheme First interpolates from the cells to the vertices and then from the face vertices to the face centres. This may be a useful scheme for calculating gradients on bad/tet meshes when used in conjunction with the Gauss gradient scheme. " Regards. Markus. 

April 7, 2011, 06:35 

#17 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Hi all, just some additonal comments:
For Andrea: as Mad said (nice to hear you too), tet meshes will never be exactly as accurate as hex ones. However, this not means that a cfd code cannot reach some remarkably good results with tetra (and, in general, unstructured) meshes, and OF indeed can (at least to my experience). Refining the mesh (whatever the type) will always improve accuracy, but of course this is not a universal solution to all problems (especially when you have a combination of stability and accuracy issues). When I suggested you to refine the unstructured mesh I was not really aware that it was so much finer than the hexa one, and of course a 9times finer tetra mesh is far enough with the respect to the hexa one...Anyway, the "trik" whit the unstructured meshes is, provided a sufficiently good mesh quality (which means no highly skewed cells neither too strong nonorthogonality), to find the most reasonably good fvSchemes/Solution setting, which unfortunately tends to vary with the type of problem that has to be solved (compressible/incompressible, laminar/turbulent, multiphase, reacting and/or burning flows, etc.). For Mad (and Andrea too): I don't think the proper fvSchemes/Solution setting could be the same for very different kind of problems (thus, of physical quantities that have to be solved). For instance, I suggested to Andrea a setting which worked very well for incompressible HighRe turbulent flows (my case), but for his multiphase problem they seemed not so good...In addition, I tried the reconCentral scheme, but honestly it does not seem a so extraordinary improvement: in my computations there were practically no differences between results obtained with the reconCentral scheme and, for instance, the linearUpwind scheme with poperly set limiters on gradients. Best Regards V. 

April 8, 2011, 03:26 

#18 
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 296
Rep Power: 8 
Hi V., mad and Markus
..some questions... 1)What do you mean for no highly skewed cells and too strong nonorthogonality? Can you give an example of high and low values? 2)I agree with V. on the fact that the fvScheme and fvSolution could hardly be the same for different problem. First of all the equations are different in multiphase problem (you have an additional term which takes into account the surface forces), the velocity normally are low, no gravity is a good approximation..etc etc.. so i guess that i have to find my "proper" schemes. For example i found that the linearUpwind for the div(rho*phi,U) would not be not a good choise but it would seem better to use Gauss limitedLinearV (i also found some Henry's comments that supports this thing). At the moment i ran some test and i could find a not too bad solution (is not a perfect circular bubble but at least is fixed at the centre of the domain and does not go away) using leastSquare limited 0.333 for the gradient and the same of V. post for all other quantities, but limited to 0.333. Now i want to check the pressure solution (and the convergence) to see if it is close to the analytical one. So my question is: How can i print out the residual for each iteration? (From the log i can see the initial and the final residual and total number of iteration for every loop) 3)I'm using OF 1.7, is there the reconCentral schemes? Because i can't find it on the manual. 4)with respect to the comment by markus, also the pointLinear schemes is not there. Can i just add the .H and .C in src/finiteVolume and then recompile the library? Thanks andrea 

April 8, 2011, 04:53 

#19  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 364
Rep Power: 11 
Quote:
2) If you mean a plot of the residuals, you can simply run the foamLog <log> utility (with <log> name of your logfile) and then plot whatever you need inside the logs folder which is created by the utility itself (I use gnuplot for this) 3) The reconCentral scheme is an additional feature of the ext release of OpenFOAM (thus it is not part of OF1.7.0/1). However, thanks to Alberto Passalacqua, you can download the source code and compile it for your OF version from here: Wing Aerodynamics Fluent OF 15 comparison Hope this helps V. 

April 8, 2011, 06:21 

#20  
Senior Member
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 182
Rep Power: 9 
Hi,
Quote:
Markus 

Thread Tools  
Display Modes  

