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

cyclic boundary condition doesn't work in a tube ??

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

Like Tree1Likes
  • 1 Post By David_010

Reply
 
LinkBack Thread Tools Display Modes
Old   April 15, 2010, 08:51
Default cyclic boundary condition doesn't work in a tube ??
  #1
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 229
Rep Power: 9
Cyp is on a distinguished road
Hi !

I try to simulate a flow through a tube (2D) with cyclic boundary condition at the inlet and the outlet.

A poiseuille's flow profile (parabolic) is expected along the tube. However, it works only after a distance from the inlet. I belevied that with the cyclic condition, I will get the same profile all along the tube. Why does my simulation is like that ??

I obtain exactly the same profile when I define my inlet and outlet as "patch".

My initial conditions are :

p :
inlet
{
type zeroGradient;
}

outlet
{
type fixedValue;
value uniform 0;


U :
outlet
{
type zeroGradient;
}

inlet
{
type fixedValue;
value uniform (4 0 0);
}


thank you for your help,
Cyp
Cyp is offline   Reply With Quote

Old   April 15, 2010, 08:55
Default
  #2
Member
 
David
Join Date: Dec 2009
Location: Spain
Posts: 58
Rep Power: 7
David_010 is on a distinguished road
Can you put here the cyclic condition you have used?

David
David_010 is offline   Reply With Quote

Old   April 15, 2010, 09:05
Default
  #3
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 229
Rep Power: 9
Cyp is on a distinguished road
In my blockMeshDict I defined :

patches
(
cyclic inlet
(
(0 4 7 3)
)
cyclic outlet
(
(1 5 6 2)
)
wall fixedWalls
(
(3 7 6 2)
(1 5 4 0)
)
empty frontAndBack
(
(0 3 2 1)
(4 5 6 7)
)
);


It is the question you asked me ?
Cyp is offline   Reply With Quote

Old   April 15, 2010, 11:00
Default
  #4
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 229
Rep Power: 9
Cyp is on a distinguished road
Indeed, I haven't well defined my cyclic condition.

It works now when I defined inout instead of inlet and outlet patches.

cyclic inout
(
(0 4 7 3)
(1 5 6 2)
)


For my pressure initial condition, I used

inout
{
type fan;
patchType cyclic;
f List<scalar> 1(-5.00); // p_OF = p_real / rho
value uniform 0;
}

Is it the single manner to do ?
Cyp is offline   Reply With Quote

Old   April 15, 2010, 12:11
Default
  #5
Member
 
David
Join Date: Dec 2009
Location: Spain
Posts: 58
Rep Power: 7
David_010 is on a distinguished road
Yes it is the correct manner. You can also use de b.c directMapped, which takes the outlet velocity and recicles it to the inlet patch.

David
David_010 is offline   Reply With Quote

Old   April 15, 2010, 14:11
Default
  #6
Cyp
Senior Member
 
Cyprien
Join Date: Feb 2010
Location: Stanford University
Posts: 229
Rep Power: 9
Cyp is on a distinguished road
I improved my case considering two-phase flow. I use the interFoam solver.

I patch well my region.

When I launch the interFoam solver, I get this error message :

Code:
Unable to set reference cell for field p
    Please supply either pRefCell or pRefPoint


file: /home/cyp/OpenFOAM/cyp-1.6/run/test07/system/fvSolution::PISO from line 55 to line 60.

    From function void Foam::setRefCell
(
    const volScalarField&,
    const dictionary&,
    label& scalar&,
    bool
)
    in file cfdTools/general/findRefCell/findRefCell.C at line 112.

FOAM exiting
Do you know what does it mean and what's wrong in my case ??
Cyp is offline   Reply With Quote

Old   April 16, 2010, 15:21
Default
  #7
Member
 
David
Join Date: Dec 2009
Location: Spain
Posts: 58
Rep Power: 7
David_010 is on a distinguished road
I don't know what is the problem. When you tried the cyclic condition and it worked, was whit this solver? If not, perhaps this solver doesn't work with the cyclic condition and this b.c for pressure. But I'm just guessing. I suppose you want to fix the mass flux? if not, you can set only:

inout
{
type cyclic;
value uniform 0;
}

for pressure. And fix the velocity in the "boundary field" (0/U). If it still doesn't work...I don't know. But may be the directMapped condition is a good alternative.

Sorry for my poor help

David
mohammad81 likes this.
David_010 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
Boundary Conditions Thomas P. Abraham Main CFD Forum 20 July 7, 2013 05:05
Cyclic Boundary Condition SG CD-adapco 0 June 1, 2008 14:56
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15
About Partial Cyclic Boundary Condition Jiaying Xu CD-adapco 2 October 31, 2002 21:12
The Boundary Condition about the Flat Plate boing Main CFD Forum 1 January 6, 2002 17:53


All times are GMT -4. The time now is 22:07.