Physiologically based pharmacokinetic models have been used to describe the distribution and elimination of ethanol after intravenous administration. These models have been used to estimate the ethanol infusion profile that is sufficient for achieving a prescribed breath ethanol concentration time course in individuals, providing a useful platform for several pharmacokinetic and pharmacodynamic investigations. Mathematical foundations of these models are examined, including the derivation of an explicit set of governing equations in the form of a system of nonlinear ordinary differential equations. These equations can then be used to formulate and refine parameter identification and control strategies. Finally, a framework in which models related to this model can be constructed and analyzed is described.