How to output an animation of the flow in paraview

[LS-Maxwell]

I want to output an animation of the flow. For example, when simulating Stokes flow, openpnm can only view the final flow result, but I want to see an animation from t=0 to steady state in paraview. I wonder if you can give me some suggestions as following code for example.

import openpnm as op
import matplotlib.pyplot as plt
import numpy as np
import math
import scipy.stats as spst
import datetime
import time

pn = op.network.Cubic(shape=[15, 15, 15], spacing=1e-6)

#------------------------------------Geometry----------------------------------#
geo = op.models.collections.geometry.spheres_and_cylinders
geo.keys()
geo['pore.diameter']
pn.add_model_collection(geo)
pn.regenerate_models()
#------------------------------------Geometry----------------------------------#

#------------------------------------Phase----------------------------------#
phase = op.phase.Phase(network=pn)
phase['pore.viscosity']=1.0
phase.add_model_collection(op.models.collections.physics.basic)
phase.regenerate_models()
#------------------------------------Phase----------------------------------#


#------------------------------------Stokes flow----------------------------------#
inlet = pn.pores('left')
outlet = pn.pores('right')
flow = op.algorithms.StokesFlow(network=pn, phase=phase)
flow.set_value_BC(pores=inlet, values=1)
flow.set_value_BC(pores=outlet, values=0)
flow.run()
phase.update(flow.soln)
#------------------------------------Stokes flow----------------------------------#

path_to_file = 'E:\code\openpnm\OpenPNM3'
op.io._vtk.project_to_vtk(pn.project, filename=path_to_file+'/Paraview_net')

[ma-sadeghi]

StokesFlow is a steady state algorithm, it's not transient.