Particle Heating Rate Research Articles

Numerical simulation of the transient temperature distribution inside moving particles

AbstractThis work is devoted to the two‐dimensional (2D) numerical simulation of heat and fluid flow by granular mixing in a horizontally rotating kiln. The heat and fluid flow in the gas phase are solved directly using a fixed Eulerian grid. At the same time, the particle dynamics and their collisions are solved on a Lagrangian grid. The no‐slip boundary condition on the particle surface is implemented using the fictitious boundary method. The heat transfer inside the particles is calculated utilising two models: the first is the direct solution of the energy conservation equation in Lagrangian and Eulerian space and the second is the so‐called linear model, which assumes a homogeneous distribution of the temperature inside each particle. Numerical simulations showed that if the thermal diffusivity of the gas phase significantly exceeds the same parameter of the particles, the linear model over‐predicts the heating rate of the particles. The analysis of the time‐averaged flow field inside the kiln showed that in the gas phase a double vortex structure is formed which increases the convective heat transfer in the upper part of the particulate bed. The influence of the particle size, the angular velocity of the drum and the fluid on the heating rates of particles is studied and discussed.

Utilisation of waste material for environmental applications: calcination of mussel shells for waste water treatment

Waste mussel shells from the New Zealand aquaculture industry have been investigated as a potential source of lime for use in waste water treatment. The calcination of raw mussel shells to lime was studied as a function of particle size, temperature, heating rate, treatment time and atmosphere. Lime formation went through a maximum with respect to particle size for all heating conditions, as a result of competition between increased surface area and a tendency for particle to form necks resulting in caking. Optimised heat treatment conditions resulted in transformations of more than 95%. Preliminary experiments to investigate removal of phosphates from model waste water solutions showed that raw shells had limited ability due to a mechanism that relies on adsorption only. For these raw materials the finer particle sizes removed more phosphate due to increased surface area. The heat treated mussel shells on the other hand were shown to remove greater than 95% of the phosphates for both fine and coarse size ranges due to a precipitation mechanism of phosphate removal. This work demonstrates the potential for using mussels as an alternative to limestone for phosphate removal, at the same time overcoming the expensive and problematic issues of waste disposal.

Direct numerical simulation of particulate flows with heat transfer in a rotating cylindrical cavity

The purpose of this work was the direct numerical simulation of heat and fluid flow by granular mixing in a horizontal rotating kiln. To model particle behaviour and the heat and fluid flow in the drum, we solve the mass conservation, momentum and energy conservation equations directly on a fixed Eulerian grid for the whole domain including particles. At the same time the particle dynamics and their collisions are solved on a Lagrangian grid for each particle. To calculate the heat transfer inside the particles we use two models: the first is the direct solution of the energy conservation equation in the Lagrangian and Eulerian space, and the second is our so-called linear model that assumes homogeneous distribution of the temperature inside each particle. Numerical simulations showed that, if the thermal diffusivity of the gas phase significantly exceeds the same parameter of the particles, the linear model overpredicts the heating rate of the particles. The influence of the particle size and the angular velocity of the drum on the heating rates of particles is studied and discussed.

A numerical simulation of pulverized coal combustion employing a tabulated-devolatilization-process model (TDP model)

A new coal devolatilization model employing a tabulated-devolatilization-process model (TDP model) is developed, and its validity is investigated by performing a numerical simulation of a pulverized coal combustion field formed by an industrial low-NO x burner in a 100 kg-coal/h test furnace. The predicted characteristics of the pulverized coal combustion field obtained from the simulation employing the TDP model are compared with those employing the conventional devolatilization model, those employing the two competing reaction rate model, and the experiments. The results show that drastic differences in the gas flow patterns and coal particle behavior appear between simulations. In particular, the recirculation flow behavior is strongly affected by the difference in the coal devolatilization model because of the difference in the volatile matter evolution rate. The TDP model captures the observed behavior of the coal particles in the experiment better than the other models. Although it is considered that by adjusting the devolatilization parameters the prediction similar to the TDP model is also possible by the other models, appropriate devolatilization parameters are automatically set to particles depending on the particle heating rate without trial–error method by employing the TDP model.

Pyrolysis of chromated copper arsenate (CCA) treated wood waste at elevated pressure: Influence of particle size, heating rate, residence time, temperature and pressure

Lab-scale pyrolysis experiments with weathered CCA treated wood chips have been performed and the influence of particle size, residence time (10–40 min), heating rate (5–20 °C/min), temperature (330–430 °C) and pressure (0 bar, 5 bar) has been investigated. Few data, covering the pyrolysis of weathered wood was found in the literature and the literature data on pyrolysis experiments with a controlled CCA wood input, showed that results were often highly affected by experimental uncertainty. In order to reduce the uncertainty on the results, a thorough characterization of the wood input has been performed and a ratio method has been proposed which allows to study the effect of particle size on arsenic and chromium volatilization. Larger wood particles show a higher arsenic and chromium retention during pyrolysis which is attributed to the higher mass transfer resistance in these particles. Residence time has a limited effect on arsenic retentions. Increasing heating rate results in a limited increase in arsenic retentions and a more profound increase in chromium retentions. The latter is attributed to a lower average particle temperature during heating caused by the thermal lag in larger particles. Elevated pressure results in a significant increase of arsenic retentions, which is probably due to higher mass transfer resistance. Increasing temperature results in a slight decrease in arsenic retentions till 390 °C, with a sharp decrease at higher temperatures. Chromium retentions are less affected by increasing temperature, especially at higher temperatures. To conclude, a mechanism is proposed for the volatilization of chromium and arsenic during low temperature pyrolysis of CCA wood. Mass transfer resistance and the formation of As 4O 6 are crucial for the control of arsenic volatilization, while heat transfer resistance and thermal lag are more important for the control of chromium volatilization.

Coal Swelling Model for High Heating Rate Pyrolysis Applications

Thermal swelling of coal during pyrolysis strongly influences combustion and gasification rates. Coal swelling is known to vary strongly with coal rank, heating rate, and total pressure. New experiments confirm and clarify previous observations that maximum swelling occurs for high-rank bituminous coals at heating rates slightly below 104 K/s. Advanced swelling models based on bubble physics yield good qualitative trends at low heating rates, but they fail to predict observed decreases in coal swelling as particle heating rates increase beyond 104 K/s. An empirical swelling correlation for computational fluid dynamics (CFD) applications that correctly describes experimentally observed trends with the heating rate is proposed. Model parameters were fit to data from the literature using a coal rank index based on the chemical structure of the coal. The correlation has a form that allows for calculation of shrinkage for lignites, increasing swelling ratios as rank increases from sub-bituminous to bituminous ranks, and decreasing swelling ratios as rank increases from low volatile bituminous to anthracite. The correlation accurately predicts swelling as a function of coal rank and heating rate at atmospheric pressure.

Validation of a numerical approach to model pyrolysis of biomass and assessment of kinetic data

The objective of this contribution is to introduce a numerical approach that predicts pyrolysis of various biomasses under various boundary conditions such as particle geometry and heating rate. Hence, the approach is applicable to numerical tools that help exploiting further the potential of biomass as an energy resource. Furthermore, it is intended to represent the general trend of state of the art in model developments. The presented model itself is based on one-dimensional and transient differential conservation equations for mass, momentum, species and energy. In order to account for a variety of reaction schemes a formulation based on the Arrhenius equation including intrinsic modelling is chosen. It is believed to represent the majority of kinetic pyrolysis data derived from experiments. Applying this single numerical approach to predict pyrolysis under varying experimental conditions for different materials always revealed good agreement with measurements.

Pyrolysis of wheat straw in a thermogravimetric analyzer: Effect of particle size and heating rate on devolatilization and estimation of global kinetics

The influence of different parameters such as particle size, initial weight of the sample, and heating rate on the devolatilization of wheat straw particles have been studied using thermogravimetric analysis. In addition, the variations in proximate analysis with different particle sizes of wheat straw have also been investigated. Results show that the curves corresponding to the third stage of pyrolysis differ for variations in particle size, initial weight, and heating rate of the pyrolysis process. A devolatilization model with three parallel nth-order reactions has been considered to determine the global kinetic parameters using thermogravimetric data. The kinetic parameters such as activation energy (kJ/mol), frequency factor (1/min), and order of the reaction for the three stages considered in devolatilization model were E 1 = 69, E 2 = 78, E 3 = 80; k 01 = 2.57 × 10 12, k 02 = 3.97 × 10 7, k 03 = 3.17 × 10 6; and n 1 = 2.3, n 2 = 0.65, n 3 = 2.7, respectively. It was noted from the order of the reaction that the second stage of the pyrolysis curve corresponds to the degradation of cellulose and hemicellulose, and the third stage to the lignin degradation.

Advanced kinetics for calorimetric techniques and thermal stability screening of sulfide minerals

Thermal methods of analysis such as differential scanning calorimetry (DSC) provide a powerful methodology for the study of solid reactions. This paper proposes an improved thermal analysis methodology for thermal stability investigation of complex solid-state reactions. The proposed methodology is based on differential iso-conversional approach and involves peak separation, individual peak analysis and combination of isothermal/non-isothermal DSC measurements for kinetic analysis and prediction. The proposed thermal analysis, which coupled with Mineral Libration Analyzer (MLA) technique was employed to investigate thermal behavior of sulfide mineral oxidation. The importance of various experimental variables such as particle size, heating rate and atmosphere were investigated and discussed. The information gained from such an advanced thermal analysis method is useful for scale-up processes with potential of significant savings in plant operations, as well as in mitigating adverse environmental and safety issues arising from handling and storage of sulfide minerals.

Dye adsorption onto char from bamboo

Experiments have been carried out to prepare char from waste bamboo scaffolding for wastewater treatment. Carbonisation parameters such as temperature, holding time, heating rate and particle size were investigated. When the material was heated to 1173 K for 2 h, surface area ( S BET- N 2 ) and total pore volume ( V total) were 327 m 2/g and 0.185 cm 3/g, respectively. Particle size and heating rate appear less influential to the yield and textural characteristics of the resultant chars. Bamboo char obtained in this study did not show very high adsorption capacities for two acid dyes (Acid Blue 25 and Acid Yellow 117) but exhibited significant uptake of Methylene Blue. The equilibrium data were modelled by three different isotherms, namely, Langmuir, Freundlich and Redlich–Peterson model. Comparing the predicted data using optimised parameters from each model based on the SSE error function, the Redlich–Peterson isotherm appears the ‘best’ model to correlate experimental data.

Tangential force and heating rate of a neutral relativistic particle mediated by equilibrium background radiation

A relativistic theory of fluctuation-electromagnetic interaction of a small neutral polarizable particle moving with respect to the equilibrium background radiation is developed. It is assumed that the particle radius is smaller then the characteristic wave length of the background radiation. General relativistic relations for the particle heating rate and tangential force are obtained and compared with the nonrelativistic and relativistic expressions of other authors.

Temperature enhancements and vertical winds in the lower thermosphere associated with auroral heating during the DELTA campaign

[1] A coordinated observation of the atmospheric response to auroral energy input in the polar lower thermosphere was conducted during the Dynamics and Energetics of the Lower Thermosphere in Aurora (DELTA) campaign. N2 rotational temperature was measured with a rocket-borne instrument launched from the Andoya Rocket Range, neutral winds were measured from auroral emissions at 557.7 nm with a Fabry-Perot Interferometer (FPI) at Skibotn and the KEOPS, and ionospheric parameters were measured with the European Incoherent Scatter (EISCAT) UHF radar at Tromso. Altitude profiles of the passive energy deposition rate and the particle heating rate were estimated using data taken with the EISCAT radar. The local temperature enhancement derived from the difference between the observed N2 rotational temperature and the MSISE-90 model neutral temperature were 70–140 K at 110–140 km altitude. The temperature increase rate derived from the estimated heating rates, however, cannot account for the temperature enhancement below 120 km, even considering the contribution of the neutral density to the estimated heating rate. The observed upward winds up to 40 m s �1 seem to respond nearly instantaneously to changes in the heating rates. Although the wind speeds cannot be explained by the estimated heating rate and the thermal expansion hypothesis, the present study suggests that the generation mechanism of the large vertical winds must be responsible for the fast response of the vertical wind to the heating event.

In situX-ray diffraction of the transformation of gibbsite to α-alumina through calcination: effect of particle size and heating rate

A study was conducted on the gibbsite to α-alumina (α-Al2O3) transformation dynamics with particular reference to the influence of particle size and heating rate. Coarse- and fine-grained gibbsites were used to examine the transformation paths for the two materials, with particular reference to the upper and lower branch sequences described by Wefers & Misra [Oxides and Hydroxides of Aluminium, (1987), Alcoa Technical Paper No. 19 Revised, Aluminium Company of America]. The main techniques used to assess gibbsite calcination were thermogravimetric differential scanning calorimetry andin situX-ray diffraction, which provided complementary information on transformation behaviour. Whilst the quantitative phase analysis results indicated that the coarse- and fine-grained gibbsites followed both the upper and the lower branch sequences, detailed analysis of the results highlighted specific differences in the transformation behaviour for the two materials. With the loss of three molecules of water as gibbsite transformed to transition aluminas and finally α-Al2O3, there was a high degree of disorder in the crystal structure, resulting in broad and diffuse reflections in the diffraction patterns.

Fundamentals, kinetics and endothermicity of the biomass pyrolysis reaction

The paper reviews the pyrolysis of biomass constituents and possible secondary reactions. Biomass pyrolysis yields mostly liquid and solid fuel, easy to store and transport. Relevant working conditions for experiments and large-scale operation are: (i) biomass particles < 200 μm, (ii) a particle heating rate of at least about 80 K min −1 and (iii) a reactor environment where the internal resistance to heat penetration is smaller than the external resistance to heat transfer (Biot-number, Bi < 1). The circumstances of TGA and DSC experiments meet these requirements and fully determine the reaction kinetics and endothermicity of the pyrolysis reaction. The reaction rate constant and the heat of reaction are essential parameters in the design of a pyrolysis reactor. For most of the biomass species tested, the first order reaction rate constant is large and >0.5 s −1. The heat of reaction ranges from 207 to 434 kJ kg −1. All results tie in with literature data, although the reader is cautioned in using literature data since experiments were not always performed under relevant testing conditions.

Electron microscopy studies of the thermal stability of gold nanoparticle arrays

A series of monolayer protected gold nanoparticle colloidal solutions have been prepared with average sizes in the 2–15nm range. If a drop of such a colloidal suspension is deposited onto a Si3N4 substrate and the solvent allowed to evaporate, the particles have a tendency to self-assemble into monolayer rafts with varying degrees of structural order depending on the initial mono-dispersity of the particles. The thermal stability of these selfassembled gold nanoparticle rafts as a function of particle size, heating method, heating rate and ligand identity have been assessed in this study. In-situ TEM studies show that sub-8nm Au nanoparticles on Si3N4 have a tendency to coarsen upon slow heating, whereas those comprised of larger particles exhibit densification. Increasing the heating rate for the smaller particles promoted densification, forcing them to form highly interconnected string-like structures. Finally, rafts of sub-4nm alkanethiol protected Au nanoparticles are shown to sinter spontaneously under ambient conditions at room temperature on the timescale of several months. This unexpected effect may have important implications for the long term structural stability of any device constructed from sub-4nm gold nanoparticles.

Thermogravimetric analysis and pyrolysis of waste mixtures of paint and tar slag

We describe thermogravimetric analyses and pyrolysis kinetic studies carried out on hazardous waste mixtures of tar slag, paint slag, paper, sodium sulfate and calcium oxide. Both thermogravimetric (TG) and differential thermogravimetric (DTG) profiles were measured by a thermogravimetric analyzer at different final temperatures, particle sizes and heating rates. Pyrolysis kinetic parameters were calculated by the Coats-Redfern method. Influences of particle size, heating rate and final temperature on pyrolysis yields and kinetic parameters are also discussed. The results show that final temperature and particle size have a great effect on pyrolysis yields. We find that with increasing temperature the activation energy initially increases to a maximum value and then decreases.

Size-dependant heating rates of iron oxide nanoparticles for magnetic fluid hyperthermia

Using the thermal decomposition of organometallics method we have synthesized high-quality, iron oxide nanoparticles of tailorable size up to ∼15 nm and transferred them to a water phase by coating with a biocompatible polymer. The magnetic behavior of these particles was measured and fit to a log-normal distribution using the Chantrell method and their polydispersity was confirmed to be very narrow. By performing calorimetry measurements with these monodisperse particles we have unambiguously demonstrated, for the first time, that at a given frequency, heating rates of superparamagnetic particles are dependent on particle size, in agreement with earlier theoretical predictions.

Effects of particle shape and size on devolatilization of biomass particle

Experimental and theoretical investigations indicate particle shape and size influence biomass particle dynamics, including drying, heating rate, and reaction rate. Experimental samples include disc/flake-like, cylindrical/cylinder-like, and equant (nearly spherical) shapes of wood particles with similar particle masses and volumes but different surface areas. Small samples (320 μm) passed through a laboratory entrained-flow reactor in a nitrogen atmosphere and a maximum reactor wall temperature of 1600 K. Large samples were suspended in the center of a single-particle reactor. Experimental data indicate that equant particles react more slowly than the other shapes, with the difference becoming more significant as particle mass or aspect ratio increases and reaching a factor of two or more for particles with sizes over 10 mm. A one-dimensional, time-dependent particle model simulates the rapid pyrolysis process of particles with different shapes. The model characterizes particles in three basic shapes (sphere, cylinder, and flat plate). With the particle geometric information (particle aspect ratio, volume, and surface area) included, this model simulates the devolatilization process of biomass particles of any shape. Model simulations of the three shapes show satisfactory agreement with the experimental data. Model predictions show that both particle shape and size affect the product yield distribution. Near-spherical particles exhibit lower volatile and higher tar yields relative to aspherical particles with the same mass under similar conditions. Volatile yields decrease with increasing particle size for particles of all shapes. Assuming spherical or isothermal conditions for biomass particles leads to large errors at most biomass particle sizes of practical interest.

Evolution of biomass char structure during oxidation in O2 as revealed with FT-Raman spectroscopy

FT-Raman spectroscopy has been used to identify structural features and evaluate the structural evolution of biomass chars during gasification with air. Chars prepared from the pyrolysis of a cane trash sample with a fast particle heating rate in a novel fluidised-bed/fixed-bed reactor at 500, 700 and 900 °C were oxidised at 400 °C in air in a TGA. The data derived from the spectral deconvolution of Fourier Transform — Raman spectra suggest that the 500 °C char showed very different structural features after pyrolysis and during oxidation from the 700 and 900 °C chars, while the differences between the latter two chars were small. Preferential consumption by O 2 of smaller aromatic rings and structures of somewhat aliphatic characteristics left the char more enriched with larger aromatic ring systems. The changes in char structure are in agreement with the observed reactivity measured in O 2 in a thermogravimetric analyser.

Study of Thermal Degradation Kinetics of Elastomeric Powder (Ground Tire Rubber)

This paper investigates the use of thermogravimetric analysis to study the kinetics of thermal degradation of ground tire rubber (GTR) obtained from reclaimed tires. We analyzed the composition of GTR and determined the content percentage of its components (volatile compounds, rubber, and ash) using thermogravimetric analysis (TGA). T 1 and T 2 temperature peaks corresponding to the maximum normalized weight loss rate (NWLR) of the two main rubber components of tires NR (natural rubber) and a BR/SBR (butadiene/styrene-butadiene rubber blend) using the Gaussian deconvolution method has been determined. The influence of GTR particle size and heating rate on reaction rates (dx/dt) and on degradation time during the degradation process has been analyzed.