Re: [Yade-dev] Capillary scripts

2016-09-21 Thread Jerome Duriez
See (for reference) 
https://github.com/yade/trunk/commit/e23a181f023956bd195be89d44632748403f1ddb 
and previous commit, thanks !


[https://avatars1.githubusercontent.com/u/5427081?v=3=200]<https://github.com/yade/trunk/commit/e23a181f023956bd195be89d44632748403f1ddb>

Capillary scripts commit (http://www.mail-archive.com/yade-dev@lists ยท 
yade/trunk@e23a181<https://github.com/yade/trunk/commit/e23a181f023956bd195be89d44632748403f1ddb>
github.com
...launchpad.net/msg12103.html). Located in a new 
examples/capillaryLaplaceYoung folder intended to illustrates capillary Law2. 
With one new simulation script example and another old one, moved here





--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367



From: Yade-dev <yade-dev-bounces+jerome.duriez=ucalgary...@lists.launchpad.net> 
on behalf of Jerome Duriez <jerome.dur...@ucalgary.ca>
Sent: September-08-16 6:02 PM
To: Yade dev
Subject: [Yade-dev] Capillary scripts


Hi there,


Attached are three .m scripts aimed at generating capillary files for use with 
Law2_ScGeom_CapillaryPhys_Capillarity. Obviously, the idea, discussed 
previously here and there on the mailing lists, is to make this part of YADE 
more "open-source". Also, there is two major differences with previous (L. 
Scholtes') work:


* these scripts may generate capillary files accommodating any contact angle 
value, instead of 0 previously


* they provide new output as an "orientation tensor" defined from the 
fluid-fluid interface of capillary bridges. These data would constitute new 
attributes of CapillaryPhys, see "nn11" and "nn33" in the attached version of 
CapillaryPhys.hpp. (Note that these nn11 and nn33 can be used to get the bridge 
fluid-fluid interface S as S = 2*nn11 + nn33). This is about post-processing 
only.


The attached README.txt details the role of the .m files and how to use them: 
interested users in the end just need to give a look to, and launch, 
writesCapFile.m.



The core of the work, as usual, is solving the Laplace-Young equation for a 
capillary bridge configuration. The corresponding, quite classical, equations 
appear in solveLiqBridge file/function. There will be a dedicated journal paper 
on this topic, currently under "Minor Revision" review.



These files are ready to be pushed to trunk, together with updated 
CapillaryPhys (see attached), and Law2_ScGeom_CapillaryPhys_Capillarity, in 
order to read the new data (nn11 and nn33) from generated capillary files.


This email is for possible discussion before I perform the commit.


Cheers,


Jerome


PS: Vaclav, I took the liberty cc-ing you since I had once the feeling you 
would be interested in such scripts. (Woo does not include any capillary 
interactions at the moment, though ?)



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367
___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-dev] Capillary scripts

2016-09-14 Thread Jerome Duriez
Yes, they were coded by myself, indeed.



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367



From: Anton Gladky <gladky.an...@gmail.com>
Sent: September-14-16 1:14 PM
To: Jerome Duriez
Cc: Yade dev
Subject: Re: [Yade-dev] Capillary scripts

Hi Jerome,

from my point of view, it is not a problem to push those scripts to
the trunk as they are just scripts and you think they can be
useful to other user too.

Of course, they need to get a description in documentation or
(and) have a good example, how to use them. Yes, they can
be pushed to the scripts/ or to examples/corresponding_subfolder.

The most important question is the license. Are those scripts created
by you?

Best regards,

Anton


2016-09-14 3:19 GMT+02:00 Jerome Duriez <jerome.dur...@ucalgary.ca>:
> Hello,
>
>
> After a couple of private emails with Bruno, I can confirm that these
> scripts work with GNU Octave as well (in addition to the less open-source
> friendly MATLAB )
>
>
> We mentionned with Bruno the possibility uploading these .m scripts on the
> wiki, however I would greatly prefer they are in trunk, in order to enjoy
> the revision tracking tools.
>
> As these scripts are not bound, from my point of view, to any non
> open-source tool, I think it should be fine. Any opinion ? (Bruno ? others
> ?)
>
>
> Another question is what would be the best trunk subfolder for these files ?
> trunk/scripts ?
___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-dev] Capillary scripts

2016-09-14 Thread Anton Gladky
Hi Jerome,

from my point of view, it is not a problem to push those scripts to
the trunk as they are just scripts and you think they can be
useful to other user too.

Of course, they need to get a description in documentation or
(and) have a good example, how to use them. Yes, they can
be pushed to the scripts/ or to examples/corresponding_subfolder.

The most important question is the license. Are those scripts created
by you?

Best regards,

Anton


2016-09-14 3:19 GMT+02:00 Jerome Duriez :
> Hello,
>
>
> After a couple of private emails with Bruno, I can confirm that these
> scripts work with GNU Octave as well (in addition to the less open-source
> friendly MATLAB )
>
>
> We mentionned with Bruno the possibility uploading these .m scripts on the
> wiki, however I would greatly prefer they are in trunk, in order to enjoy
> the revision tracking tools.
>
> As these scripts are not bound, from my point of view, to any non
> open-source tool, I think it should be fine. Any opinion ? (Bruno ? others
> ?)
>
>
> Another question is what would be the best trunk subfolder for these files ?
> trunk/scripts ?

___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-dev] Capillary scripts

2016-09-14 Thread Jerome Duriez
Hi Klaus,


You're right about Python, the reason why I used Matlab/Octave instead is quite 
stupid: I'm much more at ease with these softwares...


I will give a close look to what's inside "examples" to figure out which folder 
to use, thanks for suggestion !


Jerome



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367



From: Yade-dev <yade-dev-bounces+jerome.duriez=ucalgary...@lists.launchpad.net> 
on behalf of Klaus Thoeni <klaus.tho...@gmail.com>
Sent: September-13-16 8:58 PM
To: yade-dev@lists.launchpad.net
Subject: Re: [Yade-dev] Capillary scripts

Hi Jerome,

great addition and thanks for sharing the scripts.

One question, why did you use Matlab/Octave and not python? Python scripts
would "fit" better into the trunk ;-)

I terms of where to upload the scripts. I would suggest a folder in
"examples" with an example, your scripts and instructions. I think some people
already used this approach. But anyway, this is just a suggestion.

Cheers
Klaus

On Wed, 14 Sep 2016 01:19:16 AM Jerome Duriez wrote:
> Hello,
>
>
> After a couple of private emails with Bruno, I can confirm that these
> scripts work with GNU Octave as well (in addition to the less open-source
> friendly MATLAB )
>
>
> We mentionned with Bruno the possibility uploading these .m scripts on the
> wiki, however I would greatly prefer they are in trunk, in order to enjoy
> the revision tracking tools.
>
> As these scripts are not bound, from my point of view, to any non
> open-source tool, I think it should be fine. Any opinion ? (Bruno ? others
> ?)
>
>
> Another question is what would be the best trunk subfolder for these files ?
> trunk/scripts ?
>
>
>
> --
>
> Jerome Duriez, Research Associate
>
> University of Calgary, Dpt of Civil Engineering
>
> +1 403 220 7367
>
>
> 
> From: Yade-dev
> <yade-dev-bounces+jerome.duriez=ucalgary...@lists.launchpad.net> on behalf
> of Jerome Duriez <jerome.dur...@ucalgary.ca> Sent: September-08-16 6:02 PM
> To: Yade dev
> Subject: [Yade-dev] Capillary scripts
>
>
> Hi there,
>
>
> Attached are three .m scripts aimed at generating capillary files for use
> with Law2_ScGeom_CapillaryPhys_Capillarity. Obviously, the idea, discussed
> previously here and there on the mailing lists, is to make this part of
> YADE more "open-source". Also, there is two major differences with previous
> (L. Scholtes') work:
>
>
> * these scripts may generate capillary files accommodating any contact angle
> value, instead of 0 previously
>
>
> * they provide new output as an "orientation tensor" defined from the
> fluid-fluid interface of capillary bridges. These data would constitute new
> attributes of CapillaryPhys, see "nn11" and "nn33" in the attached version
> of CapillaryPhys.hpp. (Note that these nn11 and nn33 can be used to get the
> bridge fluid-fluid interface S as S = 2*nn11 + nn33). This is about
> post-processing only.
>
>
> The attached README.txt details the role of the .m files and how to use
> them: interested users in the end just need to give a look to, and launch,
> writesCapFile.m.
>
>
>
> The core of the work, as usual, is solving the Laplace-Young equation for a
> capillary bridge configuration. The corresponding, quite classical,
> equations appear in solveLiqBridge file/function. There will be a dedicated
> journal paper on this topic, currently under "Minor Revision" review.
>
>
>
> These files are ready to be pushed to trunk, together with updated
> CapillaryPhys (see attached), and Law2_ScGeom_CapillaryPhys_Capillarity, in
> order to read the new data (nn11 and nn33) from generated capillary files.
>
>
> This email is for possible discussion before I perform the commit.
>
>
> Cheers,
>
>
> Jerome
>
>
> PS: Vaclav, I took the liberty cc-ing you since I had once the feeling you
> would be interested in such scripts. (Woo does not include any capillary
> interactions at the moment, though ?)
>
>
>
> --
>
> Jerome Duriez, Research Associate
>
> University of Calgary, Dpt of Civil Engineering
>
> +1 403 220 7367


___
Mailing list: https://launchpad.net/~yade-dev
Yade developers in Launchpad<https://launchpad.net/~yade-dev>
launchpad.net
yade-dev@lists.launchpad.net Policy: You must be a team member to subscribe to 
the team mailing list. View public archive View subscribers


Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp
___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-dev] Capillary scripts

2016-09-13 Thread Klaus Thoeni
Hi Jerome,

great addition and thanks for sharing the scripts.

One question, why did you use Matlab/Octave and not python? Python scripts 
would "fit" better into the trunk ;-)

I terms of where to upload the scripts. I would suggest a folder in  
"examples" with an example, your scripts and instructions. I think some people 
already used this approach. But anyway, this is just a suggestion. 

Cheers
Klaus

On Wed, 14 Sep 2016 01:19:16 AM Jerome Duriez wrote:
> Hello,
> 
> 
> After a couple of private emails with Bruno, I can confirm that these
> scripts work with GNU Octave as well (in addition to the less open-source
> friendly MATLAB )
> 
> 
> We mentionned with Bruno the possibility uploading these .m scripts on the
> wiki, however I would greatly prefer they are in trunk, in order to enjoy
> the revision tracking tools.
> 
> As these scripts are not bound, from my point of view, to any non
> open-source tool, I think it should be fine. Any opinion ? (Bruno ? others
> ?)
> 
> 
> Another question is what would be the best trunk subfolder for these files ?
> trunk/scripts ?
> 
> 
> 
> --
> 
> Jerome Duriez, Research Associate
> 
> University of Calgary, Dpt of Civil Engineering
> 
> +1 403 220 7367
> 
> 
> 
> From: Yade-dev
> <yade-dev-bounces+jerome.duriez=ucalgary...@lists.launchpad.net> on behalf
> of Jerome Duriez <jerome.dur...@ucalgary.ca> Sent: September-08-16 6:02 PM
> To: Yade dev
> Subject: [Yade-dev] Capillary scripts
> 
> 
> Hi there,
> 
> 
> Attached are three .m scripts aimed at generating capillary files for use
> with Law2_ScGeom_CapillaryPhys_Capillarity. Obviously, the idea, discussed
> previously here and there on the mailing lists, is to make this part of
> YADE more "open-source". Also, there is two major differences with previous
> (L. Scholtes') work:
> 
> 
> * these scripts may generate capillary files accommodating any contact angle
> value, instead of 0 previously
> 
> 
> * they provide new output as an "orientation tensor" defined from the
> fluid-fluid interface of capillary bridges. These data would constitute new
> attributes of CapillaryPhys, see "nn11" and "nn33" in the attached version
> of CapillaryPhys.hpp. (Note that these nn11 and nn33 can be used to get the
> bridge fluid-fluid interface S as S = 2*nn11 + nn33). This is about
> post-processing only.
> 
> 
> The attached README.txt details the role of the .m files and how to use
> them: interested users in the end just need to give a look to, and launch,
> writesCapFile.m.
> 
> 
> 
> The core of the work, as usual, is solving the Laplace-Young equation for a
> capillary bridge configuration. The corresponding, quite classical,
> equations appear in solveLiqBridge file/function. There will be a dedicated
> journal paper on this topic, currently under "Minor Revision" review.
> 
> 
> 
> These files are ready to be pushed to trunk, together with updated
> CapillaryPhys (see attached), and Law2_ScGeom_CapillaryPhys_Capillarity, in
> order to read the new data (nn11 and nn33) from generated capillary files.
> 
> 
> This email is for possible discussion before I perform the commit.
> 
> 
> Cheers,
> 
> 
> Jerome
> 
> 
> PS: Vaclav, I took the liberty cc-ing you since I had once the feeling you
> would be interested in such scripts. (Woo does not include any capillary
> interactions at the moment, though ?)
> 
> 
> 
> --
> 
> Jerome Duriez, Research Associate
> 
> University of Calgary, Dpt of Civil Engineering
> 
> +1 403 220 7367


___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


Re: [Yade-dev] Capillary scripts

2016-09-13 Thread Jerome Duriez
Hello,


After a couple of private emails with Bruno, I can confirm that these scripts 
work with GNU Octave as well (in addition to the less open-source friendly 
MATLAB )


We mentionned with Bruno the possibility uploading these .m scripts on the 
wiki, however I would greatly prefer they are in trunk, in order to enjoy the 
revision tracking tools.

As these scripts are not bound, from my point of view, to any non open-source 
tool, I think it should be fine. Any opinion ? (Bruno ? others ?)


Another question is what would be the best trunk subfolder for these files ? 
trunk/scripts ?



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367



From: Yade-dev <yade-dev-bounces+jerome.duriez=ucalgary...@lists.launchpad.net> 
on behalf of Jerome Duriez <jerome.dur...@ucalgary.ca>
Sent: September-08-16 6:02 PM
To: Yade dev
Subject: [Yade-dev] Capillary scripts


Hi there,


Attached are three .m scripts aimed at generating capillary files for use with 
Law2_ScGeom_CapillaryPhys_Capillarity. Obviously, the idea, discussed 
previously here and there on the mailing lists, is to make this part of YADE 
more "open-source". Also, there is two major differences with previous (L. 
Scholtes') work:


* these scripts may generate capillary files accommodating any contact angle 
value, instead of 0 previously


* they provide new output as an "orientation tensor" defined from the 
fluid-fluid interface of capillary bridges. These data would constitute new 
attributes of CapillaryPhys, see "nn11" and "nn33" in the attached version of 
CapillaryPhys.hpp. (Note that these nn11 and nn33 can be used to get the bridge 
fluid-fluid interface S as S = 2*nn11 + nn33). This is about post-processing 
only.


The attached README.txt details the role of the .m files and how to use them: 
interested users in the end just need to give a look to, and launch, 
writesCapFile.m.



The core of the work, as usual, is solving the Laplace-Young equation for a 
capillary bridge configuration. The corresponding, quite classical, equations 
appear in solveLiqBridge file/function. There will be a dedicated journal paper 
on this topic, currently under "Minor Revision" review.



These files are ready to be pushed to trunk, together with updated 
CapillaryPhys (see attached), and Law2_ScGeom_CapillaryPhys_Capillarity, in 
order to read the new data (nn11 and nn33) from generated capillary files.


This email is for possible discussion before I perform the commit.


Cheers,


Jerome


PS: Vaclav, I took the liberty cc-ing you since I had once the feeling you 
would be interested in such scripts. (Woo does not include any capillary 
interactions at the moment, though ?)



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367
___
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help   : https://help.launchpad.net/ListHelp


[Yade-dev] Capillary scripts

2016-09-08 Thread Jerome Duriez
Hi there,


Attached are three .m scripts aimed at generating capillary files for use with 
Law2_ScGeom_CapillaryPhys_Capillarity. Obviously, the idea, discussed 
previously here and there on the mailing lists, is to make this part of YADE 
more "open-source". Also, there is two major differences with previous (L. 
Scholtes') work:


* these scripts may generate capillary files accommodating any contact angle 
value, instead of 0 previously


* they provide new output as an "orientation tensor" defined from the 
fluid-fluid interface of capillary bridges. These data would constitute new 
attributes of CapillaryPhys, see "nn11" and "nn33" in the attached version of 
CapillaryPhys.hpp. (Note that these nn11 and nn33 can be used to get the bridge 
fluid-fluid interface S as S = 2*nn11 + nn33). This is about post-processing 
only.


The attached README.txt details the role of the .m files and how to use them: 
interested users in the end just need to give a look to, and launch, 
writesCapFile.m.



The core of the work, as usual, is solving the Laplace-Young equation for a 
capillary bridge configuration. The corresponding, quite classical, equations 
appear in solveLiqBridge file/function. There will be a dedicated journal paper 
on this topic, currently under "Minor Revision" review.



These files are ready to be pushed to trunk, together with updated 
CapillaryPhys (see attached), and Law2_ScGeom_CapillaryPhys_Capillarity, in 
order to read the new data (nn11 and nn33) from generated capillary files.


This email is for possible discussion before I perform the commit.


Cheers,


Jerome


PS: Vaclav, I took the liberty cc-ing you since I had once the feeling you 
would be interested in such scripts. (Woo does not include any capillary 
interactions at the moment, though ?)



--

Jerome Duriez, Research Associate

University of Calgary, Dpt of Civil Engineering

+1 403 220 7367
/*
*  Copyright (C) 2006 by luc Scholtes*
*  luc.schol...@hmg.inpg.fr  *
**
*  This program is free software; it is licensed under the terms of the  *
*  GNU General Public License v2 or later. See file LICENSE for details. *
*/
#pragma once
#include
#include
#include

class CapillaryPhys : public FrictPhys
{
	public :
		int currentIndexes [4]; // used for faster interpolation (stores previous positions in tables)
		
		virtual ~CapillaryPhys() {};

	YADE_CLASS_BASE_DOC_ATTRS_INIT_CTOR_PY(CapillaryPhys,FrictPhys,"Physics (of interaction) for :yref:`Law2_ScGeom_CapillaryPhys_Capillarity`.",
 ((bool,meniscus,false,Attr::readonly,"True when a meniscus with a non-zero liquid volume (:yref:`vMeniscus`) has been computed for this interaction"))
 ((bool,isBroken,false,,"Might be set to true by the user to make liquid bridge inactive (capillary force is zero)"))
 ((Real,capillaryPressure,0.,,"Value of the capillary pressure Uc. Defined as Ugas-Uliquid, obtained from :yref:`corresponding Law2 parameter`"))
 ((Real,vMeniscus,0.,,"Volume of the meniscus"))
 ((Real,Delta1,0.,,"Defines the surface area wetted by the meniscus on the smallest grains of radius R1 (R1