CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Meshing & Mesh Conversion (https://www.cfd-online.com/Forums/openfoam-meshing/)
-   -   [mesh manipulation] DynamicMesh - Move internal points (https://www.cfd-online.com/Forums/openfoam-meshing/141311-dynamicmesh-move-internal-points.html)

nucat September 4, 2014 15:29

DynamicMesh - Move internal points
 
Hallo foamers,

actually I am programming an y+ adaption for OpenFoam. I need a tool, which manipulates the mesh by making the thickness of the first boundary layer smaller if there is detected y+>1.

I am able to detect that vertices with toposet:
  • select all faces with y+>1
  • select all cells with that selected faces
  • select points of the selected cells (first boundary layer)
  • deselect points of the boundary face
The result is a selection of points of the first boundary cells, but not points directly of the boundary (they have to stay fixed). Now I want to move them down and want to do it with dynamicMesh library.


Problem: I didn't find in the web a tutorial, how to move internal points (not boundary points) of the mesh. Do anybody knows how it works?


Thanks for help,
Nucat

mturcios777 September 4, 2014 18:35

If you look at the layeredEngineMesh code, you can see there is code that moves certain mesh points to a new position based on the calculated crank angle. You will need to do something similar, although I don't quite know how you would determine the new position.

It would have to be some kind of function that finds distance to the nearest point on a boundary patch, then takes a fraction of the distance along the vector to the nearest point on the patch and moves the new point to that position (like half the distance?). Or just a set distance along that line. There are some built in functions like that in meshTools and other libraries. Its been a while since I've had to do anything with them.

nucat September 6, 2014 08:54

Hey,

thanks for reply. I think I will solve that problem by myself with toposet and python ))


All times are GMT -4. The time now is 06:23.