Question about Plotting Vapor Penetration Length

asakusa · October 30, 2017, 04:31

Hello!
I'm a beginner of OpenFOAM and I'm using sprayFoam(in OF4.1). I need to plot the liquid penetration length and vapor penetration length, comparing to the standard result to check if my model is correct. Liquid penetration length can be calculated by the function penetration.(fraction)(the pretty standard of fraction for liquid penetration is 0.95), as the following code in SprayCloud.C:
template<class CloudType>
void Foam::SprayCloud<CloudType>::info()
{
CloudType::info();
scalar d32 = 1.0e+6*this->Dij(3, 2);
scalar d10 = 1.0e+6*this->Dij(1, 0);
scalar dMax = 1.0e+6*this->Dmax();
scalar pen = this->penetration(0.95);

Info << " D10, D32, Dmax (mu) = " << d10 << ", " << d32
<< ", " << dMax << nl
<< " Liquid penetration 95% mass (m) = " << pen << endl;

}
But how can I calculate the vapor penetration of every time step? Should I still use this function? Thank you in advance.

rmn · August 18, 2020, 09:24

Hi asakusa,

have you found a solution back then?
I am at the same point as you have been.

I also struggle to identify which files and entries from the time directories I have to use for the plot.

I hope you (or anybody else) will still see this.
Thank you so much in advance.

Greetings

October 30, 2017, 04:31	Question about Plotting Vapor Penetration Length	#1
asakusa New Member FengHaifeng Join Date: May 2017 Posts: 4 Rep Power: 8	Hello! I'm a beginner of OpenFOAM and I'm using sprayFoam(in OF4.1). I need to plot the liquid penetration length and vapor penetration length, comparing to the standard result to check if my model is correct. Liquid penetration length can be calculated by the function penetration.(fraction)(the pretty standard of fraction for liquid penetration is 0.95), as the following code in SprayCloud.C: template<class CloudType> void Foam::SprayCloud<CloudType>::info() { CloudType::info(); scalar d32 = 1.0e+6this->Dij(3, 2); scalar d10 = 1.0e+6this->Dij(1, 0); scalar dMax = 1.0e+6*this->Dmax(); scalar pen = this->penetration(0.95); Info << " D10, D32, Dmax (mu) = " << d10 << ", " << d32 << ", " << dMax << nl << " Liquid penetration 95% mass (m) = " << pen << endl; } But how can I calculate the vapor penetration of every time step? Should I still use this function? Thank you in advance. jonghak5724@unist.ac.kr likes this.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
SixDoFRigidBodyMotion under OF2.3 ( self oscillating cylinder)	Scabbard	OpenFOAM Running, Solving & CFD	1	July 22, 2014 04:50
Quick Reference Value Question (Length)	kyles.sverige	FLUENT	0	May 9, 2014 15:54
Integral Length Scale vs Turbulent Length Scale	Scott Nordsen	FLUENT	2	January 20, 2014 01:17
Cases with small length scale work fine on a single processor but fail in parallel	adona058	OpenFOAM Bugs	5	April 17, 2009 04:41
Characteristic Length Scale	Thomas P. Abraham	Main CFD Forum	2	April 20, 1999 17:16

August 18, 2020, 09:24		#2
rmn New Member Roman Reinhardt Join Date: Jun 2020 Location: Germany Posts: 4 Rep Power: 5	Hi asakusa, have you found a solution back then? I am at the same point as you have been. I also struggle to identify which files and entries from the time directories I have to use for the plot. I hope you (or anybody else) will still see this. Thank you so much in advance. Greetings