No description available.

There are two parts in this study. The first part is to test the validity of the assumption of thermodynamic equilibrium between fine particulate (PM2.5) nitrate and ammonium and gas-phase nitric acid (HNO3(g)) and ammonia (NH3(g)). A rough estimation of the characteristic time to achieve thermodynamic equilibrium is first carried out, which suggests that PM2.5 and gas-phase species are in thermodynamic equilibrium. Then equilibrium is tested by calculating the equilibrium concentrations of HNO3(g) and NH3(g) implied by the PM2.5 inorganic composition, temperature and relative humidity observed at the Atlanta Supersite 1999 using ISORROPIA model. The second part of this study is to analyze the ground-level ozone pollution precursor relationships in Hong Kong area. Characteristics of O3 precursors are explored. Trace gases NO and CO, VOCs, absorption coefficient, temperature and solar radiation are associated with the O3 formation. Specific VOC and VOC-sources that contribute most to the formation of photochemical smog are identified. The accuracy of pollutant emission inventories for Hong Kong and PRD region is also assessed. Combined with back trajectory information, dCO/dNOy is used to define whether O3 is locally or regionally occurred. An OBM is used to investigate the relative benefits of various emission-control strategies. Generally the formation of O3 throughout much of Hong Kong area is limited by VOC, in which reactive aromatics are dominant.

Thesis (MEng)--Stellenbosch University, 2014. / ENGLISH ABSTRACT: Thermodynamic models have been investigated extensively since Johannes van der Waals first devised a mathematical relation capable of predicting both vapour and liquid phases for a mixture at equilibrium. With the advent of modern computing power, these equations of state have gone from their humble empirical beginnings to the comprehensive and fundamentally derived models we have today. One such physically sound model is the Statistical Associating Fluid Theory (SAFT) family of equations, derived from the molecular perturbation theories of the 1980’s. The relative youth of this thermodynamic framework has meant that much work has gone into modification and optimisation of the model recently. The variants of particular interest to this work are the simplified perturbed chain SAFT equations with the Jog & Chapman (sPC-SAFTJC) and Gross & Vrabec (sPC-SAFTGV) polar terms. Each of the polar terms supports one adjustable polar parameter that relates to the quantity of polar segments in the reference fluid but not necessarily its position in the carbon chain. The strength of polar interactions is known to decrease as the functional group moves away from the terminal methyl group and the effects of steric hindrance increase. Thus, in question here is whether the models can account for the change in polar interactions associated with the changing position of the polar group, by only adjusting the values of the existing pure component parameters; that is, in lieu of a position specific parameter. The carbonyl group in ketone molecules is one such polar group, and it is this homologous series that is the focus of this study. The decrease in polar interactions as the carbonyl group in a ketone molecule shifts centrally is apparent from the lower boiling points of the isomers where the polar group is central as compared to those where the functional group is nearer the terminal methyl group. The effect of this functional group shift on binary phase behaviour has not previously been assessed for any system however, as the lack of experimental data attests. Thus, experiments had to be conducted to generate phase equilibrium data for systems comprising each structural isomer of a mid-length ketone with a common second component with no functionality. This limitation was imposed to isolate the cause of experimentally observed phenomena to the shifting polar group alone. The generated data could then be appropriately modelled using the polar sPC-SAFT variants and the capabilities of each model, as outlined above, assessed. To this end, isobaric binary vapour-liquid equilibrium data were measured for 2-, 3- & 4-heptanone with three separate normal alkanes of similar length (n-octane, n-nonane & n-decane) at 40kPa. The apparatus used was a dynamic Gillespie VLE still with temperature and pressure accuracies of 0.03°C and 1.6mbar respectively. Equipment verification was achieved through the reproduction of experimental data for the ethanol/1-butanol system at 1.013bar. The vapour and liquid samples for all nine systems were analysed by gas chromatography with a maximum compositional error of ±0.022 mole fraction. All reported data were found to be thermodynamically consistent using both the L/W and McDermott-Ellis consistency tests. When paired with a common n-alkane, all three structural heptanone isomers displayed similar qualitative trends in phase behaviour. Minimum boiling azeotropes were measured in all nine systems; in the high alkane region for n-octane systems (~98 mole%), the equal concentration region for n-nonane systems (34 mole% to 53 mole%) and in the very dilute n-alkane region for n-decane systems (~3 mole%). The n-nonane systems in particular highlighted the effect of shifting functional group, with completely separate phase envelopes away from the pure alkane composition space evident in a particularly small temperature range. Modelling was performed using in-house developed software, with pure component parameters generated for each system using five different regression procedures. The first was traditional fitting of the segment diameter (σ), segment number (m), segment energy (є/k) and the respective polar parameter (xp, np) to DIPPR correlations of pure component saturated vapour pressure, liquid density and the heat of vaporisation. The latter four procedures included the fixing of the polar parameter according to functional group correlations and the three instances of including the binary VLE data set for each of the three alkanes considered in this work. When applied to the nine binary ketone-alkane systems measured in this work, excellent predictions of the experimental data were in evidence in most cases and only small binary interaction parameters were necessary to correlate the data where pure predictions were poor. The performance of the parameter sets based on the fixing of the polar parameter and the inclusion of VLE data were consistent and of a high quality for both models, with near identical parameters generated in all four cases for each of the nine systems. The parameter sets generated in this fashion were shown to be applicable not only to the systems measured in this work, but also successfully predicted the independently measured experimental data of the n-hexane/4-heptanone system. It was thus concluded that either of these regression alternatives are viable for the generation of accurate component parameters, and the choice of VLE data set included is trivial. The pure predictions of the sPC-SAFTGV model were generally better than its sPC-SAFTJC counterpart, particularly in the case of the traditionally regressed parameter sets. sPC-SAFTGV displayed constant qualitative agreement with the experimental data for each of the heptanone isomers with a given n-alkane. The quality of the predictions of sPC-SAFTJC, however, worsened significantly as the polar interactions diminished from 2- to 4-heptanone, with no predictions even possible for the least polar isomer. This was attributed to the different perturbation theories used in the development of these terms, but a more detailed study would be necessary to confirm this. This work thus shows an apparent inability of the sPC-SAFTJC equation of state to account for the decreasing polar interactions associated with the carbonyl group in a ketone molecule shifting centrally, while sPC-SAFTGV produces qualitatively good fits for all three isomers. These flaws can be overcome through the incorporation of VLE data in the regression procedure if such data is available, or otherwise through the use of group specific correlations for fixing the polar parameter value. / AFRIKAANSE OPSOMMING: Sedert Johannes van der Waals die eerste wiskundige verhouding ontwikkel het wat beide die damp- en vloeibare fases van 'n mengsel by ewewig kon voorspel, is die veld van termodinamiese modellering al deeglik ondersoek. Na die koms van die moderne rekenaars het hierdie vergelykings van hul nederige empiriese wortels gegroei tot die omvattende, fundamentele modelle wat ons vandag het. Een so 'n fundamenteel gebaseerde familie van vergelykings is die ‘Statistical Associating Fluid Theory’ (SAFT) modelle, wat afgelei is vanaf molekulêre versteuringsteorieë, ontwikkel in die 1980s. Hierdie relatiewe jong modelle het in die afgelope ruk aansienlike aanpassing en optimering ondervind. Modelvariante van besondere belang tot hierdie werk, is die vereenvoudigde versteurde ketting of ‘simplified perturbed chain’ SAFT vergelykings, met Jog & Chapman (sPC- SAFTJC) en Gross & Vrabec (sPC- SAFTGV) polêre terme. In die sPC-SAFT toestandsvergelyking word elkeen van die polêre terme ondersteun deur een polêre veranderlike. Hierdie veranderlike is afhanklik van die aantal polêre segmente in die verwysingsvloeistof, maar nie noodwendig hul posisie in die koolstofketting nie. Daarteen is dit bekend dat die polêre interaksies tussen molekules swakker word soos die polêre groep wegbeweeg van die terminale metielgroep, en steriese hindernis ʼn groter rol begin speel. Dus is die vraag of die model die verandering in die polêre interaksie, as gevolg van veranderende posisie van die polêre groep, kan voorspel deur in plek van ʼn posisie afhanklike parameter, slegs ʼn aanpassing van die polêre waardes van die suiwer komponente te maak. Die karbonielgroep in ketoon molekules is een so 'n polêre groep, en ʼn homoloë reeks ketone word in hierdie studie ondersoek. Die afname in die polêre interaksie soos wat die karbonielgroep in 'n ketoon molekule weg skuif vanaf die terminale metiel groep is sigbaar deur die afname in kookpunt van die verskillende isomere. Hierdie effek van die funksionele groepsposisie op binêre fasegedrag is nog nie voorheen vir enige stelsels geëvalueer nie en geen eksperimentele data is vrylik beskikbaar nie. Om hierdie tekortkoming in die literatuur aan te spreek, is eksperimentele fase ewewig data gemeet. ʼn Reeks stelsels is ondersoek wat elk bestaan uit ʼn struktuurisomeer van ʼn mid-lengte ketoon en ʼn tweede komponent met geen funksionele bydrae. Eksperimente is so opgestel om die effek van die skuiwende polêre groep op die fasegedrag te isoleer en kwalitatief te ondersoek. Die gegenereerde data is dan gemodelleer met behulp van die polêre sPC- SAFT variante, soos hierbo gespesifiseer, en die vermoëns van elke model is beoordeel. Isobariese binêre fase ewewig data is by 40kPa gemeet vir damp-vloeistof stelsels bestaande uit 2, 3 & 4 heptanoon, gemeng met drie verskillende normaal alkane van vergelykbare kettinglengte (n-oktaan, n-nonaan & n-dekaan). Die apparaat wat gebruik was is 'n dinamiese Gillespie VLE eenheid met temperatuur- en drukakkuraatheid van 0,03°C en 1.6mbar, onderskeidelik. Die akkuraatheid van die toerusting is bevestig deur eksperimentele data vir ʼn etanol/1-butanol stelsel by 1.013bar te reproduseer. Die damp en vloeibare monsters vir al nege stelsels is ontleed deur gaschromatografie met 'n maksimum komposisionele fout van ± 0,022 (molfraksie). Alle data is as termodinamies konsekwent gevind deur van beide die L/W en McDermott-Ellis konsekwentheidstoetse gebruik te maak. Mengsels van die drie strukturele isomere van heptanoon met ʼn gemene n-alkaan het tydens eksperimente soortgelyke kwalitatiewe tendense in fasegedrag getoon. Gedurende eksperimente is die lae kookpunt asiotrope gemeet vir al nege stelsels. Die asiotrope verskyn in die hoë alkaan konsentrasies (~98 mol%) vir n-oktaan stelsels, medium konsentrasies (34 mol% tot 53 mol%) vir n-nonaan stelsels en baie verdunde konsentrasies (~ 3 mol%) vir n-dekaan stelsels. Die n-nonaan stelsels beeld veral die effek van die verskuiwing van die funksionele groep uit, met diskrete fasegrense wat duidelik apart staan van die suiwer alkaan ruimte, binne ʼn klein temperatuurverskil. Modellering van die stelsels is uitgevoer met behulp van sagteware wat in-huis ontwikkel is. Suiwer komponent data is gegenereer vir elke stelsel deur van vyf verskillende regressie prosedures gebruik te maak. Die eerste is die tradisionele passing van die segment deursnee (σ), segment nommer (m), segment energie (є/k) en die onderskeie polêre parameters (xp, np) op DIPPR korrelasies van die suiwer komponent versadigde dampdruk, vloeistof digtheid en die hitte van verdamping. Die oorblywende vier prosedures sluit in die bepaling van die polêre parameter deur funksionelegroep korrelasies, en drie gevalle waar die binêre VLE data vir elk van die drie alkane ingesluit is. Deur hierdie prosedures op die modellering van die nege binêre ketoon/alkaan stelsels toe te pas, is uitstekende passings van die eksperimentele data verkry met slegs baie klein binêre interaksie parameters nodig waar voorspellings minder akkuraat was. Die prestasie van die parameter stelle, gebaseer op die bepaling van die polêre parameter en die insluiting van VLE data, is konsekwent en van 'n hoë gehalte vir albei modelle, met 'n byna identiese parameters gegenereer in al vier gevalle vir elk van die nege stelsels. Die parameter stelle wat op hierdie metode gegenereer is, is nie net toepaslik gevind op eksperimentele data gemeet in hierdie werk nie, maar ook op onafhanklike data vir die n-heksaan/4-heptanoon stelsel. Daar is tot die gevolgtrekking gekom dat beide van die regressie alternatiewe lewensvatbaar is vir die generasie van akkurate suiwer komponent parameters, en dat die insluiting van die VLE data triviaal is. Die suiwer sPC - SAFTGV voorspelling was oor die algemeen beter as die suiwer sPC- SAFTJC model met die voorspelling van data, veral in die geval van passings gedoen met parameters verkry vanaf tradisionele regressie metodes. sPC- SAFTGV het ʼn voortdurende, kwalitatiewe ooreenkoms met eksperimentele data getoon vir elk van die nege stelsels. Daarteen het voorspellings deur sPC- SAFTJC beduidend verswak soos die polêre interaksies afgeneem het vanaf 2- na 4- heptanoon, met geen akkurate voorspelling moontlik vir die minste polêre isomeer nie. Die verskynsel kan toegeskryf word aan die verskil in versteuringsteorieë wat gebruik word in die ontwikkeling van die onderskeie toestandsvergelykings, maar ʼn meer in-diepte ondersoek is nodig om hierdie teorie te bevestig. Hierdie werk toon dus 'n skynbare onvermoë van die sPC - SAFTJC toestandsvergelyking om die verandering in polêre interaksie, as gevolg van die veranderende posisie van die polêre groep, vir die karbonielgroep in ʼn ketoon te voorspel, terwyl die sPC-SAFTGV toestandsvergelyking goeie kwalitatiewe passings vir al drie isomere bied. Hierdie tekortkominge kan oorkom word deur VLE data, indien beskikbaar, in die regressie prosedure in te sluit, of deur die gebruik van groep spesifieke korrelasies vir die aanpassing van die polêre parameter.

La thèse porte sur la modélisation de la transformation de l'austénite en ferrite dans les aciers en mettant l'accent sur les conditions thermodynamiques et cinétiques aux interfaces alpha/gamma en cours de croissance de la ferrite. Dans une première partie, la thèse se concentre sur la description des équilibres thermodynamiques entre alpha et gamma à l'aide de la méthode CalPhad. Nous avons développé un nouvel algorithme hybride combinant la construction d'une enveloppe convexe avec la méthode classique de Newton-Raphson. Nous montrons ses possibilités pour des aciers ternaire Fe-C-Cr et quaternaire Fe-C-Cr-Mo dans des cas particulièrement difficiles. Dans un second chapitre, un modèle à interface épaisse a été développé. Il permet de prédire l'ensemble du spectre des conditions à l'interface alpha/gamma au cours de la croissance de la ferrite, de l'équilibre complet au paraéquilibre avec des cas intermédiaires des plus intéressants. Nous montrons que de nombreux régimes cinétiques particuliers dans les systèmes Fe-C-X peuvent être prévus avec un minimum de paramètres d'ajustement, principalement le rapport entre les diffusivités de l'élément substitutionnel dans l'interface épaisse et dans le volume d'austénite. Le troisième chapitre porte sur l'étude d'un modèle de champ de phase. Une analyse approfondie des conditions à l'interface données par le modèle est réalisée en utilisant la technique des développements asymptotiques. En utilisant les connaissances fournies par cette analyse, le rôle de la mobilité intrinsèque d'interface sur la cinétique et les régimes de croissance est étudié, à la fois dans le cas simple d'alliages binaires Fe-C et dans le cas plus complexe d'alliages Fe-C-Mn / Transformation in steels focusing on the thermodynamic and kinetics conditions at the alpha/gamma interfaces during the ferrite growth. The first chapter deals with the determination of thermodynamic equilibria between alpha and gamma with CalPhad thermodynamic description. We have developed a new hybrid algorithm combining the construction of a convex hull to the more classical Newton-Raphson method to compute two phase equilibria in multicomponent alloys with two sublattices. Its capabilities are demonstrated on ternary Fe-C-Cr and quaternary Fe-C-Cr-Mo steels. In the second chapter, we present a thick interface model aiming to predict the whole spectrum of conditions at an alpha/gamma interface during ferrite growth, from full equilibrium to paraequilibrium with intermediate cases as the most interesting feature. The model, despite its numerous simplifying assumptions to facilitate its numerical implementation, allows to predict some peculiar kinetics in Fe-C-X systems with a minimum of fitting parameters, mainly the ratio between the diffusivities of the substitutional element inside the thick interface and in bulk austenite. The third chapter deals with the phase field model of austenite to ferrite transformation in steels. A thorough analysis on the conditions at the interface has been performed using the technique of matched asymptotic expansions. Special attention is given to clarify the role of the interface mobility on the growth regimes both in simple Fe-C alloys and in more complex Fe-C-Mn alloys

Removal of polynuclear aromatics from diesel fuel has become a focus of intense research due to the stringent environmental legislation associated with clean fuels. In this work, selective ring opening of model compound decalin over the set of catalysts comprising of 1) Ir-Pt supported on mesoporous Zr-MCM-41, large and medium pore zeolites like HY and H-Beta and 2) Ni-Mo/carbide on HY, H-Beta, Al-SBA-15, ¥ã- alumina and silica alumina were studied. All the catalysts were extensively characterized by BET surface area measurement, CO-chemisorption, XRD, FTIR, TPR and TPD of ammonia. Ring opening of decalin was studied on these catalysts in a trickle-bed reactor in a temperature range of 200- 400 ¡ÆC, pressure range of 2-7 MPa and LHSV of 1 to 3 h- 1. 31.7 and 65.0 wt.% of RO yield and selectivity were observed on Ir-Pt/HY catalyst at 220 ¡ÆC, whereas 34.0 and 40.0 wt.% of ring opening yield and selectivity were observed on Ni-Mo carbide/HY catalyst at 240 ¡ÆC. From the model compound studies, Ir-Pt/HY, Ni-Mo carbide/HY and Ni-Mo carbide/H-Beta were selected for study of hydrotreated light gas oil in a trickle bed reactor. Ni-Mo carbide/HY performed better over other catalysts and increased the cetane index of hydrotreated light gas oil by 12 units at 325 ¡ÆC. A first order kinetic model was fitted for the hydrotreated light gas oil study. 89, 111 and 42 KJ/gmol of activation energies was observed for dearomatization, aromatization and naphthenes cracking steps, respectively. The thermodynamic equilibrium calculations reveal that the selectivity of ring opening products of decalin can be maximized by favoring the formation of unsaturated compounds at higher operating temperatures. Energetics of dealkylation and ring opening reactions of naphthenes in gas phase and on the surface of Br©ªnsted acid sites were calculated using quantum chemical simulations. In iv gas phase, ratio of Arrhenius activation energies for forward and reverse reactions of RO and dealkylation reactions are 1.92 and 1.82 respectively. Deakylation on different level clusters revealed that surface reaction is the rate controlling.

Removal of polynuclear aromatics from diesel fuel has become a focus of intense research due to the stringent environmental legislation associated with clean fuels. In this work, selective ring opening of model compound decalin over the set of catalysts comprising of 1) Ir-Pt supported on mesoporous Zr-MCM-41, large and medium pore zeolites like HY and H-Beta and 2) Ni-Mo/carbide on HY, H-Beta, Al-SBA-15, ¥ã- alumina and silica alumina were studied. All the catalysts were extensively characterized by BET surface area measurement, CO-chemisorption, XRD, FTIR, TPR and TPD of ammonia. Ring opening of decalin was studied on these catalysts in a trickle-bed reactor in a temperature range of 200- 400 ¡ÆC, pressure range of 2-7 MPa and LHSV of 1 to 3 h- 1. 31.7 and 65.0 wt.% of RO yield and selectivity were observed on Ir-Pt/HY catalyst at 220 ¡ÆC, whereas 34.0 and 40.0 wt.% of ring opening yield and selectivity were observed on Ni-Mo carbide/HY catalyst at 240 ¡ÆC. From the model compound studies, Ir-Pt/HY, Ni-Mo carbide/HY and Ni-Mo carbide/H-Beta were selected for study of hydrotreated light gas oil in a trickle bed reactor. Ni-Mo carbide/HY performed better over other catalysts and increased the cetane index of hydrotreated light gas oil by 12 units at 325 ¡ÆC. A first order kinetic model was fitted for the hydrotreated light gas oil study. 89, 111 and 42 KJ/gmol of activation energies was observed for dearomatization, aromatization and naphthenes cracking steps, respectively. The thermodynamic equilibrium calculations reveal that the selectivity of ring opening products of decalin can be maximized by favoring the formation of unsaturated compounds at higher operating temperatures. Energetics of dealkylation and ring opening reactions of naphthenes in gas phase and on the surface of Br©ªnsted acid sites were calculated using quantum chemical simulations. In iv gas phase, ratio of Arrhenius activation energies for forward and reverse reactions of RO and dealkylation reactions are 1.92 and 1.82 respectively. Deakylation on different level clusters revealed that surface reaction is the rate controlling.

The concern about fossil fuel reserves depletion, pollution due to extensive use of conventional fuels and stricter standards for emission are the forces that drive research on alternative fuels. Amongst these are the biodiesel in which the use of vegetable oils has become an attractive to replace diesel from petroleum because of its physico-chemical characteristics and no toxic emissions. This work aimed to study the effect of diferente variables on conversion of peanut oil to biodiesel; determine the influence of temperature and pH in the purification process in order to optimize the amount of water used and study the thermodynamic equilibrium of the system biodiesel-glycerine-alcohol (decreasing the production of effluents). The process of production and purification of peanut biodiesel obtained by transesterification was studied. In this reaction was used sodium hydroxide (NaOH) as a basic catalyst, anhydrous ethanol and peanut crude oil (non-degummed) in a pilot production unit. In the purification process was used sulphuric acid diluted in distilled water and drying agent (magnesium sulphate). A 22 factorial design was performed to determine the influence of chemical variables on the amount of water required for thepurification. The variables for the production process were: temperature (30o, 70oC), stirring speed (150 rpm, 350 rpm), molar ratio oil/alcohol (1:4 and 1:10), amount of catalyst (0,5% and 1,5%) and reaction time (0,5 h, 2 h). For the purification process the variables analysed were: temperature (30o and 70oC) and pH (2 and 5). For the study of the system´s thermodynamic equilibrium was used purified biodiesel, jacketed glass reactor with circulating water, titration with biodiesel and glycerine. The evaluation parameter for the production process is the conversion to biodiesel determined by gaseous chromatography and for the purification process were the amount of water utilized and the determination of the influence of the temperature on the system´s thermodynamic equilibrium. The conversion to biodiesel varied between total saponification to a conversion of 87%. In the purification was observed that the effect of the pH variable could be ignored and the effect of the temperature and the interaction of the temperature and pH can be considered. The thermodynamic equilibrium of the system biodiesel/glycerine/alcohol and the tie lines did not have its behaviour influenced by the different temperatures, not influencing the miscibility of the mixture. / A preocupação com o esgotamento das reservas de combustíveis fósseis, a poluição devido ao grande uso de combustíveis convencionais e as normas rigorosas em relação às emissões são as forças que impulsionam as pesquisas com combustíveis alternativos. Dentre estes, temos o biodiesel no qual o uso de óleos vegetais tem se tornado uma alternativa atrativa em substituição ao diesel do petróleo devido às suas características físico-químicas e emissões não tóxicas. Esse trabalho teve como objetivos estudar o efeito de diferentes variáveis na conversão do óleo de amendoim a biodiesel; determinar a influência da temperatura e do pH no processo de purificação, a fim de otimizar a quantidade de água utilizada e, estudar o equilíbrio termodinâmico do sistema biodiesel-glicerina-álcool (diminuindo a produção de efluentes). Foi estudado o processo de produção e purificação do biodiesel de amendoim, obtido através da reação de transesterificação. Nessa reação foi utilizado hidróxido de sódio (NaOH) como catalisador básico, etanol anidro e óleo de amendoim bruto (não-degomado) em uma unidade piloto de produção. No processo de purificação foi utilizado ácido sulfúrico diluído em água destilada e agente dessecante (sulfato de magnésio). Um planejamento fatorial 2² foi realizado a fim de se determinar a influência das variáveis químicas na quantidade de água requerida para a purificação. As variáveis para o processo de produção foram: temperatura (30° e 70°C), velocidade de agitação (150 rpm e 350 rpm), razão molar óleo/álcool (1:4 e 1:10), quantidade de catalisador em massa (0,5% e 1,5%) e tempo de reação (0,5 h e 2 h). Para o processo de purificação foram analisadas as variáveis: temperatura (30° e 70°C) e pH (2 e 5). Para o estudo do equilíbrio termodinâmico foi utilizado o biodiesel purificado, reator de vidro com camisa de circulação de água, titulação com glicerina e biodiesel. O parâmetro de avaliação para a produção é o rendimento em biodiesel determinado por cromatografia gasosa e para o processo de purificação, a quantidade de água utilizada e a determinação da influência da temperatura no equilíbrio termodinâmico do sistema. A conversão a biodiesel variou de total saponificação até uma conversão de 87%. Na purificação, foi observado que o efeito da variável pH pode ser desprezado e o efeito da temperatura e da interação dessas duas variáveis pode ser considerado. O equilíbrio termodinâmico do sistema biodiesel/glicerina/álcool e suas tie lines não sofreram influência de comportamento nas diferentes temperaturas estudadas, não influenciando na miscibilidade da mistura.

The substance of present work is to provide standard thermodynamic functions (STF) of small size molecules for the calculation of the composition and thermodynamic properties of low-temperature plasma, and also method for such a calculation applying obtained STF under non-ideal plasma conditions. With a view to further application in modelling the phenomena in thermal plasma the range of pressures is limited to the region from 0.01 bar to 100 bar, and that of temperature to 298.15–50 kK. To obtain STF the method of partition function resulting from statistical mechanics was proposed. State of art in the given scientific area and theoretical basis of the statistical mechanics required for establishing of the proposed method together with discussion of partition function divergence problem have been reviewed. For the calculation of STF of diatomic molecules the method of direct summation has been employed, whereas, as for the larger size molecules, the rigid rotor and harmonic oscillator model have been generally adopted. The spectral data required for the calculations have been taken from literature, or, in selected cases, these have been computed by quantum chemistry ab initio techniques. The resulting STF have been included into already existing database system of thermodynamic properties and those can serve as input data for subsequent thermodynamic calculations. A general method has been worked out for the purpose of the computation of thermodynamic properties and composition of non-ideal homogenous plasma system in thermodynamic equilibrium. The method is based on minimizing total Gibbs energy to compute at constant pressure or Helmholtz energy to compute at constant volume. The computation algorithm was implemented into computer program and subsequently applied to the computation of the composition and thermodynamic properties of SF6 dissociation and ionization products using obtained STF.

Malgré sa disponibilité, la biomasse de troisième génération est peu utilisée pour la production d’énergie en raison notamment de sa forte contamination en métaux lourds. Afin d’améliorer la compréhension du comportement et de l’éventuel rôle catalytique ou inhibiteur de ces métaux au cours de la pyro-gazéification, et optimiser ainsi la valorisation énergétique de la biomasse contaminée, une approche intégrée expérimentale et thermodynamique a été utilisée. Pour mener à bien ces travaux, une méthode d’insertion dans le bois du métal spécifiquement étudié pour son rôle catalytique, le nickel, a été développée. Cette méthode permet de maîtriser la composition et l’homogénéité des échantillons étudiés, sans altérer la structure du bois. Les essais de pyro-gazéification de bois brut et contaminés dans un réacteur à lit fixe, montrent que, même à faibles teneurs en nickel (entre 0,016 et 0,086 mol/ kg de bois), celui-ci a une activité catalytique importante. Les analyses des phases gaz et solides confirment que l’augmentation de la teneur en Ni dans le bois favorise les réactions de pyro-gazéification qui ont ainsi lieu à des températures plus basses (d’environ 100 °C). En parallèle, une étude thermodynamique a été réalisée afin de conforter les résultats expérimentaux et d’aider à mieux comprendre le comportement des métaux lourds présents dans le bois. Les calculs thermodynamiques, basés sur la minimisation de l’enthalpie libre de Gibbs d’un système constitué de 28 éléments (C, H, O, N et 24 éléments mineurs ou traces), fournissent la spéciation et la répartition dans les différentes phases (gaz, liquides et cendres) des métaux lourds, dont le Ni. Ces calculs, et des analyses spécifiques, ont notamment permis de déterminer la spéciation chimique et physique du nickel au cours de la pyro-gazéification et fournir ainsi une base théorique à son activité catalytique en fonction de la température. / Despite its availability, contaminated biomass is not widely used for energy production due to its high contamination with heavy metals. Understanding the role and behavior of those heavy metals in the pyro-gasification process is a major scientific challenge for optimizing the thermochemical valorization of contaminated biomass. For this purpose, experimental and thermodynamic approaches were used. To carry out this work, a methodology for inserting the studied metal (nickel) in the wood matrix was developed without modifying its structure and controlling the composition and homogeneity of the studied samples. The pyro-gasification tests in a fixed bed reactor have shown that nickel has a catalytic performance even at low concentrations (between 0.016 and 0.086 mol / kg wood). The analysis of gas and solid products confirmed that the increase of Ni content in the wood generated a decrease in the samples pyro-gasification temperature by 100°C. In parallel, a thermodynamic study was conducted to confirm the experimental results and improve the understanding of heavy metals behaviors during pyro-gasification reactions. Thermodynamic calculations based on the Gibbs free energy minimization of a system consisting of 28 elements provide speciation and distribution in the different phases of heavy metals, including Ni. These calculations and specific analyses have made it possible to determine the chemical and physical speciation of nickel during pyro-gasification and thus provide a theoretical basis for its catalytic activity as a function of temperature.

[pt] Devido aos efeitos nocivos dos combustíveis de origem fóssil, cresceu o interesse no uso de fontes sustentáveis de energia, em particular os derivados da biomassa. O presente trabalho consiste no desenvolvimento de uma ferramenta de modelagem computacional da composição da mistura gasosa, chamado gás de síntese ou syngas resultante da gaseificação, baseado no equilíbrio termodinâmico. Assim, são avaliadas a flexibilidade, a confiabilidade e a estabilidade computacional das abordagens de equilíbrio estequiométrico e de minimização da energia livre de Gibbs para determinação das composições do syngas formado no processo. São investigados a influência da umidade relativa da biomassa (0, 20 e 40 por cento em base mássica), da temperatura de gaseificação (800, 1100 e 1400 K), da razão ar combustível em relação a estequiometria (0,2, 0,3 e 0,4) e da quantidade de oxigênio no agente de gaseificação (oxigênio puro, ar ou mistura equimolar dos dois) na composição do syngas, bem como nas eficiências energética e exergética, para biomassas compostas de carbono, oxigênio e hidrogênio. O melhor modelo é usado para otimizar a concentração de agente oxidante e de temperatura que maximizem a eficiência energética para biomassas com 20 por cento de umidade relativa no ar. Esses dados, estão usados para desenvolver modelos matemáticos preditivos usando o conceito de análise de misturas, cuja robusteza é avaliada por análise de variância (ANOVA). Esses resultados estão comparados com as simulações realizadas em um modelo teórico considerando a cinética das reações. Observou-se que o método usando a minimização da energia livre de Gibbs apresentou os melhores resultados para uma aproximação inicial, se tornando uma boa ferramenta para a fase de projeto. No entanto, é necessário desenvolvimento de outros modelos como cinético, CFD ou RNA para análises precisas da gaseificação. / [en] Due to the harmful effects of fossil fuels, interest in the use of sustainable energy sources, in particular those derived from biomass, has grown. The present work consists in the development of a computational modeling tool for the composition of the gas mixture, called synthesis gas or syngas resulting from gasification, based on thermodynamic equilibrium. Thus, the flexibility, reliability and computational stability of the stoichiometric equilibrium and Gibbs free energy minimization approaches are evaluated to determine the compositions of the syngas formed in the process. The influence of biomass relative humidity (0, 20 and 40 percent on mass basis), gasification temperature (800, 1100 and 1400 K), fuel air ratio related to stoichiometry (0.2, 0.3, 0.4) and the amount of oxygen in the gasification agent (pure oxygen, air or equimolar mixture of the two) are investigated in the composition of the syngas, as well as in the energy and exergy efficiencies, for biomasses composed of carbon, oxygen and hydrogen. The best model is used to optimize the oxidizing agent concentration and temperature that maximize energy efficiency for biomass with 20 percent relative humidity in the air. These data are used to develop predictive mathematical models using the concept of analysis of mixtures, whose robustness is assessed by analysis of variance (ANOVA). These results are compared with simulations performed in a theoretical model considering the reaction kinetics. the method using Gibbs free energy minimization presented the best results for an initial approximation, becoming a good tool for a design phase. However, it is necessary to develop other models such as kinetic, CFD or RNA for accurate gasification analysis.