I am solving a time-dependent fluid mechanics problem. I want to store the velocity vector field, the pressure and the time date. Then, I want to vizualize the results with Paraview.
Here is my attempt:
// (u,v) is the velocity vector
// p is the pressure
// i*dt is the time date
savevtk("results/results"+i+".vtu", Th, [i*dt, u,v,p], dataname="time u v p", order=[0,1,1,1]);