# Constant flow rate through a small area inside the fluid domain.

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 September 22, 2010, 20:55 Constant flow rate through a small area inside the fluid domain. #1 Member   Robin Gilbert Join Date: Jan 2010 Posts: 66 Rep Power: 16 Hi, I want to force a fluid at a constant flow rate through a small area inside the fluid domain. I know about the cyclic boundary condition but i do not know how to set it up so that it will give a constant flow rate. I looked into channelFoam and i know how to modify the simpleFoam solver so that it can fix the flowrate. but i just cant figure out how to make it so that the forcing is only through a small area inside the fluid domain. I considered adding a forcing term to the UEqn which is something like: (v/A-U)*alpha where v is volume flow rate A-area of the surface U- current velocity at the surface alpha- a term which has dimensions 1/time which can drive down the solution to required flow rate. is this approach right? Please give me some pointers.

 September 23, 2010, 11:37 #2 Senior Member   Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 21 The only way to do what you want, is literally to have a fixed mass flow outlet adjacent to a fixed flow inlet that maps the velocity from the adjacent outlet (and vice-versa for pressure). You can do this by starting from the low level coupled patches (cyclics), but it isn't trivial unless you know what you are doing and the resulting system probably wont be very stable. Unfortunately I don't think there is an easy solution for this particular aspect of the problem. A different approach is to put in something like an actuator disc from 1.7 wind solver and adjust the driving force such that the flow rate approximates your target value.

 September 23, 2010, 17:45 #3 Member   Robin Gilbert Join Date: Jan 2010 Posts: 66 Rep Power: 16 Thank you Eugene, I will try the actuator disk and see if that works. So according to ur first suggestion, do i need to make a hollow volume without any meshes and has an inlet (into the fluid domain) and an outlet (into the hollow volume) and map the velocity from inlet to the outlet?

September 23, 2010, 18:48
#4
Senior Member

Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 21
Quote:
 Originally Posted by robingilbert Thank you Eugene, I will try the actuator disk and see if that works. So according to ur first suggestion, do i need to make a hollow volume without any meshes and has an inlet (into the fluid domain) and an outlet (into the hollow volume) and map the velocity from inlet to the outlet?
No, no - no holes required. You take your basic cyclic boundary (which can act just like an internal boundary) and start chopping and changing. For U you want one side of the cyclic patch pair to be a scaled "zero gradient" outlet such that the total flux matches your specification. The other side is an inlet that simply gets its value from the outlet side, just like a normal cyclic. For pressure, you want the reverse, the inlet side is zero gradient and this value is mapped to the outlet side and treated as fixed. The switching is done based on flux direction. I did something like this at my previous job and remember that it was a bit tricky. Unfortunately I don't have the code any more, but it is certainly possible.

 September 23, 2010, 21:05 #5 Member   Robin Gilbert Join Date: Jan 2010 Posts: 66 Rep Power: 16 Thanks once again Eugene. one more question. forgive me if its a stupid question, can I add a source term in the momentum eqn so that it will force the fluid at a particular flow rate?

 September 24, 2010, 04:49 #6 Senior Member   Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 21 Of course. I pointed you toward the actuator disk because this is exactly the framework you need to implement the source term you need. An example usage can be found here: OpenFOAM-1.7.x/tutorials/incompressible/simpleWindFoam and the source code for the actuator disk and base class can be found here: OpenFOAM-1.7.x/src/finiteVolume/cfdTools/general/fieldSources/basicSource

 September 24, 2010, 05:35 #7 Member   Robin Gilbert Join Date: Jan 2010 Posts: 66 Rep Power: 16 Thank you soooooooooo much Eugene. I will try that out!! thank you once again.

 October 4, 2010, 16:19 #8 Senior Member     maddalena Join Date: Mar 2009 Posts: 436 Rep Power: 23 Hi Robin, I probably need something similar for this kind of problem: http://www.cfd-online.com/Forums/ope...flow-rate.html. Eugene suggested me to modify the actuator disk model as well. Have you made any progress on the subject? Shall we collaborate? Regards mad

 Tags simplefoam, source, source term, ueqn

 Thread Tools Search this Thread Search this Thread: Advanced Search 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 Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Young CFX 5 October 6, 2008 23:17 jane luo Main CFD Forum 15 April 12, 2004 17:49 Abhi Main CFD Forum 12 July 8, 2002 09:11 Eric Poindexter Main CFD Forum 2 September 22, 2000 09:21 ram Main CFD Forum 5 June 17, 2000 21:31

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

 Contact Us - CFD Online - Privacy Statement - Top