Mathematical epidemiology holds prime importance for comprehending the dynamics of infectious diseases. Consequently, mathematical model of hepatitis B with fractional-order derivative under Caputo sense is primarily focused in this research. The analysis of the required solution is qualitatively derived by applying the fixed-point theory approach. By perturbing the proposed model, the Ulam–Hyer’s stability techniques are further derived. To achieve the iterative series solution of the proposed system of hepatitis, the modified Euler method like Taylor’s series method is utilized. For validation and importance of the fractional operators, sufficient significant numerical results at various fractional orders are presented and compared them with the integer order. It is inferred from this research that, by using the fractional-order method, the transmission mechanism of hepatitis B disease can be acutely revealed. This study may provide positive theoretical support for the prevention and treatment of hepatitis B disease.