Pharyngeal bolus flow was simulated numerically using a finite element method. The bolus liquids were X-ray medium, glucose, and thickener solutions. For a lowviscosity bolus, the simulation showed a reasonable agreement of bolus velocity with X-ray measurements. The influence of bolus density on swallowing velocity was investigated numerically. Although a higher density resulted in a higher bolus velocity, the increase in velocity was modest. When the bolus viscosity was high, it was necessary to apply the slip boundary condition to obtain an agreement for bolus velocity between the simulation and X-ray measurements. The simulations also showed that the method of characteristic shear rate proposed by Zhu et al., Journal of Texture Studies, 2014, 45, 430-439 is effective for predicting the bolus velocity for shear-thinning fluids. In order to discuss the effect of saliva lubrication and the physical meaning of the characteristic shear rate, an immiscible two-layer flow of the core and wall layer was analyzed theoretically by analogy with mesopharyngeal bolus flow. The characteristic shear rate enabled us to correlate the macroscopic flow