# Open Channel Flow using InterFoam type solver

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

 April 6, 2009, 05:25 Open Channel Flow using InterFoam type solver #1 New Member   Su Xiaohui Join Date: Mar 2009 Location: Singapore Posts: 29 Rep Power: 9 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

 April 21, 2009, 02:47 #2 New Member   Su Xiaohui Join Date: Mar 2009 Location: Singapore Posts: 29 Rep Power: 9 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

 May 5, 2009, 10:28 #3 New Member   HuangWei Join Date: Mar 2009 Posts: 8 Rep Power: 9 hello, Su Xiaohui I am also studying in the interFoam and thanks for your share.

 May 5, 2009, 21:58 #4 Member   Jitao Liu Join Date: Mar 2009 Location: Jinan , China Posts: 64 Rep Power: 9 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post chegdan OpenFOAM Running, Solving & CFD 1 September 9, 2008 01:13 Mehdi Moayeri FLUENT 0 August 26, 2008 10:55 Willem Brantegem FLUENT 2 April 4, 2007 02:40 Willem Brantegem Main CFD Forum 0 April 3, 2007 09:39 yan FLUENT 0 July 4, 2005 23:36

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