  Radiative turbulent mixing layers should be ubiquitous in multi-phase gas
with shear flow. They are a potentially attractive explanation for the high
ions such as OVI seen in high velocity clouds and the circumgalactic medium
(CGM) of galaxies. We perform 3D MHD simulations with non-equilibrium (NEI) and
photoionization modeling, with an eye towards testing simple analytic models.
Even purely hydrodynamic collisional ionization equilibrium (CIE) calculations
have column densities much lower than observations. Characteristic inflow and
turbulent velocities are much less than the shear velocity, and the layer width
$h \propto t_\mathrm{cool}^{1/2}$ rather than $h \propto t_\mathrm{cool}$.
Column densities are not independent of density or metallicity as analytic
scalings predict, and show surprisingly weak dependence on shear velocity and
density contrast. Radiative cooling, rather than Kelvin-Helmholtz instability,
appears paramount in determining the saturated state. Low pressure due to fast
cooling both seeds turbulence and sets the entrainment rate of hot gas, whose
enthalpy flux, along with turbulent dissipation, energizes the layer.
Regardless of initial geometry, magnetic fields are amplified and stabilize the
mixing layer via magnetic tension, producing almost laminar flow and depressing
column densities. NEI effects can boost column densities by factors of a few.
Suppression of cooling by NEI or photoionization can in principle also increase
OVI column densities, but in practice is unimportant for CGM conditions. To
explain observations, sightlines must pierce hundreds or thousands of mixing
layers, which may be plausible if the CGM exists as a `fog' of tiny cloudlets.
