CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Open Channel Flow using InterFoam type solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 6, 2009, 05:25
Default Open Channel Flow using InterFoam type solver
  #1
Member
 
Su Xiaohui
Join Date: Mar 2009
Location: Singapore
Posts: 30
Rep Power: 17
sxhdhi is on a distinguished road
Send a message via Skype™ to sxhdhi
Dear Foamers,

Please allow me to create this thread for those of who are interesting in Open Channel Flow simulation by using InterFoam type solvers such as InterFoam, rasInterFOam, and so on.

I have spend over one month on trying to run open channel flow case using InterFOAM, but result is very disappointed. After that, I gave up for testing various BCs and instead it I start to look into the source code of interFoam.

From now on, I will write down and share my doubt, ideas and I hope Foamers who are interesting in this part can join me to fight the problems together in order to solve it.

today, I would like to write down the Governing equations used in InterFoam

(1) continuity Eq
div(U)=0
(2) momentum Eq
ddt(rho,U)+div(rho*U*U)-rho*g=laplacian(muf,U)-grad(p)-Fs
(3) free surface capturing Eq
ddt(gamma)+div(U,gamma)+div((1-gamma)*Ur,gamma)=0

where p=pd+rho*gh (4)

laplacian Eq for pd
(5)
laplacian(rUa, pd)=1/ap*[H(U)-div(Fs)-g*div(rho)-gh*div(rho)]

now I have two questions during looking into the source code:
(1) in UEqu.C where is the source term rho*g, this term should be in R.H.S.
(2) In PEqu.C still miss source term rho*g and I cannot understand the laplacian equation for pd
laplacian(rUa, pd)=div(phi)
in my view, phi is already the term after intergretion. why still have to do divergent process.

find problem---> solve it --> new problem -->solve again!!!!!

Regards

Su Xiaohui
sxhdhi is offline   Reply With Quote

Old   April 21, 2009, 02:47
Default
  #2
Member
 
Su Xiaohui
Join Date: Mar 2009
Location: Singapore
Posts: 30
Rep Power: 17
sxhdhi is on a distinguished road
Send a message via Skype™ to sxhdhi
It is so regret that no one is interested in this thread.

I am still focusing on this problem.

What I found now is p=pd+rho*gh, where rho*gh is not particularly pointed to hydrostatic pressure because in code, gh=g& mesh.C(), the real hydrostatic presure shoul be the expression rho*g*(free surface -mesh.C()).

I also found out that in Ueqn and Peqn, source term (gravity term) may not be used.

please comments if anyone is interested in this.

Best regards,

Su Xiaohui
sxhdhi is offline   Reply With Quote

Old   May 5, 2009, 10:28
Default
  #3
New Member
 
HuangWei
Join Date: Mar 2009
Posts: 8
Rep Power: 17
huoyinghw is on a distinguished road
hello, Su Xiaohui
I am also studying in the interFoam and thanks for your share.
huoyinghw is offline   Reply With Quote

Old   May 5, 2009, 21:58
Default
  #4
Member
 
Jitao Liu
Join Date: Mar 2009
Location: Jinan , China
Posts: 64
Rep Power: 17
awacs is on a distinguished road
Hi,Su Xiaohui

I am interested in multiphase flow too.I have modified the interFoam solver to simulate three dimensional,non-newtonian,nonisothermal,incompressible,laminar flow.The result seems good.
My problem is that when I use the added Cross-WLF viscosity model,the Courant Number increase enormously after a few steps.Then the calculation stopped.
I'll pay more attention to it.

Thanks for your thread.

Yours,
Jitao Liu
awacs is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
InterFoam and Channel Flow chegdan OpenFOAM Running, Solving & CFD 1 September 9, 2008 01:13
Open channel flow Mehdi Moayeri FLUENT 0 August 26, 2008 10:55
pressure outlet (open channel flow) Willem Brantegem FLUENT 2 April 4, 2007 02:40
pressure outlet (open channel flow) Willem Brantegem Main CFD Forum 0 April 3, 2007 09:39
open channel flow boundary conditions yan FLUENT 0 July 4, 2005 23:36


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