Kelvin-Helmholtz (K-H) instability in a stably stratified background is studied numerically in order to support the AirBorne Laser (ABL) project. The wind shear associated with K-H events produces irregularities in the air density which result in both refraction and scattering of laser beams used for communications and for anti-missile purposes. Optical distortion is measured from the simulations through analysis of the second order structure function of the potential temperature fluctuations. To date, direct numerical simulations (DNS) have been conducted at bulk Richardson numbers of 0.05-0.1 and Reynolds numbers of 2000-2500 using up to 1000x350x2000 Fourier modes, run on as many as 1000 processors of the Cray T3E at ERDC. The structure function coefficients, power law exponents, and turbulence inner scale from these simulations are in good agreement with both theoretical predictions and field measurements. An initial phase of a companion large-eddy simulation (LES) study has been completed where the existing DNS has been replicated at lower resolution using dynamic subgrid scale modeling ideas. The results of these simulations show that the first and second order statistics (including the structure function) can be reproduced with acceptable accurately using meshes that are coarser by factors of up to six in each spatial direction. When the effect of the increased time step is taken into account, the LES can be completed at computational cost that is 1/6⁴ = 1/1296 of that required for the DNS. We also anticipate that it should be possible to produce high-fidelity simulation results at Reynolds numbers that are one to two orders of magnitude higher using LES, which places the simulations within the range of high-altitude wind shear events.