dolfin team mailing list archive

[Bento <question153689@xxxxxxxxxxxxxxxxxxxxx>]

New question #153689 on DOLFIN:
https://answers.launchpad.net/dolfin/+question/153689

Hi,

I try to solve a time-dependent problem, however the plotting does not show the change of the function over time. The plot window always shows the function in that timestep in which the plot command was first called. What am I doing wrong? The time-dependent example look similar, but they work.
Here is the simplest example that goes wrong:

from dolfin import *

T_initial = Constant(293.0)

dt = 1.0
t = dt

mesh = Interval(200, 0, 20.0)

V = FunctionSpace(mesh, "CG", 1)

T = TrialFunction(V)
v = TestFunction(V)

def inflow(x, on_boundary):
    return on_boundary and x[0] < 1.0
    
BC = DirichletBC(V, 800.0, inflow)
T_prev = interpolate(T_initial, V)

a = (T*v + dt*T.dx(0)*v.dx(0))*dx
L = T_prev*v*dx

f = Function(V)


while t < 1000:
    A = assemble(a)
    l = assemble(L)
    BC.apply(A,l)
    solve(A,f.vector(), l)
    plot(f)
    T_prev.assign(f)
    t += dt

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.