# How to compute total liquid volume of the whole domain

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

 April 21, 2010, 15:13 How to compute total liquid volume of the whole domain #1 Member   angel Join Date: May 2009 Location: Spain Posts: 46 Rep Power: 17 I finished an interDyMFoam run. The case has variations of liquid in the domain, i would like to discover it but at the same time, I would like to compute the total liquid volume at each time step (post-porcessing). Can someone suggest the way(s) to do this? Where i have to put the code from http://www.cfd-online.com/Forums/ope...le-domain.html Info << " Liquid: " << sum(mesh.V()*gamma) << endl; Thanks!

 April 22, 2010, 03:05 #2 Senior Member     Roman Thiele Join Date: Aug 2009 Location: Eindhoven, NL Posts: 374 Rep Power: 20 You can put the code into the main function, close to the end of the loop like before the solver writes the time information. After that you have to recompile the solver. __________________ ~roman

 April 22, 2010, 15:22 #3 Member   angel Join Date: May 2009 Location: Spain Posts: 46 Rep Power: 17 Thanks a lot, I have recompile interDyMFoam, and now i can see the total volume of fluid inside my domain as: Time = 0.05 Liquid: sum((V*alpha1)) [0 3 0 0 0 0 0] 0.033936 Time = 9.08938 Liquid: sum((V*alpha1)) [0 3 0 0 0 0 0] 0.0338348 Unfortunately, i still have another problem. I have some fluid, but my simulation is a Tank (wall) excite by a sinusoidal lateral excitation using InterDyMFoam. Anybody know what could be the reason? Any ideas? Thanks,

April 23, 2010, 02:16
code correction
#4
Senior Member

Roman Thiele
Join Date: Aug 2009
Location: Eindhoven, NL
Posts: 374
Rep Power: 20
Quote:
 Originally Posted by anmartin Thanks a lot, Time = 0.05 Liquid: sum((V*alpha1)) [0 3 0 0 0 0 0] 0.033936 Time = 9.08938 Liquid: sum((V*alpha1)) [0 3 0 0 0 0 0] 0.0338348
Sorry I forgot something. Of course the code must read
Code:
`Info << " Liquid: " << sum(mesh.V()*gamma).value() << endl;`
__________________
~roman