Next Article in Journal
Response Surface Methods Used for Optimization of Abrasive Waterjet Machining of the Stainless Steel X2 CrNiMo 17-12-2
Previous Article in Journal
Low-Frequency Dielectric Relaxation in Structures Based on Macroporous Silicon with Meso-Macroporous Skin-Layer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Rheological Analysis of Biomaterial Behaviour as a Tool to Detect the Dilution of Heather Honey

1
Department of Zoology, Fisheries, Hydrobiology and Apidology, Faculty of AgriSciences, Mendel University in Brno, Zemedelska 1, 613 00 Brno, Czech Republic
2
Department of Agricultural, Food and Environmental Engineering, Faculty of AgriSciences, Mendel University in Brno, Zemedelska 1, 613 00 Brno, Czech Republic
3
Department of Food Technology, Faculty of AgriSciences, Mendel University in Brno, Zemedelska 1, 613 00 Brno, Czech Republic
4
Department of Technology and Automobile Transport (Section Physics), Faculty of AgriSciences, Mendel University in Brno, Zemedelska 1, 613 00 Brno, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 15 April 2021 / Revised: 4 May 2021 / Accepted: 6 May 2021 / Published: 11 May 2021

Abstract

:
Heather honey is a valuable and rheologically special type of honey. Its above-average selling price may motivate its intentional violation with a mixture of honey from another botanical origin, the price of which is lower on the market. This work deals with the rheological properties of such devalued heather honey in order to determine the changes in the individual rheological parameters depending on the degree of dilution of the heather honey. For this purpose, a differently diluted heather honey sample series was created and the following rheological parameters were determined: hysteresis area, n-value, yield stress (τ0), parameter B (Weltman model), parameter ϕ, or parameter C (model describing the logarithmic dependence of the complex viscosity on the angular frequency). Part of the work was research into whether the set parameters can be used as comparative parameters. It was found that the hysteresis area does not appear to be a suitable relative comparison parameter due to the high variability. The parameters that appear to be suitable are the relative parameters n-value and the parameter ϕ, which showed the greatest stability. The change in the determined rheological parameters is, depending on the degree of dilution, non-linear with a step change between the samples containing 40% (w/w) and 60% (w/w) of a heather honey.

1. Introduction

Honey is a natural food product and is considered to be the oldest sweetener. It is a very good source of many valuable nutrients that are beneficial to the human body. Honey has antioxidant, antibacterial, antifungal and antiviral properties, which are also used in the so-called apitherapy [1]. The antibacterial activity of honey is caused by its low acidity, high osmolarity and the content hydrogen peroxide and other components, such as polyphenolic compounds [2]. Ling heather honey (hereinafter “heather honey”) is a special type of honey arising from nectar of Calluna vulgaris. This honey is highly valued for its characteristic strong taste, health properties and unusual texture characteristics. It is a low-growing shrub predominantly on acid soils in the sun or semi-shade and blooms in late summer [3]. This plant species has a high beekeeping value. Honeybees from its nectar produce honey with special sensory properties, such as an amber to dark colour, a plant fragrance reminiscent of a moist ground and a persistent taste of bitterness, which are highly appreciated by consumers [4].
The price of heather honey on the market is considerably higher than common honey for reasons of it’s the exceptionality and difficulty in harvesting. This also applies to other special honeys, such as manuka honey (Leptospermum scoparium), which is one of the highly priced honeys [5]. Moreover, climatic change is bound to alter the conditions in heathlands and increase the risk of a reduced honey production [6]. It can, therefore, be assumed that these special honeys can be the target of fraud. Artificial dilutions are used for the violation of common honeys [7]. Alternatively, a common honey can be used for dilution of the special honeys. There are standards where the parameters are stated which define special honeys [8]. Various analytical methods can be used to determine these parameters and characterise the honey [9]. The determination of the rheological properties can be considered one of the indirect methods [10,11]. The fact that some special honeys show a different rheological behaviour, compared to common honeys, is used as an indicator in this area.
It is well known that some honey types in dependence on the botanical origin were classified as non-Newtonian fluids with a thixotropy behaviour [12]. Heather honey is one of them [1,13]. There is a presumption that, in the case of diluting heather honey with common honey, there will be a change in the rheological properties. The rheological behaviour of a heather honey mixed with a common honey has not been studied. Determination of the rheological properties of such mixtures can be a support method for the detection of a honey’s adulteration to protect consumers.
There are many works dealing with the problem of honey adulteration [7,14,15,16,17] and less works, where the knowledge of rheology is used for this purpose [11,18,19,20,21,22,23].
This study describes rheological properties and their changes in the heather honey in dependence of its experimental dilution by a honey with the Newtonian properties.

2. Materials and Methods

2.1. Samples

The lime (Tilia spp.) and ling heather honey (Calluna vulgaris, hereinafter “heather honey”) were used. The lime honey was produced at apiary of the first author (Příbram na Moravě, Czech Republic) and the ling heather honey was processed in the company Honningcentralen (Kløfta, Norway). Initially, both samples of original honey were de-crystallised in a thermostat and heated up to 50 °C. The warming-up lasted until the crystals in the honey were melted. Samples of the above and below defined original honeys (Table 1) as well as the mixed samples were prepared for the rheological tests. The melted honeys were mixed in a weight concentration with an accuracy 0.01 mg. The mixed samples contained 80% (S80), 60% (S60), 40% (S40), 20% (S20) and 10% (w/w) (S10) of the heather honey. The sample identification for pure lime honey is S0 and for pure heather honey is it S100.

2.2. Characterisation of Honey

The water content in the samples was measured by a refractive index at a temperature of 20–22 °C with extrapolation to 20 °C (ABBE refractometer, Carl Zeiss, type G, Jena, Germany) [24]. The electrolytic conductivity was measured with a Mitronic MVM-1 device at room temperature of 22 °C with extrapolation to 20 °C according to the correction factor described by [24].
The botanical origins of the honey samples were evaluated by melissopalynological analysis, in accordance with the harmonised methodology developed by [25] and other supporting protocols by [8,26,27]. The pollen grains were identified with the pollen collection at Mendel University in Brno, and the pollen grain catalogue by [28]. One tetrad of heather pollen grains was counted as a one pollen grain. The sensory properties of the samples (colour, odour and taste) corresponded with results of the pollen analysis. The geographical origin of the honey samples was not necessarily verified due to the honey samples coming from known apiaries. The honeydew elements were evaluated subjectively as follows: (1) The abundance of clusters, and (2) the abundance of fragments of algae and fungi. Their abundance was determined on the following subjective scale: Absent, sporadic, few, moderate and widely abundant.
The water content and the electrolytic conductivity of the honeys were determined in accordance with the principles of the harmonised methods of the European Commission for honey [9]. The standard-deviation was calculated for every honey parameter from three replicate assays.

2.3. Rheological Assays

The rheological evaluation of the honey in this paper was performed on an MCR 102 rheometer (Anton Paar, Graz, Austria) with a measuring geometry cone-plate. The gap between the cone and the plate was set to a stable value of 0.103 mm. The diameter of the cone was equal to 50 mm with an angle of 1°. The rheological tests were performed at a temperature of 25 °C. Each measurement was repeated four times.

2.3.1. Hysteresis Loop Test

The range of the shear rate was set from 1 to 100 s−1 [29]. The hysteresis area was calculated with the use the corporate software RHEOPLUS v. 3.61 (Anton Paar, Austria).

2.3.2. The Determination of Time-Dependent Behaviour

A constant value of 50 s−1 of the shear rate was set [30]. The values of the apparent viscosity were recorded for 300 s. Subsequently, the proportions of the apparent viscosity in the 300th and in the 1st second of the test were carried out. The time dependence of the shear stress versus time at a constant shear rate (50 s−1) was carried out. The dependence was described by the Weltman model [31]:
τ = A B ( ln t )
where:
  • τ —shear stress (Pa)
  • A —value of the shear stress in the first second of the measurement (Pa)
  • B —time coefficient of the thixotropic breakdown (-)
  • t —time of shearing (s)

2.3.3. Dynamic Frequency Sweep Test

The dynamic frequency sweep test was carried out in a frequency range of 0.1–10 Hz at a temperature of 25 °C and at a 1% strain. It was determined by an amplitude sweep test. The amplitude sweep test was determined at 1 Hz in a strain range of 0.1–100%.

2.3.4. Mathematical Modelling

The dependence of the shear stress on the shear rate was described by two non-Newtonian mathematical models. The mathematical models were determined with use MATLAB R2018b (MathWorks, Natick, MA, USA). The first model was the power-law fluid (also known as the Ostwald-de Waele model). The Ostwald-de Waele model is given by following equation [32]:
τ = K γ ˙ n
The second mathematical model was the Herschel-Bulkley rheological model [33]:
τ = τ 0 + K γ ˙ n
where:
  • τ —shear stress (Pa)
  • τ 0 —yield stress (Pa)
  • K —consistency factor (Pa·sn)
  • γ ˙ —shear rate (s−1)
  • n —flow index (-)

2.4. Statistical Analysis

The open source software R (v. 3.5.1) [34] was used for the statistical analysis. All the tests were performed at a significance level of α = 0.05. The Shapiro–Wilk test was used to test the hypothesis whether the data are from a normal distribution. A one-factor ANOVA (analysis of variance) was used to test the hypothesis whether there are significant differences between the individual characteristics of a data set. Tukey’s test was used for the post hoc analysis. However, when the hypothesis h0 about the homogeneity of variance was rejected, Nemenyi’s test was used. The hypothesis h0 about the homogeneity of variance was tested with use of Leven’s test.

3. Results

3.1. Hysteresis Area

In the first phase, the dependence of the shear stress on the shear rate was measured. Subsequently, the hysteresis areas were calculated. The course of curve showing the dependence of the shear stress on the shear rate for sample S80 and S100 is obvious from Figure 1. Other samples (S0, S10, S40 and S60) are given in Figure 2. The value of the shear stress at 100 s−1 is higher in sample S80 than the shear stress value of sample S100.
The exploratory data analysis of the hysteresis area values for the individual samples is shown in Table 2. The average values of the hysteresis loop are also shown in the Figure 3 for clarity. The sample S0 showed a hysteresis area value despite that it is demonstrably a pure lime nectar honey. The values of a hysteresis area are negative for the samples S0, S10. The value of a hysteresis area grows non-linearly as a proportion of the heather honey increases in the sample. A sudden increase of a hysteresis area value can be seen between samples S20 and S40 (see Figure 3). The values of hysteresis areas showed a high degree of variability, as can be seen from the Table 2. The sample S20 showed the highest value of a variation coefficient. Contrarily, the sample S100 showed the lowest value of a variation coefficient.
The datasets containing hysteresis area values of all the samples were compared with each other by use of an ANOVA analysis and Tukey’s test for the post hoc analysis. The results are shown in Appendix A (Table A1). The hypothesis that the differences between sample S60 and all other samples, sample S80 and all other samples, and S100 and all other samples are not significant was rejected.
In the next step, the mathematical description of the shear stress and shear rate dependence was carried out. For this purpose, the Ostwald-de Waele (OW) and Herschel-Bulkley (HB) models were used. The results of the mathematical modelling are given in Table 3. The flow index n , consistency factor K , possible yield stress τ0 were determined. In contrast to the hysteresis area values, the other parameters ( n , K and τ 0 ) changed adequately with the change in the honey dilution.
The flow index n determines the degree of curvature. This curvature (including the change in the curvature) can also be seen in Figure 1 and Figure 2. The factor n of sample S0 is very close to one in the case of both models. This means that the function describing the dependence of the shear stress on the shear rate is linear and this is a property of Newtonian fluids. As the amount of heather honey increases, the value of n gradually decreases. The values of n < 1 are typical for shear-thinning materials. In the case of both models, the decrease in the value of the change in index n is, at first, gradual. Between samples S40 and S60, the change increases significantly (see Figure 4). The index of change in index n between samples S40 and S60 for the OW model is the lowest ( i = 0.86), i.e., the change from the previous value is the highest in the data file. In the case of the HB model, the index of change between samples S40 and S60 is i = 0.86 and it is also the highest change in the data file. It can also be observed from Table 3 that the standard deviation of a parameter n is relatively small (despite the very high variability of the hysteresis area values). The maximum value of the coefficient of variation was 2.8% for sample S100.
The change in the consistency factor ( K ) is also nonlinear. The largest change occurs, as with index n , between samples S40 and S60. For the OW model, the index of change is i = 2.52, for the HB model, it is i = 2.26.
In the HB model, negative yield stress values ( τ 0 ) were calculated in some cases. For samples S0–S40, the values of τ 0 are around the value of 0 Pa. Again, there is a step change between samples S40 and S60.

3.2. Dependence of the Apparent Viscosity versus Time

The results of the measurement showing the dependence of the apparent viscosity versus time are graphically presented in Figure 5. In the case of samples S80 and S100, it can be observed that in the first seconds of the measurements, high apparent viscosity values are measured. Subsequently, a sharp decline occurs. Such behaviour is attributed to shear-thinning fluids with thixotropic behaviour. With that, the values of the apparent viscosity in the case of sample S80 are generally higher than the apparent viscosity values of sample S100. It is similar to the case of the values of the hysteresis area. However, if the proportion of the apparent viscosity value in the 300th second of the measurement and the apparent viscosity value in the 1st second of the measurement is set, the situation is changes. This proportion can be marked by the Greek letter ϕ . The following equation for the calculation of ϕ can be introduced:
ϕ = η A 1 η A 300
where:
  • η A 1 —apparent viscosity in the 1st second of the assay (Pa·s)
  • η A 300 —apparent viscosity in the 300th second of the assay (Pa·s)
The results of these data are evident from Table 4. As in the case of parameter n , the change of parameter ϕ is, at first, gradual, then the largest change occurs between samples S40 and S60 (see Figure 6). The change is 13% here (the second highest change is 9%). The Low data variability is also evident, where the value of the variation coefficient ranges from 0.4 to 5.8%.
The hypothesis whether there are significant differences between ϕ was tested. The results are given in Table A2 in Appendix A. A significant difference was found between the ϕ value of sample S100 and all the other samples. The exception is sample S80, where the hypothesis h0 was not rejected.
The determination of the thixotropic index, using the Weltman model, is necessary for the determination of thixotropy degree. Results are given in Table 5. Coefficient B (time index of thixotropy) is the most important parameter of the Weltman model. The negative value means that the curve grows over time, a positive value means the opposite. It is possible to observe the trend in the results of coefficient B . The significant change again occurs between samples S40 and S60 (see Figure 7). The highest coefficient B value can be seen in the case of sample S100. However, it should be noted that a high degree of variability is possible in the data sets. The greatest variability of data was found in the sample S40, where SD = 2.2 and the coefficient of variation was 96% (see Table 5). It is also clear from the results that the highest values of the coefficients of determination are achieved in sample models S0 and S100, i.e., the original unmixed honey solutions. A significant difference was calculated (ANOVA, Tukey HSD) between sample S100 and samples S0, S10, S20, S40 (see Appendix A, Table A3).

3.3. Frequency Sweep Test

The results of frequency sweep test are evident from Figure 8 and Figure 9. The first figure represents the dependence of the loss modulus ( G ) and storage modulus ( G ) on the angular frequency in samples S100 and S60. The second figure represents the dependence of the complex viscosity ( η * ) on the angular frequency of S0, S40, S60, S80 and S100.
The results of the frequency sweep test show that the curve storage modulus is gradually approaching the curve of the loss modulus as the proportion of heather honey in the sample increases. In other words, the curves are closest in the case of sample S100. On the contrary, the largest gap can be seen in the case of S0 sample. In the all the cases, the viscous module is dominant.
It is possible to describe the dependence of the complex viscosity ( η * ) on the angular frequency by a logarithmic function. This function can be represented by the following equation:
η * = C ( ln ω ) + D
where:
  • ω —the angular frequency (rad·s−1)
  • C —the degree of a declining complex viscosity
  • D —the initial value of the complex viscosity at 1 rad·s−1
It can be assumed that the higher the slope of the curve, the better the system is structured. Additionally, otherwise, the higher the negative value of the C parameter, the system is better structured. The change of parameter C depending on the sample is evident from the Figure 10. Taking the variability of the data into account, the highest negative value of parameter C is seen in the case of samples S80 and S100 (see Table 6). Significant differences were recorded between samples S100 and the other samples (ANOVA, HSD Tukey), except for samples S80 and S100, where the hypothesis h0 was not rejected.

4. Discussion

A relatively frequently used method for determining the degree of thixotropy is to determine the hysteresis area [29]. It is assumed that the value of the hysteresis area will be zero in the case of a Newtonian liquid (lime nectar honey, sample S0), or it will be close to zero. Conversely, the hysteresis area value, in the case of non-Newtonian fluid (heather honey, sample S100), will take on certain positive values (the samples show structural breakdown), as is clear from other professional work dealing with the rheological properties of heather honey e.g., [1]. It is further assumed that in the case of diluting heather honey by lime nectar honey, the hysteresis area values will gradually decrease as the proportion of lime nectar honey increases.
From the resulting values of the hysteresis area, it is clear that the previously mentioned hypothesis is not quite in accordance with the measured data. The point is, that the mean values of the hysteresis area in the case of sample S0 takes values other than zero, even became negative.
The values of a hysteresis area also often showed the very high degree of a variability at various samples (the highest variation coefficient was approx. 400%, see sample S20 in Table 2). It is, therefore, difficult to accept hysteresis area as an evaluation parameter. These findings need to be commented on in more detail.
The hysteresis loop test is still relatively popular due to its time-saving potential and simplicity. At the same time, however, this test has been criticised from several points of view. For example, the measurements are affected by the inertia effect (although this effect is already often removed using the software that is part of the measuring device). Furthermore, for example, the shear rate is dependent on both shear rate and on time for time-dependent materials. However, in the hysteresis loop test, we are not able to separate the effects of these two variables [35]. The size of the hysteresis area is also dependent on the conditions under which the test is performed, e.g., the shear history prior to the start of the experiment, the maximum shear rate and the acceleration rate [36]. However, it follows that the hysteresis area is a relative parameter dependent on factors that we cannot influence and, most importantly, we are not able to recognise. This is especially true of the sample history. We can never be sure how the material was handled before it was in the laboratory. The time required to restore the structure of the thixotropic material can be very different [37,38]. It should also be noted that honey is a viscous liquid in which solid particles (e.g., pollen grains) are dispersed. These particles settle naturally. The dynamics of these changes may be different for different honeys due to the different densities of the solutions. The dynamics of these processes are also continuously disrupted by taking a new sample when dosing honey into the rheometer. These factors can explain the high variability in the measured values of the hysteresis area. By increasing the number of repetitions, it is possible to reduce the degree of a variability. The number of repetitions may generally seem a weakness in these types of experiments. An above-standard number of repetitions was chosen in the work, i.e., four repetitions. In most works of this type, three repetitions are performed. However, even such a number of repetitions seems to be insufficient. However, the situation cannot be solved by increasing the number of repetitions, as this would be very time consuming and the test sets would not be completed in one day. As mentioned above, this type of sample changes over time, and, therefore, it is possible that over the next few days we would be working with a “different” sample, which would already have different properties than in the previous days, even if only slightly. It also does not seem practical to increase the number of repetitions when there are other (more stable) rheological parameters, which have a potential for the evaluation of the diluted honeys.
Another fact that needs to be clarified is the measured negative value of the hysteresis area in the case of sample S0 (pure lime nectar honey). This type of honey has not been reported as acting as a time dependent substance [8]. However, for example the phenomenon of a thixotropy can occur when honey contains crystals [39,40]. However, the honey samples were thermally thoroughly de-crystallised. The content of the crystals was examined under a light microscope at high magnification (600 times) and with a phase-contrast. Moreover, if there were any crystals in the original sample, they would be present in the all-tested samples and the impact of the crystals on the rheology would be similar in the all samples. In addition, the negative value of a hysteresis area shows rather to shear-thickening fluid with rheopectic behaviour. A closer view of the index n , however, we find that the curve expressing the dependence of the shear stress on the shear rate pointing up (shear rate from 0 s−1 to 100 s−1) is described by the power function with factor n u p = 0.9924   ±   0.007 (see Ostwald-de Waele model), in the case sample S0. However, the curve pointing down (shear rate from 100 s−1 to 0 s−1) is described by the power function with factor n d o w n = 0.9692   ±   0.014 . So first the curve is linear, then it arches. This is not a typical property of the shear-thickening fluid with rheopectic behaviour. Taking into account this fact, other measured parameters, and a relatively high degree of variability of the hysteresis loop in the sample S0, the sample S0 cannot be considered as shear-thickening fluid with rheopectic behaviour, but Newtonian fluid.
For the above reasons and also on the basis of information from scientific publications e.g., [36], it is necessary to approach the use of the hysteresis area as a relative parameter for evaluating the degree of the thixotropy of a honey in a somewhat reserved manner, and we are inclined rather not to use HA as a relative evaluation parameter. However, this contradicts the results of the work of [11] who, in turn, state that the area of a hysteresis loop is a suitable parameter for heather honey authentication. This assumption was probably made on the basis of the principal component analysis, where the correlations between the four rheological parameters were determined, with the strongest correlation being between hysteresis area and B (index of the Weltman model). The other parameters in the work [11] are not included in the principal component analysis. Thus, the question is: if the authors [11] had included the other parameters in the analysis, would strong correlations be found between them or not?
In the next step, the parameters of the rheological models were compared. Two rheological models were used: The Herschel–Bulkley model and the Ostwald–de Waele model. These two models are also used in other professional work on rheological properties, the Herschel–Bulkley model was used in [1,40,41] and the Ostwald–de Waele model was used in [11,41]. The results of these studies are summarised in Table 7. If we compare the measured values in sample S100 with the pure heather honey with the values derived from the literature sources, it is clear that the n -values from the literature are higher. The n -value is a relative parameter and the value of this parameter will be influenced by many other factors, as in the case of the hysteresis area. The design of the experiment is an important factor. Some of these are, for example, the maximum achieved value of the shear rate ( γ ˙ m a x ) , the temperature of the sample, etc. The same value ( γ ˙ m a x ) was set, as in our case, in the works [1,11].
Another parameter that can be compared is the yield stress ( τ 0 ). If the yield stress value ( τ 0 ), is measured in the sample, a non-Newtonian behaviour can be expected. However, [40] state that the yield stress ( τ 0 ) can be even in the case of honeys, which are considered Newtonian. The yield stress is represented as a finite stress required to achieve flow and, at lower temperatures (10, 20 °C), the yield stress occurs as an effect of the microparticles (crystals) in the honey [40]. From the measured values, it is clear that in the case of the gradual addition of heather honey to lime nectar honey, the values of τ 0 are first close to zero in the samples with a proportion of heather honey up to 40%. A higher proportion of heather honey led to a sharp increase in the parameter τ 0 . Negative yield stress values ( τ 0 ) are physically unrealistic [42]. However, the yield stress values are estimated based on a function describing the flow curve. For this reason, the yield stress should be interpreted as a model parameter rather than the actual yield stress value, which is a physical property of the fluid itself [1]. In the work [41] negative yield stress values ( τ 0 ) were even measured for the pure heather honey.
The parameters n , τ 0 showed the non-linear changes in dependence on the dilution as the essential attribute with a significant change between the 40% and 60% dilution. The same applies to other parameters, such as parameter B in the Weltman model, parameter ϕ , or parameter C in the case of the logarithmic dependence of the complex viscosity on the angular frequency.
The reliable parameters for the description of the rheological features of the diluted samples seems to be the n -value according to both models, Ostwald–de Waele and Herschel–Bulkley, and parameter ϕ , in which the most stable parameter with the lowest variability for the characterisation of the thixotropic rate of a honey was proposed.
The standard for heather honey is defined very widely, namely, as far as melissopalynological characteristics is concerned [8]. This situation makes it difficult exactly to decide on the botanical origin of a honey or its violation in spite there are also other chemical-physical markers for heather honey (e.g., specific rheological properties or high content of protein). For instance, a heather honey with a high content of pollen grains (e.g., over 10,000 PG/g and 70% share of the heather PG) can be diluted 1:1 or more while still remaining within the range for the unifloral characteristic for heather honey proposed by [8]. The rheological parameters are relative values and their share in the detection of honey botanic purity has still the characteristic of a supporting, not an arbitrary, method. In the future research, it is appropriate to pay attention to this area, where there is a step change in the measured parameters. It is necessary to examine whether this area changes depending on the choice of the diluent (e.g., different types of honey, special sugar solutions, etc.), and possibly determine how much it changes. If this area were found to be within the same range for different diluents and heather honeys from different localities, this finding could be used to detect honey defects in relation to the purity of the botanical origin. For example, if considered that a step change will be demonstrable in the specified dilution area and testing a random sample of heather honey would show that this area will be shifted towards a lower degree of dilution (i.e., a higher proportion of the random sample than the diluent that may be, e.g., lime nectar honey), this may be a signal that the heather honey has been devalued.

5. Conclusions

  • The area value of the hysteresis loop is sensitive to both external and internal factors. This is an unstable parameter and we do not recommend using it as a comparison parameter.
  • Due to the low variability of the measured values, the relative comparison n-value parameter and parameter ϕ appear to be suitable.
  • The dependence of the measured parameters n, B, ϕ and C on the degree of the dilution is non-linear and a distinct step change occurs between samples S40 and S60, i.e., the samples that contained 40% (w/w) heather honey and 60% (w/w), respectively.
The fact that the dependence of the measured parameters on the degree of the dilution is non-linear with a step change in a certain phase of the dilution is a new finding in the field of the rheological behaviour of honeys. This knowledge could be used, for example, to identify adulterated heather honeys. For this reason, more attention should be paid to the dilution interval between 40% (w/w) and 60% (w/w) of heather honey.
It is also important to consider the number of repetitions in this regard. Usually, three replicates are performed in similar experiments. Four replicates were performed in our experiment. However, the mean values calculated from the measurement results from a small number of replicates are prone to higher uncertainty. We can refine the decision-making process by increasing the number of repetitions, but this can lead to difficulties in organising and securing the experiment. The experiment may become unfeasible. Already, the total number of four repetitions was relatively time consuming and, in some cases, this number of repetitions proved to be insufficient (see, for example, the determination of the area of the hysteresis loop). A suitable strategy may be to perform pre-tests with a lower number of repetitions, but on a larger range. Then focus on a narrow area that not behaves as expected and perform an increased number of repetitions in the experiment. In our case, this may be the area defined by the dilution interval of 40% (w/w) and 60% (w/w) of heather honey.

Author Contributions

Conceptualization: A.P. and P.T.; methodology: A.P., P.T. and V.K.; validation: A.P. and P.T.; formal analysis: A.P., P.T. and V.K.; investigation: J.K., Š.N., S.O., D.T.; data curation: J.K.; resources: A.P.; writing—original draft preparation: A.P., P.T., Š.N., S.O., D.T. and V.K.; project administration: P.T.; funding acquisition: V.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Internal Grant Agency of the Faculty of AgriSciences no. AF-IGA-2020-TP006 “Modelling of rheological properties of liquid and semi-liquid food raw materials and foods showing non-Newtonian behaviour”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Appendix A

Table A1. Post-hoc analysis by Tukey’s test of the hysteresis area values.
Table A1. Post-hoc analysis by Tukey’s test of the hysteresis area values.
SampleS0S10S20S40S60S80S100
S0X------
S100.9999X-----
S200.99500.9999X----
S400.36520.55960.7345X---
S600.00020.00020.00020.0002X--
S800.00020.00020.00020.00020.0482X-
S1000.00020.00020.00020.00020.00020.0427X
Table A2. Post-hoc analysis by Tukey’s test of the φ values.
Table A2. Post-hoc analysis by Tukey’s test of the φ values.
SampleS0S10S20S40S60S80S100
S0X------
S100.27992X-----
S200.096630.99876X----
S40<0.0010.044490.14767X---
S60<0.001<0.001<0.001<0.001X--
S80<0.001<0.001<0.001<0.0010.00133X-
S100<0.001<0.001<0.001<0.001<0.0010.5623X
Table A3. Post-hoc analysis by Tukey’s test of the B values.
Table A3. Post-hoc analysis by Tukey’s test of the B values.
SampleS0S10S20S40S60S80S100
S0X------
S100.999368X-----
S200.9997711.000000X----
S400.9273890.9944050.990261X---
S600.2522950.4721610.4315020.841039X--
S800.0228380.0574870.0496940.1949410.873087X-
S1000.0008790.0021720.0018690.0089770.1374720.739904X

References

  1. Witczak, M.; Juszczak, L.; Gałkowska, D. Non-Newtonian behaviour of heather honey. J. Food Eng. 2011, 104, 532–537. [Google Scholar] [CrossRef]
  2. Dezmirean, D.S.; Mărghitaș, L.A.; Fiț, N.; Chirilă, F.; Gherman, B.; Mărgăoan, R.; Bobiș, O. Antibacterial Effect of Heather Honey (Calluna vulgaris) against Different Microorganisms of Clinical Importance. Bulletin of University of Agricultural Sciences and Veterinary Medicine Cluj-Napoca. Anim. Sci. Biotechnol. 2015, 72, 72–77. [Google Scholar] [CrossRef]
  3. Lehébel-Péron, A.; Sidawy, P.; Dounias, E.; Schatz, B. Attuning local and scientific knowledge in the context of global change: The case of heather honey production in southern France. J. Rural Stud. 2016, 44, 132–142. [Google Scholar] [CrossRef]
  4. Rodríguez-Flores, M.S.; Escuredo, O.; Seijo-Rodríguez, A.; Seijo, M.C. Characterization of the honey produced in heather communities (NW Spain). J. Apic. Res. 2019, 58, 84–91. [Google Scholar] [CrossRef]
  5. Frydman, G.H.; Olaleye, D.; Annamalai, D.; Layne, K.; Yang, I.; Kaafarani, H.M.A.; Fox, J.G. Manuka honey microneedles for enhanced wound healing and the prevention and/or treatment of Methicillin-resistant Staphylococcus aureus (MRSA) surgical site infection. Sci. Rep. 2020, 10, 13229. [Google Scholar] [CrossRef] [PubMed]
  6. Le Conte, Y.; Navajas, M. Climate change: Impact on honey bee populations and diseases. Rev. Sci. Tech. Off. Int. Epizoot. 2008, 27, 499–510. [Google Scholar]
  7. Zábrodská, B.; Vorlová, L. Adulteration of honey and available methods for detection—A review. Acta Vet. Brno 2014, 83, S85–S102. [Google Scholar] [CrossRef] [Green Version]
  8. Persano Oddo, L.; Piro, R. Main European unifloral honeys: Descriptive sheets. Apidologie 2004, 35, S38–S81. [Google Scholar] [CrossRef]
  9. Bogdanov, S.; Martin, P.; Lüllman, C. Harmonised methods of the European Honey Commission. Apidologie 1997, 28, 1–59. [Google Scholar]
  10. Louveaux, J. Essai de caractérisation des miels de Callune (Calluna vulgaris Salisb.). Annales de L’Abeille 1966, 9, 351–358. [Google Scholar]
  11. Osés, S.M.; Ruiz, M.O.; Pascual-Maté, A.; Bocos, A.; Fernández-Muiño, M.A.; Sancho, M.T. Ling heather honey authentication by thixotropic parameters. Food Bioprocess Technol. 2017, 10, 973–979. [Google Scholar] [CrossRef]
  12. Munro, J.A. The viscosity and thixotropy of honey. J. Entomol. 1943, 36, 769–777. [Google Scholar] [CrossRef]
  13. Waś, E.; Rybak-Chmielewska, H.; Szczęsna, T.; Kachaniuk, K.; Teper, D. Characteristic of polish unifloral honey. III. Heather honey (Calluna vulgaris L.). J. Apic. Sci. 2011, 55, 129–136. [Google Scholar]
  14. Aliaño-González, M.J.; Ferreiro-González, M.; Espada-Bellido, E.; Palma, M.; Barbero, G.F. A screening method based on Visible-NIR spectroscopy for the identification and quantification of different adulterants in high-quality honey. Talanta 2019, 203, 235–241. [Google Scholar] [CrossRef]
  15. Song, Y.Q.; Milne, R.I.; Zhou, H.X.; Ma, X.L.; Fang, J.Y.; Zha, H.G. Floral nectar chitinase is a potential marker for monofloral honey botanical origin authentication: A case study from loquat (Eriobotrya japonica Lindl.). Food Chem. 2019, 282, 76–83. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Q.; Zhao, H.; Zhu, M.; Zhang, J.; Cheng, N.; Cao, W. Method for identifying acacia honey adulterated by resin absorption: HPLC-ECD coupled with chemometrics. LWT Food Sci. Technol. 2020, 118, 108863. [Google Scholar] [CrossRef]
  17. Geana, E.I.; Ciucure, C.T. Establishing authenticity of honey via comprehensive Romanian honey analysis. Food Chem. 2020, 306, 125595. [Google Scholar] [CrossRef]
  18. Yilmaz, M.T.; Tatlisu, N.B.; Toker, O.S.; Karaman, S.; Dertli, E.; Sagdic, O.; Arici, M. Steady, dynamic and creep rheological analysis as a novel approach to detect honey adulteration by fructose and saccharose syrups: Correlations with HPLC-RID results. Food Res. Int. 2014, 64, 634–646. [Google Scholar] [CrossRef] [PubMed]
  19. Kamboj, U.; Mishra, S. Prediction of Adulteration in Honey Using Rheological Parameters. Int. J. Food Prop. 2015, 18, 2056–2063. [Google Scholar] [CrossRef]
  20. Oroian, M.; Ropciuc, S.; Paduret, S.; Todosi, E. Rheological analysis of honeydew honey adulterated with glucose, fructose, inverted sugar, hydrolysed inulin syrup and malt wort. LWT Food Sci. Technol. 2018, 95, 1–8. [Google Scholar] [CrossRef]
  21. Anidiobu, V.O.; Nwalor, J.U.; Babalola, F.U. Comparative Study of Rheological Characterization and Classification of Honeys with Their Physicochemical Properties. Int. J. Food Eng. 2019, 5, 268–275. [Google Scholar] [CrossRef]
  22. Nayik, G.A.; Dar, B.N.; Nanda, V. Physico-chemical, rheological and sugar profile of different unifloral honeys from Kashmir valley of India. Arab. J. Chem. 2019, 12, 3151–3162. [Google Scholar] [CrossRef] [Green Version]
  23. Kurt, A.; Palabiyik, I.; Gunes, R.; Konar, K.; Toker, O.S. Determining Honey Adulteration by Seeding Method: An Initial Study with Sunflower Honey. Food Anal. Methods 2020, 13, 952–961. [Google Scholar] [CrossRef]
  24. Szczęsna, T.; Rybak-Chmielewska, H. The temperature correction factor for electrical conductivity of honey. J. Apic. Sci. 2004, 48, 97–102. [Google Scholar]
  25. Ohe, W.; Persano Oddo, L.; Piana, M.L.; Morlot, M.; Martin, P. Harmonized methods of melissopalynology. Apidologie 2004, 35, S18–S25. [Google Scholar] [CrossRef]
  26. Maurizio, A.; Louveaux, J. Union des Groupements Apicoles Français, 1st ed.; Pollens de Plantes Melliféres d’Europe: Paris, France, 1965; p. 148. [Google Scholar]
  27. Sawyer, R. Honey Identification, 1st ed.; Cardiff Academic Press: Cardiff, UK, 1988; p. 115. [Google Scholar]
  28. Ohe, K.; Ohe, W. LAVES—Institut für Bienenkunde, 3rd ed.; Celle’s Melissopalynological Collection: Celle, Germany, 2007; p. 236. [Google Scholar]
  29. Dapčević, T.; Dokić, P.; Hadnađev, M.; Krstonošić, V. An approach in numerical evaluation of thixotropy. Food Process. Qual. Saf. 2008, 35, 33–39. [Google Scholar]
  30. Tovar, C.; Rodríguez-Flores, M.S.; Escuredo, O.; Del Carmen Seijo, M. Rheology of Honey-Advances in Rheology Research, 1st ed.; Nova Science Publisher: New York, NY, USA, 2017; pp. 175–191. [Google Scholar]
  31. Karasu, S.; Toker, O.S.; Yilmaz, M.T.; Karaman, S.; Dertli, E. Thermal loop test to determine structural changes and thermal stability of creamed honey: Rheological characterization. J. Food Eng. 2015, 150, 90–98. [Google Scholar] [CrossRef]
  32. Kumbár, V.; Strnková, J.; Nedomová, Š.; Buchar, J. Fluid dynamics of liquid egg products. J. Biol. Phys. 2015, 41, 303–311. [Google Scholar] [CrossRef] [Green Version]
  33. Kumbár, V.; Nedomová, Š.; Trnka, J.; Buchar, J.; Pytel, R. Effect of storage duration on the rheological properties of goose liquid egg products and eggshell membranes. Poult. Sci. 2016, 95, 1693–1701. [Google Scholar] [CrossRef]
  34. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria; Available online: http://www.R-project.org/ (accessed on 23 March 2020).
  35. Barnes, H.A. Thixotropy—A review. J. Non Newton. Fluid Mech. 1997, 70, 1–33. [Google Scholar] [CrossRef]
  36. Mewis, J.; Wagner, N.J. Thixotropy. Adv. Colloid Interface Sci. 2009, 147–148, 214–227. [Google Scholar] [CrossRef]
  37. Adamczyk, G.; Krystyjan, M.; Dobosz, A.; Sikora, M. Thixotropic properties of starch. Food Sci. Technol. Qual. 2013, 20, 16–31. [Google Scholar] [CrossRef]
  38. Razavi, S.M.A.; Behrouzian, F.; Alghooneh, A. Temperature dependency of the interaction between xanthan gum and sage seed gum: An interpretation of dynamic rheology and thixotropy based on creep test. J. Texture Stud. 2017, 48, 470–484. [Google Scholar] [CrossRef] [PubMed]
  39. Stelmakienė, A.; Ramanauskiene, K.; Briedis, V.; Leskauskaite, D. Examination of rheological and physicochemical characteristics in Lithuanian honey. Afr. J. Biotechnol. 2012, 11, 12406–12414. [Google Scholar] [CrossRef]
  40. Samanalieva, J.; Senge, B. Analytical and rheological investigations into selected unifloral German honey. Eur. Food Res. Technol. 2009, 229, 107–113. [Google Scholar] [CrossRef]
  41. Afonso, M.J.; Magalhães, M.; Fernandes, L.; Castro, M.; Ramalhosa, E.C.D. Temperature Effect on Rheological Behavior of Portuguese Honeys. Pol. J. Food Nutr. Sci. 2018, 68, 217–222. [Google Scholar] [CrossRef] [Green Version]
  42. Nedomova, S.; Kumbar, V.; Pytel, R.; Buchar, J. Mechanical properties of sugar beet root during storage. Int. Agrophys. 2017, 31, 507–513. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The dependence of the shear stress on the shear rate (sample: S80 and S100).
Figure 1. The dependence of the shear stress on the shear rate (sample: S80 and S100).
Materials 14 02472 g001
Figure 2. The dependence of the shear stress on the shear rate (sample: S0, S10, S40, S60).
Figure 2. The dependence of the shear stress on the shear rate (sample: S0, S10, S40, S60).
Materials 14 02472 g002
Figure 3. Hysteresis area average values for samples S0–S100.
Figure 3. Hysteresis area average values for samples S0–S100.
Materials 14 02472 g003
Figure 4. Flow index n average values for samples S0–S100.
Figure 4. Flow index n average values for samples S0–S100.
Materials 14 02472 g004
Figure 5. The dependence of the apparent viscosity on the time (samples: S0, S10, S40, S60).
Figure 5. The dependence of the apparent viscosity on the time (samples: S0, S10, S40, S60).
Materials 14 02472 g005
Figure 6. Parameter ϕ average values for samples S0–S100.
Figure 6. Parameter ϕ average values for samples S0–S100.
Materials 14 02472 g006
Figure 7. Parameter B values for samples S0–S100.
Figure 7. Parameter B values for samples S0–S100.
Materials 14 02472 g007
Figure 8. The dependence of the loss modulus ( G ) and storage modulus ( G ) on the angular frequency in the selected samples.
Figure 8. The dependence of the loss modulus ( G ) and storage modulus ( G ) on the angular frequency in the selected samples.
Materials 14 02472 g008
Figure 9. The dependence of the complex viscosity ( η * ) on the angular frequency in the selected samples.
Figure 9. The dependence of the complex viscosity ( η * ) on the angular frequency in the selected samples.
Materials 14 02472 g009
Figure 10. Parameter C values for samples S0–S100.
Figure 10. Parameter C values for samples S0–S100.
Materials 14 02472 g010
Table 1. Description of the honey samples.
Table 1. Description of the honey samples.
SampleBotanical OriginGeographical OriginPollen Grains in 1 g of HoneyPercentage of Important
Pollen Grains
Important Pollen Grains in 1 g of HoneyHoneydew ElementsWater Content
(%)
Electrolytic Conductivity
(mS·cm−1)
PA/386heather nectar honeyNorway8886 ± 46728.1 ± 1.21786 ± 21sporadic clusters, no bodies16.9 ± 0.286.6 ± 4.2
PA/416lime nectar honeyCzech Republic2252 ± 1399.6 ± 0.5216 ± 1sporadic clusters, no bodies17.8 ± 0.150.4 ± 1.9
Table 2. Exploratory data analysis—hysteresis area.
Table 2. Exploratory data analysis—hysteresis area.
SampleMean
(Pa·s−1)
Median
(Pa·s−1)
SD
(Pa·s−1)
VC
(%)
IQR
(Pa·s−1)
S-W
(-)
S0−493.87−532.96240.619−48.7208320.6310.72
S10−217.28−245.9769.333−31.909582.4030.05
S2020.390.6781.724400.7579112.6160.54
S401254.581476.30629.88950.2071769.8390.14
S607742.537386.651342.52917.33972020.3500.39
S8010,409.8310,577.502469.51023.72293373.3500.91
S10013,124.5012,859.00974.2237.42291473.0000.35
SD—standard deviation, VC—variation coefficient, IQR—interquartile range, S-W—Shapiro-Wilk test (the measurement of each sample was performed in four replicates).
Table 3. Resulting parameters of the rheological models.
Table 3. Resulting parameters of the rheological models.
Ostwald-De Waele ModelHerschel-Bulkley Model
SampleK
(Pa·sn)
n
(-)
R2
(-)
τ 0
( P a )
K
(Pa·sn)
n
(-)
R2
S07.91 ± 0.30.9924 ± 0.0071.0000−1.15 ± 2.28.0 ± 0.40.9891 ± 0.0111.0000
S108.79 ± 0.30.9864 ± 0.0031.00000.29 ± 1.78.8 ± 0.50.9872 ± 0.0081.0000
S2010.35 ± 0.10.9768 ± 0.0031.0000−0.02 ± 1.210.4 ± 0.20.9768 ± 0.0031.0000
S4017.39 ± 1.10.9265 ± 0.0031.0000−2.86 ± 2.817.8 ± 1.40.9220 ± 0.0071.0000
S6043.80 ± 2.10.7970 ± 0.0060.999817.28 ± 7.840.2 ± 2.30.8141 ± 0.0140.9998
S8067.22 ± 6.10.7634 ± 0.0170.999536.45 ± 6.758.8 ± 5.50.7898 ± 0.0190.9997
S10074.50 ± 6.20.6745 ± 0.0190.998644.94 ± 22.961.0 ± 6.40.7134 ± 0.0160.9990
±standard deviation (the measurement of each sample was performed in four replicates).
Table 4. Exploratory data analysis of the ϕ values.
Table 4. Exploratory data analysis of the ϕ values.
SampleMean
(-)
Median
(-)
SD
(-)
VC
(%)
IQR
(-)
S-W
(-)
S01.05341.02600.06100.03385.79360.01092
S101.01081.01100.00430.00620.42180.48
S201.00051.00490.01310.01431.30600.1875
S400.95140.95220.01430.02331.49840.1185
S600.82510.82580.02090.03022.53220.5364
S800.73900.73900.01840.02972.48380.1605
S1000.70540.70740.00750.00751.06500.4843
SD—standard deviation, VC—variation coefficient, IQR—interquartile range, S-W—Shapiro-Wilk test (the measurement of each sample was performed in four replicates).
Table 5. Parameters of the Weltman model.
Table 5. Parameters of the Weltman model.
SampleA (Pa)B (-)R2
S0368.5 ± 14.2−4.5 ± 2.80.83 ± 0.13
S10441.3 ± 9.4−1.7 ± 1.30.68 ± 0.17
S20515.9 ± 27.0−2.2 ± 1.60.55 ± 0.34
S40673.7 ± 22.72.3 ± 2.20.76 ± 0.19
S60864.5 ± 111.410.6 ± 7.00.75 ± 0.17
S801284.8 ± 72.218.4 ± 12.40.73 ± 0.21
S1001064.2 ± 75.628.0 ± 18.50.84 ± 0.19
±standard deviation (the measurement of each sample was performed in four replicates).
Table 6. The parameters of the model describing the logarithmic dependence of the complex viscosity ( η * ) on the angular frequency ( ω ) .
Table 6. The parameters of the model describing the logarithmic dependence of the complex viscosity ( η * ) on the angular frequency ( ω ) .
SampleC (-)R2
S0−1.93 ± 1.340.81 ± 0.2
S10−0.91 ± 0.850.66 ± 0.5
S20−0.61 ± 1.050.31 ± 0.4
S40−2.93 ± 3.220.92 ± 0.1
S60−5.05 ± 1.430.95 ± 0.01
S80−18.6 ± 3.900.95 ± 0.002
S100−15.7 ± 2.250.95 ± 0.01
±standard deviation (the measurement of each sample was performed in four replicates).
Table 7. Heather honeys’ parameters by various authors.
Table 7. Heather honeys’ parameters by various authors.
Herschel-BulkleyOstwald-de Waele
Sourcen
(-)
K
(Pa·sn)
τ 0
( P a )
n
(-)
K (Pa·sn)t
(°C)
ϕ
(%)
HA
(Pa·s−1)
γ ˙ m a x
  ( s 1 )
(A)0.70 ± 0.0150.7 ± 1.250.2 ± 1.3--2018.015,000100
0.77 ± 0.0129.1 ± 0.350.2 ± 1.9--2018.27000
0.88 ± 0.0110.8 ± 0.63.8 ± 0.2--20202000
(B)0.90113.390.15--2018.7-50
(C)0.988112−0.640.9964.713024.0-5
(D)---0.88 ± 0.123.67 ± 8.02517.56994 ± 1945100
(A)0.70 ± 0.0150.7 ± 1.250.2 ± 1.3--2018.015,000100
t —temperature, φ —humidity, HA—hysteresis area, γ ˙ m a x —maximal value of shear rate. (A)—Ref. [1], (B)—Ref. [40], (C)—Ref. [41], (D)—Ref. [11].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Přidal, A.; Trávníček, P.; Kudělka, J.; Nedomová, Š.; Ondrušíková, S.; Trost, D.; Kumbár, V. A Rheological Analysis of Biomaterial Behaviour as a Tool to Detect the Dilution of Heather Honey. Materials 2021, 14, 2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14102472

AMA Style

Přidal A, Trávníček P, Kudělka J, Nedomová Š, Ondrušíková S, Trost D, Kumbár V. A Rheological Analysis of Biomaterial Behaviour as a Tool to Detect the Dilution of Heather Honey. Materials. 2021; 14(10):2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14102472

Chicago/Turabian Style

Přidal, Antonín, Petr Trávníček, Jan Kudělka, Šárka Nedomová, Sylvie Ondrušíková, Daniel Trost, and Vojtěch Kumbár. 2021. "A Rheological Analysis of Biomaterial Behaviour as a Tool to Detect the Dilution of Heather Honey" Materials 14, no. 10: 2472. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14102472

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop