In this paper, we propose a multiscale numerical algorithm to simulate the hemolytic release of hemoglobin (Hb) from red blood cells (RBCs) flowing through sieves containing micropores with mean diameters smaller than RBCs. Analyzing the RBC damage in microfiltration is important in the sense that it can quantify the sensitivity of human erythrocytes to mechanical hemolysis while they undergo a high shear rate and high deformation. Here, the numerical simulations are carried out via the lattice Boltzmann method and spring connected network (SN) coupled with an immersed boundary method. To predict the RBC sublytic damage, a sub-cellular damage model derived from molecular dynamics simulations is incorporated in the cellular solver. In the proposed algorithm, the local RBC strain distribution calculated by the cellular solver is used to obtain the pore radius on the RBC membrane. Index of hemolysis (IH) is …