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

Annotation of /trunk/modellib/test/python/run_flow.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 938 - (hide annotations)
Thu Jan 25 03:40:54 2007 UTC (16 years ago) by gross
File MIME type: text/x-python
File size: 1484 byte(s)
the FinleyReade accepts now gmsh files (use format="gmsh")
and
Simulations are accepting Models only. Moreover, there is a test now
if a all Models targeted by a model in the simulation or a subsimulation
are included in  simulation or a subsimulation. 


1 jgs 147 # $Id$
2    
3     #
4     # small test problem fro temperture advection:
5     #
6     # p=(x0+x1)*t
7     #
8 elspeth 628 __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 gross 706 import os
15 jgs 150 from esys.escript.modelframe import Link,Simulation
16 gross 911 from esys.modellib.geometry import RectangularDomain,VectorConstrainerOverBox
17 jgs 150 from esys.modellib.input import Sequencer
18     from esys.modellib.probe import Probe,EvaluateExpression
19     from esys.modellib.flow import SteadyIncompressibleFlow
20 jgs 147
21     dom=RectangularDomain()
22     dom.order=2
23    
24 gross 911 constraints=VectorConstrainerOverBox()
25 jgs 147 constraints.domain=Link(dom)
26     constraints.left=[1,0,0]
27     constraints.right=[1,0,0]
28     constraints.top=[0,1,0]
29     constraints.bottom=[0,1,0]
30     constraints.front=[0,0,1]
31     constraints.back=[0,0,1]
32    
33     sqe=Sequencer()
34 gross 820 sqe.dt_max=0.5
35 jgs 147 sqe.t_end=1.
36    
37     source=EvaluateExpression()
38     source.domain=Link(dom)
39     source.t=Link(sqe)
40     source.expression=["t","t"]
41    
42     flow=SteadyIncompressibleFlow()
43 gross 911 flow.domain=Link(dom,"domain")
44 jgs 147 flow.internal_force=Link(source,"out")
45     flow.location_prescribed_velocity=Link(constraints,"location_of_constraint")
46     flow.prescribed_velocity=[0.,0.]
47    
48     ptest=Probe()
49     ptest.expression="(x[0]+x[1]-1.)*t"
50     ptest.t=Link(sqe)
51     ptest.value=Link(flow,"pressure")
52    
53 gross 938 s=Simulation([sqe,constraints,Simulation([flow],debug=True),ptest],debug=True)
54 jgs 147 s.writeXML()
55     s.run()

Properties

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

  ViewVC Help
Powered by ViewVC 1.1.26