ANALYTIC MHD THEORY FOR EARTHS BOW SHOCK AT LOW MACH NUMBERS

A previous MHD theory for the density jump at the Earth's bow shock, which assumed the Alfven (M(A)) and sonic (M(s)) Mach numbers are both >> 1, is reanalyzed and generalized. It is shown that the MHD jump equation can be analytically solved much more directly using perturbation theory, with the ordering determined by M(A) and M(s), and that the first order perturbation solution is identical to the solution found in the earlier theory. The second-order Perturbation solution is calculated, whereas the earlier approach cannot be used to obtain it. The second-order terms generally are important over most of the range of M(A) and M(s) in the solar wind when the angle a between the normal to the bow shock and magnetic field is not close to 0 degrees or 180 degrees (the solutions are symmetric about 90 degrees). This new perturbation solution is generally accurate under most solar wind conditions at 1 AU, with the exception of low Mach numbers when theta is close to 90 degrees. In this exceptional case the new solution does not improve on. the first-order solutions obtained earlier, and the predicted density ratio can vary by 10-20% from the exact numerical MHD solutions. For theta similar to 90 degrees another perturbation solution is derived that predicts the density ratio much more accurately. This second solution is typically accurate for quasi-perpendicular conditions. Taken together, these two analytical solutions are generally accurate for the Earth's bow shock, except in the rare circumstance that M(A) less than or equal to 2. MHD and gasdynamic simulations have produced empirical models ih which the shock's standoff distance alpha(s) is linearly related to the density jump ratio X at the subsolar point. Using an empirical relationship between alpha(s) and X obtained from MHD simulations, alpha(s) values predicted using the MHD solutions for X are compared with the predictions of phenomenological models commonly used for modeling observational data, and with the predictions of a modified phenomenological model proposed recently. The similarities and differences between these results are illustrated using plots bf X and alpha(s) predicted for the Earth's bow shock. The plots. Show that the new analytic solutions agree very well with the exact numerical MHD solutions and that these MHD solutions should replace the the corresponding phenomenological relations in comparisons with data. Furthermore; significant differences exist between the standoff distances Predicted at low M(A) using the MHD models versus those predicted by the new modified phenomenological model. These differences should be amenable to observational testing.