|
[Sponsors] |
May 28, 2012, 05:41 |
Discontinuity in Paraview
|
#1 | ||
New Member
Jonathan
Join Date: May 2012
Posts: 4
Rep Power: 13 |
Hi everyone,
First post, i hope you can help me I'm french so, sorry if i make too big mistakes^^ So, i have a discontinuity in my simulation (using twoLiquidMixingFoam) when i try to simulate an avalanche interacting with an obstacle. You can see that with attached files (bug, bug2, bug3). I hope you can tell me how to fix that. I'm using OpenFOAM 1.6. I also join blockMeshDict and boundaries. Quote:
Quote:
Jonathan Last edited by Mkoll; May 28, 2012 at 05:57. |
|||
May 28, 2012, 07:29 |
|
#2 |
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 255
Rep Power: 21 |
Hi Jonathan,
you should refine your mesh around the obstacle. Have a look at the different cell sizes there ("Surface with edges" visualization in Paraview). You can do this by using simpleGrading parameters in the blocks definition part like this: Code:
FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 1; vertices ( (0 0 0) (20 0 0) (20.4 0 0) (38 0 0) (0 0.6 0) (20 0.6 0) (20.4 0.6 0) (38 0.6 0) (0 10 0) (20 10 0) (20.4 10 0) (38 10 0) (0 0 0.1) (20 0 0.1) (20.4 0 0.1) (38 0 0.1) (0 0.6 0.1) (20 0.6 0.1) (20.4 0.6 0.1) (38 0.6 0.1) (0 10 0.1) (20 10 0.1) (20.4 10 0.1) (38 10 0.1) ); blocks ( hex (0 1 5 4 12 13 17 16) (125 10 1) simpleGrading (0.05 1 1) hex (2 3 7 6 14 15 19 18) (125 10 1) simpleGrading (20 1 1) hex (4 5 9 8 16 17 21 20) (125 50 1) simpleGrading (0.05 1 1) hex (5 6 10 9 17 18 22 21) (25 50 1) simpleGrading (1 1 1) hex (6 7 11 10 18 19 23 22) (125 50 1) simpleGrading (20 1 1) ); edges ( ); patches ( wall leftWall ( (0 12 16 4) (4 16 20 8) ) wall rightWall ( (7 19 15 3) (11 23 19 7) (8 20 21 9) (9 10 22 21) (10 11 23 22) ) wall lowerWall ( (0 1 13 12) (1 5 17 13) (5 6 18 17) (2 14 18 6) (2 3 15 14) ) //(8 20 21 9) //(9 10 22 21) //(10 11 23 22) ); mergePatchPairs ( ); |
|
May 28, 2012, 08:22 |
|
#3 |
New Member
Jonathan
Join Date: May 2012
Posts: 4
Rep Power: 13 |
Hi Martin
Thanks for answering. I tried with your blockMesh It seems to be a little better ? I attached some screenshots. Jonathan |
|
May 28, 2012, 16:55 |
|
#4 |
Senior Member
|
As already pointed out by Martin, it seems like you have a way too big discontinuity in your mesh. I would suggest you to post a picture that depicts your mesh around the obstacle.
|
|
May 28, 2012, 17:38 |
|
#5 |
New Member
Jonathan
Join Date: May 2012
Posts: 4
Rep Power: 13 |
Hi lovecraft22
Thanks for answering. I hope it's what you ask. |
|
May 28, 2012, 17:44 |
|
#6 |
Senior Member
|
Yes it is. As you can see you have a big difference in your cells, as a matter of fact it the picture you attached the cells above the obstacle cannot be seen as the blue lines are too close to each other. So you can either make a coarser mesh in that region or make a finer mesh in the surrounding regions.
|
|
May 29, 2012, 02:32 |
|
#7 |
New Member
Jonathan
Join Date: May 2012
Posts: 4
Rep Power: 13 |
Hi lovecraft22
The goal was to make finer the mesh around the obstacle. I'll try with the same precision everywhere. Thanks, Jonathan |
|
May 29, 2012, 04:11 |
|
#8 |
Senior Member
|
It can (and should) be finer around the obstacle but the transition between one mesh level and the one close to it should be smooth.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[General] Paraview nice animation | PyGloo | ParaView | 4 | June 7, 2012 13:34 |
Newbie: Install ParaView 3.81 on OF-1.6-ext/OpenSuse 11.2? | lentschi | OpenFOAM Installation | 1 | March 9, 2011 03:32 |
Paraview not found | fusij | OpenFOAM Installation | 2 | January 1, 2011 21:44 |
[OpenFOAM] Distributed ParaView and PV3FoamReader | micalil | ParaView | 4 | July 1, 2010 06:09 |
paraFoam reader for OpenFOAM 1.6 | smart | OpenFOAM Installation | 13 | November 16, 2009 22:41 |