Abstract Numerical approaches have long been used to examine material behaviors that are governed by diffusion. The surface water and groundwater become polluted with regard to time. The pollutants travel from place to place in extended time intervals interacting with the solid soil particles. The diffusion of solute with solid soil particles idealized as spheres are studied in this dissertation. The objective was to numerically simulate diffusion through a sphere. A steady-state finite difference method was used to study the diffusion in a sphere. It was assumed that a bulk liquid solution interacts with solid spheres through a diffusion limited process which is equivalent to a volume of limited bath boundary condition. A sphere of radius 0.07 cm is placed in the beaker containing the solution. The radius of the sphere is divided into 15 nodes each of radius 0.005 cm. Initially, the sphere is free from the solute. As time progresses the diffusion in the sphere takes place. The solute from outside the sphere gets infused in the sphere. By using a numerical technique, the concentration at different nodes at different time intervals is studied. The process is carried out for two different solutions with different concentration values.