The purpose of this thesis is to accurately simulate in 3D the HH901 jet in the Mystic Mountain Formation of the Carina Nebula. Astronomers present a narrow-band Wide Field Camera image of Carina and the morphology of some astrophysical jets, including HH901. The simulation attempts to replicate features of the jet, among which are pulses, bow shock, terminal Mach disk, and Kelvin-Helmholtz rollup. We use the gas dynamical equations to solve for density, velocity, and temperature. The numerical methods used to solve the equations are third-order WENO (weighted essentially non-oscillatory) and third-order Runge-Kutta. Graphs of density and radiative cooling demonstrate the effect of adding wind (nonzero ambient velocity). The paper discusses the altering of the ambient velocity and final time to fit the shape of the jet in the Hubble image. The suggested next steps are simulating the other HH901 jet and comparing the jets’ atomic makeups to advance understanding of astrophysical jets.