# flowRateInletVelocity - unexpected results

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

 February 8, 2013, 14:47 flowRateInletVelocity - unexpected results #1 Member   Chris Join Date: Aug 2012 Location: Calgary, Alberta, Canada Posts: 73 Rep Power: 4 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

February 10, 2013, 22:29
#2
Senior Member

Join Date: Nov 2009
Location: Michigan
Posts: 135
Rep Power: 7
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?

February 11, 2013, 05:04
fixedFluxPressure + flowRateInletVelocity to fix mass flow?
#3
New Member

Tatu Pinomaa
Join Date: Oct 2012
Location: Finland
Posts: 16
Rep Power: 4
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

 February 11, 2013, 06:23 #4 Senior Member   Olivier Join Date: Jun 2009 Location: France, grenoble Posts: 235 Rep Power: 9 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

 February 11, 2013, 13:05 #5 Member   Chris Join Date: Aug 2012 Location: Calgary, Alberta, Canada Posts: 73 Rep Power: 4 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.

 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 Ahmed OpenFOAM Running, Solving & CFD 9 June 22, 2011 18:59 pixie Main CFD Forum 1 August 20, 2009 08:02 Valle OpenFOAM Running, Solving & CFD 4 August 19, 2009 05:53 Cirion0000 CFX 4 July 18, 2009 12:48 shokry FLOW-3D 5 February 3, 2009 14:56

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