Citation: Nour Khlaifat, Ali Altaee, John Zhou, Yuhan Huang. A review of the key sensitive parameters on the aerodynamic performance of a horizontal wind turbine using Computational Fluid Dynamics modelling[J]. AIMS Energy, 2020, 8(3): 493-524. doi: 10.3934/energy.2020.3.493
[1] | M Toufiq Reza, Jessica Nover, Benjamin Wirth, Charles J Coronella . Correction: Hydrothermal carbonization of glucose in saline solution: sequestration of nutrients on carbonaceous materials. AIMS Energy, 2018, 6(2): 269-271. doi: 10.3934/energy.2018.2.269 |
[2] | Diakaridia Sangare, Ayoub Missaoui, Stéphane Bostyn, Verónica Belandria, Mario Moscosa-Santillan, Iskender Gökalp . Modeling of Agave Salmiana bagasse conversion by hydrothermal carbonization (HTC) for solid fuel combustion using surface response methodology. AIMS Energy, 2020, 8(4): 538-562. doi: 10.3934/energy.2020.4.538 |
[3] | Vo Thanh Phuoc, Kunio Yoshikawa . Comparison between direct transesterification of microalgae and hydrochar. AIMS Energy, 2017, 5(4): 652-666. doi: 10.3934/energy.2017.4.652 |
[4] | Vo Thanh Phuoc, Kunio Yoshikawa . Effect of the storage condition of microalgae on hydrochar lipids and direct esterification-transesterification of hydrochar lipids for biodiesel production. AIMS Energy, 2017, 5(1): 39-53. doi: 10.3934/energy.2017.1.39 |
[5] | Thananat Lungkadee, Tossapon Katongtung, Pich-ramon Pokkanta, Tossaporn Onsree, Chawannat Jaroenkhasemmeesuk, Nakorn Tippayawong . Current technology landscape of biochar in carbon capture and storage research via bibliometric analysis. AIMS Energy, 2024, 12(1): 277-303. doi: 10.3934/energy.2024014 |
[6] | Arief Heru Kuncoro, Afri Dwijatmiko, Noer’aida, Vetri Nurliyanti, Agus Sugiyono, Widhiatmaka, Andri Subandriya, Nurry Widya Hesty, Cuk Supriyadi Ali Nandar, Irhan Febijanto, La Ode Muhammad Abdul Wahid, Paul Butarbutar . Mapping the landscape of carbon trading & carbon offset research: A global and Indonesian perspective. AIMS Energy, 2025, 13(1): 86-124. doi: 10.3934/energy.2025004 |
[7] | Ankita Juneja, Ganti S. Murthy . Evaluating the potential of renewable diesel production from algae cultured on wastewater: techno-economic analysis and life cycle assessment. AIMS Energy, 2017, 5(2): 239-257. doi: 10.3934/energy.2017.2.239 |
[8] | Vui Thi Hoang, Diem Huong Hoang, Ngoc Duc Pham, Hanh My Tran, Ha Thi Viet Bui, Tien Anh Ngo . Hydrogen production by newly isolated Clostridium species from cow rumen in pure- and co-cultures on a broad range of carbon sources. AIMS Energy, 2018, 6(5): 846-865. doi: 10.3934/energy.2018.5.846 |
[9] | Nirwan Syarif, Dedi Rohendi, Ade Dwi Nanda, M. Try Sandi, Delima Sukma Wati Br Sihombing . Gas diffusion layer from Binchotan carbon and its electrochemical properties for supporting electrocatalyst in fuel cell. AIMS Energy, 2022, 10(2): 292-305. doi: 10.3934/energy.2022016 |
[10] | Eric Danso-Boateng, Osei-Wusu Achaw . Bioenergy and biofuel production from biomass using thermochemical conversions technologies—a review. AIMS Energy, 2022, 10(4): 585-647. doi: 10.3934/energy.2022030 |
Hydrothermal carbonization (HTC) has gained attentions from various scientific fields from solid fuel to carbon based catalyst and from waste treatment to low cost adsorbents over the last several years [1,2,3,4,5,6,7,8,9]. Because of its very high ionic product [10], when biomass is treated with subcritical water (200-280 °C), organic components of the biomass degrade and re-polymerize into solid hydrochar. Hydrochar is hydrophobic, energy dense, carbon-rich amorphous solid, although the quality depends on HTC reaction parameters as well as raw feedstock [11,12,13]. In addition to fuel application, hydrochar’s potential for soil amendment, slow release fertilizer, and carbon sink have been evaluated by several research groups around the world [8,14].
Several publications have already reported the fate of elemental nitrogen, phosphorus, and potassium during HTC from various feedstocks including lignocellulosic or waste biomass [2,15,16]. In fact, a recent article discussed about the fate of plant available nutrients in the form of NH4+ and PO43− during HTC of digestate [14]. During HTC, nitrogen could be recovered either as HTC process liquor or distributed evenly in hydrochar and liquid phase [16,17]. The difference in feedstock as well as the different reaction condition makes the comparison unrealistic. Phosphorus has been recovered primarily in hydrochar, although there are several cases, where majority of the phosphorus was found in liquid phase [2,15]. In contrary, most literatures agree with the fate of potassium in the process liquid [2,14,16]. It is impossible to discuss the opposing reports because different feedstocks, experimental setups, and reaction conditions have been applied. However, all of these literatures failed to explain the nutrient recovery phenomena.
The HTC reaction chemistry has been evaluated very extensively studying both model compounds and various biomass feedstocks [12,18,19,20]. There are evidences of chemisorption sites on the hydrochar [21], which makes hydrochar effective for chemisorption of selective substances like U(VI), Cu, Cd, triclosan etc. [22,23,24,25]. However, the sorption phenomena for plant available nutrients (N, P, K) need further investigation. It is still unknown whether the nutrients in the solid hydrochar result from nutrient adsorption or nutrient remaining due to the incomplete degradation of fiber components.
The main goal of this study was to find out whether nutrients are sequestered during HTC. The term sequestration can be defined for this study as the nutrients uptake on the hydrochar during HTC. To eliminate nutrients sorption from the feedstock nutrients, glucose solution has been used as feedstock. Various salts containing targeted nutrients (NPK) are used as nutrient sources. Investigations of specific nutrient contents as well as their positions on hydrochar are explained here. Moreover, a simple reaction mechanism is proposed and discussed in this article.
All the salts (potassium chloride, ammonium hydrophosphate, potassium sulfate, and anhydrous ferric chloride) and sugar (glucose) were purchased from VWR, Germany. The salts and sugar were more than 98% pure. Rest 2% of the salts is impurity including copper, chromium, calcium, aluminum, and cobalt. Individual 1 M solutions of all sugar and salts are produced using de-ionized water and are stored in a sealed glass bottle until HTC.
In a HTC experiment, 60 mL of 1M glucose solution, 10 mL of 1M KCl, 10 mL of 1M NH4H2PO4, 10 mL of 1M K2SO4, and 10 mL of 1M FeCl3.6H2O were poured in a 200 mL Parr pressure bomb (Parr model 4570, Moline, IL). The experimental condition was set to 180,230 or 300 °C for holding time 16 h with 3 K min−1 heating rate. According to the previous studies, HTC reactions, especially in the liquid state, requires several hours to complete depending on the reaction parameters [8,19,31]. To ensure the completion of the HTC reaction, a very long reaction time (16 h) was chosen for this study. After the end of reaction period, the heater was turned off and let the reactor cooled down naturally. It took 3-4 h to cool down from 300 °C to 25 °C (about 30-45 min from 300 °C to 180 °C). The gaseous product was purged in the fume hood; hydrochar was filtered using a folded filter paper (ROTH Type 113 P filter) for 20 min and washed with excess de-ionized water for another 10 min to remove physically adhered nutrients on the surface. Hydrochar was dried in a heating oven at 105 °C for overnight. Dried hydrochar was placed into a zip-lock bag and stored for further physico-chemical analyses. All the experiments were at least duplicated for this study. Hydrochars were named as HTC-T, where HTC stands for hydrothermally carbonized char and T is the HTC temperature, respectively.
Elemental carbon, hydrogen, nitrogen, and sulfur (CHNS) were measured using a Vario El elemental analyzer (Elementar Analysesysteme Hanau, Germany). Sulphonic acid was used in this elemental analysis as reference and two ovens were set at 1150 °C, and 850 °C, respectively. Each sample was analyzed three times, and average elemental compositions are reported.
A Nicolet iS5 ATR-FTIR (Thermo Scientific, Madison, WI, USA) with mid- and far-IR capabilities was used on the raw and hydrothermally treated biomass. IR spectra of all solid samples were recorded at 30 °C using ATR-FTIR. All samples were milled into fine powder for homogenize and dried in 105 °C for 24 h in an oven prior to FTIR. Only 5-10 mg of dry sample was placed in the FTIR for this analysis and pressed against the instrument's diamond surface with its metal rod. All spectra were obtained using 64 scans for the background (air) and 64 scans for the samples, which were scanned from 500-4000 cm−1.
The morphology of samples HTC-180, HTC-230, and HTC-300 were analyzed using SEM-EDX and nitrogen adsorption. For SEM, hydrochar samples were dispersed on a conductive pad and introduced into the SEM equipment. The investigation was carried out using a Cambridge S200 scanning electron microscope (SEM) at an excitation voltage of 20kV and a sample distance of 18mm. The SEM samples were prepared by depositing about 50 mg of dry powdered hydrochar on an aluminum stud covered with carbon conduction tapes, and then coated with Pt for 1 min to prevent charging during observation.
The specific surface area was determined based on nitrogen adsorption at 77 K and the BET method according to ISO 9277. Barrett-Joyner-Halenda (BJH), Dubinin-Radushkevich (DR), Horvath-Kawazoe (HK), and Non-Local Density Functional Theory (NLDFT) methods were used to determine the volume of meso- and micropores according to DIN 66134 and DIN 66135, respectively. Hitachi scanning electron microscope (SEM) S-2700 was used to investigate the influence of the process setting on hydrochars’ microscopic structure.
Targeted nutrients (potassium and phosphorus) and trace elements (aluminum, calcium, cobalt, chromium, copper, iron, magnesium, manganese, sodium, nickel, sulfur, and zinc) were determined by an induced coupled plasma optical emission spectrometry (ICP-OES). The hydrochar samples were digested by nitric acid (ACS grade, 65% pure) and hydrogen peroxide (10% pure) at 250 °C for 20 min in a MWS 640 microwave (MLS GmbH, Germany). For ICP-OES, an iCAP 6000 analyzer from Thermo Fisher Scientific, Inc. USA was used for inorganic determination.
A series of HTC experiments were conducted on glucose solution in presence of various salts and varying HTC temperature. The salt and sugar concentration in the feed was maintained similar to find out the HTC temperature effect on nutrient sequestration. As HTC of glucose as model compound was reported from 180-250°C [24,26,27,28,29], two different HTC temperatures (180 and 230 °C) and one higher temperature (300 °C) were chosen for this study. The long reaction time (16h) was chosen to ensure the completion of HTC reactions even at the lowest temperature (180 °C). It was previously reported that the presence of iron catalyze in favor of hydrochar during HTC process[21,30]. Moreover, iron itself is an essential micronutrient for plant growth. Hence, iron salt was added with others for nutrient sequestration in hydrochar. As nutrient sequestration is the main objective of this study, the process liquid and the gaseous product were not further studied. Solid production (weight of solid hydrochar product from 100 mL of feedstock), ultimate analysis (CHNS), and selective nutrient elements are reported in Table 1.
Hydrochar | Solid production1 (g per 100 mL) | Mass yield (g g-1dry reactants) | Elemental analysis2 | ICP-OES3 | |||||||
Nitrogen (%) | Carbon (%) | Sulfur (%) | Hydrogen (%) | Oxygen (%)4 | Iron (mg kg-1) | Potassium (mg kg-1) | Phosphorus (mg kg-1) | Sulfur (mg kg-1) | |||
Dry reactants | NA | 1.0 | 1.4 | 23.3 | 2.3 | 4.2 | 68.8 | 55,800 | 117,194 | 31,000 | 23,593 |
HTC-180 | 9.1 ± 0.5 | 0.53 ± 0.04 | 0.9 | 62.7 | 0.2 | 4.5 | 31.7 | 5271 | 3800 | 3424 | 533 |
HTC-230 | 9.5 ± 0.4 | 0.56 ± 0.04 | 1.0 | 64.0 | 0.2 | 4.5 | 30.3 | 8920 | 6389 | 5346 | 790 |
HTC-300 | 11.1 ± 0.5 | 0.65 ± 0.06 | 0.9 | 63.1 | 0.3 | 3.9 | 32.4 | 15,621 | 6309 | 6412 | 1579 |
1 @Solid \; production = \frac{{mass \; of \; hydrochar \; dry \; \left( g \right)}}{{Total \; volume \; of \; feed \; \left( {ml} \right)}} \times 100@ | |||||||||||
2 Precision of the elemental analysis is ±0.5 wt% | |||||||||||
3 Precision of the ICP-OES is 100 mg kg−1 | |||||||||||
4 % O = 100 - (%C + %H + %N + %S); NA: Not applicable |
Nutrient sequestration on the hydrochar during HTC can be seen from Table 1. All the nutrients (N, K, and P) at different concentrations have been found in the hydrochar. As the starting HTC feedstock was liquid solution of sugar and salts, formation of hydrochar during HTC with nutrients proves that nutrients from the liquid streams can be sequestered in the solid hydrochar during HTC. It can also be found from Table 1 that solid recovery was increased with HTC temperature. According to the previous studies [26,28,31,32], HTC temperature above 180 °C has minimal effects on glucose-derived-hydrochar production for an extended HTC reaction time. In other words, carbon in hydrochar, once produced from glucose, is fairly stable under subcritical water. In fact, the elemental carbon content from Table 1 implies the consistent result in this study. Now, an increase of mass yield indicates either an increase of nutrients in hydrochar with the increase of HTC temperature, or the salts catalyze HTC reaction and catalytic activity increases with HTC temperature, or both.
If the mass yield is expressed as mass of hydrochar produced from mass of salt and sugar input, then the mass yield become 53.3, 55.6, and 64.5% in HTC-180, HTC-230, and HTC-300, respectively. Moreover, the carbon yield (ratio of carbon (g) in product and carbon(g) in feedstock) is very similar to the overall mass yield. In fact, it was reported earlier that only mass yield or solid production is around 50-53% of the overall mass when glucose was hydrothermally carbonized from 180-250 °C for 1-16 h [33,39]. Unlike the previous study, the mass yield increases with HTC temperature. The carbon concentration remains similar, which means that the overall mass of the carbon along with the nutrient concentrations increase in hydrochar with HTC temperature.
Unlike carbon, nitrogen concentration in the hydrochar was increasing with HTC temperature as well. In terms of mass recovery of nitrogen, around 61.5-71 wt% of nitrogen was sequestered on hydrochar and mass of nitrogen in the hydrochar increases with HTC temperature. Other nutrients e.g., iron, potassium, sulfur, and phosphorus show an increasing trend of their concentration on hydrochar with the increase of HTC temperature. Similar increase of nutrients can be observed from previous study [40]. The nutrients did not follow a linear relationship; however, overall they increase with the increase of HTC temperature from 180 °C to 230 °C. The increment of temperature above
230 °C further sequester several nutrients, while the others like potassium remains similar. Potassium may have reached the maximum sequestration at around 230 °C. Among these nutrients, iron and sulfur concentration increases about two times from HTC-180 to HTC-300. In terms of mass recovery, mass of iron sequestration increases from 8.6 wt% to 31.5 wt% in HTC-180 to HTC-300 hydrochar. Phosphorus concentration also increases about 87 % from HTC-180 to HTC-300. Phosphorus sequestration on hydrochar increases from 10.1 wt% to 23 wt% in HTC-180 to HTC-300. Unlike nitrogen, iron, and phosphorus, potassium sequestration in hydrochar was very low. In fact, only 3.0-5.2 wt% of potassium was recovered in hydrochar and the highest recovery was observed in HTC-230. Potassium, available as potassium ion (monovalent cation) in liquid, has the highest dissociation constant among the other nutrients. As a result, potassium has the tendency to remain in the liquid even during HTC.
N2 adsorption-desorption isotherms for HTC-180, HTC-230, and HTC-300 are shown in Figure 1. The adsorption isotherms of the hydrochars, regardless of the activation conditions, showed type IV isotherms, indicating presence of large fractions of mesopores [31]. HTC-180 and HTC-230 show more mesopores than HTC-300. The N2 adsorption isotherms in Figure 1 showed that there is a significant downward shift in N2 adsorption for all samples, due to the effect of heat shrinking and oxidation and decomposition of the products. The trend is similar as found in the desorption isotherms. The increase of nutrient concentration in the hydrochar may decrease the gas uptake even at higher pressure. The similar phenomenon was observed and reported in a previous study [29], where zinc chloride and lithium chloride eutectic salts are used during HTC of glucose at 180 °C.
Lower gas uptake might also be caused by the smaller particle size, which might indicate that the average particle size reduces with HTC temperature. Surface area, pore size, and pore volume were calculated in various methods and the results are listed in Table 2. All the hydrochars have very similar multipoint BET surface area in the range of 1.0-3.5 m2 g−1. However, the BJH and DH surface areas are highest for HTC-180 (BJH and DH adsorption surface area 8.0 and 8.1 m2 g−1, respectively) and lowest for HTC-300 (BJH and DH adsorption surface area 3.8 and 3.9 m2 g−1, respectively). This might be the reason for less pronounced gas uptake by HTC-300 than HTC-180.
Morphology | Method | HTC-180 | HTC-240 | H-300 |
Surface area | Multi point BET | 1.85 | 2.05 | 0.98 |
(m2 g-1) | BJH method cumulative adsorption surface area | 7.96 | 8.26 | 3.80 |
BJH method cumulative desorption surface area | 5.47 | 6.44 | 2.96 | |
DH method cumulative adsorption surface area | 8.08 | 8.40 | 3.87 | |
DH method cumulative desorption surface area | 5.65 | 6.68 | 3.07 | |
t-method external surface area | 1.85 | 2.05 | 0.98 | |
t-method micropore surface area | 1.81 | ND | ND | |
NLDFT cumulative surface area | 0.00 | 2.01 | 0.89 | |
Pore volume | BJH method cumulative adsorption pore volume | 1.10×10-2 | 1.11×10-2 | 5.70×10-3 |
(mL g-1) | BJH method cumulative desorption pore volume | 1.11×10-2 | 1.08×10-2 | 5.54×10-3 |
DH method cumulative adsorption pore volume | 1.06×10-2 | 1.08×10-2 | 5.52×10-3 | |
DH method cumulative desorption pore volume | 6.89×10-4 | 1.05×10-2 | 5.38×10-3 | |
t-method micropore volume | ND | ND | ND | |
HK method cumulative pore volume | 6.89×10-4 | 7.65×10-4 | 3.47×10-4 | |
SF method cumulative pore volume | 7.11×10-4 | 7.89×10-4 | 3.60×10-4 | |
NLDFT method cumulative pore volume | 5.45×10-3 | 5.73×10-3 | 1.94×10-3 | |
Pore size | Average pore Radius | 10.71 | 9.23 | 10.05 |
(nm) | BJH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 |
BJH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
DH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 | |
DH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
HK method pore Radius (Mode) | 0.34 | 0.18 | 0.18 | |
SF method pore Radius (Mode) | 0.18 | 0.18 | 0.23 | |
NLDFT pore Radius (Mode) | 20.09 | 20.09 | 0.75 | |
ND: not determined |
The pore volumes are presented in the Table 2 as well. With the lower BJH surface area, it is expected that the pore volume would be lower which is observed in Table 2. In fact, the BJH and DH pore volume of HTC-300 were around 50% lower than those of HTC-180. Average pore size was very similar for all the hydrochars. In fact, pore size of hydrochars determined by BJH, DH, HK, and SF were around 0.18-0.56 nm. HTC-300 was predominantly mesoporous, as can be found from the NLDFT pore size for HTC-180 and HTC-230 were both 20.09 nm.
From the porosity analysis, it can be concluded that all the hydrochars are mesoporous, which is different than pyrolysis biochar [34]. In biochar, physisorption is more dominant, as the large microporous surface area. However, hydrochar seems to have very low BET surface area and it is mostly mesoporous. Therefore, in the case of nutrient sequestration, it should be chemisorbed with the chemically active sites on the hydrochar surface. To find out the chemisorption sites on the hydrochar surface and possible nutrient bonding with hydrochar carbon, FTIR was carried out for all the hydrochars and the results are discussed below.
FTIR spectroscopy has been extensively used in biomass research, as it shows the bond energies of characteristic groups in biomass, and can indicate changes in molecular formation resulting from various treatments [12]. The FTIR spectra of raw glucose, and different hydrochars were performed and the spectra are shown in Figure 2. Raw glucose shows strong bonds at 896, 1036, 1186, 1244, 1367, 1426, 1618, 2845, and 3300 cm−1 in the spectrum. These correspond to carbon hydrogen bond (C-H), alcohol group (C-OH), aryl-alkyl ether (C-O-C), organic acid (-COOH), carbon skeleton (C=C), carbon hydrogen bond (C-H), and hydroxyl bond (-OH), respectively [12,32]. With HTC, oxygen containing bonds in hydrochars are deemed mainly because the volatile oxygen-rich components degrade into liquid and gaseous products [12].
The elemental analysis, presented in Table 1, shows an increase of elemental carbon, also in complement with the FTIR spectra. Two new peaks at 1600 and 1700 cm−1 corresponding to C=C bonds are produced and sharpen for all hydrochar spectra. This might be another confirmation of carbonization during HTC as previously addressed by various researchers[12,35]. Besides C=C bonds, there is also evidence of organometallic bonds e.g. C-N or C-P bonds in 1030 cm−1, C=S bond in 1170 cm−1, -NH2 or -CH2 twist in 1350 cm−1 etc. In addition, two shifted peaks observed at 745 and 830 cm−1 further confirmed the formation of C-O-Fe bonds due to the interaction between iron and the carbon matrix. They were related to Fe-O stretching modes of the tetrahedral and octahedral sites in its spinel structure [36]. Supporting the ICP-OES results in Table 1, the iron and phosphorous bonds were sharpening with HTC temperature. Thus, combining the elemental analysis and FTIR indicates the formation of nutrient containing carbon spheres during HTC of glucose with salt addition. SEM-EDX may reveal the sites and distribution of nutrients on the hydrochar surface.
It is important to know whether the nutrients chemisorb on the hydrochar or they formed inorganic salts during the cooling phase. One can only speculate from porosity analysis that the nutrients are probably chemisorbed on the hydrochar as the physisorbed salts may have been washed away during washing after hydrochar production. However, to prove this as well as to find the effect of various nutrient concentrations in hydrochar formation, SEM-EDX was performed in all three hydrochars. SEM images are reported in Figure 3, whereas elemental scanning for HTC-180, HTC-230, and HTC 300 is presented in Figures 4, 5, and 6, respectively.
From Figure 3, carbon spheres formed during HTC can be observed. The starting feedstock was glucose dissolved in saline solution. The spheres were nucleated and developed during HTC. The similar carbon spheres produced during HTC of saccharides were previously reported [25,26,27,31]. It has also been reported that the size of carbon spheres are around 1 ± 0.8 µm, when glucose was hydrothermally carbonized from 170-240 °C for 1-15h[31]. The carbon spheres produced in this study are clearly larger compared to that previous study. In fact, a higher degree of aromatization was reported earlier in the presence of eutectic salt solution during HTC of polysaccharides [29]. The presence of salts during HTC yields higher energy value, according to another study [37]. In fact, iron in the solution might catalyze the HTC and probably nucleate the formation of carbon spheres. As a result, the carbon spheres are larger and contain nutrients in their structure, which can be observed from the EDX.
Elemental scanning or EDX was performed to find out the nutrients distribution on the hydrochar’s surfaces. Elements of interest for this study are carbon, oxygen, phosphorus, potassium, and iron along with other possible inorganics (due to their presence in the salts) distribution in HTC-180, HTC-230, and HTC-300 are presented in Figures 4, 5, and 6, respectively. As carbon tape was used to attach the hydrochar in the stud during SEM-EDX, carbon can be observed on the other sites besides the samples as well. From Figure 4-6, carbon and oxygen were dominating in the hydrochars, regardless of HTC temperature. Phosphorus and iron were very widely distributed in all of the hydrochars. As discussed earlier (section 3.1), potassium is unlikely to chemisorb in the solid hydrochar during HTC, there is no visible evidence of potassium in the EDX of hydrochars at any HTC temperature. Chromium, the impurity of iron salt was also visible in the EDX in all hydrochars as well. According to Table 1, sulfur concentration increases around three times in HTC-300 than HTC-180. As a result, there is some evidence of sulfur in the EDX of HTC-300 hydrochar.
A number of studies has been reported the HTC reaction pathways to produce hydrochar from model compounds (e.g., sugars, cellulose) to real biomass (e.g., wood) [18,20,31,38]. Regardless of the feedstock type, four major HTC reactions occur in this following sequence: (i) dehydration of sugars into furan compounds, (ii) polymerization or condensation of furan-derivatives, (iii) aromatization of polymerized products, and finally (iv) solidification of macromolecules into macromolecules [33,38]. Based on the FTIR results previously explained in Section 3.3, it is evident that similar reaction mechanism takes place for sugar molecules with the presence of salinity as well. In addition, nutrient sequestration on the surface take place via chemisorption. The complete reaction mechanism is shown in Figure 7. Previously, it was reported that hydrochar spheres consist of a hydrophobic core and hydrophilic surface with various oxygen functional groups [18,31]. Due to their charge density, these functional groups (esters, acids, ketones etc.) are behaving as chemisorption sites (section 3.3) [22,23]. Charged nutrient molecules (NH4+, PO43−) are attracted by the surface sites and adsorb as shown in Figure 7. These functional groups on the surface chemisorb nutrient ions, which are also found from FTIR spectra. The longer reaction temperature results more functional groups on the solid surface. Therefore, more nutrients are sequestered on the hydrochar surface.
In presence of saline solution, hydrochar yield from glucose increases with HTC temperature, although carbon concentration in the hydrochar remains similar. Maximum 71, 31, and 23.1 wt% nitrogen, iron, and phosphorus respectively were recovered in solid hydrochar at highest HTC temperature (300 °C). Potassium was unlikely to be recovered in solid during HTC in the temperature range of 180-300 °C. The hydrochars produced during HTC are spherical and their surface area varies from 1.0-3.5 m2 g−1, which is very low compared to charcoal or pyrolysis biochar. Therefore, physisorption is highly unlikely but the sequestration of nutrients indicates chemisorption. FTIR and SEM-DEX provide supporting information on chemisorption. However, the optimum sequestration and kinetics of sequestration in terms of specific nutrients needs to be further investigated.
This research is supported by funds from the Bioenergy 2021 program delegated from the German Federal Ministry of Research and Education (BMBF) to Project Management Juelich (PtJ). Additionally, the authors acknowledge Western Sun Grant Initiative (grant no. 2010-38502-21839) for the financial support. The authors also acknowledge the analytical support of Zeta Partikel Analytik, Germany, especially Dr. Andreas Hahn for SEM and nitrogen sorption analyses.
All authors declare no conflicts of interest in this paper.
[1] | Conti J, Holtberg P, Diefenderfer J, et al. (2016) International energy outlook 2016 with projections to 2040. USDOE Energy Information Administration (EIA), Washington, DC (United States), 2016. |
[2] |
Sims RE, Rogner HH, Gregory K (2003) Carbon emission and mitigation cost comparisons between fossil fuel, nuclear and renewable energy resources for electricity generation. Energy Policy 31: 1315-1326. doi: 10.1016/S0301-4215(02)00192-1
![]() |
[3] | W. W. E. A. (February 25, 2019) Wind Power Capacity Worldwide Reaches 597 GW, 50,1 GW added in 2018. |
[4] |
Kaviani H, Nejat A (2017) Aeroacoustic and aerodynamic optimization of a MW class HAWT using MOPSO algorithm. Energy 140: 1198-1215. doi: 10.1016/j.energy.2017.08.011
![]() |
[5] |
Ashrafi ZN, Ghaderi M, Sedaghat A (2015) Parametric study on off-design aerodynamic performance of a horizontal axis wind turbine blade and proposed pitch control. Energy Convers Manage 93: 349-356. doi: 10.1016/j.enconman.2015.01.048
![]() |
[6] |
Lydia M, Kumar SS, Selvakumar AI, et al. (2014) A comprehensive review on wind turbine power curve modeling techniques. Renewable Sustainable Energy Rev 30: 452-460. doi: 10.1016/j.rser.2013.10.030
![]() |
[7] |
Tchakoua P, Wamkeue R, Ouhrouche M, et al. (2014) Wind turbine condition monitoring: State-of-the-art review, new trends, and future challenges. Energies 7: 2595-2630. doi: 10.3390/en7042595
![]() |
[8] |
Sanderse B, van der Pijl SP, Koren B (2011). Review of computational fluid dynamics for wind turbine wake aerodynamics. Wind Energy 14: 799-819. doi: 10.1002/we.458
![]() |
[9] | Infield D, Freris L (2020). Renewable energy in power systems, John Wiley & Sons. |
[10] |
Bai CJ, Chen PW, Wang WC (2016) Aerodynamic design and analysis of a 10 kW horizontal-axis wind turbine for Tainan, Taiwan. Clean Technol Environ Policy 18: 1151-1166. doi: 10.1007/s10098-016-1109-z
![]() |
[11] | Manwell JF, McGowan JG, Rogers AL (2010) Wind energy explained: theory, design and application, John Wiley & Sons. |
[12] |
Katsigiannis YA, Stavrakakis GS (2014). Estimation of wind energy production in various sites in Australia for different wind turbine classes: A comparative technical and economic assessment. Renewable Energy 67: 230-236. doi: 10.1016/j.renene.2013.11.051
![]() |
[13] | Masters GM (2013) Renewable and efficient electric power systems, John Wiley & Sons. |
[14] |
Rocha PAC, de Sousa RC, de Andrade CF, et al. (2012) Comparison of seven numerical methods for determining Weibull parameters for wind energy generation in the northeast region of Brazil. Appl Energy 89: 395-400. doi: 10.1016/j.apenergy.2011.08.003
![]() |
[15] |
Seguro J, Lambert T (2000) Modern estimation of the parameters of the Weibull wind speed distribution for wind energy analysis. J Wind Eng Ind Aerodyn 85: 75-84. doi: 10.1016/S0167-6105(99)00122-1
![]() |
[16] | Laban ON, Maghanga CM, Joash K (2019) Determination of the surface roughness parameter and wind shear exponent of Kisii Region from the On-Site measurement of wind profiles. J Energy 2019. |
[17] |
Wang L, Cholette ME, Tan AC, et al. (2017) A computationally-efficient layout optimization method for real wind farms considering altitude variations. Energy 132: 147-159. doi: 10.1016/j.energy.2017.05.076
![]() |
[18] |
Islam M, Saidur R, Rahim N (2011) Assessment of wind energy potentiality at Kudat and Labuan, Malaysia using Weibull distribution function. Energy 36: 985-992. doi: 10.1016/j.energy.2010.12.011
![]() |
[19] | Krenn A, Winkelmeier H, Cattin R, et al. (2010) Austrian wind atlas and wind potential analysis. DEWEK, Available from: www. windatlas. at/downloads/20101117_Paper_Dewek. pdf. |
[20] |
Celik AN (2004) A statistical analysis of wind power density based on the Weibull and Rayleigh models at the southern region of Turkey. Renewable Energy 29: 593-604. doi: 10.1016/j.renene.2003.07.002
![]() |
[21] |
Mentis D, Hermann S, Howells M, et al. (2015) Assessing the technical wind energy potential in Africa a GIS-based approach. Renewable Energy 83: 110-125. doi: 10.1016/j.renene.2015.03.072
![]() |
[22] |
Ma PC, Zhang Y (2014) Perspectives of carbon nanotubes/polymer nanocomposites for wind blade materials. Renewable Sustainable Energy Rev 30: 651-660. doi: 10.1016/j.rser.2013.11.008
![]() |
[23] | Perez-Blanco H (2011) Optimization of wind energy capture using BET. Proc., ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition, American Society of Mechanical Engineers Digital Collection, 879-887. |
[24] |
Duquette MM, Visser KD (2003) Numerical implications of solidity and blade number on rotor performance of horizontal-axis wind turbines. J Sol Energy Eng 125: 425-432. doi: 10.1115/1.1629751
![]() |
[25] |
McKenna R, vd Leye PO, Fichtner W (2016) Key challenges and prospects for large wind turbines. Renewable Sustainable Energy Rev 53: 1212-1221. doi: 10.1016/j.rser.2015.09.080
![]() |
[26] | Letcher TM (2017) Wind energy engineering: A handbook for onshore and offshore wind turbines, Academic Press. |
[27] | Heier S (2014) Grid integration of wind energy: Onshore and offshore conversion systems, John Wiley & Sons. |
[28] | Boldea I, Tutelea LN (2009) Electric machines: steady state, transients, and design with MATLAB, CRC press. |
[29] | Orabi M, El-Sousy F, Godah H, et al. (2004) High-performance induction generator-wind turbine connected to utility grid. Proc., INTELEC 2004. 26th Annual International Telecommunications Energy Conference, IEEE, 697-704. |
[30] | Mishnaevsky Jr L, Favorsky O (2011) Composite materials in wind energy technology. Thermal to Mechanical Energy Conversion: Engines and Requirements, EOLSS Publishers: Oxford, UK. |
[31] | Ashwill T Materials and innovations for large blade structures: research opportunities in wind energy technology. Proc., 50th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 17th AIAA/ASME/AHS Adaptive Structures Conference 11th AIAA No, 2407. |
[32] | Thresher RW, Dodge DM (1998) Trends in the evolution of wind turbine generator configurations and systems. Wind Energy: Int J Prog Appl Wind Power Convers Technol 1: 70-86. |
[33] |
Sieros G, Chaviaropoulos P, Sørensen JD, et al. (2012) Upscaling wind turbines: theoretical and practical aspects and their impact on the cost of energy. Wind Energy 15: 3-17. doi: 10.1002/we.527
![]() |
[34] | Al-Abadi A, Ertunç Ö, Beyer F, et al. Torque-matched aerodynamic shape optimization of HAWT rotor. Proc, J Physics: Conference Series, IOP Publishing, 012003. |
[35] |
Al-Abadi A, Ertunç Ö, Weber H, et al. (2015) A design and optimization method for matching the torque of the wind turbines. J Renewable Sustainable Energy 7: 023129. doi: 10.1063/1.4918748
![]() |
[36] |
Chattot JJ (2003) Optimization of wind turbines using helicoidal vortex model. J Sol Energy Eng 125: 418-424. doi: 10.1115/1.1621675
![]() |
[37] |
Jureczko M, Pawlak M, Mężyk A (2005) Optimisation of wind turbine blades. J Mater Process Technol 167: 463-471. doi: 10.1016/j.jmatprotec.2005.06.055
![]() |
[38] |
Henriques J, Da Silva FM, Estanqueiro, et al. (2009) Design of a new urban wind turbine airfoil using a pressure-load inverse method. Renewable Energy 34: 2728-2734. doi: 10.1016/j.renene.2009.05.011
![]() |
[39] |
Vardar A, Alibas I (2008) Research on wind turbine rotor models using NACA profiles. Renewable Energy 33: 1721-1732. doi: 10.1016/j.renene.2007.07.009
![]() |
[40] |
Ameku K, Nagai BM, Roy JN (2008). Design of a 3 kW wind turbine generator with thin airfoil blades. Exp Therm Fluid Sci 32: 1723-1730. doi: 10.1016/j.expthermflusci.2008.06.008
![]() |
[41] | Leung D, Deng Y, Leung M (2010) Design optimization of a cost-effective micro wind turbine. Proc, WCE 2010-World Congress on Engineering 2010, International Association of Engineers. |
[42] |
Hirahara H, Hossain MZ, Kawahashi M, et al. (2005) Testing basic performance of a very small wind turbine designed for multi-purposes. Renewable Energy 30: 1279-1297. doi: 10.1016/j.renene.2004.10.009
![]() |
[43] |
Liu X, Wang L, Tang X (2013) Optimized linearization of chord and twist angle profiles for fixed-pitch fixed-speed wind turbine blades. Renewable Energy 57: 111-119. doi: 10.1016/j.renene.2013.01.036
![]() |
[44] |
Sedaghat A, Mirhosseini M (2012) Aerodynamic design of a 300 kW horizontal axis wind turbine for province of Semnan. Energy Convers Manage 63: 87-94. doi: 10.1016/j.enconman.2012.01.033
![]() |
[45] |
Darwish AS, Shaaban S, Marsillac E, et al. (2019) A methodology for improving wind energy production in low wind speed regions, with a case study application in Iraq. Comput Ind Eng 127: 89-102. doi: 10.1016/j.cie.2018.11.049
![]() |
[46] |
Derakhshan S, Tavaziani A, Kasaeian N (2015) Numerical shape optimization of a wind turbine blades using artificial bee colony algorithm. J Energy Resour Technol 137: 051210. doi: 10.1115/1.4031043
![]() |
[47] | Alpman E (2014) Effect of selection of design parameters on the optimization of a horizontal axis wind turbine via genetic algorithm. Institute of Physics Publishing. |
[48] |
Kong C, Bang J, Sugiyama Y (2005) Structural investigation of composite wind turbine blade considering various load cases and fatigue life. Energy 30: 2101-2114. doi: 10.1016/j.energy.2004.08.016
![]() |
[49] |
Seo S, Oh SD, Kwak HY (2019) Wind turbine power curve modeling using maximum likelihood estimation method. Renewable Energy 136: 1164-1169. doi: 10.1016/j.renene.2018.09.087
![]() |
[50] |
Kishore RA, Priya S (2013) Design and experimental verification of a high efficiency small wind energy portable turbine (SWEPT). J Wind Eng Ind Aerodyn 118: 12-19. doi: 10.1016/j.jweia.2013.04.009
![]() |
[51] |
Singh RK, Ahmed MR, Zullah MA, et al. (2012) Design of a low Reynolds number airfoil for small horizontal axis wind turbines. Renewable Energy 42: 66-76. doi: 10.1016/j.renene.2011.09.014
![]() |
[52] | Gasch R, Twele J (2011) Wind power plants: fundamentals, design, construction and operation, Springer Science & Business Media. |
[53] | Scholbrock A, Fleming P, Fingersh L, et al., Field testing LIDAR-based feed-forward controls on the NREL controls advanced research turbine. Proc., 51st AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, 818. |
[54] |
Oerlemans S, Sijtsma P, Méndez López B (2007) Location and quantification of noise sources on a wind turbine. J Sound Vib 299: 869-883. doi: 10.1016/j.jsv.2006.07.032
![]() |
[55] |
Deshmukh AP, Allison JT (2016) Multidisciplinary dynamic optimization of horizontal axis wind turbine design. Struct Multidiscip Optim 53: 15-27. doi: 10.1007/s00158-015-1308-y
![]() |
[56] | Tong W (2010) Wind power generation and wind turbine design, WIT press. |
[57] |
Sedaghat A, Mirhosseini M (2012) Aerodynamic design of a 300 kW horizontal axis wind turbine for province of Semnan. Energy Convers Manage 63: 87-94. doi: 10.1016/j.enconman.2012.01.033
![]() |
[58] | Hau E, von Renouard H (2003) Wind turbines: fundamentals, technologies, application, economics, Springer. |
[59] |
Devinant P, Laverne T, Hureau J (2002) Experimental study of wind-turbine airfoil aerodynamics in high turbulence. J Wind Eng Ind Aerodyn 90: 689-707. doi: 10.1016/S0167-6105(02)00162-9
![]() |
[60] | Bertagnolio F, Sørensen NN, Johansen J, et al. (2001) Wind turbine airfoil catalogue. |
[61] |
Fupeng H, Yuhong L, Zuoyi C (2001) Suggestions for improving wind turbines blade characteristics. Wind Eng 25: 105-113. doi: 10.1260/0309524011495908
![]() |
[62] | Lawson MJ, Li Y, Sale DC (2011) Development and verification of a computational fluid dynamics model of a horizontal-axis tidal current turbine. Proc, ASME 2011 30th international conference on ocean, offshore and arctic engineering, American Society of Mechanical Engineers Digital Collection, 711-720. |
[63] | Yılmaz M, Köten H, Çetinkaya E, et al. (2018) A comparative CFD analysis of NACA0012 and NACA4412 airfoils. J Energy Syst 2: 145-159. |
[64] |
Fuglsang P, Madsen HA (1999) Optimization method for wind turbine rotors. J Wind Eng Ind Aerodyn 80: 191-206. doi: 10.1016/S0167-6105(98)00191-3
![]() |
[65] | Fuglsang P, Bak C (2004) Development of the Risø wind turbine airfoils. Wind Energy: Int J Prog Appl Wind Power Convers Technol 7: 145-162. |
[66] | Dahl KS, Fuglsang P (1998) Design of the wind turbine airfoil family RISØ-A-XX. |
[67] | Timmer W, Rooij R, Summary of the Delft University wind turbine dedicated airfoils. Proc, 41st aerospace sciences meeting and exhibit, 352. |
[68] | Tangler JL, Somers DM (1995) NREL airfoil families for HAWTs. National Renewable Energy Lab., Golden, CO (United States). |
[69] | Somers DM (2005) S833, S834, and S835 Airfoils: November 2001--November 2002. National Renewable Energy Lab.(NREL), Golden, CO (United States). |
[70] |
Sagol E, Reggio M, Ilinca A (2013) Issues concerning roughness on wind turbine blades. Renewable Sustainable Energy Rev 23: 514-525. doi: 10.1016/j.rser.2013.02.034
![]() |
[71] | Somers DM, Tangler JL (1995) Wind-tunnel test of the S814 thick root airfoil. National Renewable Energy Lab., Golden, CO (United States). |
[72] |
Van Rooij R, Timmer W (2003) Roughness sensitivity considerations for thick rotor blade airfoils. J Sol Energy Eng 125: 468-478. doi: 10.1115/1.1624614
![]() |
[73] |
Lanzafame R, Messina M (2009) Design and performance of a double-pitch wind turbine with non-twisted blades. Renewable Energy 34: 1413-1420. doi: 10.1016/j.renene.2008.09.004
![]() |
[74] | Laursen J, Enevoldsen P, Hjort S, 3D CFD rotor computations of a multi-megawatt HAWT rotor. Proc., Proceedings of the European Wind Energy Conference, Milan, Italy. |
[75] |
Jeong JH, Kim SH (2018) CFD investigation on the flatback airfoil effect of 10 MW wind turbine blade. J Mech Sci Technol 32: 2089-2097. doi: 10.1007/s12206-018-0418-z
![]() |
[76] |
Ahmed MR, Narayan S, Zullah MA, et al. (2011) Experimental and numerical studies on a low Reynolds number airfoil for wind turbine blades. J Fluid Sci Technol 6: 357-371. doi: 10.1299/jfst.6.357
![]() |
[77] |
Sayed MA, Kandil HA, Shaltot A (2012) Aerodynamic analysis of different wind-turbine-blade profiles using finite-volume method. Energy Convers Manage 64: 541-550. doi: 10.1016/j.enconman.2012.05.030
![]() |
[78] |
Sicot C, Devinant P, Laverne T, et al. (2006) Experimental study of the effect of turbulence on horizontal axis wind turbine aerodynamics. Wind Energy 9: 361-370. doi: 10.1002/we.184
![]() |
[79] | Delnero J, Marañon di Leo J, Bacchi F, et al. (2005) Experimental determination of the influence of turbulent scale on the lift and drag coefficients of low Reynolds number airfoils. Lat Am Appl Res 35: 183-188. |
[80] | Swalwell KE, Sheridan J, Melbourne W, The effect of turbulence intensity on stall of the NACA 0021 aerofoil. Proc, 14th Australasian Fluid Mechanics Conference, 10-14. |
[81] |
Larsen JW, Nielsen SR, Krenk S (2007) Dynamic stall model for wind turbine airfoils. J Fluids Struct 23: 959-982. doi: 10.1016/j.jfluidstructs.2007.02.005
![]() |
[82] |
Hoffmann JA (1991) Effects of freestream turbulence on the performance characteristics of an airfoil. AIAA J 29: 1353-1354. doi: 10.2514/3.10745
![]() |
[83] |
Kamada Y, Maeda T, Murata J, et al. (2011) Effects of turbulence intensity on dynamic characteristics of wind turbine airfoil. J Fluid Sci Technol 6: 333-341. doi: 10.1299/jfst.6.333
![]() |
[84] | Eke G, Onyewudiala J (2010) Optimization of wind turbine blades using genetic algorithm. Global J Res Eng 10. |
[85] | Berg D, Zayas J, Aerodynamic and aeroacoustic properties of flatback airfoils. Proc, 46th AIAA Aerospace Sciences Meeting and Exhibit, 1455. |
[86] | McWilliam M, Crawford C, Manufacturing defect effects on Bend-Twist coupled wind turbine blades. Proc, 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, 1062. |
[87] |
Lee SG, Park SJ, Lee KS, et al. (2012) Performance prediction of NREL (National Renewable Energy Laboratory) Phase VI blade adopting blunt trailing edge airfoil. Energy 47: 47-61. doi: 10.1016/j.energy.2012.08.007
![]() |
[88] |
Kim SH, Bang HJ, Shin HK, et al. (2014) Composite structural analysis of flat-back shaped blade for multi-MW class wind turbine. Appl Compos Mater 21: 525-539. doi: 10.1007/s10443-013-9362-3
![]() |
[89] |
Fischer GR, Kipouros T, Savill AM (2014) Multi-objective optimisation of horizontal axis wind turbine structure and energy production using aerofoil and blade properties as design variables. Renewable Energy 62: 506-515. doi: 10.1016/j.renene.2013.08.009
![]() |
[90] |
Lago LI, Ponta FL, Otero AD (2013) Analysis of alternative adaptive geometrical configurations for the NREL-5 MW wind turbine blade. Renewable Energy 59: 13-22. doi: 10.1016/j.renene.2013.03.007
![]() |
[91] | Kahn DL, van Dam C, Berg DE (2008) Trailing edge modifications for flatback airfoils. Sandia National Laboratories. |
[92] | Baker JP, Van Dam C (2009) Drag reduction of a blunt trailing-edge airfoil, University of California, Davis. |
[93] | Cooperman A, McLennan A, Baker J, et al., Aerodynamic performance of thick blunt trailing edge airfoils. Proc, 28th AIAA Applied Aerodynamics Conference, 4228. |
[94] |
Chow R, Van Dam C (2013) Computational investigations of blunt trailing‐edge and twist modifications to the inboard region of the NREL 5 MW rotor. Wind Energy 16: 445-458. doi: 10.1002/we.1505
![]() |
[95] | Schreck S, Fingersh L, Siegel K, et al., Rotational Augmentation on a 2.3 MW Rotor Blade with Thick Flatback Airfoil Cross Sections. Proc, 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, 915. |
[96] | Murcia JP, Pinilla Á (2011) CFD analysis of blunt trailing edge airfoils obtained with several modification methods. Revista de Ingeniería, 14-24. |
[97] |
Standish K, Van Dam C (2003) Aerodynamic analysis of blunt trailing edge airfoils. J Sol Energy Eng 125: 479-487. doi: 10.1115/1.1629103
![]() |
[98] | Law S, Gregorek G (1987) Wind tunnel evaluation of a truncated NACA 64-621 airfoil for wind turbine applications. |
[99] | Homsrivaranon K (2016) Investigation of Active Flow Control on an Extremely Thick Wind Turbine Airfoil. University of Kansas. |
[100] |
Zhang Tt, Huang W, Wang ZG, et al. (2016) A study of airfoil parameterization, modeling, and optimization based on the computational fluid dynamics method. J Zhejiang University-SCIENCE A 17: 632-645. doi: 10.1631/jzus.A1500308
![]() |
[101] |
Baker J, Mayda E, Van Dam C (2006) Experimental analysis of thick blunt trailing-edge wind turbine airfoils. J sol Energy Eng 128: 422-431. doi: 10.1115/1.2346701
![]() |
[102] |
Göçmen T, Özerdem B (2012) Airfoil optimization for noise emission problem and aerodynamic performance criterion on small scale wind turbines. Energy 46: 62-71. doi: 10.1016/j.energy.2012.05.036
![]() |
[103] |
Kim T, Jeon M, Lee S, et al. (2014) Numerical simulation of flatback airfoil aerodynamic noise. Renewable Energy 65: 192-201. doi: 10.1016/j.renene.2013.08.036
![]() |
[104] | Velte CM, Hansen MOL, Meyer KE, at al., Evaluation of the performance of vortex generators on the DU 91-W2-250 profile using stereoscopic PIV. Proc, International Symposium on Energy, Informatics and Cybernetics: Focus Symposium in the 12th World Multiconference on Systemics, Cybernetics and Informatics (WMSCI 2008). |
[105] |
Chamorro LP, Arndt R, Sotiropoulos F (2013) Drag reduction of large wind turbine blades through riblets: Evaluation of riblet geometry and application strategies. Renewable Energy 50: 1095-1105. doi: 10.1016/j.renene.2012.09.001
![]() |
[106] | Ceyhan O, Timmer W, Experimental evaluation of a non-conventional flat back thick airfoil concept for large offshore wind turbines. Proc, 2018 Applied Aerodynamics Conference, 3827. |
[107] |
Kim T, Lee S (2014) Aeroacoustic simulations of a blunt trailing-edge wind turbine airfoil. J Mech Sci Technol 28: 1241-1249. doi: 10.1007/s12206-014-0114-6
![]() |
[108] | Fuglsang P, Sangill O, Hansen P (2002) Design of a 21 m blade with risø-a1 airfoils for active stall controlled wind turbines. |
[109] | Johansen J, Madsen HA, Gaunaa M, at al. (2009) Design of a wind turbine rotor for maximum aerodynamic efficiency. Wind Energy: Int J Prog Appl Wind Power Convers Technol 12: 261-273. |
[110] |
Pinto RN, Afzal A, D'Souza LV, et al. (2017) Computational fluid dynamics in turbomachinery: A review of state of the art. Arch Comput Methods Eng 24: 467-479. doi: 10.1007/s11831-016-9175-2
![]() |
[111] | Mikkelsen R, Sørensen JN, Øye S, et al., Analysis of power enhancement for a row of wind turbines using the actuator line technique. Proc, J Physics: Conference Series, IOP Publishing, 012044. |
[112] | Snel H (2003) Review of aerodynamics for wind turbines. Wind Energy: Int J Prog Appl Wind Power Convers Technol 6: 203-211. |
[113] | Zhou P (2017) CFD Simulation of the Wind Turbine Wake Under Different Atmospheric Boundary Conditions. Purdue University. |
[114] | Schmidt S, McIver D, Blackburn HM, et al., Spectral element based simulation of turbulent pipe flow. Proc, 14th A/Asian Fluid Mech. Conf. |
[115] | Sargsyan A (2010) Simulation and modeling of flow field around a horizontal axis wind turbine (HAWT) using RANS method, Florida Atlantic University. |
[116] | Launder BE, Sharma B (1974) Application of the energy-dissipation model of turbulence to the calculation of flow near a spinning disc. Lett Heat Mass Transfer 1: 131-137. |
[117] | Shih TH, Liou WW, Shabbir A, et al. (1994) A new k-epsilon eddy viscosity model for high Reynolds number turbulent flows: Model development and validation. |
[118] |
Yakhot V, Orszag S, Thangam S, et al. (1992) Development of turbulence models for shear flows by a double expansion technique. Phys Fluids A: Fluid Dyn 4: 1510-1520. doi: 10.1063/1.858424
![]() |
[119] | Yakhot V, Orszag SA (1986) Renormalization group analysis of turbulence. I. Basic theory. J Sci Comput 1: 3-51. |
[120] |
Mohamed M (2012) Performance investigation of H-rotor Darrieus turbine with new airfoil shapes. Energy 47: 522-530. doi: 10.1016/j.energy.2012.08.044
![]() |
[121] | Mielke A, Naumann J (2018) On the existence of global-in-time weak solutions and scaling laws for Kolmogorov's two-equation model of turbulence. arXiv preprint arXiv:1801.02039. |
[122] |
Wilcox DC (2008) Formulation of the kw turbulence model revisited. AIAA J 46: 2823-2838. doi: 10.2514/1.36541
![]() |
[123] |
Menter FR (2009) Review of the shear-stress transport turbulence model experience from an industrial perspective. Int J Comput Fluid Dyn 23: 305-316. doi: 10.1080/10618560902773387
![]() |
[124] | Smirnov PE, Menter FR (2009) Sensitization of the SST turbulence model to rotation and curvature by applying the Spalart-Shur correction term. J Turbomachinery 131. |
[125] |
Langtry RB, Menter FR (2009) Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes. AIAA J 47: 2894-2906. doi: 10.2514/1.42362
![]() |
[126] | Spalart P, Allmaras SA, One-equation turbulence model for aerodynamic flows. Proc, 30th Aerospace Sciences Meeting And Exhibit, 439. |
[127] | Bouhelal A, Smaïli A, Masson C, et al., Effects of surface roughness on aerodynamic performance of horizontal axis wind turbines. Proc, The 25th Annual Conference of the Computational Fluid Dynamics Society of Canada, CFD2017-337, University of Windsor, 18-21. |
[128] |
Gatski TB, Rumsey CL, Manceau R (2007) Current trends in modelling research for turbulent aerodynamic flows. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 365: 2389-2418. doi: 10.1098/rsta.2007.2015
![]() |
[129] |
Spalart P, Shur M (1997) On the sensitization of turbulence models to rotation and curvature. Aerosp Sci Technol 1: 297-302. doi: 10.1016/S1270-9638(97)90051-1
![]() |
[130] |
Rahman M, Siikonen T, Agarwal R (2011) Improved low-Reynolds-number one-equation turbulence model. AIAA J 49: 735-747. doi: 10.2514/1.J050651
![]() |
[131] |
Deardorff JW (1970) Numerical study of three-dimensional turbulent channel flow at large Reynolds numbers. J Fluid Mech 41: 453-480. doi: 10.1017/S0022112070000691
![]() |
[132] |
Touil H, Hussaini M, Gotoh T, et al. (2007) Development of stochastic models for turbulence. New J Phys 9: 215. doi: 10.1088/1367-2630/9/7/215
![]() |
[133] |
Kraichnan RH (1976) Eddy viscosity in two and three dimensions. J Atmos Sci 33: 1521-1536. doi: 10.1175/1520-0469(1976)033<1521:EVITAT>2.0.CO;2
![]() |
[134] |
Chasnov JR (1991) Simulation of the Kolmogorov inertial subrange using an improved subgrid model. Phys Fluids A: Fluid Dyn 3: 188-200. doi: 10.1063/1.857878
![]() |
[135] |
Piomelli U (1993) High Reynolds number calculations using the dynamic subgrid‐scale stress model. Phys Fluids A: Fluid Dyn 5: 1484-1490. doi: 10.1063/1.858586
![]() |
[136] |
Fröhlich J, Mellen CP, Rodi W, et al. (2005) Highly resolved large-eddy simulation of separated flow in a channel with streamwise periodic constrictions. J Fluid Mech 526: 19-66. doi: 10.1017/S0022112004002812
![]() |
[137] | Spalart PR, Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. Proc, Proceedings of first AFOSR international conference on DNS/LES, Greyden Press. |
[138] |
Spalart PR (2009) Detached-eddy simulation. Annu Rev Fluid Mech 41: 181-202. doi: 10.1146/annurev.fluid.010908.165130
![]() |
[139] | Verhoeven O (2011) Trailing Edge Noise Simulations: Using IDDES in OpenFOAM. |
[140] |
Travin A, Shur M, Strelets M, et al. (2000) Detached-eddy simulations past a circular cylinder. Flow, Turbul Combust 63: 293-313. doi: 10.1023/A:1009901401183
![]() |
[141] |
Sørensen JN (2011) Aerodynamic aspects of wind energy conversion. Annu Rev Fluid Mech 43: 427-448. doi: 10.1146/annurev-fluid-122109-160801
![]() |
[142] |
Li Y, Paik KJ, Xing T, et al. (2012) Dynamic overset CFD simulations of wind turbine aerodynamics. Renewable Energy 37: 285-298. doi: 10.1016/j.renene.2011.06.029
![]() |
[143] |
Lanzafame R, Mauro S, Messina M (2013) Wind turbine CFD modeling using a correlation-based transitional model. Renewable Energy 52: 31-39. doi: 10.1016/j.renene.2012.10.007
![]() |
[144] | Potsdam M, Mavriplis D, Unstructured mesh CFD aerodynamic analysis of the NREL Phase VI rotor. Proc, 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, 1221. |
[145] | Rajvanshi D, Baig R, Pandya R, et al., Wind turbine blade aerodynamics and performance analysis using numerical simulations. Proc, 11th Asian International Conference on Fluid Machinery. |
[146] | Moshfeghi M, Song YJ, Xie YH (2012) Effects of near-wall grid spacing on SST-K-ω model using NREL Phase VI horizontal axis wind turbine. J Wind Eng Ind Aerodyn 107: 94-105. |
[147] |
El Kasmi A, Masson C (2008) An extended k-ε model for turbulent flow through horizontal-axis wind turbines. J Wind Eng Ind Aerodyn 96: 103-122. doi: 10.1016/j.jweia.2007.03.007
![]() |
[148] |
AbdelSalam AM, Ramalingam V (2014) Wake prediction of horizontal-axis wind turbine using full-rotor modeling. J Wind Eng Ind Aerodyn 124: 7-19. doi: 10.1016/j.jweia.2013.11.005
![]() |
[149] | Rütten M, Penneçot J, Wagner C (2009) Unsteady Numerical Simulation of the Turbulent Flow around a Wind Turbine. Progress in Turbulence III, Springer, 103-106. |
[150] |
Abdulqadir SA, Iacovides H, Nasser A (2017) The physical modelling and aerodynamics of turbulent flows around horizontal axis wind turbines. Energy 119: 767-799. doi: 10.1016/j.energy.2016.11.060
![]() |
[151] |
You JY, Yu DO, Kwon OJ (2013) Effect of turbulence models on predicting HAWT rotor blade performances. J Mech Sci Technol 27: 3703-3711. doi: 10.1007/s12206-013-0916-y
![]() |
[152] |
Bouhelal A, Smaili A, Guerri O, et al. (2018) Numerical investigation of turbulent flow around a recent horizontal axis wind Turbine using low and high Reynolds models. J Appl Fluid Mech 11: 151-164. doi: 10.29252/jafm.11.01.28074
![]() |
[153] | Sørensen NN, Michelsen J, Schreck S (2002) Navier-Stokes predictions of the NREL phase VI rotor in the NASA Ames 80 ft×120 ft wind tunnel. Wind Energy: Int J Prog Appl Wind Power Convers Technol 5: 151-169. |
[154] | Johansen J, Sørensen NN, Michelsen J, et al. (2002) Detached‐eddy simulation of flow around the NREL Phase VI blade. Wind Energy: Int J Prog Appl Wind Power Convers Technol 5: 185-197. |
[155] |
Duque EP, Burklund MD, Johnson W (2003) Navier-Stokes and comprehensive analysis performance predictions of the NREL phase VI experiment. J Sol Energy Eng 125: 457-467. doi: 10.1115/1.1624088
![]() |
[156] | Johansen J, Sørensen NN (2004) Aerofoil characteristics from 3D CFD rotor computations. Wind Energy: Int J Prog Appl Wind Power Convers Technol 7: 283-294. |
[157] | Mandas N, Cambuli F, Carcangiu CE (2006) Numerical prediction of horizontal axis wind turbine flow. University of Caglairi, EWEC. |
[158] | Sezer-Uzol N, Long L, 3-D time-accurate CFD simulations of wind turbine rotor flow fields. Proc, 44th AIAA Aerospace Sciences Meeting and Exhibit, 394. |
[159] |
Hu D, Hua O, Du Z (2006) A study on stall-delay for horizontal axis wind turbine. Renewable Energy 31: 821-836. doi: 10.1016/j.renene.2005.05.002
![]() |
[160] | Simms D, Robinson M, Hand M, et al. (1995) A comparison of baseline aerodynamic performance of optimally-twisted versus non-twisted HAWT blades. National Renewable Energy Lab., Golden, CO (United States). |
[161] | Wußow S, Sitzki L, Hahm T, 3D-simulation of the turbulent wake behind a wind turbine. Proc, J Physics: Conference Series, IOP Publishing, 012033. |
[162] |
Thumthae C, Chitsomboon T (2009) Optimal angle of attack for untwisted blade wind turbine. Renewable Energy 34: 1279-1284. doi: 10.1016/j.renene.2008.09.017
![]() |
[163] | Fletcher TM, Brown R, Kim DH, et al. (2009) Predicting wind turbine blade loads using vorticity transport and RANS methodologies. Proc, European Wind Energy Conference and Exhibition, EWEC 2009. |
[164] | Sørensen NN (2009) CFD modelling of laminar‐turbulent transition for airfoils and rotors using the γ− model. Wind Energy: Int J Prog Appl Wind Power Convers Technol 12: 715-733. |
[165] | Gomez-Iradi S, Steijl R, Barakos G (2009) Development and validation of a CFD technique for the aerodynamic analysis of HAWT. J Sol Energy Eng, 131. |
[166] |
Tachos N, Filios A, Margaris D (2010) A comparative numerical study of four turbulence models for the prediction of horizontal axis wind turbine flow. Proc Inst Mech Eng Part C: J Mech Eng Sci 224: 1973-1979. doi: 10.1243/09544062JMES1901
![]() |
[167] |
Fu P, Farzaneh M (2010) A CFD approach for modeling the rime-ice accretion process on a horizontal-axis wind turbine. J Wind Eng Ind Aerodyn 98: 181-188. doi: 10.1016/j.jweia.2009.10.014
![]() |
[168] |
Bechmann A, Sørensen NN, Zahle F (2011) CFD simulations of the MEXICO rotor. Wind Energy 14: 677-689. doi: 10.1002/we.450
![]() |
[169] |
Elfarra MA, Sezer‐Uzol N, Akmandor IS (2014) NREL VI rotor blade: numerical investigation and winglet design and optimization using CFD. Wind Energy 17: 605-626. doi: 10.1002/we.1593
![]() |
[170] |
Abdelsalam AM, Boopathi K, Gomathinayagam S, et al. (2014) Experimental and numerical studies on the wake behavior of a horizontal axis wind turbine. J Wind Eng Ind Aerodyn 128: 54-65. doi: 10.1016/j.jweia.2014.03.002
![]() |
[171] |
Song Y, Perot JB (2015) Cfd simulation of the nrel phase vi rotor. Wind Eng 39: 299-309. doi: 10.1260/0309-524X.39.3.299
![]() |
[172] |
Derakhshan S, Tavaziani A (2015) Study of wind turbine aerodynamic performance using numerical methods. J Clean Energy Technol 3: 83-90. doi: 10.7763/JOCET.2015.V3.174
![]() |
[173] | Sørensen NN, Zahle F, Boorsma K, et al. CFD computations of the second round of MEXICO rotor measurements. Proc, J Physics: Conference Series, IOP Publishing, 022054. |
[174] |
Wang L, Quant R, Kolios A (2016) Fluid structure interaction modelling of horizontal-axis wind turbine blades based on CFD and FEA. J Wind Eng Ind Aerodyn 158: 11-25. doi: 10.1016/j.jweia.2016.09.006
![]() |
[175] |
Menegozzo L, Dal Monte A, Benini E, et al. (2018) Small wind turbines: A numerical study for aerodynamic performance assessment under gust conditions. Renewable Energy 121: 123-132. doi: 10.1016/j.renene.2017.12.086
![]() |
[176] |
Bolinger M, Wiser R (2009) Wind power price trends in the United States: struggling to remain competitive in the face of strong growth. Energy Policy 37: 1061-1071. doi: 10.1016/j.enpol.2008.10.053
![]() |
[177] | Giebel G, Brownsword R, Kariniotakis G, et al. (2011) The state-of-the-art in short-term prediction of wind power: A literature overview. |
[178] |
Albadi M, El-Saadany E (2010) Overview of wind power intermittency impacts on power systems. Electr Power Syst Res 80: 627-632. doi: 10.1016/j.epsr.2009.10.035
![]() |
[179] |
Cavallo A (2007) Controllable and affordable utility-scale electricity from intermittent wind resources and compressed air energy storage (CAES). Energy 32: 120-127. doi: 10.1016/j.energy.2006.03.018
![]() |
[180] |
Sovacool BK (2009) The intermittency of wind, solar, and renewable electricity generators: Technical barrier or rhetorical excuse? Utilities Policy 17: 288-296. doi: 10.1016/j.jup.2008.07.001
![]() |
[181] |
Spinato F, Tavner PJ, Van Bussel GJ, et al. (2009) Reliability of wind turbine subassemblies. IET Renewable Power Gener 3: 387-401. doi: 10.1049/iet-rpg.2008.0060
![]() |
[182] | Walford CA (2006) Wind turbine reliability: understanding and minimizing wind turbine operation and maintenance costs. Sandia National Laboratories. |
[183] |
Faulstich S, Hahn B, Tavner PJ (2011) Wind turbine downtime and its importance for offshore deployment. Wind Energy 14: 327-337. doi: 10.1002/we.421
![]() |
[184] | Boyle G (2009) Renewable electricity and the grid: the challenge of variability, Routledge. |
[185] | Piwko R, MillerN, Sanchez-Gasca J, et al. (2006) Integrating large wind farms into weak power grids with long transmission lines. Proc, 2006 CES/IEEE 5th International Power Electronics and Motion Control Conference, IEEE, 1-7. |
[186] |
Perveen R, Kishor N, Mohanty SR (2014) Off-shore wind farm development: Present status and challenges. Renewable Sustainable Energy Rev 29: 780-792. doi: 10.1016/j.rser.2013.08.108
![]() |
[187] |
Barlas TK, van Kuik GA (2010) Review of state of the art in smart rotor control research for wind turbines. Prog Aerosp Sci 46: 1-27. doi: 10.1016/j.paerosci.2009.08.002
![]() |
[188] |
Bossanyi E, Savini B, Iribas M, et al. (2012) Advanced controller research for multi‐MW wind turbines in the UPWIND project. Wind Energy 15: 119-145. doi: 10.1002/we.523
![]() |
[189] |
Márquez FPG, Tobias AM, Pérez JMP, et al. (2012) Condition monitoring of wind turbines: Techniques and methods. Renewable Energy 46: 169-178. doi: 10.1016/j.renene.2012.03.003
![]() |
[190] |
Hameed Z, Hong Y, Cho Y, et al. (2009) Condition monitoring and fault detection of wind turbines and related algorithms: A review. Renewable Sustainable Energy Rev 13: 1-39. doi: 10.1016/j.rser.2007.05.008
![]() |
[191] |
Blaabjerg F, Ma K (2013) Future on power electronics for wind turbine systems. IEEE J Emerging Selected Topics Power Electron 1: 139-152. doi: 10.1109/JESTPE.2013.2275978
![]() |
1. | Nepu Saha, Kyle McGaughy, Michael A. Held, M. Toufiq Reza, Hydrothermal degradation of β-estradiol and oxytetracycline at selective reaction severities, 2020, 2, 2523-3963, 10.1007/s42452-020-03436-0 | |
2. | Dennis Jung, Michael Zimmermann, Andrea Kruse, Hydrothermal Carbonization of Fructose: Growth Mechanism and Kinetic Model, 2018, 6, 2168-0485, 13877, 10.1021/acssuschemeng.8b02118 | |
3. | Camila Centurion Silva, Camila Almeida Melo, Francisco Holanda Soares Junior, Altair Benedito Moreira, Odair Pastor Ferreira, Márcia Cristina Bisinoti, Effect of the reaction medium on the immobilization of nutrients in hydrochars obtained using sugarcane industry residues, 2017, 237, 09608524, 213, 10.1016/j.biortech.2017.04.004 | |
4. | Piyaporn Wataniyakul, Panatpong Boonnoun, Armando T. Quitain, Mitsuru Sasaki, Tetsuya Kida, Navadol Laosiripojana, Artiwan Shotipruk, Preparation of hydrothermal carbon as catalyst support for conversion of biomass to 5-hydroxymethylfurfural, 2018, 104, 15667367, 41, 10.1016/j.catcom.2017.10.014 | |
5. | M Toufiq Reza, Jessica Nover, Benjamin Wirth, Charles J Coronella, Correction: Hydrothermal carbonization of glucose in saline solution: sequestration of nutrients on carbonaceous materials, 2018, 6, 2333-8334, 269, 10.3934/energy.2018.2.269 | |
6. | Kyle McGaughy, M. Toufiq Reza, Recovery of Macro and Micro-Nutrients by Hydrothermal Carbonization of Septage, 2018, 66, 0021-8561, 1854, 10.1021/acs.jafc.7b05667 | |
7. | Sunday E. Elaigwu, Gillian M. Greenway, Characterization of Energy-Rich Hydrochars from Microwave-Assisted Hydrothermal Carbonization of Coconut Shell, 2019, 10, 1877-2641, 1979, 10.1007/s12649-018-0209-x | |
8. | María E. Chiosso, Ileana D. Lick, Mónica L. Casella, Andrea B. Merlo, Acid functionalized carbons as catalyst for glycerol etherification with benzyl alcohol, 2020, 37, 0104-6632, 129, 10.1007/s43153-019-00002-z | |
9. | Magdalini Tsarpali, Jessica Martin, John Kuhn, George P. Philippidis, Valorization of Brassica carinata biomass through conversion to hydrolysate and hydrochar, 2022, 2190-6815, 10.1007/s13399-022-02578-y | |
10. | Erika Kutalkova, Tomas Plachy, Semi-conducting microspheres formed from glucose for semi-active electric field-responsive electrorheological systems, 2022, 18, 1744-683X, 9037, 10.1039/D2SM01145J | |
11. | Meng Xiao, Wenqian Qi, Shuanghui Jia, Mengting Pang, Fengcheng Shi, Hui Mao, High-performance Removal of Tetracycline Enabled by Fe0 Nanoparticles Supported on Carbon@ZIF-8, 2022, 38, 1005-9040, 1349, 10.1007/s40242-022-2255-y | |
12. | Paulo Tarso Ferreira Sales, Katia Maria Souza, Alyne Gonçalves Bezerra, Satu Anneli Ojala, Sérgio Botelho Oliveira, Pierre Alexandre Santos, Maria Teresa Freitas Bara, Chemical Reaction: Understanding the Key to the Formation of Carbonaceous Materials from Sucralose, 2021, 6, 2365-6549, 11846, 10.1002/slct.202102243 | |
13. | Yen-Ju Lai, Jo-Shu Chang, Duu-Jong Lee, Synthesis of a novel solid mediator Z-scheme heterojunction photocatalysis Fe3O4/C/uio66-nh2: Used for oxidation of Rh6G in water, 2023, 231, 00139351, 116264, 10.1016/j.envres.2023.116264 |
Hydrochar | Solid production1 (g per 100 mL) | Mass yield (g g-1dry reactants) | Elemental analysis2 | ICP-OES3 | |||||||
Nitrogen (%) | Carbon (%) | Sulfur (%) | Hydrogen (%) | Oxygen (%)4 | Iron (mg kg-1) | Potassium (mg kg-1) | Phosphorus (mg kg-1) | Sulfur (mg kg-1) | |||
Dry reactants | NA | 1.0 | 1.4 | 23.3 | 2.3 | 4.2 | 68.8 | 55,800 | 117,194 | 31,000 | 23,593 |
HTC-180 | 9.1 ± 0.5 | 0.53 ± 0.04 | 0.9 | 62.7 | 0.2 | 4.5 | 31.7 | 5271 | 3800 | 3424 | 533 |
HTC-230 | 9.5 ± 0.4 | 0.56 ± 0.04 | 1.0 | 64.0 | 0.2 | 4.5 | 30.3 | 8920 | 6389 | 5346 | 790 |
HTC-300 | 11.1 ± 0.5 | 0.65 ± 0.06 | 0.9 | 63.1 | 0.3 | 3.9 | 32.4 | 15,621 | 6309 | 6412 | 1579 |
1 @Solid \; production = \frac{{mass \; of \; hydrochar \; dry \; \left( g \right)}}{{Total \; volume \; of \; feed \; \left( {ml} \right)}} \times 100@ | |||||||||||
2 Precision of the elemental analysis is ±0.5 wt% | |||||||||||
3 Precision of the ICP-OES is 100 mg kg−1 | |||||||||||
4 % O = 100 - (%C + %H + %N + %S); NA: Not applicable |
Morphology | Method | HTC-180 | HTC-240 | H-300 |
Surface area | Multi point BET | 1.85 | 2.05 | 0.98 |
(m2 g-1) | BJH method cumulative adsorption surface area | 7.96 | 8.26 | 3.80 |
BJH method cumulative desorption surface area | 5.47 | 6.44 | 2.96 | |
DH method cumulative adsorption surface area | 8.08 | 8.40 | 3.87 | |
DH method cumulative desorption surface area | 5.65 | 6.68 | 3.07 | |
t-method external surface area | 1.85 | 2.05 | 0.98 | |
t-method micropore surface area | 1.81 | ND | ND | |
NLDFT cumulative surface area | 0.00 | 2.01 | 0.89 | |
Pore volume | BJH method cumulative adsorption pore volume | 1.10×10-2 | 1.11×10-2 | 5.70×10-3 |
(mL g-1) | BJH method cumulative desorption pore volume | 1.11×10-2 | 1.08×10-2 | 5.54×10-3 |
DH method cumulative adsorption pore volume | 1.06×10-2 | 1.08×10-2 | 5.52×10-3 | |
DH method cumulative desorption pore volume | 6.89×10-4 | 1.05×10-2 | 5.38×10-3 | |
t-method micropore volume | ND | ND | ND | |
HK method cumulative pore volume | 6.89×10-4 | 7.65×10-4 | 3.47×10-4 | |
SF method cumulative pore volume | 7.11×10-4 | 7.89×10-4 | 3.60×10-4 | |
NLDFT method cumulative pore volume | 5.45×10-3 | 5.73×10-3 | 1.94×10-3 | |
Pore size | Average pore Radius | 10.71 | 9.23 | 10.05 |
(nm) | BJH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 |
BJH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
DH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 | |
DH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
HK method pore Radius (Mode) | 0.34 | 0.18 | 0.18 | |
SF method pore Radius (Mode) | 0.18 | 0.18 | 0.23 | |
NLDFT pore Radius (Mode) | 20.09 | 20.09 | 0.75 | |
ND: not determined |
Hydrochar | Solid production1 (g per 100 mL) | Mass yield (g g-1dry reactants) | Elemental analysis2 | ICP-OES3 | |||||||
Nitrogen (%) | Carbon (%) | Sulfur (%) | Hydrogen (%) | Oxygen (%)4 | Iron (mg kg-1) | Potassium (mg kg-1) | Phosphorus (mg kg-1) | Sulfur (mg kg-1) | |||
Dry reactants | NA | 1.0 | 1.4 | 23.3 | 2.3 | 4.2 | 68.8 | 55,800 | 117,194 | 31,000 | 23,593 |
HTC-180 | 9.1 ± 0.5 | 0.53 ± 0.04 | 0.9 | 62.7 | 0.2 | 4.5 | 31.7 | 5271 | 3800 | 3424 | 533 |
HTC-230 | 9.5 ± 0.4 | 0.56 ± 0.04 | 1.0 | 64.0 | 0.2 | 4.5 | 30.3 | 8920 | 6389 | 5346 | 790 |
HTC-300 | 11.1 ± 0.5 | 0.65 ± 0.06 | 0.9 | 63.1 | 0.3 | 3.9 | 32.4 | 15,621 | 6309 | 6412 | 1579 |
1 @Solid \; production = \frac{{mass \; of \; hydrochar \; dry \; \left( g \right)}}{{Total \; volume \; of \; feed \; \left( {ml} \right)}} \times 100@ | |||||||||||
2 Precision of the elemental analysis is ±0.5 wt% | |||||||||||
3 Precision of the ICP-OES is 100 mg kg−1 | |||||||||||
4 % O = 100 - (%C + %H + %N + %S); NA: Not applicable |
Morphology | Method | HTC-180 | HTC-240 | H-300 |
Surface area | Multi point BET | 1.85 | 2.05 | 0.98 |
(m2 g-1) | BJH method cumulative adsorption surface area | 7.96 | 8.26 | 3.80 |
BJH method cumulative desorption surface area | 5.47 | 6.44 | 2.96 | |
DH method cumulative adsorption surface area | 8.08 | 8.40 | 3.87 | |
DH method cumulative desorption surface area | 5.65 | 6.68 | 3.07 | |
t-method external surface area | 1.85 | 2.05 | 0.98 | |
t-method micropore surface area | 1.81 | ND | ND | |
NLDFT cumulative surface area | 0.00 | 2.01 | 0.89 | |
Pore volume | BJH method cumulative adsorption pore volume | 1.10×10-2 | 1.11×10-2 | 5.70×10-3 |
(mL g-1) | BJH method cumulative desorption pore volume | 1.11×10-2 | 1.08×10-2 | 5.54×10-3 |
DH method cumulative adsorption pore volume | 1.06×10-2 | 1.08×10-2 | 5.52×10-3 | |
DH method cumulative desorption pore volume | 6.89×10-4 | 1.05×10-2 | 5.38×10-3 | |
t-method micropore volume | ND | ND | ND | |
HK method cumulative pore volume | 6.89×10-4 | 7.65×10-4 | 3.47×10-4 | |
SF method cumulative pore volume | 7.11×10-4 | 7.89×10-4 | 3.60×10-4 | |
NLDFT method cumulative pore volume | 5.45×10-3 | 5.73×10-3 | 1.94×10-3 | |
Pore size | Average pore Radius | 10.71 | 9.23 | 10.05 |
(nm) | BJH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 |
BJH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
DH method adsorption pore Radius (Mode dV(r)) | 0.45 | 0.46 | 0.45 | |
DH method desorption pore Radius (Mode dV(r)) | 0.55 | 0.55 | 0.56 | |
HK method pore Radius (Mode) | 0.34 | 0.18 | 0.18 | |
SF method pore Radius (Mode) | 0.18 | 0.18 | 0.23 | |
NLDFT pore Radius (Mode) | 20.09 | 20.09 | 0.75 | |
ND: not determined |