doc/cookbook/sandbox.py

from esys.escript import *
from esys.escript.linearPDEs import LinearPDE
from esys.finley import Rectangle
from esys.pyvisi import Scene, DataCollector, Map, Camera 
from esys.pyvisi.constant import *
import os

#... set some parameters ...
xc=[0.0,0.0]
qc=50.e6
Tref=0.
rhocp=1
eta=10
kappa=1
tend=5.
# ... time, time step size and counter ...
t=0
h=0.1
i=0

#... generate domain ...
mydomain = Rectangle(l0=5,l1=0.00,n0=500, n1=1)
#... open PDE ...
mypde=LinearPDE(mydomain)
mypde.setSymmetryOn()
mypde.setValue(A=kappa*kronecker(mydomain),D=rhocp/h,d=eta,y=eta*Tref)
# ... set heat source: ....
x=mydomain.getX()
print x
#qH=qc*whereNegative(length(x-xc)-r)
1	from esys.escript import *
2	from esys.escript.linearPDEs import LinearPDE
3	from esys.finley import Rectangle
4	from esys.pyvisi import Scene, DataCollector, Map, Camera
5	from esys.pyvisi.constant import *
6	import os
7
8	#... set some parameters ...
9	xc=[0.0,0.0]
10	qc=50.e6
11	Tref=0.
12	rhocp=1
13	eta=10
14	kappa=1
15	tend=5.
16	# ... time, time step size and counter ...
17	t=0
18	h=0.1
19	i=0
20
21	#... generate domain ...
22	mydomain = Rectangle(l0=5,l1=0.00,n0=500, n1=1)
23	#... open PDE ...
24	mypde=LinearPDE(mydomain)
25	mypde.setSymmetryOn()
26	mypde.setValue(A=kappakronecker(mydomain),D=rhocp/h,d=eta,y=etaTref)
27	# ... set heat source: ....
28	x=mydomain.getX()
29	print x
30	#qH=qc*whereNegative(length(x-xc)-r)