/[escript]/trunk/doc/examples/diffusion.py
ViewVC logotype

Annotation of /trunk/doc/examples/diffusion.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 102 - (hide annotations)
Wed Dec 15 07:08:39 2004 UTC (15 years, 11 months ago) by jgs
Original Path: trunk/esys2/doc/user/examples/diffusion.py
File MIME type: text/x-python
File size: 720 byte(s)
*** empty log message ***

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/h*T,eta=eta,g=eta*Tref)
33     T=mypde.getSolution()
34     T.saveDX("T%d.dx"%i)

Properties

Name Value
svn:eol-style native
svn:keywords Author Date Id Revision

  ViewVC Help
Powered by ViewVC 1.1.26