|
[Sponsors] |
April 9, 2014, 10:14 |
Periodic BC rotating channel
|
#1 |
Senior Member
Jiri
Join Date: Mar 2014
Posts: 221
Rep Power: 13 |
Hi there,
I have ascii mesh .msh of rotor segment and another .msh of draft tube. Desired BC conditions are velocity at the inlet of the rotor and static pressure at the draft tube outlet and mixing plane between rotor and drafut tube. I merge these meshes by using mergeMeshes but I cant figure out the correct setting of periodical conditions for the rotor faces (rotor mesh is channel with one blade - but there are five blades really). Moreover I d like to know the settings for mixing plane between this rotor segment and the draft tube. I am using OF 2.3.0 I struggle a lot but cannot sort it out. I am grateful for any help! |
|
April 11, 2014, 06:15 |
|
#2 |
Senior Member
Jiri
Join Date: Mar 2014
Posts: 221
Rep Power: 13 |
I would like to prepare calculation in the following way:
- Having two separate mesh files (e.g. rotor.msh and draft_tube.msh) - by using fluentMeshToFoam rotor.msh -writeSets a convert .msh to OF format. The same I do with draft_tube.msh. - then merge them in OF (as I said I use mergeMeshes as I do not want to merge them in another program) - set boundary conditions, but still I cannot find out how. Thank you for any advice |
|
April 15, 2014, 04:45 |
|
#3 |
Senior Member
Jiri
Join Date: Mar 2014
Posts: 221
Rep Power: 13 |
I successfully merged two regions -rotor and draft tube, however I am struggling with BC.
Rotor domain is just a rotating segment (just 90° to save mesh size). Rotor domain is composed of: inflow outflow blade interface on hub (gap between hub and blade) interface on shroud(gap between shroud and blade) per1 - first side face of Rotor domain per2 - second side face of Rotor domain Draft tube domain (is stationary) is composed of: wall-4 ...inlet (outlet from rotor domain) wall-5 ...outlet (pressure outlet) wall-3, wall-6, wall-7 are walls. Could anyone help me please with setting BC for rotating Rotor domain? I guess the main trouble is in setting per1 and per2. I use OpenFoam 2.3.0 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.3.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class polyBoundaryMesh; location "polyMesh"; object boundary; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 ( hub { type wall; inGroups 1(wall); nFaces 3825; startFace 2246933; } shroud { type wall; inGroups 1(wall); nFaces 3098; startFace 2250758; } inflow { type patch; nFaces 3128; startFace 2253856; } outflow { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch wall-4; zone outflow_ZONE; transform noOrdering; nFaces 3332; startFace 2256984; } blade { type wall; inGroups 1(wall); nFaces 6254; startFace 2260316; } hub-tip-ggi-side-1 { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch hub-tip-ggi-side-2; zone hub-tip-ggi-side-1_ZONE; transform noOrdering; nFaces 816; startFace 2266570; } hub-tip-ggi-side-2 { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch hub-tip-ggi-side-1; zone hub-tip-ggi-side-2_ZONE; transform noOrdering; nFaces 876; startFace 2267386; } shroud-tip-ggi-side-1 { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch shroud-tip-ggi-side-2; zone shroud-tip-ggi-side-1_ZONE; transform noOrdering; nFaces 930; startFace 2268262; } shroud-tip-ggi-side-2 { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch shroud-tip-ggi-side-1; zone shroud-tip-ggi-side-2_ZONE; transform noOrdering; nFaces 915; startFace 2269192; } per1 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 3468; startFace 2270107; transform rotational; neighbourPatch per2; rotationAxis (0 0 1); rotationCentre (0 0 0); rotationAngle 90; // degrees } per2 { type cyclicAMI; inGroups 1(cyclicAMI); nFaces 3468; startFace 2273575; transform rotational; neighbourPatch per1; rotationAxis (0 0 1); rotationCentre (0 0 0); rotationAngle 90; // degrees } wall-3 { type wall; inGroups 1(wall); nFaces 3272; startFace 2277043; } wall-4 { type cyclicAMI; inGroups 1(cyclicAMI); neighbourPatch outflow; zone wall-4_ZONE; transform noOrdering; nFaces 5580; startFace 2280315; } wall-5 { type patch; nFaces 2390; startFace 2285895; } wall-6 { type wall; inGroups 1(wall); nFaces 6118; startFace 2288285; } wall-7 { type wall; inGroups 1(wall); nFaces 22156; startFace 2294403; } ) // ************************************************** *********************** // |
|
April 23, 2014, 06:00 |
|
#4 |
Senior Member
Jiri
Join Date: Mar 2014
Posts: 221
Rep Power: 13 |
any idea please?
|
|
April 25, 2014, 08:11 |
|
#5 |
Senior Member
Jiri
Join Date: Mar 2014
Posts: 221
Rep Power: 13 |
Please have a look at my files, I uploaded it here:
http://leteckaposta.cz/632363549 I would like to set periodic boundaries on per1 and per2 which are the side faces of the blade segment. When I set per1 and per2 as cyclic I cannot visualise the mesh in openfoam. |
|
May 23, 2014, 14:16 |
|
#6 | |
Member
Sami
Join Date: Nov 2012
Location: Cap Town, South Africa
Posts: 87
Rep Power: 14 |
Hi dear Jiricbeng,
To set periodic BC, you need to have the same mesh in the two faces. You can mesh a face and then copy to another before setting the cyclic BC. Mehrez Quote:
|
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
3D Rotating Channel | Curtin University | FLUENT | 5 | November 9, 2018 03:25 |
question on simulating rotating channel with a pressure gradient | huangxianbei | CFX | 10 | January 15, 2014 00:00 |
Periodic channel flow with time dependent mass flow rate | QBeast | FLUENT | 3 | May 10, 2013 14:14 |
rotating channel | emilio | Siemens | 1 | April 6, 2012 01:21 |
Setting up for Rotating 2D straight channel | pkgupta | FLUENT | 0 | April 2, 2010 06:56 |