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

Revision 819 - (show annotations)
Sun Aug 27 23:07:51 2006 UTC (15 years, 11 months ago) by gross
File MIME type: text/x-python
File size: 1453 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 # 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,dom,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

 ViewVC Help Powered by ViewVC 1.1.26