we have tc and en:
(*) tc = pi 2 m_eff / sqrt( 4 kn m_eff - cn^2 ) (note your delta in
brackets)
en = exp( -cn tc / 2 / meff )
where we have kn and cn:
kn = m_eff / tc^2 * ( pi^2 + ln^2(en) )
cn = -2 * m_eff / tc * ln(en)
Then, for pair spheres with mass m1 we have
m1_eff = m1/2, en1, tc1 => kn1 and cn1.
For pair spheres with mass m2 we have
m2_eff = m2/2, en2, tc2 => kn2 and cn2.
For pair spheres with mass m1 and m2 we have
1) 1/m12_eff = 1/m1+1/m2, en12, tc12 => kn12, cn12,
OR we have
2) 1/kn12 = 1/kn1 + 1/kn2 and 1/cn12 = 1/cn1 + 1/cn2.
From (*) cn^2 can't be more than cn_crit^2 = (4 kn m_eff)...
OK, now I believe that I understand.
Does _OR_ means that 1) is equivalent to 2)?
Can you tell me more about the time step and the physical meaning of tc?
Thanks,
VR
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : [email protected]
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
