Nomenclature
- Symbol
Unit
Quantity
-
${c_{pg}}$
-
${{kJ}}/{{kg}} \cdot {{K}}$
Constant-pressure specific heat
-
${D_0}$
mm
Initial diameter of droplet
-
${k_g}$
-
${{W}}/{{m}} \cdot {{K}}$
Thermal conductivity
-
${X_F}$
pixel
Size of the SiC fibre in the image
-
${X_{im}}$
pixel
Corresponding value in the image
-
${\rho _l}$
-
${{kg}}/{{{m}}^3}$
Density
- B
–
Spalding transfer number
- D
mm
Diameter of droplet
- DCN
–
Derived cetane number
- K
-
${{m}}{{{m}}^2}/{{s}}$
Burning rate constant
- t
s
Time
1.0 Introduction
There are four processes of liquid fuel combustion: fuel atomises forming droplets; the initiation of droplet evaporation; formation of a combustible mixture from the evaporated droplet vapor; and the combustion of said mixture [Reference Warnatz, Maas and Dibble1]. Liquid fuel is required to be atomised for fuel to burn efficiency in a combustion chamber. The process of atomisation of the liquid fuel is generally accomplished by injecting the fuel or by mixing the fuel with high pressure air or gas [Reference Chigier2]. Atomised fuel combustion has an exchange of heat, mass, momentum and chemical reactions simultaneously. The obvious complexity of these processes and interactions increases the difficulty of studying atomization combustion [Reference Xin, Ricart and Reitz3–Reference Shen, Li, Liu, Wang and Shen5]. Therefore, droplet combustion is a good way to study liquid fuel combustion [Reference Kotake and Okazaki6, Reference Han, Yang, Zhao, Fu, Ma and Song7]. Rasid and Zhang [Reference Rasid and Zhang8] mentioned that during the evaporation of a neat fuel droplet there are three common combustion phases which are the hearing phase, boiling phase and a disruptive phase. Nomura et al. [Reference Nomura, Takahashi, Suganuma and Kikuchi9] experimentally investigated the mechanism of flame spread along fuel droplet array where they noted that for the multicomponent fuel droplets, the liquid phase process is extremely important to determine the behaviour of the entire fuel droplet’s combustion process [Reference Faik and Zhang10]. There are no chemical reactions that occur in the multicomponent fuel, so each component maintains their own physical and chemical characteristics, such as boiling points and evaporation rates. Both interior heat and mass transfer occur in combustion and evaporation of multicomponent fuel droplets. These characteristics introduce more complexity and uncertainty during the multicomponent fuel droplet combustion processes, as compared to the single component fuel droplet combustion processes.
It has been shown that expansion, ejection and micro-explosion behaviours consistently happen during multicomponent droplet combustion [Reference Lasheras and Dryer11, Reference Chausalkar, Kweon and Michael12]. The occurrence of puffing and micro-explosion directly impacts the overall combustion behaviour of the fuel droplets [Reference Rao and Karmakar13]. Faik and Zhang [Reference Faik and Zhang10] explained different behaviours, such as generation of bubbles, micro-explosion and atomisation of different diesel base multicomponent droplets during the combustion process. Antonov et al. [Reference Antonov, Fedorenko, Strizhak, Nissar and Sazhin14] investigated the expansion and micro-explosion of rapeseed oil/water droplets placed in ethanol as well as a propane/butane mixture flame. They concluded that the amount of time to puffing/micro-explosion of propane/butane mixture flame was always shorter than with an ethanol flame. Faik and Zhang [Reference Faik and Zhang15] studied binary fuel droplet combustion and showed that interacting fuel droplets have higher nucleation rates than with corresponding isolated fuel droplets. Also, the study found that the relationship between nucleation rate and additive concentration is inversely proportional. Meng et al. [Reference Meng, Bao, Li, Wang and Lin16] pointed out the relationship among micro-explosion, temperature and oxygen, as well as combustion efficiency and micro-explosion intensity of different droplets. Due to the water-fuel mixture allowing the reduction of harmful emissions and improving fuel efficiency, many researchers have conducted research on water fuel emulsions [Reference Watanabe, Suzuki, Harada, Matsushita, Aoki and Miura17–Reference Faik, Theeb and Zhang21]. Different models on puffing/micro-explosion, heating and evaporation of water-fuel emulsion droplets have been developed [Reference Nissar, Rybdylova, Sazhin, Heikal, Aziz and Ismael22–Reference Sazhin, Rybdylova, Crua, Heikal, Ismael, Nissar and Aziz24].
Conventional jet fuels have between 20% and 25% concentration of aromatics by mass. Various engine parts require aromatics in the jet fuels to achieve seals compatibility [Reference Almohammadi, Singh, Sharma, Kumar and Khandelwal25, Reference Undavalli, Gbadamosi Olatunde, Boylu, Wei, Haeker, Hamilton and Khandelwal26], and required energy density per unit volume. Although aromatics are an important part of the fuel system, they are also responsible for majority of the PM emissions from liquid fuel combustion [Reference DeWitt, Corporan, Graham and Minus27, Reference Zheng, Singh, Cronly, Ubogu, Ahmed, Ling, Zhang and Khandelwal28]. Different aromatic species typically produce different amounts of PM emissions, so it is important to complete a comprehensive analysis of different aromatics to achieve both low PM emissions and fuel system component compatibility [Reference Sharma, Singh, Bhardwaj, Khandelwal and Kumar29]. Research has been conducted that specifically focused on the aromatic fuel combustion and sprays. Singh et al. [Reference Singh, Sharma, Almohammadi, Khandelwal and Kumar30] studied the emissions of dearomatised hydrocarbon solvent mixed with particular type of aromatic species. Sabourin et al. [Reference Sabourin, Boteler and Hallett31] experimentally studied auto ignition of n-paraffins and n-alkyl aromatics mixtures, then designed a continuous mixture droplet ignition model based on an existing droplet evaporation model. Wijesinghe and Khandelwal [Reference Wijesinghe and Khandelwal32] investigated spray and atomisation characteristics of different aromatic species mixed with a base blend. Moriue et al. [Reference Moriue, Eigenbrod, Rath, Sato, Okai, Tsue and Kono33] experimentally studied the spontaneous ignition of n-decane dilution by aromatic hydrocarbon droplets and suggested volatility of each component dependence should be considered when applied to spray ignition problems. It can be observed from the above discussion that there are studies in literature that have investigated the impact of difference in fuel properties on droplet combustion. Though how different aromatic species burn in droplet combustion is still an area which needs further investigation.
Contaminations have significant impact on the fuel droplet combustion. The use of bio-oil as an additive helps to generate slightly decomposed substances that initiate polymerisation and form residues [Reference Wu and Yang34]. Rasid and Zhang [Reference Rasid and Zhang35] conducted a study on combustion and ignition of soot contamination on the surface of a diesel fuel droplet and concluded that soot particles on the droplet surface can significantly reduce the evaporation rate and disrupt droplet behaviour throughout its lifetime. Following the investigation on combustion characteristics of soot contamination on the surface of diesel fuel droplet, Rasid and Zhang [Reference Abdul Rasid and Zhang36] also studied the combustion behaviour of homogeneous distribution of soot particles in the diesel fuel droplet. They found that the burning rate of contaminated droplets is smaller than neat droplets, and that burning rates can be slightly improve when the loading of soot particles is 0.2% by mass. Polymer addition to fuel has been studied by different researchers, and its impact on droplet combustion has been preliminarily explored. Wang et al. [Reference Wang, Liu and Law37] found that polymer addition in a fuel changed the fuel’s characteristic to a multicomponent fuel. Singh et al. [Reference Singh, Esmaeilpour and Ratner38] experimentally evaluated the impact of carbon-based nanoparticles in crude oil droplets. Results show that the addition of carbon-based nanoparticles can decrease the vaper pressure of crude oil droplets but increase bulk heat conductivity and radiation absorption. Some researchers have shown that the addition of polymers to hydrocarbon fuels imparts non-Newtonian properties to the fuel, thereby suppressing splash behaviour when spilling over surfaces [Reference Ghamari and Ratner39, Reference David, Wei, Liu, Bathel, Plog, Ratner and Kornfield40]. Hazart et al. [Reference Hazrat, Rasul, Khan, Ashwath and Rufford41] experimentally studied the emissions produced by combustion of polymer additive biodiesel fuel and it was shown that CO emissions were reduced when polymer additive biofuel was used as compared to no additive case. It was also said that polymer blended fuels might be used as a new clean fuel source. David et al. [Reference David, Wei, Liu, Bathel, Plog, Ratner and Kornfield42] investigated the effects of mist control of aviation fuel with polymeric additive and found that associative polymers could improve antimisting. Ghamari and Ranter [Reference Ghamari and Ratner39] found that at early stages of fuel droplet combustion, polymer additive could reduce evaporation rate, increase swelling behaviours, which eventually reduced the burning time of diesel and jet fuel droplets. Investigations into fuels with polymer addition has been done, although there is a lack of understanding how these fuels will perform under droplet combustion.
Liquid fuel is commonly stored in different types of containers. Plastic drums are suitable for a wide range of materials and are low cost, so they are one of the most versatile containers used for industrial storage and transportation needs. One of the biggest problems is that if liquid fuel is stored in plastic drums, plastic could react with fuel and break down. The polymer-based plastic in the plastic drums readily interact with the fuel. So, fuel storage in plastic or steel drums with polymer lining may change its characteristics by adding polymeric content from the containers into the fuel. All of the above discussed studies on polymer additive have focused on the intentional addition of polymers to fuels. There is no research in literature about impact of polymer contaminated fuel by storage in plastic drums. One of the main objectives of present study is to investigate the impact of polymer contaminated alternative fuels, including pure aromatics and other various multicomponent fuels on droplet evaporation during the combustion.
2.0 Methodology
2.1 Droplet lifetime and combustion phases
The droplet evaporation and combustion model used in this analysis is the classical theory (also known as D2-law) of droplet combustion designed by Godsave [Reference Godsave43] and Spalding [Reference Spalding44] in 1953 [Reference Choi, Dryer and Ross45]. The theory exemplifies the essential physics and gives an elementary speculation on droplet regression rate [Reference Law46] and is still widely used in various research on droplet combustion [Reference Choi, Dryer and Ross45].
In this study, measurement of the burning rate of a fuel droplet is in accordance with D 2-law, which stated that a squared diameter of fuel droplet through its lifetime produced a linear slope that is measured as burning rate constant, K (
${\rm{m}}{{\rm{m}}^2}/{\rm{s}}$
). K is calculated by the terms of fuel properties and the transfer number, expressed as [Reference Law46]:
\begin{align}K = \frac{{8{k_g}}}{{{\rho _l}{c_{pg}}}}\ln \!\left[ {1 + B} \right]\end{align}
Where,
${k_g}$
is Thermal conductivity (
${\rm{W}}/{\rm{m}} \cdot {\rm{k}}$
);
${\rho _l}$
is Density (
${\rm{kg}}/{{\rm{m}}^3}$
);
${c_{pg}}$
is Constant-pressure specific heat (
${\rm{kJ}}/{\rm{kg}} \cdot {\rm{K}}$
); B is the Spalding transfer number.
For the steady burning state, K is a constant to represent the slope of the droplet size variation with time. Hence, the general equation of the temporal droplet size variation (D 2-law) is expressed as:
\begin{align}{D^2}\!\left( t \right) = {D_0}^2 - Kt\end{align}
Where, D is the diameter of droplet,
${D_0}$
is the initial diameter of droplet; t is time.
Figure 1(a) shows the original image of droplet and fibre while Fig. 1(b) shows the morphological binary image of the droplet. The image has been improved to make it easier determine boundary line and extract the size features of droplets more accurately. Transferring the color image into the morphological mode has been done through a series of the image processing steps. The calibration of mapping between digital and real dimensions has been done before the experiment. The real size of any droplet is calibrated in accordance with the hanging fibre size, shown as Equation (3).
\begin{align}D = 100 \kern1pt \ast {10^{ - 3}}\frac{{{X_{im}}}}{{{X_F}}}\end{align}
(a) The original image of droplet and fibre with fixed size; (b) Morphological binary image of droplet.
Where: D is the real value of the droplet (mm);
${X_{im}}$
is the corresponding value in the image (pixel);
${X_F}$
is the size of the SiC fibre in the image (pixel).
The actual droplet diameter has been converted as per Equation (3) with the reference to the scale. SiC fibre of known size (
${10^{ - 3}}{\rm{mm}}$
) is used as reference scale to determine the value of
${X_F}$
, shown in Fig. 1(a). In this experiment, the average value of the measured transverse, and longitudinal data is taken as the calculated diameter of the droplet (
${X_{im}}$
). For example, the resolution of the droplet is shown in Fig. 1(b), there are 164.50 pixels in the horizontal direction and 164.13 pixels in the vertical direction, the average pixels number
${X_{im}}$
is 164.32. By calculation, spatial resolution is 10
${\rm{\mu m}}/{\rm{pixel}}$
. Therefore, the actual size of the droplet (D) in this example is 1.6432mm. All the camera setting and spatial resolution data has been listed in Table 1.
Several combustion phases were observed throughout the lifetime of the combusted fuel droplets. Specifically, a binary fuel mixture will have a few additional phases of combustion behaviour as shown in Fig. 2. To facilitate the comparison of different droplet sizes, based on the D 2-law analysis, the droplet size is normalised for the initial droplet size
${(D/{D_0})^2}$
and the time is normalised to the initial droplet size squared
$t/{({D_0})^2}$
.
It can be observed in Fig. 2 that the squared droplet diameter is continuously reduced through the passing of time while the slope of the plot changes. Therefore, according to the droplet burning stage, after ignition, the droplet combustion state is divided into five phases (P1 to P5). The specific combustion states and the associated behaviours are shown in Table 2.
Camera setting and spatial resolution data
Binary fuel droplet phases throughout its lifetime during combustion process based on D2-law analysis.
Description of combustion phases in the combustion process of binary fuel droplet
2.2 Experiment setup
A Photron-SA4 high-speed color camera with Nikon AF Micro NIKKOR 60mm f/2.8D lens was used for the tracking of the lifetime and dynamics of the droplet. Backlighting was incorporated in this experiment to track the size and sharp edge data of the transparent droplets. The backlit imaging was performed by placing a high-intensity light source (IDT 19-LED) with a light diffuser behind the droplet. The images for this study were recorded with a frame rate of 40,000fps. Heated thermal Kanthal wire was placed 1mm below the surface of the liquid drop, and the surrounding ambient air is rapidly heated to 800°C. At this time the droplet exceeds the ignition temperature where it can be ignited without smoke contamination. When the droplet flame is formed, the hot filament is removed quickly to avoid affecting the normal combustion state of the droplet. The overall experimental arrangement is shown in Fig. 3.
Experimental setup in high-speed imaging of fuel droplet combustion.
For each test, the initial diameters of droplets are kept at 1mm with a measured standard deviation of 0.05mm. The images (in TIFF format) were recorded between the activation of the ignition and the extinguishing of the flame. MATLAB was used to process these images.
2.3 Selected fuels and properties
To determine the impact of polymer contamination on droplet combustion, a number of fuels were selected, and the fuel details are listed in Table 3. Fuels are separated into two categories based on their storage types: plastic or steel. Both sets of drums were stored indoors in dry and cool conditions with no exposure to sunlight for three months prior to testing. Actual pictures of fuel drums storage are shown in Fig. 4. The plastic drums used in this study were being sold on the market as suitable for fuel storage. The material of the plastic drums is HDPE, and it was rated for UN Rated 3H1. The steel drums used in this study were lined and approved for fuel storage.
Fuel stored in steel drum and fuel stored in plastic drum
Test fuel batches: plastic (left), steel (right).
Table 4 shows the selected aromatics used in this experiment. The aim is to determine the burning rate constant of selected aromatics droplet. Table 5 shows various multicomponent and cetane number of Jet A fuels used in this experiment, in purpose of measuring their burning rates.
3.0 Results and discussion
3.1 Droplet combustion characteristics of fuel stored in steel and plastic drum
Figure 5 shows the evolution of diameter squared of fuel stored in steel drum and fuel stored in plastic drum. In each multicomponent and contaminated fuel, there is a distinct phase transition from high volatility (P2) to low volatility (P4) phase. This effect is clearly illustrated in the contaminated Jet A fuel sample, which produced two slopes of combustion phases – with a high volatility and low volatility combustion phases. In P2 stage, there is no obvious difference between the fuels stored in the two containers; however, in P4 stage, the burning lifetime of the fuels stored in the plastic containers can be prolonged significantly. Especially in AR-series fuel products, the higher the mixing degree of Toluene, the more obvious the difference of droplet size. However, Jet A did not show significant difference when stored in plastic drums as compared to steel drums. It indicates that some fuel components are more sensitive to contamination in plastics.
Aromatics fuel samples
Various multicomponent fuels and Jet A with different cetane number
Evolution of diameter square of fuel stored in steel drum and fuel stored in plastic drum for (a) AR-1A, (b) AR-1B, (c) AR-1C and (d) A2.
Figures 6 and 7 shows the burning rate constant of tested fuels in high volatility (P2) and low volatility (P4) combustion phases, respectively. It can be observed from Fig. 6 that burning rate constant has changed in high volatility (P2) phase when the fuel is stored in plastic drums as compared to steel drums. For the steel drums, with increasing content of Toluene the burning rate constant at high volatility (P2) phase has decreased, and for plastic drum case it has increased. This change in P2 stage could be due to the fact that Toluene is an aromatic hydrocarbon, and its molecular bond is relatively stable compared with Bannersol during combustion. Therefore, theoretically, with the increase of T content, the combustion rate should be lower. The combustion rate of AR series fuels stored in steel drums shows a decreasing trend, which basically conforms to the normal theoretical cognition. The fuel in the plastic drums shows a completely opposite trend to that expected. Therefore, it indicates that impurities in the plastic drum may change its chemical properties and cause this unexplained result.
Burning rate constant, K in high volatility combustion phase (P2) of fuel stored in steel drum and fuel stored in plastic drum. The error bars show standard deviation.
Burning rate constant, K in low volatility combustion phase (P4) of fuel stored in steel drum and fuel stored in plastic drum. The error bars show standard deviation.
Liquid phases of fuel stored in steel drum and fuel stored in plastic drum in low volatility combustion phase.
Evolution of diameter square of pure aromatics and bannersol.
Different aromatic species droplet phases throughout its lifetime during combustion process based on D2-law analysis.
Burning rate of pure aromatics and bannersol. Error bars show standard deviation.
Evolution of squared diameter of various multicomponent fuels and cetane number of Jet A.
Burning rate of various multicomponent and cetane number of Jet A. The error bars show standard deviation.
At the P4 stage, in general, burning rate constant in P4 is higher than in P2. Regardless of the term ‘low volatility’, higher burning rate in these phases were due to its higher boiling point of the components, thus increased the droplet temperature and evaporation rate in the process. All fuels follow the rule that the combustion rate decreases with the increase of aromatic hydrocarbons. In addition, the uncertainty is also reduced. This indicates that with the increase of droplet temperature, the impact of burning rate caused by pollution from plastic drums is alleviated. This means that if burning in the engine, the impact of plastic drum pollution on the fuel may not be obvious at high temperature or rich fuel combustion stage, but it will have a greater impact at low temperature or lean fuel burning and ignition.
Test done by Ghamari and Ratner [Reference Ghamari and Ratner39] showed that polymer contamination in Jet-A has low effect on the burning rate reduction in volatile phase combustion compared to diesel fuel. This shows that changes in the burning rate of the fuel is highly dependent on the type of fuel contaminated by polymer additives. From Figs. 6 and 7, it can be observed that Toluene has a stronger dependence on polymer contamination, as fuel with higher Toluene composition has increased the burning rate in high volatility phase.
This justified the result obtained in this test, which has proven that there was an increase in burning rate for contaminated fuels of AR-1A, AR-1b and A2 fuel, higher than their respective fine fuels. During this phase, the composition of polymer contamination was denser, thus increased the evaporation rate due to improved heat conductivity within the fuel droplet.
Polymer has a higher density than both diesel and jet fuel. It also has very low volatility because of a very high boiling point. As droplets continue to evaporate, the mass fraction of polymer increases, which lead to increasing droplet surface temperature to exceed boiling temperature of fuel. Under this condition, polymer will cause the droplets to swell and bubble. An interesting phenomenon revealed by Yang et al. [Reference Yang, Hamins and Donnelly47] is, during combustion, most of the polymer sphere is made by gaseous bubbles, then bubbles will grow and escape from the polymer, causing sputter.
Figure 8 shows the images of liquid phases of fuel stored in steel drum and fuel stored in plastic drum in low volatility combustion phase. Referring to liquid state of fuel stored in plastic drum in Fig. 8, it was observed that polymer presence in the fuel played as the nucleation sites that created bubbles emerging from inside the droplet. More active swelling of the droplet was recorded for each contaminated fuel, as higher droplet temperature tends to increase the core temperature and thus increased the fuel temperature, exceeding its boiling temperature and resulting a puffing condition. Referring to the contaminated fuels liquid state in Fig. 8, it was observed that polymer presence in the fuel played as the nucleation sites that created bubbles emerging from inside the droplet.
As reported by Ghamari and Ratner [Reference Ghamari and Ratner39], soot produced during combustion might be captured by the polymer strings in polymer additive fuel and explains the agglomerations shown in contaminated fuels in Fig. 8. Higher amounts of soot trapped inside the droplet further increased the boiling temperature, thus elevated the fuel temperature beyond its boiling temperature, causing micro shootings of the soot particle and fuel sub-droplets.
3.2 Droplet combustion characteristics of selected aromatics
Figure 9 shows the structure of different aromatic species selected in this study for investigation on droplet combustion characteristics. Figure 10 shows different aromatic species droplet phases throughout its lifetime during combustion process based on D2-law analysis. It can be observed from Figs. 9 and 10 that the biphenyl benzene ring with larger molecules, such as the droplet of Tetralin, is the most difficult to decompose and often burns the longest. On the contrary, some fuels with smaller molecules (branched molecule), such as Toluene and Terbutylbenzene, rapidly decrease in droplet size and burn out. This is an indication that in combustion of aromatics droplet, low volatility phase burns with higher rate due to the breakup of strong molecular chains, thus evaporated the fuel even faster.
Figure 11 compares the combustion rates of different fuels in different combustion stages (P2 and P4). As expected, the burning rate constant, from high to low-rate, correlates with aromatics that possess the lowest boiling point to the highest boiling point. As the fuel with the smallest molecule, Toluene’s combustion rate is obviously higher than other fuels in both P2 and P4 stages. The reducing trend of burning rates are obvious in low volatility combustion phase (P4, as they go through chemical compositions). O-Xylene has shown significantly different combustion properties. In P2 stage, O-Xylene combustion rate is the lowest, which even lower than that of bicyclic aromatic hydrocarbons (Tetralin). However, in the P4 stage, its combustion rate is only second to toluene. In terms of the complexity of the molecular chain, its combustion performance should be in the middle position, so this study guesses that in addition to the factors of volatility and molecular arrangement, there are other physical and chemical characteristics that need to be considered, and further research needs to be carried out.
3.3 Droplet combustion characteristics of various multicomponent fuels
Figure 12 shows droplet evolution of various fuel mixtures. From the chemical composition table, C8 contained aromatics blend. Obvious irregular change in droplet diameter was observed in its combustion behaviour; having a chemical decomposition phase as illustrated in the Fig. 12. Theoretically, for Jet A with distinct DCN values, a rise in cetane number should correspond to an increased combustion rate. However, in this particular experiment, Fig. 13 reveals that no clear trend changes were detected in relation to shifts in DCN.
4.0 Conclusion
Droplet evaporation during the combustion of different alternative and conventional fuels with a range properties and storage conditions have been experimentally investigated in this paper. Selected fuels are divided into three groups. In the first group there are four different fuels. They are stored in plastic drums (approved for fuel storage) and steel drums respectively to determine the impact of storage drums. The second group is made up of seven different pure aromatic species. The third group contains various multicomponent refined fuels and Jet A with different cetane number. For fuels stored in plastic and steel drums, it was observed that combustion properties change when fuel was stored in plastic drums. Plastic drums were found to be contaminating the fuel and it is found that they are more likely to produce unstable phenomena such as soot, jet and bubble, which leads to fuel consumption but not complete combustion. It was also found that the impact of plastic drums contamination on fuel may not be obvious at high temperature combustion or during the rich fuel combustion stage, but it will have greater impact at low temperature or lean fuel combustion and ignition stages. For the pure aromatics, molecular structure, atomic bond stability and volatility jointly impact the combustion rate of droplets. It was observed that smaller molecules and highly volatile fuels burn faster. However, O-Xylene showed different burn rates trends as compared to other aromatic species. Therefore, the influence of other chemical and physical factors is not excluded. For the various multicomponent fuels and Jet A with different DCN, theoretically, with the increase of cetane number, the combustion rate should show an upward trend. However, in this experiment, no obvious trend change with the change of DCN has been observed.
Open access
