A generalised dynamic model for char particle gasification with structure evolution and peripheral fragmentation

A generalised model for the prediction of single char particle gasification dynamics, accounting for multi-component mass transfer with chemical reaction, heat transfer, as well as structure evolution and peripheral fragmentation is developed in this paper. Maxwell-Stefan analysis is uniquely applied to both micro and macropores within the framework of the dusty-gas model to account for the bidisperse nature of the char, which differs significantly from the conventional models that are based on a single pore type. The peripheral fragmentation and random-pore correlation incorporated into the model enable prediction of structure/reactivity relationships. The occurrence of chemical reaction within the boundary layer reported by Biggs and Agarwal (Chem. Eng. Sci. 52 (1997) 941) has been confirmed through an analysis of CO/CO2 product ratio obtained from model simulations. However, it is also quantitatively observed that the significance of boundary layer reaction reduces notably with the reduction of oxygen concentration in the flue gas, operational pressure and film thickness. Computations have also shown that in the presence of diffusional gradients peripheral fragmentation occurs in the early stages on the surface, after which conversion quickens significantly due to small particle size. Results of the early commencement of peripheral fragmentation at relatively low overall conversion obtained from a large number of simulations agree well with experimental observations reported by Feng and Bhatia (Energy & Fuels 14 (2000) 297). Comprehensive analysis of simulation results is carried out based on well accepted physical principles to rationalise model prediction. (C) 2001 Elsevier Science Ltd. AH rights reserved.

Propagation of nonstationary curved and stretched premixed flames with multistep reaction mechanisms

The propagation speed of a thin premixed flame disturbed by an unsteady fluid flow of a larger scale is considered. The flame may also have a general shape but the reaction zone is assumed to be thin compared to the flame thickness. Unlike in preceding publications, the presented asymptotic analysis is performed for a general multistep reaction mechanism and, at the same time, the flame front is curved by the fluid flow. The resulting equations define the propagation speed of disturbed flames in terms of the properties of undisturbed planar flames and the flame stretch. Special attention is paid to the near-equidiffusion limit. In this case, the flame propagation speed is shown to depend on the effective Zeldovich number Z(f) , and the flame stretch. Unlike the conventional Zeldovich number, the effective Zeldovich number is not necessarily linked directly to the activation energies of the reactions. Several examples of determining the effective Zeldovich number for reduced combustion mechanisms are given while, for realistic reactions, the effective Zeldovich number is determined from experiments. Another feature of the present approach is represented by the relatively simple asymptotic technique based on the adaptive generalized curvilinear system of coordinates attached to the flame (i.e., intrinsic disturbed flame equations [IDFE]).

Recent work on gaseous detonations

The paper reviews recent progress in the field of gaseous detonations, with sections on shock diffraction and reflection, the transition to detonation, hybrid, spherically-imploding, and galloping and stuttering fronts, their structure, their transmission and quenching by additives, the critical energy for initiation and detonation of more unusual fuels. The final section points out areas where our understanding is still far from being complete and contains some suggestions of ways in which progress might be made.

Efficient time-independent wave packet scattering calculations within a Lanczos subspace: H+O-2 (J=0) state-to-state reaction probabilities

An efficient Lanczos subspace method has been devised for calculating state-to-state reaction probabilities. The method recasts the time-independent wave packet Lippmann-Schwinger equation [Kouri , Chem. Phys. Lett. 203, 166 (1993)] inside a tridiagonal (Lanczos) representation in which action of the causal Green's operator is affected easily with a QR algorithm. The method is designed to yield all state-to-state reaction probabilities from a given reactant-channel wave packet using a single Lanczos subspace; the spectral properties of the tridiagonal Hamiltonian allow calculations to be undertaken at arbitrary energies within the spectral range of the initial wave packet. The method is applied to a H+O-2 system (J=0), and the results indicate the approach is accurate and stable. (C) 2002 American Institute of Physics.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

Upstream porthole injection in a 2-D scramjet model

Injection from portholes upstream of the combustion chamber was investigated as a method of delivering fuel into a scramjet. This method reduces the viscous drag on a model by allowing a reduction in the length of the combustion chamber. At experimental enthalpies of 3.0 MJ/kg in the T4 shock tunnel, there was no evidence of combustion in the intake, either by shadowgraph or pressure measurements. Combustion was observed in the combustion chamber. A theoretical extension of these results is made to a hot wall scenario.

Role of oxygen in the nitrous oxide/carbon reaction

An investigation of the role of oxygen in the nitrous oxide/carbon reaction was carried out on various carbon samples (both graphitic and nongraphitic) over a range of temperatures and partial pressures. Previous work reported that oxygen strongly inhibited the nitrous oxide/carbon reaction. Large ratios of O-2/N2O were used in all previous work. In this work, the O-2/N2O ratio was kept below 1, and we found that oxygen did not inhibit the rate of the C + N2O reaction. Instead, the rate of the reaction in the presence of oxygen was essentially that predicted by the two independent reactions, nitrous oxide/carbon and oxygen/carbon, occurring simultaneously. A simple theoretical explanation is given for the observations, both past and present, on the basis of competitive chemisorption of nitrous oxide and oxygen on active sites.

Assessing the self-heating behaviour of Callide coal using a 2-metre column

A 2-m, adiabatic column has been successfully refurbished and recommissioned for coal self-heating research at The University of Queensland. Subbituminous coal from the Callide Coalfields reached thermal runaway in just under 19 days from a starting temperature of 20-22 degreesC. The coal was loaded as two layers, with an R-70 index of 2.73 degreesC h(-1) and 5.90 degreesC h(-1) for the upper layer and lower layer respectively. Initially, a hotspot developed in the upper layer between 120 and 140 cm from the air inlet due to moisture adsorption. After 7 days, self-heating in the lower half of the column began to take over, consistent with the higher R-70 index of this coal. The location of the final hotspot was approximately 60 cm from the air inlet. Further tests on Australian coals, with the column, will enable a better understanding of coal self-heating under conditions closely resembling mining, transport and storage of coal. The results from the column will also provide industry with the information needed to manage the coal self-heating hazard. (C) 2002 Elsevier Science Ltd. All rights reserved.

Structural ordering of coal char during heat treatment and its impact on reactivity

The effect of heat treatment on the structure of an Australian semi-anthracite char was studied in detail in the 850-1150degreesC temperature range using XRD, HRTEM, and electrical resistivity techniques. It was found that the carbon crystallite size in the char does not change significantly during heat treatment in the temperature range studied, for both the raw coal and its ash-free derivative obtained by acid treatment. However, the fraction of the organized carbon in the raw coal chars, determined by XRD, increased with increase of heat treatment time and temperature, while that for the ash-free coal chars remained almost unchanged. This suggests the occurrence of catalytic ordering during heat treatment, supported by the observation that the electrical resistivity of the raw coal chars decreased with heat treatment, while that of the ash-free coal chars did not vary significantly. Further confirmatory evidence was provided by high resolution transmission electron micrographs depicting well-organized carbon layers surrounding iron particles. It is also found that the fraction of organized carbon does not reach unity, but attains an apparent equilibrium value that increases with increase in temperature, providing an apparent heat of ordering of 71.7 kJ mol(-1) in the temperature range studied. Good temperature-independent correlation was found between the electrical resistivity and the organized carbon fraction, indicating that electrical resistivity is indeed structure sensitive. Good correlation was also found between the electrical resistivity and the reactivity of coal char. All these results strongly suggest that the thermal deactivation is the result of a crystallite-perfecting process, which is effectively catalyzed by the inorganic matter in the coal char. Based on kinetic interpretation of the data it is concluded that the process is diffusion controlled, most likely involving transport of iron in the inter-crystallite nanospaces in the temperature range studied. The activation energy of this transport process is found to be very low, at about 11.8 kJ mol(-1), which is corroborated by model-free correlation of the temporal variation of organized carbon fraction as well as electrical resistivity data using the superposition method, and is suggestive of surface transport of iron. (C) 2002 Elsevier Science Ltd. All rights reserved.

gamma-alumina nanofibers prepared from aluminum hydrate with poly(ethylene oxide) surfactant

Introducing poly(ethylene oxide) surfactant to aluminum hydrate colloids can effectively direct the crystal growth of boehmite and the crystal morphology of final gamma-alumina crystallites. Fibrous crystallites of gamma-alumina about 3-4 nm thick and 30-60 nm long are obtained. They stack randomly, resulting in a structure with a low contact area between the fibers but with a very large porosity. Such a structure exhibits strong resistance to sintering when heated to high temperatures. A sample retains a BET surface area of 68 m(2)/g, after being heated to 1473 K. The surfactant molecules form micelles that interact with the colloid particles of aluminum hydroxide through hydrogen bonding. This interaction is not sufficient to change the intrinsic crystal structure of boehmite, but induces profound changes in the morphology of boehmite crystallites and their growth. The surfactant-induced fiber formation (SIFF) process has distinct features from templated synthesis but shows similarities in some respects to biomineralization processes in which inorganic crystals with complex morphological shapes can be formed in biological systems. SIFF offers an effective approach to create new nanostructures of inorganic oxide from aqueous media.

Blast design using measurement while drilling parameters

Measurement while drilling (MWD) techniques can provide a useful tool to aid drill and blast engineers in open cut mining. By avoiding time consuming tasks such as scan-lines and rock sample collection for laboratory tests, MWD techniques can not only save time but also improve the reliability of the blast design by providing the drill and blast engineer with the information specially tailored for use. While most mines use a standard blast pattern and charge per blasthole, based on a single rock factor for the entire bench or blast region, information derived from the MWD parameters can improve the blast design by providing more accurate rock properties for each individual blasthole. From this, decisions can be made on the most appropriate type and amount of explosive charge to place in a per blasthole or to optimise the inter-hole timing detonation time of different decks and blastholes. Where real-time calculations are feasible, the system could extend the present blast design even be used to determine the placement of subsequent holes towards a more appropriate blasthole pattern design like asymmetrical blasting.

Asymmetric blasting: a rock mass dependent blast design method

Blasting has been the most frequently used method for rock breakage since black powder was first used to fragment rocks, more than two hundred years ago. This paper is an attempt to reassess standard design techniques used in blasting by providing an alternative approach to blast design. The new approach has been termed asymmetric blasting. Based on providing real time rock recognition through the capacity of measurement while drilling (MWD) techniques, asymmetric blasting is an approach to deal with rock properties as they occur in nature, i.e., randomly and asymmetrically spatially distributed. It is well accepted that performance of basic mining operations, such as excavation and crushing rely on a broken rock mass which has been pre conditioned by the blast. By pre-conditioned we mean well fragmented, sufficiently loose and with adequate muckpile profile. These muckpile characteristics affect loading and hauling [1]. The influence of blasting does not end there. Under the Mine to Mill paradigm, blasting has a significant leverage on downstream operations such as crushing and milling. There is a body of evidence that blasting affects mineral liberation [2]. Thus, the importance of blasting has increased from simply fragmenting and loosing the rock mass, to a broader role that encompasses many aspects of mining, which affects the cost of the end product. A new approach is proposed in this paper which facilitates this trend 'to treat non-homogeneous media (rock mass) in a non-homogeneous manner (an asymmetrical pattern) in order to achieve an optimal result (in terms of muckpile size distribution).' It is postulated there are no logical reasons (besides the current lack of means to infer rock mass properties in the blind zones of the bench and onsite precedents) for drilling a regular blast pattern over a rock mass that is inherently heterogeneous. Real and theoretical examples of such a method are presented.

Critical comparison of Julius Kruttschnitt Mineral Research Centre (JKMRC) blast fragmentation models

Blast fragmentation can have a significant impact on the profitability of a mine. An optimum run of mine (ROM) size distribution is required to maximise the performance of downstream processes. If this fragmentation size distribution can be modelled and controlled, the operation will have made a significant advancement towards improving its performance. Blast fragmentation modelling is an important step in Mine to Mill™ optimisation. It allows the estimation of blast fragmentation distributions for a number of different rock mass, blast geometry, and explosive parameters. These distributions can then be modelled in downstream mining and milling processes to determine the optimum blast design. When a blast hole is detonated rock breakage occurs in two different stress regions - compressive and tensile. In the-first region, compressive stress waves form a 'crushed zone' directly adjacent to the blast hole. The second region, termed the 'cracked zone', occurs outside the crush one. The widely used Kuz-Ram model does not recognise these two blast regions. In the Kuz-Ram model the mean fragment size from the blast is approximated and is then used to estimate the remaining size distribution. Experience has shown that this model predicts the coarse end reasonably accurately, but it can significantly underestimate the amount of fines generated. As part of the Australian Mineral Industries Research Association (AMIRA) P483A Mine to Mill™ project, the Two-Component Model (TCM) and Crush Zone Model (CZM), developed by the Julius Kruttschnitt Mineral Research Centre (JKMRC), were compared and evaluated to measured ROM fragmentation distributions. An important criteria for this comparison was the variation of model results from measured ROM in the-fine to intermediate section (1-100 mm) of the fragmentation curve. This region of the distribution is important for Mine to Mill™ optimisation. The comparison of modelled and Split ROM fragmentation distributions has been conducted in harder ores (UCS greater than 80 MPa). Further work involves modelling softer ores. The comparisons will be continued with future site surveys to increase confidence in the comparison of the CZM and TCM to Split results. Stochastic fragmentation modelling will then be conducted to take into account variation of input parameters. A window of possible fragmentation distributions can be compared to those obtained by Split . Following this work, an improved fragmentation model will be developed in response to these findings.

The modeling of turbulent reactive flows based on multiple mapping conditioning

A new modeling approach-multiple mapping conditioning (MMC)-is introduced to treat mixing and reaction in turbulent flows. The model combines the advantages of the probability density function and the conditional moment closure methods and is based on a certain generalization of the mapping closure concept. An equivalent stochastic formulation of the MMC model is given. The validity of the closuring hypothesis of the model is demonstrated by a comparison with direct numerical simulation results for the three-stream mixing problem. (C) 2003 American Institute of Physics.

Kinetic parameters associated with self-heating of New Zealand coals under adiabatic conditions

Adiabatic self-heating tests were carried out on five New Zealand coal samples ranging in rank from lignite to high-volatile bituminous. Kinetic parameters of oxidation were obtained front the self-heating curves assuming Arrhenius behaviour. The activation energy E (kJ mol(-1)) and the pre-exponential factor A (s(-1)) were determined in the temperature range of 70-140 degreesC. The activation energy exhibited a definite rank relationship with a minimum E of 55 kJ mol(-1) occurring at a Suggate rank of similar to6.2 corresponding to subbituminous C. Either side of this rank there was a noticeable increase in the activation energy indicating lower reactivity of the coal. A similar rank trend was also observed in the R-70 self-heating rate index values that were taken from the initial portion of the self-heating curve front 40 to 70 degreesC. From these results it is clear that the adiabatic method is capable of providing reliable kinetic parameters of coal oxidation.