MODELAGEM E SIMULAÇÃO DA DESTERPENAÇÃO DO ÓLEO DA CASCA DE LARANJA COM CO2 SUPERCRÍTICO EM MODO SEMI-CONTÍNUO

A desterpenação do óleo da casca de laranja com CO2 supercrítico foi investigada através da modelagem e simulação da separação de uma mistura sintética de limoneno (90 % em peso) e linalol (10 %), em um extrator operando em modo semi-contínuo. A modelagem matemática da extração supercrítica foi realizada por analogia com a destilação de uma mistura binária, em batelada, expressando-se a composição das fases em equilíbrio numa base molar livre de CO2. O cálculo das variáveis do processo foi feito por integração numérica da equação de Rayleigh empregando-se o método de Runge-Kutta de quarta ordem. Para a determinação da relação de equilíbrio entre as fases, adotou-se a equação de Peng-Robinson modificada, com os parâmetros de interação obtidos de dados de ELV dos sistemas binários CO2+limoneno e CO2+linalol e do ternário CO2+limoneno+linalol.

Modelagem e simulação dos reatores de polimerização em fase líquida do processo spheripol

O processo Spheripol é, atualmente, um dos mais utilizados na produção industrial de polipropileno e seus copolímeros. É composto por dois reatores tipo loop em série e, em alguns casos, um ou dois reatores tipo fase gás também em série. Neste trabalho é apresentado um modelo para os reatores loop do processo Spheripol. Os reatores loop são modelados como CSTR não ideais. O esquema cinético implementado é relativamente amplo, sendo capaz de englobar copolimerizações multisítio. Uma análise qualitativa do modelo é feita utilizando parâmetros cinéticos presentes na literatura e uma análise quantitativa do modelo é feita com base em dados reais de planta industrial. Ainda com base em dados de planta é feita uma estimação de parâmetros e a validação do modelo. Finalmente, o modelo com parâmetros estimados é aplicado num estudo de caso no qual também é criado um modelo de equilíbrio termodinâmico para a mistura reacional baseado na lei de Henry e na equação de estado de Peng-Robinson.

Recuperação e purificação de quitosanases usando adsorção em leito expandido com streamline DEAE com modelagem e simulação usando redes neurais

Expanded Bed Adsorption (EBA) is an integrative process that combines concepts of chromatography and fluidization of solids. The many parameters involved and their synergistic effects complicate the optimization of the process. Fortunately, some mathematical tools have been developed in order to guide the investigation of the EBA system. In this work the application of experimental design, phenomenological modeling and artificial neural networks (ANN) in understanding chitosanases adsorption on ion exchange resin Streamline® DEAE have been investigated. The strain Paenibacillus ehimensis NRRL B-23118 was used for chitosanase production. EBA experiments were carried out using a column of 2.6 cm inner diameter with 30.0 cm in height that was coupled to a peristaltic pump. At the bottom of the column there was a distributor of glass beads having a height of 3.0 cm. Assays for residence time distribution (RTD) revelead a high degree of mixing, however, the Richardson-Zaki coefficients showed that the column was on the threshold of stability. Isotherm models fitted the adsorption equilibrium data in the presence of lyotropic salts. The results of experiment design indicated that the ionic strength and superficial velocity are important to the recovery and purity of chitosanases. The molecular mass of the two chitosanases were approximately 23 kDa and 52 kDa as estimated by SDS-PAGE. The phenomenological modeling was aimed to describe the operations in batch and column chromatography. The simulations were performed in Microsoft Visual Studio. The kinetic rate constant model set to kinetic curves efficiently under conditions of initial enzyme activity 0.232, 0.142 e 0.079 UA/mL. The simulated breakthrough curves showed some differences with experimental data, especially regarding the slope. Sensitivity tests of the model on the surface velocity, axial dispersion and initial concentration showed agreement with the literature. The neural network was constructed in MATLAB and Neural Network Toolbox. The cross-validation was used to improve the ability of generalization. The parameters of ANN were improved to obtain the settings 6-6 (enzyme activity) and 9-6 (total protein), as well as tansig transfer function and Levenberg-Marquardt training algorithm. The neural Carlos Eduardo de Araújo Padilha dezembro/2013 9 networks simulations, including all the steps of cycle, showed good agreement with experimental data, with a correlation coefficient of approximately 0.974. The effects of input variables on profiles of the stages of loading, washing and elution were consistent with the literature

Determinação de Dados de Equilíbrio Líquido-Vapor para Sistemas Hidrocarbonetos e Desenvolvimento de uma Nova Célula Dinâmica

In the present work are established initially the fundamental relationships of thermodynamics that govern the equilibrium between phases, the models used for the description of the behavior non ideal of the liquid and vapor phases in conditions of low pressures. This work seeks the determination of vapor-liquid equilibrium (VLE) data for a series of multicomponents mixtures of saturated aliphatic hydrocarbons, prepared synthetically starting from substances with analytical degree and the development of a new dynamic cell with circulation of the vapor phase. The apparatus and experimental procedures developed are described and applied for the determination of VLE data. VLE isobarics data were obtained through a Fischer's ebulliometer of circulation of both phases, for the systems pentane + dodecane, heptane + dodecane and decane + dodecane. Using the two new dynamic cells especially projected, of easy operation and low cost, with circulation of the vapor phase, data for the systems heptane + decane + dodecane, acetone + water, tween 20 + dodecane, phenol + water and distillation curves of a gasoline without addictive were measured. Compositions of the equilibrium phases were found by densimetry, chromatography, and total organic carbon analyzer. Calibration curves of density versus composition were prepared from synthetic mixtures and the behavior excess volumes were evaluated. The VLE data obtained experimentally for the hydrocarbon and aqueous systems were submitted to the test of thermodynamic consistency, as well as the obtained from the literature data for another binary systems, mainly in the bank DDB (Dortmund Data Bank), where the Gibbs-Duhem equation is used obtaining a satisfactory data base. The results of the thermodynamic consistency tests for the binary and ternary systems were evaluated in terms of deviations for applications such as model development. Later, those groups of data (tested and approved) were used in the KijPoly program for the determination of the binary kij parameters of the cubic equations of state original Peng-Robinson and with the expanded alpha function. These obtained parameters can be applied for simulation of the reservoirs petroleum conditions and of the several distillation processes found in the petrochemistry industry, through simulators. The two designed dynamic cells used equipments of national technology for the determination Humberto Neves Maia de Oliveira Tese de Doutorado PPGEQ/PRH-ANP 14/UFRN of VLE data were well succeed, demonstrating efficiency and low cost. Multicomponents systems, mixtures of components of different molecular weights and also diluted solutions may be studied in these developed VLE cells

Simulação numérica de escoamentos gás-sólido em leito fluidizado borbulhante utilizando a teoria cinética dos escoamentos granulares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Efeitos numéricos na simulação de escoamentos gás-sólido em leito fuidizado borbulhante utilizando a teoria cinética dos escoamentos granulares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Modelagem e simulação de transferência de calor em leito de partículas estacionárias sujeito a agitação intermitente

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Modelagem do absorvedor e do gerador de ciclos de refrigeração por absorção de calor com o par amônia/água baseados na tecnologia de filme descendente sobre placas inclinadas.

Esse trabalho constitui o desenvolvimento da modelagem térmica e simulação por métodos numéricos de dois componentes fundamentais do ciclo de refrigeração por absorção de calor com o par amônia/água: o absorvedor e o gerador. A função do absorvedor é produzir mistura líquida com alta fração mássica de amônia a partir de mistura líquida com baixa fração mássica de amônia e mistura vapor mediante retirada de calor. A função do gerador é produzir mistura líquido/vapor a partir de mistura líquida mediante o fornecimento de calor. É proposto o uso da tecnologia de filmes descendentes sobre placas inclinadas e o método de diferenças finitas para dividir o comprimento da placa em volumes de controle discretos e realizar os balanços de massa, espécie de amônia e energia juntamente com as equações de transferência de calor e massa para o filme descendente. O objetivo desse trabalho é obter um modelo matemático simplificado para ser utilizado em controle e otimização. Esse modelo foi utilizado para calcular as trocas de calor e massa no absorvedor e gerador para diversas condições a partir de dados operacionais, tais como: dimensões desses componentes, ângulo de inclinação da placa, temperatura de superfície e condições de entrada da fase líquida e vapor. Esses resultados foram utilizados para estabelecer relações de causa e efeito entre as variáveis e parâmetros do problema. Os resultados mostraram que o ângulo de inclinação da placa ótimo tanto para o absorvedor como para o gerador é a posição vertical, ou 90°. A posição vertical proporciona o menor comprimento de equilíbrio (0,85 m para o absorvedor e 1,27 m para o gerador com as condições testadas) e se mostrou estável, pois até 75° não foram verificadas variações no funcionamento do absorvedor e gerador. Dentre as condições testadas para uma placa de 0,5 m verificou-se que as maiores efetividades térmicas no absorvedor e gerador foram respectivamente 0,9 e 0,7 e as maiores efetividades mássicas no absorvedor e gerador foram respectivamente 0,6 e 0,5. É esperado que os dados obtidos sejam utilizados em trabalhos futuros para a construção de um protótipo laboratorial e na validação do modelo.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

This work approaches the forced air cooling of strawberry by numerical simulation. The mathematical model that was used describes the process of heat transfer, based on the Fourier's law, in spherical coordinates and simplified to describe the one-dimensional process. For the resolution of the equation expressed for the mathematical model, an algorithm was developed based on the explicit scheme of the numerical method of the finite differences and implemented in the scientific computation program MATLAB 6.1. The validation of the mathematical model was made by the comparison between theoretical and experimental data, where strawberries had been cooled with forced air. The results showed to be possible the determination of the convective heat transfer coefficient by fitting the numerical and experimental data. The methodology of the numerical simulations was showed like a promising tool in the support of the decision to use or to develop equipment in the area of cooling process with forced air of spherical fruits.

Influência da função distribuição de quebra na modelagem e simulação do circuito de moagem de Carajás - Parte 1 - ajuste de modelos

The perfect mixing model (PMM) is based on parameters derived from the equipment characteristics as well as ore breakage characteristics. Ore characteristics are represented through the appearance function. This function may be determined using JKMRC laboratorial methods or by standard functions. This work describes the model fitting process of the Carajas grinding circuit, using the JKSimMet simulator Two scenarios were used in model fitting exercises: 1) standard appearance function; and 2) appearance fund ion based on testing carried out on samples taken at circuit feed. From this assessment, the appearance function`s influence in the PMM,fit and it`s relation with the breakage rate were determined. The influence of the appearance function on the respective breakage rate distribution was assessed.

Influência da função distribuição de quebra na modelagem e simulação do circuito de moagem de Carajás - parte 2 - simulação

Using the previously described method for appearance function determination, described in Part I of this article, the breakage characterization of the main Carajas ore types was carried out. Based on such characteristics, the ball mill circuit performance was evaluated through simulations. The model described in the first part was used. The results were assessed by comparing ball mill products and cyclone overflow size distribution, as well as simulated recirculating loads. The simulations indicated the potential for processing such ore types at the Carajas grinding circuit, which until now was unknown.

Modelagem de dados de entrada para simulação estocástica de lavra

Mine simulation depends on data that is both coherent and representative of the mining operation. This paper describes a methodology for modeling operational data which has been developed for mine simulation. The methodology has been applied to a case study of an open-pit mine, where the cycle times of the truck fleet have been modeled for mine simulation purposes. The results obtained have shown that once the operational data has been treated using the proposed methodology, the system variables have proven to be adherent to theoretical distributions. The research indicated the need jar tracking the origin of data inconsistencies through the development of a process to manage inconsistent data from the mining operation.

Modelagem e simulação de sistema logístico de distribuição de carne de frango

O sistema logístico para distribuição de produtos acabados caracteriza-se pela integração dos serviços de comunicação, transporte e financeiros com a finalidade de atender às demandas do consumidor final. Estima-se que no estado do Espírito Santo, o consumo de carne de frango seja de 44,4 quilos per capita por ano. Para atender a esta demanda, o estado conta com matadouros-frigoríficos distribuídos pelo seu território, bem como, com a participação de outras empresas localizadas no país. Em sistemas de transportes, são característicos Problemas de Roteamento de Veículos (VRP), que precisam ser estudados, caracterizados e otimizados, normalmente, através de rotinas computacionais, que permitem avaliar maior quantidade de variáveis. O presente trabalho teve por objetivo caracterizar um VRP de um matadouro-frigorífico da região do Sul do Espírito Santo e desenvolver um aplicativo computacional que seja suporte para os gestores de logística, servindo para avaliar e propor rotas, e analisar parâmetros logísticos do processo de distribuição de produtos. No desenvolvimento do aplicativo computacional foi necessário caracterizar o sistema logístico da empresa, coletar e analisar os dados das operações logísticas, desenvolver as rotinas computacionais que representassem o sistema em estudo, verificar a confiabilidade dos resultados fornecidos pelo aplicativo, validá-lo e então, poder realizar as experimentações. O aplicativo desenvolvido permitiu reproduzir dados do sistema estudado e avaliar rotas segundo parâmetros logísticos. Pode-se concluir que o aplicativo computacional desenvolvido é útil aos gestores de logística, permitindo a avaliação das rotas praticadas e de novas configurações de rotas.

Estudo do equilíbrio líquido – líquido em sistemas aquosos bifásicos formados por polietilenoglicol (1500, 3350 e 6000), sais de fosfato e água em diferentes temperaturas

No presente trabalho, foram estudados os dados de equilíbrio de fases dos sistemas aquosos bifásicos formados por polietilenoglicol (1500, 3350 e 6000) + fosfato monobásico e dibásico de sódio (pH 7) + água. Estudou-se o efeito da variação da temperatura (10, 25 e 40°C) bem como da massa molar sobre os dados de equilíbrio. Para o PEG 3350 observou-se um aumento da área bifásica com diminuição da temperatura, mostrando que a formação do sistema aquoso bifásico é exotérmico. Para o PEG 1500 e 6000 houve aumento da área bifásica à 10 e 40°C se comparado à temperatura de 25°C. Em todas as temperaturas em estudo, o aumento da massa molar contribuiu para o aumento da área bifásica. Fato este que foi explicado pelo aumento do caráter hidrofóbico com o aumento da massa do polietilenoglicol.

Modelagem e simulação da dinâmica do arrastamento por ar

Os compostos orgânicos voláteis constituem uma fonte vulgar de contaminação da água subterrânea, a qual pode ser eliminada pela tecnologia do arrastamento por ar (air stripping) em colunas com enchimento desordenado e utilizando fluxos das fases em contra-corrente. Propõe-se neste trabalho uma nova metodologia de dimensionamento destas colunas, para qualquer tipo de enchimento e de contaminante, onde não há necessidade de se arbitrar nenhum diâmetro, onde se evita o recurso a ábacos experimentais e onde o regime hidráulico conveniente é seleccionado à partida. O procedimento proposto foi algoritmizado e convertido num programa em linguagem C++. Para verificar e testar não só o dimensionamento mas também o comportamento teórico estacionário e dinâmico construiu-se de raiz uma coluna experimental. Seleccionou-se como contaminante uma solução de clorofórmio em água destilada. A experimentação permite, ainda, corrigir o coeficiente de transferência de massa global teórico estimado pelas correlações de Onda e que depende de inúmeros parâmetros nem sempre controláveis experimentalmente. Apresenta-se, em seguida, um modelo original de simulação dinâmica do comportamento da coluna e que é constituído por um sistema de equações diferenciais não lineares (parâmetros distribuidos). No entanto, se os débitos forem arbitrados como constantes, o sistema passa a ser linear apesar de não possuir solução analítica evidente (p.e. por transformações integrais). A discretização por diferenças finitas permitiu superar estas dificuldades. Existe uma notável concordância entre os valores experimentais e os previstos no modelo.