search

VasiliBaranov / packing-generation

Hard-sphere packing generation in C++ with the Lubachevsky–Stillinger, Jodrey–Tory, and force-biased algorithms and packing post-processing.
MIT License
106 stars 43 forks source link

meaning of insertion_radii.txt #13

Open est522 opened 4 years ago

est522 commented 4 years ago

Hi, I generated a monodisperse packing with your code and I'm interested in calculating the mean pore size within the packing. Does the insertion_radii.txt file contain the values of each pore size? Why are there negative values? I expected to find radii size (positive values). Best wishes!

VasiliBaranov commented 4 years ago

Hi,

Does the insertion_radii.txt file contain the values of each pore size?

Under some definition of pore sizes, yes. I.e. the program generates "pore centers" in the packing randomly and uniformly and calculates distances to the nearest sphere surfaces. The resulting distribution also sometimes called "nearest surface distribution".

Why are there negative values?

I write negative values for "pore centers" that appeared inside hard particles. I.e. these are the distances to the nearest surface from inside the hard particle (taken with the minus sign). So you can just filter out the negative values.

There was actually another project in this research group where we were calculating pore size distributions with different methods, so you may find it interesting: https://github.com/VasiliBaranov/PorousMediaAnalysis . The related paper is here: https://pubs.rsc.org/en/content/articlelanding/2016/nj/c5nj02814k#!divAbstract

Hope this helps.

Best, Vasili

est522 commented 4 years ago

Thank you very much!