A simple model of nonlinear diffusive shock acceleration

We present a simple model of nonlinear diffusive shock acceleration (also called first-order Fermi shock acceleration) that determines the shock modification, spectrum, and efficiency of the process in the plane-wave, steady state approximation as a function of an arbitrary injection parameter, eta. The model, which uses a three-power-law form for the accelerated particle spectrum and contains only simple algebraic equations, includes the essential elements of the full nonlinear model and has been tested against Monte Carlo and numerical kinetic shock models. We include both adiabatic and Alfven wave heating of the upstream precursor. The simplicity and ease of calculation make this model useful for studying the basic properties of nonlinear shock acceleration, as well as providing results accurate enough for many astrophysical applications. It is shown that the shock properties depend upon the shock speed u(0) with respect to a critical value u(0)* proportional to eta p(max)(1/4), which is a function of the injection rate eta and maximum accelerated particle momentum p(max). For u(0) < U-0*, acceleration is efficient and the shock is strongly modified by the back pressure of the energetic particles. In this case, the overall compression ratio is given by r(tot) approximate to 1.3M(S0)(3/4) if M-S0(2) > M-A0, Or by r(tot) approximate to 1.5M(A0)(3/8) in the opposite case (M-S0 is the sonic Mach number and M-A0 is the Alfven Mach number). If u(0) > u(0)*, the shock, although still strong, becomes almost unmodified and accelerated particle production decreases inversely proportional to u(0).