Hydric Properties Evolution of Spirulina platensis during Drying: Experimental Analysis and Modeling

Eloi Salmwendé Tiendrebeogo^{1,2,3}^{*},
Guy Christian Tubreoumya^{1},
A. O. Dissa^{1},
A. Compaoré^{1},
Jean Koulidiati^{1},
F. Cherblanc^{2},
J.-C. Bénet^{2},
I. Youm^{3}

Show more

1. Introduction

In many countries of sub-Saharan Africa, food imports absorbe for a large proportion of national income. For example, in Burkina Faso, per capita food production is now lower than it was 40 years ago. Post-harvest losses related to the degradation of agro-food products contribute to the decline in agricultural productivity. Spirulina platensis, a nutritional and therapeutic seaweed, is used today for complementary food requirements [1] [2] [3] . Cultivated all over the world under several techniques, it is highly porous and deteriorates rapidly in the presence of water. Methods have been proposed in previous studies to produce healthy and nutritious foods, reduce losses associated with the development of severe microorganisms and improve their quality [4] [5] [6] [7] .

Drying can make a significant contribution to sustainable agricultural development and improved food security by reducing post-harvest losses and improving the shape of dried produce and their nutritional quality. This process, which aims to strengthen the bonds of water molecules by reducing the latter, is one of the most used techniques to protect the product from degrading agents and to extend its conservation. In addition to stresses and cracks, the drying process often induces water density gradients in the solid matrix. The distribution of water in the solid matrix during drying differs from one product to another. Water is often trapped inside the structure by a crusting phenomenon that accelerates the product degradation. The hydrics properties evolution during drying differs from one product to another and has been the subject of various studies because of its crucial importance in the modeling of the drying process [8] - [13] . Some food products such as the Spirulina platensis continue to be victim mechanicals and biologicals stress.

The diffusion of water in Spirulina platensis had already been approached experimentally by Desmorieux and Dissa [14] [15] in thin layer. This method does not allow to take into account the contribution of the internal phenomena and neglects the influence of the product shrinkage. The complexity of material shrinkage has led some authors [16] - [24] to develop numerical methods in support of experimental data to predict the withdrawal of food products such as potato, okra, banana, Spirulina, cereals, onion etc. The results they have obtained show that the influence of the solid contraction on the transport of water is negligible. These earlier studies are generally carried out under very simplifying assumptions that lead to limited results. The variation of the hydric characteristic parameters in the solid matrix and in time during the drying of the Spirulina platensis has not known an advanced understanding. In this article, the overall objective is to evaluate the evolution of the hydric characteristic parameters of the Spirulina platensis.

The present study is divided into three parts. First, we will determine experimentally the water content profiles of the Spirulina platensis. Then we will evaluate the evolution of the mass flow and concentration gradient and establish a relationship between them. Finally, we will estimate the transport coefficient.

2. Method and Materials

2.1. Origin, Materials Pretreatment

The materials used in this study are Spirulina platensis samples from two farm sites located in two different areas:

“Loumbila” farm situated at 15 km in the north of Ouagadougou, Burkina Faso. This material sample will be named sample-B,

“La Fon del Cardaire” farm situated at Gignac in the south of France. This material will be named sample-F.

Their mode of production is not different from one farm to another, from one country to another. These are the methods that change, adapted to the size of the farm, the workforce and the materials available on site. Indeed, the “sample-B” is grown in open ponds and the “sample-F” in glass ponds. However, if these two products were grown in localities with different climates, their seeding is done by strains of the same type “platensis” in basins whose salinity and alkalinity are provided by sodium bicarbonate (NaHCO_{3}). Pretreatments were carried out on each biomass harvested to remove a good part of the culture water until reaching a given water content. The methods pretreatment used is different from one farm to another depending on the available materials. Thus, the water content of the sample-B is reduced to about 2.73 kg_{w}/kg_{db} by manual pressing while that of the sample-F is reduced to about 3.12 kg_{w}/kg_{db} by a spin under the action of vacuum. These water contents obtained after pretreatment are considered as the water contents initial of the samples. The initial characteristic parameters recorded in Table 1 have been determined using standard methods. The volume of a substance corresponds to the space occupied by this substance. For measurements, a container of cylindrical shape with known dimensions is previously weighed empty so that its mass can be measured and subtracted from that of the container loaded mass m. It is then in an oven maintained at 70˚C for 48 hours to determine the anhydrous mass m_{s}.

2.2. Sampling

First, the filtering of biomass through one or two devices superimposed. In general the first consists of a thin mesh of the order of 300 microns which retains any impurities (lumps, insects, leaves). The second is a filter in nylon is of the

Table 1. Initial characteristics of the materials.

order of 50 μm. After filtering, fresh Spirulina platensis is processed into test-ready samples. A cylindrical shape commonly used in farms for the drying is proposed (Figure 1). To avoid variability and ensure identical initial conditions between samples, the following design process has been adopted: the fresh biomass is moulded in a cylindrical device of 20 mm diameter and 40 mm length composed of a tube and a piston made in Teflon; the cylindrical surface of samples are wrapped with plastic film to prevent adhesion between samples. Both end faces are free to ensure a one-dimensional water transport along sample axis; a sample is kept if its initial mass is m_{0} = 0.013 ± 0.0006 kg.

It is assumed that these conditions are met and the samples are uniform in the initial state.

To understand experimentally the distribution of water in the product during drying, we chose to follow the variation of the water content in the product by cutting samples. The water content in a sample slice is determined by considering the mass of this slice at a given instant and that of the solid phase according to the following Equation (1):

$w=\frac{m\left(x,t\right)-{m}_{s}}{{m}_{s}}$ (1)

2.3. Experimental Device and Procedure

The experimental setup is constituted of a desiccator where the relative humidity is regulated at 6% by a potassium hydroxide KOH solution (Figure 2). The cylindrical samples (Figure 2(a)) are put above the solution (Figure 2(b)) while the desiccator is placed in a thermo-regulated chamber (Memmert UFP 600) where the temperature is maintained at 50˚C. At specific time steps of the drying process, one of the samples is removed from the desiccator and immediately ox introduced in a cutting device to obtain slices of 2 mm thickness. The water content of each slice is determined by differential weighing after 48 hours drying at 70˚C for complete dehydration. In order for the desired water content profiles to be characteristic of a given single sample, all the samples should be in the same state at the process beginning. Thus, they are prepared under the same

Figure 1. Study sample cylinder.

Figure 2. Experimental Principe setup for the characterization of water transport coefficient in Spirulina samples.

conditions as indicated above. The material being symmetrical, the water contents averages of symmetrical points with respect to the median plane are considered. This approach provides to access to the moisture content profiles $w\left(x,t\right)$ corresponding to the drying times (Figure 2(d)). Since we aim to characterize water transport phenomena in a single sample from measurements obtained on various samples, it is essential to have the same initial conditions in all samples. Although special attention is paid to sample preparation, this methodology may have some differences. For this, the experimental water content profiles are approximated by simple polynomial equations (Equation (2)) introducing three parameters a, b and c at each time step.

$w\left(x,t\right)=a{x}^{4}+b{x}^{2}+c$ (2)

This expression ensures to have a symmetrical evolution with a no flux condition at the symmetry plane. The Experiments are carried out and lasted at most 120 hours (5 days). As it for a drying period long the material begins to stiffen and the cutting deviate delicate. For low temperature (less than 25˚C), the drying is too slow for reach values with low water content and the range of water content investigated is too narrow. It appears that for longer time period, the material starts to denature and the development of micro-organisms are observed. This could change the process of water transport and should be avoided.

2.4. Hypothesis and Theoretical Modeling

Simplifying assumptions are considered to develop the descriptive equations of water transport: the material is assumed to remain diphasic along the experiment, the initial water content in the product is homogeneous, both solid and liquid phases are incompressible, water transport is one-dimensional along the cylindrical sample axis x.

The real mass densities are determined over intrinsic volume of water and solid phases:

${\rho}_{w}^{*}=\frac{{m}_{w}}{{V}_{w}}$ ; ${\rho}_{s}^{*}=\frac{{m}_{s}}{{V}_{s}}$ (3)

The material being biphasic, the apparent densities are expressed by [25] :

${\rho}_{w}={\rho}_{w}^{*}\frac{w}{\alpha +w}$ (4)

${\rho}_{s}={\rho}_{s}^{*}\frac{1}{\alpha +w}$ (5)

$\alpha ={\rho}_{w}^{*}/{\rho}_{s}^{*}$ is the ratio of real densities. Considering ${u}_{w}$ and ${u}_{s}$ as the water and solid phase velocities along x, mass balance equations are written [26] [27] :

$\frac{\partial {\rho}_{w}}{\partial t}+\frac{\partial {\rho}_{w}{u}_{w}}{\partial x}=0$ (6)

$\frac{\partial {\rho}_{s}}{\partial t}+\frac{\partial {\rho}_{s}{u}_{s}}{\partial x}=0$ (7)

The mass flow of water with respect to the solid phase is written as follows:

${F}_{w}={\rho}_{w}\left({u}_{w}-{u}_{s}\right)={\rho}_{w}{u}_{w}-w{\rho}_{s}{u}_{s}$ (8)

For reasons of symmetry, ${u}_{w}\left(0,t\right)={u}_{s}\left(0,t\right)=0$ and from balance Equations (6)-(7), we obtain by integration the Equations (9) and (10):

${\rho}_{w}{u}_{w}\left(x,t\right)-{\rho}_{w}{u}_{w}\left(0,t\right)={\rho}_{w}{u}_{w}\left(x,t\right)=-{\displaystyle \underset{0}{\overset{x}{\int}}\frac{\partial {\rho}_{w}}{\partial t}\text{d}x}$ (9)

${\rho}_{s}{u}_{s}\left(x,t\right)-{\rho}_{s}{u}_{s}\left(0,t\right)={\rho}_{s}{u}_{s}\left(x,t\right)=-{\displaystyle \underset{0}{\overset{x}{\int}}\frac{\partial {\rho}_{s}}{\partial t}\text{d}x}$ (10)

By replacing the expressions 9 and 10 in Equation (8) and considering the Equations (8) and (9), the relation of the water flow becomes:

${F}_{w}={\rho}_{w}^{*}\left(w{\displaystyle \underset{0}{\overset{x}{\int}}\frac{\partial}{\partial t}}\frac{1}{\alpha +w}\text{d}x-{\displaystyle \underset{0}{\overset{x}{\int}}\frac{\partial}{\partial t}}\frac{w}{\alpha +w}\text{d}x\right)$ (11)

Deriving the Equation (4), the density gradient of the water is written:

${G}_{w}={\rho}_{w}^{*}\frac{\partial}{\partial x}\left(\frac{w}{\alpha +w}\right)$ (12)

Fick’s law (Equation (13)) is a relationship that links the water flow to the mass concentration gradient through a phenomenological coefficient [28] .

${F}_{w}=-{D}_{w}\frac{\partial {\rho}_{w}}{\partial x}$ (13)

In the case of geometry and boundary conditions used in experiments, radial water transport can be neglected so that water content is considered uniform in a cross-section at abscissa x. From relations 11 and 13, the water transport coefficient is accessed as a function of the water content.

${D}_{w}=-\frac{{F}_{w}}{{G}_{w}}$ (14)

The experimental determination of the water transport coefficient from the exploitation of the water content profiles is carried out as follows. For a particular time step ${t}_{j}$ and a given water content ${w}_{k}$ , the fitting of the associated water content profile (Equation (2)) can be inverted to give the corresponding abscissa ${x}_{i}$ . Knowing this abscissa, the values of the water flow (Equation (11)) and those of the water density gradient (Equation (12)) can be calculated. In fact, in Figure 2(d), the horizontal line associated with the given water content ${w}_{k}$ can cut several water content profiles providing a set of values ( ${x}_{i}$ , ${t}_{i}$ ) leading to a set of couples $\left\{{F}_{w}\left({x}_{i},{t}_{j}\right),{G}_{w}\left({x}_{i},{t}_{j}\right)\right\}$ . This set of couples is associated with a single value of the moisture content ${w}_{k}$ and must correspond to a single value of the coefficient of transport ${G}_{w}$ . In this way, by scrolling through a range of water content values between $0.7<{w}_{k}<2.5$ , this approach makes it possible to determine the evolution of the transport coefficient ${D}_{w}\left(w\right)$ as a function of the water content. The precision on the determination of the experimental results is estimated by statistical parameters: the coefficient ${R}^{2}$ (Equation (15)): and the root mean square error (RMSE) (Equation (16)): the RMSE squared error and the coefficient (Equation (16)).

${R}^{2}=1-\frac{{\displaystyle {\sum}_{z=1}^{N}{\left({Y}_{exp,z}-{Y}_{pre,z}\right)}^{2}}}{{\displaystyle {\sum}_{z=1}^{N}{\left({\stackrel{\xaf}{Y}}_{exp,z}-{Y}_{pre,z}\right)}^{2}}}$ (15)

$RMSE=\sqrt{\frac{{\displaystyle {\sum}_{z=1}^{N}{\left({Y}_{exp,z}-{Y}_{pre,z}\right)}^{2}}}{{\displaystyle {\sum}_{z=1}^{N}{\left({\stackrel{\xaf}{Y}}_{exp,z}-{Y}_{pre,z}\right)}^{2}}}}$ (16)

3. Results and Discussion

3.1. Water Content Profiles

The experimental approach presented above made it possible to obtain a representative evolution of the water distribution inside of the two Spirulina platensis samples types studied (Figure 3 and Figure 4). These profiles are accurately estimated by a polynomial expression (Equation (2)) whose parametric coefficients at each time t of measurement are recorded in Table 2 and Table 3. Data for both samples were recorded differently. It appears that the recording time has no influence on the evolution of the parameters. This function constitutes a base of the calculation program for the coefficient of transport. The water content profiles (axial and temporal variation of the water content in the product) show that the drying process is more important at the ends than inside of the samples. A heterogeneous water distribution is observed along the cylindrical axis. Ce qui montre que la perte d’eau par évaporation à la surface d’échange n’est pas complètement compensée par le transport interne. Considerable shrinkage of the solid matrix is observed and the thickness of the samples is decreased by about 0.6 mm. Trends similar to that of axial shrinkage

Figure 3. Impact of the thickness on the water content profiles of sample-B at 50˚C.

Figure 4. Impact of the thickness on the water content profiles of sample-F at 50˚C.

Table 2. Parameters of the descriptive equation of the sample-B Profiles.

Table 3. Parameters of the descriptive equation of the sample-F Profiles.

of the spirulina have been theoretically demonstrated for the alumina gel by Chemkhi (2008) [29] . The profiles make it possible to establish the drying kinetics for various abscissas of x up to a value of the order of 5 mm. These water content profiles will serve as a basis for determining the transport coefficient and its evolution throughout the drying process.

3.2. Drying Kinetics

In Figure 5 and Figure 6 are carried the variations water content at various abscissa as a function of time for both types of samples. The results are classic: the surfaces dry faster than the heart of the product. But, the differences between the curves are quite considerable. In fact, the water discharged from the zones close to the surface is not compensated by the water flowing from the heart towards the surface under the effect of the large gradients of the liquid phase, the motor terms of the flow. The study is carried out for abscissa between 5 mm and 16 mm. Because, out of this interval, the variations of water content and concentration gradient are too small and the experiments could not be exploited.

3.3. Water Mass Flow and Apparent Concentration Gradient

From these profiles above, we have represented on the one hand the evolution of the water transport flow and that of density gradient at different thicknesses in the initial configuration of the sample ( $x\le 16$ mm) as function of time (Figure 7(a) and Figure 7(b)) and secondly their variation as function of the water content (Figure 7(c) and Figure 7(d)). On the sample-B and sample-F, we note that the transport flow and the gradient have a similar appearance at each abscissa. They increase that the further one moves away from the median plane (x = 0). The gradients grow from the drying beginning and decrease towards the process end. For a given abscissa, the water transport flow becomes constant and has some inflection points that reflect the influence of an instantaneous contraction of the material during drying which contributes to migrate the water.

Figure 5. Thickness Influence on the sample-B water content.

Figure 6. Thickness Influence on the sample-F water content.

Decreasing curves indicate that the flow of water carried by filtration from the interior to the exchange surface is not compensated for in the pores. The main interest of the representation of the gradient as a function of time and as a function of the water content lies in the identification for each slice, a critical period and a critical water content (Figures 8(a)-(d)) to from which the contribution of the solid contraction to the transport of water becomes very negligible. We have just chosen to draw the curves of Figure 7(a), Figure 7(b), Figure 8(a) and Figure 8(b) on abscissa less than 13 mm. Because beyond these values, the rigidity of the products deviates more important and the exploitation of the data becomes delicate.

3.4. Relation between Water Flux and Concentration Gradient

In Figure 9 and Figure 10, a relation between the mass flux and the density gradient for each material studied is shown. This function is quasi-linear for each slice considered. What approaches us to Fick’s law which states that in a

Figure 7. (a) Variation of the water flow in sample-B in function of time for various abscissas at 50˚C; (b) Variation of the water flow in sample-F in function of time for various abscissas at 50˚C; (c) Variation of the water flow in sample-B in function of the water content for various abscissas at 50˚C; (d) Variation of the water flow in sample-F in function of the water content for various abscissas at 50˚C.

Figure 8. (a) Water concentration gradient of the sample-B in function of time for various abscissas at 50˚C; (b) Water concentration gradient of the sample-F in function of time for various abscissas at 50˚C; (c) Water concentration gradient in sample-B depending of the water content for various abscissas at 50˚C; (d) Water concentration gradient in sample-F depending of the water content for various abscissas at 50˚C.

Figure 9. Water flow depending of the water concentration gradient in sample-B at 50˚C.

Figure 10. Water flow depending of the water concentration gradient in sample-F at 50˚C.

homogeneous medium, the flow of the liquid phase can be described by a linear relation (Equation (13)).

3.5. Water Transport Coefficient

From the water content profiles, density gradients and fluxes presented above, the range of water content that can be studied is 0.7 < w_{w} < 2.5 kg_{w}/kg_{dm}. That is, from a horizontal line intersecting at least two water content profiles. In this range, the water flow and the apparent density gradient of water are determined by Equations (11) and (12). For each water content value, a relatively linear correlation is observed between the water flow and the gradient. Therefore, the water transport coefficient is determined from the slope of each linear approximation. For both strains (sample-B and sample-F), the dependencies of the water transport coefficients with respect to the water content are presented in Figure 11. It is noted that the experimental results are practically concordant. We can

Figure 11. Water transport coefficient depending of the water content at 50°C: comparison of the experimental results and those of the model.

therefore consider that the two results are representative of a single evolution of the water transport coefficient of Spirulina type platensis. Experimental data are accurately described by an exponential equation (Equation (17)).

${D}_{w}\left(w\right)=2.78\times {10}^{-10}+3.33\times {10}^{-12}\mathrm{exp}\left(3.12w\right)$ (17)

An optimization algorithm based on a nonlinear regression procedure made it possible to obtain a good correlation between the experimental data and those of the theory with a correlation coefficient R^{2} > 97% and a standard deviation RMSE < 1.16 × 10^{−10} m^{2}∙s^{−1}.

There is a sharp decrease in the transport coefficient for water contents above 2 kg_{w}/kg_{dm} and a low decrease for low water contents. This trend has already been observed with rubber [30] , agar gel [31] [32] , clay [33] , coffee [34] and wood [35] . W = 2 kg_{w}/kg_{dm} represents a critical mean value from which the contribution of solid shrinkage to water transport is negligible. In the case of Spiruline platenis, previous work [36] [37] has shown that the solid acquires a mechanical strength and becomes progressively rigid during the last stages of drying. It should be remembered that the proposed destructive method satisfactorily assesses the water transport coefficient of Spirulina platensis over a wide range of water content.

4. Conclusions

This work enabled us to experimentally establish water content profiles during the one-dimensional drying process of Spirulina platensis by cutting samples. These profiles are then used to determine the influence of the water content on the mass flow, the water density gradient and the water transport coefficient that governs the internal phenomena. The distribution of water in the biomass is evaluated with the inclusion of the shrinkage along the cylindrical axis. The determined drying kinetics have similar profiles either in Sample-B or in Sample-F. It is shown that the proposed destructive method satisfactorily evaluates the evolution of the water content of Spirulina platensis over a wide range of water content.

The work carried out therefore provides a prospective tool that may be of interest for the local analysis of the water transport coefficient in other deformable materials. The development of a numerical model is necessary to compare the results and to predict the internal mechanisms during drying.

Nomenclature

RH: relative humidity (%)

db: dry matter (kg_{w}/kg_{dm})

d: cylinder diameter

w: water content (kg/kg)

${a}_{w}$ : water activity (/)

${D}_{w}$ : transport coefficient (m2/s)

${F}_{w}$ : mass water flux (kg∙m^{−2}∙s^{−1})

${G}_{w}$ : gradient of water density (kg/m4)

${m}_{eq}$ : equilibrium mass of sample (kg)

${m}_{i}$ : sample initial mass (kg)

x: sample thickness (m)

${m}_{s}$ : dry mass of sample (kg)

${u}_{s}$ : velocity of the solid phase (m/s)

${u}_{w}$ : velocity of the liquid phase (m/s)

N: number of observations

${R}^{2}$ : correlation coefficient (%)

RMSE: root mean square error (kg/kg)

${w}_{eq}$ : equilibrium moisture content (kg/kg)

Y: given drying parameter

$m\left(t\right)$ : sample mass at a instant (kg)

${Y}_{pre,j}$ : the ith predicted value of Y

${\stackrel{\xaf}{Y}}_{exp,j}$ : mean experimental value of Y

${Y}_{exp,j}$ : given drying parameter Y

Greek Letters

$\alpha $ : real densities ratio (/)

${\rho}_{s}$ : solid apparent density (kg/m3)

${\rho}_{w}$ : water apparent density (kg/m3)

${\rho}_{s}^{*}$ : solid real density (kg/m3)

${\rho}_{w}^{*}$ : water real density (kg/m3)

$\varphi $ : porosity (/)

Subscripts

0: initial

Z: number counter of parameters

i, j: respectively abscissa and time counters

w, db: water, dry matter

eq: equilibrium

Exp: experimental

Pre: predicted

s: solid

1D: unidimensionnel

References

[1] Bouraout, M., Richard, P. and Durance, T. (1994) Microwave and Convective Drying of Potato Slices. Journal of Food Process Engineering, 17, 353-363.

https://doi.org/10.1111/j.1745-4530.1994.tb00343.x

[2] Yongsawatdigul, J. and Gunasekaran, S. (1996) Microwave Vaccum Drying of Cranberries Part II. Journal of Food Processing and Preservation, 20, 145-156.

https://doi.org/10.1111/j.1745-4549.1996.tb00851.x

[3] Feng, H. and Tang, J. (1998) Microwave Finish Drying of Diced Apples in a Spouted Bed. Journal of Food Science, 63, 679-683.

https://doi.org/10.1111/j.1365-2621.1998.tb15811.x

[4] Yang, F., Zhang, M., Mujumdar, A.S., Zhong, Q. and Wang, Z. (2018) Enhancing Drying Efficiency and Product Quality Using Advanced Pretreatments and Analytical Tools: An Overview. Drying Technology, 36, 1824-1838.

https://doi.org/10.1080/07373937.2018.1431658

[5] Feng, Y., Zhang, M., Fan, K. and Mujumdar, A.S. (2018) Effects of Drying Methods on Quality of Fermented Plant Extract Powder. Drying Technology, 36, 1913-1919.

https://doi.org/10.1080/07373937.2018.1461649

[6] Hiep Le, K., Hampel, N., Kharaghani, A., Bück, A. and Tsotsas, E. (2018) Superheated Steam Drying of Single Wood Particles: A Characteristic Drying Curve Model Deduced from Continuum Model Simulations and Assessed by Experiments. Drying Technology, 36, 1866-1881.

https://doi.org/10.1080/07373937.2018.1444633

[7] R1etz, S., Porley, V.E., von Gersdorff, G., Hensel, O., Crichton, S. and Sturm, B. (2017) Effect of Maturation and Freezing on Quality and Drying Kinetics of Beef. Drying Technology, 35, 2002-2014.

https://doi.org/10.1080/07373937.2017.1295051

[8] Tham, T.C., Ng, M.X., Gan, S.H., Chua, L.S., Aziz, R., Chuah, L.A., Hii, C.L., Ong, S.P., Chin, N.L. and Law, C.L. (2017) Effect of Ambient Conditions on Drying of Herbs in Solar Greenhouse Dryer with Integrated Heat Pump. Drying Technology, 35, 1721-1732.

https://doi.org/10.1080/07373937.2016.1271984

[9] Bornhorst, M., Walzel, P., Rahimi, A., Kharaghani, A., Tsotsas, E., Nestle, N., Besser, A., Jager, F.K. and Metzger, T. (2016) Influence of Pore Structure and Impregnation—Drying Conditions on the Solid Distribution in Porous Support Materials. Drying Technology, 34, 1964-1978.

https://doi.org/10.1080/07373937.2016.1147048

[10] Tsuruta, T., Tanigawa, H. and Sashi, H. (2015) Study on Shrinkage Deformation of Food in Microwave-Vacuum Drying. Drying Technology, 33, 1830-1836.

https://doi.org/10.1080/07373937.2015.1036286

[11] Miri Tari, S.M., Habibzade, S. and Taghiyari, H.R. (2015) Effects of Drying Schedules on Physical and Mechanical Properties in Paulownia Wood. Drying Technology, 33, 1981-1990.

https://doi.org/10.1080/07373937.2014.948553

[12] Bénet, J.-C., Ramirez-Martinez, A., Ouedraogo, F. and Cherblanc, F. (2012) Measurement of the Chemichal Potential of a Liquid in Porous Media. Journal of Porous Media, 15, 1019-1029.

https://doi.org/10.1615/JPorMedia.v15.i11.30

[13] Henrikson, R. (2000) Earthfood Spirulina. Earthrise, Ronore Enterprises.

[14] Maskan, M. (2001) Kinetics of Colour Change of Kiwifruits during Hot Air and Microwave Drying. Journal of Food Engineering, 48, 169-175.

https://doi.org/10.1016/S0260-8774(00)00154-0

[15] Gowen, A.A., Abu-Ghannam, N., Frias, J. and Oliveira, J. (2008) Modeling Dehydratation and Rehydratation of Cooked Soybeans Subjected to Combined Microwave-Hot-Air Drying. Innovative Food Science and Emerging Technologies, 9, 129-137.

https://doi.org/10.1016/j.ifset.2007.06.009

[16] Bourassa, J., Ramachandran, R.P., Paliwal, J. and Cenkowski, S. (2015) Drying Characteristics and Moisture Diffusivity of Distillers’ Spent Grains Dried in Superheated Steam. Drying Technology, 33, 2012-2018.

https://doi.org/10.1080/07373937.2015.1040883

[17] Compaore, A., Dissa, A.O., Rogaume, Y., Putranto, A., Chen, X.-D., Mangindaan, D., Zoulalian, A., Rémond, R. and Tiendrebeogo, E. (2016) Application of the Reaction Engineering Approach (REA) for Modeling of the Convective Drying of Onion. Drying Technology, 35, 500-508.

https://doi.org/10.1080/07373937.2016.1192189

[18] Tsarahevitra, J., Loc, C. and Marie-Jose, L. (2004) Adaptation des souches de spiruline du sud de madagascar à la culture en eau de mer. In: Colloque international sur les cyanobactéries pour la santé, la science et le développement, Ile des Embiez.

[19] Enzing, C., Ploeg, M., Barbosa, M. and Sijtsma, L. (2014) Microalgae Based Products for the Food and Feed Sector: An Outlook for Europe. IPTS Institute for Prospective Technological Studies, JRC.

[20] Dissa, A., Desmorieux, H., Savadogo, P., Segda, B.G. and Koulidiati, J. (2010) Shrinkage, Porosity and Density Behaviour during Convective Drying of Spirulina. Journal of Food Engineering, 97, 410-418.

https://doi.org/10.1016/j.jfoodeng.2009.10.036

[21] Bonazzi, C., Dumoulin, E. and Bimbenet, J.-J. (2008) Le séchage des produits alimentaires. Industrie Alimentaire Agricole, 125, 12-22.

[22] Desmorieux, H. and Hernandez, F. (2004). Biochemical and Physical Criteria of Spirulina after Different Drying Processes. In: Proceedings of the 14th International Drying Symposium (IDS 2004), Sào Paulo, Brazil.

[23] Desmorieux, H. and Decaen, N. (2005) Convective Drying of Spirulina in Thin Layer. Journal of Food Engenering, 66, 497-503.

https://doi.org/10.1016/j.jfoodeng.2004.04.021

[24] Tiendrebeogo, E.S., Dissa, A.O., Cherblanc, F., Youm, I., Bénet, J.-C., Compaoré, A. and Koulidiati, J. (2015) Characterization of Two Different Stumps of Spirulina platensis Drying: Assessment of Water Transport Coefficient. Food and Nutrition Sciences, 6, 1437-1449.

https://doi.org/10.4236/fns.2015.615148

[25] Mihoubi, D., Zagrouba, F., Vaxelaire, J., Bellagi, A. and Roques, M. (2004) Transfer Phenomena during the Drying of a Shrinkable Product, Modeling and Simulation. Drying Technology, 22, 91-109.

https://doi.org/10.1081/DRT-120028216

[26] Hawladsr, M., Ho, J. and Qing, Z. (1999) A Mathematical Model for Drying of Shrinking Materials. Drying Technology, 17, 27-47.

https://doi.org/10.1080/07373939908917517

[27] Couture, F., Laurent, S. and Roques, M.A. (2007) Drying of Two-Phase Media: Simulation with Liquid Pressure as Driven Force. AIChE Journal, 53, 1703-1717.

https://doi.org/10.1002/aic.11213

[28] Crank, J. (1975) The Mathematics of Diffusion. Oxford University Press, London.

[29] Chemkhi, S., Zagrouba, F. and Bellagi, A. (2004) Mathematical Model for Drying of Highly Shrinkable Media. Drying Technology, 22, 1023-1039.

https://doi.org/10.1081/DRT-120038578

[30] Auria, R. and Bénet, J.-C. (1990) Transport de l’eau dans une feuille de caoutchouc naturel pendant la période de séchage à vitesse décroissante. International Journal of Heat and Mass Transfer, 33, 1885-1894.

https://doi.org/10.1016/0017-9310(90)90220-O

[31] Mrani, I., Bénet, J. and Fras, G. (1995) Transport of Water in a Biconstituent Elastic Medium. Applied Mechanics Review, 48, 717-721.

https://doi.org/10.1115/1.3005053

[32] Mrani, I. (1993) Transport de matière en milieu biconstituant élastique: application au séchage du gel d’agar. Ph.D. Thesis, Université Montpellier, Montpellier, France.

[33] Ketelaars, A., Pel, L., Coumans, W. and Kerkhof, P. (1995) Drying Kinetics: A Comparison of Diffusion Coefficients from Moisture Concentration Profiles and Drying Curves. Chemical Enaineerir Science, 50, 1187-1191.

https://doi.org/10.1016/0009-2509(94)00494-C

[34] Ramirez-Martinez, A., Salgado-Cervantes, M., Rodriguez-Jimenes, G., Garcia-Alvarado, M., Cherblanc, F. and Bénet, J.-C. (2013) Water Transport in Parchment and Endosperm of Coffee Bean. Journal of Food Engineering, 114, 375-383.

https://doi.org/10.1016/j.jfoodeng.2012.08.028

[35] Anoua, M., Ramirez-Martinez, A., Cherblanc, F. and Bénet, J.-C. (2014) The Use of Chemical Potential to Describe Water Transfer in Complex Media with Strong Solid-Liquid Bonding. Transport in Porous Media, 102, 111-122.

https://doi.org/10.1007/s11242-013-0265-4

[36] Dissa, A., Compaore, A., Tiendrebeogo, E. and Koulidiati, J. (2014) An Effective Moisture Diffusivity Model Deduced from Experiment and Numerical Solution of Mass Transfer Equations for a Shrinkable Drying Slab of Microalgae Spirulina. Drying Technology, 32, 1231-1244.

https://doi.org/10.1080/07373937.2014.897234

[37] Tiendrebeogo, S.E. and Dissa, A.O. (2015) Guy Christian Tubreoumya, Kayaba Haro, Fabien Cherblanc, J Koulidiati, Couplage des mécanismes de déshydratation et du retrait de la Spiruline platensis “Arthrospira platensis” lors d’un séchage isotherme. Journal de la Société Ouest-Africaine de Chimie, 040, 17-23.