The interactions between boundary layer turbulence, including Langmuir turbulence, and submesoscale processes in the oceanic mixed layer are described using large-eddy simulations of the spindown of a temperature front in the presence of submesoscale eddies, winds, and waves. The simulations solve the surface-wave-averaged Boussinesq equations with Stokes drift wave forcing at a resolution that is sufficiently fine to capture small-scale Langmuir turbulence. A simulation without Stokes drift forcing is also performed for comparison. Spatial and spectral properties of temperature, velocity, and vorticity fields are described, and these fields are scale decomposed in order to examine multiscale fluxes of momentum and buoyancy. Buoyancy flux results indicate that Langmuir turbulence counters the restratifying effects of submesoscale eddies, leading to small-scale vertical transport and mixing that is 4 times greater than in the simulations without Stokes drift forcing. The observed fluxes are also shown to be in good agreement with results from an asymptotic analysis of the surface-wave-averaged, or Craik–Leibovich, equations. Regions of potential instability in the flow are identified using Richardson and Rossby numbers, and it is found that mixed gravitational/symmetric instabilities are nearly twice as prevalent when Langmuir turbulence is present, in contrast to simulations without Stokes drift forcing, which are dominated by symmetric instabilities. Mixed layer depth calculations based on potential vorticity and temperature show that the mixed layer is up to 2 times deeper in the presence of Langmuir turbulence. Differences between measures of the mixed layer depth based on potential vorticity and temperature are smaller in the simulations with Stokes drift forcing, indicating a reduced incidence of symmetric instabilities in the presence of Langmuir turbulence.