CODEWITHSUNDEEP
pdefipy
Baudelaire
Baudelaire

Reputation: 1

fit pde solution to data

I have solved a pde (second Fick's law with non-uniform diffusivity) using fipy.

I now would like to estimate the diffusion coefficient in one of the layers by fitting solutions of the pde to the measurements.

Do you have any suggestions on how to proceed?

Many thanks in advance!

Upvotes: 0

Views: 228

Answers (1)

jeguyer
jeguyer

Reputation: 2484

Interesting problem. As @wd15 suggests, scipy.optimize has useful routines.

I put together a notebook that combines the example code from scipy.optimize.curve_fit with our steady-state diffusion example with non-uniform-diffusivity in examples.diffusion.mesh1D.

Your GitHub issue may be better for any extended discussion on this problem.

Upvotes: 1

Related Questions