/[escript]/trunk/modellib/test/python/run_convection.py
ViewVC logotype

Contents of /trunk/modellib/test/python/run_convection.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 700 - (show annotations)
Thu Apr 6 00:13:40 2006 UTC (14 years, 3 months ago) by gross
File MIME type: text/x-python
File size: 3057 byte(s)
A few changes in the build mechanism and the file structure so scons can build release tar files:

  * paso/src/Solver has been moved to paso/src 
  * all test_.py are now run_.py files and are assumed to be passing python tests. they can run by 
    scons py_tests and are part of the release test set
  * escript/py_src/test_ are moved to escript/test/python and are installed in to the build directory 
    (rather then the PYTHONPATH).
  * all py files in test/python which don't start with run_ or test_ are now 'local_py_tests'. they are installed i
    by not run automatically.
  * CppUnitTest is now treated as a escript module (against previous decisions).
  * scons realse builds nor tar/zip files with relvant source code (src and tests in seperate files)

the python tests don't pass yet due to path problems.


1 # $Id$
2
3 #
4 #
5 # a very simple convection model in model frame:
6 #
7
8 __copyright__=""" Copyright (c) 2006 by ACcESS MNRF
9 http://www.access.edu.au
10 Primary Business: Queensland, Australia"""
11 __license__="""Licensed under the Open Software License version 3.0
12 http://www.opensource.org/licenses/osl-3.0.php"""
13
14 from esys.escript.modelframe import Link,Simulation
15 from esys.modellib.geometry import RectangularDomain,ScalarConstrainer,VectorConstrainer
16 from esys.modellib.input import Sequencer,InterpolateOverBox,GaussianProfile,LinearCombination
17 from esys.modellib.flow import SteadyIncompressibleFlow
18 from esys.modellib.temperature import TemperatureAdvection
19 from esys.modellib.materials import SimpleEarthModel,GravityForce
20 from esys.modellib.visualization import WriteVTK
21
22 dom=RectangularDomain()
23 dom.order=2
24
25 temp_constraints=ScalarConstrainer()
26 temp_constraints.domain=Link(dom)
27 temp_constraints.top=1
28 temp_constraints.bottom=1
29
30 vel_constraints=VectorConstrainer()
31 vel_constraints.domain=Link(dom)
32 vel_constraints.left=[1,0,0]
33 vel_constraints.right=[1,0,0]
34 vel_constraints.top=[0,1,0]
35 vel_constraints.bottom=[0,1,0]
36 vel_constraints.front=[0,0,1]
37 vel_constraints.back=[0,0,1]
38
39
40 temp_val=InterpolateOverBox()
41 temp_val.domain=Link(dom,"domain")
42 temp_val.right_top_back=Link(dom,"l")
43 temp_val.value_left_bottom_front=1.
44 temp_val.value_right_bottom_front=1.
45 temp_val.value_left_top_front=0.
46 temp_val.value_right_top_front=0.
47 temp_val.value_left_bottom_back=1.
48 temp_val.value_right_bottom_back=1.
49 temp_val.value_left_top_back=0.
50 temp_val.value_right_top_back=0.
51
52 mat=SimpleEarthModel()
53 mat.density0=1.
54 mat.viscocity0=1.
55 mat.rayleigh_number=10000.
56 mat.alpha=0.001
57
58 temp=TemperatureAdvection(debug=True)
59 temp.domain=Link(dom)
60 temp.density=Link(mat,"density0")
61 temp.heat_capacity=Link(mat,"heat_capacity")
62 temp.location_fixed_temperature=Link(temp_constraints,"location_of_constraint")
63 temp.fixed_temperature=Link(temp_val,"out")
64 temp.safety_factor=0.01
65 mat.temperature=Link(temp,"temperature")
66
67
68 grav=GravityForce()
69 grav.domain=Link(dom,"domain")
70 grav.direction=[0.,-1.,0.]
71 grav.density=Link(mat,"density")
72 grav.gravity=Link(mat,"gravity")
73
74
75 vel=SteadyIncompressibleFlow(debug=True)
76 vel.domain=Link(dom)
77 vel.internal_force=Link(grav,"gravity_force")
78 vel.viscosity=Link(mat,"viscosity")
79 vel.location_prescribed_velocity=Link(vel_constraints,"location_of_constraint")
80 vel.rel_tol=1.e-6
81 temp.velocity=Link(vel,"velocity")
82
83 sq=Sequencer()
84 sq.t_end=0.005
85
86 vis=WriteVTK()
87 vis.t=Link(sq)
88 vis.scalar=Link(temp,"temperature")
89 vis.vector=Link(vel,"velocity")
90 vis.stride=5
91
92 per=GaussianProfile()
93 per.domain=Link(dom)
94 per.x_c=[0.5,0.5,0.5]
95 per.A=0.0001
96 per.width=0.01
97 per.r=0
98
99 lc=LinearCombination()
100 lc.f0=1.
101 lc.v0=Link(per,"out")
102 lc.f1=1.
103 lc.v1=Link(temp_val,"out")
104 temp.temperature=Link(lc,"out")
105
106 # s=Simulation([dom,sq,temp_constraints,vel_constraints,temp_val,per,lc,mat,grav,Simulation([vel],debug=True),temp,vis],debug=True)
107 s=Simulation([dom,sq,Simulation([vel],debug=True),temp,vis],debug=True)
108 s.writeXML()
109 s.run()

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26