Two-level network design with intermediate facilities: An application to electrical distribution systems

We consider the two-level network design problem with intermediate facilities. This problem consists of designing a minimum cost network respecting some requirements, usually described in terms of the network topology or in terms of a desired flow of commodities between source and destination vertices. Each selected link must receive one of two types of edge facilities and the connection of different edge facilities requires a costly and capacitated vertex facility. We propose a hybrid decomposition approach which heuristically obtains tentative solutions for the vertex facilities number and location and use these solutions to limit the computational burden of a branch-and-cut algorithm. We test our method on instances of the power system secondary distribution network design problem. The results show that the method is efficient both in terms of solution quality and computational times. © 2010 Elsevier Ltd.

Mathematical decomposition technique applied to the probabilistic power flow problem

In this paper a framework based on the decomposition of the first-order optimality conditions is described and applied to solve the Probabilistic Power Flow (PPF) problem in a coordinated but decentralized way in the context of multi-area power systems. The purpose of the decomposition framework is to solve the problem through a process of solving smaller subproblems, associated with each area of the power system, iteratively. This strategy allows the probabilistic analysis of the variables of interest, in a particular area, without explicit knowledge of network data of the other interconnected areas, being only necessary to exchange border information related to the tie-lines between areas. An efficient method for probabilistic analysis, considering uncertainty in n system loads, is applied. The proposal is to use a particular case of the point estimate method, known as Two-Point Estimate Method (TPM), rather than the traditional approach based on Monte Carlo simulation. The main feature of the TPM is that it only requires resolve 2n power flows for to obtain the behavior of any random variable. An iterative coordination algorithm between areas is also presented. This algorithm solves the Multi-Area PPF problem in a decentralized way, ensures the independent operation of each area and integrates the decomposition framework and the TPM appropriately. The IEEE RTS-96 system is used in order to show the operation and effectiveness of the proposed approach and the Monte Carlo simulations are used to validation of the results. © 2011 IEEE.

Probabilistic analysis of the distributed power generation in weakly meshed distribution systems

This paper presents an approach for probabilistic analysis of unbalanced three-phase weakly meshed distribution systems considering uncertainty in load demand. In order to achieve high computational efficiency this approach uses both an efficient method for probabilistic analysis and a radial power flow. The probabilistic approach used is the well-known Two-Point Estimate Method. Meanwhile, the compensation-based radial power flow is used in order to extract benefits from the topological characteristics of the distribution systems. The generation model proposed allows modeling either PQ or PV bus on the connection point between the network and the distributed generator. In addition allows control of the generator operating conditions, such as the field current and the power delivery at terminals. Results on test with IEEE 37 bus system is given to illustrate the operation and effectiveness of the proposed approach. A Monte Carlo Simulations method is used to validate the results. © 2011 IEEE.

Distributed generation modelling for unbalanced three-phase power flow calculations in smart grids

Distributed Generators (DG) are generally modeled as PQ or PV buses in power flow studies. But in order to integrate DG units into the distribution systems and control the reactive power injection it is necessary to know the operation mode and the type of connection to the system. This paper presents a single-phase and a three-phase mathematical model to integrate DG in power flow calculations in distribution systems, especially suited for Smart Grid calculations. If the DG is in PV mode, each step of the power flow algorithm calculates the reactive power injection from the DG to the system to keep the voltage in the bus in a predefined level, if the DG is in PQ mode, the power injection is considered as a negative load. The method is tested on two well known test system, presenting single-phase results on 85 bus system, and three-phase results in the IEEE 34 bus test system. © 2011 IEEE.

Pricing of reactive power support provided by distributed generators in transmission systems

Distributed Generation, microgrid technologies, two-way communication systems, and demand response programs are issues that are being studied in recent years within the concept of smart grids. At some level of enough penetration, the Distributed Generators (DGs) can provide benefits for sub-transmission and transmission systems through the so-called ancillary services. This work is focused on the ancillary service of reactive power support provided by DGs, specifically Wind Turbine Generators (WTGs), with high level of impact on transmission systems. The main objective of this work is to propose an optimization methodology to price this service by determining the costs in which a DG incurs when it loses sales opportunity of active power, i.e, by determining the Loss of Opportunity Costs (LOC). LOC occur when more reactive power is required than available, and the active power generation has to be reduced in order to increase the reactive power capacity. In the optimization process, three objectives are considered: active power generation costs of DGs, voltage stability margin of the system, and losses in the lines of the network. Uncertainties of WTGs are reduced solving multi-objective optimal power flows in multiple probabilistic scenarios constructed by Monte Carlo simulations, and modeling the time series associated with the active power generation of each WTG via Fuzzy Logic and Markov Chains. The proposed methodology was tested using the IEEE 14 bus test system with two WTGs installed. © 2011 IEEE.

Harmonic propagation analysis in electric energy distribution systems

An important alteration of the equivalent loads profile has been observed in the electrical energy distribution systems, for the last years. Such fact is due to the significant increment of the electronic processors of electric energy that, in general, behave as nonlinear loads, generating harmonic distortions in the currents and voltages along the electric network. The effects of these nonlinear loads, even if they are concentrated in specific sections of the network, are present along the branch circuits, affecting the behavior of the entire electric network. For the evaluation of this phenomenon it is necessary the analysis of the harmonic currents flow and the understanding of the causes and effects of the consequent voltage harmonic distortions. The usual tools for calculation the harmonic flow consider one-line equivalent networks, balanced and symmetrical systems. Therefore, they are not tools appropriate for analysis of the operation and the influence/interaction of mitigation elements. In this context, this work proposes the development of a computational tool for the analysis of the three-phase harmonic propagation using Norton modified models and considering the real nature of unbalanced electric systems operation. © 2011 IEEE.

Simulation of the flow through automatic valves of hermetic compressors by the immersed boundary method approach

The main goal of the present work is to verify the applicability of the Immersed Boundary Method together with the Virtual Physical Model to solve the flow through automatic valves of hermetic compressors. The valve was simplified to a two-dimensional radial diffuser, with diameter ratio of D/d = 1.5, and simulated for a one cycle of opening and closing process with a imposed velocity of 3.0 cm/s for the reed, dimensionless gap between disks in the range of 0.07 < s/d < 0.10, and inlet Reynolds number equal to 1500. The good results obtained showed that the methodology has great potential as project tool for this type of valve systems. © The Authors, 2011.

Parameterization technique for the continuation power flow using the trivial and tangent predictor

This paper presents efficient geometric parameterization techniques using the tangent and the trivial predictors for the continuation power flow, developed from observation of the trajectories of the load flow solution. The parameterization technique eliminates the Jacobian matrix singularity of load flow, and therefore all the consequent problems of ill-conditioning, by the addition of the line equations which pass through the points in the plane determined by the variables loading factor and the real power generated by the slack bus, two parameters with clear physical meaning. This paper also provides an automatic step size control around the maximum loading point. Thus, the resulting method enables not only the calculation of the maximum loading point, but also the complete tracing of P-V curves of electric power systems. The technique combines robustness with ease of understanding. The results to the IEEE 300-bus system and of large real systems show the effectiveness of the proposed method. © 2012 IEEE.

Model to increase the lateral line length of drip irrigation systems

Non-pressure compensating drip hose is widely used for irrigation of vegetables and orchards. One limitation is that the lateral line length must be short to maintain uniformity due to head loss and slope. Any procedure to increase the length is appropriate because it represents low initial cost of the irrigation system. The hypothesis of this research is that it is possible to increase the lateral line length combining two points: using a larger spacing between emitters at the beginning of the lateral line and a smaller one after a certain distance; and allowing a higher pressure variation along the lateral line under an acceptable value of distribution uniformity. To evaluate this hypothesis, a nonlinear programming model (NLP) was developed. The input data are: diameter, roughness coefficient, pressure variation, emitter operational pressure, relationship between emitter discharge and pressure. The output data are: line length, discharge and length of the each section with different spacing between drippers, total discharge in the lateral line, multiple outlet adjustment coefficient, head losses, localized head loss, pressure variation, number of emitters, spacing between emitters, discharge in each emitter, and discharge per linear meter. The mathematical model developed was compared with the lateral line length obtained with the algebraic solution generated by the Darcy-Weisbach equation. The NLP model showed the best results since it generated the greater gain in the lateral line length, maintaining the uniformity and the flow variation under acceptable standards. It had also the lower flow variation, so its adoption is feasible and recommended.

A comparative study of aquifer systems occurring at the Paraná sedimentary basin, Brazil: U-isotopes contribution

The radioactivity due to 238U and 234U in three aquifer systems occurring within the Paraná sedimentary basin, South America, has been investigated. Uranium is much less dissolved from fractured igneous rocks than from the porous sedimentary rocks as indicated by the U-mobility coefficients between 7. 6 × 10-6 and 1. 2 × 10-3 g cm-3. These values are also compatible with the U preference ratios relative to Na, K, Ca, Mg and SiO2, which showed that U is never preferentially mobilized in the liquid phase during the flow occurring in cracks, fissures, fractures and faults of the igneous basaltic rocks. Experimental dissolution of diabase grains on a time-scale laboratory has demonstrated that the U dissolution appeared to be a two-stage process characterized by linear and second-order kinetics. The U dissolution rate was 8 × 10-16 mol m-2 s-1 that is within the range of 4 × 10-16-3 × 10-14 mol m-2 s-1 estimated for other rock types. The 234U/238U activity ratio of dissolved U in solutions was higher than unity, a typical result expected during the water-rock interactions when preferential 234U-leach from the rock surfaces takes place. Some U-isotopes data allowed estimating 320 ka for the groundwater residence time in a sector of a transect in São Paulo State. A modeling has been also realized considering all U-isotopes data obtained in Bauru (35 samples), Serra Geral (16 samples) and Guarani (29 samples) aquifers. The results indicated that the Bauru aquifer waters may result from the admixture of waters from Guarani (1. 5 %) and Serra Geral (98. 5 %) aquifers. © 2012 Springer-Verlag.

Improving the grid operation and reliability cost of distribution systems with dispersed generation

In this work, a mathematical model to analyze the impact of the installation and operation of dispersed generation units in power distribution systems is proposed. The main focus is to determine the trade-off between the reliability and operational costs of distribution networks when the operation of isolated areas is allowed. In order to increase the system operator revenue, an optimal power flow makes use of the different energy prices offered by the dispersed generation connected to the grid. Simultaneously, the type and location of the protective devices initially installed on the protection system are reconfigured in order to minimize the interruption and expenditure of adjusting the protection system to conditions imposed by the operation of dispersed units. The interruption cost regards the unsupplied energy to customers in secure systems but affected by the normal tripping of protective devices. Therefore, the tripping of fuses, reclosers, and overcurrent relays aims to protect the system against both temporary and permanent fault types. Additionally, in order to reduce the average duration of the system interruption experienced by customers, the isolated operation of dispersed generation is allowed by installing directional overcurrent relays with synchronized reclose capabilities. A 135-bus real distribution system is used in order to show the advantages of using the mathematical model proposed. © 1969-2012 IEEE.

Transdentinal cytotoxicity of experimental adhesive systems of different hydrophilicity applied to ethanol-saturated dentin

The aim of this study was to evaluate the transdentinal cytotoxicity of experimental adhesive systems (EASs) with different hydrophilicity and dentin saturation solutions on odontoblast-like cells. One hundred 0.4-mm-thick dentin discs were mounted in in vitro pulp chambers and assigned to 10 groups. MDPC-23 cells were seeded onto the pulpal side of the discs, incubated for 48 h. The EASs with increasing hydrophilicity (R1, R2, R3 and R4) were applied to the occlusal side after etching and saturation of etched dentin with water or ethanol. R0 (no adhesive) served as controls. R1 is a non-solvated hydrophobic blend, R2 is similar to a simplified etch-and-rinse adhesive system and R3 and R4 are similar to self-etching adhesives. After 24 h, cell metabolism was evaluated by MTT assay (n = 8 discs) and cell morphology was examined by SEM (n = 2 discs). Type of cell death was identified by flow cytometry and the degree of monomer conversion (%DC) was determined by infrared spectroscopy (FTIR) after 10 s or 20 s of photoactivation. Data were analyzed by the Kruskal-Wallis and Mann-Whitney tests (α = 0.05). Dentin saturation with ethanol resulted in higher necrotic cell death ratios for R2, R3 and R4 compared with water saturation, although R2 and R3 induced higher SDH production. Photoactivation for 20 s significantly improved the %DC of all EASs compared with 10 s. A significant positive correlation was observed between the degree of hydrophilicity and %DC. In conclusion, except for R1, dentin saturation with ethanol increased the cytotoxicity of EASs, as expressed by the induction of necrotic cell death. © 2013 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Time-resolved random laser spectroscopy of inhomogeneously broadened systems

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

A two-fluid model for refrigerant flow through adiabatic capillary tubes

This work presents a numerical model to simulate refrigerant flow through capillary tubes, commonly used as expansion devices in refrigeration systems. The flow is divided in a single-phase region, where the refrigerant is in the subcooled liquid state, and a region of two-phase flow. The capillary tube is considered straight and horizontal. The flow is taken as one-dimensional and adiabatic. Steady-state condition is also assumed and the metastable flow phenomena are neglected. The two-fluid model, considering the hydrodynamic and thermal non-equilibrium between the liquid and vapor phases, is applied to the two-phase flow region. Comparisons are made with experimental measurements of the mass flow rate and pressure distribution along two capillary tubes working with refrigerant R-134a in different operating conditions. The results indicate that the present model provides a better estimation than the commonly employed homogeneous model. Some computational results referring to the quality, void fraction, velocities, and temperatures of each phase are presented and discussed.

Monitoring of diclofenac with biomimetic sensor in batch and FIA systems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)