search

KikeM / romtime

Unsteady Finite Element Reduced Order Model for Time Moving Domains
1 stars 1 forks source link

BUG : Wrong interpolation by fenics for probes #34

Closed KikeM closed 3 years ago

KikeM commented 3 years ago

I am implementing a linear interpolation between nodal values to compare FOM and ROM solutions.

Since I had already implemented a probe interpolation using the fenics __call__ method for CoefficientFunctions, I expected my interpolation to match these values.

To my surprise, they don't. In fact, they do at the outflow (which is reassuring although trivial. Reassuring because the mass conservation calculation uses this value; trivial because it is a nodal value, so no interpolation takes place.), but they don't at the mid-point x=0.5.

I think there might be a problem with the moving mesh and the interpolation routine used inside fenics. I am going to stick to my interpolation and use it to compare FOM and ROM solutions.

legend:

half_0

outflow

KikeM commented 3 years ago

Closed with https://github.com/KikeM/romtime/pull/35