2 Attachment(s)
Thank you Bastian and Bruno for your quick replies to my problem!
@Bastian: I am not completely sure what you mean by Quote:
@Bruno: Just out of curiosity would checking the .eMesh file give me the same information or is there a difference one should have in mind between .obj and eMesh features? According to your suggestion I checked the .obj files which seems to be right. But not all of the lines lay inside my blockMesh block. Actually this could be the reason for the bad edges I get. I think I have to change my procedure as I used to have stl geoemtries overlapping my blockmesh (see images of my first post). Now I seem to need the stl geometry beeing completely surrounded by a blockmesh. In other words all edges have to be inside the blockmesh. Can you confirm that according to your experience? Today I only had time for a quick test which actually showed some improvements. Tomorrow I will have some more time for further test and will let you know about my progress. Kind regards Sebastian Attachment 10417 Attachment 10418 |
Quote:
You need to generate feature lines at the intersections in your emesh-File. Regards Bastian |
1 Attachment(s)
This whole topic seems to be more complicated than expected.
Bastian, I read your answer several times and thought about it. But I'm still not completely sure if I got you right.The white lines in the picture below are my feature lines (eMesh-File). They intersect the block. I think this is what you mean by: Quote:
I tried some workarounds for this problem: 1.) Having the geometry completely surrounded by the blockMesh box. As front and back are now inside the box, there are now intersection lines. But this didn't solve the problem. 2.) Paradoxically if only the front of the prism lies inside whereas the back overlaps the box I get relatively nice results. Unfortunately I don't understand why as I would expect my first workaround to be better than the second. As the second workaround doesn't make any sense to me, although I now get a nice prism, I still don't have the impression that I really solved the problem. Was that a lucky strike? cheers Sebastian Attachment 10451 |
Quote:
Regards Bastian |
Hi all,
sorry that I havn't replied until now. I didn't have the time to follow up on this until this week. @bastil: thank you very much for your input and tips on this!! I tried with several different modifications of the featureEdges - unfortunately without success. For all who might have some similar issues: unfortunately I couldn't find a clean way to solve my problem. But I now have kind of a smooth solution with which I can deal. I used the described workarround - one side of the geometry overlapping and one side inside the blockMesh block. Even after hours of trying and thinking about, I really don't now why this is working. My best guess is that the featureEdge funtion is not fully developped until now. With the workaround I used one problem remains: at the edges my wall patch (prism sides) owns one line of cells of the inlet patch. But this I could solve with a combination of autoPatch and createPatch. Tip: sometimes it helps to run snappyHexMesh a second time. Thanks for your help! Sebastian |
Quote:
I tried to run the command : surfaceFeatureExtract -includedAngle 150 icube.stl icube but the error showed up such as : --> FOAM FATAL ERROR: Cannnot read "icube.stl" From function triSurface::read(const fileName&, const word&, const bool) in file triSurface/triSurface.C at line 370. Could you please tell me how to solve this? my icube.stl file was created by ProE, if i just perform SHM, it worked fine. But why couldn't I extract the feature? Thank you so much! Aqua |
Hey aqua
Are you sure your STL is located in the right directory? It should be in constant/triSurface |
Quote:
Thank you so much for your reply. But yes, I have the stl file under triSurface. Actually I successfully creat mesh by snappyHexMesh, which should means the stl file is right. But just couldnt perform surfaceFeatureExtract... Regards! Yebo |
But in this case, shouldn't you also adjust the directory in the surfaceFeatureExtract command as well?
For example, I use surfaceFeatureExtract -includedAngle 150 -writeObj constant/triSurface/myGeometry.stl myGeometry |
Greetings to all!
Quote:
@Aqua: you had quoted my post #2... on post #4 was the answer to your question!! ;) Best regards, Bruno |
Quote:
Yes, I got it! Thank you so much! Cheers! Aqua |
Quote:
You are right! I forgot the directory.... Thank you so much for your help~ Cheers~ Aqua |
Hey, Bruno,
I am using foamToCGNS, the mesh is created in OpenFoam by SnappyHexMesh, but when run foamToCGNS, error happens like : Wrong number of vertices in cell expected 4,5,6, or 8, found 0 From function foamToCGNS in file writeCGNS.H at line 219. Does CGNS deal with polyhedral cell created by snappyHexMesh? If not, is there some other way to convert OF mesh to CFX? Thank you so much! Aqua |
Hi Aqua,
Quote:
Code:
checkMesh Bruno |
Quote:
Thank you so much for your reply! I did checkMesh, but the solution is mesh ok! " Min volume = 0.0373576. Max volume = 936694. Total volume = 1.02747e+10. Cell volumes OK. Mesh non-orthogonality Max: 54.4594 average: 10.6017 Threshold = 70 Non-orthogonality check OK. Face pyramids OK. Max skewness = 1.18109 OK. Mesh OK. End" Do you know some way to convert the polyhedral cell into hex or tetra? as long as no polyhedral is fine for me.. Thank you so much! Aqua |
Quote:
|
Quote:
That was me, too... but no one answers me, so I asked many times... I will check netgen, thank you so much for your suggestion! Aqua |
Hi Aqua,
Quote:
And there seems to be some hope: Quote:
Bruno |
Hi there,
Can anyone please guide me on how to control the mesh density at the tip of the cylinder as show in the post picture? Which parameter and how to realize? Thanks in advance. Quote:
|
I just set 2 levels for the mesh of the cylinder in snappy and then I used feature edge to have a sharp edge. Read the whole discussion, I think you should be able to sort that out and if not just get back and explain where you're having a hard time.
|
All times are GMT -4. The time now is 10:12. |