CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Meshing & Mesh Conversion

Cyclic mesh for fully-developed profile

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

Like Tree1Likes
  • 1 Post By stevenvanharen

Reply
 
LinkBack Thread Tools Display Modes
Old   February 13, 2012, 17:01
Default Cyclic mesh for fully-developed profile
  #1
Member
 
Join Date: Nov 2011
Posts: 44
Rep Power: 5
fferroni is on a distinguished road
Hello

I am trying to find the fully developed velocity profile in a straight duct.

Rather than compute a 3D duct, and look at its profile towards the end, I would like to apply periodic boundary conditions to make it more efficient. I have never done this on openFOAM and would like some help.

I have a inlet and an outlet boundary. So:


boundary
(
inlet
{
type cyclic;
neighbourPatch outlet;
faces
(
(0 4 7 3)
);
}

outlet
{
type cyclic;
neighbourPatch inlet;
faces
(
(2 6 5 1)
);


However, now I am wondering what kind of boundary field definitions I need.
Previously I had zeroGradient or fixedValue fields (depends on variable) for both the inlet and outlet. What do I need to set for the periodic case?

Thank you very much.

Kindest Regards,

F
fferroni is offline   Reply With Quote

Old   February 14, 2012, 07:17
Default
  #2
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
Take a look at the channelFoam tutorial. This will help you I think.

Please keep in mind that the pressure field is never periodic. Therefore you will have to split the pressure field in a periodic part and a pressure drop. This can be seen in the channelFoam solver.
Chandsome likes this.
stevenvanharen is offline   Reply With Quote

Old   February 14, 2012, 08:27
Default
  #3
Member
 
Join Date: Nov 2011
Posts: 44
Rep Power: 5
fferroni is on a distinguished road
I see..

So, it means that I need to modify my solver?

Looking at channelFoam.c, I would need to slightly change the first part before the PISO loop and add the section correcting the driving force for a constant mass flow rate. I would not need to change any field definitions, apart from creating createGradP.H, is that correct?

In the first section, what is "sgsModel->correct(); " ?

Kindest regards,

Francesco
fferroni is offline   Reply With Quote

Old   February 14, 2012, 08:40
Default
  #4
Senior Member
 
Steven van Haren
Join Date: Aug 2010
Location: The Netherlands
Posts: 149
Rep Power: 6
stevenvanharen is on a distinguished road
That you can do indeed, or you can use a jumpcyclic. This has different names for different OF versions I think. Browse around the OF source to see the difference.

src » finiteVolume » fields » fvPatchFields


The same holds for sgsModel->correct(); Browse around the turbulence models to see what it does.

src » turbulenceModels
stevenvanharen is offline   Reply With Quote

Old   February 15, 2012, 03:58
Default
  #5
Member
 
Join Date: Nov 2011
Posts: 44
Rep Power: 5
fferroni is on a distinguished road
Hmm, I don't quite understand how jumpCyclic works (especially trying to read the source code).. Is there a tutorial on how to use it?

There is a fan tutorial in OF that apparently uses it but I'm not sure how that is relevant to what I'm trying to do...

Kind Regards,

F
fferroni is offline   Reply With Quote

Old   July 11, 2012, 10:49
Default
  #6
Senior Member
 
Jian Zhong
Join Date: Feb 2012
Location: Birmingham
Posts: 108
Rep Power: 5
zxj160 is on a distinguished road
Quote:
Originally Posted by fferroni View Post
Hello

I am trying to find the fully developed velocity profile in a straight duct.

Rather than compute a 3D duct, and look at its profile towards the end, I would like to apply periodic boundary conditions to make it more efficient. I have never done this on openFOAM and would like some help.

I have a inlet and an outlet boundary. So:


boundary
(
inlet
{
type cyclic;
neighbourPatch outlet;
faces
(
(0 4 7 3)
);
}

outlet
{
type cyclic;
neighbourPatch inlet;
faces
(
(2 6 5 1)
);


However, now I am wondering what kind of boundary field definitions I need.
Previously I had zeroGradient or fixedValue fields (depends on variable) for both the inlet and outlet. What do I need to set for the periodic case?

Thank you very much.

Kindest Regards,

F
Dear F,

I also met the same problem that I want to set zeroGradient or fixedValue fields (depends on variable) for both the cyclic inlet and outlet. Did you find a way to achieve it?

My best regards,
Jian
zxj160 is offline   Reply With Quote

Reply

Tags
cyclic boundaries, openfoam, periodic boundary

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
Mesh motion with Translation & Rotation Doginal CFX 2 January 12, 2014 07:21
No layers in a small gap bobburnquist OpenFOAM Native Meshers: snappyHexMesh and Others 2 November 25, 2012 09:54
Interpolation Error on FAM Mesh with Cyclic BCs ngj OpenFOAM Bugs 1 August 9, 2011 06:12
problem when I import mesh with cyclic graduated BC Cyp OpenFOAM 0 March 3, 2011 11:38
inlet velocity profile for fully developed turbule pankaj saha Main CFD Forum 4 August 24, 2007 16:53


All times are GMT -4. The time now is 13:55.