# Contents of /trunk/modellib/test/python/run_flow.py

Revision 728 - (show annotations)
Mon May 15 02:01:18 2006 UTC (15 years ago) by gross
File MIME type: text/x-python
File size: 1449 byte(s)
```modelframe test examples work now
```
 1 # \$Id\$ 2 3 # 4 # small test problem fro temperture advection: 5 # 6 # p=(x0+x1)*t 7 # 8 __copyright__=""" Copyright (c) 2006 by ACcESS MNRF 9 10 Primary Business: Queensland, Australia""" 11 __license__="""Licensed under the Open Software License version 3.0 12 13 14 import os 15 from esys.escript.modelframe import Link,Simulation 16 from esys.modellib.geometry import RectangularDomain,VectorConstrainer 17 from esys.modellib.input import Sequencer 18 from esys.modellib.probe import Probe,EvaluateExpression 19 from esys.modellib.flow import SteadyIncompressibleFlow 20 21 dom=RectangularDomain() 22 dom.order=2 23 24 constraints=VectorConstrainer() 25 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 sqe.dt_max=0.3 35 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 flow.domain=Link(dom) 44 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 s=Simulation([sqe,Simulation([flow],debug=True),ptest],debug=True) 54 s.writeXML() 55 s.run()

## Properties

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