Heat transmission in systems or thermal devices relies on the thermal conductivity of nanofluids. Thermal conductivity (32%) of the hybrid nanofluid (i.e. ethelene glycol conveying zinc and titania nanoparticles) is higher than that (13%) of the ordinary nanofluid (i.e. ethelene glycol conveying zinc nanoparticles), according to experimental data. When Lorentz force, thermal radiation, viscous dissipation, heat flux causes to concentration gradient as well as mass flux causes to thermal gradient are significant, nothing is known on the flow of chemically reactive Ethelene glycol conveying zinc and titania nanoparticles due to heat flux due to thermal inertia through porous medium on sphere. Keller box scheme with the combination of finite difference procedure is opted to solve the governing equations by transforming them into a nonlinear and non-dimensional system of partial differential equations. Outcomes (such as the variation in the temperature profile due to Dufour number) are presented in the form of bar plots and plots in two cases i.e., 2 EG+ ZnO+TiO and EG+ ZnO. An in-built function, ‘nntool’ (ANN) in MATLAB is used to validate the data of the engineering parameters including shear stress. An important observation is that, with respect to magnetic field, the shear stress drops to -0.146 (HNF-Hybrid Nano Fluid) and -0.08887 (NF-Nano Fluid), whereas the skin friction factor increases by 2.764055 (HNF), 2.028481 (NF) with an increment in volume fraction of nanoparticle. The reported decreases in heat transfer rate for Dufour number are - 0.15607 (HNF) and - 0.15815 (NF). Also it was revealed that when a growth in Soret number causes to lessen the mass transfer rate and that the decrement amounts for HNF and NF in relation to the Soret number are - 0.03187 and - 0.03904. And also, it is noticed that the greater the thermal relaxation parameter, the lower the temperature of the fluid.