CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   ANSYS Meshing & Geometry (https://www.cfd-online.com/Forums/ansys-meshing/)
-   -   [GAMBIT] periodic boundary condition for a rectangular duct domain (https://www.cfd-online.com/Forums/ansys-meshing/115703-periodic-boundary-condition-rectangular-duct-domain.html)

Das April 5, 2013 01:01

periodic boundary condition for a rectangular duct domain
 
Hi all,

I am new to CFD and i m doing a work on the film cooling on a flat plate, and i used gambit for modelling a 3D rectangular geometry as a domain. The bottom face of the rectangular duct is the taken as the test plate.

I gave the following boundary conditions:-

1. duct inlet > Velocity inlet
2. duct outlet> pressure outlet
3. test plate > wall
4. Top side> wall( opp face of test plate)
5. right side face> symmetry
6. left side face> Periodic

After that i did export > mesh, but the file is not exporting successfully.:confused:

It is showing some errors with the periodic boundary condition,
warn: boundary entity periodic side shadow does not contain any valid entity and is not written.
error: periodic boundary entities are not specified on the same boundary entity.

please help... what is the error and how can I resolve it...

thanks in advance

-mAx- April 5, 2013 08:44

can you attach a sketch of your geometry with your BC.

Das April 5, 2013 12:47

2 Attachment(s)
thanks for showing interest to help me, sir
here i m attaching its sketch...

i did this in gambit...
for the left face of the duct i want to giv periodic boundary condition.:confused:

how can I resolve it...sir

-mAx- April 5, 2013 13:10

can show where is the periodic BC on the picture?
Periodicity has some rules to be well defined
Anyway your splits around your cylinder are wrong, since it will generate skewed cells.

Das April 7, 2013 05:14

1 Attachment(s)
i gave periodic boundary condition for the flat side(opposite to the face having semicircle groove) of the rectangular domain.

sir, here i m supposed to model a film cooling over a falt plate.

in actual experimental setup consist of a rectangular duct with a test plate at the base portion of the rectangle, where the test plate contains certain number of holes inclined at 35 deg to the streamwise direction. hot stream of air is fed through the inlet pace of the rectangle while coolant air is blown through inclined holes.

for reducing the iteration, i took one half of the coolant hole for modelling...

here i m attaching the top view of the test plate of actual experimental setup.

that section shown by dotted line is taken in to consideration for the simulation.

-mAx- April 8, 2013 01:48

ok
Then you shouldn't use the symetry, and you have to model the reflected zone.
Thus, as I said, periodic BC has some rules: mesh from both periodic surfaces has to be linked, and you should give at least 2 surfaces (not coplanar) in the BC set
check this tutorial: http://aerojet.engr.ucdavis.edu/flue.../tg/node26.htm

Das April 9, 2013 05:08

thank you...sir:)

niloogh May 23, 2013 07:21

1 Attachment(s)
Quote:

Originally Posted by -mAx- (Post 418972)
ok
Then you shouldn't use the symetry, and you have to model the reflected zone.
Thus, as I said, periodic BC has some rules: mesh from both periodic surfaces has to be linked, and you should give at least 2 surfaces (not coplanar) in the BC set
check this tutorial: http://aerojet.engr.ucdavis.edu/flue.../tg/node26.htm

hi max
would you please take a look at my geometry.
i don't know to use symmetry BC or periodic for the surfaces which i had labeled in the photo:confused:

it is a 120 degree model of a wind turbine.

-mAx- May 23, 2013 08:01

periodic: all surfaces which belongs to plane (yz) are supposed to be mesh-linked with the surfaces which belongs to 120°-rotated (yz)
That means that both suface-groups have exactly same mesh, which allows you to use periodicity.
Now I don't remember how to define the periodic BC: all surfaces as one BC, or each plane as a BC, etc...
But I think, you only need one BC

niloogh May 23, 2013 08:30

Quote:

Originally Posted by -mAx- (Post 429496)
periodic: all surfaces which belongs to plane (yz) are supposed to be mesh-linked with the surfaces which belongs to 120°-rotated (yz)
That means that both suface-groups have exactly same mesh, which allows you to use periodicity.
Now I don't remember how to define the periodic BC: all surfaces as one BC, or each plane as a BC, etc...
But I think, you only need one BC



you mean,after meshing all the volume,i must go to "link face meshes" panel and choose the faces which are going to be linked??

and i don't know to choose what for vertices which are needed in this panel???:confused:

-mAx- May 23, 2013 09:14

1 Attachment(s)
no before meshing the volume. Thus you will have to modify the surface mesh.
I confirm you need to set all your periodic surfaces as one periodic BC (in my case the 4 color surfaces)
Attachment 22045

niloogh May 23, 2013 16:23

Quote:

Originally Posted by -mAx- (Post 429513)
no before meshing the volume. Thus you will have to modify the surface mesh.
I confirm you need to set all your periodic surfaces as one periodic BC (in my case the 4 color surfaces)
Attachment 22045

tnx for your help.
i had done it successfully but

after using this mesh in fluent and starting itarating i recived this message:
Error: divergence detected in AMG solver: pressure correction
Error: divergence detected in AMG solver: x-momentum
Error: divergence detected in AMG solver

why is that?
is that related to my mesh???

-mAx- May 24, 2013 01:13

could be many things
Check all your BC (set up).
If you don't see anything, then you may try to set your periodic BC as wall, and iterate (just to see if it comes from periodicity).
And off course a check mesh

niloogh May 24, 2013 04:34

Quote:

Originally Posted by -mAx- (Post 429677)
could be many things
Check all your BC (set up).
If you don't see anything, then you may try to set your periodic BC as wall, and iterate (just to see if it comes from periodicity).
And off course a check mesh

i,m sure that it comes from my periodic boundary.
because when i used symmetry instead of periodic.i did,t have this problem.
but i don't know what i must change about this periodic bcs:(

-mAx- May 24, 2013 04:52

how did you set it?

niloogh May 24, 2013 05:24

Quote:

Originally Posted by -mAx- (Post 429716)
how did you set it?

i linked the faces as you said just without ticking the periodic in "link face meshes" panel.
after that i meshed all of my volume with tetra.and set the periodic for every both linked faces and saved it as mesh and sent it to fluent.

-mAx- May 24, 2013 05:31

I meant how did you set it in fluent?
http://aerojet.engr.ucdavis.edu/flue...4.htm#Periodic

niloogh May 24, 2013 05:43

Quote:

Originally Posted by -mAx- (Post 429725)
I meant how did you set it in fluent?
http://aerojet.engr.ucdavis.edu/flue...4.htm#Periodic

I set it as rotational.and i didn't do any thing else in fluent for this bc.

-mAx- May 24, 2013 06:05

what gives your a check mesh?

niloogh May 24, 2013 06:18

Quote:

Originally Posted by -mAx- (Post 429732)
what gives your a check mesh?

Grid Check

Domain Extents:
x-coordinate: min (m) = -1.194447e+005, max (m) = -1.000000e+002
y-coordinate: min (m) = -1.680000e+005, max (m) = 6.120000e+005
z-coordinate: min (m) = -1.200000e+005, max (m) = 6.000000e+004
Volume statistics:
minimum volume (m3): 1.474610e+002
maximum volume (m3): 2.305678e+012
total volume (m3): 1.166610e+016
Face area statistics:
minimum face area (m2): 5.263843e+000
maximum face area (m2): 3.819772e+008
Checking number of nodes per cell.
Checking number of faces per cell.
Checking thread pointers.
Checking number of cells per face.
Checking face cells.
Checking bridge faces.
Checking right-handed cells.
Checking face handedness.
Checking face node order.
Checking element type consistency.
Checking boundary types:
Checking face pairs.
Checking periodic boundaries.
Zone 6: x-translation: min (m) = 0.000000e+000, max (m) = 1.029000e+005
y-translation: min (m) = 0.000000e+000, max (m) = 0.000000e+000
z-translation: min (m) = -1.800000e+005, max (m) = 0.000000e+000
Zone 8: x-translation: min (m) = -1.029000e+005, max (m) = 0.000000e+000
y-translation: min (m) = -3.955078e-002, max (m) = 3.878906e+000
z-translation: min (m) = 0.000000e+000, max (m) = 1.800000e+005
Zone 10: x-translation: min (m) = -1.029000e+005, max (m) = 0.000000e+000
y-translation: min (m) = 0.000000e+000, max (m) = 0.000000e+000
z-translation: min (m) = 0.000000e+000, max (m) = 1.800000e+005
Checking node count.
Checking nosolve cell count.
Checking nosolve face count.
Checking face children.
Checking cell children.
Checking storage.
Done.

-mAx- May 24, 2013 06:29

you can provide your mesh file, I can check but only on Monday.

Are you sure you don't have this problem with symetry instead of periodic?
What about wall instead of periodic?

Give all the setup, density, viscosity etc...

You are handling a domain with huge dimensions (order of 100 km), is it right?
Seems to be quite big for a impeller if I refered to your picture, isn't it?

niloogh May 24, 2013 07:29

Quote:

Originally Posted by -mAx- (Post 429747)
you can provide your mesh file, I can check but only on Monday.

Are you sure you don't have this problem with symetry instead of periodic?
What about wall instead of periodic?

Give all the setup, density, viscosity etc...

You are handling a domain with huge dimensions (order of 100 km), is it right?
Seems to be quite big for a impeller if I refered to your picture, isn't it?

yes max
i,m sure that with symmetry and with wall in stead of periodic there is no problem.
i,m simulating a wind turbine which the radios of it,s rotor is 63 meters
the radius of outer volume is 180 meters and the height of that is about 900 meters.
the other bcs and situation are as:
velocity inlet 11 m/s
pressure outlet 0 gauge
tow faces as interior
wind turbine wall
i,m starting with simple situation so i set my turbine as stationary wall and laminar situation.

-mAx- May 24, 2013 07:43

2 points:

*
Domain Extents:
x-coordinate: min (m) = -1.194447e+005, max (m) = -1.000000e+002

--> DX=119347m = 119 km
y-coordinate: min (m) = -1.680000e+005, max (m) = 6.120000e+005

--> DY=780000m=780 km
z-coordinate: min (m) = -1.200000e+005, max (m) = 6.000000e+004

--> DZ=180000m=180 km

In Gambit how big is your outer volume? (Gambit has no unit, if you only import the mesh, then it is automatically scaled in m)

*Did you specify your 2 faces as interior? (In Gambit, interior faces are automatically fused in interior domain). Do they appear as interior or wall in fluent?

niloogh May 24, 2013 07:55

Quote:

Originally Posted by -mAx- (Post 429767)
2 points:

*
Domain Extents:
x-coordinate: min (m) = -1.194447e+005, max (m) = -1.000000e+002

--> DX=119347m = 119 km
y-coordinate: min (m) = -1.680000e+005, max (m) = 6.120000e+005

--> DY=780000m=780 km
z-coordinate: min (m) = -1.200000e+005, max (m) = 6.000000e+004

--> DZ=180000m=180 km

In Gambit how big is your outer volume? (Gambit has no unit, if you only import the mesh, then it is automatically scaled in m)

*Did you specify your 2 faces as interior? (In Gambit, interior faces are automatically fused in interior domain). Do they appear as interior or wall in fluent?

i bring my geometry from catia into gambit,in catia they were in mm
so i interred my amounts in gambit with mm

i set that tow faces as interior in fluent too.

-mAx- May 24, 2013 08:05

*then you have to scale your model
*you cannot set your faces as interior in fluent, since they are not supposed to exist, except if you set them in Gambit. Did you set them in Gambit as interior?
If not, then fluent shouldn't create them, and then something went wrong.

niloogh May 24, 2013 18:09

Quote:

Originally Posted by -mAx- (Post 429772)
*then you have to scale your model
*you cannot set your faces as interior in fluent, since they are not supposed to exist, except if you set them in Gambit. Did you set them in Gambit as interior?
If not, then fluent shouldn't create them, and then something went wrong.

dear max tnx really for all of your helps
i set the correct amounts for one of my walls and now it is working:D

Avistar June 5, 2013 07:48

Select boundaries as " wall "type only then give the names to the boundaries such that you can easily identify the which boundary is periodic and his shadow.
then In Fluent use command
mesh<enter>modify zones<enter>make periodic
give number of periodic and shadow boundaries <enter>
use yes for rotational no for translation


All times are GMT -4. The time now is 00:08.