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

Revision 628 - (show annotations)
Thu Mar 23 02:27:57 2006 UTC (15 years, 2 months ago) by elspeth
File MIME type: text/x-python
File size: 1443 byte(s)
```Copyright information added.

```
 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 from esys.escript.modelframe import Link,Simulation 15 from esys.modellib.geometry import RectangularDomain,VectorConstrainer 16 from esys.modellib.input import Sequencer 17 from esys.modellib.probe import Probe,EvaluateExpression 18 from esys.modellib.flow import SteadyIncompressibleFlow 19 20 dom=RectangularDomain() 21 dom.order=2 22 23 constraints=VectorConstrainer() 24 constraints.domain=Link(dom) 25 constraints.left=[1,0,0] 26 constraints.right=[1,0,0] 27 constraints.top=[0,1,0] 28 constraints.bottom=[0,1,0] 29 constraints.front=[0,0,1] 30 constraints.back=[0,0,1] 31 32 sqe=Sequencer() 33 sqe.dt_max=0.3 34 sqe.t_end=1. 35 36 source=EvaluateExpression() 37 source.domain=Link(dom) 38 source.t=Link(sqe) 39 source.expression=["t","t"] 40 41 flow=SteadyIncompressibleFlow() 42 flow.domain=Link(dom) 43 flow.internal_force=Link(source,"out") 44 flow.location_prescribed_velocity=Link(constraints,"location_of_constraint") 45 flow.prescribed_velocity=[0.,0.] 46 47 ptest=Probe() 48 ptest.expression="(x[0]+x[1]-1.)*t" 49 ptest.t=Link(sqe) 50 ptest.value=Link(flow,"pressure") 51 52 s=Simulation([dom,sqe,Simulation([flow],debug=True),ptest],debug=True) 53 s.writeXML() 54 s.run()

Properties

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