Abstract
The extensive use of traditional cooking and heating stoves to meet domestic requirements creates a serious problem of indoor and outdoor air pollution. This study reports the impacts of two fuel feeding methods – front-loading and top-loading on the thermal and emissions performance of a modern coal-fired water-heating and cooking stove using a contextual test sequence that replicates typical patterns of domestic use. Known as a low-pressure boiler, when this stove was fueled with raw coal, the findings indicate that front-loading the fuel, which devolatilizes the new fuel gradually, produced consistently higher space heating efficiency and lower emission factors than top-loading the same stove, which devolatilizes new fuel all at once. Comparing the performance at both high and low power gave the similar results: front-loading with raw coal produced consistently better results than top-loading. The average water heating efficiency when front-loading was (58.6±2.3 )% and (53.4±1.8 )% for top-loading. Over the sixteen-hour test sequence, front-loading produced 22% lower emissions of PM2.5 (3.9±0.6 ) mg/MJNET than top-loading (4.7±0.9 ) mg/MJNET. The same pattern was observed for carbon monoxide and the CO/CO2 ratio. CO was reduced from (5.0±0.4 ) g/MJNET to (4.1±0.5 ) g/MJNET. The combustion efficiency (CO/CO2 ratio) improved from (8.2±0.8 )% to (6.6±0.6 )%. Briquetted semi-coked coal briquettes are promoted as a raw coal substitute, and the tests were replicated using this fuel. Again, the same pattern of improved performance was observed. Front loading produced 3.5% higher heating efficiency, 10% lower CO and a 0.9% lower CO/CO2 ratio. It is concluded that, compared with top loading, the manufacturers recommended front-loading refueling behavior delivered better thermal, emissions and combustion performance under all test conditions with those two fuels.
1. Introduction
Coal is a major energy carrier in most countries. Important domestic uses of such energy are cooking and space heating[1,2]. It is estimated that some three billion people use solid fuels in the form of coal, crop residues, charcoal, wood, peat, and coke[3]. Chinese residential consumption was approximately 90 Mt of coal in 2011[4]. More than 50% of the energy in urban households and 22% in rural households was provided by coal burned in stoves and small low-pressure boilers (LPB)[5,6].
Consumption of coal and coal-based fuels for domestic cooking and space heating is a major source of greenhouse gases (GHGs) and atmospheric pollutants such as particular matter with a diameter of less than 2.5 micrometers (PM2.5), carbon monoxide (CO), nitrogen oxides (NOX), black carbon, and hydrogen sulfide (H2S). Some of those pollutants cause significant human health problems even in low concentrations. Household air pollution (HAP) is considered among the top environmental risks in China as is the case in many other developing countries[7]. According to the World Health Organization (WHO), an estimated 4.3 million people a year die prematurely from illness attributable to HAP created by the inefficient combustion of solid fuels[8]. It has been reported that low-quality coal is a significant cause of air pollution in China[9-12]. For example, raw coal combusted by rural households was reported to be responsible for 26% of urban and 32% of rural ambient PM2.5 during the winter heating season in Xi’an[13]. Those pollutants do not necessarily originate from inherent coal properties but are more often from the incomplete combustion of fuels in appliances not well-matched to the fuel properties. User behaviors can exacerbate under-performance arising from the mismatch between stoves and fuels.
There have been several projects to study and distribute modern stoves and fuels in rural China attempting to mitigate these pollution problems[14,15]. However, to be widely adopted, the features of a stove design must match the cultural and functional needs of the users. Thus, new standards have been implemented that consider thermal efficiency, cooking efficiency, quantification of emissions and a ten-hour fire endurance requirement[16,17]. Meeting these cultural demands and the new standards requires that operators understand the correct usage of stove features, including refueling techniques. Traditional domestic stoves in this sector have a hole on top that functions both as a cooking station and a refueling port. Some modern stoves that were presented on the market after 2000 feature a separate upper front door for refueling and increasingly common feature. Such a front-loading stove might (correctly) be refueled through the upper door on the front, or (incorrectly) through the top cooking hole – a well-known cultural behavior.
This study investigates the influence of these two alternative operators refueling behaviors, front loading and top loading, that might affect the emissions and thermal efficiency of such a modern domestic coal stove. The fire burns in a fundamentally different manner with these two methods. Top-loading places cold fuel on a hot coke base, leading to rapid devolatilization of the whole added batch, and tends to lead to anoxic, flameless gasification of at least part of the added fuel. Front-loading tends to feed the fire from one side which devolatilizes the new fuel gradually. This tends not to create a flameless condition, usually with better combustion conditions releasing a greater fraction of the embedded energy. The evaluation was performed using two common fuels, raw coal, and semi-coked coal briquettes. The assessments were performed using an online, real-time stove performance testing system developed at the China Agricultural University’s Biomass Energy Stove Testing (BEST) Laboratory. The stove selected was an 11.7 kW-rated coal-fired low-pressure boiler (LPB) with integrated cooking, suited to the requirements of typical rural households in China. The results of those tests can be used to inform policy and product selection for the residences in the agricultural sector of China.
2. Materials and Methods
2.1 Fuel preparation
Two fuels, raw coal passing through a 40 mm screen, and 50 mm diameter semi-coked coal briquettes, were used during the tests (Figure 1). The fuels were stored in an environment with stable humidity.
The ultimate and proximate analysis of the fuels was performed by the North China Electric Power University using relevant Chinese standards: (NY/T 1881.1-2010)[18] – Industrial Standard of the People’s Republic of China, (GB/T 28731-2012)[19] and (GB/T 28734-2012)[20] – National Standards of the People's Republic of China.
2.2 Testing procedures
The tests were performed at the BEST Laboratory at China Agricultural University (CAU), Beijing. The Hebei Farmer contextual testing sequence was used. It was derived from field observations and used in the Hebei Clean Heating Project in 2017[21]. The testing sequence is sixteen hours long and represents the typical cooking and heating behaviors of rural farmers in Hebei Province, China (Figure 2).
Figure 1 Raw coal (left) and semi-coked coal briquettes (right)
Figure 2 16-hour contextual test sequence
The poorest combustion conditions are most common when igniting, refueling or when ramping from low power heating to high power cooking. Banking the fire from high to low power sometimes produces high emissions. Increasing the number of power changes or the number of cooking events will affect the emitted mass. To establish “comparative performance”, a full range of typical behaviors should be replicated during a standard test sequence. If the ignition is once per heating season, the ignition emissions should be ignored. The performance parameters reported are the thermal efficiency, combustion efficiency and emission factors for PM2.5 and CO expressed in unit mass per Megajoule of energy delivered into the circulating water. Measurements were automatically logged at least every 10 s and stored for analysis. The schematic diagram of the stove used for two fuel feeding option shown in Figure 3.
Figure 3 Stove with two fuel feeding options
2.3 Test apparatus
Temperatures were measured with K-type thermocouples logged once per second. An ultrasonic flowmeter was used to determine the water flow rate through the heat exchanger [L/s]. A gas analyzer (MRU MGA5/Vario Plus, Germany) was used to measure flue gas concentrations. It has a non-dispersive infrared (NDIR) sensor for measuring carbon dioxide (CO2) and electrochemical cells for trace gases: O2, CO, NOx, and SO2. The particulate matter was measured in a diluted air stream using a DustTrak DRX aerosol monitor (Model 8533) which consecutively reports PM1, PM2.5, PM4, and PM10. The dilution system was fashioned after that of the SeTAR Centre (University of Johannesburg, South Africa)[22]. A sample drawn continuously from the flue is diluted by CO- and CO2-free dry air. The manually variable dilution ratio is provided by comparing simultaneous measurements of CO2 in the flue gas and the diluted sample stream.
The stove testing system is shown in Figure 4. Before each test, the function of the instruments was checked. A zero calibration of the Dusttrak particle counter was performed before each test. The gas analyzers were zero calibrated before every test. A total of twelve tests were conducted - three replications of each combination of fuel and loading behavior.
The level of dilution by the particle diluter can be adjusted at the control panel by altering the flow rate of the dilutant gas, permitting the measurement of a wide range of smoke concentrations. The dilution ratio can be varied from 3:1 to 100:1 as required, to keep the concentration within the operating limits of the particle counter (0.001-150 mg/m3). The dilutant is air scrubbed in an absorber to produce the zero air containing low concentrations of CO (< 1 ppm), CO2 (< 1 ppm) and water vapor (< 1 ppm). The dilutant dehydrates the particles and enhances the condensation of particles from the gaseous phase by cooling. The integrated mass from the Dusttrack continuous recordings is compared with the mass collected on a 37 mm filter attached to the outlet to derive the particle density factor for the particles (conversion from particle number to mass). The particle detector has an optical measurement system. Above 50% humidity, the presence of water vapor affects the optical diameter of smoke particles, so it is necessary to monitor the relative humidity. The dilution ratio is adjusted during a test to maintain the relative humidity within the range of 15%-30%. The dilution level is thus determined by two considerations: keeping the maximum PM concentration under the 150 mg/m3 limit or optimizing the detection limit at low concentrations, and simultaneously keeping the relative humidity under 30%.
Figure 4 Online testing system
The dilution ratio for each interval is required for the calculation of the emitted mass of PM. The CO2 in the dilution tunnel is compared with the flue gas concentration and the dilution level obtained per recording interval.
At the start of a test, the stove was placed on the electronic mass balance, together with a pre-weighed load of kindling and fuel. The fire was ignited with paper and small pieces of coal from the pre-weighed load. Measurements were recorded every 10 s. Parameters recorded included mass of the stove system, space heating heat exchanger inlet and outlet water temperatures, temperature of the cooking pot contents, inlet and outlet temperatures of the cooking pot heat exchanger, flow rates of the coolant through the heat exchangers, all PM channels and trace gas emissions from the flue and from the diluter, ambient and diluted gas humidity, and the ambient air temperature and pressure. During the high-power portion of the test sequence, the flow of coolant through the space heating heat exchanger was adjusted as needed to keep the outlet temperature under 70°C . The concentrations of trace gases H2, CxHy, SO2, NO, CO, CO2, O2 were measured directly from flue while in the diluter measurements included CO2, H2O, H2S. The calculated outputs include the
2.4 Calculated values
The performance parameters were calculated using the equations listed below.
The heat gained by the LPB heat exchanger (QHEx) and the water contained within it is calculated as:
where, CpHEx is the specific heat capacity of the heat exchanger; CpW is the specific heat capacity of water, J/g; mHEx is the mass of the heat exchanger, g; mWHEx is the mass of the water in the heat
exchanger, g; ΔTHEx is the change in average the temperature of the water in the heat exchanger, °C, which is calculated using Equation (3):
where, the item with subscript f is the final time; the item with subscript i is the initial time; T1 is inlet temperature, °C; T2 is the outlet temperature, °C; Heat gained by the water flowing through the heat exchanger per recording interval (QjWHEx) is calculated as:
where, (t2–t1) is the recording interval, s; FRWHEx is the water flow rate through the heat exchanger, g/s.
The total heat gained by the water flowing through the heat exchanger during the test (QWHEx) is calculated as:
where, Fr is the flow rate of water through cooking pot heat exchanger g/s; n and j–1 from measurement 1 to measurement n is the total heat gained.
The total heat gained by the heat exchanger system (QTH) is calculated as:
The heat gained by the cooking pot heat exchanger and the water contained within (Qc) it is calculated as:
where, CpCEx is the specific heat capacity of the cooking pot heat exchanger; mCEx is the mass of the pot heat exchanger, g; mWCEx is the mass of the water in the pot heat exchanger, g; ΔTCEx is the change in the average temperature of the water in the pot heat exchanger, °C, calculated using Equation (8):
where, the item with subscript f is the final time; the item with subscript i is the initial time.
Heat gained by the water flowing through the cooking pot heat exchanger per recording interval (QjWCEx) is calculated as:
where, T3 is the cooking pot heat exchanger water inlet temperature; T4 is the cooking pot heat exchanger water outlet temperature; FRWCEx is the water flow rate through the cooking pot heat
exchanger, g/s.
The total heat gained by the water flowing through the cooking pot heat exchanger (QWHEx) is calculated as:
The heat gained by the cooking pot and its uncirculated contents per recording interval (QjPot) is calculated using Equation (11):
where, CpPot is the specific heat capacity of the cooking pot, J/g; mPot is the mass of the pot, g.
The total heat gained by the cooking pot and its uncirculated contents (Qp) is calculated as:
where, j is heat gained by uncirculated water in cooking pot J/g. The total energy gained by the cooking system (QTC) is calculated as:
The total energy gained by the heating and cooking system (QT) is calculated as:
The total chemical energy available (UT) from the fuel fed into the stove during the test (Equation (15)):
where, BK is the mass of kindling materials, kg; LHVK is the lower heating value of the kindling materials as received, MJ/kg; BF is the mass of fuel, kg; LHVF is the LHV of the fuel as received, MJ/kg. The fractional energy utilization is calculated as the ratio of UT – the useful heat delivered by the stove - to QT - the total heat available from the fuel fed. The pre-weighed fuel was placed on the scale before starting the test. The test continued for sixteen hours from ignition to completion. The one-hour cooking phase commenced twelve hours after ignition. The cooking efficiency is based on readings from that period. Changes in mass over selected intervals were determined from the mass balance data and were interpreted to be the mass of fuel consumed during the respective space heating or cooking phases.
The cooking efficiency (ηc) was calculated as:
The space heating efficiency (ηh) is calculated as:
The system (overall) (ηt) efficiency was calculated using Equation (18):
The total volume of diluted gases emitted (VDil~1) is calculated as:
where, VSto is the stoichiometric volume of gases produced by the combustion of 1 kg of fuel (AR), m3; mfuel is the mass of fuel fed, kg; λ is the total air demand factor providing both used and unused
oxygen.
The dilution ratio (DR) between the flue gases and the diluted gas sample is calculated as:
The total mass of PM reported during the sampling interval (TPM2.5) is calculated using Equation (21):
where, PM2.5 is the total measured mass of particles reported by the instrument.
The PM2.5 emission factor (EFPM₂.₅) is calculated as:
where, QT is the total useful heat energy, MJ, from Equation (14).
The average mass concentration of CO is calculated as:
where, ppm(v) is the average CO concentration expressed on a volumetric basis, ppm; and 0.0012334 is the conversion constant for CO ppm(v) to CO(m), g/m3.
The total mass of CO emitted is calculated as:
The CO emission factor (EFCO) is calculated as:
The percentage increase in thermal efficiency (ηTI) is calculated as:
where, ηTL is the thermal efficiency of top-load, %; ηFL is the thermal efficiency of top-load %.
The percentage decrease in CO, CO/CO2, and PM2.5 is calculated using Equations (27)-(29), respectively.
where, COD is the decrease in CO, %; COTL is the CO emission of top-load g/MJNET; COFL is the CO emission of front-load g/MJNET .
where, (CO/CO2)D is the decrease in CO/CO2 ratio, %; CO/CO2TL is the CO/CO2 emission of top-load %; CO/CO2FL is the CO/CO2 emission of front-load %.
where, PM2.5D is the decrease in PM2.5 emission %; PM2.5TL is the PM2.5 emission from top-load mg/MJNET; PM2.5FL is the PM2.5 emission from front-load mg/MJNET.
2.5 Equipment Specifications
A catalogue of testing indicators, instruments and testing principles is given in Table 1. Sixteen-hour duration tests were conducted to evaluate the performance of the stoves, combining the relevant Chinese standard methods with modifications from observed local behaviors. The testing sequence comprises various stages informed by observing the cooking and heating practices of rural Chinese users.
cooking and space heating efficiency, several emission factors and the combined cooking and heating efficiency according to Beijing Municipal Standard (DB11/T540–2008) General Technical Specification of Domestic Biomass Stove/Boiler (Beijing Municipal Bureau of Quality and Technical Supervision, 2008)[16].
The fire was set and ignited according to the manufacturer’s instructions. The quantity of fuel adequate for the test is estimated using the following equation:
where, Mf is the mass of the fuel set aside, kg; 3.6 converts J/s to kJ/h; Dt is the testing duration, h; PrNET is rated heating power, kWNET; LHVf is the lower heating value of the fuel as received (AR), MJ/kg; ηr is the rated high power thermal efficiency, %; 1.1 was added as a safety margin so as not to run out before the test ended.
2.4 Calculated values
The performance parameters were calculated using the equations listed below.
The heat gained by the LPB heat exchanger (QHEx) and the water contained within it is calculated as:
exchanger, g; ΔTHEx is the change in average the temperature of the water in the heat exchanger, °C, which is calculated using Equation (3):
The total heat gained by the water flowing through the heat exchanger during the test (QWHEx) is calculated as:
The total heat gained by the heat exchanger system (QTH) is calculated as:
Heat gained by the water flowing through the cooking pot heat exchanger per recording interval (QjWCEx) is calculated as:
exchanger, g/s.
The total heat gained by the water flowing through the cooking pot heat exchanger (QWHEx) is calculated as:
The total heat gained by the cooking pot and its uncirculated contents (Qp) is calculated as:
The cooking efficiency (ηc) was calculated as:
oxygen.
The dilution ratio (DR) between the flue gases and the diluted gas sample is calculated as:
The PM2.5 emission factor (EFPM₂.₅) is calculated as:
The average mass concentration of CO is calculated as:
The total mass of CO emitted is calculated as:
The percentage decrease in CO, CO/CO2, and PM2.5 is calculated using Equations (27)-(29), respectively.
2.5 Equipment Specifications
A catalogue of testing indicators, instruments and testing principles is given in Table 1. Sixteen-hour duration tests were conducted to evaluate the performance of the stoves, combining the relevant Chinese standard methods with modifications from observed local behaviors. The testing sequence comprises various stages informed by observing the cooking and heating practices of rural Chinese users.
Table 1 Description of test indicators, equipments, and testing principles
3. Results and Discussion
3.1 Ultimate analysis of the fuels
Both fuels were analyzed for their elemental composition, moisture content, volatile matter fraction and higher heating value (HHV) from which the lower heating value (LHV) was calculated. The moisture content of the raw coal (9.04%) was higher than that of the semi-coked coal briquettes (3.5%), determined on a wet weight basis. The volatiles content of the briquettes was low (Table 2). The nitrogen content of the briquettes was higher than that of the raw coal. The efficiency of the stove when top-loaded with raw coal (54.7%±1.8%) was significantly lower than when top-loaded with briquettes (62.7%±2.7%) . This demonstrates that the fuel can strongly affect stove performance. Fuels analyses are provided in Table 2.
Table 2 Proximate analysis as received and ultimate analysis after drying
Table 3 Increase in thermal efficiency and decrease in CO, CO/CO2, and PM2.5 of front load method as compared to top load (%)
3.2 Thermal performance
The test results showed that front-loading the fuel produced a higher thermal efficiency in all cases, even when combustion efficiency was good. Thermal efficiencies are presented in Figure 5. The front-loading thermal efficiency burning raw coal was (54.7±1.8 )% and (62.7±2.7 )% for briquettes. The top-loading thermal efficiency burning raw coal was (49.2±1.8 )% and (58.1±2.1 )% for briquettes. Regarding fuel savings, front-loading the stove saved (7.3±3.4 )% for briquettes and (10.1±2.5 )% for raw coal. Increase in thermal efficiency of front load method as compared to top load are shown in Table 3.
Figure 5 Thermal performance of the front-loading versus top-loading for low and high power and two types of fuel
3.3 Emission measurement
3.3.1 Particular matter (PM2.5)
The average PM2.5 emitted when top-loaded was 4.6±0.8 mg/MJNET burning raw coal while when burning briquettes, it was 2.9±0.7 mg/MJNET. Front-loading was better for burning both raw coal (2.8±0.6 mg/MJNET) and briquettes (2.5±0.7 mg/MJNET) (Figure 6). Analysis of the real-time testing data that the PM2.5 emission numbers were higher for top-loading than front-loading because the fuel dropped directly on fire and disturbed fuel already in the combustion zone.
In a previous study, Li et al.[23] found PM2.5 was 53± 6 mg/MJNET when top-loading raw coal and 23±5 mg/MJNET top-loading briquettes in a modern stove, a reduction of 57%. This study showed a reduction of (26.2±6.8 )% achieved at high power when substituting semi-coked coal briquettes for raw coal. The decrease in PM2.5 of front load method as compared to top load are shown in Table 3.
Figure 6 PM2.5 mass of per Megajoule delivered, front-loaded versus top-loaded for the two fuels at low and high power
3.3.2 CO/CO2 ratio
The lower the CO/CO2 ratio, the more complete the combustion of carbon. The CO/CO2 ratios during high power for front-loaded raw coal was lower than for top-loaded raw coal: (5.7±0.7 )% versus (7.6±0.8 )% (Figure 7). The CO/CO2 ratios for both feeding methods were found to be lower during high power combustion. The CO/CO2 ratio for top-loading was lower during high power (7.63±0.81 )% and (6.80±0.72 )%, and higher during low power, (8.7±0.9 )% and (7.2±0.8 )% burning raw coal and briquettes respectively. The CO/CO2 ratios for front-loading was lower than when top-loaded for both fuels and all operating conditions. The CO/CO2 ratios for front-loading at high power for raw coal and briquettes was, respectively, (5.7±0.6 )% and (5.3±0.5 )% while at low power it was (7.6±0.6 )% and (6.7±0.6 )%. It was observed that for both refueling behaviors the stove operating at high power had relatively good CO/CO2 ratios, but the combustion quality declined slightly as the burn rate decreased due to an increasing level of excess air. The results for decrease in CO/CO2 of front load method as compared to top load are shown in Table 3.
Figure 7 CO/CO2 ratios of front-loaded versus top-loaded for low and high power and two fuels
The results supported previous findings that the CO/CO2 ratio is affected by several factors including air-fuel ratio, burn rate, combustion temperature, combustion efficiency, thermal efficiency, residence time in the combustion chamber and flame turbulence[24-27].
Makonese[28] observed that the low combustion efficiency during ignition and in the smouldering phase, the CO/CO2 ratio was poor at 15%, indicating a need for further design changes to improve the combustion efficiency.
3.3.3 CO emission factors
The CO emission factor for front-loading found to be lower than the top-loading. The minimum emission factors were 2.0±0.2 g/MJNET at low power and 2.9±0.2 g/MJNET at high power when front-loading the briquettes. When top-loaded, this increased to 3.6±0.4 g/MJNET at high power and 2.4±0.3 g/MJNET at low power. Both loading methods produced a higher CO emission factor burning raw coal, but still, front-loading was significantly lower than the top-loading as shown in Figure 8. The decreases in CO of front load method as compared to top load are shown in Table 3.
Figure 8 CO emission factors for front-loading versus top-loading refueling for low and high power and two fuel types
4. Conclusions
This study investigated the effects of fuel feeding methods on thermal efficiencies and emissions of a modern cooking and heating stove burning semi-coked coal briquettes and raw coal. Sixteen-hour contextual tests representing typical behavior patterns were conducted to evaluate stove performance using a real-time detection system. The results showed that when the fuel was fed into the front of the stove, it invariably produced better in performance than when top-loaded. When burning semi-coked coal briquettes, the front-loading method delivered a higher average thermal efficiency (66.7±3.8 )% and the lowest emissions of PM2.5 2.9±1.0 mg/MJNET and CO 2.4±0.3 g/MJNET.
On the other hand, it was observed that fuel type also influenced the performance. It is concluded that 50 mm semi-coked coal briquettes gave better performance than randomly sized raw coal. The front-loading method delivered a higher thermal efficiency, saving fuel, and contributing to a cleaner environment.
Additional experiments were conducted using a similar stove designed for top loading only. At all power levels and for both fuels, all metrics show front loading to be superior refueling behavior. We conclude that to maximize the benefits of modern stoves, user training is essential because something as elementary as how the stove is refueled can make a considerable difference to the efficiency and emissions performance.
Source: Riaz Ahmad1,2,3, Yuguang Zhou1,2,3*, Nan Zhao1,2,3, Crispin Pemberton-Pigott3,4,
Harold John Annegarn3,4, Muhammad Sultan5, Renjie Dong1,2,3,6, Xinxin Ju7
1. Bioenergy and Environmental Science and Technology Laboratory, College of Engineering, China Agricultural University
2.Key Laboratory of Clean Renewable Energy Utilization Technology, Ministry of Agriculture and Rural Affairs, China;
3.National Center for International Research of BioEnergy Science and Technology, Ministry of Science and Technology, Beijing
4.School of Geo- and Spatial Sciences, North-West University, South Africa;
5.Department of Agricultural Engineering, Bahauddin Zakariya University
6.Yantai Institute, China Agricultural University
7.Shandong Sino-March Environmental Technology Co., Ltd.
The 10 largest coal producers and exporters in Indonesia:Harold John Annegarn3,4, Muhammad Sultan5, Renjie Dong1,2,3,6, Xinxin Ju7
1. Bioenergy and Environmental Science and Technology Laboratory, College of Engineering, China Agricultural University
2.Key Laboratory of Clean Renewable Energy Utilization Technology, Ministry of Agriculture and Rural Affairs, China;
3.National Center for International Research of BioEnergy Science and Technology, Ministry of Science and Technology, Beijing
4.School of Geo- and Spatial Sciences, North-West University, South Africa;
5.Department of Agricultural Engineering, Bahauddin Zakariya University
6.Yantai Institute, China Agricultural University
7.Shandong Sino-March Environmental Technology Co., Ltd.
- Bumi Resouces (BUMI)
- Adaro Energy (ADRO)
- Indo Tambangraya Megah (ITMG)
- Bukit Asam (PTBA)
- Baramulti Sukses Sarana (BSSR)
- Harum Energy (HRUM)
- Mitrabara Adiperdana (MBAP)
- Samindo Resources (MYOH)
- United Tractors (UNTR)
- Berau Coal
