Error estimation using a function of your solution vector

[Bruno Blais]

Hello everyone,

I was wondering what was the simplest way to carry out error estimation using a function of your solution vector. For example, to carry out mesh adaptation using  a Kelly Error estimator but using the vorticity as the input variable instead of the velocity vector.

Is there a way to do it in a similar fashion as the way we can post-process with DataPostprocessorScalar<dim> or DataPostprocessorVector<dim>? (which is what we use for the vorticity and the q-criterion post-processing)

I have tried looking through the documentation of the Kelly error estimator or the steps, but I have not found a similar example.

If such postprocessor cannot be used, what would you suggest to be the best way?
I am afraid my knowledge on error estimator is really limited.

Best regards and thank you for your time.

Bruno

[Daniel Arndt]

Bruno

Hello everyone,

I was wondering what was the simplest way to carry out error estimation using a function of your solution vector. For example, to carry out mesh adaptation using  a Kelly Error estimator but using the vorticity as the input variable instead of the velocity vector.

Is there a way to do it in a similar fashion as the way we can post-process with DataPostprocessorScalar<dim> or DataPostprocessorVector<dim>? (which is what we use for the vorticity and the q-criterion post-processing)

DataPostprocessor computes values in certain evaluation points that are not (necessarily) related to any given finite element space.

On the other hand, KellyerrorEstimator expects a finite element vector. What you need to do is to represent the vorticity in some finite element space.

You would typically either interpolate or project with respect to this new space.

Best,

Daniel