|
[Sponsors] |
[ICEM] Negative volume error while meshing concentric spheres |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
December 21, 2020, 13:20 |
Negative volume error while meshing concentric spheres
|
#1 |
Member
Join Date: Oct 2012
Posts: 46
Rep Power: 14 |
Hello! I need help with meshing in ICEM after banging my head against the wall for long.
I am trying to mesh two concentric spheres in ICEM using blocking. I used two methods: First: I created a 3D sphere, by making a quarter-circle and then doing the steps (2d-Blocking, mirror, Extrusion by rotate) shown in images titled Method1_step Second: I tweaked the steps of the sphere-cube tutorial to mesh the two spheres. I don't do any association, except I snap vertices for the 3D blocks for outer and inner spheres. Also, I create an O-grid around the inner-sphere. The images have been attached as Method2. In both methods, I encounter the negative volume error. For method 2, the error is shown in the 3x3x3 determinant quality. When I try to smoothen pre-mesh using the Pre-Mesh Smooth option, and although it smoothens the quality of determinant to above 0, the change is not reflected in the pre-mesh and if I do any other operation on the mesh, it reverts to previous quality. Can someone please help me with this issue? I can try some other method, if anyone can suggest. Thanks ! |
|
December 21, 2020, 13:25 |
|
#2 |
Member
Join Date: Oct 2012
Posts: 46
Rep Power: 14 |
These are the images for method 2 I've used. I can post the project files/, more images and explanation if required. The unstruct mesh check gives NO error/possible problems, and yet the solver tells me that there are negative volumes in mesh. This is one of the worst problems I've faced in ICEM meshing, please help!
|
|
December 21, 2020, 15:06 |
|
#3 |
Senior Member
Sebastian Engel
Join Date: Jun 2011
Location: Germany
Posts: 567
Rep Power: 21 |
Hi samurai_01,
Have a look at the blocks whose corners show low quality elements. Just the respective block. You might see that that corner with the bad element is almost squashed to a flat plate... Imagine a dice. Look at one corner where you can see all three attached edges. Now try to imagine how the cube/dice would have to look when all three edges are found in same flat/2d plane. What does this mean for the other faces? And what volume would this object have? Now, when we investigate the surface of the sphere, we can assume that it is locally flat, since the curvature is low at a large radius. Another observation is, that you created an o-grid around the inner sphere. But now, the block containing that o-grid is a cuboid again. And that cuboid is now inside a sphere again. This leaves you with a case which is similar but in a way inverse to the o-grid around the inner sphere. This requires the same solution strategy, an ogrid. In method 1 you correctly extended the o-layer to the outer sphere, but method 2 does not have this merit Bottom line is, you need to fix your blocking in a way that keeps every block as close as possible to a perfectly rectangular cuboid. I see two solutions straight forward solutions. The first, create an ogrid with a single layer just as in the method1 where the single layer spans across the radius from the inner to the outer sphere. The advantage is, the block structure is very simple. The disadvantage is the huge difference in mesh density between regions of small radii and regions of large radii. The second approach is to use the blocking you have from method2 and just introduce another o-grid all around. That way you can reshape the degenerated blocks at those 'corners'. I imagine just selecting all blocks and hit create in the o-grid feature should do the trick. Though, you might need to add the block in the void of the inner sphere, or select the boundary faces on the inner sphere. The advantage is you can continue with what you have. Moreover, the mesh density distribution is narrower. The disadvantage is, you will have to do some tedious manual vertex moving. Best, Sebastian |
|
December 22, 2020, 01:04 |
|
#4 | |
Member
Join Date: Oct 2012
Posts: 46
Rep Power: 14 |
Dear Sabastian,
Thank you for the reply, let me just quickly attempt these and get back to you. I Quote:
|
||
December 28, 2020, 04:35 |
Update, Problem solved!
|
#5 |
Member
Join Date: Oct 2012
Posts: 46
Rep Power: 14 |
Dear Sebastian,
It worked!! I got rid of all the negative volumes by method2 you suggested. I have attached the image of the quality check and a cut plane! Thanks a ton!! I can't thank you enough |
|
Tags |
blocking, icem, mesh, negative, volume |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to use "translation" in solidBodyMotionFunction in OpenFOAM | rupesh_w | OpenFOAM Running, Solving & CFD | 5 | August 16, 2016 05:27 |
blockMesh error ... | balkrishna | OpenFOAM Pre-Processing | 0 | August 17, 2010 03:39 |
[blockMesh] BlockMesh FOAM warning | gaottino | OpenFOAM Meshing & Mesh Conversion | 7 | July 19, 2010 15:11 |
On the damBreak4phaseFine cases | paean | OpenFOAM Running, Solving & CFD | 0 | November 14, 2008 22:14 |
[blockMesh] Axisymmetrical mesh | Rasmus Gjesing (Gjesing) | OpenFOAM Meshing & Mesh Conversion | 10 | April 2, 2007 15:00 |