CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

InterFoam and Channel Flow

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 8, 2008, 19:17
Default I'm attempting to adapt the in
  #1
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 525
Rep Power: 17
chegdan will become famous soon enough
I'm attempting to adapt the interFoam solver to incorporate cyclic boundary conditions. I have had a look at the ChannelOodles solver and added several files (mainly the createGrad.H and writeGradP.H). The solver was able to compile and I have run the damBreak case with no problems, receiving similar result to the interFoam solver. However, when I apply cyclic BC's to another case, then there are problems. I must choose a ridiculously small courant number (1e-5 or less) to even have continuity converge and even with that the pressure gradient still diverges. I was unable to attach such large files (a mere 250 kb) so here is a link where they can be downloaded. Please take a look at it, see if there are any issues and maybe make some comments.

Thanks for any help that can be given. Also, I will be sure to share any knowledge gained with the community.

Dan
__________________
Dan

Find me on twitter @dancombest
chegdan is offline   Reply With Quote

Old   September 11, 2008, 21:17
Default Hello again, I seem to be
  #2
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 525
Rep Power: 17
chegdan will become famous soon enough
Hello again,

I seem to be posting a lot lately, so I'll just continue to do so until I figure out my problem. Let me get to it. I'm trying to make a multiphase channel flow solver called rasInterChannelFoam (it borrows from rasInterFoam and channelOodles). So far I have been able to run cases where all of the fields agree with the cyclic boundary called "inOut" except for the gamma field. I have a simple case that as a block of water approaches the lower boundary, it is supposed to past through and fall from the top boundary. However, when the brick contacts the lower boundary, it compresses and then only passes a small fraction through the interface to the top boundary. Here are two pictures taken during the hit and just after the contact with the boundary.





I have been able to simulate single phase flow (with and without turbulence) that shows vortex shedding and field gradients (U, k, epsilon, etc.) pass through to through the cyclic boundary. However, as you can see there is no gamma passed through. Does anyone have a suggestion as to how I can correct this? Do I need to change the gamma equation to no include the compression term (just a thought although I'm not really sure how that will affect anything) or do I need to define another boundary type/condition to pass gamma? Lastly, is this a bug? I am using OF-1.5.


Procedure for generating mesh in Gambit and importing into OF:

1.Draw geometry (a box)
2.link top and bottom edges
3.mesh the object
4.define the top and bottom edge as one boundary called "inOut" of type "wall" and the sides as one boundary called "wall of type "wall".
5.Export the 2D mesh as .msh format.
6.Use fluentMeshToFoam utility to import the mesh
7.go into the polymesh folder and change the "inOut" to base type "cyclic".
8.After that I make sure that Ubar in the transport properties file has a low value that matches the initial uniform field of U.
9.Run my case.
10.Watch the crazy brick of liquid bounce back.

Here is a link to the video on youtube so you can see for yourself. Any Help is much appreciated. Here is a link to my zipped solver and test case.

Dan
__________________
Dan

Find me on twitter @dancombest
chegdan is offline   Reply With Quote

Old   December 14, 2011, 11:26
Default
  #3
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 525
Rep Power: 17
chegdan will become famous soon enough
I finally found the solution to this one, but never replied again to my own thread.

interFoam and cyclic boundary condition issue

and

Cyclic boundaries importing Gambit msh file and interFoam

better late than never i guess. I know there are other answers out there, but this is just closing this thread.
__________________
Dan

Find me on twitter @dancombest
chegdan 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
maintain a constant flow rate in channel flow ? Lewis Main CFD Forum 2 September 28, 2010 12:35
InterFoam and Channel Flow chegdan OpenFOAM Running, Solving & CFD 1 September 9, 2008 02:13
InterFoam with vacuum in the flow domain possible hsieh OpenFOAM Running, Solving & CFD 0 May 1, 2008 11:02
Outlet boundary condition in channel flow solving with interFoam liuzhw OpenFOAM Running, Solving & CFD 0 November 29, 2005 22:07


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