* PhiView 3.1 - force analysis for a dielectric solution
* Field Precision LLC, Copyright 2010
* Variables: &Ex, &Ey, &Ez, &EMag, &Q[n]
* Stored quantities
*    Phi: electrostatic potential (V)
*    EpsiR: relative permitivity
*    Rho: space-charge density (C/m3)
* OPERATOR SET
*   &Phi, &Ex, &Ey, &Ez, &EMag, &Q[N]

HEADER
  ProgramName = PhiView 3.1
  ManualName = hiphi.pdf
  WebSiteURL = http://www.fieldp.com
  DataSuffix = HOU
  About = PhiView 3.0;Copyright 2009;Field Precision LLC;Albuquerque, NM USA
  NRecord = 3
  SliceDisplay = 5
  ScanDisplay = 4
  VectorDisplay = 1
  PlaneDisplay = 5
  SurfaceDisplay = 4
END

STOREDQUANT
  Phi N
  Epsi E
  Rho E
END

REGIONQUANT
END

PROGPARAM
  $Epsi0 = 8.854187E-12
* Origin for torque calculations, specify in meters
  $xt = 0.00
  $yt = 0.00
  $zt = 0.00
END

RUNPARAM
END

REALTIMEPARAM
END

INTERPOLATION
  Ex = &Ex
  Ey = &Ey
  Ez = &Ez
  |E| = &EMag
  Phi = &Phi
  EpsiR = &Q[2]
  Rho = &Q[3]
END

VECTOR
  EVector = &Ex;&Ey;&Ez
END

* Volumetric force components from a space-charge density
VOLUME
  FxVol = &Ex &Q[3] *
  FyVol = &Ey &Q[3] *
  FzVol = &Ez &Q[3] *
  TxVol = &Ez $y $yt - * &Ey $z $zt - * - &Q[3] *
  TyVol = &Ex $z $zt - * &Ez $x $xt - * - &Q[3] *
  TzVol = &Ey $x $xt - * &Ex $y $yt - * - &Q[3] *
END

* Maxwell stress tensor
SURFACE
*  FxSurf = s11 $Epsi0 *;s12 $Epsi0 *;s13 $Epsi0 *
*  FySurf = S21 $Epsi0 *;S22 $Epsi0 *;s23 $Epsi0 *
*  FzSurf = s31 $Epsi0 *;s32 $Epsi0 *;s33 $Epsi0 *
  FxSurf = &Ex 2 ^ &EMag 2 ^ 2.0 / - $Epsi0 *;&Ex &Ey * $Epsi0 *;&Ex &Ez * $Epsi0 *
  FySurf = &Ey &Ex * $Epsi0 *;&Ey 2 ^ &EMag 2 ^ 2.0 / - $Epsi0 *;&Ey &Ez * $Epsi0 *
  FzSurf = &Ez &Ex * $Epsi0 *;&Ez &Ey * $Epsi0 *;&Ez 2 ^ &EMag 2 ^ 2.0 / - $Epsi0 *
*  TxSurf =
*      $y $yt - S31 * $z $zt - S21 * - $Epsi0 *;
*      $y $yt - S32 * $z $zt - S22 * - $Epsi0 *;
*      $y $yt - S33 * $z $zt - S23 * - $Epsi0 *;
*  TySurf =
*      $z $zt - S11 * $x $xt - S31 * - $Epsi0 *;
*      $z $zt - S12 * $x $xt - S32 * - $Epsi0 *;
*      $z $zt - S13 * $x $xt - S33 * - $Epsi0 *;
*  TzSurf =
*      $x $xt - S21 * $y $yt - S11 * - $Epsi0 *;
*      $x $xt - S22 * $y $yt - S12 * - $Epsi0 *;
*      $x $xt - S23 * $y $yt - S13 * - $Epsi0 *;
  TxSurf = $y $yt - &Ez &Ex * * $z $zt - &Ey &Ex * * - $Epsi0 *;$y $yt - &Ez &Ey * * $z $zt - &Ey 2 ^ &EMag 2 ^ 2.0 / - * - $Epsi0 *;$y $yt - &Ez 2 ^ &EMag 2 ^ 2.0 / - * $z $zt - &Ey &Ez * * - $Epsi0 *
  TySurf = $z $zt - &Ex 2 ^ &EMag 2 ^ 2.0 / - * $x $xt - &Ez &Ex * * - $Epsi0 *;$z $zt - &Ex &Ey * * $x $xt - &Ez &Ey * * - $Epsi0 *;$z $zt - &Ex &Ez * * $x $xt - &Ez 2 ^ &EMag 2 ^ 2.0 / - * - $Epsi0 *
  TzSurf = $x $xt - &Ey &Ex * * $y $yt - &Ex 2 ^ &EMag 2 ^ 2.0 / - * - $Epsi0 *;$x $xt - &Ey 2 ^ &EMag 2 ^ 2.0 / - * $y $yt - &Ex &Ey * * - $Epsi0 *;$x $xt - &Ey &Ez * * $y $yt - &Ex &Ez * * - $Epsi0 *
END

ENDFILE

-------------------------------------------------------------------------------
                             STRESS TENSOR COMPONENTS
-------------------------------------------------------------------------------

s11: &Ex 2 ^ &EMag 2 ^ 2.0 / -
s12: &Ex &Ey *
s13: &Ex &Ez *

s21: &Ey &Ex *
s22: &Ey 2 ^ &EMag 2 ^ 2.0 / -
s23: &Ey &Ez *

s31: &Ez &Ex *
s32: &Ez &Ey *
s33: &Ez 2 ^ &EMag 2 ^ 2.0 / -