I am trying to solve the problem of 2d transient diffusion equation.
I created block like this figure, 60×6.
All of edges (E1,E2,E3,E4) is Dirichlet Boundary, 0.
where S is the concentration of each species(A,B,C,D), Ds is the diffusion constant for each species, and Rs is the reaction term for each species.
Diffusion constant :
What should I write program?
Can anybody help me?