Avaliação da capacidade de adsorção do CO2 em zeólita 12X com gases sintéticos e originados da pirólise do lodo de esgoto

This work depicts a study of the adsorption of carbon dioxide on zeolite 13X. The activities were divided into four stages: study batch adsorption capacity of the adsorbent with synthetic CO2 (4%), fixed bed dynamic evaluation with the commercial mixture of gases (4% CO2, 1.11% CO, 1 2% H2, 0.233% CH4, 0.1% C3, 0.0233% C4 argon as inert closing balance), fixed bed dynamic modeling and evaluation of the breakthrough curve of CO2 originated from the pyrolysis of sewage sludge. The sewage sludge and the adsorbent were characterized by analysis TG / DTA, SEM, XRF and BET. Adsorption studies were carried out under the following operating conditions: temperature 40 °C (for the pyrolysis of the sludge T = 600 °C), pressures of 0.55 to 5.05 bar (batch process), flow rate of the gaseous mixture between 50 - 72 ml/min and the adsorbent masses of 10, 15 and 20 g (fixed bed process). The time for the adsorption batch was 7 h and on the fixed bed was around 180 min. The results of this study showed that in batch adsorption process step with zeolite 13X is efficient and the mass of adsorbed CO2 increases with the increases pressure, decreases with temperature increases and rises due the increase of activation temperature adsorbent. In the batch process were evaluated the breakthrough curves, which were compared with adsorption isotherms represented by the models of Langmuir, Freündlich and Toth. All models well adjusted to the experimental points, but the Langmuir model was chosen in view of its use in the dynamic model does not have implications for adsorption (indeterminacy and larger number of parameters such as occurred with others) in solving the equation. In the fixed bed dynamic study with the synthetic gas mixture, 20 g of mass adsorbent showed the maximum adsorption percentage 46.7% at 40 °C temperature and 50 mL/min of flow rate. The model was satisfactorily fitted to the three breakthrough curves and the parameters were: axial dispersion coefficient (0.0165 dm2/min), effective diffusivity inside the particle (dm2/min 0.0884) and external transfer coefficient mass (0.45 dm/min). The breakthrough curve for CO2 in the process of pyrolysis of the sludge showed a fast saturation with traces of aerosols presents in the gas phase into the fixed bed under the reaction process

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

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

Evolução da biodegradabilidade da matéria orgânica em um sistema de lagoas de estabilização

Waste stabilization ponds (WSP) have been widely used for sewage treatment in hot climate regions because they are economic and environmentally sustainable. In the present study a WSP complex comprising a primary facultative pond (PFP) followed by two maturation ponds (MP-1 and MP-2) was studied, in the city of Natal-RN. The main objective was to study the bio-degradability of organic matter through the determination of the kinetic constant k throughout the system. The work was carried out in two phases. In the first, the variability in BOD, COD and TOC concentrations and an analysis of the relations between these parameters, in the influent raw sewage, pond effluents and in specific areas inside the ponds was studied. In the second stage, the decay rate for organic matter (k) was determined throughout the system based on BOD tests on the influent sewage, pond effluents and water column samples taken from fixed locations within the ponds, using the mathematical methods of Least Squares and the Thomas equation. Subsequently k was estimated as a function of a hydrodynamic model determined from the dispersion number (d), using empirical methods and a Partial Hydrodynamic Evaluation (PHE), obtained from tracer studies in a section of the primary facultative pond corresponding to 10% of its total length. The concentrations of biodegradable organic matter, measured as BOD and COD, gradually reduced through the series of ponds, giving overall removal efficiencies of 71.95% for BOD and of 52.45% for COD. Determining the values for k, in the influent and effluent samples of the ponds using the mathematical method of Least Squares, gave the following values respectively: primary facultative pond (0,23 day-1 and 0,09 day-1), maturation 1 (0,04 day-1 and 0,03 day-1) and maturation 2 (0,03 day-1 and 0,08 day-1). When using the Thomas method, the values of k in the influents and effluents of the ponds were: primary facultative pond (0,17 day-1 and 0,07 day-1), maturation 1 (0,02 day-1 and 0,01 day-1) and maturation 2 (0,01 day-1 and 0,02 day-1). From the Partial Hydrodynamic Evaluation, in the first section of the facultative pond corresponding to 10% of its total length, it can be concluded from the dispersion number obtained of d = 0.04, that the hydraulic regime is one of dispersed flow with a kinetic constant value of 0.20 day-1

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Estudo da perda de carga e calor no poço injetor no processo De drenagem gravitacional assistida por vapor (SAGD)

The oil companies in the area in general are looking for new technologies that can increase the recovery factor of oil contained in reservoirs. These investments are mainly aimed at reducing the costs of projects which are high. Steam injection is one of these special methods of recovery in which steam is injected into the reservoir in order to reduce the viscosity of the oil and make it more mobile. The process assisted gravity drainage steam (SAGD) using steam injection in its mechanism, as well as two parallel horizontal wells. In this process steam is injected through the horizontal injection well, then a vapor chamber is formed by heating the oil in the reservoir and, by the action of gravitational forces, this oil is drained down to where the production well. This study aims to analyze the influence of pressure drop and heat along the injection well in the SAGD process. Numerical simulations were performed using the thermal simulator STARS of CMG (Computer Modeling Group). The parameters studied were the thermal conductivity of the formation, the flow of steam injection, the inner diameter of the column, the steam quality and temperature. A factorial design was used to verify the influence of the parameters studied in the recovery factor. We also analyzed different injection flow rates for the model with pressure drop and no pressure drop, as well as different maximum flow rates of oil production. Finally, we performed an economic analysis of the two models in order to check the profitability of the projects studied. The results showed that the pressure drop in injection well have a significant influence on the SAGD process.

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developing the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. It s important to point out that, in spite of the loads being normally connected to the transformer s secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity

Uso de técnicas de otimização baseadas em derivadas como suporte do planejamento operacional de redes de distribuição de energia elétrica

The usual programs for load flow calculation were in general developped aiming the simulation of electric energy transmission, subtransmission and distribution systems. However, the mathematical methods and algorithms used by the formulations were based, in majority, just on the characteristics of the transmittion systems, which were the main concern focus of engineers and researchers. Though, the physical characteristics of these systems are quite different from the distribution ones. In the transmission systems, the voltage levels are high and the lines are generally very long. These aspects contribute the capacitive and inductive effects that appear in the system to have a considerable influence in the values of the interest quantities, reason why they should be taken into consideration. Still in the transmission systems, the loads have a macro nature, as for example, cities, neiborhoods, or big industries. These loads are, generally, practically balanced, what reduces the necessity of utilization of three-phase methodology for the load flow calculation. Distribution systems, on the other hand, present different characteristics: the voltage levels are small in comparison to the transmission ones. This almost annul the capacitive effects of the lines. The loads are, in this case, transformers, in whose secondaries are connected small consumers, in a sort of times, mono-phase ones, so that the probability of finding an unbalanced circuit is high. This way, the utilization of three-phase methodologies assumes an important dimension. Besides, equipments like voltage regulators, that use simultaneously the concepts of phase and line voltage in their functioning, need a three-phase methodology, in order to allow the simulation of their real behavior. For the exposed reasons, initially was developped, in the scope of this work, a method for three-phase load flow calculation in order to simulate the steady-state behaviour of distribution systems. Aiming to achieve this goal, the Power Summation Algorithm was used, as a base for developping the three phase method. This algorithm was already widely tested and approved by researchers and engineers in the simulation of radial electric energy distribution systems, mainly for single-phase representation. By our formulation, lines are modeled in three-phase circuits, considering the magnetic coupling between the phases; but the earth effect is considered through the Carson reduction. Its important to point out that, in spite of the loads being normally connected to the transformers secondaries, was considered the hypothesis of existence of star or delta loads connected to the primary circuit. To perform the simulation of voltage regulators, a new model was utilized, allowing the simulation of various types of configurations, according to their real functioning. Finally, was considered the possibility of representation of switches with current measuring in various points of the feeder. The loads are adjusted during the iteractive process, in order to match the current in each switch, converging to the measured value specified by the input data. In a second stage of the work, sensibility parameters were derived taking as base the described load flow, with the objective of suporting further optimization processes. This parameters are found by calculating of the partial derivatives of a variable in respect to another, in general, voltages, losses and reactive powers. After describing the calculation of the sensibility parameters, the Gradient Method was presented, using these parameters to optimize an objective function, that will be defined for each type of study. The first one refers to the reduction of technical losses in a medium voltage feeder, through the installation of capacitor banks; the second one refers to the problem of correction of voltage profile, through the instalation of capacitor banks or voltage regulators. In case of the losses reduction will be considered, as objective function, the sum of the losses in all the parts of the system. To the correction of the voltage profile, the objective function will be the sum of the square voltage deviations in each node, in respect to the rated voltage. In the end of the work, results of application of the described methods in some feeders are presented, aiming to give insight about their performance and acuity

Estudo numérico de escoamento turbulento em padrão anular gás-líquido em dutos verticais

Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm

Decaimento bacteriano em lagoas de estabilização no Nordeste brasileiro

Stabilization pond system consisting in more sewage treatment used in Rio Grande do Norte (RN), Brazil, representing about 90% of all systems. Fecal bacteria are removed mainly facultative ponds and in maturation ponds. Many factors influence bacterial decay, such as the levels of pH and DO, temperature, light intensity, HDT and nutrient availability. The bacterial decay rate (Kb) is calculated considering many variables, but the hydraulic regime is a significant influence for microorganisms removal, and the dispersed flow which best characterizes a stabilization pond. However, some authors developed equations for the Kb accordant plug flow and complete mixing. This research study aimed to evaluate the bacterial decay of fecal coliform and Enterococcus sp. in stabilization ponds designed to treat domestic sewage, full-scale, in RN. All systems have assessed pretreatment, a facultative pond (LF) followed by two maturation (LM1 and LM2). The parameters availed were: temperature, pH, DO, BOD5, COD, fecal coliform, Enterococcus sp., Chlorophyll a, total suspended solids, fixed and volatile. In general, there were not significant differences for pH, DO and temperature in the ponds, except for the new systems, since they have low flow and hydraulic loads. The removal of organic matter in the ponds was low, about 70%, and nearly all are overloaded organic and operational problems. The bacterial removals were low, with average 96% for LF for fecal coliform, and 98% for Enterococcus sp.; LM1 were in itself a removal for fecal coliform about 71%, and 81% for Enterococcus sp.; LM2 have efficiency of 69% for fecal coliform, and 68% for Enterococcus sp. The equation proposed by Von Sperling (1999), according to the dispersed flow regime, generated empirical values of Kb more approximate to calculated values of Kb. On average, the calculated Kb to coliforms in the LF was 0.31 d-1, and for both maturation ponds were 0.35 d-1. For Enterococcus sp. the average was 0.40 d-1 for LF, 0.55 d-1 for LM1, and 0.58 d-1 for LM2. These results also showed that the Kb obtained in full-scale systems are smaller than those found in pilot-scale ponds. Moreover, one can say that the equation proposed by Marais (1974), according to the complete-mix regime, overestimates Kb. Actual results of Kb indicated that fecal coliforms are more resistant to adverse conditions present in stabilization ponds than Enterococcus sp., therefore, an indicator of microbiological safety and efficiency. The factors significant interventions in the rate of bacterial decay were concentrations of COD, the organic loading and HDT. The few Kb relationship between pH, DO and temperature were not significant. Finally, we conclude that it s essential to correct operation and maintenance, for not performing these activities is one of the main factors contributing to low rates of bacterial decay.

Estudo de manejo de águas pluviais urbanas na cidade Natal-Rio Grande do Norte

Urban stormwater can be considered as potential water resources as well as problems for the proper functioning of the manifold activities of the city, resulting from inappropriate use and occupation of the soil, usually due to poor planning of the occupation of the development areas, with little care for the environmental aspects of the drainage of surface runoff. As a basic premise, we must seek mechanisms to preserve the natural flow in all stages of development of an urban area, preserving the soil infiltration capacity in the scale of the urban area, comprising the mechanisms of natural drainage, and noting preserving natural areas of dynamic water courses, both in the main channel and in the secondary. They are challenges for a sustainable urban development in a harmonious coexistence of modern developmental, which are consistent with the authoritative economic environmental and social quality. Integrated studies involving the quantity and quality of rainwater are absolutely necessary to achieve understanding and obtaining appropriate technologies, involving both aspects of the drainage problems and aspects of use of water when subjected to an adequate management of surface runoff , for example, the accumulation of these reservoirs in detention with the possibility of use for other purposes. The purpose of this study aims to develop a computer model, adjusted to prevailing conditions of an experimental urban watershed in order to enable the implementation of management practices for water resources, hydrological simulations of quantity and, in a preliminary way, the quality of stormwater that flow to a pond located at the downstream end of the basin. To this end, we used in parallel with the distributed model SWMM data raised the basin with the highest possible resolution to allow the simulation of diffuse loads, heterogeneous characteristics of the basin both in terms of hydrological and hydraulic parameters on the use and occupation soil. The parallel work should improve the degree of understanding of the phenomena simulated in the basin as well as the activity of the calibration models, and this is supported by monitoring data acquired during the duration of the project MAPLU (Urban Stormwater Management) belonging to the network PROSAB (Research Program in Basic Sanitation) in the years 2006 to 2008

Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico

LINS, Filipe C. A. et al. Modelagem dinâmica e simulação computacional de poços de petróleo verticais e direcionais com elevação por bombeio mecânico. In: CONGRESSO BRASILEIRO DE PESQUISA E DESENVOLVIMENTO EM PETRÓLEO E GÁS, 5. 2009, Fortaleza, CE. Anais... Fortaleza: CBPDPetro, 2009.

Estudo da perda de carga e calor no poço injetor no processo De drenagem gravitacional assistida por vapor (SAGD)

The oil companies in the area in general are looking for new technologies that can increase the recovery factor of oil contained in reservoirs. These investments are mainly aimed at reducing the costs of projects which are high. Steam injection is one of these special methods of recovery in which steam is injected into the reservoir in order to reduce the viscosity of the oil and make it more mobile. The process assisted gravity drainage steam (SAGD) using steam injection in its mechanism, as well as two parallel horizontal wells. In this process steam is injected through the horizontal injection well, then a vapor chamber is formed by heating the oil in the reservoir and, by the action of gravitational forces, this oil is drained down to where the production well. This study aims to analyze the influence of pressure drop and heat along the injection well in the SAGD process. Numerical simulations were performed using the thermal simulator STARS of CMG (Computer Modeling Group). The parameters studied were the thermal conductivity of the formation, the flow of steam injection, the inner diameter of the column, the steam quality and temperature. A factorial design was used to verify the influence of the parameters studied in the recovery factor. We also analyzed different injection flow rates for the model with pressure drop and no pressure drop, as well as different maximum flow rates of oil production. Finally, we performed an economic analysis of the two models in order to check the profitability of the projects studied. The results showed that the pressure drop in injection well have a significant influence on the SAGD process.