User:Tom Allen/TensorExperimentsR14

From J Wiki
Jump to: navigation, search


GTER1400C.jpg


NB. ... script experimentalR.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

rfEEMcn0=:-@FEMcvcv30
rfEEMcn1=:-@FEMcvcv31
rfEEMcn2=:-@FEMcvcv32

rfEEMcn=:(0|:])@(rfEEMcn0,rfEEMcn1,:rfEEMcn2)


NB. ... script experimentalR.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

rfPEMcn0=:(rfEEMcn1*HEMcn2)+-@(rfEEMcn2*HEMcn1)
rfPEMcn1=:(rfEEMcn2*HEMcn0)+-@(rfEEMcn0*HEMcn2)
rfPEMcn2=:(rfEEMcn0*HEMcn1)+-@(rfEEMcn1*HEMcn0)

rfPEMcn=:(0|:])@(rfPEMcn0,rfPEMcn1,:rfPEMcn2)


NB. ... script experimentalR.ijs ...

NB. ... these verbs are valid only in Cartesian coordinates ...

rfJEMcn  =:+/"1@((<1 2)|:])"3@(-@FEMcvcvdn)
rfJEMcndn=:+/"1@((<1 2)|:])"4@(-@FEMcvcvdndn)


NB. ... execute (ijx) ...

   gXRAll=:,:@(1 1 2 8 1 1 1 _1 0 0"_)

   xCpts=:_ 0.7 1.3 26,_ 0.4p1 0.6p1 26,_ 0 2p1 26,:_ 0 2p1 26"_

NB. ... for Maxwell's equations ...

   *./^:3"3 (0=]) ((2^_44)qtsz]) (gXR([(]+0|:])"2@:((0 1|:])"3)@(+/"1@((<0 1)|:])"5@R)xRef@hkxR03)]) ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

   *./^:3"3 (0=]) ((2^_44)qtsz]) (]+(0|:])+0 1|:])"3 (gXR([+/"1@((<0 1)|:])"5@B2kdc yRef@hkxR03)])   ((0;50"_)pTRandom xCpts)''
1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd (((0.8723"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160

pd (((1.0912"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160

pd (((1.3202"_)@]"0);]"0) pTsteps _ _0.00004 0.00004 160

pd 'pdf'
pd 'show'


JTER1400F.jpg
... electric field ...


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.00p1 0.00p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.25p1 0.25p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.40p1 0.40p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.50p1 0.50p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.60p1 0.60p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 0.75p1 0.75p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 0.4 2.0 100,_ 0.5p1 0.5p1 0,_ 0.25p1 0.25p1 0,:(_ 1.00p1 1.00p1 0"_)*10"_

pd (gXR(cx0"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd (((0.8723"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160

pd (((1.0912"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160

pd (((1.3202"_)@]"0);]"0) pTsteps _ _0.0002 0.0002 160

pd 'pdf'
pd 'show'


JTER1401E.jpg
... charge distribution ...


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfEEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER1402D.jpg
... components of electric field ...


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 HEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER1403D.jpg
... components of magnetic field ...


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(0{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(1{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 500,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(2{])"1@(}:"1@}:"2@xdy smx"2 1 rfPEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER1404D.jpg
... components of Poynting vector ...


load 'plot'

gXRAll=:,:@(0.05 1 2 8 1 1 1 _0.1 0 0"_)

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 0p1 0p1 0"_

pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 2p1 2p1 0"_

pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

xCpts=:_ 1.0912 1.0912 0,_ 0.5p1 0.5p1 0,_ 0.25p1 2.25p1 100,:_ 4p1 4p1 0"_

pd (gXR(cx2"_ 1;[(3{])"1@(xdy smx"2 1 rfJEMcn)yRef@hkxR03)]) pTtotal xCpts''

pd 'pdf'
pd 'show'


JTER1405D.jpg
... charge distribution at times x3=0p1 x3=2p1 x3=4p1 ...



Download: File:LTER1400C.txt