I'm trying to solve the heat equation in FiPy and right now it works, but I have one problem:
The material I have to simulate has different diffusion coefficients in the x and y direction (due to its fibre layerd structure). I know, that you can make your diffusion coefficient dependent to a variable in FiPy.
My grid is rectengular, but in the future is has to be able to be unstructured. This would mean, the diffusion coefficient is different for every face.
Has anybody a clue on how to implement this type of dependency?