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

ChannelOodles Correcting driving force for a constant mass flow rate

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 20, 2005, 10:39
Default channelOodles: I'm trying t
  #1
Senior Member
 
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 18
maka is on a distinguished road
channelOodles:

I'm trying to edit the channelOodles to use it for channelDNS. I found some part of the code that is about correcting driving force, this part does not exist in icoFoam, which means: it is either included because of the channel or because of LES? I'm a ware that the driving force in the channel enters in the non-periodic part of pressure gradient, but I do not know about such correction. Can any one explain why we need this part in the channel code or give any reference? Many thanks!

---------------------------------------------
// Correct driving force for a constant mass flow rate

// Extract the velocity in the flow direction
dimensionedScalar magUbarStar =
(flowDirection & U)().weightedAverage(mesh.V());

// Calculate the pressure gradient increment needed to
// adjust the average flow-rate to the correct value
dimensionedScalar gragPplus =
(magUbar - magUbarStar)/rUA.weightedAverage(mesh.V());

U += flowDirection*rUA*gragPplus;

gradP += gragPplus;

Info<< "Uncorrected Ubar = " << magUbarStar.value() << tab
<< "pressure gradient = " << gradP.value() << endl;


/Maka
maka is offline   Reply With Quote

Old   September 20, 2005, 12:58
Default All to do with the 'channel' b
  #2
Senior Member
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,419
Rep Power: 26
mattijs is on a distinguished road
All to do with the 'channel' bit, not with LES.

One specifies an average flow-rate (Ubar in constant/transportProperties) and the pressure gradient is (hopefully slightly) adjusted to keep the wanted flow-rate.
mattijs is offline   Reply With Quote

Old   September 21, 2005, 04:56
Default In my case, I'm interested
  #3
Senior Member
 
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 18
maka is on a distinguished road
In my case,

I'm interested to keep the dimensionless driving pressure gradient (non-periodic part) (scaled by u_tau and half channel width) to be -1 while specifying Re_tau for the simulation. The Ubar can come as a result of the simulation. I think in this case, I do not need the correction step, am I right? On the other hand I was thinking to add this deriving pressure gradient as a source term in Su, what do you think? Many thanks for your help.

Best regards,
Maka.
maka is offline   Reply With Quote

Old   September 21, 2005, 05:18
Default Yes, you are correct, I have d
  #4
Senior Member
 
Eugene de Villiers
Join Date: Mar 2009
Posts: 725
Rep Power: 21
eugene is on a distinguished road
Yes, you are correct, I have done this myself.

The pressure gradient was excluded from the momentum matrix for a reason, I just cant remember why. It should be simple to try both options and see if there is a difference.
eugene 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
periodic, unsteady, non-constant mass flow rate Lee FLUENT 42 October 5, 2011 15:20
maintain a constant flow rate in channel flow ? Lewis Main CFD Forum 2 September 28, 2010 11:35
Constant Volumetric Flow Rate Boundary Condition Hadley Main CFD Forum 0 October 20, 2005 23:13
Possible?: Periodic conditions with non-constant mass flow rate. Ray FLUENT 0 April 10, 2000 06:10
Constant mass flow rate / choking outlet bounardy condition Min Zhu Main CFD Forum 1 September 29, 1998 15:33


All times are GMT -4. The time now is 14:09.