Computational fluid dynamics simulations using the WENO-LF method are applied to high Mach number nonrelativistic astrophysical jets, including the effects of radiative cooling. Our numerical methods have allowed us to simulate astrophysical jets at much higher Mach numbers than have been attained (Mach 20) in the literature. Our simulations of the HH 1-2 astrophysical jets are at Mach 80. Simulations at high Mach numbers and with radiative cooling are essential for achieving detailed agreement with the astrophysical images.