extrudeMesh limits Adaptive/Dynamic Mesh Refinenemt
Dear Foamers,
I have an issue using the utility 'extrudeMesh' (foam version 1606+ or 4.1): Using Snappy works fine for dynamic mesh refinement (or adaptive mesh refinement, AMR), like here http://www.openfoam.com/version-v1606+/meshing.php The problem occurs if 'extrudeMesh' is used after using Snappy - all cells with a hanging node will become protected from refinement. I prepared a simple case to demonstrate this problem. This case is based on the motor bike tutorial: $FOAM_TUTORIALS/multiphase/interDyMFoam/ras/motorBike in version 1606 (see link above). It also works with version 4.1. This case can be downloaded here: https://www.dropbox.com/s/3g40075603...rBike.zip?dl=0 Two figures illustrate this problem. The grey colored cells show the 'protectedCells' which cannot be refined. Firstly, this result is gained for the standard tutorial case: https://s13.postimg.org/r9ai68zv7/pr...trude_Mesh.png Secondly, the following figure shows the result after 'extrudeMesh' was used (one cell layer is extruded from the top boundary): https://s29.postimg.org/o49ig5js3/pr...trude_Mesh.png Three scripts lead to different results: - 'runscript_likeTutorial': creates more or less the tutorial case (corresponds to the first figure). - 'runscript_extrude1/2': uses 'extrudeMesh' in two different ways (lead to second figure). Question: How can you use 'extrudeMesh' and dynamic mesh refinement together? Best Regards, Peter Details: 1.) One problem is that the fields cellLevel and pointLevel in the directory constant/polyMesh need to be deleted for 'interDyMFoam'. In 'runscript_extrude1' they are deleted, in 'runscript_extrude2' the cellLevel field is kind of reconstructed with the 'refinementLevel'-utility. Both do not work (-> many protected cells). Remarks: same question asked here but not answered yet unfortunately: https://www.cfd-online.com/Forums/op...tml#post647491 and here: https://www.cfd-online.com/Forums/op...tml#post647510 I also asked this question on the ESI Foam help page https://develop.openfoam.com/Develop...lus/issues/471 Please let me know if you have an idea - I am stuck here... |
Problem solved ?
Did you get any solution to your problem ? I am facing the same problem.
Moreover when using interDyMFoam (and without extrudeMesh), non pure hexahedron cells are not refined (cells at junction between 2 different level). I tried to use the option mergePatchFaces of SHM without any change on the results. The debugging option of dynamicRefineFvmesh gave me : Quote:
Faust |
1 Attachment(s)
Quote:
EDIT : I find the problem. It comes from how the script (./Allrun) is written. I will come back when I would have clearly find out the mistake. |
Dear Paul,
probably this issue has been solved in v1712, see e.g.: https://www.cfd-online.com/Forums/op...tml#post674315 |
All times are GMT -4. The time now is 16:16. |