user/examples/diffusion.py

# $Id$
from mytools import *
from esys.escript import *
import esys.finley
#... set some parameters ...
x_c=[0.02,0.002]
r=0.001
q0=50.e6
Tref=0.
rhocp=2.6e6
eta=75.
kappa=240.
t_end=5.
# ...time step size and counter ...
h=0.1
i=0
t=0
#... generate domain ...
mydomain = esys.finley.Rectangle(l0=0.05,l1=0.01,n0=250, n1=50)
#... open PDE ...
mypde=Helmholtz(mydomain)
# ... set heat source: ....
x=mydomain.getX()
q=q0*(length(x-x_c)-r).whereNegative()
# ... set initial temperature ....
T=Tref
# ... start iteration:
while t<t_end:
      i+=1
      t+=h
      print "time step :",t
      mypde.setValue(kappa=kappa,omega=rhocp/h,f=q+rhocp/h*T,eta=eta,g=eta*Tref)
      T=mypde.getSolution()
      T.saveDX("T%d.dx"%i)
1	# $Id$
2	from mytools import *
3	from esys.escript import *
4	import esys.finley
5	#... set some parameters ...
6	x_c=[0.02,0.002]
7	r=0.001
8	q0=50.e6
9	Tref=0.
10	rhocp=2.6e6
11	eta=75.
12	kappa=240.
13	t_end=5.
14	# ...time step size and counter ...
15	h=0.1
16	i=0
17	t=0
18	#... generate domain ...
19	mydomain = esys.finley.Rectangle(l0=0.05,l1=0.01,n0=250, n1=50)
20	#... open PDE ...
21	mypde=Helmholtz(mydomain)
22	# ... set heat source: ....
23	x=mydomain.getX()
24	q=q0*(length(x-x_c)-r).whereNegative()
25	# ... set initial temperature ....
26	T=Tref
27	# ... start iteration:
28	while t<t_end:
29	i+=1
30	t+=h
31	print "time step :",t
32	mypde.setValue(kappa=kappa,omega=rhocp/h,f=q+rhocp/hT,eta=eta,g=etaTref)
33	T=mypde.getSolution()
34	T.saveDX("T%d.dx"%i)
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision