An analytical bulk mobility model for hydrodynamic transport equations is developed from a microscopic level and designed for silicon device simulation. Applying Kohler's variational method extended to the regime of nonlinear transport yields the general expression for the mobility as function of carrier temperature, lattice temperature, and doping. Assuming a nonparabolic, isotropic band model and a heated Maxwellian allows for the analytical calculation of the collision integrals. A nonelastic approximation for intravalley acoustic-phonon scattering is proposed, which improves the model in the low-temperature range. Intervalley scattering is treated in a one-mode, equipartition approximation. Here, an accurate analytical approach for all carrier temperatures is derived. For impurity scattering the Brooks-Herring theory is used including Fermi statistics and the effect of dispersive screening. The influence of other effects like anisotropic valleys and perturbation of the density of states by heavy doping are discussed quantitatively. Despite the oversimplified band structure, all essential features of the measured mobility in silicon can be reproduced except in the heavy doping range. The adjusted deformation potentials coincide with estimated sums of the corresponding sets used in full-band Monte Carlo simulation. The method has the potential of an extension to the SiϪSiO 2 system.