Mass flow rate: calculation v/s computation
In my simulations density (Rho) is constant. I have computed mass flow rate (M) through an interior face using standard mass flow rate calculation option from fluent. Then I calculated mass averaged velocity (Vm), area weighted averaged velocity (Va) and facet averaged velocity (Vf). To recheck, I calculated the mass flow rate as M' = Rho*A*V. By using all the three above mentioned velocities in the above formula I am not able to get both mass flow rates the same (i.e. M' and M are not equal). Which is the apropriate value of velocity that one should consider to re-calculate accurately the mass flow rate?
I guess, the mass flow rate is computed in fluent as summation of (Rho*A*V) values at all the facets. Could any body throw some light on this issue?
Hi beguxa, I believe the velocity you want to use is actually the axial velocity. When running my current simulation of flow rate through a nozzle I used surface integrals in the report menu to get an accurate reading but when calculating the flow through use of a custom field function the velocity which returned the same value for flow rate was the axial velocity.
you can use standard massflow option in fluent. but to check this use velocity normal to the surface in your calculation. I think both should agree.
beguxa, There are subtle differences in all four values.
For correct mass flow rate, we want addition of all the mass flow rates going into the every surface element of your surface.
Let's compare all four approaches with this bold statement.
Area averaged velocity is calculated as . This is a good representation of average velocity over the surface. But depending on velocity distribution and mesh, smaller mesh elements may have higher velocities and vice versa, so this doesn't equate to bold sentence at the beginning.
Facet average is calculated as , where n is no. of facets in the surface. If your mesh is not uniform and if there are more mesh elements in regions of high velocity, your facet average velocity will be unnecessarily higher than the average velocity. You can't multiply this with area and still get the definition of bold statement at the beginning.
Mass weighted average is . We are slowly getting closer now, since there is a slight reference of actual mass going into every element - - in the denominator. But again since it is "weighted" it bears same limitations as above.
Mass flow rate is calculated as: . This definition is closest to the bold statement at the beginning. And hence is recommended.
Although, if mesh is of high quality, structured and if solution is fully converged, the differences between these should become smaller. But in real world, with often-used unstructured meshes, this is not always possible.
If you define a variable in the Define> Custom Field Functions... as Density*Velocity, and then integrated this variable in the report>surface integral using the "integral" in the "Report type" you can observe that this value is equal to the mass flow rate in your desired surface, which is reported by FLUENT post processing from Report> FLuxes
|All times are GMT -4. The time now is 11:36.|