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

flowRateInletVelocity - unexpected results

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 8, 2013, 14:47
Default flowRateInletVelocity - unexpected results
  #1
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 73
Rep Power: 4
ChrisA is on a distinguished road
I'm trying to use flowRateInletVelocity to model a jet-in cross flow using a modified rhoCentralFoam solver (to support multiple gas species). I have a 2d model with a jet width of 0.81mm. I'm trying to specify a mass flow rate of 0.0337kg/(m^2*s), my 2d model has an "empty" depth of 0.01m, so I calculate out an equivalent kg/s (since the units for flowRateInletVelocity are kg/s) to be 2.7317e-7 kg/s and put that into my boundary condition.

I also specify the temperature and total pressure at the jet outlet.

My jet velocity BC looks like:

jetIn
{
type flowRateInletVelocity;
flowRate 2.7317e-7; // Volumetric/mass flow rate [m3/s or kg/s]
value uniform (0 0 0); // placeholder
}

The issue I'm having is that when I look at my results in paraview the mass flow rate is much lower than expected. If I take my species mass fraction * rho * Uy at the boundary over the jet I get a peak value of about 0.0115 kg/(m^2*s) and an average value of about 0.009. This is nowhere near the 0.0337 that I calculated the boundary conditions from. Does anyone have any experience with this BC and able to offer some input/advice? I've attached a plot of my mass flow rate in kg/(m^2*s) in paraView across the jet boundary.

mdot.jpg
ChrisA is offline   Reply With Quote

Old   February 10, 2013, 22:29
Default
  #2
Senior Member
 
Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 7
doubtsincfd is on a distinguished road
Hi Chris,

Quote:
If I take my species mass fraction * rho * Uy at the boundary over the jet
Can you please explain the "boundary over the jet"? Is it the inlet face of the jet?

If possible, can you also post a sketch of your domain, and where you are calculating the values?
doubtsincfd is offline   Reply With Quote

Old   February 11, 2013, 05:04
Default fixedFluxPressure + flowRateInletVelocity to fix mass flow?
  #3
New Member
 
Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 4
tatu is on a distinguished road
Hi Chris!

I am in a very similar situation: I aswell have to specify subsonic mass flow at the inlet, and I have modified rhoCentralFoam to support multiple species. However, at the moment I don't have a working KNP "upwinding" predictor for the species transport equation.

Quote:
I also specify the temperature and total pressure at the jet outlet.
When I used fixedFluxPressure and flowRateInletVelocity for a subsonic inlet (using OF 2.1.1), I got precisely the right mass flow through the inlet. In other words,
Code:
 patchIntegrate phi inletName
outputs the flowRate I specify in flowRateInletVelocity.

However, I get an "unnatural" pressure distribution for the cells near the inlet, so I'm not sure if fixedFluxPressure+flowRateInletVelocity is a valid way to fix inlet BC for subsonic inlet in a compressible flow. I think that you should fix the pressure at the inlet, and let the velocity "float" as a zeroGradient. Can anyone verify this? If I understood correctly, fixedFluxPressure + flowRateInletVelocity fixes a nonuniform pressure gradient, and a nonuniform velocity, at the inlet patch to match the flowRate (phi).

Tatu
tatu is offline   Reply With Quote

Old   February 11, 2013, 06:23
Default
  #4
Senior Member
 
Olivier
Join Date: Jun 2009
Location: France, grenoble
Posts: 235
Rep Power: 9
olivierG is on a distinguished road
hello chris,

You are setting a total mass flow rate at inlet, and you evaluate the mass flow of a specie in paraview ... try rho*Uy, not mass_fraction*rho*Uy.

regards,
olivier
olivierG is offline   Reply With Quote

Old   February 11, 2013, 13:05
Default
  #5
Member
 
Chris
Join Date: Aug 2012
Location: Calgary, Alberta, Canada
Posts: 73
Rep Power: 4
ChrisA is on a distinguished road
Oliver,

This is true, but the jet is blowing a pure gas so along the boundary of interest the mass fraction of interest is equal to 1. I'm interested in controlling the mass flow of said pure gas.

Omkar,

That is correct, I was looking at the inlet face of the jet. My domain is essentially a flat plate (0.2m long) with a small jet blowing perpendicular to the plate at x=0.029m. I'm calculating the values at y=0 in the figure shown (surface of the plate).


Tatu,

Thanks for the input, I'll check out what you're suggesting, although it seems like fixing the pressure with a "floating" velocity would make it difficult to control the mass flow rate.
ChrisA 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
OpenFOAM - Validation of Results Ahmed OpenFOAM Running, Solving & CFD 9 June 22, 2011 18:59
wind tunnel results vs fluent pixie Main CFD Forum 1 August 20, 2009 08:02
Velocity spots in openFoam results Valle OpenFOAM Running, Solving & CFD 4 August 19, 2009 05:53
How to plot a function over a time period? Cirion0000 CFX 4 July 18, 2009 12:48
how to view the complete results shokry FLOW-3D 5 February 3, 2009 14:56


All times are GMT -4. The time now is 15:19.