CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Problem with symmetry boundary condition

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 30, 2017, 07:21
Default Problem with symmetry boundary condition
  #1
New Member
 
Somashekhar Kulkarni
Join Date: Oct 2017
Posts: 14
Rep Power: 9
somas123 is on a distinguished road
The actual geometry is a hollow cylinder and I'll be using the axisymmetric condition, hence I've taken only half part of the cylinder and rotated it by 5 degrees. I've highlighted the problem areas in the attached pic. The initial problem was that at these highlighted corners, 7-node hexa elements were created. So CFX wasn't accepting the mesh. Hence I used the divide 7 node hexas into pyramids option. Now I'm able to import the mesh but CFX generates extra surfaces (which are the faces of the mesh element) at those corner nodes. Because of these extra surfaces I'm unable to use the symmetry BC in my simulation. I even tried assigning the symmetry BC to these extra surfaces but the solver still gave an error. What am I supposed to do?
Attached Images
File Type: jpg geometry.JPG (25.8 KB, 12 views)
File Type: jpg geom 2.JPG (28.7 KB, 17 views)
somas123 is offline   Reply With Quote

Old   October 30, 2017, 09:01
Default
  #2
Senior Member
 
JuPa's Avatar
 
Mr CFD
Join Date: Jun 2012
Location: Britain
Posts: 361
Rep Power: 15
JuPa is on a distinguished road
A nice trick is to replace the axisymmetric "pizza slice" with a "wedge". The advantage of this is you completely avoid pyramid and triangular cells. You can get a high quality hexahedral mesh.

So instead of the geometry finishing to a point like it does in your example, make it finish to a tiny flat surface.

This paper did exactly what I'm referring to: http://www.tandfonline.com/doi/abs/1...32.2015.994459
JuPa is offline   Reply With Quote

Old   October 30, 2017, 18:48
Default
  #3
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
What meshing software did you use to generate this mesh?

I do not understand the original problem. What is a 7 node hexa? Do you mean a 8 node hex with one degenerate edge where two nodes are coincident?

If you take a quad meshed surface and rotate it about one edge you get 8 node hex for most of the geometry and 6 node wedges along the rotation edge.

You have to be careful with Mr CFD's suggestion. While that approach does eliminate the problem at the rotation axis, it does change the geometry. Whether it is significant or not depends on what you are doing.
ghorrocks is offline   Reply With Quote

Old   October 31, 2017, 00:57
Default
  #4
New Member
 
Somashekhar Kulkarni
Join Date: Oct 2017
Posts: 14
Rep Power: 9
somas123 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
What meshing software did you use to generate this mesh?

I do not understand the original problem. What is a 7 node hexa? Do you mean a 8 node hex with one degenerate edge where two nodes are coincident?

If you take a quad meshed surface and rotate it about one edge you get 8 node hex for most of the geometry and 6 node wedges along the rotation edge.

You have to be careful with Mr CFD's suggestion. While that approach does eliminate the problem at the rotation axis, it does change the geometry. Whether it is significant or not depends on what you are doing.
I used ICEM CFD (ANSYS 15 version). In the 2nd image you can see that the extra surfaces formed are basically 2d mesh elements. That element's one side is on the corner (located on the axis of symmetry) of the geometry whereas the other edges are inside. Therefore a total of 7 nodes are available for the element instead of six or eight. I am getting such elements at the corners because of the O-grid that I have used to specify first cell height. CFX wasn't accepting the mesh because of these 7 node elements. Hence I had to divide them into pyramids and now it is accepting the mesh but it shows these extra surfaces which are clubbed under primitive 2D.
somas123 is offline   Reply With Quote

Old   October 31, 2017, 01:26
Default
  #5
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
In that case, in ICEM try defining the elements which are at the end of the sweep as parts, and make sure the correct element faces are in the parts and it should flow through to CFX.

In ICEM if you do a rotate mesh operation it should generate correct wedge elements on the rotation axis. This would avoid this problem entirely.
ghorrocks is offline   Reply With Quote

Old   October 31, 2017, 01:44
Default
  #6
New Member
 
Somashekhar Kulkarni
Join Date: Oct 2017
Posts: 14
Rep Power: 9
somas123 is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
In that case, in ICEM try defining the elements which are at the end of the sweep as parts, and make sure the correct element faces are in the parts and it should flow through to CFX.

In ICEM if you do a rotate mesh operation it should generate correct wedge elements on the rotation axis. This would avoid this problem entirely.
I have generated my mesh using a similar method. First I created a 2d block by using sweep method. Then I meshed only the surface, to which the 2d block was associated, using pre-mesh. After this, I converted the 2d block to 3D using rotate option and then updated the mesh. The main reason why I can't get a wedge at that corner is because the O-grid is dividing the block at 45 deg. at that corner node.
Attached Images
File Type: jpg o-grid.JPG (19.5 KB, 9 views)
somas123 is offline   Reply With Quote

Old   November 1, 2017, 02:39
Default
  #7
New Member
 
Somashekhar Kulkarni
Join Date: Oct 2017
Posts: 14
Rep Power: 9
somas123 is on a distinguished road
Is there any method by which I can convert the mesh in those problematic areas to tet mesh and then merge it with the hex mesh in other areas?
somas123 is offline   Reply With Quote

Old   November 1, 2017, 06:17
Default
  #8
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,871
Rep Power: 144
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
ICEM can handle doing that, and much more. There is not much mesh editting you cannot do in ICEM. But ICEM is a complex package and takes a while to get used to. Make sure you have done the tutorial examples (from the ANSYS Community website).
ghorrocks 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
Wind turbine simulation Saturn CFX 60 July 17, 2024 06:45
Symmetry boundary condition problem jhtcten CFX 4 April 20, 2017 01:00
problem of implementing a new solver with a special boundary condition huyidao OpenFOAM Programming & Development 0 April 27, 2015 09:24
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00
symmetry boundary condition icemaniac178 CFX 3 March 13, 2011 06:40


All times are GMT -4. The time now is 21:53.