/[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 819 - (show annotations)
Sun Aug 27 23:07:51 2006 UTC (14 years, 9 months ago) by gross
File MIME type: text/x-python
File size: 3014 byte(s)
new version of FinleyReader etc. They can noe update node coordinates but have to be part of the Simulation
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 import os
15 from esys.escript.modelframe import Link,Simulation
16 from esys.modellib.geometry import RectangularDomain,ScalarConstrainer,VectorConstrainer
17 from esys.modellib.input import Sequencer,InterpolateOverBox,GaussianProfile,LinearCombination
18 from esys.modellib.flow import SteadyIncompressibleFlow
19 from esys.modellib.temperature import TemperatureAdvection
20 from esys.modellib.materials import SimpleEarthModel,GravityForce
21 from esys.modellib.visualization import WriteVTK
22
23 try:
24 WORKDIR=os.environ['MODELLIB_WORKDIR']
25 except KeyError:
26 WORKDIR='.'
27
28
29 dom=RectangularDomain()
30 dom.order=2
31
32 temp_constraints=ScalarConstrainer()
33 temp_constraints.domain=Link(dom)
34 temp_constraints.top=1
35 temp_constraints.bottom=1
36
37 vel_constraints=VectorConstrainer()
38 vel_constraints.domain=Link(dom)
39 vel_constraints.left=[1,0,0]
40 vel_constraints.right=[1,0,0]
41 vel_constraints.top=[0,1,0]
42 vel_constraints.bottom=[0,1,0]
43 vel_constraints.front=[0,0,1]
44 vel_constraints.back=[0,0,1]
45
46
47 temp_val=InterpolateOverBox()
48 temp_val.domain=Link(dom,"domain")
49 temp_val.value_left_bottom_front=1.
50 temp_val.value_right_bottom_front=1.
51 temp_val.value_left_top_front=0.
52 temp_val.value_right_top_front=0.
53 temp_val.value_left_bottom_back=1.
54 temp_val.value_right_bottom_back=1.
55 temp_val.value_left_top_back=0.
56 temp_val.value_right_top_back=0.
57
58 mat=SimpleEarthModel()
59 mat.density0=1.
60 mat.viscocity0=1.
61 mat.rayleigh_number=10000.
62 mat.alpha=0.001
63
64 temp=TemperatureAdvection(debug=True)
65 temp.domain=Link(dom)
66 temp.density=Link(mat,"density0")
67 temp.heat_capacity=Link(mat,"heat_capacity")
68 temp.location_fixed_temperature=Link(temp_constraints,"location_of_constraint")
69 temp.fixed_temperature=Link(temp_val,"out")
70 temp.safety_factor=0.01
71 mat.temperature=Link(temp,"temperature")
72
73
74 grav=GravityForce()
75 grav.domain=Link(dom,"domain")
76 grav.direction=[0.,-1.,0.]
77 grav.density=Link(mat,"density")
78 grav.gravity=Link(mat,"gravity")
79
80
81 vel=SteadyIncompressibleFlow(debug=True)
82 vel.domain=Link(dom)
83 vel.internal_force=Link(grav,"gravity_force")
84 vel.viscosity=Link(mat,"viscosity")
85 vel.location_prescribed_velocity=Link(vel_constraints,"location_of_constraint")
86 vel.rel_tol=1.e-6
87 temp.velocity=Link(vel,"velocity")
88
89 sq=Sequencer()
90 sq.t_end=0.005
91
92 vis=WriteVTK()
93 vis.t=Link(sq)
94 vis.scalar=Link(temp,"temperature")
95 vis.vector=Link(vel,"velocity")
96 vis.dt=0.0005
97 vis.filename=WORKDIR+"/temp.xml"
98
99 per=GaussianProfile()
100 per.domain=Link(dom)
101 per.x_c=[0.5,0.5,0.5]
102 per.A=0.0001
103 per.width=0.01
104 per.r=0
105
106 lc=LinearCombination()
107 lc.f0=1.
108 lc.v0=Link(per,"out")
109 lc.f1=1.
110 lc.v1=Link(temp_val,"out")
111 temp.temperature=Link(lc,"out")
112
113 s=Simulation([sq,dom,Simulation([vel],debug=True),temp,vis],debug=True)
114 s.writeXML()
115 s.run()

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26