Electrostatics Model

The Electrostatics contact model allows you to add electrostatic forces to the existing kinematic forces in the Particle Body Force section.

Electrostatic force is a force between particles exerted due to their charge, and can be either attractive or repulsive. Since the EDEM Electrostatic model uses Electrostatic screening, you must specify the size of a screen when using this model. Particles outside the screen will be removed from electrostatic force calculations on the central particle as shown:


The screened Coulomb force Fs is defined as: 

F s = q 1 q 2 4 π 0 k r + 1 r 2 e k r MathType@MTEF@5@5@+= feaahyart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOramaaBa aaleaacaWGZbaabeaakiabg2da9maalaaabaGaamyCamaaBaaaleaa caaIXaaabeaakiaadghadaWgaaWcbaGaaGOmaaqabaaakeaacaaI0a GaeqiWdaNaeyicI48aaSbaaSqaaiaaicdaaeqaaaaakmaabmaabaWa aSaaaeaacaWGRbaabaGaamOCaaaacqGHRaWkdaWcaaqaaiaaigdaae aacaWGYbWaaWbaaSqabeaacaaIYaaaaaaaaOGaayjkaiaawMcaaiaa dwgadaahaaWcbeqaaiabgkHiTiaadUgacaWGYbaaaaaa@4CC4@

Where q1 and q2 are the charges of two particles, r is the distance between their centers, and ε0 is the permittivity of free space (8.854 X 10-12 Fm-1). κ is an inverse of the Debye length, λD. 

k= q e 1 0 K B T i n i z i 2 MathType@MTEF@5@5@+= feaahyart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4Aaiabg2 da9iaadghadaWgaaWcbaGaamyzaaqabaGcdaGcaaqaamaabmaabaWa aSaaaeaacaaIXaaabaGaeyicI4SaeyicI48aaSbaaSqaaiaaicdaae qaaOGaam4samaaBaaaleaacaWGcbaabeaakiaadsfaaaWaaabuaeaa caWGUbWaaSbaaSqaaiaadMgaaeqaaOGaamOEamaaBaaaleaacaWGPb aabeaakmaaCaaaleqabaGaaGOmaaaaaeaacaWGPbaabeqdcqGHris5 aaGccaGLOaGaayzkaaaaleqaaaaa@4B41@
Where
  • qe is the charge of an electron (1.602 X 10-19C)
  • ε is the relative permittivity of the medium. EDEM uses a value of 1, which is a relative permittivity of a vacuum.
  • Kb is the Boltzmann’s constant equal to 1.38 X 10-23 JK-1.
  • T is the temperature in Kelvin, EDEM uses a value of 295.4K.
  • The last term is the local charge concentration (where n is the number of particles of charge z).