* EStat 7.0 - electrostatic force in dielectric solutions
* Field Precision LLC, Copyright 2010
* Variables: &Phi, &Exz, &Eyr, &EMag, &Exz0, &Eyr0, &EpsSig1, &EpsSig2, &Rho
* Stored quantities
*    Phi: electrostatic potential (V)
*    EpsiRUp: relative permitivity in upper triangle
*    EpsiRDn: relative permitivity in lower triangle
*    RhoUp: space-charge density in upper triangle (C/m3)
*    RhoDn: space-charge density in lower triangle (C/m3)
* Region quantities
*  $RegEps1
*  $RegRho
*  $RegPot
*  $RegEps2
*  $RegAng
*  $RegSinAng
*  $RegCosAng

PROGPARAM
  $Epsi0 = 8.854187E-12
END

RUNPARAM
  $CondFlag = 0.0
END

REALTIMEPARAM
* Origin for torque calculations (in meters)
*  $xt = 0.00
*  $yt = 0.00
END

INTERPOLATION
  Phi = &Phi
  Ex/Ez = &Exz
  Ey/Er = &Eyr
  |E| = &EMag
END

VECTOR
  EVect = &Exz;&Eyr
  DVect = &Exz0 $Epsi0 *;&Eyr0 $Epsi0 *
END

VOLUME
* Volumetric force integrals for solutions with space charge
  FxVol/FzVol = &Exz &Rho *
  FyVol/FrVol = &Eyr &Rho *

* Torque is meaningful only in planar solutions
*  TzVol = $x $xt - &Eyr * $y $yt - &Exz * - &Rho *
END

SURFACE
*  FxSurf = s11 $Epsi0 *;s12 $Epsi0 *;s13 $Epsi0 *
*  FySurf = S21 $Epsi0 *;S22 $Epsi0 *;s23 $Epsi0 *
  FxSurf = &Exz 2 ^ &Eyr 2 ^ - $Epsi0 * 0.5 *;&Exz &Eyr * $Epsi0 *
  FySurf = &Eyr &Exz * $Epsi0 *;&Eyr 2 ^ &Exz 2 ^ - $Epsi0 * 0.5 *
* Torque is meaningful only in planar solutions
*  TzSurf =
*      $x $xt - S21 * $y $yt - S11 * - $Epsi0 *;
*      $x $xt - S22 * $y $yt - S12 * - $Epsi0 *;
*      $x $xt - S23 * $y $yt - S13 * - $Epsi0 *;
*  TzSurf = $x $xt - &Eyr &Exz * * $y $yt - &Exz 2 ^ &EMag 2 ^ 2.0 / - * - $Epsi0 *;$x $xt - &Eyr 2 ^ &EMag 2 ^ 2.0 / - * $y $yt - &Exz &Eyr * * - $Epsi0 *
END

ENDFILE