Solute diffusion through solid-solid interfaces is paramount to many physical processes. From a modeling point of view, the discontinuities in the energy landscape at a sharp interface represent difficulties in predicting solute diffusion that, to date, have not been solved in a consistent manner across length scales. Using an explicit finite volume method, this work is the first to derive numerical solutions to the diffusion equations at a continuum level while including discrete variations in the energy landscape at a bicrystal interface. An atomic jump equation consistent with atomistic descriptions is derived and scaled up into a compendium of model interfaces: monolayer energy barriers, monolayer interfacial traps, multilayered traps, and heterogeneous interfaces. These can track solute segregation behavior and long-range diffusion effects. We perform simulations with data for hydrogen diffusion in structural metals, of relevance to the assessment of the hydrogen embrittlement phenomenon, and point defects in electronic devices. The approach developed represents an advancement in the mathematical treatment of solute diffusion through solid-solid interfaces and an important bridge between the atomistic and macroscopic modeling of diffusion, with potential applications in a variety of fields in materials science and physics.