CFD Online Discussion Forums - Smooth reconstruction of flowsolution

there are variaty of methods for this purpose,

your field is formally C0 space, because in formal numerical method we just attend to continuty of field not its derivative (it is very obvious in FEM)

a very general method (not computationally cheap) is to assume a desired smooth space (C^n, n>=1) with its basis (e.g. high order spline basis) and then solve a least square problem to find unknown interpolation coefficient, this system is not essentially wel-posed (and could be over/under specify) and maybe methods like SVD are needed.

if you search computer graphics literature you find several alternatives, your problem is common there.

if you look for a toolbox have a look at this: http://www.farfieldtechnology.com/

as a very simple and easy to implement solution:

assume your filed is member of L2 space (space filled by square integrable function, Hilbert space, L2-norm = \int (.)^2 ), so there is not garantee to have continuous derivative.

then you want to have a filedwith continuty of first derivative, so we should map from L2 to H1 (H1 is Sobolev space filled with function has continuous first derivative, H1-norm = \int (.)^ + \nabla(.)^2)

assume g is original and G is smoothed filed, for this mapping you should solve this poisson equation (weak solution):

- \nabla \cdot (\nabla G) + G = g

to control lenght scale of smoothing, it could be modified to:

- \nabla \cdot (D \nabla G) + G = g

where D is diffusion coefficent control smoothing, D = 0, could give original field.

Hope this helps.