Decentralized AC power flow for real-time multi-TSO power system operation

This paper adjusts decentralized OPF optimization to the AC power flow problem in power systems with interconnected areas operated by diferent transmission system operators (TSO). The proposed methodology allows finding the operation point of 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. The methodology is based on the decomposition of the first-order optimality conditions of the AC power flow, which is formulated as a nonlinear programming problem. To allow better visualization of the concept of independent operation of each TSO, an artificial neural network have been used for computing border information of the interconnected TSOs. A multi-area Power Flow tool can be seen as a basic building block able to address a large number of problems under a multi-TSO competitive market philosophy. The IEEE RTS-96 power system is used in order to show the operation and effectiveness of the decentralized AC Power Flow. ©2010 IEEE.

Obtenção de materiais magnéticos para o estudo da propriedade magnetocalórica

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

Photoelectrocatalytic/photoelectro-Fenton coupling system using a nanostructured photoanode for the oxidation of a textile dye: kinetics study and oxidation pathway

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

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth's magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.

Higher order ionospheric effects in GNSS positioning in the European region

After removal of the Selective Availability in 2000, the ionosphere became the dominant error source for Global Navigation Satellite Systems (GNSS), especially for the high-accuracy (cm-mm) demanding applications like the Precise Point Positioning (PPP) and Real Time Kinematic (RTK) positioning.The common practice of eliminating the ionospheric error, e. g. by the ionosphere free (IF) observable, which is a linear combination of observables on two frequencies such as GPS L1 and L2, accounts for about 99% of the total ionospheric effect, known as the first order ionospheric effect (Ion1). The remaining 1% residual range errors (RREs) in the IF observable are due to the higher - second and third, order ionospheric effects, Ion2 and Ion3, respectively. Both terms are related with the electron content along the signal path; moreover Ion2 term is associated with the influence of the geomagnetic field on the ionospheric refractive index and Ion3 with the ray bending effect of the ionosphere, which can cause significant deviation in the ray trajectory (due to strong electron density gradients in the ionosphere) such that the error contribution of Ion3 can exceed that of Ion2 (Kim and Tinin, 2007).The higher order error terms do not cancel out in the (first order) ionospherically corrected observable and as such, when not accounted for, they can degrade the accuracy of GNSS positioning, depending on the level of the solar activity and geomagnetic and ionospheric conditions (Hoque and Jakowski, 2007). Simulation results from early 1990s show that Ion2 and Ion3 would contribute to the ionospheric error budget by less than 1% of the Ion1 term at GPS frequencies (Datta-Barua et al., 2008). Although the IF observable may provide sufficient accuracy for most GNSS applications, Ion2 and Ion3 need to be considered for higher accuracy demanding applications especially at times of higher solar activity.This paper investigates the higher order ionospheric effects (Ion2 and Ion3, however excluding the ray bending effects associated with Ion3) in the European region in the GNSS positioning considering the precise point positioning (PPP) method. For this purpose observations from four European stations were considered. These observations were taken in four time intervals corresponding to various geophysical conditions: the active and quiet periods of the solar cycle, 2001 and 2006, respectively, excluding the effects of disturbances in the geomagnetic field (i.e. geomagnetic storms), as well as the years of 2001 and 2003, this time including the impact of geomagnetic disturbances. The program RINEX_HO (Marques et al., 2011) was used to calculate the magnitudes of Ion2 and Ion3 on the range measurements as well as the total electron content (TEC) observed on each receiver-satellite link. The program also corrects the GPS observation files for Ion2 and Ion3; thereafter it is possible to perform PPP with both the original and corrected GPS observation files to analyze the impact of the higher order ionospheric error terms excluding the ray bending effect which may become significant especially at low elevation angles (Ioannides and Strangeways, 2002) on the estimated station coordinates.

Efeitos de segunda e terceira ordem da ionosfera no posicionamento GNSS no Brasil

O Sistema de Posicionamento Global (GPS) transmite seus sinais em duas freqüências, o que permite eliminar matematicamente os efeitos de primeira ordem da ionosfera através da combinação linear ionosphere free. Porém, restam os efeitos de segunda e terceira ordem, os quais podem provocar erros da ordem de centímetros nas medidas GPS. Esses efeitos, geralmente, são negligenciados no processamento dos dados GPS. Os efeitos ionosféricos de primeira, segunda e terceira ordem são diretamente proporcionais ao TEC presente na ionosfera, porém, no caso dos efeitos de segunda e terceira ordem, comparecem também o campo magnético da Terra e a máxima densidade de elétrons, respectivamente. Nesse artigo, os efeitos de segunda e terceira ordem da ionosfera são investigados, sendo que foram levados em consideração no processamento de dados GPS na região brasileira para fins de posicionamento. Serão apresentados os modelos matemáticos associados a esses efeitos, as transformações envolvendo o campo magnético da Terra e a utilização do TEC advindo dos Mapas Globais da Ionosfera ou calculados a partir das observações GPS de pseudodistância. O processamento dos dados GPS foi realizado considerando o método relativo estático e cinemático e o posicionamento por ponto preciso (PPP). Os efeitos de segunda e terceira ordem foram analisados considerando períodos de alta e baixa atividade ionosférica. Os resultados mostraram que a não consideração desses efeitos no posicionamento por ponto preciso e no posicionamento relativo para linhas de base longas pode introduzir variações da ordem de poucos milímetros nas coordenadas das estações, além de variações diurnas em altitude da ordem de centímetros.

Resonant excitation of Mn local vibrational modes in the higher order Raman spectra of nanocrystalline Ga(1-x)Mn(x)N films

The effect of manganese on the vibrational properties of Ga(1-x)Mn(x)N (0 <= x <= 0.18) films has been investigated by Raman scattering using 488.0 and 632.8 nm photon excitations. The first-order transverse and longitudinal optical GaN vibrational bands were observed in the whole composition range using both excitations, while the corresponding overtones, as well as a prominent peak located in 1238 cm(-1) (153.5 meV) were only observed in the Mn-containing films under 488.0 nm excitation. We propose that the peak observed at 1238 cm(-1) is due to resonant Mn local vibrational modes, the excitation process being related to electronic transitions involving the Mn acceptor band.

Deposition and characterization of BiVO4 thin films and evaluation as photoanodes for methylene blue degradation

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

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.

Um método primal-dual aplicado na resolução do problema de fluxo de potência ótimo

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

Fluxo de potência ótimo reativo via método da função lagrangiana barreira modificada

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

Loss minimization by the predictor-corrector modified barrier approach

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

LIMIT CYCLES of DISCONTINUOUS PIECEWISE LINEAR DIFFERENTIAL SYSTEMS

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

Optimal attitude corrections for cylindrical spacecraft

A first order analytical model for optimal small amplitude attitude maneuvers of spacecraft with cylindrical symmetry in an elliptical orbits is presented. The optimization problem is formulated as a Mayer problem with the control torques provided by a power limited propulsion system. The state is defined by Seffet-Andoyer's variables and the control by the components of the propulsive torques. The Pontryagin Maximum Principle is applied to the problem and the optimal torques are given explicitly in Serret-Andoyer's variables and their adjoints. For small amplitude attitude maneuvers, the optimal Hamiltonian function is linearized around a reference attitude. A complete first order analytical solution is obtained by simple quadrature and is expressed through a linear algebraic system involving the initial values of the adjoint variables. A numerical solution is obtained by taking the Euler angles formulation of the problem, solving the two-point boundary problem through the shooting method, and, then, determining the Serret-Andoyer variables through Serret-Andoyer transformation. Numerical results show that the first order solution provides a good approximation to the optimal control law and also that is possible to establish an optimal control law for the artificial satellite's attitude. (C) 2003 COSPAR. Published by Elsevier B.V. Ltd. All rights reserved.

Robust broadband vibration control of a flexible structure using an electrical dynamic absorber

This paper presents a simple but practical feedback control method to suppress the vibration of a flexible structure in the frequency range between 10 Hz and 1 kHz. A dynamic vibration absorber is designed for this, which has a natural frequency of 100 Hz and a normalized bandwidth (twice the damping ratio) of 9.9. The absorber is realized electrically by feeding back the structural acceleration at one position on the host structure to a collocated piezoceramic patch actuator via an analog controller consisting of a second-order lowpass filter. This absorber is equivalent to a single degree-of-freedom mechanical oscillator consisting of a serially connected mass-spring-damper system. A first-order lowpass filter is additionally used to improve stability at very high frequencies. Experiments were conducted on a free-free beam embedded with a piezoceramic patch actuator and an accelerometer at its center. It is demonstrated that the single absorber can simultaneously suppress multiple vibration modes within the control bandwidth. It is further shown that the control system is robust to slight changes in the plant. The method described can be applied to many other practical structures, after retuning the absorber parameters for the structure under control.