A carbon activation model with application to longan seed char gasification

In this paper a new structural model is presented to describe the evolution of porosity of char during the gasification process. The model assumes the char structure to be composed of bundles of parallel graphite layers, and the reactivities of each layer with the gasification agent are assumed to be different to represent the different degree of heterogeneity of each layer (i.e. each layer will react with the gasification agent at a different rate). It is this difference in the reactivity that allows micropores to be created during the course of gasification. This simple structural model enables the evolution of pore volume, pore geometrical surface area and the pore size distribution to be described with respect to the extent of char burn-off. The model is tested against the experimental data of gasification of longan seed-derived char with carbon dioxide and it is found that the agreement between the model and the data is reasonably satisfactory, especially the evolution of surface area and pore volume with burn-off.

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.

Ab initio modelling of basal plane oxidation of graphenes and implications for modelling char combustion

Ab initio calculations have been performed to determine the energetics of oxygen atoms adsorbed onto graphene planes and the possible reaction path extracting carbon atorns in the form of carbon monoxide. Front the energetics it is confirmed that this reaction path will not significantly contribute to the gasification of well ordered carbonaceous chars. Modelling results which explore this limit Lire presented. (C) 2002 Elsevier Science Ltd, All rights reserved.

On the validity of thermogravimetric determination of carbon gasification kinetics

Thermogravimetric analysis has been widely applied in kinetic studies of carbon gasification, with the associated temporal weight change profiles being used to extract kinetic information and to validate gasification models. However the weight change profiles are not always governed by the intrinsic gasification activity because of the effect of chemisorption and its dynamics. In the present work we theoretically determine the criteria under which weight change profiles can be used to determine intrinsic kinetics for CO2 and O2 gasification by examining the region in which the chemisorption dynamics can be assumed pseudo-steady. It is found that the validity of the pseudo-steady assumption depends on the experimental conditions as well as on the initial surface area of carbon. Based on known mechanisms and rate constants an active surface area region is identified within which the steady state assumption is valid and the effect of chemisorption dynamics is negligible. The size of the permissible region is sensitive to the reaction temperature and gas pressure. The results indicate that in some cases the thermogravimetric data should be used with caution in kinetic studies. A large amount of literature on thermogravimetric analyzer determined char gasification kinetics is examined and the importance of chemisorption dynamics for the data assessed.

Biomass Gasification - Process analysis and dimensioning aspects for downdraft units and gas cleaning lines

In such territories where food production is mostly scattered in several small / medium size or even domestic farms, a lot of heterogeneous residues are produced yearly, since farmers usually carry out different activities in their properties. The amount and composition of farm residues, therefore, widely change during year, according to the single production process periodically achieved. Coupling high efficiency micro-cogeneration energy units with easy handling biomass conversion equipments, suitable to treat different materials, would provide many important advantages to the farmers and to the community as well, so that the increase in feedstock flexibility of gasification units is nowadays seen as a further paramount step towards their wide spreading in rural areas and as a real necessity for their utilization at small scale. Two main research topics were thought to be of main concern at this purpose, and they were therefore discussed in this work: the investigation of fuels properties impact on gasification process development and the technical feasibility of small scale gasification units integration with cogeneration systems. According to these two main aspects, the present work was thus divided in two main parts. The first one is focused on the biomass gasification process, that was investigated in its theoretical aspects and then analytically modelled in order to simulate thermo-chemical conversion of different biomass fuels, such as wood (park waste wood and softwood), wheat straw, sewage sludge and refuse derived fuels. The main idea is to correlate the results of reactor design procedures with the physical properties of biomasses and the corresponding working conditions of gasifiers (temperature profile, above all), in order to point out the main differences which prevent the use of the same conversion unit for different materials. At this scope, a gasification kinetic free model was initially developed in Excel sheets, considering different values of air to biomass ratio and the downdraft gasification technology as particular examined application. The differences in syngas production and working conditions (process temperatures, above all) among the considered fuels were tried to be connected to some biomass properties, such elementary composition, ash and water contents. The novelty of this analytical approach was the use of kinetic constants ratio in order to determine oxygen distribution among the different oxidation reactions (regarding volatile matter only) while equilibrium of water gas shift reaction was considered in gasification zone, by which the energy and mass balances involved in the process algorithm were linked together, as well. Moreover, the main advantage of this analytical tool is the easiness by which the input data corresponding to the particular biomass materials can be inserted into the model, so that a rapid evaluation on their own thermo-chemical conversion properties is possible to be obtained, mainly based on their chemical composition A good conformity of the model results with the other literature and experimental data was detected for almost all the considered materials (except for refuse derived fuels, because of their unfitting chemical composition with the model assumptions). Successively, a dimensioning procedure for open core downdraft gasifiers was set up, by the analysis on the fundamental thermo-physical and thermo-chemical mechanisms which are supposed to regulate the main solid conversion steps involved in the gasification process. Gasification units were schematically subdivided in four reaction zones, respectively corresponding to biomass heating, solids drying, pyrolysis and char gasification processes, and the time required for the full development of each of these steps was correlated to the kinetics rates (for pyrolysis and char gasification processes only) and to the heat and mass transfer phenomena from gas to solid phase. On the basis of this analysis and according to the kinetic free model results and biomass physical properties (particles size, above all) it was achieved that for all the considered materials char gasification step is kinetically limited and therefore temperature is the main working parameter controlling this step. Solids drying is mainly regulated by heat transfer from bulk gas to the inner layers of particles and the corresponding time especially depends on particle size. Biomass heating is almost totally achieved by the radiative heat transfer from the hot walls of reactor to the bed of material. For pyrolysis, instead, working temperature, particles size and the same nature of biomass (through its own pyrolysis heat) have all comparable weights on the process development, so that the corresponding time can be differently depending on one of these factors according to the particular fuel is gasified and the particular conditions are established inside the gasifier. The same analysis also led to the estimation of reaction zone volumes for each biomass fuel, so as a comparison among the dimensions of the differently fed gasification units was finally accomplished. Each biomass material showed a different volumes distribution, so that any dimensioned gasification unit does not seem to be suitable for more than one biomass species. Nevertheless, since reactors diameters were found out quite similar for all the examined materials, it could be envisaged to design a single units for all of them by adopting the largest diameter and by combining together the maximum heights of each reaction zone, as they were calculated for the different biomasses. A total height of gasifier as around 2400mm would be obtained in this case. Besides, by arranging air injecting nozzles at different levels along the reactor, gasification zone could be properly set up according to the particular material is in turn gasified. Finally, since gasification and pyrolysis times were found to considerably change according to even short temperature variations, it could be also envisaged to regulate air feeding rate for each gasified material (which process temperatures depend on), so as the available reactor volumes would be suitable for the complete development of solid conversion in each case, without even changing fluid dynamics behaviour of the unit as well as air/biomass ratio in noticeable measure. The second part of this work dealt with the gas cleaning systems to be adopted downstream the gasifiers in order to run high efficiency CHP units (i.e. internal engines and micro-turbines). Especially in the case multi–fuel gasifiers are assumed to be used, weightier gas cleaning lines need to be envisaged in order to reach the standard gas quality degree required to fuel cogeneration units. Indeed, as the more heterogeneous feed to the gasification unit, several contaminant species can simultaneously be present in the exit gas stream and, as a consequence, suitable gas cleaning systems have to be designed. In this work, an overall study on gas cleaning lines assessment is carried out. Differently from the other research efforts carried out in the same field, the main scope is to define general arrangements for gas cleaning lines suitable to remove several contaminants from the gas stream, independently on the feedstock material and the energy plant size The gas contaminant species taken into account in this analysis were: particulate, tars, sulphur (in H2S form), alkali metals, nitrogen (in NH3 form) and acid gases (in HCl form). For each of these species, alternative cleaning devices were designed according to three different plant sizes, respectively corresponding with 8Nm3/h, 125Nm3/h and 350Nm3/h gas flows. Their performances were examined on the basis of their optimal working conditions (efficiency, temperature and pressure drops, above all) and their own consumption of energy and materials. Successively, the designed units were combined together in different overall gas cleaning line arrangements, paths, by following some technical constraints which were mainly determined from the same performance analysis on the cleaning units and from the presumable synergic effects by contaminants on the right working of some of them (filters clogging, catalysts deactivation, etc.). One of the main issues to be stated in paths design accomplishment was the tars removal from the gas stream, preventing filters plugging and/or line pipes clogging At this scope, a catalytic tars cracking unit was envisaged as the only solution to be adopted, and, therefore, a catalytic material which is able to work at relatively low temperatures was chosen. Nevertheless, a rapid drop in tars cracking efficiency was also estimated for this same material, so that an high frequency of catalysts regeneration and a consequent relevant air consumption for this operation were calculated in all of the cases. Other difficulties had to be overcome in the abatement of alkali metals, which condense at temperatures lower than tars, but they also need to be removed in the first sections of gas cleaning line in order to avoid corrosion of materials. In this case a dry scrubber technology was envisaged, by using the same fine particles filter units and by choosing for them corrosion resistant materials, like ceramic ones. Besides these two solutions which seem to be unavoidable in gas cleaning line design, high temperature gas cleaning lines were not possible to be achieved for the two larger plant sizes, as well. Indeed, as the use of temperature control devices was precluded in the adopted design procedure, ammonia partial oxidation units (as the only considered methods for the abatement of ammonia at high temperature) were not suitable for the large scale units, because of the high increase of reactors temperature by the exothermic reactions involved in the process. In spite of these limitations, yet, overall arrangements for each considered plant size were finally designed, so that the possibility to clean the gas up to the required standard degree was technically demonstrated, even in the case several contaminants are simultaneously present in the gas stream. Moreover, all the possible paths defined for the different plant sizes were compared each others on the basis of some defined operational parameters, among which total pressure drops, total energy losses, number of units and secondary materials consumption. On the basis of this analysis, dry gas cleaning methods proved preferable to the ones including water scrubber technology in al of the cases, especially because of the high water consumption provided by water scrubber units in ammonia adsorption process. This result is yet connected to the possibility to use activated carbon units for ammonia removal and Nahcolite adsorber for chloride acid. The very high efficiency of this latter material is also remarkable. Finally, as an estimation of the overall energy loss pertaining the gas cleaning process, the total enthalpy losses estimated for the three plant sizes were compared with the respective gas streams energy contents, these latter obtained on the basis of low heating value of gas only. This overall study on gas cleaning systems is thus proposed as an analytical tool by which different gas cleaning line configurations can be evaluated, according to the particular practical application they are adopted for and the size of cogeneration unit they are connected to.

Downdraft gasification of biomass

The objectives of this research were to investigate the parameters affecting the gasification process within downdraft gasifiers using biomass feedstocks. In addition to investigations with an open-core gasifier, a novel open-topped throated gasifier was designed and used. A sampling system was designed and installed to determine the water, tar and particular content of the raw product gas. This permitted evaluation of the effects of process parameters and reactor design on tar and particular production, although a large variation was found for the particulate measurements due to the capture of large particles. For both gasifiers, the gasification process was studied in order to identify and compare the mechanisms controlling the position and shape of the reaction zones. The stability of the reaction zone was found to be governed by the superficial gas velocity within the reactor. A superficial gas velocity below 0.2 Nms-1 resulted in a rising reaction zone in both gasifiers. Turndown is achieved when the rate of char production by flaming pyrolysis equals the rate of char gasification over a range of throughputs. A turndown ratio of 2:1 was achieved for the hybrid-throated gasifier, compared to 1.3:1 for the open-core. It is hypothesized that pyrolysis is a surface area phenomenon, and that in the hybrid gasifier the pyrolysis front can expand to form a dome-shape. The rate of char gasification is believed to increase as the depth of the gasification zone increases. Vibration of the open-core reactor bed decreased the bed pressure drop, reduced the voidage, aided solids flow and gave a minor improvement in the product gas energy content. Insulation improved the performance of both reactors by reducing heat losses resulting in a reduced air to feed ratio requirement. The hybrid gasifier gave a higher energy conversion efficiency, a higher product gas heating value, and a lower tar content than the open-core gasifier due to efficient gas mixing in a high temperature tar cracking region below the throat and reduced heat losses.

Numerical simulation of group combustion of pulverized coal

A mathematical model for the group combustion of pulverized coal particles was developed in a previous work. It includes the Lagrangian description of the dehumidification, devolatilization and char gasification reactions of the coal particles in the homogenized gaseous environment resulting from the three fuels, CO, H2 and volatiles, supplied by the gasification of the particles and their simultaneous group combustion by the gas phase oxidation reactions, which are considered to be very fast. This model is complemented here with an analysis of the particle dynamics, determined principally by the effects of aerodynamic drag and gravity, and its dispersion based on a stochastic model. It is also extended to include two other simpler models for the gasification of the particles: the first one for particles small enough to extinguish the surrounding diffusion flames, and a second one for particles with small ash content when the porous shell of ashes remaining after gasification of the char, non structurally stable, is disrupted. As an example of the applicability of the models, they are used in the numerical simulation of an experiment of a non-swirling pulverized coal jet with a nearly stagnant air at ambient temperature, with an initial region of interaction with a small annular methane flame. Computational algorithms for solving the different stages undergone by a coal particle during its combustion are proposed. For the partial differential equations modeling the gas phase, a second order finite element method combined with a semi-Lagrangian characteristics method are used. The results obtained with the three versions of the model are compared among them and show how the first of the simpler models fits better the experimental results.

Electronic structure methods applied to gas-carbon reactions

A review is given on the fundamental studies of gas-carbon reactions using electronic structure methods in the last several decades. The three types of electronic structure methods including semi-empirical, ab initio and density functional theory, methods are briefly introduced first, followed by the studies on carbon reactions with hydrogen and oxygen-containing gases (non-catalysed and catalysed). The problems yet to solve and possible promising directions are discussed. (C) 2002 Elsevier Science Ltd. All rights reserved.

Percolative fragmentation of char particles during gasification

Percolative fragmentation was confirmed to occur during gasification of three microporous coal chars. Indirect evidence obtained by the variation of electrical resistivity (ER) with conversion was supported by direct observation of numerous fragments during gasification. The resistivity increases slowly at low conversions and then sharply after a certain conversion value, which is a typical percolation phenomenon suggesting the occurrence of internal fragmentation at high conversion. Two percolation models are applied to interpret the experimental data and determine the percolation threshold. A percolation threshold of 0.02-0.07 was found, corresponding to a critical conversion of 92-96% for fragmentation. The electrical resistivity variation at high conversions is found to be very sensitive to diffusional effects during gasification. Partially burnt samples with a narrow initial particle size range were also observed microscopically, and found to yield a large number of small fragments even when the particles showed no disintegration and chemical control prevailed. It is proposed that this is due to the separation of isolated clusters from the particle surface. The particle size distribution of the fragments was essentially independent of the reaction conditions and the char type, and supported the prediction by percolation theory that the number fraction distribution varies linearly with mass in a log-log plot. The results imply that perimeter fragmentation would occur in practical combustion systems in which the reactions are strongly diffusion affected.

Variation of the crystalline structure of coal char during gasification

The variation of the crystallite structure of several coal chars during gasification in air and carbon dioxide was studied by high-resolution transmission electron microscopy (HRTEM) and X-ray diffraction (XRD) techniques. The XRD analysis of the partially gasified coal chars, based on two approaches, Scherrer's equation and Alexander and Sommer's method, shows a contradictory trend of the variation of the crystallite height with carbon conversion, despite giving a similar trend for the crystallite width change. The HRTEM fringe images of the partially gasified coal chars indicate that large and highly ordered crystallites exist at conversion levels as high as 86%. It is also demonstrated that the crystalline structure of chars can be very different although their pore structures are similar, suggesting a combination of crystalline structure analysis with pore structure analysis in studies of carbon gasification.

Thermogravimetric determination of the carbon dioxide reactivity of char from some New Zealand coals and its association with the inorganic geochemistry of the parent coal

Thermogravimetrically-determined carbon dioxide reactivities of chars formed from New Zealand coals, ranging in rank from lignite to high volatile bituminous, vary from 0.12 to 10.63 mg/h/mg on a dry, ash-free basis. The lowest rank subbituminous coal chars have similar reactivities to the lignite coal chars. Calcium content of the char shows the strongest correlation with reactivity, which increases as the calcium content increases. High calcium per se does not directly imply a high char reactivity. Organically-bound calcium catalyses the conversion of carbon to carbon monoxide in the presence of carbon dioxide, whereas calcium present as discrete minerals in the coal matrix, e.g., calcite, fails to significantly affect reactivity. Catalytic effects of magnesium, iron, sodium and phosphorous are not as obvious, but can be recognised for individual chars. The thermogravimetric technique provides a fast, reliable analysis that is able to distinguish char reactivity differences between coals, which may be due to any of the above effects. Published by Elsevier Science B.V.

Structural and surface property changes of macadamia nut-shell char upon activation and high temperature treatment

Structural and surface property changes of macadamia nut-shell (MNS) char upon activation and high temperature treatment (HTT) were studied by high-resolution nitrogen adsorption, diffuse reflectance infra-red Fourier transform spectroscopy, X-ray photoelectron spectroscopy, and temperature-programmed desorption. It is found that activation of MNS char can be divided into the low extent activation which may involve the reactions of internal oxygen-containing groups and leads to the formation of comparatively uniform micropores, and the high extent activation which induces reactions between carbon and activating gas and produces a large amount of micropores. The surface functional groups (SFGs) basically increase with the increase of activation extent, but high extent activation preferentially increases the amount of -C-O and -C=O. HTT in air for a short tithe at a high temperature (1173 K) greatly increases the micropore volume and the amounts of SFGs. By appropriately choosing the activation and HTT conditions, it is possible to control both the textural structure and the type and amounts of SFG. (C) 2002 Published by Elsevier Science Ltd.

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.

Sewage sludge gasification. Dolomite performance under different operating conditions

Gasification is a technology that can replace traditional management alternatives used up to date to deal with this waste (landfilling, composting and incineration) and which fulfils the social, environmental and legislative requirements. The main products of sewage sludge gasification are permanent gases (useful to generate energy or to be used as raw material in chemical synthesis processes), liquids (tars) and char. One of the main problems to be solved in gasification is tar production. Tars are organic impurities which can condense at relatively high temperatures making impossible to use the produced gases for most applications. This work deals with the effect of some primary tar removal processes (performed inside the gasifier) on sewage sludge gasification products. For this purpose, analysis of the gas composition, tar production, cold gas efficiency and carbon conversion were carried out. The tests were performed with air in a laboratory scale plant consisting mainly of a bubbling bed gasifier. No catalyzed and catalyzed (10% wt of dolomite in the bed and in the feeding) tests were carried out at different temperatures (750ºC, 800ºC and 850ºC) in order to know the effect of these parameters in the gasification products. As far as tars were concerned, qualitative and quantitative tar composition was determined. In all tests the Equivalence Ratio (ER) was kept at 0.3. Temperature is one of the most influential variables in sewage sludge gasification. Higher temperatures favoured hydrogen and CO production while CO2 content decreased, which might be partially explained by the effect of the cracking, Boudouard and CO2 reforming reactions. At 850ºC, cold gas efficiency and carbon conversion reached 49% and 76%, respectively. The presence of dolomite as catalyst increased the production of H2 reaching contents of 15.5% by volume at 850 °C. Similar behaviour was found for CO whereas CO2 and CnHm (light hydrocarbons) production decreased. In the presence of dolomite, a tar reduction of up to 51% was reached in comparison with no catalyzed tests, as well as improvements on cold gas efficiency and carbon conversion. Several assays were developed in order to test catalyst performance under more rough gasification conditions. For this purpose, the throughput value (TR), defined as kg sludge “as received” fed to the gasifier per hour and per m2 of cross sectional area of the gasifier, was modified. Specifically, the TR values used were 110 (reference value), 215 and 322 kg/h·m2. When TR increased, the H2, CO and CH4 production decreased while the CO2 and the CnHm production increased. Tar production increased drastically with TR during no catalysed tests what is related to the lower residence time of the gas inside the reactor. Nevertheless, even at TR=322 kg/h·m2, tar production decreased by nearly 50% with in-bed use of dolomite in comparison with no catalyzed assays under the same operating conditions. Regarding relative tar composition, there was an increase in benzene and naphthalene content when temperature increased while the content of the rest of compounds decreased. The dolomite seemed to be effective all over the range of molecular weight studied showing tar removal efficiencies between 35-55% in most cases. High values of the TR caused a significant increase in tar production but a slight effect on tar composition.

Characterization of sugar cane waste biomass derived chars from pressurized gasification

There is interest in the use of sugar cane waste biomass for electricity cogeneration, by integrated gasification combined cycle (IGCC) processes. This paper describes one aspect of an overall investigation into the reactivity of cane wastes under pressurized IGGC conditions, for input into process design. There is currently a gap in understanding the morphological transformations experienced by cane waste biomass undergoing conversion to char during pressurized gasification, which is addressed by this work. Char residuals remaining after pressurized pyrolysis and carbon dioxide gasification were analysed by optical microscope, nitrogen (BET) adsorption analysis, SEM/EDS, TEM/EDS and XPS techniques. The amorphous cane plant silica structures were found to remain physically intact during entrained flow gasification, but chemically altered in the presence of other inorganic species. The resulting crystalline silicates were mesoporous (with surface areas of the order of 20 m(2) g(-1)) and contributed to much of the otherwise limited pore volume present in the residual chars. Coke deposition and intimate blending of the carbonaceous and inorganic species was identified. Progressive sintering of the silicates appeared to trap coke deposits in the pore network. As a result ash residuals showed significant organic contents, even after extensive additional oxidation in air. The implications of the findings are that full conversion of cane trash materials under pressurized IGCC conditions may be significantly hampered by the silica structures inherent in these biomass materials and that further research of the contributing phenomena is recommended.