Para abordagem do controle não-linear geométrico, a síntese do controle e elaborada diretamente a partir da descrição do processo com a dinâmica não-linear em espaço de estados. O presente trabalho trata da aplicação dos principais conceitos e formalismos do controle não-linear geométrico para os processos multivariaveis típicos da engenharia química: o controle não-linear continuo da coluna de destilação e o controle não-linear discreto da unidade de craqueamento catalítico em leito fluidizado. A síntese e o projeto do controlador não-linear são enfocados separadamente. O projeto do controlador tem importância pratica para as aplicações industriais. O presente trabalho apresenta metodologias para a abordagem dos seguintes aspectos da aplicação multivariavel do controle geométrico não-linear: (a) como relaxar a sintonia do controlador interno de desacoplamento não-linear; (b) como definir o controlador externo como controle linear de alocação de pólos com coeficientes de hurwitz; (c) neste controlador externo, como incluir a ação integral com prevenção da saturação; e (d) como definir a dinâmica dossetpoints externos. / For the geometric nonlinear control approach, the controller synthesis is elaborated directly from the nonlinear dynamics state space description of the process. This work concerns the application of the main concepts and formalisms of the geometric nonlinear control theory to typical multivariable (MIMO) chemical engineering process as illustrative case studies: the continuous nonlinear control of the distillation column and the discrete nonlinear control of the fluid catalytic cracking unit. The synthesis and the project issues of the nonlinear controller are focused separately. The controller project has the practical importance for the industrial controller applications. This work applies the methodologies to approach the following issues for the MIMO applications of the geometric nonlinear control: (a) to detune the internal nonlinear decoupling controller; (b) to define the external controller as linear pole-placement controllers with Hurwitz coefficients; (c) to include the integral action with anti-reset windup on this external controllers and (d) to define the dynamics of the external setpoints.

A grande dificuldade na simulação de processos envolvendo substâncias naturais é o reduzido número de propriedades termo físicas e dados experimentais disponível na literatura. Devido à diversidade e complexidade das moléculas presentes, o estudo do processamento de compostos naturais ainda é pequeno. Neste trabalho, são investigados modelos para predição de tais propriedades termodinâmicas, fundamentais à modelagem e simulação do processo de destilação em batelada de óleos essenciais. Um método de contribuição de grupos é utilizado para predição da pressão de vapor de substâncias puras (CSGC-PVR), quando dados experimentais não são disponíveis. Em relação ao cálculo do coeficiente de atividade, um modelo totalmente preditivo baseado em química quântica computacional é utilizado (COSMO-SAC). De modo geral, este trabalho tem como objetivo a aplicação de propriedades termodinâmicas preditas em um modelo dinâmico, capaz de descrever o processo de retificação através da destilação em batelada. As simulações foram realizadas no simulador orientado a equações EMSO, a fim de demonstrar a viabilidade do método. Os óleos essenciais de Eucalyptus globulus Labill, Eucalyptus citriodora Hook e Cymbopogon winterianus Jowitt foram selecionados devido ao potencial econômico de suas frações. Inicialmente, um conjunto de dados experimentais é comparado às simulações para validação do método. Em um segundo momento, as simulações são realizadas para correlacionar as condições operacionais e cortes de destilado para as misturas dadas. Os perfis de composição e a determinação dos cortes majoritários são apresentados. Duas alternativas são avaliadas para obtenção de maior pureza nos cortes; o uso de bateladas de reciclo e o controle da razão de refluxo através de uma ação proporcional. A influência de diferentes configurações de coluna também é avaliada, variando o número de estágios teóricos, regime de refluxo e acúmulo de coluna. Os métodos utilizados se mostram capazes de predizer o comportamento do processo de fracionamento e podem ser aplicados a um grande número de óleos essenciais. / The great obstacle in simulating processes involving natural products is the small number of thermo-physical properties and experimental data available in the literature. Due to the diversity and complexity of these molecules, studies concerning essential oils processing is still limited. In this work, thermodynamic models are investigated in order to predict such properties, which are requisites to the modelling and simulation of the essential oil batch distillation process. In this work a group contribution method is used to predict the vapor pressure (CSGC-PVR) when experimental data was unavailable. Regarding the activity coefficients, a fully predictive model based on quantum calculations (COSMO-SAC) is used. Moreover, this work also aims the use of those predicted properties in a dynamic model, capable of describing the fractionation process by batch distillation. The simulations were performed in the equation-oriented simulator EMSO, in order to demonstrate the feasibility of the method. Oils of Eucalyptus globulus Labill, Eucalyptus citriodora Hook and Cymbopogon winterianus Jowitt were selected due to their high production and economic potential of its fractions. Initially, a set of experimental data is compared to the simulations in order to validate the method. Then, simulations are performed to relate operational conditions and distillate cuts for the given mixtures. The distillate composition profiles and the determination of the major cuts are presented. Two alternatives are evaluated in order to obtain cuts of higher purity; the use of a recycle batch and a reflux ratio proportional control based on the distillate composition. The influence of different column configurations is also evaluated, varying its number of theoretical stages, reflux ratio regime and column holdup. The methods used are shown able to predict a reasonable behavior for the fractionation process of the essential oils tested and the method can be extended to most of the essential oils.

Neste trabalho foram determinadas as razões de reatividade para a copolimerização entre estireno e VeoVa-10 através de estudo experimental de copolimerização em solução em tolueno. Os parâmetros foram obtidos através do ajuste aos dados experimentais dos modelos de Mayo-Lewis e Meyer-Lowry, utilizando o critério dos erros-nas-variáveis (método dos mínimos quadrados não-linear que considera a existência de erros em todas as variáveis medidas). A metodologia experimental desenvolvida, de polimerização em solução em ampolas preparadas dentro de uma glove box e submetidas a temperatura de 80°C, resultou em baixa conversão dos monômeros, adequada para a estimação das razões de reatividade pelos modelos. As amostras foram analisadas quantitativamente por espectroscopia de infravermelho para determinar, ao longo da reação, a conversão e as frações molares de cada um dos monômeros na mistura, sendo a composição do copolímero obtida pelo consumo dos monômeros. Foram obtidos valores de razões de reatividade bastante diferentes (restireno = 49.926 e rVeoVa10 = 1,062), evidenciando que a copolimerização destes monômeros em batelada leva a grandes desvios de composição, com tendência à formação de dois homopolímeros. O mesmo sistema de copolimerização foi também estudado no processo de polimerização em emulsão, em processos batelada alimentada e batelada intermitente. A conversão global dos monômeros e a distribuição de diâmetro médio das partículas foram analisadas durante a reação através das amostras coletadas ao longo do processo. Os efeitos das estratégias de alimentação avaliados foram apenas preliminares, pois não apresentaram resultados conclusivos. / In this work the reactivity ratios for the copolymerization between styrene and VeoVa-10 were determined through an experimental study of solution copolymerization in toluen. The parameters were obtained by fitting the experimental data with the Mayo-Lewis and Meyer-Lowry models using the error-in-variables criterion (non-linear least squares method that considers the existence of errors in all measured variables). The experimental methodology developed, solution polymerization in ampoules prepared in a glove box and reacted at 80 °C, resulted in low conversion of the monomers, suitable for estimation of the reactivity ratios by the models employed. The samples were analyzed quantitatively by infrared spectroscopy to determine, during the reaction, the conversion and the mole fractions of each monomer in the mixture, the composition of the copolymer being obtained by the consumption of the monomers. Very different reactivity ratios were obtained (rstyrene = 49.926 and rVeoVa10 = 1.062), showing that the batch copolymerization of these monomers leads to large composition drift, with a tendency to form two homopolymers. The same copolymerization system was also tested in emulsion polymerization process, in semi-continuous and intermittent batch modes. The overall conversion of the monomers and the mean particle diameter distribution were analyzed during the reaction through the samples collected throughout the process. The effects of the feeding strategies evaluated were only preliminary, as they did not present conclusive results.

O propósito deste trabalho foi o desenvolvimento de um procedimento simulador de processo reproduzindo a etapa de destilação extrativa de uma unidade de extração de butadieno a partir de uma corrente de hidrocarbonetos na faixa de quatro átomos de carbono, através da adição do solvente n-metil-2- pirrolidona (NMP). Os resultados obtidos foram comparados e validados com dados de processo obtidos por uma unidade industrial de extração de butadieno. O aprofundamento nos conceitos do processo de separação através de uma ferramenta em simulador de processo capaz de predizer condições de operação permitiu avaliações de aumento de capacidade. A capacidade dos elementos internos dos equipamentos envolvidos na separação pode ser avaliada e a identificação do ponto inicial de engargalamento da unidade foi possível. O procedimento proposto também permite reduzir incertezas para identificação de novos pontos de engargalamento a partir de uma nova configuração dos elementos internos identificados como ineficientes com a elevação de carga processada. / The purpose of this study was the development of a model using a commercial process simulator software in order to reproduce the extractive distillation step of a butadiene extraction unit from a hydrocarbon stream by the addition of the solvent n- methyl-2-pyrrolidone (NMP). The results obtained from the simulations were compared with process data from an operational butadiene extraction unit. The knowledge of the separation process due the model should predict operating conditions allowing revamps on the actual asset. The internals capacity for the equipment involved in the separation could be evaluated and units bottlenecks points should be identified.

A bactéria Zymomonas mobilis atraiu considerável interesse nas últimas décadas devido ao seu metabolismo único e eficientes características fermentativas na produção de etanol a partir de açúcares simples. No entanto, apesar das aparentes vantagens na conversão e taxas específicas quando comparada com as leveduras, estas ainda dominam o mercado produtivo de etanol. Dentre os diversos modelos encontrados na literatura para representar o processo fermentativo de glicose através de Z. mobilis, o modelo proposto por Jöbses et al. (1986) aparenta ser bem estruturado e ideal para a criação de estratégias de controle e otimização de processos fermentativos com este microrganismo, viabilizando sua utilização em escala industrial. Fermentações de Z. mobilis em regime contínuo apresentam comportamento oscilatório para baixas taxas de diluição, fenômeno este comprovado experimentalmente. Experimentos laboratoriais foram utilizados para a estimação de parâmetros do modelo de Jöbses, através da técnica de minimização da derivada do erro da função objetivo. O novo modelo ajustado apresenta dinâmica equivalente ao modelo original de Jöbses, fato este comprovado através da construção dos diagramas de bifurcação. Os diagramas de fase dos modelos apresentam algumas diferenças estruturais entre si, entretanto, a ocorrência de multiplicidade de estados estacionários em baixos valores de taxa de diluição está presente em todos. Em bioprocessos, a obtenção de informações sobre o sistema é um tanto complexa, devido a estes meios serem muito suscetíveis a distúrbios, contaminações, além de que medidas de concentrações são relativamente dispendiosas, associadas com um alto tempo morto, o que impossibilita, muitas vezes, a implementação de sistemas de controle eficientes. Para contribuir com a resolução deste problema foi desenvolvida uma nova inferência das variáveis de estado através da informação de consumo de hidróxido de sódio, utilizado no controle de pH do meio. Inferidores de variáveis de estado que utilizam a variação no consumo de hidróxido de sódio foram ajustados e apresentaram resultados promissores, mostrando a viabilidade do desenvolvimento de metodologias para este tipo de análise, tornando a caracterização de sistemas fermentativos mais rápida e acessível, não apenas em nível acadêmico, que a desenvolvida atualmente, principalmente devido ao baixo custo associado, e a dinâmica rápida deste tipo de sensor. / The bacterium Zymomonas mobilis has attracted considerable interest in recent decades due to their unique metabolism and efficient characteristics in the production of ethanol by simple sugars. However, despite the apparent advantages in the conversion and specific rates, when compared with yeast, there is no industrial-scale fermentations with these bacteria. Among the various models found in literature to represent glucose fermentation with Z. mobilis, the model proposed by Jöbses et al. (1986) appears to be well structured and ideal for the creation of control strategies and optimization methods for increase productions. Fermentation in continuous mode with this organism exhibit oscillatory behavior at low dilution rates, phenomenon showed or demonstrated by experiments performed on laboratory scale. A laboratorial experiment also was used for estimation of new parameters group of Jöbses models, using the technique of minimization of error derivative for objective function. The new adjusted model set presents dynamics similar to the original model, fact confirmed by bifurcation analyses of both models. The diagram for these models show structural differences, but both presents steady states multiplicity and Hopf bifurcations. The obtained of fermentation characteristics is very complex, due to these reactive means being very susceptible to disturbances, contamination, among others, and concentrations measures are relatively expensive and whit a high dead time associated, which often prevents the implementation of control tools. The proposed inference of state variables by consumption rate of hydroxide sodium, used for maintenance pH medium, was tested and confirmed. Virtual sensors using the information of consumption rate of hydroxide were adjusted and have shown promising results, demonstrating the feasibility of developing methods of analysis based on this methodology, making the characterization of fermentative systems faster and cheaper them currently developed, mainly due to low cost associated, and the dynamics of this fast type of sensor.

O presente trabalho investiga a aplicação do processo de extração com fluido supercrítico (EFS) para a obtenção de compostos fenólicos a partir de Butia capitata, coletado na Região Sul do Brasil. A planta, chamada popularmente de butiá, é uma espécie nativa da Argentina, Uruguai, Paraguai e Brasil. Os frutos do butiá são comestíveis e bastante apreciados pelas populações locais, que os consomem ao natural ou na forma de licores, geléias ou doces. Estudos recentes revelam que compostos fenólicos estão presentes nos frutos do butiá. Estes compostos merecem atenção por apresentarem atividade antioxidante, por prevenirem doenças cardiovasculares e reduzirem o risco de câncer. Atualmente existem diferentes formas para a obtenção de compostos fenólicos presentes em matrizes vegetais, porém muitas delas utilizam solventes orgânicos a altas temperaturas, o que deprecia o valor comercial dos extratos. Neste estudo, é utilizada uma tecnologia limpa, que tem como solvente o dióxido de carbono (CO2) em condições supercríticas. O CO2 é não-tóxico e não-inflamável, o que faz dele uma opção mais segura, quando comparado a outros solventes. Além disso, no estado supercrítico, o CO2 apresenta condições bastante favoráveis para a extração, como altas difusividades e massas específicas elevadas. Uma característica que deve ser ressaltada é a possibilidade de utilizar baixas temperaturas na extração, comparativamente aos processos que utilizam solventes líquidos, uma vez que os compostos fenólicos são termolábeis. Nos experimentos, foi utilizada apenas a polpa desidratada. O planejamento experimental foi elaborado avaliando a influência de três variáveis de processo em três níveis: pressão (150, 200 e 250 bar), temperatura (40, 50 e 60°C) e teor de cossolvente (1, 3 e 5% de etanol). Os extratos obtidos foram analisados em cromatografia líquida de alta eficiência (CLAE), sendo os rendimentos expressos em mg do composto identificado por CLAE /100g de polpa de butiá desidratada (PBD), e pelo método Folin-Ciocalteau com os rendimentos em compostos fenólicos totais (CFT) expressos em mg de rutina equivalente (RE)/100g PBD. Foram identificados cinco compostos majoritários nos extratos nas análises em CLAE. As maiores concentrações foram encontradas nas condições de 250 bar, 60°C e 5% de etanol. Os compostos identificados foram o ácido clorogênico (0,46mg/100g PBD), o ácido ferúlico (1,45mg/100g PBD), a rutina (3,47mg/100g PBD), a quercetina (0,90mg/100g PBD) e o canferol (1,93mg/100g PBD). Em termos de CFT, a melhor condição de extração foi de 200 bar, 50°C e 5% de etanol, resultando em 23,09mg RE/100g PBD. Os extratos também foram avaliados quanto ao poder de redução do radical DPPH•, pelo método IC50. O método fornece a quantidade necessária de extrato pra varrer 50% do radical. Como esperado, o extrato obtido na corrida experimental de maior rendimento em CFT apresentou maior atividade antioxidante (17,4 mg/mL), seguido da corrida experimental de maior rendimento em compostos fenólicos isolados (56,1 mg/mL) e, com menos atividade, o extrato obtido por solvente orgânico (80,0 mg/mL). / In this work it is investigated the application of supercritical fluid extraction to obtain phenolic compounds from Butia capitata, collected in southern Brazil. This plant, known popularly as butiá, is a native species in Argentina, Uruguay, Paraguay and Brazil. The edible fruits are popular and eaten raw or as liqueurs, jellies or jams. In recent studies, it is shown that phenolic compounds are present in the fruits of butiá. These compounds deserve attention because they have antioxidant activity, preventing cardiovascular disease and reducing the risk of cancer. There are different processes to obtain phenolic compounds from plant matrices, but many of them make use of organic solvents at high temperatures, which depreciates the extract value. In this study, it is used a clean technology, with carbon dioxide (CO2) as solvent in supercritical conditions. The CO2 is nontoxic and nonflammable that makes it a safer option than other solvents. Furthermore, in supercritical conditions the CO2 has favorable conditions for extraction, as high diffusivities and densities. A feature that should be emphasized is the possibility of using low temperatures in the extract, when compared to processes using liquid solvents, since the phenolic compounds are decomposed by high temperature. In the experiments only the dried pulp was used. The experimental design was developed to assess the influence of three process variables at three levels: pressure (150, 200 e 250 bar), temperature (40, 50 e 60°C) and cosolvent content (1, 3 e 5% of ethanol). The extracts were analyzed by high performance liquid chromatography (HPLC), in which the yields were expressed in terms of mg of the identified phenolic compounds by HPLC/100g of dried pulp and by the Folin-Ciocalteau method, with yields, as total phenolic compounds, expressed in mg of rutin equivalent/100g of dried pulp. It were identified five majority phenolic compounds by HPLC. Its highest concentrations were found in the conditions of 250 bar, 60°C and 5% of cosolvent. The identified compounds were chlorogenic acid (0,46mg/100g of dried pulp), ferulic acid (1,45mg/100g of dried pulp), rutin (3,47mg/100g of dried pulp), quercetin (0,90mg/100g of dried pulp) and kaempferol (1,93mg/100g of dried pulp). In terms of total phenolic compounds, the best conditions were 200 bar, 50°C and 5% of cosolvent, which resulted in 23,09mg of rutin equivalent/100g of dried pulp. Also, the antioxidant activity of the extracts was measured by scavenging of the DPPH• radical, by the IC50 method. The method provides the necessary amount of extract to scavenging 50% of the radial. As expected, the extract obtained in experimental run with the highest yield in CFT showed a higher antioxidant activity (17.4 mg / mL), followed by the experimental run with the highest yield in phenolic compounds isolated (56.1 mg / mL), and, with less activity, the extract obtained by organic solvent (80.0 mg / mL).

Durante a operação de um controlador MPC, a planta pode ir para outro ponto de operação principalmente pela decisão operacional ou pela presença de perturbações medidas/não-medidas. Assim, o modelo do controlador deve ser adaptado para a nova condição de operação favorecendo o controle sob as novas condições. Desta forma, as condições ótimas de controle podem ser alcançadas com a maior quantidade de modelos identificados e com um controlador adaptativo que seja capaz de selecionar o melhor modelo. Neste trabalho é apresentada uma metodologia de controle adaptativo com identificação on-line do melhor modelo o qual pertence a um conjunto previamente levantado. A metodologia proposta considera um controlador em duas camadas e a excitação do processo através de um sinal GBN na camada de otimização com o controlador em malha fechada. Está sendo considerada a validação deste controlador adaptativo através da comparação dos resultados com duas diferentes técnicas Controlador MMPC e Identificação ARX, para a comprovação dos bons resultados desta metodologia. / During the operation of a MPC, the plant can change the operation point mainly due to management decision or due to the presence of measured or unmeasured disturbances. Thus, the model of the controller must be adapted to improve the control in the new operation conditions. In such a way, a better control policy can be achieved if a large number of models are identified at the possible operation points and it is available an adaptive controller that is capable of selecting the best model. In this work is presented a methodology of adaptive control with on-line identification of the most adequate model which belongs to a set of models previously obtained. The proposed methodology considers a two-layer controller and process excitation by a GBN signal in the LP optimization layer with the controller in closed loop mode. It is also presented the adaptive controller validation by comparing the proposed approach with two different techniques - MMPC and ARX Identification, to confirm the good results with this new methodology to the adaptive controller.

Dissertação apresentada para obtenção do grau de Mestre em Engenharia dos Processos Químicos, na Faculdade de Engenharia da Universidade do Porto, sob a orientação do Prof. Doutor José Inácio Martins

A construção de um analisador virtual é sustentada basicamente por três pilares: o modelo, as variáveis que integram o modelo e a estratégia de correção/atualização do modelo. Os modelos matemáticos são classificados quanto ao nível de conhecimento do processo contido nele, indo de modelos complexos baseados em relações fundamentais e leis físico-químicas, denominados white-box, até modelos obtidos através de técnicas de análise multivariável, como técnicas de regressão multiváriavel e redes neurais, referenciados como black box. O presente trabalho objetiva uma análise de dois dos pilares: os modelos, focando em modelos obtidos através das técnicas de redução de dimensionalidade do tipo PLS, e metodologias de seleção de variáveis para a construção dessa classe de modelos. Primeiramente é realizada uma revisão das principais variantes lineares e não lineares da metodologia PLS, compreendendo desde o seu desenvolvimento até a sua combinação com redes neurais. Posteriormente são apresentadas algumas das técnicas popularmente utilizadas para a seleção de variáveis em modelos do tipo black-box, técnicas de validação cruzada e técnicas de seleção de dados para calibração e validação de modelos. São propostas novas abordagens para os procedimentos de seleção de variáveis, originadas da combinação das técnicas de seleção de dados com duas metodologias de seleção de variáveis. Os resultados produzidos por essas novas abordagens são comparados com o método clássico através de casos lineares e não lineares. A viabilidade das técnicas analisadas e desenvolvidas é verificada através da aplicação das mesmas no desenvolvimento de um analisador virtual para uma coluna de destilação simulada através do simulador dinâmico Aspen Dynamics®. Por fim são apresentadas as etapas e desafios da implementação de um analisador virtual baseados em técnicas PLS em uma Torre Depropanizadora de uma central de matérias primas de um pólo petroquímico. / The construction of a virtual analyzer is sustained basically by three pillars: the model, the variables that integrate the model and the updating strategy of the model. The mathematical models are classified with relationship at the level of the process knowledge within it, going from complex models, based on fundamental relationships and physical-chemistries laws, called white-box, until models obtained through multivariable analysis techniques, as multiple linear regression and neural networks, also called as black box. The focus of the present work is the analysis of two of the pillars: the models, specially the ones obtained by dimension reduction techniques, like PLS, and methodologies used in the development of this class of models. Initially, a revision of the main linear and non linear variants of the PLS methodology is done, embracing since its development to its combination with neural networks. Later on, some popularly variables selection techniques for black-box models are explained, as well as some cross validation techniques and strategies for data selection for calibration and validation of models. New approaches for variables selection procedures are proposed, originated by the combination of data selection strategies and two variables selection techniques. The results produced by those new approaches are compared with the classic method through linear and non linear case studies. The viability of the analyzed and developed techniques is verified through the application of the same ones in the development of a virtual analyzer for a distillation column, simulated by the dynamic simulator Aspen Dynamics®. The steps and challenges faced in the implementation of a virtual analyzer based on PLS technical for a Depropanizer Unit are finally presented.

A presente dissertação aborda a otimização de processos industriais através da utilização de projeto de experimentos. Em experimentos planejados, variáveis de respostas são usualmente consideradas como normalmente distribuídas. No entanto, em algumas situações, tal suposição é violada; por exemplo, quando respostas expressam contagens, podendo assumir somente valores inteiros e não-negativos. Nesses casos, é mais provável que as respostas sigam uma distribuição de Poisson e, em sua modelagem, deve-se utilizar os modelos lineares generalizados (MLG), adequados para respostas da chamada família exponencial, à qual pertence a distribuição de Poisson. Se ainda persiste a dúvida quanto ao comportamento das respostas, o modelo de quase-verossimilhança também é uma alternativa possível. Esta dissertação apresenta a reanálise de um experimento, apresentado em Arriba (2005), onde algumas respostas são categóricas. Na análise original do experimento, respostas categóricas foram modeladas através de regressão dos mínimos quadrados ordinários, desconsiderando a violação do pressuposto de normalidade das respostas. Na reanálise aqui apresentada, as variáveis são corretamente abordadas usando-se os modelos lineares generalizados. Como o objetivo do trabalho de Arriba (2005) era a otimização de um processo descrito por múltiplas respostas, comparam-se os resultados da otimização mediante as duas estratégias de modelagem das respostas, além de se propor um método alternativo, mais simplificado, de otimização experimental. / This dissertation deals with the optimization of industrial processes using Design of Experiments. In designed experiments, response variables are often considered as normally distributed. However, in some situations, such assumption is violated; for example, when responses express counts, and only non-negative integers numbers may come up as outcomes. In these cases, it is likely that responses follow a Poisson distribution which is then modeled by generalized linear models (GLM), since such distribution is a member of the exponential family. If a question still holds on the responses behavior, their modeling through the quasi-likelihood method is another option that should be considered. This dissertation analyses an experiment performed by Arriba (2005), where some responses are of categorical nature. In the original analysis presented by Arriba (2005), categorical responses were modeled using ordinary least squares regression, violating the normality assumption associated with that method. In the analysis presented here, variables are appropriately modeled using GLM. Since the objective in Arriba (2005) was optimizing a multiresponse process, results from the two optimization processes are compared. In addition, an alternative experimental optimization method, simpler than the one in Arriba (2005), is also presented.