CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > ANSYS Meshing & Geometry

Meshing Sudden Contraction and Expansion

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 23, 2010, 01:01
Default Meshing Sudden Contraction and Expansion
  #1
New Member
 
Robert Manning
Join Date: Mar 2010
Posts: 2
Rep Power: 0
carldaru is on a distinguished road
Hello,

I am trying to mesh a 2D sudden contraction and expansion using Gambit. However, I am not sure what type of meshing to use. Every type except "Quad Pave" keeps kicking up some sort of error for me. Either it's too skewed, or if I try using triangles, the triangles are too big, or I get a negative volume, etc. When I used "Quad Pave," the analysis that came out was not what was expected. Is it possible to mesh 2D figures that have a conical contraction and sudden expansion with plain cartesian meshes that literally go parallel and perpendicular to the length of the pipe in a perfect grid pattern (i.e. a 10 x 100 grid). If the perfect cartesian grid mesh is not possible, what type of meshing do you guys suggest I use? Also, do I just mesh the face or just the edges, or both?

Also, when I want to turn this into 3D, do I just do a revolve and then use the same type of mesh? Would that work, or is something else needed?

Thank you.
carldaru is offline   Reply With Quote

Old   March 23, 2010, 01:56
Default
  #2
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
you can control the mesh refinment with size functions.
In 2d I assume you can handle your problem with quads. (maybe you need to split your domain to achieve it)
you may post a picture.
Regarding 3d model, then you can revolve your faces, hex-cooper meshing schema should be enabled (just take care with face connectivity)
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   March 24, 2010, 13:28
Default
  #3
New Member
 
Robert Manning
Join Date: Mar 2010
Posts: 2
Rep Power: 0
carldaru is on a distinguished road
Quote:
Originally Posted by -mAx- View Post
you can control the mesh refinment with size functions.
In 2d I assume you can handle your problem with quads. (maybe you need to split your domain to achieve it)
you may post a picture.
Regarding 3d model, then you can revolve your faces, hex-cooper meshing schema should be enabled (just take care with face connectivity)
So by "size functions", I assume you are referring to the interval count/interval size controls? If I wanted a grid of say 10 x 100 (100 lines perpendicular to the length of the pipe), then would I select the edges appropriately and just type in 10 and 100?

When you say "split your domain", what are you referring to, and how does one go about doing that?

Once I get the 2D meshing down I will try the 3d model with the hex-cooper meshing. I assume face connectivity has to do with the splitting of the domain as referenced above? Am I correct in saying this?

A picture is posted with this message.

Thank you.
Attached Images
File Type: jpg Nozzle.jpg (15.9 KB, 23 views)
carldaru is offline   Reply With Quote

Old   March 25, 2010, 02:55
Default
  #4
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
In your case, just split your domain into logical and mappable surfaces.
See the picture (I just meshed one half (upper) side)
Sans titre.jpg
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   April 13, 2010, 13:02
Default need help
  #5
New Member
 
hassan
Join Date: Jun 2009
Posts: 14
Rep Power: 9
hassan79 is on a distinguished road
which turbulent model used can i have your mobile no
hassan79 is offline   Reply With Quote

Old   January 2, 2017, 13:03
Default Difficulty in cooper mesh hexahedral in sudden contraction
  #6
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
Hi everyone,
I am trying to model oil-water flow through two ducts, encountered a sudden contraction. I attached you my geometry. My intention is to inject oil in face 2 and water in face 1. It is very important to me to have a refined mesh near boundary. if i mesh the upstream pipe with tetrahedral and downstream with hexahedral mesh, everything goes well. however, when I want to use both pipe hexahedral mesh with defining boundary layer, Gambit cannot make it an error occurs regarding to source face. first, I defined two cylindrical geometry, then I split downstream pipe from face3. I meshed face 1 with refining and also face2.then I tried to sweep to these face (as sources) to the upstream pipe, but face 3 does not allow. any help would be appreciated.
Attached Images
File Type: png PICTURE.png (21.2 KB, 7 views)
mk091088 is offline   Reply With Quote

Old   January 2, 2017, 13:10
Default difficulty in creating cooper mesh in sudden contraction
  #7
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
I also attached you the error associated with cooper mesh.
Attached Images
File Type: jpg PICTURE2.jpg (73.8 KB, 5 views)
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 02:14
Default
  #8
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
you need to split your domain,
Split it at face 3
Then split the cylinder which contains faces 1 & 2 in 2 concentric cylinders.
The common diameter from those 2 cylinders will be non-constant (use the intermediate edges from face 1 & 2 to create the surface split)
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 4, 2017, 07:41
Default
  #9
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
Thank you Max for your response.
I did not understand what you said about concentric cylinders. Let's consider upstream pipe as volume 1, and downstream as volume 2. As you said, I split volume 1 from face 3 to create face 4 so that I am able to define face 3 and 4 as interior and wall, respectively.
You said ''Then split the cylinder which contains faces 1 & 2 in 2 concentric cylinders''. The faces 2 and 3 do not have the same diameters. The faces 2 and 3 have the diameters 15 mm and 21 mm, respectively. Honestly, I did not get what you said about spliting volume 1 into two concentric cylinders.
Attached Images
File Type: png picture 3.png (29.0 KB, 8 views)
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 08:13
Default
  #10
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
Quote:
Originally Posted by mk091088 View Post
I am able to define face 3 as interior
you don't need to set it as interior since it will automatically done

Regarding the split I mentionned, create an edge from both vertices as mentionned in the picture
Capture.JPG

Once it is done, split each circular edge in its middle and create another straight edge.

Now you can create the conical surface with help from both straight edges and 2 semi-circles.

Split the volume with the conical surface

It should do the job
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 4, 2017, 15:13
Default
  #11
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
As I mentioned before, I want to model two phase flow in contraction pipe. Thus, the transport equation in volume 1 must be solved for mixture (1 fluid). The method that you suggested is perfect, but I am curious because creating the conical surface would create an extra surface inside the volume 1 (there are two volumes in volume 1 which included the annular between wall and conical surface and another inside the conical surface), and definitely it should not be wall and it must be interior, according to my physical problem. I focused first on volume 1, I was properly able to separate face 1 from 2 and face 3 from face 4. Furthermore, I created two semi-circles on face 1 and face 3 and two straight edges as you said. but I do not know how to make a conical surface. It must be two semi-conical surfaces and then should be united, right?
Sorry to ask you two many questions, I am new to Gambit.
Attached Images
File Type: png image 5.png (25.3 KB, 6 views)
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 15:18
Default
  #12
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
If you split the volume as I described, don't set anything on the conical surface.
Fluent will treat it as interior.
So it will understand you only have 1 fluid volume.
(the split is only here for enabling you a full hex mesh)

In your picture delete your 2 straight lines, and create 2 new ones using the vertex (already existing)
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 4, 2017, 15:34
Default
  #13
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
Thank you. I do know if it is true. Please verify. I deleted previous edges, and created two straight edges with vertexes and I used two semi-circle on two different edges and two straight lines to create one conical surface as you said. I used Geometry/Face/Wireframe.
Attached Images
File Type: png image 6.png (18.6 KB, 5 views)
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 15:37
Default
  #14
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
yes
now check if topology is proper by meshing this surface (map)
it will give you a good idea (sometimes Gambit creates disturbed topology )

If it is correct, do the same with the other surface. And finally split the volume with those 2 surfaces
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 4, 2017, 16:13
Default
  #15
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
I meshed the conical surface with mapped mesh, it is ok. Thus, I start meshing the domain. first I created a boundary layer for annulus in face 1, then I used a mappable quad for face 1. I used a paved quad mesh for face 2. For volume mesh, I still have problem because when I swept all the surface meshes face 1 and 2 through the domain by cooper mesh, the edges (which was used to create conical surfaces) does mot allow to mesh the volume properly. Moreover, the pannel related to source faces is deactivated and I needed to define manually which I do not know which source face I must use. It is worth noting that I did not unit two semi-conical surfaces before meshing.
Attached Images
File Type: jpg testing the mesh surface.jpg (93.5 KB, 4 views)
File Type: jpg meshing the inlet surface.jpg (62.8 KB, 4 views)
File Type: jpg cooper mesh.jpg (95.1 KB, 6 views)
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 16:21
Default
  #16
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
first don't use (yet) any BL

delete all mesh
Mesh one annulus face with quad
Mesh the conical face with quad
Mesh the volume (annulus), enforce the cooper scheme if necessary by picking both annulus surfaces as source.
If you still have problems, I will check it when I go back to the office (monday)
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 4, 2017, 17:37
Default
  #17
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
It is done. some important points to be addressed: I did what you said with interval count=42 for faces and volumes. here are the results of meshes:

For volume 1 (annulus)=21,168
For volume 1 (conical volume)=26,586
For volume 2=322,197

I do not know why there is a significant difference between upstream and downstream mesh in terms of the quantity with the same interval count.
Most importantly, I still need to apply boundary layers to my walls in volumes 1 and 2.
mk091088 is offline   Reply With Quote

Old   January 4, 2017, 19:06
Default
  #18
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
If your volumes are well meshed (good quality); save your *.dbs file and add the BL.


Sent from my iPhone using CFD Online Forum mobile app
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   January 5, 2017, 08:25
Default
  #19
New Member
 
parham
Join Date: Oct 2015
Posts: 17
Rep Power: 3
mk091088 is on a distinguished road
These are final questions, Max:

1)Boundary layers must be added to faces 1 and face 3 at the inlet and the point of singularity, right? In this case how boundary layer can be swept all along the wall inside the volumes? I mean the boundary layer can be only applied to specific faces not the whole volume.

2)After adding boundary, should I export mesh, right?

3)The second question is regarding the difference between the number of elements in volume 1 and 2. Why there is so much difference. According to your method (Upstream pipe, volume 1 is so coarse and downstream is fine). The mesh quality, however, is good.
Attached Images
File Type: jpg final mesh.jpg (24.7 KB, 6 views)
File Type: png boundary layer.png (19.3 KB, 4 views)
mk091088 is offline   Reply With Quote

Old   January 5, 2017, 08:42
Default
  #20
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,172
Rep Power: 32
-mAx- will become famous soon enough
1) No. BL have to be attached to the walls
2)Once BL have been created, remesh the volumes to be sure that BL are included in the volume mesh.
Once it is checked (examine mesh), you can export the mesh
3)check the node distribution along the axis of the volume with higher cells number.
The cells distribution in the section is given because of the first volume (annular and conical)
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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



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