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

simulate cross section of 3D fully Devolped poiseuille rectangular closed channel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 2, 2012, 17:57
Red face simulate cross section of 3D fully Devolped poiseuille rectangular closed channel
  #1
Member
 
Amin Shariat KHah
Join Date: Apr 2011
Location: Shiraz
Posts: 86
Rep Power: 6
amin144 is on a distinguished road
Hi Dear FOAMers

I want to simulate one cross section (2D) of a 3D rectangular closed channel in fully developed poiseuille flow.
I know should use cyclic BC but I have no idea to use it.

Anybody can upload me a case for study and learning ?
amin144 is offline   Reply With Quote

Old   July 2, 2012, 18:11
Default
  #2
Member
 
Amin Shariat KHah
Join Date: Apr 2011
Location: Shiraz
Posts: 86
Rep Power: 6
amin144 is on a distinguished road
I don't know if cyclic BC is useable for poiseuille flow?
How should I implement pressure gradient?
Pressure is not cyclic in this flow!
amin144 is offline   Reply With Quote

Old   July 3, 2012, 03:28
Default
  #3
Member
 
Amin Shariat KHah
Join Date: Apr 2011
Location: Shiraz
Posts: 86
Rep Power: 6
amin144 is on a distinguished road
NoBody has done some case like this?
amin144 is offline   Reply With Quote

Old   July 3, 2012, 04:20
Default
  #4
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 72
Rep Power: 7
Ohbuchi is on a distinguished road
Hi, Amin.

Please check "channelFoam" source code and "channel395" tutorial.
The channelFoam can handle pressure gradient and cyclic boundary in streamwise direction.
If you want to simulate duct flow, the boundary condition for span direction should be changed to non-slip wall from cyclic boundary.
Ohbuchi is offline   Reply With Quote

Old   July 3, 2012, 07:16
Default
  #5
Member
 
Amin Shariat KHah
Join Date: Apr 2011
Location: Shiraz
Posts: 86
Rep Power: 6
amin144 is on a distinguished road
Hello Dear Masashi

Tanx for your quick answer
What should I do with pressure BC?

Thank you very much...

Last edited by amin144; July 3, 2012 at 09:20.
amin144 is offline   Reply With Quote

Old   July 3, 2012, 20:38
Default
  #6
Member
 
Masashi Ohbuchi
Join Date: Oct 2009
Posts: 72
Rep Power: 7
Ohbuchi is on a distinguished road
Hi, Amin,

Pressure gradient is coded in solver source code as follows.

channelFoam.C
00063 fvVectorMatrix UEqn
00064 (
00065 fvm::ddt(U)
00066 + fvm::div(phi, U)
00067 + sgsModel->divDevBeff(U)
00068 ==
00069 flowDirection*gradP <-- here
00070 );
00071
00072 if (momentumPredictor)
00073 {
00074 solve(UEqn == -fvc::grad(p));
00075 }

So the pressure boundary condition at inlet-outlet should be cyclic.
Ohbuchi is offline   Reply With Quote

Old   July 4, 2012, 03:54
Default
  #7
Member
 
Amin Shariat KHah
Join Date: Apr 2011
Location: Shiraz
Posts: 86
Rep Power: 6
amin144 is on a distinguished road
Hi Dear Masashi

Tanx for your answer
It seems this way doesn't work for another solver? am I right?

So what should I do if I want simulate one cross secttion normal to basis axis of a fully developed rectangular duct with 2d mesh? (my flow is pressure driven)

Best wishes for you
amin144 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
Neumann pressure BC and velocity field Antech Main CFD Forum 0 April 25, 2006 02:15
pressure driven flow by pressure correction method justentered Main CFD Forum 0 December 30, 2003 00:52
Gas pressure question Dan Moskal Main CFD Forum 0 October 24, 2002 22:02
pressure gradient term in low speed flow Atit Koonsrisuk Main CFD Forum 2 January 10, 2002 11:52
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 03:25.