CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

Periodic BC rotating channel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 9, 2014, 10:14
Default Periodic BC rotating channel
  #1
Member
 
Jiri
Join Date: Mar 2014
Location: Czech Republic
Posts: 38
Rep Power: 2
Jiricbeng is on a distinguished road
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!
Jiricbeng is offline   Reply With Quote

Old   April 11, 2014, 06:15
Default
  #2
Member
 
Jiri
Join Date: Mar 2014
Location: Czech Republic
Posts: 38
Rep Power: 2
Jiricbeng is on a distinguished road
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
Jiricbeng is offline   Reply With Quote

Old   April 15, 2014, 04:45
Default
  #3
Member
 
Jiri
Join Date: Mar 2014
Location: Czech Republic
Posts: 38
Rep Power: 2
Jiricbeng is on a distinguished road
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;
}
)

// ************************************************** *********************** //
Jiricbeng is offline   Reply With Quote

Old   April 23, 2014, 06:00
Default
  #4
Member
 
Jiri
Join Date: Mar 2014
Location: Czech Republic
Posts: 38
Rep Power: 2
Jiricbeng is on a distinguished road
any idea please?
Jiricbeng is offline   Reply With Quote

Old   April 25, 2014, 08:11
Default
  #5
Member
 
Jiri
Join Date: Mar 2014
Location: Czech Republic
Posts: 38
Rep Power: 2
Jiricbeng is on a distinguished road
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.
Jiricbeng is offline   Reply With Quote

Old   May 23, 2014, 14:16
Arrow
  #6
Member
 
Mehrez
Join Date: Nov 2012
Location: Bordeaux, France
Posts: 43
Rep Power: 4
Mehrez is on a distinguished road
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:
Originally Posted by Jiricbeng View Post
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.
Mehrez is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
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
3D Rotating Channel Curtin University FLUENT 3 April 10, 2012 04:24
rotating channel emilio CD-adapco 1 April 6, 2012 01:21
Setting up for Rotating 2D straight channel pkgupta FLUENT 0 April 2, 2010 06:56


All times are GMT -4. The time now is 01:33.