yade-users team mailing list archive

[Jan Stránský <question230574@xxxxxxxxxxxxxxxxxxxxx>]

Question #230574 on Yade changed:
https://answers.launchpad.net/yade/+question/230574

    Status: Open => Answered

Jan Stránský proposed the following answer:
Hello Nolan,

for that purpose, you would need to create new Material class, IPhys class,
Ip2 functor (combining two Materials into IPhys) and Law2 functor (actual
constitutive material law). It should not be very difficult, but I
understand you don't want to hide deep in the source code :-) What is the
time horizont you will need it (days, week, months)? If you are ok to have
in a few weeks time (let's say 3 :-), I can do it, but I don't have time
now..
cheers
Jan



2013/6/11 Nolan Dyck <question230574@xxxxxxxxxxxxxxxxxxxxx>

> New question #230574 on Yade:
> https://answers.launchpad.net/yade/+question/230574
>
> Hi,
>
> This is a follow up to a question I asked a few weeks ago. I am using yade
> to model randomly generated periodic domains of spherical void phase porous
> media (I simply interpret the solid spheres in the yade simulation as
> revolved cuts in CAD). So far I have been quite successful in modifying
> some of the python scripts to produce very nice looking periodic geometries
> that I have been able to conduct pore level CFD simulations on using
> commercial software. Now in order to move forward and start generating more
> physically realistic representations of graphite foam I need to be able to
> correctly model the force displacement relationships between bubbles.
>
> Unfortunately the force displacement relationship between 2 interacting
> bubbles is very non-linear. It is also unfortunate that none of the other
> interaction models that yade currently supports (such as the Hertz-Midlin
> contact theory) is applicable to bubbles. As explained in a paper called
> Film drainage and coalescence between deformable drops and bubbles by Chan
> et. al (2010) two interacting bubbles have an approximate
> force-displacement relationship of:
>
> D = (F/[2*PI*S])*log(F/[4*PI*S*{R1+R2}])
>
> Where D is the separation between the surfaces of the two bubbles, F is
> the normal reaction force between the bubbles, S is the surface tension of
> the fluid (a material constant), and R1 and R2 are the radii of the
> interacting bubbles. This equation is a highly simplified version of the
> one in the paper for my application. A big problem is that I need to solve
> for the force when given the displacement. This means the solution must be
> iterative.
>
> This post has become very long but I now arrive at my problem. I am still
> very new to Yade and python (and new to Unix for that matter) and I find
> the prospect of modifying Yade itself to implement the contact law very
> daunting. Can someone point me in the right direction on how to do this?
> What classes do I need to create? Which existing contact law should I use
> as a template? Is there a simple trick I can use in python to make this
> happen without diving into the c++ code? I'm just a master's student and
> this work is not really the core of my thesis so I can't really commit to
> working on this for very long. I would feel more inclined to do this if the
> community felt that this is a worthwhile addition to yade, and if someone
> would be willing to maintain the code after I am finished.
>
> Thanks,
> Nolan
>
> --
> You received this question notification because you are a member of
> yade-users, which is an answer contact for Yade.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~yade-users
> Post to     : yade-users@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~yade-users
> More help   : https://help.launchpad.net/ListHelp
>

-- 
You received this question notification because you are a member of
yade-users, which is an answer contact for Yade.