|
[Sponsors] |
2D-axisymmetric case - influence wedge angle on accurancy |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
F.F.
Join Date: Dec 2011
Posts: 14
Rep Power: 14 ![]() |
Hi,
assuming a simple 2D-axisymmetric case (side = symmetry planes, no swirl) what effect has the chosen wedge-angle on the results? Referring to the documentation 1 to 5 degree is preferable. From OpenFOAM I know 5 deg is recommended. Regarding my case: transient simulation of thermal stratification in a cylindrical tank. I do get "better" results (less convection, less movement in the layer) if my model is only 1deg. Mesh is quite fine and timesteps as well (1e-2s). From commen sense it's somehow clear that, e.g. a 1 deg model introduces a smaller error than a 5deg model. But what's the detailled explanation and why is 5 deg recommended? Shouldn't the recommended angle rather be a function of mesh-density, transport variables and numerical error? Looking forward to get some answers ![]() |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,816
Rep Power: 144 ![]() ![]() ![]() ![]() |
When you make the wedge angle smaller it means you more closely approximate the full circular shape and this reduces error. But you also make the quality of the mesh worse - as you now have a small angled element in the central axis. Different simulations have different sensitivities to poor mesh quality, so some simulations will run OK with a 1° wedge angle, but more sensitive ones will never converge.
So the compromise is better resolution of the full geometry against decreasing numerical stability. |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 1,173
Rep Power: 23 ![]() |
Axis-Symmetric is one major fallback of CFX.
One thing I do when meshing, if it is appropriate: 1.) Leave out a very small portion of the center, so that you have an annular space instead of a wedge. Place a symmetry plane on that small surface. Now you can sweep quads 1 element through the thickness. 2.) Another thing I do if I can't leave out the center, is use a mesh method similar to this: it requires more than 1 element through the thickness unfortunately, but maybe only 2,4,or 6. I haven't perfected this yet obviously, but it is better than doing a full 90 degrees. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
F.F.
Join Date: Dec 2011
Posts: 14
Rep Power: 14 ![]() |
Thanks for the input! I wasn't aware of that issue.
So the recommendation for 2d-axisymmetric cases should rather be: 1) Use a 5°-Mesh and make sure the case is converging (model is allright). 2) Verify results from 1) with a 1°-Mesh. 3) If case is diverging, rsp. mesh quality is too bad, improve meshing near axis. Luckily, the center of my geometry is already empty. So I will just change all my cases from 5 to 1 degree. Erik, your mesh from 2) looks interesting! I keep that in mind. For another case I had quite bad elements in the center, but CFX was alright with that. Coarsen the mesh along the axis helped. Smoothing hexa-meshes made it worse (icem). Also deleting degenerating elements didn't help. Anyway, I had good experiences with the collapseEdges utility in OpenFOAM after creating axisymmetric meshes from a quasi-2D/3D-icem Mesh. But I couldn't reimport it to CFX.. |
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
Join Date: May 2018
Posts: 1
Rep Power: 0 ![]() |
To aid in the convergence, you could use the 5deg mesh case as initial condition for the 1deg mesh case.
|
|
![]() |
![]() |
![]() |
Tags |
2d axisymmetric, angle, wedge |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
a serious problem on a case with attack angle | anne | CFX | 0 | November 18, 2008 05:28 |
Heat flux BC in axisymmetric case | kevin | FLUENT | 0 | February 7, 2008 18:07 |
2D case, UDF for finding y1/2 in axisymmetric jet | sagar | FLUENT | 0 | November 1, 2006 04:41 |
Free surface boudary conditions with SOLA-VOF | Fan | Main CFD Forum | 10 | September 9, 2006 12:24 |
B C for 2D axisymmetric FV case | Amith | Main CFD Forum | 7 | June 20, 2002 18:06 |