The stability of the flow under the magnetic force is one of the classical problems in fluid mechanics. In this paper, the flow in a rectangular duct with different Hartmann (Ha) number is simulated. The finite volume method and the SIMPLE algorithm are used to solve a system of equations and the energy gradient theory is then used to study the (associated) stability of magnetohydrodynamics (MHD). The flow stability of MHD flow for different Hartmann (Ha) number, from Ha=1 to 40, at the fixed Reynolds number, Re=190 are investigated. The simulation is validated firstly against the simulation in literature. The results show that, with the increasing Ha number, the centerline velocity of the rectangular duct with MHD flow decreases and the absolute value of the gradient of total mechanical energy along the streamwise direction increases. The maximum of K appears near the wall in both coordinate axis of the duct. According to the energy gradient theory, this position of the maximum of K would initiate flow instability (if any) than the other positions. The higher the Hartmann number is, the smaller the K value becomes, which means that the fluid becomes more stable in the presence of higher magnetic force. As the Hartmann number increases, the K value in the parallel layer decreases more significantly than in the Hartmann layer. The most dangerous position of instability tends to migrate towards wall of the duct as the Hartmann number increases. Thus, with the energy gradient theory, the stability or instability in the rectangular duct can be controlled by modulating the magnetic force.

In this paper, the stability of inviscid parallel flow between two parallel walls is studied. Firstly, it is obtained that the profile of the base flow for this classical problem is a uniform flow. Secondly, it is shown that the solution of the disturbance equation is cr=U and ci=0, i.e., the propagation speed of the disturbance equals the flow velocity and the disturbance in this flow is neutral. Finally, it is suggested that the classical Rayleigh Theorem on inflectional velocity instability is incorrect which states that the necessary condition for instability of inviscid parallel flow is the existence of an inflection point on the velocity profile.

The distribution of energy loss due to viscosity friction in plane Couette flow and Taylor-Couette Flow between concentric rotating cylinders are studied in detail for various flow conditions. The energy loss is related to the industrial processes in some fluid delivery devices and has significant influence on the flow efficiency, flow stability, turbulent transition, mixing, and heat transfer behaviours, etc. Therefore, it is important to know about the energy loss distribution in the flow domain and to know its influence on the flow for better understanding of the flow physics. The calculation or methodology of calculating the energy loss distribution in the Taylor-Couette flow between concentric rotating cylinders is not readily found in the open literature. In this paper, the principle and the calculation are given for single cylinder rotation of either the inner or outer cylinder, and counter and same direction rotation of two cylinders. For comparison, the distribution of energy loss in a plane Couette flow is also derived for various flow conditions. Discussions of the effect of energy loss on the flow behaviour are carried out from which some findings are suggested.

The energy gradient theory is used to study the instability of Taylor-Couette flow between concentric rotating cylinders. This theory has been proposed in our previous works. In our previous studies, the energy gradient theory was demonstrated to be applicable for wall-bounded parallel flows. It was found that the critical value of the energy gradient parameter Kmax at turbulent transition is about 370-389 for wall-bounded parallel flows (which include plane Poiseuille flow, pipe Poiseuille flow and plane Couette flow) below which no turbulence occurs. In this paper, the detailed derivation for the calculation of the energy gradient parameter in the flow between concentric rotating cylinders is provided. The calculated results for the critical condition of primary instability (with semi-empirical treatment) are found to be in very good agreement with the experiments in the literature. A possible mechanism of spiral turbulence generation observed for counter-rotation of two cylinders can also be explained using the energy gradient theory. The energy gradient theory can serve to relate the condition of transition in Taylor-Couette flow to that in plane Couette flow. The latter reasonably becomes the limiting case of the former when the radii of cylinders tend to infinity. It is our contention that the energy gradient theory is possibly fairly universal for analysis of flow instability and turbulent transition, and is found valid for both pressure and shear driven flows in parallel and rotating flow configurations.

The effect of an initial disturbance on the detonation front structure in a narrow duct is studied by three-dimensional numerical simulation. The numerical method used includes a high resolution fifth-order weighted essentially non-oscillatory scheme for spatial discretization, coupled with a third order total variation diminishing Runge-Kutta time stepping method. Two types of disturbances are used for the initial perturbation. One is a random disturbance which is imposed on the whole area of the detonation front, and the other is a symmetrical disturbance imposed within a band along the diagonal direction on the front. The results show that the two types of disturbances lead to different processes. For the random disturbance, the detonation front evolves into a stable spinning detonation. For the symmetrical diagonal disturbance, the detonation front displays a diagonal pattern at an early stage, but this pattern is unstable. It breaks down after a short while and it finally evolves into a spinning detonation. The spinning detonation structure ultimately formed due to the two types of disturbances is the same. This means that spinning detonation is the most stable mode for the simulated narrow duct. Therefore, in a narrow duct, triggering a spinning detonation can be an effective way to produce a stable detonation as well as to speed up the deflagration to detonation transition process.