Why is tribocharging model not matching the reference paper?

satri
satri Altair Community Member
edited August 14 in Community Q&A

I am trying to find out why the charge transfer in Altair does not match the equations outlined in the reference paper given?

I am trying to model 1 particle falling on a plane due to gravity

the particle and geometry got charged to the same level which was expected. However, the value does not match the paper given in the documentation of Altair. please tell me why?

image

 

Tagged:

Answers

  • satri
    satri Altair Community Member
    edited July 12

    This is the calculation based on the equation given in the paper

    image

    EDEM gives a charge of 2.37E-15nc on particles where, as the equation gives 3.7e-13nc. Why is there such a huge difference? I am only giving Hertz Mindlin contact model, no rolling friction too

  • satri
    satri Altair Community Member
    edited July 12

    image

    I try to stop the simulation aat the time it contacts the ground. Even then the value computed from equation is atleast 100 times more than Altair

  • Stephen Cole
    Stephen Cole
    Altair Employee
    edited August 14
    Shyam Prasad V Atri said:

    image

    I try to stop the simulation aat the time it contacts the ground. Even then the value computed from equation is atleast 100 times more than Altair

    Hi Shyam,

    The code that EDEM uses is:

        double area = 4.0 * PI * (nPhysicalCurvature1 * nPhysicalCurvature1);     double q_s = 2.66e-5 * area;      *chargeMovedToSurface1 = m_modelData.getAlpha() * (q_s - fabs(charge1)) * nTimestep;

    Where nPhysicalCurvature1 is the radius of the smallest element in contact.

    Regards

    Stephen

  • satri
    satri Altair Community Member
    edited August 14
    Stephen Cole_21117 said:

    Hi Shyam,

    The code that EDEM uses is:

        double area = 4.0 * PI * (nPhysicalCurvature1 * nPhysicalCurvature1);     double q_s = 2.66e-5 * area;      *chargeMovedToSurface1 = m_modelData.getAlpha() * (q_s - fabs(charge1)) * nTimestep;

    Where nPhysicalCurvature1 is the radius of the smallest element in contact.

    Regards

    Stephen

    Hey, yes, I want the information of the sphere to be used, which is correct. I am using the same information to compute in Matlab, too, where I use the equation. So why are the information not matching?

    MATLAB evaluation and EDEM evaluation should be same. i am trying to understand where the difference is coming from.

Categories