Time-domain models for power system stability and unbalance

It is an important and difficult challenge to protect modern interconnected power system from blackouts. Applying advanced power system protection techniques and increasing power system stability are ways to improve the reliability and security of power systems. Phasor-domain software packages such as Power System Simulator for Engineers (PSS/E) can be used to study large power systems but cannot be used for transient analysis. In order to observe both power system stability and transient behavior of the system during disturbances, modeling has to be done in the time-domain. This work focuses on modeling of power systems and various control systems in the Alternative Transients Program (ATP). ATP is a time-domain power system modeling software in which all the power system components can be modeled in detail. Models are implemented with attention to component representation and parameters. The synchronous machine model includes the saturation characteristics and control interface. Transient Analysis Control System is used to model the excitation control system, power system stabilizer and the turbine governor system of the synchronous machine. Several base cases of a single machine system are modeled and benchmarked against PSS/E. A two area system is modeled and inter-area and intra-area oscillations are observed. The two area system is reduced to a two machine system using reduced dynamic equivalencing. The original and the reduced systems are benchmarked against PSS/E. This work also includes the simulation of single-pole tripping using one of the base case models. Advantages of single-pole tripping and comparison of system behavior against three-pole tripping are studied. Results indicate that the built-in control system models in PSS/E can be effectively reproduced in ATP. The benchmarked models correctly simulate the power system dynamics. The successful implementation of a dynamically reduced system in ATP shows promise for studying a small sub-system of a large system without losing the dynamic behaviors. Other aspects such as relaying can be investigated using the benchmarked models. It is expected that this work will provide guidance in modeling different control systems for the synchronous machine and in representing dynamic equivalents of large power systems.

Methodology for Dynamic Stability and Robustness Analysis of Commercial-Power-Module-Based DC-Distributed Systems

El propósito de esta tesis es presentar una metodología para realizar análisis de la dinámica en pequeña señal y el comportamiento de sistemas de alimentación distribuidos de corriente continua (CC), formados por módulos comerciales. Para ello se hace uso de un método sencillo que indica los márgenes de estabilidad menos conservadores posibles mediante un solo número. Este índice es calculado en cada una de las interfaces que componen el sistema y puede usarse para obtener un índice global que indica la estabilidad del sistema global. De esta manera se posibilita la comparación de sistemas de alimentación distribuidos en términos de robustez. La interconexión de convertidores CC-CC entre ellos y con los filtros EMI necesarios puede originar interacciones no deseadas que dan lugar a la degradación del comportamiento de los convertidores, haciendo el sistema más propenso a inestabilidades. Esta diferencia en el comportamiento se debe a interacciones entre las impedancias de los diversos elementos del sistema. En la mayoría de los casos, los sistemas de alimentación distribuida están formados por módulos comerciales cuya estructura interna es desconocida. Por ello los análisis presentados en esta tesis se basan en medidas de la respuesta en frecuencia del convertidor que pueden realizarse desde los terminales de entrada y salida del mismo. Utilizando las medidas de las impedancias de entrada y salida de los elementos del sistema, se puede construir una función de sensibilidad que proporciona los márgenes de estabilidad de las diferentes interfaces. En esta tesis se utiliza el concepto del valor máximo de la función de sensibilidad (MPC por sus siglas en inglés) para indicar los márgenes de estabilidad como un único número. Una vez que la estabilidad de todas las interfaces del sistema se han evaluado individualmente, los índices obtenidos pueden combinarse para obtener un único número con el que comparar la estabilidad de diferentes sistemas. Igualmente se han analizado las posibles interacciones en la entrada y la salida de los convertidores CC-CC, obteniéndose expresiones analíticas con las que describir en detalle los acoplamientos generados en el sistema. Los estudios analíticos realizados se han validado experimentalmente a lo largo de la tesis. El análisis presentado en esta tesis se culmina con la obtención de un índice que condensa los márgenes de estabilidad menos conservativos. También se demuestra que la robustez del sistema está asegurada si las impedancias utilizadas en la función de sensibilidad se obtienen justamente en la entrada o la salida del subsistema que está siendo analizado. Por otra parte, la tesis presenta un conjunto de parámetros internos asimilados a impedancias, junto con sus expresiones analíticas, que permiten una explicación detallada de las interacciones en el sistema. Dichas expresiones analíticas pueden obtenerse bien mediante las funciones de transferencia analíticas si se conoce la estructura interna, o utilizando medidas en frecuencia o identificación de las mismas a través de la respuesta temporal del convertidor. De acuerdo a las metodologías presentadas en esta tesis se puede predecir la estabilidad y el comportamiento de sistemas compuestos básicamente por convertidores CC-CC y filtros, cuya estructura interna es desconocida. La predicción se basa en un índice que condensa la información de los márgenes de estabilidad y que permite la obtención de un indicador de la estabilidad global de todo el sistema, permitiendo la comparación de la estabilidad de diferentes arquitecturas de sistemas de alimentación distribuidos. ABSTRACT The purpose of this thesis is to present dynamic small-signal stability and performance analysis methodology for dc-distributed systems consisting of commercial power modules. Furthermore, the objective is to introduce simple method to state the least conservative margins for robust stability as a single number. In addition, an index characterizing the overall system stability is obtained, based on which different dc-distributed systems can be compared in terms of robustness. The interconnected systems are prone to impedance-based interactions which might lead to transient-performance degradation or even instability. These systems typically are constructed using commercial converters with unknown internal structure. Therefore, the analysis presented throughout this thesis is based on frequency responses measurable from the input and output terminals. The stability margins are stated utilizing a concept of maximum peak criteria, derived from the behavior of impedance-based sensitivity function that provides a single number to state robust stability. Using this concept, the stability information at every system interface is combined to a meaningful number to state the average robustness of the system. In addition, theoretical formulas are extracted to assess source and load side interactions in order to describe detailed couplings within the system. The presented theoretical analysis methodologies are experimentally validated throughout the thesis. In this thesis, according to the presented analysis, the least conservative stability margins are provided as a single number guaranteeing robustness. It is also shown that within the interconnected system the robust stability is ensured only if the impedance-based minor-loop gain is determined at the very input or output of each subsystem. Moreover, a complete set of impedance-type internal parameters as well as the formulas according to which the interaction sensitivity can be fully explained and analyzed, is provided. The given formulation can be utilized equally either based on measured frequency responses, time-domain identified internal parameters or extracted analytic transfer functions. Based on the analysis methodologies presented in this thesis, the stability and performance of interconnected systems consisting of converters with unknown internal structure, can be predicted. Moreover, the provided concept to assess the least conservative stability margins enables to obtain an index to state the overall robust stability of distributed power architecture and thus to compare different systems in terms of stability.

Astrophysical parameters and orbital solution of the peculiar X-ray transient IGR J00370+6122

Context. BD + 60° 73 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. Aims. We aim to characterise the binary system IGR J00370+6122 by deriving its orbital and physical parameters. Methods. We obtained high-resolution spectra of BD + 60° 73 at different epochs. We used the fastwind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. Results. BD + 60° 73 is a BN0.7 Ib low-luminosity supergiant located at a distance ~3.1 kpc, in the Cas OB4 association. We derive Teff = 24 000 K and log gc = 3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1-σ level with a mass M∗ ≈ 15 M⊙. The recurrence time of the X-ray flares is the orbital period of the system. The neutron star is in a high-eccentricity (e = 0.56 ± 0.07) orbit, and the X-ray emission is strongly peaked around orbital phase φ = 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. Conclusions. The X-ray behaviour of IGR J00370+6122 is reminiscent of “intermediate” supergiant X-ray transients, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which when combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will very likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Short-term stability in refractive status despite large fluctuations in glucose levels in diabetes mellitus type 1 and 2

Purpose: This work investigates how short-term changes in blood glucose concentration affect the refractive components of the diabetic eye in patients with long-term Type 1 and Type 2 diabetes. Methods: Blood glucose concentration, refractive error components (mean spherical equivalent MSE, J0, J45), central corneal thickness (CCT), anterior chamber depth (ACD), crystalline lens thickness (LT), axial length (AL) and ocular aberrations were monitored at two-hourly intervals over a 12-hour period in: 20 T1DM patients (mean age ± SD) 38±14 years, baseline HbA1c 8.6±1.9%; 21 T2DM patients (mean age ± SD) 56±11 years, HbA1c 7.5±1.8%; and in 20 control subjects (mean age ± SD) 49±23 years, HbA1c 5.5±0.5%. The refractive and biometric results were compared with the corresponding changes in blood glucose concentration. Results: Blood glucose concentration at different times was found to vary significantly within (p<0.0005) and between groups (p<0.0005). However, the refractive error components and ocular aberrations were not found to alter significantly over the day in either the diabetic patients or the control subjects (p>0.05). Minor changes of marginal statistical or optical significance were observed in some biometric parameters. Similarly there were some marginally significant differences between the baseline biometric parameters of well-controlled and poorly-controlled diabetic subjects. Conclusion: This work suggests that normal, short-term fluctuations (of up to about 6 mM/l on a timescale of a few hours) in the blood glucose levels of diabetics are not usually associated with acute changes in refractive error or ocular wavefront aberrations. It is therefore possible that factors other than refractive error fluctuations are sometimes responsible for the transient visual problems often reported by diabetic patients. © 2012 Huntjens et al.

Purity evaluation and influence of carbon nanotube on carbon nanotube/graphite thermal stability

High resolution thermogravimetry has been used to evaluate the carbonaceous content in a commercial sample of single-walled carbon nanotube (SWNT). The content of SWNTs in the sample was found to be at least 77mass% which was supported by images obtained with scanning and transmission electron microscopies (SEM and TEM). Furthermore, the influence of SWNT addition on the thermal stability of graphite in mixtures of SWNT/graphite at different proportions was investigated. The graphite stability decreased with the increased of SWNT content in the overall range of composition. This behavior could be due to the close contact between these carbonaceous species as determined by SEM analysis.

Development of diagnostic and prognostic algorithms for SF6 puffer circuit breakers from transient waveforms : a validation proposal

- This paper presents a validation proposal for development of diagnostic and prognostic algorithms for SF6 puffer circuit-breakers reproduced from actual site waveforms. The re-ignition/restriking rates are duplicated in given circuits and the cumulative energy dissipated in interrupters by the restriking currents. The targeted objective is to provide a simulated database for diagnosis of re-ignition/restrikes relating to the phase to earth voltage and the number of re-ignition/restrikes as well as estimating the remaining life of SF6 circuit-breakers. The model-based diagnosis of a tool will be useful in monitoring re-ignition/restrikes as well as predicting a nozzle’s lifetime. This will help ATP users with practical study cases and component data compilation for shunt reactor switching and capacitor switching. This method can be easily applied with different data for the different dielectric curves of circuit breakers and networks. This paper presents modelling details and some of the available cases, required project support, the validation proposal, the specific plan for implementation and the propsed main contributions.