I am reading some research papers on Taylor-Sedov type self-similar solution to the blast wave problems...
It is written that for the existence of similarity solutions Mach number should be constant...please explain this?
and Please explain when it is possible to find self-similar flows with variable Mach numbers?
Also Please tell me range of Mach and Knudsen numbers which are applicable in the interstellar medium ?
Please suggests data or reference....
According to Wikipedia, a similarity solution is a form of scale invariance such that, in the case of fluid flow, the flow "looks" the same independent of time or length scale.
It is written that for the existence of similarity solutions Mach number should be constant...please explain this?
In the case of Taylor-Sedov-von Neumann like solutions [e.g., see pages 192-196 in Whitham, 1999], the relevant parameter is the position of the shock wave at $r = R\left( t \right)$, given by: $$ R\left( t \right) = k \left( \frac{ E }{ \rho_{up} } \right)^{1/5} \ t^{2/5} \tag{1} $$ where $t$ is time from the initial release of energy $E$ from a point source, $\rho_{up}$ is the ambient gas mass density, and $k$ is a dimensionless parameter used for scaling. These solutions are founded upon two assumptions given as follows:
The second assumption implies the strong shock limit, namely that the downstream parameters are given by (in shock frame): $$ \begin{align} U_{dn} & = \left( \frac{ 2 }{ \gamma + 1 } \right) U_{up} \tag{2a} \\ \rho_{dn} & = \left( \frac{ \gamma + 1 }{ \gamma - 1 } \right) \rho_{up} \tag{2b} \\ P_{dn} & = \left( \frac{ 2 }{ \gamma + 1 } \right) \rho_{up} \ U_{up}^{2} \tag{2c} \end{align} $$ where subscript $up$($dn$) corresponds to upstream(downstream) averages, $\gamma$ is the ratio of specific heats, $U_{j}$ is bulk flow speed in region $j$ (i.e., $up$ or $dn$), and $P_{j}$ is the pressure (here just using dynamic or ram pressure) in region $j$.
It can be seen from Equation 1 the only parameter of the ambient air/gas relevant is the density, $\rho_{up}$. Note that the shock speed, $U_{up}$, will then be given by $dR/dt$ (i.e., the time derivative of Equation 1), or: $$ \begin{align} U_{up}\left( t \right) & = \frac{ 2 \ k }{ 5 } \left( \frac{ E }{ \rho_{up} } \right)^{1/5} \ t^{-3/5} \tag{3a} \\ & = \frac{ 2 \ k^{5/2} }{ 5 } \sqrt{ \frac{ E }{ \rho_{up} } } \ R^{-3/2} \tag{3b} \end{align} $$ Since the upstream density and pressure are assumed constant, then the upstream sound speed, $C_{s,up}$, must be constant as well.
Answer 1: The Mach number would be proportionally related to time as $M \propto t^{-3/5} \propto R^{-3/2}$, which is not constant but can be expressed independently of temporal or spatial scales. Over very small radial distances one can, of course, approximate the Mach number as being roughly constant but this is not a good approximation in many cases unless the shock is "old."
and Please explain when it is possible to find self-similar flows with variable Mach numbers?
The self-similarity of the above relations arise because there are no independent length or time scales in the solutions. This is because we can rewrite Equations 2a and 2c solely in terms of $E$ and $R$ such that they do not explicitly depend upon length or time. We can solve Equation 1 for $t$ and use it to invert our pressure and flow speed expressions. These expressions are given as: $$ \begin{align} U_{dn} & = \left( \frac{ 4 \ k }{ 5 \left( \gamma + 1 \right) } \right) \left( \frac{ E }{ \rho_{up} } \right)^{1/5} \ t^{-3/5} \tag{4a} \\ & = \left( \frac{ 4 \ k^{5/2} }{ 5 \left( \gamma + 1 \right) } \right) \sqrt{ \frac{ E }{ \rho_{up} } } \ R^{-3/2} \tag{4b} \\ P_{dn} & = \left( \frac{ 8 \ \rho_{up} \ k^{2} }{ 25 \left( \gamma + 1 \right) } \right) \left( \frac{ E }{ \rho_{up} } \right)^{2/5} \ t^{-6/5} \tag{4c} \\ & = \left( \frac{ 8 \ E \ k^{2} }{ 25 \left( \gamma + 1 \right) } \right) R^{-3} \tag{4d} \end{align} $$
Thus you can see that $U_{dn}$ and $P_{dn}$ are only functions (explicitly) of $E$ and $R$. Since $E/\rho_{up}$ has units of $L^{5}/T^{2}$, we can define a dimensionless parameter $\zeta = E t^{2}/\rho_{up} r^{5}$ that all relevant parameters (e.g., $U_{dn}$) must depend upon in some way, where $r$ is just a radial position. We can further define the dimensionless parameter $\xi = r/R\left( t \right) \propto \zeta^{-1/5}$, upon which all of our characteristics must depend.
Answer 2: The Mach number can be variable (i.e., it depends upon radial position), but one can construct a set of equations whereby it does not explicitly depend upon a length or time scale. I think you may have been confusing constant with independence from time and length scales.
Also Please tell me range of Mach and Knudsen numbers which are applicable in the interstellar medium ?
Wikipedia has a list of relevant parameters for the interstellar medium (ISM). You can use the expressions found in this answer to determine relevant sound speeds. To determine the Knudsen number, $K_{T}$, you can estimate the Coulomb collision mean free path and use the $L_{T} = \lvert d \ln{T_{e}}/dx \rvert^{-1}$ as the characteristic scale length, where $T_{e}$ is the electron temperature (you could use a total temperature too) and $d/dx$ is just a one-dimensional spatial gradient.
Answer 3: In the solar wind $K_{T}$ ~ 0.01-10 [e.g., Bale et al., 2013; Horaites et al., 2015; Landi et al., 2014] and $K_{T}$ ~ 1 for the local interstellar medium or LISM [e.g., Baranov, 2009]. For other ISM regions, the uncertainty in $K_{T}$ will be large since it must be inferred from electromagnetic radiation measurements alone, which require numerous assumptions in order to invert the spectrum into a particle velocity distribution.
As for the Mach number in the ISM, that depends upon with respect to what object or reference frame. If you are looking for the Mach number of, say, the heliosphere with respect to the LISM, there are numerous papers on the heliospheric bow shock (it's actually under debate at the moment whether the flow of the LISM is large enough to produce a bow shock here).
