Orientadores: João Alexandre Ferreira da Rocha Pereira, Roger Josef Zem / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T18:36:44Z (GMT). No. of bitstreams: 1 BounousRodriguez_Gaston_M.pdf: 3883781 bytes, checksum: 1b217dd247c73d9b0be9947f109696c4 (MD5) Previous issue date: 2003 / Resumo: O etileno é um dos principais intermediários da indústria petroquímica, sendo que sua principal rota de obtenção industrial é através da pirólise de etano. A reação de pirólise é realizada em um fomo que consiste num longo reator tubular colocado no interior de uma fornalha, sendo que a reação é rápida, altamente endotérmica e a distribuição de produtos é fortemente dependente das condições operacionais. O presente trabalho apresenta a modelagem fenomenológica do fomo de pirólise de etano e uma estratégia para o controle do mesmo, utilizando o simulador desenvolvido. A abordagem empregada consiste na resolução dos balanços transientes de massa, energia e queda de pressão no interior do reator tubular em conjunto com o balanço transiente de energia que descreve a transferência de calor da fornalha para o mesmo. O conjunto de equações diferenciais parciais obtido é resolvido mediante a aplicação do Método das Diferenças Finitas e de Runge-Kutta-Fehlberg em um programa computacional desenvolvido em linguagem Fortran. A simulação do modelo proposto permite obter e analisar a dinâmica do fomo de pirólise de etano para diferentes perturbações, assim como os perfis de temperatura, conversão de etano, pressão e distribuição dos produtos ao longo do reator tubular. A seletividade do processo para o etileno, variável importante para a economia do mesmo, é analisada para diferentes condições operacionais, observando-se a dependência da mesma com a carga térmica fornecida a mistura reacional. A estratégia de controle proposta para o fomo de pirólise de etano, consiste em controlar a temperatura dos produtos efluentes do fomo de pirólise, que na industria é um importante indicador do seu funcionamento, atuando na vazão de combustível na fornalha. Inicialmente é utilizada uma técnica de controle convencional empregando-se um controlador PI digital. Em face aos resultados é proposta a aplicação da técnica de controle preditivo por matriz dinâmica ("Dynamic Matrix Contra/" - DMC). A influência dos parâmetros de sintonia no desempenho de ambos controladores é analisada, sintonizando-se os mesmo pelo método da integral do erro absoluto (IAE). Os comportamentos da variável controlada e da variável manipulada obtidos com os controladores PI e DMC propostos são comparados. Observa-se um melhor desempenho do controlador DMC, especialmente devido ao fato de implementar ações de controle menos severas que o controlador PI / Abstract: Ethylene is one of the main intermidiates of the petrochemical industry, being the ethane pyrolisis the main industrial route for its obtention. The pyrolisis reaction takes place in a furnace that consists of a long tubular reactor enclosed inside a combustion chamber, being fast, highly exothermic and whose products distrubition is strongly dependent of the opperational conditions. In the present work is developed a phenomenological modelling of the pyrolisis furnace and a control strategy for it, using the developed simulator. The approach employed consists in the resolution of the transient mass, energy and pressure drop balances inside the tubular reactor along with the transient energy balance which describes the heat transfer between the combustion chamber and the reactor. The obtained set of partial differential equations is solved by the apllication of the Finite Differences Method and Runge-Kutta-FeWberg using a developed computer program in Fortran. The simulation of the proposed model allows to obtain and analyse the furnace dynamic response to different loads, as well as the temperature, ethane conversion, pressure and products distribution profiles along the tubular reactor. The selectivity of the process to ethylene, which is an important variable for the process economy, is analized for differents operational conditions, observing its dependence with the thermal load supplied to the reactional mixture. The proposed control strategy for the pyrolisis furnace consists in controlling the temperature of the eff1uent stream from the furnace, witch is an important industrial indicator of its performance, manipulating the fuel flow to the combustion chamber. First1y, a conventional control technique is applied using a digital PI controller. Then, it is applied the predictive control technique known as Dynamic Matrix Control (DMC). The influence of the tunning parameters in the performance of both controllers is analyzed, tunning them using the integral of the absolute value of the error method (IAE). The behaviours of the controlled and manipulated variables obtained with the PI and DMC controllers are compared. A better performance of the DMC controller can be observed, especially due to the fact that it applies control actions less severe than the PI controller / Mestrado / Sistemas de Processos Quimicos e Informatica / Mestre em Engenharia Química

Orientador: Liliane M. Ferrareso Lona / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T18:36:21Z (GMT). No. of bitstreams: 1 RoseroMesa_PaolaViviana_M.pdf: 5952280 bytes, checksum: 09b01c247830a67ec0b6df187b843097 (MD5) Previous issue date: 2003 / Resumo: Neste trabalho foi desenvolvido um modelo matemático para a polimerização via radical livre controlada em presença de radicais nitróxidos (NMRP), que inclui a obtenção dos pesos moleculares médios e o índice de polidispersidade para polímeros de estrutura linear. Por ser a polimerização via radical livre (FRP) uma das polimerizações mais estudadas, este trabalho contemplou o estudo do processo FRP em sistemas homogêneos como base para a modelagem da polimerização via radical livre controlada (LFRP). A LFRP se destaca por ser capaz de obter polímeros com polidispersidade perto de 1,0 e distribuições de peso molecular estreitas. Existem diferentes mecanismos LFRP, dos quais se destaca o mecanismo NMRP, o qual consegue obter polímero com peso molecular controlado mediante a adição de um radical nitróxido à polimerização via radical livre, gerando uma reação reversível capaz de bloquear as cadeias em crescimento. Para obtenção dos pesos moleculares médios e a distribuição do peso molecular (MWD) no processo FRP, foram utilizados o método instantâneo e o método dos momentos; já no novo processo (NMRP) foi usado o método dos momentos para a obtenção dos pesos moleculares médios e o índice de polidispersidade. O modelo matemático foi baseado nas constantes cinéticas de taxa e na teoria do volume livre para o processo FRP e nas constantes cinéticas no processo NMRP. Os resultados foram comparados com dados experimentais para avaliar o modelo. Como casos estudo foram analisadas as polimerizações de estireno e de metilmetacrilato para o processo FRP e as polimerizações de estireno no processo NMRP. Os resultados de validação obtidos mostraram que o modelo sugerido é capaz de modelar o processo monomolecular e bimolecular NMRP a diferentes condições de temperatura e concentração de espécies. Os resultados do modelo também mostraram como no processo NMRP, a espécie polimérica morta surge de radicais poliméricos que não conseguem ser capturados no processo. Foi mostrado como radicais nitróxidos menos estáveis geram resultados mais favoráveis em termos de produtividade e qualidade de produto / Abstract: A practical methodology for the computer modeling of number-and weight-averages molecular weight and molecular weight distribution of polymer with linear structure obtained from nitroxide mediated radical polymerization (NMRP) has been developed. Since free-radical polymerization (FRP) is one of the most important chain-growth polymerization, this work embraces a free-radical polymerization in homogeneous systems as reference point to model the "Iiving" free-radical polymerization. The "living" free-radical polymerization (LFRP) emphasizes the production of polymers with highly controlled structures, narrows MWD and polydispersity index closed to 1.0. Different mechanisms have been proposed to carry out LFRP, among them, Nitroxide-mediated radical po/ymerization (NMRP), which attains polymers with controlled molecular weight by mediated nitroxide radicais acting as trapping agent, in which the polymer chain that is growing are reversibly trapped. To obtain number- and weight-averages molecular weight and molecular weight distributions (MWD) for FRP process the instantaneous and moments methods were used. On the other hand, for the process NMRP the method of moments was used to obtain the number- and weight-averages molecular weight and polydispersity index. The mathematical model was based on kinetics constant rates and free-volume theory for FRP and kinetics constant rates for NMRP. Simulations were compared with experimental data to validate the model. Styrene and methylmethacrilate was analyzed as cases studies for FRP process and styrene polymerization for NMRP process. The model was able to modeling monomolecular and bimolecular NMRP process to different conditions of temperature and species concentrations. It showed the dead polymer chain appear by live polymeric chain not captured by nitroxide radicais in the NMRP processo. Then, NMRP can obtain "Iiving" characteristic if exists enough quantity of nitroxide radical. This permits the radicais created to be captured. Nitroxide radicais more unstable generate best results in relation to productivity and product quality / Mestrado / Mestre em Engenharia Química

Orientador: Sergio Tonini Button / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T19:44:10Z (GMT). No. of bitstreams: 1 Yamakami_WyserJose_D.pdf: 20475770 bytes, checksum: 04886753497f892b0d4f876e5e7de310 (MD5) Previous issue date: 2003 / Doutorado

Orientador: Franco Giuseppe Dedini / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T20:38:36Z (GMT). No. of bitstreams: 1 Forti_AntonioWagner_D.pdf: 2288454 bytes, checksum: edbd3279e8fdb2ec693dc723ba141b2b (MD5) Previous issue date: 2003 / Doutorado

Orientadores: Denis Jose Schiozer, Eliana Luci Ligero / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica, Instituto de Geociencias / Made available in DSpace on 2018-08-03T22:12:10Z (GMT). No. of bitstreams: 1 Sanjombi_BernabeAlfredo_M.pdf: 7416551 bytes, checksum: 9fa54684138d500059c70748ad19b2a8 (MD5) Previous issue date: 2004 / Resumo: Técnicas de transferência de escala permitem que propriedades de uma malha refinada sejam adaptadas para uma malha de simulação mais grosseira. Em casos com fortes heterogeneidades, a passagem de uma escala refinada para outra mais grosseira pode fazer com que resultados bem diferentes sejam obtidos. Outro ponto importante é que a transferência de escala de blocos finos com poços pode ter grande impacto nos resultados e, em muitos casos, deve ter tratamento diferenciado e especial. Deste modo, para uma melhor aproximação entre as malhas refinada e grosseira, surge a necessidade do emprego de tratamentos especiais aos blocos grosseiros com fortes heterogeneidades e com poços. O objetivo deste trabalho é estudar a transferência de escala da permeabilidade absoluta de casos com canais e barreiras e aproximar os resultados da simulação numérica de fluxo da malha grosseira ao da malha refinada utilizando tratamentos de poços nos blocos grosseiros. A partir dos resultados obtidos foi possível selecionar alguns tratamentos de poços e generaliza-los. Como os resultados dependem tanto do tipo de tratamento de poço, quanto do tipo de heterogeneidade e tipo de fluxo, elaborou-se uma etapa onde a malha grosseira foi gerada pela combinação de tipos de fluxo com tipos de tratamentos de poços de acordo com a localização do poço na malha refinada. Este processo mostrou ser adequado para todos os casos estudados, apesar de exigir um elevado tempo na sua elaboração. Com a finalidade de redução deste tempo, foi realizada uma generalização de técnicas de transferência de escala e de tipos de tratamentos de poços para cada condição de fluxo e heterogeneidade. Esta regra geral foi aplicada com sucesso nos casos estudados. Considerando-se todos os resultados gerados pode-se concluir quais os tipos de tratamentos de poços mais adequados para malhas refinadas com canais e barreiras. Também foi possível indicar tipos de tratamentos de poços inadequados aos casos estudados. Quanto a técnica de transferência de escala mais adequada para cada condição de fluxo e tipo de heterogeneidade pode-se observar que as mesmas dependeram. das características dos poços. Nos casos com canais as técnicas de transferência de escala mais adequadas dependeram. do fato do poço produtor estar alinhado ou desalinhado aos canais. Quanto aos casos com barreiras as técnicas indicadas estão associadas ao número de camadas com barreiras nas quais os poços produtores estão perfurados / Abstract: UpscaIing techniques allow the adaptation of fine grid properties to coarse simulation grids. In cases with strong heterogeneities, the flow simulation of coarse grids could produce very different results from those obtained with fine grids. An important point is the upscaling of fine blocks with wells, which yields a great impact on simulation results ando in many cases, demands a differentiated and special treatment. To obtain a better approximation between coarse and fine grids, it is necessary to employ special treatments to coarse blocks with strong heterogeneity and with well. The goal of this work is to study absolute permeability upscaling in cases with channels and barriers and to approximate the flow simulation results of coarse grid to those obtained for fine grids by using special well block treatment It was possible to select, from the obtained results, some well treatments and generalize them. As the results depend on the type of well treatment, of heterogeneity and of flow, a coarse grid was built via combination of flow types and with types of well treatments according to the welllocalization in the fine grid This process was appropriate for all studied cases, in spite of demanding a long time to elaborate the coarse grid In order to reduce the process time, a genera1ization of upscaling techniques and of types of well treatments according to flow and heterogeneity conditions was realized. This general role was applied to studied cases successfully. Considering all the generated results, it was possible to obtain the most adequate types of well treatments to fine grids with channel and baniers. Also, it was possible to select the inadequate well treatments to the studied cases. It was observed the most adequate upscaling techniques for each type of flow and of heterogeneity depended on the well characteristics. In cases with channels, the upscaling technique depended on the welllocalization being aligned or not to the channels. In the case with baniers, the indicated upscaling techniques were associated to the number of layers with baniers where the producer well is drilled. / Mestrado / Explotação / Mestre em Ciências e Engenharia de Petróleo

Orientador : Kamal A. R. Ismail / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T22:46:33Z (GMT). No. of bitstreams: 1 SalinasSedano_CarlosTeofilo_D.pdf: 9310832 bytes, checksum: 9f979c379774f591e2629bbb81d97e7b (MD5) Previous issue date: 2003 / Resumo: A radiação tem um papel muito importante na transferência de energia em variadas aplicações científicas e da engenharia. A determinação dos fluxos radiativos nas superfícies e os termos radiativos volumétricos em meios contendo gases reais é dificil devido à forte dependência do coeficiente de absorção destes gases no número de onda. Um dos objetivos deste trabalho é desenvolver e aplicar um método para a solução espectral da equação de transferência de calor por radiação e condução combinadas em misturas de gases não cinzentos em casos homogêneos e não homogêneos e ou outro é sua aplicação na solução de problemas práticos de aplicação na pesquisa científica e na engenharia. Aplica-se especificamente neste trabalho na modelagem térmica de janelas térmicas com gases absorventes de radiação inftavermelha. Na solução da equação de transporte radiativo é utilizado o método das ordenadas discretas. Proporse um método eficiente e capaz de eliminar o efeito raio em complexas situações bidimensionais, validando o algoritmo em variadas situações bidimensionais. Depois de alcançar esta meta, este procedimento e algoritmo serão utilizados no cálculo de transferência de calor por condução e radiação acopladas, validando os algoritmos por comparação com soluções de casos teste padrão obtidos e utilizados na literatura. O passo seguinte é a utilização desses algoritmos numéricos em cada etapa da solução numérica dos casos complexos de transferência de calor por radiação em gases reais ou não cinzentos e em misturas destes gases onde tem que se fazer a integração espectral das equações de transporte radiativo. Adota-se um modelo de integração espectral para gases reais baseado no conceito de correlação de espectro local que pode tratar de maneira geral, misturas de gases reais, não isotérmicos, não homogêneos e também com particulados tais como fuligem. O modelo de gás real e o algoritmo numérico são validados pela aplicação do modelo em casos teste na literatura. Os métodos de solução geral da equação de transporte radiativo acoplado ao mecanismo de transferência de calor por condução e dos modelos de integração espectral para misturas de gases reais que são desenvolvidas neste trabalho mostram-se úteis para o tratamento de inúmeros problemas que vão desde os problemas de combustão com gases homogêneos e não homogêneos e foram utilizados aqui para a validação dos modelos de gases reais, problemas a temperaturas ambientais e outros que não foram tratados aqui. Depois de ter-se validados o modelo de integração espectral da equação de transporte radiativo aplica-se o modelo na modelagem numérica da transferência de calor por radiação e condução combinadas em janelas térmicas com gases absorventes de radiação inftavermelha para casos bidimensionais. O desempenho térmico da janela é avaliado / Abstract: Determination of superficial radiative heat flux and volumetric source tenn in real gas media is difficult due the strong dependency of absorption coefficient to wavenumber. One objective of this work is to develop and to apply one method for spectrally radiation-conduction heat transfer solution in homogeneous and non-homogeneous real gases mixtures. Another objective is to apply the method for practical problems in engineering and scientific research. In this work, the method was applied to thermal modeling of a glass window filled with infrared absorbing gases. Radiative transport equation solution is obtained using the discrete ordinates method. An eflicient method for eliminate the. ray effect in complex two-dimensional situations is proposed and the algorithm is validated against different two-dimensional cases. The procedure and algorithm are used in coupled radiation-conduction heat transfer problems. Comparison of results with literature benchmark test cases is done for validation. The next step is to use the numerical -algorithms in the solution for complex cases of radiative heat transfer in non-gray gases and in non-gray gases mixtures where spectrally integration of radiative transport equations is needed. The model of local spectrum correlation is adopted for spectral integration of radiative transport equation. This model supports non-isothermal and non-homogeneous real gases mixtures. Also, it supports gases mixtures with addition of soot. The real gas model and the numeric algorithm are validated for comparison with results taken trom the literature. Finally, the development method is applied for to numerical modeling of the coupled radiation-conduction heat transfers in glazing filled with infrared absorbing gases mixtures. The thermal performance of the glazing is calculated. / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Made available in DSpace on 2014-06-12T17:43:01Z (GMT). No. of bitstreams: 2 arquivo9603_1.pdf: 12819797 bytes, checksum: 7102fe5a2da0f64aa3f3570246985c2a (MD5) license.txt: 1748 bytes, checksum: 8a4605be74aa9ea9d79846c1fba20a33 (MD5) Previous issue date: 2012 / A Indústria Naval no Brasil adormeceu por um período de vinte anos, com a descoberta do Pré-Sal e a necessidade de desenvolvimento do setor de Óleo e Gás no país, a atividade tem crescido vertiginosamente. Devido à escassez de dados históricos sobre a construção naval, este trabalho se propõe a apresentar um modelo para análise de risco financeiro através de um projeto de construção de um navio petroleiro do tipo Suezmax. O modelo tem natureza probabilística, utilizando a metodologia de Simulação de Monte Carlo, pois leva em consideração a incerteza e complexidade do ambiente observando diferentes graus de risco através de variáveis de saída: Valor Presente Líquido (VPL) e Taxa Interna de Retorno (TIR). As variáveis de saída são geradas a partir de dados que compõe o projeto da embarcação e, em especial, variáveis aleatórias que representam impacto significativo no resultado final da análise de risco. Para aplicar a Simulação de Monte Carlo, foi necessária a utilização dos softwares Microsoft Excel® e Crystal Ball®. Os resultados do VPL e TIR mostram valores bastante próximos, confirmando uma forte correlação entre os dados estatísticos analisados. Em todas as distribuições probabilísticas realizadas, não geraram VPL negativo, porém, em um dos cenários, houve probabilidade de haver TIR negativa. A utilização de Simulação de Monte Carlo em análise de risco financeiro de projetos da Indústria Naval pode ser uma ferramenta útil para gerar valor agregado nas projeções nos projetos de investimento do setor, dando maior conforto e segurança ao decisor

Made available in DSpace on 2014-06-12T18:02:11Z (GMT). No. of bitstreams: 2 arquivo124_1.pdf: 2840973 bytes, checksum: 6824639178df3434cffa33b46a79eae2 (MD5) license.txt: 1748 bytes, checksum: 8a4605be74aa9ea9d79846c1fba20a33 (MD5) Previous issue date: 2011 / Conselho Nacional de Desenvolvimento Científico e Tecnológico / Nos últimos anos a indústria sucroalcooleira tem apresentado um crescimento significativo. Em parte, esse fenômeno se deve a grande valorização dos biocombustíveis, no caso em questão o etanol. Dentro desse panorama a otimização e integração energéticas se tornam cada vez mais importantes nessas indústrias. Os combustíveis de segunda geração, produzidos através do processo de hidrólise dos materiais lignocelulósicos, apresentam-se como uma grande alternativa para o uso energético da biomassa. Neste trabalho foi proposto um estudo de caso, a fim de construir uma planta virtual do processo de produção de etanol a partir da hidrólise do bagaço de cana. O processo de transformação de biomassa em combustível consiste das seguintes etapas: Pré-tratamento, deslignificação, hidrólise, fermentação e destilação. A partir das simulações, realizadas no software HYSYS®, foi possível determinar o consumo de insumos e de energia necessários na planta. Para avaliar o aumento da produção de etanol dentro de uma usina foi realizada a integração do processo de segunda geração (processo de hidrólise do bagaço) e o processo de primeira geração (a partir do processamento da cana-de-açúcar), através da adição da corrente de mosto à corrente do licor de hexoses gerado no processo de hidrólise. Dos resultados obtidos verificou-se que ao se utilizar 80% de bagaço (cogerado nas moendas) como matéria-prima no processo de hidrólise foi possível aumentar a produção de etanol de cerca de 26% (em volume). Também foi realizada a análise exergética na unidade de produção de etanol com a finalidade de localizar as irreversibilidades dos sistemas, identificando as exergias destruídas e assim avaliando as eficiências. A partir dos resultados dessa análise foram propostas alternativas para aperfeiçoar o desempenho dos sistemas e de melhor aproveitar os recursos utilizados

Made available in DSpace on 2014-06-12T18:03:07Z (GMT). No. of bitstreams: 2 arquivo172_1.pdf: 3067231 bytes, checksum: 51cdafaa87f5a655f55051d85f5aa95b (MD5) license.txt: 1748 bytes, checksum: 8a4605be74aa9ea9d79846c1fba20a33 (MD5) Previous issue date: 2011 / Agência Nacional do Petróleo, Gás Natural e Biocombustíveis / Dentro do panorama atual, há uma demanda crescente pelo processamento de óleos pesados, devido principalmente a escassez de novas reservas de óleos leves. Sob esta ótica o processamento térmico de frações pesadas tem sido uma das formas encontradas para obter compostos mais leves que, em geral, possuem uma maior demanda tanto para combustíveis como para a petroquímica, assim como um maior valor agregado. Sob esses aspectos torna-se interessante o desenvolvimento de modelos matemáticos e simulação computacional para a melhor compreensão de como se dá o processo de craqueamento térmico em uma fornalha. Através do desenvolvimento de um código em linguagem Matlab® computacional foi possível prever a conversão da carga alimentada em diversas frações de petróleo usando pseudocomponentes e lumps para tratar as misturas respectivamente do ponto de vista termodinâmico e cinético. No presente trabalho foram usadas equações diferenciais ordinárias para representar o balanço de massa de um processo de craqueamento térmico, que foram resolvidas através da técnica das diferenças finitas. As diversas propriedades necessárias para modelar a perda de carga, hold up e equilíbrio termodinâmico foram obtidas através de correlações empíricas presentes na literatura ou nos procedimentos da API. A partir dos resultados das simulações foram obtidos perfis de concentração e de queda de pressão ao longo do equipamento

Made available in DSpace on 2014-06-12T18:04:23Z (GMT). No. of bitstreams: 2 arquivo213_1.pdf: 4056030 bytes, checksum: c7d78e098b753c3c2b6ce96fc5f9588a (MD5) license.txt: 1748 bytes, checksum: 8a4605be74aa9ea9d79846c1fba20a33 (MD5) Previous issue date: 2011 / Universidade Federal de Pernambuco / Nos últimos anos o biodiesel se tornou uma alternativa para a demanda crescente de combustível. O próximo passo é conseguir produzir um biodiesel economicamente competitivo com o diesel fóssil em um processo em nível industrial. Este trabalho visa estudar do ponto de vista computacional, uma planta piloto de biodiesel, simulando o processo desde a reação de transesterificação de óleos vegetais, chegando até a purificação do biodiesel, utilizando o APEN/HYSYS. Além disso, foi feito o estudo da viabilidade econômica preliminar da mesma, fazendo-se uso do custo anualizado total unitário CATU. Os resultados das simulações foram comparados com os resultados obtidos numa planta piloto montada em Pernambuco/Brasil. Em seguida, foi comparada a viabilidade econômica da planta piloto, com uma planta operando com uma coluna de destilação reativa para produção do referido combustível. Os resultados mostraram que a destilação reativa é um processo mais econômico para a produção do biodiesel do que em um processo em batelada