### saving time-varying Expressions

83

views

0

Hi,

I'm trying to save my exact solution u which was defined as an Expression with varying parameter t; unfortunately, only the first step is saved. Below is the basic code.

I'm trying to save my exact solution u which was defined as an Expression with varying parameter t; unfortunately, only the first step is saved. Below is the basic code.

```
from fenics import *
num_steps = 20
dt = 0.1
# Create mesh and define function space
nx = ny = 30
mesh = RectangleMesh(Point(-2, -2), Point(2, 2), nx, ny)
V = FunctionSpace(mesh, 'P', 1)
import sympy as sym
x, y, t = sym.symbols('x[0], x[1], t')
u_e = x*x*y + t
#f_e = -(sym.diff(sym.diff(u_e,x),x) + sym.diff(sym.diff(u_e,y),y)) + sym.diff(u_e,t)
u_code = sym.printing.ccode(u_e)
#f_code = sym.printing.ccode(f_e)
u_0 = Expression(u_code, degree = 2, t = 0)
t = 0
vtkfile = File('save_f_results/u_exact.pvd')
for n in range(num_steps):
u_0.t = t
u_exact = interpolate(u_0, V)
# Update current time
t += dt
# Save to file and plot solution
vtkfile << (u_exact, float(t))
```

My goal is to be able to save the exact u and exact f as they vary in time. This has troubled me for some time so if anyone can offer suggestions, I appreciate any help.

Community: FEniCS Project

### 1 Answer

Please login to add an answer/comment or follow this question.