Towards a comprehensive model for cooking processes

Abstract

The objective of this study was to develop mathematical models that describe the course of delignification and carbohydrate degradation during cooking processes as a function of cooking temperature and cooking chemicals concentrations. The cooking processes investigated in this work were the kraft and the autohydrolysis processes. In the case of kraft cooking different existing delignification models were fitted to experimental data from both soft- and hardwood cooks. The fact that data with different temperatures were present enabled the estimation of the activation energies for the reactions associated with the models. For the available data it turned out that a model containing two lignin components of differing reactivities was sufficient for a satisfactory description of the observed delignification. Furthermore, two modifications aimed at augmenting the authenticity of the model were conducted. First, Donnan theory, accounting for the time development of the chemical concentrations at the actual location of the delignification reactions, was incorporated to the model. Second, the rate constants in the model were calculated according to the real temperature during heating and cooling instead of resorting to the computation simplifying concept of isothermal time. Moreover, the square root relationship between the reciprocal lignin content and the anthraquinone charge earlier established for soda-anthraquinone cooks was confirmed in the kraft-anthraquinone context. For carbohydrate degradation during kraft cooking a model based on the endwise reactions peeling, stopping and alkaline hydrolysis was developed. The model provided satisfactory fits to the yield data of all the three carbohydrates (cellulose, glucomannan and xylan) under scrutiny. The portions of material degraded by the different reactions present in the model were calculated based on the parameter estimates obtained from the fits. Furthermore, the impact of dissolution was compared to that of the endwise reactions; for glucomannan and cellulose the endwise reactions degrade more material than does the dissolution, whereas the opposite is true for xylan. The yield loss model for the carbohydrates was also combined with a model for the decrease in the degree of polymerization. By performing a simultaneous fit of the two models an enhanced estimate of the model parameters was obtained. For the autohydrolysis cooks reaction schemes for the delignification and the degradation of xylan and glucan were developed and the corresponding kinetic models were fitted to the data resulting in estimates for the kinetic parameters. The kinetic models were combined with the diffusion equation resulting in a simulation program for the effect of cooking on the wood components on chip level. Furthermore, the kinetic models together with the Ergun equation for the forced flow of the heated cooking liquor through the packed bed led to a simulation program of the spatial developments of delignification and carbohydrate degradation on reactor level.

Tämän työn tavoitteena oli kehittää matemaattinen malli kuvaamaan keiton aikana tapahtuvaa puun ligniinin ja hiilihydraattien hajoamista sekä hajoamisnopeuksien riippuvuutta keittolämpötilasta ja keittokemikaalien konsentraatioista. Sulfaattimenetelmä ja kuumavesiuutto olivat tutkimuksen kohteina. Sulfaattikeiton tapauksessa olemassa olevia ligniinin hajoamisen kineettisiä malleja sovitettiin sekä lehtipuu- että havupuukeitoista saatuun kokeelliseen dataan. Koska keittoja oli suoritettu eri lämpötiloilla, oli mahdollista arvioida malliin liittyvien reaktioiden aktivointienergiat. Saatavilla olevan datan tapauksessa malli, jossa oli kaksi eri reaktionopeudet omaavaa ligniinin lajia, osoittautui riittäväksi kuvaamaan tyydyttävästi havaittua ligniinin määrän vähenemistä. Tämän lisäksi kaksi eri tapaa lisätä mallin todenmukaisuutta otettiin käyttöön. Ensinnäkin Donnanin teoria, jonka avulla voidaan laskea kemikaalien konsentraatiot siinä kohdin missä varsinainen reaktio tapahtuu, liitettiin malliin. Toiseksi mallin reaktiokertoimet laskettiin käyttäen kuumennus-ja jäähdytysaikajaksojen todellisia lämpötiloja laskuja yksinkertaistavaan isotermiseen aikaan liittyvän kuvitteellisen vakiolämpötilan asemesta. Lisäksi vahvistettiin sooda-AQ keitosta tutun ligniinimäärän käänteisarvon ja antrakinonimäärän neliöjuuren välisen riippuvuuden pätevyys myös sulfaattikeiton tapauksessa. Hiilihydraattien hajoamiselle sulfaattikeitossa kehitettiin polymeeriketjujen reunayksikköihin kohdistuviin reaktioihin (peeling, stopping, alkaline hydrolysis) perustuva malli. Malli sopi hyvin yhteen kaikkien kolmen puun tärkeimpien hiilihydraattien selluloosan, glukomannaanin ja ksylaanin kokeellisten saantoarvojen kanssa. Mallin eri reaktioiden kautta hajoavan materiaalin osuus laskettiin sovituksista saatujen parametriestimaattien avulla. Tämän lisäksi vertailtiin polymeerien liukenemisen merkitys suhteessa polymeeriketjujen reunayksikköihin kohdistuvien reaktioiden merkitykseen; glukomannaanin ja selluloosan tapauksessa päätyreaktiot hävittävät enemmän ainetta kun liukeneminen, kun taas ksylaanin tapauksessa tilanne on käänteinen. Hiilihydraattien saantohävikin malli yhdistettiin myös polymerisaatioasteen vähenemistä kuvaavaan malliin. Mainittujen kahden mallien samanaikainen sovitus vahvisti kineettisten parametrien estimaatit. Kuumavesiuutolle kehitettiin ligniinin, ksylaanin ja glukaanin hajoamisen reaktiokaaviot ja niihin liittyvät kineettiset mallit sovitettiin dataan tuloksena kineettisten parametrien estimaatit. Kineettiset mallit yhdistettynä diffuusioyhtälön kanssa mahdollistivat keiton simulointiohjelman hakepalasen tasolla. Lisäksi kineettiset mallit yhdessä hakepedin läpi kulkevan keittonesteen virtauksen kuvaavan Ergunin yhtälön kanssa johtivat ligniinin ja hiilihydraattien hajoamisen reaktoritason simulointiohjelmaan.