How can I solve 2D transient diffusion equation with nonlinear source term?

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.

Initial condition:

Diffusion constant :

What should I write program?

Can anybody help me?

Answer 1

[Ravi Kumar]

I don't see any interaction between species being defined by the stated PDEs. You can define coefficients usng PDEModel workflow: d = 1, c = Ds, and f = Rs.

Regards,

Ravi