CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Visualization & Post-Processing Software > ParaView

[OpenFOAM] Species mass flow rate

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

Like Tree2Likes
  • 1 Post By tomf
  • 1 Post By tomf

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 19, 2011, 16:03
Default Species mass flow rate
  #1
Member
 
N. A.
Join Date: May 2010
Posts: 64
Rep Power: 15
N. A. is on a distinguished road
Hello All,

I am trying to calculate the average mass flow rate of a particular species passing through a plane in the geometry, such as mass flow rate of CO2 passing through a plane.

Can you please let me know if know how we can extract this information in paraview.

Is there a way to calcualte this value directly in OpenFOAM

Thanks,
Nir
N. A. is offline   Reply With Quote

Old   November 14, 2012, 19:16
Default
  #2
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

your post is old and you just made four posts so I think none will answer me.

Do you find out how that works?

If someone can explain it, it would be very good for me and all the others


Tobi
Tobi is offline   Reply With Quote

Old   November 15, 2012, 05:16
Default
  #3
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hello Tobi,

How are you?

I think I may be able to help you on this one. The most easy way to do this with paraView would be to calculate the total mass flow through a plane using Y_CO2*rho*U_n where U_n is the velocity component normal to the plane of interest and Y_CO2 is the mass fraction of CO2. If you have a plane not aligned with x, y or z, you can use the calculator filter or the surface flow filter to get the velocity components perpendicular to that plane. Make sure the calculator for rho*U_n works on the cell data. Than go to filters>integrate Variables and you will get a table with the integral Y_CO2*rho*U*dA over Area A. This is the mass of CO2 passing through the area.

I do not have an example for doing this during the calculation, but I think something can be done with swak4Foam.

Regards,
Tom
kostnermo likes this.
tomf is offline   Reply With Quote

Old   November 15, 2012, 12:48
Default
  #4
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by tomf View Post
Hello Tobi,

How are you?

I think I may be able to help you on this one. The most easy way to do this with paraView would be to calculate the total mass flow through a plane using Y_CO2*rho*U_n where U_n is the velocity component normal to the plane of interest and Y_CO2 is the mass fraction of CO2. If you have a plane not aligned with x, y or z, you can use the calculator filter or the surface flow filter to get the velocity components perpendicular to that plane. Make sure the calculator for rho*U_n works on the cell data. Than go to filters>integrate Variables and you will get a table with the integral Y_CO2*rho*U*dA over Area A. This is the mass of CO2 passing through the area.

I do not have an example for doing this during the calculation, but I think something can be done with swak4Foam.

Regards,
Tom
Hi Tom,

I am fine. I have two weeks left to finish my masterthesis. Did you follow the libOpenSmoke thread? I find some fatal bugs in the tool and will bring out my updated model after I finished my thesis. So I think in 4 Weeks everything is avaiable (maybe + validation and description). To tool is very nice but there are some problems with the enthalpie. But I solved it


Thanks for your replay. I am too stupid
But I test it in the thermolibrary. therefor I ll create a new scalarField and will calculate the CO-emission in PPM directly. I think it should be possible and then I can compare it.

Tobi
Tobi is offline   Reply With Quote

Old   November 16, 2012, 04:32
Default
  #5
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Hi Tobi,

Yes I did follow the libOpenSmoke thread, looks like you have some interesting results!

Quote:
But I test it in the thermolibrary. therefor I ll create a new scalarField and will calculate the CO-emission in PPM directly. I think it should be possible and then I can compare it.
Sounds like a plan,

good luck
Tom
tomf is offline   Reply With Quote

Old   November 16, 2012, 07:46
Default
  #6
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Tom,

one question to you.

If I want the massweighted species in the whole fluid domain, is it possible to do like you said with mag(U) then?
Tobi is offline   Reply With Quote

Old   November 16, 2012, 10:54
Default
  #7
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi Tom,

just one question, if I make everything correct.


I wanna calculate CO in ppm on a patch.

Therefor I calculate the Mean molmass with:


M_aver = sum(Y_i / M_i) of CO, O2, H2O, N2, CO2, H2, OH

After that I calculate CO in ppm like:


COppm = M_aver/M_CO * Y_CO * 1.000.000

Well the temperature profil looks better with OpenFOAM compared with ansys but if I look to the CO emissions there is so a huge difference:


Measurement 80 ppm
Ansys 60 ppm
OF 0.0013 ppm


I am not sure if I make everything correct!
Tobi is offline   Reply With Quote

Old   November 19, 2012, 06:53
Default
  #8
Senior Member
 
Tom Fahner
Join Date: Mar 2009
Location: Breda, Netherlands
Posts: 634
Rep Power: 32
tomf will become famous soon enoughtomf will become famous soon enough
Send a message via MSN to tomf Send a message via Skype™ to tomf
Quote:
Originally Posted by Tobi View Post
Hi Tom,

just one question, if I make everything correct.


I wanna calculate CO in ppm on a patch.

Therefor I calculate the Mean molmass with:


M_aver = sum(Y_i / M_i) of CO, O2, H2O, N2, CO2, H2, OH

After that I calculate CO in ppm like:


COppm = M_aver/M_CO * Y_CO * 1.000.000

Well the temperature profil looks better with OpenFOAM compared with ansys but if I look to the CO emissions there is so a huge difference:


Measurement 80 ppm
Ansys 60 ppm
OF 0.0013 ppm


I am not sure if I make everything correct!
Hi Tobi,

I would think you have:

M_aver = sum(Y_i * M_i)/sum (Y_i) = sum(Y_i * M_i)

and I agree with your calculation of COppm.
Regards,
Tom
rksopenfoam likes this.

Last edited by tomf; November 19, 2012 at 06:54. Reason: typo
tomf is offline   Reply With Quote

Old   November 19, 2012, 08:11
Default
  #9
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by tomf View Post
Hi Tobi,

I would think you have:

M_aver = sum(Y_i * M_i)/sum (Y_i) = sum(Y_i * M_i)

and I agree with your calculation of COppm.
Regards,
Tom
Hi Tom,

thanks for your replay.
I made a validation of my COppm calculation like I mentioned in the last post and its correct, like you told me

In that case I realized a other bug in the thermodynamics
but in my thesis there is no time for it. Maybe its a problem with the flameletgeneration.

Tobi
Tobi is offline   Reply With Quote

Old   December 3, 2012, 18:45
Default
  #10
Senior Member
 
Wouter van der Meer
Join Date: May 2009
Location: Elahuizen, Netherlands
Posts: 203
Rep Power: 17
wouter is on a distinguished road
Hello,

Maybe a stupid thing to ask someone that is working on a thesis, but is ppm for gasses not volume based instead of mass based?

best regards
Wouter
wouter is offline   Reply With Quote

Old   December 3, 2012, 19:29
Default
  #11
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Quote:
Originally Posted by wouter View Post
Hello,

Maybe a stupid thing to ask someone that is working on a thesis, but is ppm for gasses not volume based instead of mass based?

best regards
Wouter
Hi,

ppm on (maybe) the species CO is based on the mol fraction [X_CO] * 1.000.000 = CO_ppm.

If you wanna get an information on a slice its better to make a mass weighted calculation to get the better result.

If you wanna get an information on a specified volume you should do a volume weighted calculation.

If I am wrong please correct me.

Greetings Tobi
Tobi is offline   Reply With Quote

Old   January 20, 2016, 13:03
Default
  #12
Senior Member
 
Thomas Oliveira
Join Date: Apr 2015
Posts: 114
Rep Power: 11
t.oliveira is on a distinguished road
Quote:
Originally Posted by N. A. View Post
calculate the average mass flow rate of a particular species passing through a plane in the geometry [...] Is there a way to calculate this value directly in OpenFOAM
I had a similar problem (integrate a variable over an internal plane) and managed to do so. Please refer to http://www.cfd-online.com/Forums/ope...tml#post581774 . You may have to use areaNormalAverage instead of areaNormalIntegrate.

Best wishes,

Thomas
t.oliveira 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
Match Pressure Inlet/Outlet Boundary Condition Mass Flow Rate MSchneid Fluent UDF and Scheme Programming 3 February 23, 2019 07:00
Shape Optimization Problem of Mass Flow Rate ggc1991 SU2 Shape Design 1 December 1, 2016 21:24
Multiphase flow - incorrect velocity on inlet Mike_Tom CFX 6 September 29, 2016 02:27
DPM total mass flow rate yow FLUENT 0 July 17, 2015 13:01
Periodic channel flow with time dependent mass flow rate QBeast FLUENT 3 May 10, 2013 14:14


All times are GMT -4. The time now is 08:41.