{Epot2D0.5ElTh.PDE  - Electric POTential  D.Bakewell 13:45 Tu/6/2/2001  }
title 'Epot2D0.5ElTh'
SELECT
  regrid=on  {enable adaptive grid - warning can take too long}
 errlim=1e-6 {8e-6 errorlimit is much faster }
COORDINATES 
  cartesian(x,y)
VARIABLES
  Phi   	{Electric potential}  Ex   {Flux in X direction}  Ey	{Flux in X direction}
  dxE2   {dy(E^2)}	  dyE2   { mgrdE2} {magnitude grad(E^2)}
  EMod
DEFINITIONS
  om=10^(0)
  W=10*om    {electrode width/gap} He=0.1*om {electrode boundary Height}
  x0=0*om    y0=0*om   { horizontal & vertical origin} H=40*om {External height}   
 {a=arbitrary}  b=y0+He  c=H
  V=1
EQUATIONS
 div[grad(Phi)]=0  		 {2D Laplace's equation }
 dx(Phi)=-Ex		 {for monitoring Electric field in X direction}
 dy(Phi)=-Ey		 {for monitoring Electric field in Y direction}
 dx(Ex*Ex + Ey*Ey)=dxE2
 dy(Ex*Ex + Ey*Ey)=dyE2 
 sqrt(Ex^2 + Ey^2)=EMod
 BOUNDARIES
 Region 1            { define region 2 RECTANGLE boundary on horizontal plane }
 start (x0,b) Value(Phi)=V line to (W/2,b) to (W/2,y0)
 Neumann(Phi)=0 line to (W,y0)
 Value(Phi)=0 line to (W,c)
 Neumann(Phi)=0 line to (x0,c)
 Neumann(Phi)=0 line to finish
 FEATURE
  {start (4.95,0.2) line to (5.05,0.2)} start (4.9,0.3) line to (5.1,0.3)
 start (4.75,0.6) line to (5.25,0.6) start (4,1) line to (6,1)
MONITORS 
for cycle=1  contour(Phi)  contour(Ex)  contour(Ey)
elevation(Phi) from (x0,y0) to (W,y0)
elevation(Phi) from (x0,c) to (W,c)
contour(EMod)
elevation(dxE2) from (x0,0.2) to (W,0.2)
elevation(dyE2) from (x0,0.2) to (W,0.2)
elevation(dxE2) from (4.9,0.2) to (4.9,c)
elevation(dyE2) from (4.9,0.2) to (4.9,c)
elevation(LOG10(-dyE2)) from (4.9,0.2) to (4.9,c)
HISTORIES
    {history(Phi) at (20*om,20*om)}
PLOTS
contour(Phi)  contour(Ex)  contour(Ey)
elevation(Phi) from (x0,y0) to (W,y0)
 {export format {"#t#b#"} "#x#b#y#b#" file="c:\windows\desktop\Phi.dat"}
elevation(Phi) from (x0,c) to (W,c)
contour(EMod)
elevation(dxE2) from (x0,0.2) to (W,0.2)
 {export format {"#t#b#"} "#x#b#y#b#" file="c:\windows\desktop\dxE2.dat"}
elevation(dyE2) from (x0,0.2) to (W,0.2)
 {export format {"#t#b#"} "#x#b#y#b#" file="c:\windows\desktop\dyE2.dat"}
elevation(dxE2) from (4.9,0.2) to (4.9,c)
 export format {"#t#b#"} "#x#b#y#b#" file="c:\windows\desktop\dxE2y.dat"
elevation(dyE2) from (4.9,0.2) to (4.9,c)
 export format {"#t#b#"} "#x#b#y#b#" file="c:\windows\desktop\dyE2y.dat"
elevation(LOG10(-dyE2)) from (4.9,0.2) to (4.9,c)
END  8745