|
[Sponsors] |
[blockMesh] Need help defining internal face as boundary |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Member
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6 ![]() |
I'm new to OpenFOAM and I'm trying to create a wedge geometry using Block mesh. I'm facing a couple of issues with the mesh creation. I'm sharing an image showing the boundaries I'll be needing for the simulation.
I'd be really grateful if someone could point me in the right direction. Here's my blockmesh file. Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // scale 1; vertices ( (0 0 0) (0.0105 0 0) (0.0105 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0.0105 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) (0 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) (0.001 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) (0.001 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) (0.001 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0.001 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) (0.0105 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) (0.0105 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0.001 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0.001 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0.001 0 0) (0 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) ); blocks ( //Orifice Block// hex (0 16 9 6 0 16 8 7) (25 10 1) simpleGrading (1 1 1) //Block to the left of orifice (Pair 1)// hex (6 9 10 17 7 8 11 18) (25 90 1) simpleGrading (1 1 1) //Block to the far left of orifice (Pair 2)// hex (17 10 14 3 18 11 15 4) (25 100 1) simpleGrading (1 1.5 1) //Block above orifice (Pair 1)// hex (16 1 13 10 16 1 12 11) (300 100 1) simpleGrading (1.5 1.5 1) //Block above orifice and to the left (Pair 2)// hex (10 13 2 14 11 12 5 15) (300 100 1) simpleGrading (1.5 1.5 1) ); edges ( ); boundary ( upperwall { type wall; faces ( (3 4 15 14) (14 15 5 2) ); } lowerwall { type wall; faces ( (6 7 18 17) (17 18 4 3) ); } pressure_outlet { type patch; faces ( (1 12 13 1) (13 12 5 2) ); } orificewall { type wall; faces ( (6 7 8 9) ); } airinlet { type patch; faces ( (0 7 6 0) ); } axis { type empty; faces ( (0 16 16 0) (16 1 1 16) ); } frontPlane { type wedge; faces ( (0 16 8 7) (16 1 12 11) (11 12 5 15) (7 8 11 18) (18 11 15 4) ); } backPlane { type wedge; faces ( (0 16 9 6) (16 1 13 10) (10 13 2 14) (6 9 10 17) (17 10 14 3) ); } ); mergePatchPairs ( ); // ************************************************************************* // |
|
![]() |
![]() |
![]() |
![]() |
#2 | |
Member
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6 ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 ![]() |
You cannot change two sides of a patch into a boundary. If you really need a boundary in the very inner of your domain, you have to create two patches.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Member
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6 ![]() |
Can you please explain how to do that? Like have duplicate vertices for the same coordinates and then define both sides of the patches as walls (with the different coordinates?)
Last edited by Venky_94; August 10, 2021 at 09:46. |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 ![]() |
> Like have duplicate vertices for the same coordinates and then define both sides of the patches as walls (with the different coordinates?
That is what I recommend.
__________________
Uwe Pilz -- Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950) |
|
![]() |
![]() |
![]() |
![]() |
#6 | |
Member
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6 ![]() |
Quote:
What would you suggest for the second part of my issue though? The interface between the blocks 0,3 and 1,3 is being detected as a wall. How can I change that? |
||
![]() |
![]() |
![]() |
![]() |
#7 | ||
Member
Venkat Ganesh
Join Date: May 2020
Location: Cincinnati, Ohio
Posts: 49
Rep Power: 6 ![]() |
Quote:
Quote:
Code:
/*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v2012 | | \\ / A nd | Website: www.openfoam.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // scale 1; vertices ( (0 0 0) (0.0105 0 0) (0.0105 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0.0105 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) (0 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) (0.001 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) (0.001 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) (0.001 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0.001 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) (0.0105 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) (0.0105 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0.001 0.0099904822158185784680695107340398 -0.00043619387365336002686244665227377) (0.001 0.0099904822158185784680695107340398 0.00043619387365336002686244665227377) (0.001 0 0) (0 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) (0 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) (0 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) //Copy of pt. 6// (0 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) //Copy of pt. 7// (0.001 0.00040011881274353401985455214173726 0.000017469564639817068040427913699908) //Copy of pt. 8// (0.001 0.00040011881274353401985455214173726 -0.000017469564639817068040427913699908) //Copy of pt. 9// (0.001 0 0) //Copy of pt. 16// (0.001 0.0039961928863274315607001518912966 -0.00017447754946134399448194607362694) //Copy of pt. 10// (0.001 0.0039961928863274315607001518912966 0.00017447754946134399448194607362694) //Copy of pt. 11// ); blocks ( //Orifice Block// hex (0 16 22 19 0 16 21 20) (25 10 1) simpleGrading (1 1 1) //Block to the left of orifice// hex (6 9 24 17 7 8 25 18) (25 90 1) simpleGrading (1 1 1) //Block to the far left of orifice// hex (17 10 14 3 18 11 15 4) (25 100 1) simpleGrading (1 1.5 1) //Block above orifice// hex (23 1 13 10 23 1 12 11) (300 100 1) simpleGrading (1.5 1.5 1) //Block above orifice and to the left// hex (10 13 2 14 11 12 5 15) (300 100 1) simpleGrading (1.5 1.5 1) ); edges ( ); boundary ( upperwall { type wall; faces ( (3 4 15 14) (14 15 5 2) ); } lowerwall { type wall; faces ( (6 7 18 17) (17 18 4 3) ); } pressure_outlet { type patch; faces ( (1 12 13 1) (13 12 5 2) ); } orificewall { type wall; faces ( (6 7 8 9) (19 20 21 22) ); } airinlet { type patch; faces ( (0 20 19 0) ); } axis { type empty; faces ( (0 16 16 0) (23 1 1 23) ); } frontPlane { type wedge; faces ( (0 16 21 20) (23 1 12 11) (11 12 5 15) (7 8 25 18) (18 11 15 4) ); } backPlane { type wedge; faces ( (0 16 22 19) (23 1 13 10) (10 13 2 14) (6 9 24 17) (17 10 14 3) ); } mpatch { type patch; faces ( (23 11 10 23) ); } spatch1 //Outlet of orifice// { type patch; faces ( (16 21 22 16) ); } spatch2 //Patch to left of orifice outlet// { type patch; faces ( (8 25 24 9) ); } ); mergePatchPairs ( (mpatch spatch1) (mpatch spatch2) ); // ************************************************************************* // |
|||
![]() |
![]() |
![]() |
Tags |
blockmesh, interface, internal face |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
[snappyHexMesh] Invalid Normals for source face to target face while making AMI? | Sorabh | OpenFOAM Meshing & Mesh Conversion | 1 | August 3, 2021 06:35 |
Error - Solar absorber - Solar Thermal Radiation | MichaelK | CFX | 12 | September 1, 2016 05:15 |
[snappyHexMesh] sHM layer process keeps getting killed | MBttR | OpenFOAM Meshing & Mesh Conversion | 4 | August 15, 2016 03:21 |
[snappyHexMesh] No layers in a small gap | bobburnquist | OpenFOAM Meshing & Mesh Conversion | 6 | August 26, 2015 09:38 |
snappyhexmesh remove blockmesh geometry | philipp1 | OpenFOAM Running, Solving & CFD | 2 | December 12, 2014 10:58 |