Distributed Load Control in Multiphase Radial Networks

The current power grid is on the cusp of modernization due to the emergence of distributed generation and controllable loads, as well as renewable energy. On one hand, distributed and renewable generation is volatile and difficult to dispatch. On the other hand, controllable loads provide significant potential for compensating for the uncertainties. In a future grid where there are thousands or millions of controllable loads and a large portion of the generation comes from volatile sources like wind and solar, distributed control that shifts or reduces the power consumption of electric loads in a reliable and economic way would be highly valuable.

Load control needs to be conducted with network awareness. Otherwise, voltage violations and overloading of circuit devices are likely. To model these effects, network power flows and voltages have to be considered explicitly. However, the physical laws that determine power flows and voltages are nonlinear. Furthermore, while distributed generation and controllable loads are mostly located in distribution networks that are multiphase and radial, most of the power flow studies focus on single-phase networks.

This thesis focuses on distributed load control in multiphase radial distribution networks. In particular, we first study distributed load control without considering network constraints, and then consider network-aware distributed load control.

Distributed implementation of load control is the main challenge if network constraints can be ignored. In this case, we first ignore the uncertainties in renewable generation and load arrivals, and propose a distributed load control algorithm, Algorithm 1, that optimally schedules the deferrable loads to shape the net electricity demand. Deferrable loads refer to loads whose total energy consumption is fixed, but energy usage can be shifted over time in response to network conditions. Algorithm 1 is a distributed gradient decent algorithm, and empirically converges to optimal deferrable load schedules within 15 iterations.

We then extend Algorithm 1 to a real-time setup where deferrable loads arrive over time, and only imprecise predictions about future renewable generation and load are available at the time of decision making. The real-time algorithm Algorithm 2 is based on model-predictive control: Algorithm 2 uses updated predictions on renewable generation as the true values, and computes a pseudo load to simulate future deferrable load. The pseudo load consumes 0 power at the current time step, and its total energy consumption equals the expectation of future deferrable load total energy request.

Network constraints, e.g., transformer loading constraints and voltage regulation constraints, bring significant challenge to the load control problem since power flows and voltages are governed by nonlinear physical laws. Remarkably, distribution networks are usually multiphase and radial. Two approaches are explored to overcome this challenge: one based on convex relaxation and the other that seeks a locally optimal load schedule.

To explore the convex relaxation approach, a novel but equivalent power flow model, the branch flow model, is developed, and a semidefinite programming relaxation, called BFM-SDP, is obtained using the branch flow model. BFM-SDP is mathematically equivalent to a standard convex relaxation proposed in the literature, but numerically is much more stable. Empirical studies show that BFM-SDP is numerically exact for the IEEE 13-, 34-, 37-, 123-bus networks and a real-world 2065-bus network, while the standard convex relaxation is numerically exact for only two of these networks.

Theoretical guarantees on the exactness of convex relaxations are provided for two types of networks: single-phase radial alternative-current (AC) networks, and single-phase mesh direct-current (DC) networks. In particular, for single-phase radial AC networks, we prove that a second-order cone program (SOCP) relaxation is exact if voltage upper bounds are not binding; we also modify the optimal load control problem so that its SOCP relaxation is always exact. For single-phase mesh DC networks, we prove that an SOCP relaxation is exact if 1) voltage upper bounds are not binding, or 2) voltage upper bounds are uniform and power injection lower bounds are strictly negative; we also modify the optimal load control problem so that its SOCP relaxation is always exact.

To seek a locally optimal load schedule, a distributed gradient-decent algorithm, Algorithm 9, is proposed. The suboptimality gap of the algorithm is rigorously characterized and close to 0 for practical networks. Furthermore, unlike the convex relaxation approach, Algorithm 9 ensures a feasible solution. The gradients used in Algorithm 9 are estimated based on a linear approximation of the power flow, which is derived with the following assumptions: 1) line losses are negligible; and 2) voltages are reasonably balanced. Both assumptions are satisfied in practical distribution networks. Empirical results show that Algorithm 9 obtains 70+ times speed up over the convex relaxation approach, at the cost of a suboptimality within numerical precision.

Detection of Brillouin scattering temperature signal in Brillouin optical time-domain reflectometer sensing system based on instantaneous frequency measurement technology

We present a simple and practical method for the single-ended distributed fiber temperature measurements using microwave (11-GHz) coherent detection and the instantaneous frequency measurement (IFM) technique to detect spontaneous Brillouin backscattered signal in which a specially designed rf bandpass filter at 11 GHz is used as a frequency discriminator to transform frequency shift to intensity fluctuation. A Brillouin temperature signal can be obtained at 11 GHz over a sensing length of 10 km. The power sensitivity dependence on temperature induced by frequency shift is measured as 2.66%/K. (c) 2007 Society of Photo-Optical Instrumentation Engineers.

Detection of Brillouin scattering temperature signal in Brillouin optical time-domain reflectometer sensing system based on instantaneous frequency measurement technology

abstract {We present a simple and practical method for the single-ended distributed fiber temperature measurements using microwave (11-GHz) coherent detection and the instantaneous frequency measurement (IFM) technique to detect spontaneous Brillouin backscattered signal in which a specially designed rf bandpass filter at 11 GHz is used as a frequency discriminator to transform frequency shift to intensity fluctuation. A Brillouin temperature signal can be obtained at 11 GHz over a sensing length of 10 km. The power sensitivity dependence on temperature induced by frequency shift is measured as 2.66%/K. © 2007 Society of Photo-Optical Instrumentation Engineers.}

Rheological destructuring of aqueous gels composed of cellulose ethers following storage in the presence of redox agents

Despite their widespread use, there is a paucity of information concerning the effect of storage on the rheological properties of pharmaceutical gels that contain organic and inorganic additives. Therefore, this study examined the effect of storage (1 month at either 4 or 37 degrees C) on the rheological and mechanical properties of gels composed of either hydroxypropylmethylcellulose (3-5% w/w, HPMC) or hydroxyethylcellulose (3-5% w/w, HEC) and containing or devoid of dispersed organic (tetracycline hydrochloride 2% w/w) or inorganic (iron oxide 0.1% w/w) agents. The mechanical properties were measured using texture profile analysis whereas the rheological properties were analyzed using continuous shear rheometry and modeled using the Power Law model. All formulations exhibited pseudoplastic flow with minimal thixotropy. Increasing polymer concentration (3-5% w/w) significantly increased the consistency, hardness, compressibility, and adhesiveness of the formulations due to increased polymer chain entanglement. Following storage (I month at 4 and 37 degrees C) the consistency and mechanical properties of additive free HPMC gets (but not HEC gels) increased, due to the time-dependent development of polymer chain entanglements. Incorporation of tetracycline hydrochloride significantly decreased and increased the rheological and mechanical properties of HPMC and HEC gels, respectively. Conversely, the incorporation of iron oxide did not affect these properties. Following storage, the rheological and mechanical properties of HPMC and HEC formulations were markedly compromised. This effect was greater following storage at 37 than at 4 degrees C and, additionally, greater in the presence of tetracycline hydrochloride than iron oxide. It is suggested that the loss of rheological/mechanical structure was due to chain depolymerization, facilitated by the redox properties of tetracycline hydrochloride and iron oxide. These observations have direct implications for the design and formulation of gels containing an active pharmaceutical ingredient. (c) 2005 Wiley Periodicals, Inc.

Tripping Control for Transient Stability in Coordinated Hydro and Wind Generation

Over-frequency generator tripping (OFGT) is used to cut off extra generation to balance power and loads in an isolated system. In this paper the impact of OGFT as a consequence of grid-connected wind farms and under-frequency load shedding (UFLS) is analysed. The paper uses a power system model to demonstrate that wind power fluctuations can readily render OFGT and UFLS maloperation. Using combined hydro and wind generation, the paper proposes a coordinated strategy which resolves problems associated with OFGT and UFLS and preserves system stability.

A Wavelet Approach to the Selective Computation of Domninant Eigenvalues

This paper introduces an algorithm that calculates the dominant eigenvalues (in terms of system stability) of a linear model and neglects the exact computation of the non-dominant eigenvalues. The method estimates all of the eigenvalues using wavelet based compression techniques. These estimates are used to find a suitable invariant subspace such that projection by this subspace will provide one containing the eigenvalues of interest. The proposed algorithm is exemplified by application to a power system model.

Fractional-order impulse response of the respiratory system

This paper presents the measurement, frequency-response modeling and identification, and the corresponding impulse time response of the human respiratory impedance and admittance. The investigated adult patient groups were healthy, diagnosed with chronic obstructive pulmonary disease and kyphoscoliosis, respectively. The investigated children patient groups were healthy, diagnosed with asthma and cystic fibrosis, respectively. Fractional order (FO) models are identified on the measured impedance to quantify the respiratory mechanical properties. Two methods are presented for obtaining and simulating the time-domain impulse response from FO models of the respiratory admittance: (i) the classical pole-zero interpolation proposed by Oustaloup in the early 90s, and (ii) the inverse discrete Fourier Transform (DFT). The results of the identified FO models for the respiratory admittance are presented by means of their average values for each group of patients. Consequently, the impulse time response calculated from the frequency response of the averaged FO models is given by means of the two methods mentioned above. Our results indicate that both methods provide similar impulse response data. However, we suggest that the inverse DFT is a more suitable alternative to the high order transfer functions obtained using the classical Oustaloup filter. Additionally, a power law model is fitted on the impulse response data, emphasizing the intrinsic fractal dynamics of the respiratory system.

La sous-traitance « boomerang » : l’influence des relations industrielles sur l’internalisation d’activités manufacturières confiées en sous-traitance

Au cours de la dernière décennie, le secteur manufacturier québécois a subi d’importantes transformations. La mondialisation ayant amené les entreprises à restructurer leurs activités, 96% des établissements manufacturiers québécois ont choisi d’avoir recours à la sous-traitance (Jalette, 2004). Or, nous savons que l’octroi de contrats en sous-traitance est un processus qui se révèle parfois décevant pour les entreprises. En raison de l’évolution de certains facteurs organisationnels internes ou externes qui contribue à changer le contexte au sein duquel les décisions de sous-traitance avaient été prises initialement (ex. : augmentation des coûts de main-d’œuvre et de transport, changements technologiques), il arrive que les entreprises renversent les décisions de sous-traitance par une décision d’internalisation (OCDE, 2013). Le concept d’internalisation réfère ici à l’action d’une entreprise de ramener vers un établissement des activités de production qui avaient été antérieurement transférées en dehors de celui-ci, et ce, que ce soit en raison d’une sous-traitance d’activités à une entreprise indépendante ou d’une délocalisation au sein d’un autre établissement de l’entreprise (Jalette & Chevance, 2008). La présente étude a pour objectif d’expliquer la décision managériale d’internalisation par le biais des dynamiques sociales inhérentes aux relations industrielles, en allant au-delà des seuls facteurs économiques. L’effet sur l’internalisation de quatre variables indépendantes liées aux relations industrielles est étudié, soit la capacité stratégique du syndicat local, la solidarité interne de l’acteur syndical local et sa solidarité externe, ainsi que le climat des relations patronales-syndicales. La typologie des positions syndicales de Jalette (2005), le modèle des ressources de pouvoir syndicales de Lévesque et Murray (2003) et les travaux de Dastmalchian (2008) traitant du climat des relations industrielles constituent le cadre théorique de cette recherche. La présente étude utilise des données recueillies dans le cadre d’une enquête par questionnaire réalisée en 2005 auprès de présidents de syndicats locaux affiliés à la CSD et à la FTQ qui œuvrent dans le secteur manufacturier québécois. L’étude révèle que 24% des syndicats locaux sondés ont vécu l’internalisation d’activités antérieurement sous-traitées dans les deux années précédentes au sondage. De plus, les résultats sont à l’effet que les chances d’internalisation d’activités augmentent lorsque plusieurs stratégies et actions visant l’internalisation sont adoptées par le syndicat local, lorsque la vie syndicale d’un syndicat local s’améliore et lorsque le temps de libération syndicale d’un syndicat local augmente. Toutefois, les chances d’internalisation d’activités diminuent lorsque le rapport de force d’un syndical local augmente. L’internalisation d’activités diminue également lorsque le climat des relations du travail se détériore. Bref, les résultats de cette recherche montrent que les syndicats n’ont pas qu’à subir la mondialisation et que, par leurs actions, ils sont à même d’en influencer le cours, notamment en favorisant l’internalisation d’activités antérieurement confiées en sous-traitance.

Modeling of complex-valued Wiener systems using B-spline neural network

In this brief, a new complex-valued B-spline neural network is introduced in order to model the complex-valued Wiener system using observational input/output data. The complex-valued nonlinear static function in the Wiener system is represented using the tensor product from two univariate B-spline neural networks, using the real and imaginary parts of the system input. Following the use of a simple least squares parameter initialization scheme, the Gauss-Newton algorithm is applied for the parameter estimation, which incorporates the De Boor algorithm, including both the B-spline curve and the first-order derivatives recursion. Numerical examples, including a nonlinear high-power amplifier model in communication systems, are used to demonstrate the efficacy of the proposed approaches.

Comments about some species abundance patterns: classic, neutral, and niche partitioning models

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

Influência das cargas selecionadas na determinação da potência crítica determinada no ergômetro de braço em dois modelos lineares

O objetivo deste estudo foi verificar o efeito da seleção das cargas e do modelo utilizado para a determinação da PC no ergômetro de braço. Participaram do estudo oito voluntários do sexo masculino, que praticavam atividade física regularmente e eram aparentemente saudáveis. Os sujeitos realizaram quatro testes com cargas constantes mantidas até a exaustão voluntária no ergômetro de braço UBE 2462-Cybex. As cargas foram individualmente selecionadas para induzir a fadiga entre 1 e 15 minutos. Para cada sujeito, a determinação da PC foi realizada através de dois modelos lineares: potência-1/tempo e trabalho-tempo. em cada um dos modelos, foram utilizadas todas as potências (1), as três maiores (2) e as três menores (3). As PC encontradas no modelo potência-1/tempo e trabalho-tempo para a condição 3 (177,5 + 29,5; 173,9 + 33,3, respectivamente) foram significantemente menores do que as da condição 2 (190,5 + 23,2; 183,4 + 22,3, respectivamente), não existindo diferenças destas com as da condição 1 (184,2 + 25,4; 176,4 + 28,8, respectivamente). As PC determinadas no modelo potência-1/tempo para as condições 1 e 2 foram significantemente maiores do que as determinadas no modelo trabalho-tempo, não existindo diferenças para a condição 3. Pode-se concluir que as cargas selecionadas e o modelo utilizado interferem na determinação da PC encontrada no ergômetro de braço, podendo interferir no tempo de exaustão durante o exercício submáximo realizado em cargas relativas a este índice.

Rheological and thermophysical properties of blackberry juice

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

Density and rheological parameters of goat milk

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

Non-newtonian flow and pressure drop of pineapple juice in a plate heat exchanger

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

Friction losses in valves and fittings for liquid food products

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