CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[mesh manipulation] Mesh improvement using mesquiteMotionSolver

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 3, 2011, 10:50
Default Mesh improvement using mesquiteMotionSolver
  #1
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 552
Rep Power: 25
philippose will become famous soon enough
Hello and a Good Afternoon,

I use pure tetrahedral meshes for my Simulations. However, there are situations when I need to merge and then stitch two tetrahedral meshes into one complete mesh for simulating cases in which my mesh consists of multiple parts.

A typical example in my case would be the simulation of multiple positions of a hydraulic spool valve. Here, I mesh the valve housing and the spool separately. Then use transformPoints to move the spool to a specific position, followed by mergeMeshes and stitchMesh to create a single mesh for simulation.

When using stitchMesh (which uses the sliding interface features of OpenFOAM internally), checkMesh usually reports many non-orthogonal and highly skewed faces created as a result of the stitching operation.

I was wondering whether it would be possible to use the mesquiteMotionSolver for improving the mesh after running stitchMesh, before starting the simulation?

stitchMesh generates a mesh which contains not just tetrahedra, but also polyhedra, wedges, prisms, etc.... can the mesquiteMotionSolver handle these cell types, and has anyone already tried this?


Have a great weekend ahead!

Philippose

(Theoretically, I could just use mesquiteMotionSolver to move the spool into a new position without creating the mesh as multiple pieces in the first place, but for some reason, mesquiteMotionSolver kept crashing while trying to do this.... hence the more complicated work flow. )
philippose is offline   Reply With Quote

Old   June 3, 2011, 11:19
Default
  #2
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25
deepsterblue will become famous soon enough
Oliver Borm and I worked on hex/prism support for mesquite, and there's a branch with updates on the extend git repository. We haven't worked on polyhedra yet, and I may look into that when I have some time. The stuff should work in parallel as well.

You can play around with the termination criteria to dictate the amount of smoothing that you would need.

Mesh smoothing can only get you so far - at some point, extreme distortion would require topo-changes.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   June 3, 2011, 12:30
Default
  #3
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 552
Rep Power: 25
philippose will become famous soon enough
Hello Sandeep,

A Good Day to you!

Thanks a lot for your reply..... I had seen the Git commits regarding Hex/Prism support for the mesquiteMotionSolver, but was not sure if I could use this in my particular case, since it also includes cell types other than Hex and Prism.

How does the mesquiteMotionSolver handle cells which it cannot modify? Is it possible to run the solver on cases which contain even polyhedral cells, or will it abort with errors?


Out of interest..... how did Oliver and you go about handling these other cell types? Do you internally split the cells temporarily into tets and then apply the mesquite functions to them, or does the mesquite library inherently in some form or the other already support arbitrary cell types?

Have a nice day ahead!

Philippose
philippose is offline   Reply With Quote

Old   June 3, 2011, 12:35
Default
  #4
Senior Member
 
Sandeep Menon
Join Date: Mar 2009
Location: Amherst, MA
Posts: 403
Rep Power: 25
deepsterblue will become famous soon enough
Unfortunately, mesquite only supports standard cell types at the moment. Polyhedra are tricky, since there's no clear way to define a good "cell-quality" for them. My attempt is to decompose into tets, which should get the job done in theory, but there's always robustness/efficiency concerns.

If you have cell-types that don't match a standard type, it will quit.
__________________
Sandeep Menon
University of Massachusetts Amherst
https://github.com/smenon
deepsterblue is offline   Reply With Quote

Old   June 3, 2011, 12:43
Default
  #5
Senior Member
 
Philippose Rajan
Join Date: Mar 2009
Location: Germany
Posts: 552
Rep Power: 25
philippose will become famous soon enough
Hi again,

Thanks once again for the information.....

Hmmm... that would be a problem in my case. I cannot control the kind of cells stitchMesh creates, and I have seen that it does create the whole range of cell-types (as reported by checkMesh).

I guess I will have to try to use the mesquiteMotionSolver directly, by creating one tet mesh with the housing and the spool, and trying to move the spool by "x" units using moveDynamicMesh before starting a simulation.

I may have to get back to you regarding this, because the last time I tried, mesquiteMotionSolver quit after a couple of steps with negative volume errors (if I remember right).

Philippose
philippose is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Add Mesh Layers doesnt work on the whole surface Kryo OpenFOAM Meshing & Mesh Conversion 13 February 17, 2022 08:34
[ICEM] surface mesh merging problem everest ANSYS Meshing & Geometry 44 April 14, 2016 07:41
[snappyHexMesh] SnappyHexMesh no layers and no decent mesh for complex geometry pizzaspinate OpenFOAM Meshing & Mesh Conversion 1 February 25, 2015 08:05
Fine/Marine mesh improvement issue inarus Fidelity CFD 1 February 16, 2015 09:03
[snappyHexMesh] snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Meshing & Mesh Conversion 2 March 27, 2011 22:11


All times are GMT -4. The time now is 05:43.