|
[Sponsors] |
July 20, 2011, 00:24 |
mergePatchPairs Misbehaving
|
#1 |
Senior Member
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 156
Rep Power: 17 |
Hello all,
I am having trouble using blockMesh & mergePatchPairs to generate a C-grid for an airfoil. I was unable to use the grading to get sufficient resolution (y+ = 30), therefore I created a "layer" of blocks around the airfoil to get the proper number of cells and expansion ratio. That led to a problem, which was cells that were too skinny in the wake. I addressed that by using mergePatchPairs to terminate the layer blocks at the trailing edge, however mergePatchPairs appears to join the wake and forward blocks improperly, leading to additional orthogonality errors and poor convergence. I have attached my blockMeshDict and an image of the discrepant cells. What is shown is the block interface along the line that runs vertically from the trailing edge. I would appreciate any assistance at all. Thanks in advance, Dan |
|
July 20, 2011, 00:49 |
|
#2 |
Senior Member
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 156
Rep Power: 17 |
To follow up, I commented out the patches in mergePatchPairs, then ran stitchMesh internalRight internalLeft. I am left with this error:
Code:
--> FOAM FATAL IO ERROR: keyword internalRight is undefined in dictionary "/home/dan/OpenFOAM/dan-1.7.0/run/thesis/cases/VISN017fix2/0/p::boundaryField" file: /home/dan/OpenFOAM/dan-1.7.0/run/thesis/cases/VISN017fix2/0/p::boundaryField from line 27 to line 43. From function dictionary::subDict(const word& keyword) const in file db/dictionary/dictionary.C at line 461. FOAM exiting Thanks, Dan |
|
July 20, 2011, 01:08 |
|
#3 |
Senior Member
Daniel
Join Date: Jul 2009
Location: Montreal, Canada
Posts: 156
Rep Power: 17 |
Last one for now... I got around the stitchMesh issue by temporarily emptying the /0 directory. Now stitchMesh runs, but it produces the same nonOrthogonality issue that blockMesh did.
Another way forward would be to split the wake from two blocks (top and bottom) to four (the thickness of the layer). To make this work, I would have to reduce the number of nodes in the y-direction, meaning it would not match the neighbouring block around the airfoil. I duplicated the vertices at the interfaces in question, and defined the left and right blocks that meet at the interface according to the separate set of duplicate vertices. It looks like the mergePatchPairs is completed successfully on the top, but not on the bottom of the wake (see pic of trailing edge). Any ideas on how I could fix this? New BlockMeshDict attached. Dan Last edited by dancfd; July 20, 2011 at 02:11. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[blockMesh] mergePatchPairs problem when migrating from OFext1.6 to OFext3.1 | ma-tri-x | OpenFOAM Meshing & Mesh Conversion | 2 | January 12, 2021 05:23 |
[blockMesh] mergePatchPairs reducing a face to less than 3 points | aow | OpenFOAM Meshing & Mesh Conversion | 2 | June 1, 2018 18:37 |
[mesh manipulation] mergePatchPairs with multiple areas | bjoern1 | OpenFOAM Meshing & Mesh Conversion | 1 | October 1, 2015 08:07 |
[blockMesh] 2:1 refined 2D blockMeshDict with mergePatchPairs? | chrisb2244 | OpenFOAM Meshing & Mesh Conversion | 0 | December 2, 2014 00:52 |
[Other] mergePatchPairs for a meshing with fixed interface | Cyp | OpenFOAM Meshing & Mesh Conversion | 17 | May 22, 2014 10:34 |