Hi, I take advantage of Christian's recent commit (*) to re-discuss a little this validity range of capillary files. In the end, I see only one validity issue that comes from the precision with which the maximum suction value (the one above which there is no more solution for Laplace's equation) could be approximated during the construction of this files.
E.g. for contacting monosized spheres the maximum succion value was found indeed to be 1493 (~ the 1500 value that was previously spotted by trial and error) : see l. 686 of M(r=1) file, (suction is the second column).... Whereas for contacting spheres with r=2, it has been found as 2849 (see l.616 of M(r=2) file..) This corresponds in both cases to smallest filling angle values around 2 degrees (two last columns), whereas probably the theoretical maximum value of suction here is infinite, and the theoretical minimum filling angle should be 0. If I am right, I guess the same issue arises for non contacting spheres. For monosized spheres at an adimensionned distance of 0.1, the maximum suction value in the capillary file is 8.7, corresponding to a lowest filling angle value of 10.6 deg (l. 12876 of M(r=1) file...). I guess these values only approximate the exact limit values of maximum suction and minimum filling angle, but the exact suction limit value is probably far from 1500 in this distant case. Do you share my understanding ? If yes, there would be only one problem in our simulations : that the code can not handle correctly the cases (existing in reality) falling between the maximum suction value written in the files (that depends on the radii ratio and distance), and the exact maximum suction value. Is it really a big deal ? (If not, it is finally maybe not catastrophic to run simulations with uc* > e.g. 1500, and Christian was right to lower the "noise" of the error message) Jerome (*) https://github.com/yade/trunk/commit/569aa859ab7e9ca8e9ac1a99b445884ad1f3f232 ________________________________ From: Yade-dev [[email protected]] on behalf of Bruno Chareyre [[email protected]] Sent: November-12-14 6:43 AM To: [email protected] Subject: Re: [Yade-dev] Some cleaning in capillary law : data files validity On 10/11/14 20:00, Jerome Duriez wrote: Ok. I'm getting zero water volume for contacting spheres for an adimensionned capillary pressure Pc* = 1500 +/- 5 (the minimum adimensionned water volume seems to be ~ 2.6e-6) Great. It means that 1500 was probably the max value when generating the data files. There could be a warning in the code if Pc*>=1500, but it is only relevant for a specific set of data files... The max should be defined automatically when reading the data. B
_______________________________________________ Mailing list: https://launchpad.net/~yade-dev Post to : [email protected] Unsubscribe : https://launchpad.net/~yade-dev More help : https://help.launchpad.net/ListHelp
