Optimization and Control of Power Flow in Distribution Networks

Climate change is arguably the most critical issue facing our generation and the next. As we move towards a sustainable future, the grid is rapidly evolving with the integration of more and more renewable energy resources and the emergence of electric vehicles. In particular, large scale adoption of residential and commercial solar photovoltaics (PV) plants is completely changing the traditional slowly-varying unidirectional power flow nature of distribution systems. High share of intermittent renewables pose several technical challenges, including voltage and frequency control. But along with these challenges, renewable generators also bring with them millions of new DC-AC inverter controllers each year. These fast power electronic devices can provide an unprecedented opportunity to increase energy efficiency and improve power quality, if combined with well-designed inverter control algorithms. The main goal of this dissertation is to develop scalable power flow optimization and control methods that achieve system-wide efficiency, reliability, and robustness for power distribution networks of future with high penetration of distributed inverter-based renewable generators.

Proposed solutions to power flow control problems in the literature range from fully centralized to fully local ones. In this thesis, we will focus on the two ends of this spectrum. In the first half of this thesis (chapters 2 and 3), we seek optimal solutions to voltage control problems provided a centralized architecture with complete information. These solutions are particularly important for better understanding the overall system behavior and can serve as a benchmark to compare the performance of other control methods against. To this end, we first propose a branch flow model (BFM) for the analysis and optimization of radial and meshed networks. This model leads to a new approach to solve optimal power flow (OPF) problems using a two step relaxation procedure, which has proven to be both reliable and computationally efficient in dealing with the non-convexity of power flow equations in radial and weakly-meshed distribution networks. We will then apply the results to fast time- scale inverter var control problem and evaluate the performance on real-world circuits in Southern California Edison’s service territory.

The second half (chapters 4 and 5), however, is dedicated to study local control approaches, as they are the only options available for immediate implementation on today’s distribution networks that lack sufficient monitoring and communication infrastructure. In particular, we will follow a reverse and forward engineering approach to study the recently proposed piecewise linear volt/var control curves. It is the aim of this dissertation to tackle some key problems in these two areas and contribute by providing rigorous theoretical basis for future work.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

I. Configurational stability and redistribution equilibria in organomagnesium compounds. II. Redistribution equilibria in organocadmium compounds

I. CONFIGURATIONAL STABILITY AND REDISTRIBUTION EQUILIBRIA IN ORGANOMAGNESIUM COMPOUNDS

The dependence of the rate of inversion of a dialkylmagnesium compound on the solvent has been studied.

Examination of the temperature dependence of the nuclear magnetic resonance spectrum of 1-phenyl-2-propylmagnesium bromide in diethyl ether solution indicates that inversion of configuration at the methylene group of this Grignard reagent occurs with an approximate rate of 2 sec^-1 at room temperature. This is the first example of a rapid inversion rate in a secondary Grignard reagent.

The rates of exchange of alkyl groups between dineopentylmagnesium and di-s-butylmagnesium, bis-(2-methylbutyl)-magnesium and bis-(4, 4-dimethyl-2-pentyl)-magnesium respectively in diethyl ether solution were found to be fast on the nmr time scale. However, the alkyl group exchange rate was found to be slow in a diethyl ether solution of dineopentylmagnesium and bis-(2-methylbutyl)-magnesium containing N, N, N', N'-tetramethylethylenediamine. The unsymmetrical species neopentyl-2-methylbutyl-magnesium was observed at room temperature in the nmr spectrum of the solution containing the diamine.

II. REDISTRIBUTION EQUILIBRIA IN ORGANOCADMIUM COMPOUNDS

The exchange of methyl groups in dimethylcadmium has been studied by nuclear magnetic resonance spectroscopy. Activation parameters for the methyl group exchange have been measured for a neat sample and for a solution in tetrahydrofuran. The exchange is faster in the basic solvent tetrahydrofuran relative to the neat sample and in tetrahydrofuran solution is retarded by the solvating agent N, N, N’, N’-tetramethylethylenediamine and greatly increased by cadmium bromide. The addition of methanol to a solution of dimethylcadmium in tetrahydrofuran appears to have very little effect on the rate of exchange. The exchange was found to proceed with retention of configuration. The rate-limiting step for the exchange of methyl groups in a basic solvent appears to be the dissociation of coordinating solvent from dimethylcadmium.

The equilibrium between methylcadmium bromide, dimethylcadmium and cadmium bromide in tetrahydrofuran solution has also been studied. At room temperature the interconversion of the species is very fast on the nmr time scale but at -100° distinct absorptions for methylcadmium bromide and imethylcadmium are observed.

The species ethylmethylcadmium has been observed in the nmr spectrum.

The rate of exchange of vinyl groups in a solution of divinylcadmium in tetrahydrofuran has been found to be fast on the nmr time scale.