Implement a caching strategy for the Turbine friction fields.
The interpolation of the turbine expressions onto the discrete function spaces
takes a lot of time so it makes sense to cache it whenever possible