A computational study is performed to develop a capability of analyzing the effect damage has on the aerodynamic performance of airfoils. A Cartesian immersed boundary method is implemented in 2D simulations of the compressible Navier-Stokes equations. These equations are discretized using a Weighted Essentially Non-Oscillatory (WENO) scheme for spatial derivatives and a 4th order Runge-Kutta scheme for temporal derivatives. Results from time-accurate, parallel computations are presented for a NACA 0009 airfoil in both undamaged and damaged states. Time histories of lift, drag, and moment coefficients are shown, along with Mach contours and turbulent kinetic energy contours at three times throughout the simulations. Based on these initial results, it was determined that the boundary conditions used in the simulations were causing reflections of propagating disturbances which contaminate the solution. Further studies must be completed with the current methodology using alternate boundary conditions to validate its results prior to continuing its development for this application in aircraft survivability.