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

[ICEM] Meshing from porous media out of CT data

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 23, 2010, 13:02
Default Meshing from porous media out of CT data
  #1
New Member
 
Join Date: Apr 2010
Posts: 17
Rep Power: 6
shohin is on a distinguished road
Hi all,
i have some delicate problem to solve. The attached pictures show the fluid fraction of a porous media which was scanned by a computer tomograph. The dimension are: Lenght 8cm and diameter 1cm.

As you can guess, my task is to prepare a mesh for CFD simulations.

My procedure until now was the following:
- I segmented the fluid part from the solid and exported it as STL.
- I imported the STL in 3matic and made some fixing procedures (overlapping elements, inverted normals, etc..) and some first remeshing of the surfaces.
- Afterwards i defined the inlet and outlet surfaces
- Sinde 3matic is only available for me as 32bit version i cannot generate volume meshes
- So i exported the surface mesh (2Mio elemnts) as NASTRAN using surface splitting to have the bocos already defined in ICEM
- Next import in ICEM
Now, the problems are arising:
- I imported as Geometry => correct or is import as Mesh better?
- remeshing the surface mesh using "patch independent" and curvature/proximity based refinement
- a mesh check, produced lots of errors, which i cannot repair easily.

So, since i am new to ICEM i wanted to ask for some hints or tips about the best procedure to solve this task.
Also, i am interested in your opinions about the overall feasiblity.

Thanks for your comments in advance.

all the best,
flo
Attached Images
File Type: jpg porous_media_fluid2.jpg (46.5 KB, 63 views)
File Type: jpg porous_media_fluid3.jpg (51.2 KB, 59 views)
File Type: jpg porous_media_fluid.jpg (11.8 KB, 49 views)
shohin is offline   Reply With Quote

Old   July 24, 2010, 10:09
Default Can you start from the imported surface mesh?
  #2
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,651
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
This probably has all sorts of separate regions... Very interesting.

Is the Nastran mesh good enough quality to use as a surface mesh? (it looks good in the pics)

If so, you could read it in as mesh or geometry and then use edit (menu) => Mesh to facets or facets to mesh to get the other. It is good to have both so the smoother can work better.

Then go to Mesh (tab) => Compute Mesh and choose Delaunay (with the TGlib - AFoption if you have the memory for it).

This will generate the volume mesh from your surface mesh...

If that won't work for you, we can talk about how to get a mesh with the Patch independent mesh (hint, you don't need a surface mesh first)...
PSYMN is offline   Reply With Quote

Old   July 24, 2010, 10:27
Default
  #3
New Member
 
Join Date: Apr 2010
Posts: 17
Rep Power: 6
shohin is on a distinguished road
Hi PSYMN,

first, thanks for your reply.

1) Yes, i think the NASTRAN mesh is suitable for further work

So, i will give it a try. I already played around with Delauney and created volume meshes. Unfortunately, i got errors when reading it in into FLUENT.
But, i will give it a try with the method you proposed.

I will report the results and are looking forward for your comments.

best,
flo
shohin is offline   Reply With Quote

Old   July 28, 2010, 10:58
Default
  #4
New Member
 
Join Date: Apr 2010
Posts: 17
Rep Power: 6
shohin is on a distinguished road
hey,

i used the delauney mesher with your options proposed and got a volume mesh. the mesh check in FLUENT brought me ~700 degenerated contact points as only errors. A first test simluation was promising though.

Since i am doing my preliminary test only with a small part of the whole model,
i will not be able to use 3matic for the whole model.
Here i will have to use ICEM for importing the "unmeshed" NASTRAN (with bocos) and further meshing. I also have to do "facets to mesh", right before surface meshing? I already made a surface mesh using "patch independent" algorithm which seems quite nice.
For creating the volume mesh i used delauney with the proposed settings but failed due to lack of memory (i have 16GB available at the moment). When playing around with the max element size i always got the same amount of volume elements when using delauney - what is the reason for this.

Is there any memory saving technique or way to use?

Thanks a lot for your kind help.

best
flo
shohin is offline   Reply With Quote

Old   July 29, 2010, 10:42
Default Answers...
  #5
Senior Member
 
PSYMN's Avatar
 
Simon Pereira
Join Date: Mar 2009
Location: Ann Arbor, MI
Posts: 2,651
Blog Entries: 1
Rep Power: 34
PSYMN has a spectacular aura aboutPSYMN has a spectacular aura about
Shohin,

Instead of generating the surface mesh with patch independent, you can skip that step and just generate the volume mesh with the Tetra/Mixed => Octree method. This will give you the PI surface and Octree volume mesh in one go and using the same memory and time (pretty much) that it took just to generate the PI surface mesh.

The 700 degenerated contact points are just "non-manifold verts". You can check for those in ICEM CFD. You could repair them or at least some of them. Depending on why they formed, you may even be able to avoid them with a smaller size or with "thin cuts" or a smaller "Edge Criterion", but I think it is reasonable to expect these sorts of things in this sort of geometry.

The "Facets to Mesh" step is only if you wanted to use those facets as a surface mesh. if you just want to patch independent mesh the geometry, you can skip it. If you had a lot of sharp features you wanted to capture then I would have suggested building diagnostic topology to create feature curves and break up the model. But since you probably just have a few or none, I would suggest drawing your features on the model with spline curves thru faceted points... I can discuss that more if necessary.

When using Delaunay, the mesh starts with the surface size and grows at the growth rate toward the max size. There is no guarantee that it will reach the max size. It may have not reached the max size in the first case, so increasing the max size even further would not have made any difference.

It is possible to save more memory by running in batch... How big of a mesh are you trying to generate? On what hardware?

Simon
PSYMN is offline   Reply With Quote

Old   July 29, 2010, 14:03
Default
  #6
New Member
 
Join Date: Apr 2010
Posts: 17
Rep Power: 6
shohin is on a distinguished road
Hi PSYMN,

thanks for your answer.

i generated a mesh w/o a surface mesh first using the octree method and ended up with ~22Mio elements. It was imported in fluent w/o error. But, a mesh check revealed 40 cells having only wall faces.
I will do mesh check and some smoothing if required in ICEM?

It would be very nice if you could comment a little bit on the method with the "spline curves thru faceted points" you talked about.

If you do not mind, i would ask two more questions:

1) Consider i did not have the bocos defined already in the NASTRAN file and i would like to import the raw STL directly into ICEM. Since both the inlet and outlet are plane surfaces, is there an easy method to select them and throw them in a new part?

2) Since it seems that my memory is too low at the moment to do the calculations, how could i split up the mesh (lets say 2 or 3 parts) and prepare them for FLUENT?


Again, i really want to thank you, since your work here in this forum is exceptional and especially the newbies like me learn a lot!!

ps. i work on a dual-quad xeon workstation with 16GB RAM at the moment (will upgrade to 24GB in the next days )
best
flo
shohin 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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Porous Media coupled with internal flow Samuel Andrade FLUENT 2 August 26, 2012 10:43
UDS porous media Nagendra FLUENT 10 February 1, 2011 21:57
combustion in porous media burner zuby ANSYS 0 July 27, 2009 03:26
DPM on porous media magnounibo FLUENT 0 April 23, 2009 07:35
Porous Flow and Data CFX 4.4 Robin CFX 0 April 14, 2005 10:56


All times are GMT -4. The time now is 17:23.