Non-linear dynamics of a two phase flow system in an evaporator: The effects of (i) a time varying pressure drop (ii) an axially varying heat flux Original Research Article

In this paper we study the phenomena of density wave oscillations (DWO) in a vertical heated channel. The homogeneous equilibrium model is used to simulate the flow in the two-phase region. The equations are solved numerically using a `shooting-methodʹ technique. This in its turn employs an implicit backward finite difference scheme. The scheme can incorporate the movement of the interface. It is very elegant and does not involve storage of variables in large N×N matrices. This scheme is sufficiently general and can be used to simulate the dynamic behaviour when: (i) the heat flux imposed at the surface is non-constant, i.e. exhibits an axial variation; and (ii) the imposed pressure drop is varied periodically at a fixed frequency. A possible explanation for the conflicting reports of the effect of a periodic variation in heat flux is provided using a linear stability analysis and the D-partition method. The interaction of the natural frequency of the DWO and the fixed forcing frequency of the imposed pressure drop gives rise to various phenomena viz relaxation oscillations, sub-harmonic oscillations, quasi-periodic and chaotic solutions. To aid the experimentalist describe this infinite-dimensional system on the basis of his experimental results we discuss the characterisation using only the velocity time series data. This is done employing the method of delay coordinate embedding. The phase portraits, stroboscopic map and correlation dimension of the actual attractor are compared with that of the reconstructed attractor from the velocity time series.