Spray simulation post process
 

Dear All,

may i know how to get spray tip penetration and cone angle and no of droplet distribution from paraview.

Thanks in advance

Balaji

Post-processing sprays in paraview gets a little tricky. I don't think there is a built in filter to get what you want; you'll have to write one yourself. Though isn't the cone angle an input parameter from the sprayProperties file?

For future reference, the KinematicCloud class (upon which the sprays are based) have a member function called penetration(scalar prc) that returns the penetration length based on the percentage prc of the current mass in the system. If you add this to the solver you can get the penetration at runtime, even write it to a logfile.

Dear Marco,

Since i am new to openfoam and paraview, could you please brief me how to write the code in order to get the spray penetration length and the droplet size distribution.

thanks
Balaji

Hi Balaji,

For the penetration length, you'll need to create/modify the solver you are using. Look at the tutorial on adding temperature to icoFoam for a quick overview of the process:

http://openfoamwiki.net/index.php/Ho...ure_to_icoFoam

In the main solver loop, you will need to add something like this:


You'll need to find the name of the spray object created in your solver, and you can use whatever percentage you want. Make the change, compile your solver and run it. You should see a message output every iteration. You will need to save the output to a logfile and extract the penetration info from there.

If you are willing to study more C++ and OF, you can make it so that the results are written to a file in 2 column format. As a lot of my textbooks used to say "This is left as an excercise for the reader".

Good luck!

Hi Balaji,

Were you able to get this to work? I am interested in adding spray penetration length code to the .../multiphase/interFoam/les/nozzleFlow2D tutorial but have no idea how to do it. Even if the spray penetration data is output to the log file that would be okay. Any advice/help? Thanks in advance!

Scott

Dear M. Turcios

First of all, I know this is a very old Thread, and OpenFOAM classes have changed since then, so I was wondering if you could help me. I am using OF301, and I would like to calculate Liquid Penetration Length in InterFoam.

Having implemented your penetration length explanation above, I found one message:

KinematicCloudI.H:404:37: error: ‘const parcelType’ has no member named ‘position0’
scalar d = mag(p.position() - p.position0());

I searched for this function position0(), and its declaration exists in SprayParcel.H, its definition in SprayParcelI.H. The latter uses one variable as a return value (position0_) which is declared and manipulated in SprayParcel.C. The problem with this file (SprayParcel.C) is that it includes unnecessary reacting and compressible features which I cannot get rid of without messing up position0_ calculation.

Any ideas would be highly appreciated
Thanks in advance

JAP