CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

Mass flow rate history over solution step- rhoSimpleFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   December 8, 2019, 10:20
Default Mass flow rate history over solution step- rhoSimpleFoam
Join Date: Sep 2017
Posts: 50
Rep Power: 8
gian93 is on a distinguished road
Hi Guys,

I'm simulating a simple steady state channel flow with a butterfly valve . I want to control the evolution of the computed mass flow rate at the outlet patch to evaluate the solution convergence (i also use the residual plot).

I used rhoSimpleFoam as solver.

Theese are the BC:

U >>> flowRateInletVelocity (I impose a mass flow rate of 0.5 kg/s)

p >>> zeroGradient


U >>> inletOutlet (in my case is the same as pressureNormalInletOutletVelocity)
P >>> fixedvalue (I impose a downstream pressure of 350000 Pa)

to calculate the mass flow rate at outlet patch I' ve added to my controlDict the functionObject for mass flow calculation :

type faceSource;
functionObjectLibs ("");
enabled true;
outputControl timeStep;
log true;
valueOutput false;

sourceName Outlet;

operation sum;



With this function object I would graph the mass flow rate at the outlet patch over each step iteration. I think that after few iteration the solver will find the solution and the mass flow rate at outlet will be equal (for mass conservation) to the inlet one (asymptotic trend). For the early time steps i would expect an oscillatory value of the variable phi (in this case mass Flow rate) that tends to stabilize when solution reach the convergence.

The problem is that this function object gives me a value of the outlet mass flow rate always equal (exactly equal) to the same value of the inlet patch without the expected initial oscillations...

Why i can't see the oscillations of the solver initial time step ?? There is something that i've missed from a phisical point of view in my BC or in the functionObject itself?
As it concerns inlet patch if i follow the Pressure evolution by an appropriate functionObject i can appreciate the "convergence trend".. but why it do not happen for mass flow rate at outlet?

gian93 is offline   Reply With Quote


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
Mass flow rate: calculation v/s computation beguxa FLUENT 5 December 2, 2018 21:02
Shape Optimization Problem of Mass Flow Rate ggc1991 SU2 Shape Design 1 December 1, 2016 20:24
Uniform mass flow rate at different outlets. mehcaecfd STAR-CCM+ 6 May 18, 2016 02:34
Periodic channel flow with time dependent mass flow rate QBeast FLUENT 3 May 10, 2013 13:14
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03

All times are GMT -4. The time now is 13:27.