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

Artificial neural networks for load flow and external equivalents studies

In this paper an artificial neural network (ANN) based methodology is proposed for (a) solving the basic load flow, (b) solving the load flow considering the reactive power limits of generation (PV) buses, (c) determining a good quality load flow starting point for ill-conditioned systems, and (d) computing static external equivalent circuits. An analysis of the input data required as well as the ANN architecture is presented. A multilayer perceptron trained with the Levenberg-Marquardt second order method is used. The proposed methodology was tested with the IEEE 30- and 57-bus, and an ill-conditioned 11-bus system. Normal operating conditions (base case) and several contingency situations including different load and generation scenarios have been considered. Simulation results show the excellent performance of the ANN for solving problems (a)-(d). (C) 2010 Elsevier B.V. All rights reserved.

Correlation between the intensity of venous reflux in the saphenofemoral junction and morphological changes of the great saphenous vein by duplex scanning in patients with primary varicosis

Aim. One of the major causes of chronic venous disease is venous reflux, the identification and quantification of which are important for diagnosis. Duplex scanning allows for the detection and quantification of reflux in individual veins. Evaluation of the great saphenous vein in primary varicosis is necessary for its preservation. Objective of the study is to evaluate a possible correlation between the intensity of reflux at the saphenofemoral junction, diameter alterations of the incompetent great saphenous vein and the practical effect of such correlation. Also to compare the clinical severity of the CEAP classification with such parameters.Methods. Three hundred limbs were submitted to duplex evaluation of their insufficient saphenous veins. Vein diameter was measured on five different points. Velocity and flow at reflux peak and reflux time were determined. The saphenous vein's diameters were correlated with velocity, flow and time. The three latter parameters and diameters were compared with clinical severity according to CEAP.Results. Correlation was found between the saphenous vein's diameters, velocity and flow. No correlation was observed between time and diameter in the thigh's upper and middle thirds. When comparing diameter, velocity and flow with CEAP clinical severity classification, an association was observed. The correlation between reflux time with clinical severity was weak.Conclusion. Reflux time is a good parameter for identifying the presence of reflux, but not for quantifying it. Velocity and peak flow were better parameters for evaluating reflux intensity as they were correlated with great saphenous vein alterations, and were associated with the disease's clinical severity. [Int Angiol 2010;29:323-30]

Measurement of the t(t)over-bar production cross section in p(p)over-bar collisions at root s=1.96 TeV using lepton plus jets events with lifetime b-tagging

We present a measurement of the top quark pair (t (t) over bar) production cross section (sigma(t (t) over bar)) in pp collisions at root s = 1.96 TeV using 230 pb(-1) of data collected by the DO experiment at the Fermilab Tevatron Collider. We select events with one charged lepton (electron or muon), missing transverse energy, and jets in the final state. We employ lifetime-based b-jet identification techniques to further enhance the t F purity of the selected sample. For a top quark mass of 175 GeV, we measure sigma(t (t) over bar) 8.6(-1.5)(+1.6) (stat. + syst.) +/- 0.6(lumi.) pb, in agreement with the standard model expectation. (c) 2005 Published by Elsevier B.V.

On negative flows of the AKNS hierarchy and a class of deformations of a bihamiltonian structure of hydrodynamic type

A deformation parameter of a bihamiltonian structure of hydrodynamic type is shown to parametrize different extensions of the AKNS hierarchy to include negative flows. This construction establishes a purely algebraic link between, on the one hand, two realizations of the first negative flow of the AKNS model and, on the other, two-component generalizations of Camassa-Holmand Dym-type equations. The two-component generalizations of Camassa-Holm- and Dym-type equations can be obtained from the negative-order Hamiltonians constructed from the Lenard relations recursively applied on the Casimir of the first Poisson bracket of hydrodynamic type. The positive-order Hamiltonians, which follow froth the Lenard scheme applied on the Casimir of the second Poisson bracket of hydrodynamic type, are shown to coincide with the Hamiltonians of the AKNS model. The AKNS Hamiltonians give rise to charges conserved with respect to equations of motion of two-component Camassa-Holm- and two-component Dym-type equations.

Measurement of the CP-violation parameter of B-0 mixing and decay with p(p)over-bar ->mu mu X data

We measure the dimuon charge asymmetry A in p (p) over bar collisions at a center of mass energy root s=1960 GeV. The data was recorded with the D0 detector and corresponds to an integrated luminosity of approximately 1.0 fb(-1). Assuming that the asymmetry A is due to asymmetric B-0 <->(B) over bar (0) mixing and decay, we extract the CP-violation parameter of B-0 mixing and decay: ((epsilon B0))/(1+vertical bar epsilon B0 vertical bar 2)=(AB0)/(4)= -0.0023 +/- 0.0011(stat)+/- 0.0008(syst).A(B)(0) is the dimuon charge asymmetry from decays of B-0(B) over bar (0) pairs. The general case, with CP violation in both B-0 and B-s(0) systems, is also considered. Finally we obtain the forward-backward asymmetry that quantifies the tendency of mu(+) to go in the proton direction and mu(-) to go in the antiproton direction. The results are consistent with the standard model and constrain new physics.

Search for supersymmetry in pp collisions at 7 TeV in events with jets and missing transverse energy

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

Simultaneous measurement of the ratio R = B(t -> Wb)/B(t -> Wq) and the top-quark pair production cross section with the D0 detector at root s=1.96 TeV

We present the first simultaneous measurement of the ratio of branching fractions, R = B(t --> Wb)/B(t --> Wq), with q being a d, s, or b quark, and the top-quark pair production cross section sigma(t (t) over bar) in the lepton plus jets channel using 0.9 fb(-1) of p (p) over bar collision data at root s = 1.96 TeV collected with the D0 detector. We extract R and sigma(t (t) over bar) by analyzing samples of events with 0, 1, and >= 2 identified b jets. We measure R = 0.97(-0.08)(+0.09) (stat + syst) and sigma(t (t) over bar) = 8.18(-0.84)(+0.90) (stat + syst) +/- 0.50(lumi) pb, in agreement with the standard model prediction.

Model-independent measurement of the W-boson helicity in top-quark decays at D0

We present the first model-independent measurement of the helicity of W bosons produced in top quark decays, based on a 1 fb(-1) sample of candidate t (t) over bar events in the dilepton and lepton plus jets channels collected by the D0 detector at the Fermilab Tevatron p (p) over bar Collider. We reconstruct the angle theta(*) between the momenta of the down-type fermion and the top quark in the W boson rest frame for each top quark decay. A fit of the resulting cos theta(*) distribution finds that the fraction of longitudinal W bosons f(0)=0.425 +/- 0.166(stat)+/- 0.102(syst) and the fraction of right-handed W bosons f(+)=0.119 +/- 0.090(stat)+/- 0.053(syst), which is consistent at the 30% C.L. with the standard model.

Evidence of WW and WZ Production with lepton plus jets Final States in pp Collisions at root s=1.96 TeV

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

Study of the (D)over-barN Interaction in a QCD Coulomb Gauge Quark Model

We study the (D) over barN interaction at low energies with a quark model inspired in the QCD Hamiltonian in Coulomb gauge. The model Hamiltonian incorporates a confining Coulomb potential extracted from a self-consistent quasiparticle method for the gluon degrees of freedom, and transverse-gluon hyperfine interaction consistent with a finite gluon propagator in the infrared. Initially a constituent-quark mass function is obtained by solving a gap equation and baryon and meson bound-states are obtained in Fock space using a variational calculation. Next, having obtained the constituent-quark masses and the hadron waves functions, an effective meson-nucleon interaction is derived from a quark-interchange mechanism. This leads to a short range meson-baryon interaction and to describe long-distance physics vector- and scalar-meson exchanges described by effective Lagrangians are incorporated. The derived effective (D) over barN potential is used in a Lippmann-Schwinger equation to obtain phase shifts. The results are compared with a recent similar calculation using the nonrelativistic quark model.

Search for a Vectorlike Quark with Charge 2/3 in t plus Z Events from pp Collisions at root s 7 TeV

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