def run(dom, stream):
constraints=VectorConstrainerOverBox()
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.5
sqe.t_end=1.
source=EvaluateExpression()
source.domain=Link(dom)
source.t=Link(sqe)
source.expression=["t","t"]
flow=SteadyIncompressibleFlow()
flow.domain=Link(dom,"domain")
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([sqe,constraints,Simulation([flow],debug=True),ptest],debug=True)
s.writeXML(stream)
s.run()
def run(dom, stream):
constraints=VectorConstrainerOverBox()
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.5
sqe.t_end=1.
source=EvaluateExpression()
source.domain=Link(dom)
source.t=Link(sqe)
source.expression=["t","t"]
flow=SteadyIncompressibleFlow()
flow.domain=Link(dom,"domain")
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([sqe,constraints,Simulation([flow],debug=True),ptest],debug=True)
s.writeXML(stream)
s.run()
except KeyError:
WORKDIR = '.'
dbg = True
dom = RectangularDomain(debug=dbg)
dom.dim = 3
dom.l = [0.5, 1., 1.]
dom.n = [30, 6, 6]
dom.order = 1
dom.integrationOrder = -1
sq = Sequencer(debug=dbg)
sq.t = 0
sq.t_end = 0.8
sq.dt_max = 100.
iob = InterpolateOverBox(debug=dbg)
iob.domain = Link(dom, "domain")
iob.value_left_bottom_front = [-1., 0., 0.]
iob.value_right_bottom_front = [0., 0., 0.]
iob.value_left_bottom_back = [-1., 0., 0.]
iob.value_right_bottom_back = [0., 0., 0.]
iob.value_left_top_front = [-1., 0., 0.]
iob.value_right_top_front = [0., 0., 0.]
iob.value_left_top_back = [-1., 0., 0.]
iob.value_right_top_back = [0., 0., 0.]
# iob.value_left_bottom_front=[-1.,0.]
# iob.value_right_bottom_front=[0.,0.]
# iob.value_left_bottom_back=[-1.,0.]
dbg=True dom=RectangularDomain(debug=dbg) dom.dim=3 dom.l=[0.5,1.,1.] dom.n=[30,6,6] dom.order=1 dom.integrationOrder=-1 sq=Sequencer(debug=dbg) sq.t=0 sq.t_end=0.8 sq.dt_max=100. iob=InterpolateOverBox(debug=dbg) iob.domain=Link(dom,"domain") iob.value_left_bottom_front=[-1.,0.,0.] iob.value_right_bottom_front=[0.,0.,0.] iob.value_left_bottom_back=[-1.,0.,0.] iob.value_right_bottom_back=[0.,0.,0.] iob.value_left_top_front=[-1.,0.,0.] iob.value_right_top_front=[0.,0.,0.] iob.value_left_top_back=[-1.,0.,0.] iob.value_right_top_back=[0.,0.,0.] # iob.value_left_bottom_front=[-1.,0.] # iob.value_right_bottom_front=[0.,0.] # iob.value_left_bottom_back=[-1.,0.]
