<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">EPE</journal-id><journal-title-group><journal-title>Energy and Power Engineering</journal-title></journal-title-group><issn pub-type="epub">1949-243X</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/epe.2024.162004</article-id><article-id pub-id-type="publisher-id">EPE-131339</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Optimization of the Pretreatment of the Mixture of Cassava Peelings and Pineapple Fibers Using Response Surface Methodology and a Process Simulator for the Bioethanol Production
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Paul</surname><given-names>Nestor Djomou Djonga</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>George</surname><given-names>Elambo Nkeng</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Madjoyogo</surname><given-names>Hervé Sirma</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ahmat</surname><given-names>Tom</given-names></name><xref ref-type="aff" rid="aff4"><sup>4</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Thierry</surname><given-names>Tchamba Tchuidjang</given-names></name><xref ref-type="aff" rid="aff5"><sup>5</sup></xref></contrib></contrib-group><aff id="aff3"><addr-line>Institut du Génie des Systèmes Industriels et textiles, Ecole Polytechnique de Ouagadougou, Ouagadougou, Burkina-Faso</addr-line></aff><aff id="aff1"><addr-line>Department of Textile and Leather Engineering, National Advanced School of Engineering of Maroua, Maroua, Cameroon</addr-line></aff><aff id="aff5"><addr-line>Department of Process Engineering, National Advanced School of Agro-Industrial Sciences, University of Ngaoundere, Ngaoundere, Cameroon</addr-line></aff><aff id="aff4"><addr-line>Department of Energy Engineering, University Institute of Technology, University of Ngaoundere, Ngaoundere, Cameroon</addr-line></aff><aff id="aff2"><addr-line>Department of Civil, Environmental and Architectural Engineering, University of Padova, Padova, Italy</addr-line></aff><pub-date pub-type="epub"><day>26</day><month>02</month><year>2024</year></pub-date><volume>16</volume><issue>02</issue><fpage>79</fpage><lpage>96</lpage><history><date date-type="received"><day>8,</day>	<month>September</month>	<year>2023</year></date><date date-type="rev-recd"><day>24,</day>	<month>February</month>	<year>2024</year>	</date><date date-type="accepted"><day>27,</day>	<month>February</month>	<year>2024</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution-NonCommercial International License (CC BY-NC).http://creativecommons.org/licenses/by-nc/4.0/</license-p></license></permissions><abstract><p>
 
 
  The increase in oil prices and greenhouse gas emissions has led to the search for substitutes for fossil fuels. In Cameroon, the abundance of lignocellulosic resources is inherent to agricultural activity. Production of bioethanol remains a challenge given the crystallinity of cellulose and the presence of the complex. The pretreatment aimed to solubilize the lignin fraction and to make cellulose more accessible to the hydrolytic enzymes, was done using the organosolv process. A mathematical modeling was performed to point out the effect of the temperature on the kinetics of the release of the reducing sugars during the pretreatment. Two mathematical model was used, SAEMAN’s model and Response surface methodology. The first show that the kinetic parameters of the hydrolysis of the cellulose and reducing sugar are: 0.05089 min
  <sup>-1</sup>, 5358.1461 J
  &#183;mol
  <sup>-1</sup>, 1383.03691 min
  <sup>-1</sup>, 51577.6100 J
  &#183;mol
  <sup>-1</sup> respectively. The second model was used. Temperature is the factor having the most positive influence whereas, ethanol concentration is not an essential factor. To release the maximum, an organosolv pre-treatment of this sub-strate should be carried out at 209.08&#176;C for 47.60 min with an ethanol-water ratio of 24.02%. Organosolv pre-treatment is an effective process for delignification of the lignocellulosic structure.
 
</p></abstract><kwd-group><kwd>Bioethanol</kwd><kwd> Cassava Peeling</kwd><kwd> Pineapple Fibers</kwd><kwd> Organosolv Process and Optimization</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Energy is an engine of growth for any society. In Third World Countries, particularly in the Sub-Saharan zone, household energy coverage remains a luxury [<xref ref-type="bibr" rid="scirp.131339-ref1">1</xref>] . The biggest concern of households in our different municipalities remains access to stable, quality energy, respectful of the environment and above all at lower cost. To survive, many households in developing countries still rely on non-renewable sources of energy, the consumption of which increases over time. In Cameroon for example, it is estimated that wood consumption for energy supply will increase at a rate of around 3% per year. Most of the increase in wood is represented by firewood (91.2%) and charcoal (1%), the predominant energy sources in rural areas [<xref ref-type="bibr" rid="scirp.131339-ref2">2</xref>] . Overall, around 80% of Cameroonians and almost 100% of rural households depend on woody biomass as the main and only source of energy [<xref ref-type="bibr" rid="scirp.131339-ref3">3</xref>] . It’s skyrocketing Deforestation. In Cameroonian communities where forests are rapidly disappearing, the use of firewood as a major energy the source has become expensive. The Cameroonian economy is essentially based on agriculture, which produces nearly $60 million in tons of different crops per year. This represents an enormous volume of agricultural residues which could be exploited to meet the country’s energy challenges. Using fuelwood as a key energy source is time consuming and labor intensive and harmful to the environment [<xref ref-type="bibr" rid="scirp.131339-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref3">3</xref>] . Biogas, which is a clean technology, has been considered in recent years as a promising source of alternative energy suitable for fuelwood, for in many other African countries in general and for rural households in Cameroon in particular. The world’s stock of economically exploitable fossil Energy (oil, gas, and coal) are limited and the number of people in the world is estimated at 7.5 billion today and 10 billion in 2050. In these circumstances, it is essential to prepare the succession of the combustible fossil towards renewable forms of energy based on organic raw material. The development of second-generation biofuels based on the conversion of cellulosic biomass into bioethanol by enzymatic hydrolysis or into biodiesel by thermochemical gasification or liquefaction by catalytic conversion [<xref ref-type="bibr" rid="scirp.131339-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref4">4</xref>] . Generally, agro food industries generate significant quantities of lignocellulosic waste which is poorly exploited. Also, the increase in the price of oil and the increase in greenhouse gas emissions justify the search for raw materials and alternative technologies capable of reducing dependence on these fossil fuels and protecting the environment. Bioethanol is already used to substitute a fraction of the volume of gasoline in many countries. Several authors have worked on the conversion of biomass into bioethanol using the organosolv processes. It emerges from this work that the organosolv process is one of the best because it makes it possible to delignify the biomass without denaturing the other compounds, while ethanol and lignin can be recovered on the other hand [<xref ref-type="bibr" rid="scirp.131339-ref4">4</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref6">6</xref>] . But to the best of our knowledge, there is no work on the production of bioethanol from association of cassava waste and pineapple fibers waste pre-treated by the organosolv process.</p><p>The high cellulose and lignin contents are a hindrance to the biochemical conversion of lignocellulosic materials into liquid and gaseous bioenergy (bioethanol and biogas), due to the crystalline structure of cellulose, and the rigidity of lignin, which act as obstacles to the enzymatic hydrolysis of fermentable sugars [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] . A pretreatment is therefore necessary to separate the lignin from the polysaccharides, to increase the specific surface area and the porosity of the cellulose, to decrease the degree of polymerization of the cellulose [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] . There are several types of pretreatment: physical, chemical, biological and physico-chemical, but all of them release more or less inhibitory compounds for the subsequent biochemical steps [<xref ref-type="bibr" rid="scirp.131339-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] . The organosolv process is the extraction of lignin using an organic solvent [<xref ref-type="bibr" rid="scirp.131339-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref10">10</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref12">12</xref>] . It gives the best yield of lignin solubilization [<xref ref-type="bibr" rid="scirp.131339-ref9">9</xref>] . Several works have focused on the non-mathematical modeling of pretreatments, in particular the organosolv process [<xref ref-type="bibr" rid="scirp.131339-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref16">16</xref>] , but very few have focused on mathematical modeling and optimization, hence the objective of our work which is to mathematically model and to optimize the pretreatment of cocoa shells by the organosolv process. More specifically, we will look at the influence of temperature on the kinetics of pretreatment and determine the optimal temperature profile as a function of time.</p></sec><sec id="s2"><title>2. Materials and Methods</title><sec id="s2_1"><title>2.1. Substrate Sampling and Characterization</title><p>The cassava and pineapple skins used were collected at the Yaound&#233; market landfill, then they were transported to the laboratory, washed with distilled water, dried in an oven and ground until a fine consistency was obtained particle size of less than 80 microns. Once in the laboratory, the mixture 50:50 of raw material underwent a physico-chemical characterization through the determination of its dry matter (MS) and its mineral matter (MM) according to the calcination method. The lignin content was determined, the hemicellulose content and the cellulose content.</p></sec><sec id="s2_2"><title>2.2. Substrate Pre-Treatments</title><p>Organosolv pre-treatment of sawdust is necessary to release the cellulose from the lignin-hemicellulose complex. It also makes it possible to reduce cellulosic crystallinity and makes the cellulose accessible to enzymes during the hydrolysis step. In order to better control this pre-treatment step, modelling by the response surface methodology was used. This is how the centred composite plan was used to determine the optimal levels of the various factors influencing treatment [<xref ref-type="bibr" rid="scirp.131339-ref6">6</xref>] .</p></sec><sec id="s2_3"><title>2.3. Choice of Factors and Operating Mode</title><p>The parameters of the organosolv treatment are the temperature, the treatment time, the liquid/solid ratio, the concentration of the solvent (ethanol) and this concentration is expressed in terms of the liquid/solid ratio generally fixed at 10:1 [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] . The optimization factors retained for the experimental device of the organosolv treatment are: the temperature (˚C), the treatment time (min).</p><p>The organosolv pretreatment consists in introducing 3 g of shell into 60 ml of an ethanol solution at a concentration indicated by the experimental matrix, the whole is brought to the set temperature and for a period indicated by the experimental matrix. The warm-up time was between 15 - 20 minutes. The treatment was stopped by quenching and the liquid fraction was recovered using a filtration device. The content of total cellulose and phenolic compounds were determined in the liquid fraction.</p><p>The optimization carried out as part of this work will be carried out by two approaches. The first will consist of developing a mathematical model and the second an experimental model based on the experimental plans.</p></sec><sec id="s2_4"><title>2.4. Mathematical Modeling and Optimization</title><p>➢ Mathematical modeling</p><p>The model chosen is that of SAEMAN which states that the hydrolysis reactions of cellulose and glucose are of first order and the rate constants are a function of the acid concentration and temperature [<xref ref-type="bibr" rid="scirp.131339-ref17">17</xref>] .</p><p>d C i d t = ( − 1 0 1 − 1 ) ⋅ ( k 1 ⋅ C c e l l u l o s e k 2 ⋅ C s u c r e s r e d ) (1)</p><p>With i the compound (cellulose, glucose, hydroxy-methyl furfural), k<sub>1</sub> and k<sub>2</sub> the hydrolysis constants of cellulose and glucose respectively.</p><p>The rate constants are related to the temperature by the Arrhenius relation:</p><p>k 1 = A 1 &#215; exp ( − E 1 R &#215; T ) (2)</p><p>k 2 = A 2 &#215; exp ( − E 2 R &#215; T ) (3)</p><p>&#183; Parameter estimation</p><p>The parameters A<sub>1</sub>, A<sub>2</sub>, E<sub>1</sub>, E<sub>2</sub> are estimated by nonlinear regression with the experimental data using the Levenberg-Marquartd algorithm to solve the least squares problem.</p><p>&#183; Model validation</p><p>The model goodness of fit was expressed by the correlation coefficient (R<sup>2</sup>).</p><p>➢ Optimization</p><p>The chosen optimization technique used is optimal control which is an application of dynamic continuous-time optimization [<xref ref-type="bibr" rid="scirp.131339-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref21">21</xref>] .</p><p>The state variables are the cellulose and glucose concentrations, the control variable is the temperature of the reactor, the adjunct variables p<sub>1</sub> and p<sub>2</sub>. Thus, we have the problem formulated as follows:</p><p>Concentration rate: d C i d t = ( − 1 0 1 − 1 ) &#215; ( k 1 &#215; C c e l l u l o s e k 2 &#215; C s u c r e s r e d )</p><p>Objective function or cost: d C i d t = ( − 1 0 1 − 1 ) &#215; ( k 1 &#215; C c e l l u l o s e k 2 &#215; C s u c r e s r e d )</p><p>Hamiltonian: H = k 2 &#215; C s − k 1 &#215; C c − p 1 &#215; k 1 &#215; C c + p 2 &#215; ( − k 2 &#215; C s + k 1 &#215; C c )</p><p>The model used in this work is similar to those used by [<xref ref-type="bibr" rid="scirp.131339-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref19">19</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref20">20</xref>] .</p></sec><sec id="s2_5"><title>2.5. Choice of Factors and Experience Matrix</title><p>The temperature ranges during organosolv treatment are between 171 and 207˚C for pretreatment Tchuidjjang et al. [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] and Carrasco [<xref ref-type="bibr" rid="scirp.131339-ref22">22</xref>] . The time range (15 - 45 min) was chosen based on the work of Kabir [<xref ref-type="bibr" rid="scirp.131339-ref17">17</xref>] and Kabel et al. [<xref ref-type="bibr" rid="scirp.131339-ref23">23</xref>] . The ethanol concentration between 30 and 70% v/v was set based on the work of Klason P [<xref ref-type="bibr" rid="scirp.131339-ref24">24</xref>] and N’Diaye [<xref ref-type="bibr" rid="scirp.131339-ref25">25</xref>] . <xref ref-type="table" rid="table1">Table 1</xref> presents the experimental matrix for the pretreatment of the organosolv method.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Matrix of real and coded variables of the centred composite plane of organosolv pre-treatment</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Experience</th><th align="center" valign="middle"  colspan="3"  >Coded variables</th><th align="center" valign="middle"  colspan="3"  >Real variables</th></tr></thead><tr><td align="center" valign="middle" >Temperature</td><td align="center" valign="middle" >Time</td><td align="center" valign="middle" >Concentration of ethanol</td><td align="center" valign="middle" >Temperature</td><td align="center" valign="middle" >Time</td><td align="center" valign="middle" >Concentration of ethanol</td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >30</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >70</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >70</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >−1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >70</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >70</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >1.41</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >218.2</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >−1.41</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >161.8</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1.41</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >161.8</td><td align="center" valign="middle" >56.15</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >−1.41</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >13.85</td><td align="center" valign="middle" >50</td></tr><tr><td align="center" valign="middle" >17</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >1.41</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >78.2</td></tr><tr><td align="center" valign="middle" >18</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >0</td><td align="center" valign="middle" >−1.41</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >21.8</td></tr></tbody></table></table-wrap><p>The mathematical model postulated is of the second-degree polynomial type with interactions. Its expression is in the form of Equation (1).</p><p>Y = a 0 + ∑ i = 1 k a i x i + ∑ i 2 = 1 k a i i x i 2 + ∑ i = 1 k ∑ j = 1 k a i j x i x j + E (4)</p><p>Y = concentration of reducing sugars, pentoses, phenolic compounds, dependent variable (response)</p><p>x<sub>i</sub> et x<sub>j</sub> = dependent variable; a<sub>0</sub> = constante; E = error</p><p>a<sub>i</sub> = linear coefficient; a<sub>ii</sub> = quadratic coefficient; a<sub>ij</sub> = interaction coefficient</p><p>The quality of fit of the model was expressed by the coefficients of determination R<sup>2</sup> and adjusted R<sup>2</sup> which must be less than 80 [<xref ref-type="bibr" rid="scirp.131339-ref26">26</xref>] . As this mathematical tool was not sufficient to validate a model, other tools and methods were used to judge the adequate quality of the models. To this end, the absolute mean deviation analysis (AADM) between the calculated and observed values [<xref ref-type="bibr" rid="scirp.131339-ref27">27</xref>] , the bias factor (Bf) [<xref ref-type="bibr" rid="scirp.131339-ref28">28</xref>] and the accuracy factor (Af) [<xref ref-type="bibr" rid="scirp.131339-ref29">29</xref>] were carried out. They were calculated from Equations (2)-(4).</p><p>➢ Bias Factor (Bf) Bf = 10<sup>B</sup> avec</p><p>B = 1 n ∑ ​ log Y t h e Y o b s (5)</p><p>➢ Accuracy Factor (Af) Af = 10<sup>A</sup> avec</p><p>A = 1 n ∑ ​ | log Y t h e Y o b s | (6)</p><p>➢ Absolute Mean Deviation Analysis (AADM)</p><p>AADM = ∑ i = 1 n | Y o b s − Y t h e Y o b s | n (7)</p><p>Y<sub>the</sub> and Y<sub>obs</sub> are respectively the experimental response and calculated from the model for an experiment i; n being the total number of experiments.</p></sec><sec id="s2_6"><title>2.6. Pre-Treatment Operation</title><p>The organosolv a pretreatment consists of taking 3 g of mixture of pineapple fiber powder and cassava skin in a ratio of 50:50 and introducing into 75 ml of an ethanol solution as indicated in the test matrix. The assembly was brought to the set temperature and for the duration indicated by the methods described in <xref ref-type="table" rid="table1">Table 1</xref>. The duration of the pretreatment operation is taken when the temperature sensor displays the set point imposed by the matrix. Warm-up time is 10 - 15 minutes. The reaction ends by immersing the water at a temperature of 10˚C and the filtrate will be recovered for subsequent analyses.</p></sec><sec id="s2_7"><title>2.7. Characterization of the Pre-Treated Substrate</title><p>The concentration of reducing sugars was determined by the spectroscopic method of Godin B [<xref ref-type="bibr" rid="scirp.131339-ref28">28</xref>] , Fisher and Stein [<xref ref-type="bibr" rid="scirp.131339-ref30">30</xref>] . Total phenolic compounds were determined by the Folin-Ciocalteu Hashimoto et al. [<xref ref-type="bibr" rid="scirp.131339-ref31">31</xref>] method with some modifications. The pentose concentration of the pretreatment hydrolysates was determined by the method of Akpinar O. [<xref ref-type="bibr" rid="scirp.131339-ref32">32</xref>] .</p></sec><sec id="s2_8"><title>2.8. Optimization</title><p>The optimization which consisted in seeking the optimal conditions allowing to have the maximum of cellulose and glucose was carried out by the software Minitap 17.</p></sec></sec><sec id="s3"><title>3. Results and Discussion</title><sec id="s3_1"><title>3.1. Mathematical Modeling of the Pretreatment of Raw Material Mixture by Simulation</title><p>Modeling made it possible to study the influence of temperature on the kinetics of the pretreatment, for a powder with a particle size of 1 millimeter, a concentration of sulfuric acid of 1% or approximately 0.0914M and a pH = 1.039, a concentration of ethanol fixed at 75% (v/v) [<xref ref-type="bibr" rid="scirp.131339-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref33">33</xref>] . Figures 1-3 present the result of Pretreatment kinetics processing at 170˚C, 185˚C and 200˚C.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Results of modeling and statistical validation criteria</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Compounds</th><th align="center" valign="middle" >Kinetic models</th><th align="center" valign="middle" >T = 170˚C</th><th align="center" valign="middle" >T = 185˚C</th><th align="center" valign="middle" >T = 200˚C</th></tr></thead><tr><td align="center" valign="middle" >Cellulose (C)</td><td align="center" valign="middle" >d C d t = − k 1 &#215; C</td><td align="center" valign="middle" >k<sub>1</sub> = 0.0119 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9992 RMSE = 0.0025</td><td align="center" valign="middle" >k<sub>1</sub> = 0.0125 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9988 RMSE = 0.0013</td><td align="center" valign="middle" >k<sub>1</sub> = 0.0130 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9965 RMSE = 0.0002</td></tr><tr><td align="center" valign="middle" >Reducing sugars (G)</td><td align="center" valign="middle" >d G d t = k 1 &#215; C − k 2 &#215; G</td><td align="center" valign="middle" >k<sub>2</sub> = 0.0011 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9933 RMSE = 0.0043</td><td align="center" valign="middle" >k<sub>2</sub> = 0.0019 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9954 RMSE = 0.0036</td><td align="center" valign="middle" >k<sub>2</sub> = 0.0028 min<sup>−</sup><sup>1</sup> R<sup>2</sup>adj = 0.9963 RMSE = 0.0003</td></tr></tbody></table></table-wrap><p>With RMSE and R<sup>2</sup>adj: square root of the mean of the square of the error and Adjusted coefficient of determination respectively.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Pre-exponential factors and activation energies</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Pre-Exponential Factor (min<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >Activation Energy (J∙mol<sup>−</sup><sup>1</sup>)</th><th align="center" valign="middle" >Radj<sup>2</sup></th><th align="center" valign="middle" >P value</th></tr></thead><tr><td align="center" valign="middle" >A<sub>1</sub> = 0.05089</td><td align="center" valign="middle" >E<sub>1</sub> = 5358.1461</td><td align="center" valign="middle" >1</td><td align="center" valign="middle" >0.000526</td></tr><tr><td align="center" valign="middle" >A<sub>2</sub> = 1383.03691</td><td align="center" valign="middle" >E<sub>2</sub> = 51577.6100</td><td align="center" valign="middle" >0.972</td><td align="center" valign="middle" >0.0403</td></tr></tbody></table></table-wrap><p>The rate constants for the hydrolysis of cellulose and glucose are given in <xref ref-type="table" rid="table2">Table 2</xref> above.</p><p>The results of the modeling show that the kinetics of the pre-treatment by the organosolv process are of order 1, the increase in temperature leads to that of the rate constants, which is in accordance with the work of Thierry Tchamba Tchuidjang, Eric Noubissi&#233;, Ahmed Ali [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] .</p><p>The kinetic parameters of the kinetic constants, Equations (2) and (3) are presented in <xref ref-type="table" rid="table3">Table 3</xref>.</p><p>The results of the modeling show that the kinetics of the pretreatment by the organosolv process is of order 1, and that the influence of temperature on the kinetics obeys Arrhenius’ law, which is in accordance with the results of Blaise (2013) Kupiainen L et al., and Kanchanalai P. et al. [<xref ref-type="bibr" rid="scirp.131339-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref20">20</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref33">33</xref>] . These results are consistent with those obtained by authors such as Park N. et al. and Fisher E.H., Stein E.A. [<xref ref-type="bibr" rid="scirp.131339-ref30">30</xref>] [<xref ref-type="bibr" rid="scirp.131339-ref34">34</xref>] in that the pre-exponential factors and the activation energies of the hydrolysis constant of reducing sugars are greater than those of the constant of hydrolysis of cellulose. However, Kanchanalai P. et al. [<xref ref-type="bibr" rid="scirp.131339-ref33">33</xref>] have found value pre-exponential factor and much higher activation energy for the constants of cellulose hydrolysis rates and reducing sugars for other agricultural residues lignocellulosic, these differences may be due to the absence of agitation in the pretreatment reactor, or to the nature of the lignocellulosic material. <xref ref-type="fig" rid="fig4">Figure 4</xref> presents the result of Simulation of the SAEMAN model of the pretreatment by the organosolv process.</p><sec id="s3_1_1"><title>3.1.1. Mathematical Simulation</title><p>According to <xref ref-type="fig" rid="fig4">Figure 4</xref>, the concentration of cellulose decreases with time and temperature while that of reducing sugar (glucose) first increases whatever the temperature this is due to the fact that the rate of formation of reducing sugar (glucose) is greater than that of its hydrolysis, then it reaches a maximum, which implies that the rates of formation and hydrolysis are equal, and finally it decreases with time and temperature, this may be due to Maillard reactions between proteins and glucose at high temperatures. The work of Kenmogne Yony Blaise [<xref ref-type="bibr" rid="scirp.131339-ref18">18</xref>] led to the same observations and attributed them to the increase in the rate of degradation of glucose by the increase in temperature. It should be noted that Hashimoto S. et al. [<xref ref-type="bibr" rid="scirp.131339-ref31">31</xref>] has found result diametrically opposite, it has been found that the activation energy and the pre-exponential factor of the hydrolysis of cellulose is higher than that of reducing sugars, and thus increasing the temperature increased the maximum concentration of reducing sugars.</p></sec><sec id="s3_1_2"><title>3.1.2. Optimization of the Pretreatment of Raw Material</title><p>The optimization carried out with the optimal control, gave the following results:</p><p>This <xref ref-type="fig" rid="fig5">Figure 5</xref> shows that the optimal temperature drops as the reaction progresses and reaches a plateau after 200 minutes. <xref ref-type="fig" rid="fig6">Figure 6</xref> shows the evolution of cellulose model and reducing sugar during processing. This confirms the results obtained by numerical simulation which shows that the maximum glucose concentration decreases with temperature. This could be because the kinetic parameters of cellulose hydrolysis rate constant are much lower than those of glucose, so the rate of degradation of reducing sugars increases faster than that of hydrolysis of cellulose at high temperatures. The work of Hashimoto S. [<xref ref-type="bibr" rid="scirp.131339-ref31">31</xref>] led to the same observation.</p><p><xref ref-type="fig" rid="fig7">Figure 7</xref> below shows the evolution of the objective function as a function of the number of iterations. It appears that after each iteration, the objective function is improved, and this for decreasing temperatures, which validates the observations</p><p>made during the simulation of the SAEMAN model with the estimated kinetic parameters.</p></sec></sec><sec id="s3_2"><title>3.2. Mathematical Models of the Different Experimental Responses</title><p>The experimental design responses are presented in <xref ref-type="table" rid="table4">Table 4</xref>.</p><p>Three second-degree polynomial mathematical models describing the extraction of the different responses during the organosolv pre-treatment of mixture raw material were generated by this centred composite plane. The content of Cellulose (Y<sub>C</sub>) materialized by Equation (5) is a measure which reflects the solubilization rate of the hemi cellulosic fraction. The content of Reducing Sugars (Y<sub>RS</sub>) materialized by Equation (6) reflects the solubilization of celluloses. The validation indicators for these three mathematical models are presented in <xref ref-type="table" rid="table5">Table 5</xref>.</p><p>Y<sub>C</sub> = content of celluloses; Y<sub>RS</sub> = content of reducing sugars;</p><p>Y C = 2.925 + 0.338 X 1 − 0.337 X 2 − 0.227 X 3 − 0.094 X 1 2 − 0.248 X 2 2     + 0.969 X 3 2 + 0.263 X 1 X 2 − 0.564 X 1 X 3 + 0.270 X 2 X 3 (8)</p><p>Y R S = 2.049 + 0.2062 X 1 + 0.4094 X 2 + 0.1901 X 3 − 0.1747 X 1 2 + 0.4271 X 2 2   − 0.1611 X 3 2 + 0.0440 X 1 X 2 − 0.1546 X 1 X 3 − 0.1952 X 2 X 3 (9)</p><p>The two mathematical models are valid under the basis of the intervals provided by the validation criteria presented in <xref ref-type="table" rid="table3">Table 3</xref>.</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref> shows the response of the Cellulose extraction as a function of factors taken two by two.</p><p>It appears from this work that increasing the temperature favors the yield of cellulose while a long treatment time (36 minutes) will favor a low yield of obtaining cellulose (<xref ref-type="table" rid="table3">Table 3</xref>). The response surfaces perfectly illustrate these observations in <xref ref-type="fig" rid="fig8">Figure 8</xref>(a) and <xref ref-type="fig" rid="fig8">Figure 8</xref>(b). It also appears that temperature is</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Experimental and calculated responses of the centred composite plan</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle"  colspan="3"  >Real variables</th><th align="center" valign="middle"  colspan="2"  >Cellulose (mg/mL)</th><th align="center" valign="middle"  colspan="2"  >Reducing Sugars (mg/mL)</th></tr></thead><tr><td align="center" valign="middle" ></td><td align="center" valign="middle" >X<sub>1</sub></td><td align="center" valign="middle" >X<sub>2</sub></td><td align="center" valign="middle" >X<sub>3</sub></td><td align="center" valign="middle" >Y<sub>ExpP</sub></td><td align="center" valign="middle" >Y<sub>CalP</sub></td><td align="center" valign="middle" >Y<sub>ExpSR</sub></td><td align="center" valign="middle" >Y<sub>CalSR</sub></td></tr><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >3.750</td><td align="center" valign="middle" >2.925</td><td align="center" valign="middle" >0.094</td><td align="center" valign="middle" >2.049</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >2.700</td><td align="center" valign="middle" >2.925</td><td align="center" valign="middle" >0.067</td><td align="center" valign="middle" >2.049</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >2.938</td><td align="center" valign="middle" >2.925</td><td align="center" valign="middle" >0.073</td><td align="center" valign="middle" >2.049</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >3.869</td><td align="center" valign="middle" >2.925</td><td align="center" valign="middle" >0.097</td><td align="center" valign="middle" >2.049</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >4.450</td><td align="center" valign="middle" >3.747</td><td align="center" valign="middle" >0.111</td><td align="center" valign="middle" >1.029</td></tr><tr><td align="center" valign="middle" >6</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >5.080</td><td align="center" valign="middle" >5.024</td><td align="center" valign="middle" >0.127</td><td align="center" valign="middle" >1.662</td></tr><tr><td align="center" valign="middle" >7</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >2.045</td><td align="center" valign="middle" >2.007</td><td align="center" valign="middle" >0.051</td><td align="center" valign="middle" >2.150</td></tr><tr><td align="center" valign="middle" >8</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >30</td><td align="center" valign="middle" >4.641</td><td align="center" valign="middle" >4.337</td><td align="center" valign="middle" >0.116</td><td align="center" valign="middle" >2.960</td></tr><tr><td align="center" valign="middle" >9</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >70</td><td align="center" valign="middle" >2.356</td><td align="center" valign="middle" >3.881</td><td align="center" valign="middle" >0.059</td><td align="center" valign="middle" >2.109</td></tr><tr><td align="center" valign="middle" >10</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >20</td><td align="center" valign="middle" >70</td><td align="center" valign="middle" >2.942</td><td align="center" valign="middle" >2.902</td><td align="center" valign="middle" >0.074</td><td align="center" valign="middle" >2.124</td></tr><tr><td align="center" valign="middle" >11</td><td align="center" valign="middle" >170</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >70</td><td align="center" valign="middle" >3.642</td><td align="center" valign="middle" >3.219</td><td align="center" valign="middle" >0.091</td><td align="center" valign="middle" >2.449</td></tr><tr><td align="center" valign="middle" >12</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >70</td><td align="center" valign="middle" >2.668</td><td align="center" valign="middle" >3.293</td><td align="center" valign="middle" >0.067</td><td align="center" valign="middle" >2.640</td></tr><tr><td align="center" valign="middle" >13</td><td align="center" valign="middle" >218.2</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >3.407</td><td align="center" valign="middle" >3.214</td><td align="center" valign="middle" >0.085</td><td align="center" valign="middle" >1.992</td></tr><tr><td align="center" valign="middle" >14</td><td align="center" valign="middle" >161.8</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >1.911</td><td align="center" valign="middle" >2.261</td><td align="center" valign="middle" >0.048</td><td align="center" valign="middle" >1.411</td></tr><tr><td align="center" valign="middle" >15</td><td align="center" valign="middle" >161.8</td><td align="center" valign="middle" >56.15</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >2.173</td><td align="center" valign="middle" >1.955</td><td align="center" valign="middle" >0.054</td><td align="center" valign="middle" >3.475</td></tr><tr><td align="center" valign="middle" >16</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >13.85</td><td align="center" valign="middle" >50</td><td align="center" valign="middle" >2.232</td><td align="center" valign="middle" >2.906</td><td align="center" valign="middle" >0.056</td><td align="center" valign="middle" >2.321</td></tr><tr><td align="center" valign="middle" >17</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >78.2</td><td align="center" valign="middle" >5.122</td><td align="center" valign="middle" >4.531</td><td align="center" valign="middle" >0.128</td><td align="center" valign="middle" >1.997</td></tr><tr><td align="center" valign="middle" >18</td><td align="center" valign="middle" >190</td><td align="center" valign="middle" >35</td><td align="center" valign="middle" >21.8</td><td align="center" valign="middle" >4.424</td><td align="center" valign="middle" >5.172</td><td align="center" valign="middle" >0.111</td><td align="center" valign="middle" >1.461</td></tr></tbody></table></table-wrap><p>X<sub>1</sub> = temperature; X<sub>2</sub> = time; X<sub>3</sub> = ratio ethanol-water; Y<sub>exp</sub> = experimental response; Y<sub>cal</sub> = calculated response RS = Reducing Sugars.</p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Validation indicators for mathematical models expressing pentoses, reducing sugars derived from the pre-treatment by the organosolv process</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Validation Indicator</th><th align="center" valign="middle" >Y<sub>C</sub></th><th align="center" valign="middle" >Y<sub>RS</sub></th></tr></thead><tr><td align="center" valign="middle" >R<sup>2</sup></td><td align="center" valign="middle" >85.66%</td><td align="center" valign="middle" >93.74%</td></tr><tr><td align="center" valign="middle" >Adjusted R<sup>2</sup></td><td align="center" valign="middle" >74.56%</td><td align="center" valign="middle" >94.23%</td></tr><tr><td align="center" valign="middle" >AADM</td><td align="center" valign="middle" >0.11</td><td align="center" valign="middle" >0,09</td></tr><tr><td align="center" valign="middle" >Bias factor</td><td align="center" valign="middle" >1.00</td><td align="center" valign="middle" >1</td></tr><tr><td align="center" valign="middle" >Accuracy factor</td><td align="center" valign="middle" >1.03</td><td align="center" valign="middle" >1.05</td></tr></tbody></table></table-wrap><p>a very important factor during cellulose extraction. This factor is followed by the preprocessing time and the E-W ratio, respectively. <xref ref-type="fig" rid="fig8">Figure 8</xref>(c) better illustrates this information where we clearly observe that it is easier to extract celluloses by neglecting the E-W factor ratio, rather than by neglecting the pretreatment time.</p><p>Reducing sugars</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref> shows the response of the RS extraction as a function of factors taken two by two.</p><p>The increase in temperature contributes favorably and significantly to more than 95% during the extraction of reducing sugars. On the other hand, the time will contribute to reducing the yield beyond 35 minutes of treatment. The analysis of <xref ref-type="fig" rid="fig9">Figure 9</xref> illustrates that the E-W ratio will promote optimal extraction of RS in association with temperature (<xref ref-type="fig" rid="fig9">Figure 9</xref>(a)) and time (<xref ref-type="fig" rid="fig9">Figure 9</xref>(b)). Considering all these factors, temperature is the factor that has a great influence on extraction yield.</p></sec><sec id="s3_3"><title>3.3. Optimal Points for the Extraction of Cellulose and Reducing Sugars</title><p>At a temperature of 218.2˚C, a treatment time of 24.65 minutes and a concentration</p><p>of 21.8% ethanol, a yield of 6.70 mg/ml of cellulose is obtained using Minitab software Version 2017. The experimental verification in these conditions allowed us to obtain a value of 7.09 &#177; 0.8 mg/mL, which is not significantly different at the 95% threshold from the expected value (6.70 mg/mL) predicted by the simulator. As for the reducing sugars, the following conditions: 206.24˚C, 56.15 min and 33.19% ethanol, allowed us to obtain a concentration of 3.66 mg/mL. The experimental test gave 4.09 &#177; 0.8 mg/mL, which is very close to the expected theoretical optimal value (3.66 mg/mL).</p><p>Effect of factors on the extraction</p><p>The mathematical answers show that the extraction of Cellulose and reducing sugars illustrates that temperature is the most important factor during this pretreatment operation, which is in agreement with the work of Nanssou et al. [<xref ref-type="bibr" rid="scirp.131339-ref6">6</xref>] ; Kanchanalai et al. [<xref ref-type="bibr" rid="scirp.131339-ref33">33</xref>] and Tchuidjang [<xref ref-type="bibr" rid="scirp.131339-ref7">7</xref>] . Indeed, temperature is the only one of the three factors whose increase contributes favorably to the extraction of RS and cellulose. The model equation reveals that this positive contribution is significant on the RS and not significant on the cellulose. This could be explained by the fact that an increase in temperature would weaken the structure of the fibber contained in the mixture and then the lignocellulosic structure of the substrate, which would cause an increase in the kinetic energy of the water molecules which are the constituents. Under these conditions, the rise in temperature would favour diffusion phenomena and self-ionization releases the hydroxide ions. Time here is the second factor that affects extraction. After a certain value, the amount of cellulose and reducing sugar decreases in the pretreatment reactor. This decrease is due to the destruction of molecules due to an increase in temperature because these compounds are heat-sensitive. They will therefore be denatured when they are exposed for longer than expected to a temperature higher than the temperature necessary for the hydrolysis of their respective substrates.</p><p>The analysis of the two models studied shows that the second model using the experimental design method provides more precision on the results and the applicability of the theory compared to the second model. The continuation of this work will involve implementing this second model in order to size the bioethanol production units while having properties on the raw materials. It emerges from this work that knowledge of the characteristics of precursors is a good indicator for mastering simulation operations because the kinetics of sugar bio decomposition depends on the states and their bioavailability in the basic raw materials.</p></sec></sec><sec id="s4"><title>4. Conclusion</title><p>The objective of this work was to optimize the pretreatment of the mixture of pineapple fibers and cassava peels in a 50:50 ratio by the organosolv process in order to produce bioethanol. Two models were used. The SAEMAN model, which suggests two consecutive first-order reactions, with a dependence of rate constants on temperature; the kinetic parameters of cellulose and glucose hydrolysis are respectively A<sub>1</sub> = 0.05089 min<sup>−1</sup>; A<sub>2</sub> = 1383.03691 min<sup>−1</sup>; E<sub>1</sub> = 5358.1461 J∙mol<sup>−1</sup>; E<sub>2</sub> = 51577.6100 J∙mol<sup>−1</sup>. The other model uses an expert design MINITAB version 2017 by Methodology Response Surfaces was also used in order to compare the results of the numerical simulation and the experimental results. It appears from both models that the optimal pretreatment time and the maximum glucose concentration decrease with temperature, the optimal temperature profile decreases, which attests to the fact that the activation energy of the cellulose hydrolysis reaction is lower than that of glucose. The formation of the latter is therefore favored at low temperatures and limited at high temperatures. This work shows that organosolv pretreatment is an effective method for the delignification of the lignocellulosic structure of precursors. The optimization method developed will make it possible to properly control production performance. To improve this work, it would be wise to carry out a technical-economic analysis of the process and make an energy assessment of the installation.</p></sec><sec id="s5"><title>Acknowledgments</title><p>The authors of this work address their gratitude to Prof. George Elambo Nkeng, Director of National Advanced School of Public Work, of Yaound&#233;, Cameroon for the financial support during this work. Our thanks also go to Tchuidjang Tchamba Thierry, doctoral thesis student in Department of process engineering for his participation in the review of this article.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The authors declare no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s7"><title>Cite this paper</title><p>Djonga, P.N.D., Nkeng, G.E., Sirma, M.H., Tom, A. and Tchuidjang, T.T. (2024) Optimization of the Pretreatment of the Mixture of Cassava Peelings and Pineapple Fibers Using Response Surface Methodology and a Process Simulator for the Bioethanol Production. Energy and Power Engineering, 16, 79-96. https://doi.org/10.4236/epe.2024.162004</p></sec></body><back><ref-list><title>References</title><ref id="scirp.131339-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Kamkuimo, P., Ayouba, S. and Mbarga, A. (2018) étude sur la situation de référence de valorisation des rebuts de l’exploitation forestière et de scierie dans la région de l’Est-Cameroun.</mixed-citation></ref><ref id="scirp.131339-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Sharma, P., Kumar, U.G. and Amiya, K.R. (2016) Dynamic Model for Anaerobic Digestion of Sewage Sludge Considering Haldane Kinetics and Wall Growth Factor. American Journal of Environmental Engineering, 6, 129-139.</mixed-citation></ref><ref id="scirp.131339-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Samomssa, I., Nono, Y.J., Tsamo, C., Dinica, M.R. and Kamga, R. (2019) Influence of Physico-Chemical Parameters on Fuel Briquettes Properties Formulated Withmixture of Biomasses. Journal of Environmental Science and Pollution Research, 5, 338-341. https://doi.org/10.30799/jespr.165.19050202</mixed-citation></ref><ref id="scirp.131339-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">WEO (World Energy Outlook) (2019) Rapport de l’Agence Internationale de l’énergie (AIE) 2019. https://www.iea.org/reports/world-energy-outlook-2019</mixed-citation></ref><ref id="scirp.131339-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Nanssou, P.A.K., Nono, Y.J. and Kapseu, C. (2016) Pretreatment of Cassava Stems and Peelings by Thermohydrolysis to Enhance Hydrolysis Yield of Cellulose in Bioethanol Production Process. Renewable Energy, 97, 252-265. https://doi.org/10.1016/j.renene.2016.05.050</mixed-citation></ref><ref id="scirp.131339-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Tchuidjang, T.T., Noubissié, E. and Ali, A. (2021) Optimization of the Pre-Treatment of White Sawdust (Triplochiton scleroxylon) by the Organosolv Process for the Production of Bioethanol. Oil &amp; Gas Science and Technology, 76, Article No. 23. https://doi.org/10.2516/ogst/2021004</mixed-citation></ref><ref id="scirp.131339-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Stevens, C., et al. (2013) Aqueous Pretreatment for Plant Biomass for Biological and Chemical Conversions to Fuels and Chemicals. John Wiley and Sons, New York.</mixed-citation></ref><ref id="scirp.131339-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Wertz, J.L. and Bédué, O. (2013) Lignocellulosic Biorefineries. EPLF Press, New York. https://doi.org/10.1201/b15443</mixed-citation></ref><ref id="scirp.131339-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Mussato, S.I. (2016) Biomass Fractionation Technologies for Lignocellulosic Feedstock Based Biorefinery. Elsevier, Amsterdam.</mixed-citation></ref><ref id="scirp.131339-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Waldron, K. (2010) Bioalcohol Production: Biochemical Conversion of Lignocellulosic Biomass. Woodhead Publishing, Cambridge.</mixed-citation></ref><ref id="scirp.131339-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Moulijn, J., Makkee, M. and Van Diepen, A. (2013) Chemical Process Technology. Wiley, New York.</mixed-citation></ref><ref id="scirp.131339-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Adabe, K.E. and Ngo-Samnick, E.L. (2014) Cocoa Production and Processing. Agricultural and Food Sciences, Douala.</mixed-citation></ref><ref id="scirp.131339-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Aregheore, E.M. (2002) Chemical Evaluation and Digestibility of Cocoa (Theobroma cacao) by-Products Fed to Goats. Tropical Animal Health and Production, 34, 339-348. https://doi.org/10.1023/A:1015638903740</mixed-citation></ref><ref id="scirp.131339-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Cateto, C.A., Filomena Barreiro, M., Egídio Rodrigues, A. and Naceur Belgacem, M. (2011) Kinetic Study of the Formation of Lignin-Based Polyurethanes in Bulk. Reactive and Functional Polymers, 71, 863-869. https://doi.org/10.1016/j.reactfunctpolym.2011.05.007</mixed-citation></ref><ref id="scirp.131339-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Dien, B.S., Jung, H.J.G., Vogel, K.P., Casler, M.D., Lamb, J.F. and Iten, L. (2006) Chemical Composition and Response to Dilute-Acid Pretreatment and Enzymatic Saccharification of Alfalfa, Reed Canarygrass, and Switchgrass. Biomass and Bioenergy, 30, 880-891. https://doi.org/10.1016/j.biombioe.2006.02.004</mixed-citation></ref><ref id="scirp.131339-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Kabir, M.M., Rajendran, K., Taherzade, M.J. and Horváth, I.S. (2014) Experimental and Economical Evaluation of Bioconversion of Forest Residues to Biogas Using Organosolv Pre-treatment. Bioresource Technology, 178, 201-208. https://doi.org/10.1016/j.biortech.2014.07.064</mixed-citation></ref><ref id="scirp.131339-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Blaise, K.Y. (2013) Optimization of Sugars in Cocoa Shells (Théobroma cacao) in Order to Produce Bioethanol. Master’s Thesis, University of Ngaoundéré, Ngaoundéré.</mixed-citation></ref><ref id="scirp.131339-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Kim, D. and Pan, X. (2010) Preliminary Study on Converting Hybrid Poplar to High-Value Chemicals and Lignin Using Organosolv Ethanol Process. Industrial &amp; Engineering Chemistry Research, 49, 12156-12163. https://doi.org/10.1021/ie101671r</mixed-citation></ref><ref id="scirp.131339-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Kupiainen, L., Ahola, J. and Tanskanen, J. (2014) Kinetics of Formic Acid-Catalyzed Cellulose Hydrolysis. BioRes, 9, 2645-2658.</mixed-citation></ref><ref id="scirp.131339-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Mosier, N., Wyman, C., et al. (2004) Features of Promising Technologies for Pretreatment of Lignocellulosic Biomass. Elsevier, Amsterdam.</mixed-citation></ref><ref id="scirp.131339-ref21"><label>21</label><mixed-citation publication-type="book" xlink:type="simple">Carrasco, R.C. and Oncina, J. (1994) Learning Stochastic Regular Grammars by Means of a State Merging Method. In: Carrasco, R.C. and Oncina, J., Eds., ICGI 1994: Grammatical Inference and Applications, Springer, Berlin, 139-152. https://doi.org/10.1007/3-540-58473-0_144</mixed-citation></ref><ref id="scirp.131339-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Kabel, M.A., Bos, G., Zeevalking, J., Voragen, A.G. and Schols, H.A. (2007) Effect of Pretreatment Severity on Xylan Solubility and Enzymatic Breakdown of the Remaining Cellulose from Wheat Straw. Bioresource Technology, 98, 2034-2042. https://doi.org/10.1016/j.biortech.2006.08.006</mixed-citation></ref><ref id="scirp.131339-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Ashraf Hussain, M., Emdadul Huq, M., Rahman, S.M. and Ahmed, Z. (2002) Estimation of Lignin in Jute by Titration Method. Pakistan Journal of Biological Sciences, 5, 521-522. https://doi.org/10.3923/pjbs.2002.521.522</mixed-citation></ref><ref id="scirp.131339-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">N’Diaye, S., Rigal, L., Larocque, P. and Vidal, P.F. (1996) Extraction of Hemicelluloses from Poplar, Populus tremuloides, Using an Extruder-Type Twin-Screw Reactor: A Feasibility Study. Bioresource Technology, 57, 61-67. https://doi.org/10.1016/0960-8524(96)00041-7</mixed-citation></ref><ref id="scirp.131339-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Joglekar, A.M. and May, A.T. (1987) Product Excellence through Design of Experiments. Cereal Foods World, 32, 857-868.</mixed-citation></ref><ref id="scirp.131339-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Bas, D. and Boyac, I.H. (2007) Modeling and Optimization I. Journal of Food Engineering, 78, 836-845. https://doi.org/10.1016/j.jfoodeng.2005.11.024</mixed-citation></ref><ref id="scirp.131339-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Godin, B., Ghysel, F., Agneessens, R., Schmit, T., Gofflot, S., Lamaudière, S. and Delcarte, J. (2010) Détermination de la cellulose, des hémicelluloses, de la lignine et des cendres dans diverses cultures lignocellulosiques dédiées à la production de bioéthanol de deuxième generation [Cellulose, Hemicelluloses, Lignin, and Ash Contents in Various Lignocellulosic Crops for Second Generation Bioethanol Production]. Biotechnologie, Agronomie, Société et Environnement, 14, 549. https://popups.uliege.be/1780-4507/index.php?id=6186</mixed-citation></ref><ref id="scirp.131339-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Tijskens, L.M., Hertog, M.L. and Nicolai, B.M. (2001) Food Process Modelling. CRC Press, Boca Raton. https://doi.org/10.1201/9781439823064</mixed-citation></ref><ref id="scirp.131339-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Park, N., Kim, Y.H., Koo, B.W., Yeo, H and Choi, I.G. (2010) Organosolv Pretreatment with Various Catalysts for Enhancing Enzymatic Hydrolysis of Pitch Pine (Pinus rigida). Bioresource Technology, 101, 7046-7053. https://doi.org/10.1016/j.biortech.2010.04.020</mixed-citation></ref><ref id="scirp.131339-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Hashimoto, S., Shogren, M.D. and Pomeranz, Y. (1987) Cereal Pentosans: Their Estimation and Significance. I. Pentosans in Wheat and Milled Wheat Products. Cereal Chemistry, 64, 30-34.</mixed-citation></ref><ref id="scirp.131339-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Akpinar, O., Serdal, S., Okan, L. and Abdulvahit, S. (2012) Evaluation of Antioxidant Activity of Dilute Acid Hydrolysate of Wheat Straw during Xylose Production. Industrial Crops and Products, 40, 39-44. https://doi.org/10.1016/j.indcrop.2012.02.035</mixed-citation></ref><ref id="scirp.131339-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Kanchanalai, P., Temani, G., Kawajiri, Y. and Realff, M.J. (2016) Reaction Kinetics of Concentrated-Acid Hydrolysis for Cellulose and Hemicellulose and Effect of Crystallinity. BioRe-sources, 11, 1672-1689. https://doi.org/10.15376/biores.11.1.1672-1689https://ojs.cnr.ncsu.edu/index.php/BioRes/article/view/BioRes</mixed-citation></ref><ref id="scirp.131339-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Mesa, L., González, E., Cara, C., González, M., Castro, E. and Mussatto, S.I. (2011) The Effect of Organosolv Pretreatment Variables on Enzymatic Hydrolysis of Sugarcane Bagasse. Chemical Engineering Journal, 168, 1157-1162. https://doi.org/10.1016/j.cej.2011.02.003</mixed-citation></ref><ref id="scirp.131339-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Fisher, E.H. and Stein, E.A. (1961) DNS Colorimetric Determination of Available Carbohydrates in Foods. Bio-chemical Preparation, 8, 30-37.</mixed-citation></ref></ref-list></back></article>