SnappyHexMesh..
Dear All,
I have a question: is it possible to generate a volume mesh inside a certain surface? I mean, I have a non-simple stl geometry (the reconstruction of nasal cavity) and I wanna mesh the inside volume. If my file is called nasal_cavity.stl, how can I set a good snappyHexMeshDict file in order to generate my mesh inside the surface?? It is a closed surface. The only problem is that is presents some loops. I hope this is not a problem. Here is my snappyHexMeshDict file: Code:
geometry Code:
samiam@penrose: ~/Desktop/naso_ok/test > snappyHexMesh snappyHexMeshDict Samuele |
volume mesh inside a region
Hi samiam,
Did you use the locationInMesh to specify the portion of your geometry you wanted to keep? You can find the details about its usage in the User guide. |
What does this mean? Also, could you link me the SnappyHexMesh user guide, please??
Do you mean the OpenFOAM u'sg?? Thanks a lot, Samuele |
Hi
two things: 1) as mentioned in the previous post the location in mesh entry is a point in the part of the mesh that will be kept 2) the correct command is >this/is/your/case/directory>snappyHexMesh and maybe one of the arguments specified in your failure message. You don't have to tell sHM that your are using a Dict-file. The stl you probably have here: >this/is/your/case/directory/constant/triSurface and the sHMD file here >this/is/your/case/directory/system then everything should be alright hope that helps regards Colin Edit: yes he is refering to the OF userguide, where sHM is explained Edit2: Have you heard about the motorbike tutorial? maybe you should try that first and then adapt your case to that one. |
Dear Collin,
thanks for your answer. Well, just 1 more question: 1. where should I add the line with a point that I wanna include in the meshed region? Thanks. Samuele |
in your snappyHexMeshDict file line 160
locationInMesh (X Y Z); this point must be in the area you want to keep (in your pipe what ever) and must not coincidence with a face! |
I think I did a mistake then..
My snappyHexMeshDict file is no bigger than 20 lines.. Code:
geometry |
Oh that is a problem.
I assumed its not bigger than 20 lines but I really hoped that you just quoted parts of your sHMD. Well then I highly recommend like I did it in one of my previous posts the motorbike- tutorial and the user guide. the motorbike tutorial you find here: openFoam171/tutorials/incompressible/simpleFoam/motorBike have a look in the system folder there you find an example file for the sHMD the user guide you find here good look with it and if you have further questions just ask regards |
I'll work hard on it.
Thanks a lot. Just one more question before beginning: in the motorbike tutorial (the Dict file!) I can find different "zones" to describe the surface.. But I do not need them. Can I delete that part, right?? |
yes you can delete them
|
Samuele, you might want to read up on sHM in the user's manual. There you'll find a nicely commented simple case.
Quote:
|
Dear All,
it's again from me. I read the tutorials and I followed your advice. However, there is a point I can not find anywhere. I have a complicated geometry and I have its internal mesh. The point is that I don't know how to set which is the inlet and which is the outlet. Of course, I can not write the cell numbers, manually. Is there an automatic tool to find the inlet and the outlet section(s)?? Thanks, Samuele |
Hi
ususally this is done by the blockMesh background mesh. Otherwise you have to specify different regions in your stl file so you get different patches and then assign in the files listed in the 0 folder different BC eg inlet outlet to these patches. check for this the different layers of the motorBike tutorial, but in this case to all of the layers the BC describing them as a wall are assigned regards Colin |
Dear Colin,
thanks for answering. I will try to have a look. I guess that your second solution (layers!) is the best one. Bu I am not sure to be able to define different regions in the stl files.. But I will try to do my best. Thanks again, Samuele |
All times are GMT -4. The time now is 17:22. |