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	jgs	102	# $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