test/python/flow_test.py

# $Id$

#
#   small test problem fro temperture advection:
#
#   p=(x0+x1)*t
#
__copyright__="""  Copyright (c) 2006 by ACcESS MNRF
                    http://www.access.edu.au
                Primary Business: Queensland, Australia"""
__license__="""Licensed under the Open Software License version 3.0
             http://www.opensource.org/licenses/osl-3.0.php"""

from esys.escript.modelframe import Link,Simulation
from esys.modellib.geometry import RectangularDomain,VectorConstrainer
from esys.modellib.input import Sequencer
from esys.modellib.probe import Probe,EvaluateExpression
from esys.modellib.flow import SteadyIncompressibleFlow

dom=RectangularDomain()
dom.order=2

constraints=VectorConstrainer()
constraints.domain=Link(dom)
constraints.left=[1,0,0]
constraints.right=[1,0,0]
constraints.top=[0,1,0]
constraints.bottom=[0,1,0]
constraints.front=[0,0,1]
constraints.back=[0,0,1]

sqe=Sequencer()
sqe.dt_max=0.3
sqe.t_end=1.

source=EvaluateExpression()
source.domain=Link(dom)
source.t=Link(sqe)
source.expression=["t","t"]

flow=SteadyIncompressibleFlow()
flow.domain=Link(dom)
flow.internal_force=Link(source,"out")
flow.location_prescribed_velocity=Link(constraints,"location_of_constraint")
flow.prescribed_velocity=[0.,0.]

ptest=Probe()
ptest.expression="(x[0]+x[1]-1.)*t"
ptest.t=Link(sqe)
ptest.value=Link(flow,"pressure")

s=Simulation([dom,sqe,Simulation([flow],debug=True),ptest],debug=True)
s.writeXML()
s.run()
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	jgs	150	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	jgs	147
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	jgs	148	s=Simulation([dom,sqe,Simulation([flow],debug=True),ptest],debug=True)
53	jgs	147	s.writeXML()
54			s.run()
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision