This paper reports an analytical investigation of magnetohydrodynamic phenomena in electromagnetic levitation processes. The flow is treated as a Stokes flow, and the turbulence in the system is accounted for by using a constant eddy viscosity model, which may be derived from the Prandtl mixing length theory. The curl of the electromagnetic force field for flow calculations and the global lifting force for earthbound levitation are derived. The stream function is introduced, and the resultant scalar equation for fluid flow is solved by the method of separation of variables. The analytical solutions obtained are applicable to all applied frequencies and all conducting materials. Both detailed calculations and asymptotic analyses are presented. Calculated results illustrate that the flow field is characterized by two toroidal recirculating loops, and is strongly correlated to the distribution of the curl of the force field, which in turn depends on the coil placement. Compared with sophisticated numerical computations, the analytical solutions predict well the flow pattern, but underestimate the velocity magnitude as a result of overestimating the eddy sizes by the eddy viscosity model. The asymptotic analyses show that, for a high frequency of applied field or a perfect conductor, the global levitational force is proportional to the square of the input current and increases with increasing conductivity or frequency. The characteristic velocity is proportional to the input current, but inversely proportional to the square root of the liquid density, for a turbulent flow, while it is proportional to the square of the input current, but inversely proportional to the molecular viscosity, for a Stokes flow.