CFD modelling of the fast pyrolysis of an in-flight cellulosic particle subjected to convective heat transfer

The pyrolysis of a freely moving cellulosic particle inside a 41.7mgs -1 source continuously fed fluid bed reactor subjected to convective heat transfer is modelled. The Lagrangian approach is adopted for the particle tracking inside the reactor, while the flow of the inert gas is treated with the standard Eulerian method for gases. The model incorporates the thermal degradation of cellulose to char with simultaneous evolution of gases and vapours from discrete cellulosic particles. The reaction kinetics is represented according to the Broido–Shafizadeh scheme. The convective heat transfer to the surface of the particle is solved by two means, namely the Ranz–Marshall correlation and the limit case of infinitely fast external heat transfer rates. The results from both approaches are compared and discussed. The effect of the different heat transfer rates on the discrete phase trajectory is also considered.

Mechanisms of convective drying in thick beds of solids

The literature on heat and mass transfer mechanisms in the convective drying of thick beds of solids has been critically reviewed. Related mathematical models of heat transfer are also considered. Experimental and theoretical studies were made of the temperature distribution within beds, and of drying rates, with various materials undergoing convective drying. The experimental work covered thick beds of hygroscopic and non-hygroscopic materials (glass beads of different diameters, polystyrene pellets, activated alumina and wood powder) at air temperatures of 54°C to 84°C. Tests were carried out in a laboratory drying apparatus comprising a wind tunnel through which the air, of controlled temperature and humidity, was passed over a sample suspended from a balance. Thermocouples were inserted at different depths within the sample bed. The temperature distribution profiles for both hygroscopic and non-hygroscopic beds exhibited a clear difference between the temperatures at the surface and bottom during the constant rate period. An effective method was introduced for predicting the critical moisture content. During the falling rate the profiles showed the existence of a receding evaporation plane; this divided the system into a hotter dry zone in the upper section and a wet zone near the bottom. A graphical procedure was established to predict accurately the position of the receding evaporation front at any time. A new mathematical model, based on the receding evaporation front phenomenon, was proposed to predict temperature distributions throughout a bed during drying. Good agreement was obtained when the model was validated by comparing its predictions with experimental data. The model was also able to predict the duration of each drying stage. In experiments using sample trays of different diameters, the drying rate was found to increase with a decrease in the effective length of the bed surface. During the constant rate period with trays of a small effective length, i.e. less than 0.08 m, an 'inversion' in temperature distribution occurred in the bed; the bottom temperature increased and became greater than that of the surface. Experimental measurements were verified in several ways to ensure this phenomenon was real. Theoretical explanations are given for both the effective length and temperature inversion phenomena.

Transition in convective flows heated internally

The stability of internally heated convective flows in a vertical channel under the influence of a pressure gradient and in the limit of small Prandtl number is examined numerically. In each of the cases studied the basic flow, which can have two inflection points, loses stability at the critical point identified by the corresponding linear analysis to two-dimensional states in a Hopf bifurcation. These marginal points determine the linear stability curve that identifies the minimum Grashof number (based on the strength of the homogeneous heat source), at which the two-dimensional periodic flow can bifurcate. The range of stability of the finite amplitude secondary flow is determined by its (linear) stability against three-dimensional infinitesimal disturbances. By first examining the behavior of the eigenvalues as functions of the Floquet parameters in the streamwise and spanwise directions we show that the secondary flow loses stability also in a Hopf bifurcation as the Grashof number increases, indicating that the tertiary flow is quasi-periodic. Secondly the Eckhaus marginal stability curve, that bounds the domain of stable transverse vortices towards smaller and larger wavenumbers, but does not cause a transition as the Grashof number increases, is also given for the cases studied in this work.