Numerical solution to hyperbolic stiff simultaneous PDEs

To develop a code to solve hyperbolic stiff simultaneous PDEs numerically (not analytically) mimicking the reaction between a fluid and solid. The problem contains four equations with four dependent variables (concentration of fluid and solid, temperature of fluid and solid) and two independent variables (time and distance). Please read the attached document.