Computation of bifurcation boundaries for power systems: A new Delta-plane method

This paper is devoted to the problems of finding the load flow feasibility, saddle node, and Hopf bifurcation boundaries in the space of power system parameters. The first part contains a review of the existing relevant approaches including not-so-well-known contributions from Russia. The second part presents a new robust method for finding the power system load flow feasibility boundary on the plane defined by any three vectors of dependent variables (nodal voltages), called the Delta plane. The method exploits some quadratic and linear properties of the load now equations and state matrices written in rectangular coordinates. An advantage of the method is that it does not require an iterative solution of nonlinear equations (except the eigenvalue problem). In addition to benefits for visualization, the method is a useful tool for topological studies of power system multiple solution structures and stability domains. Although the power system application is developed, the method can be equally efficient for any quadratic algebraic problem.

Sensitivity to measurement perturbation of single-atom dynamics in cavity QED

We consider continuous observation of the nonlinear dynamics of single atom trapped in an optical cavity by a standing wave with intensity modulation. The motion of the atom changes the phase of the field which is then monitored by homodyne detection of the output field. We show that the conditional Hilbert space dynamics of this system, subject to measurement-induced perturbations, depends strongly on whether the corresponding classical dynamics is regular or chaotic. If the classical dynamics is chaotic, the distribution of conditional Hilbert space vectors corresponding to different observation records tends to be orthogonal. This is a characteristic feature of hypersensitivity to perturbation for quantum chaotic systems.

Simultaneous steady state and dynamic design of process systems using structural indicators

The simultaneous design of the steady-state and dynamic performance of a process has the ability to satisfy much more demanding dynamic performance criteria than the design of dynamics only by the connection of a control system. A method for designing process dynamics based on the use of a linearised systems' eigenvalues has been developed. The eigenvalues are associated with system states using the unit perturbation spectral resolution (UPSR), characterising the dynamics of each state. The design method uses a homotopy approach to determine a final design which satisfies both steady-state and dynamic performance criteria. A highly interacting single stage forced circulation evaporator system, including control loops, was designed by this method with the goal of reducing the time taken for the liquid composition to reach steady-state. Initially the system was successfully redesigned to speed up the eigenvalue associated with the liquid composition state, but this did not result in an improved startup performance. Further analysis showed that the integral action of the composition controller was the source of the limiting eigenvalue. Design changes made to speed up this eigenvalue did result in an improved startup performance. The proposed approach provides a structured way to address the design-control interface, giving significant insight into the dynamic behaviour of the system such that a systematic design or redesign of an existing system can be undertaken with confidence.

An L-system-based plant modeling language

Cpfg is a program for simulating and visualizing plant development, based on the theory of L-systems. A special-purpose programming language, used to specify plant models, is an essential feature of cpfg. We review postulates of L-system theory that have influenced the design of this language. We then present the main constructs of this language, and evaluate it from a user's perspective.

L-studio/cpfg: A software system for modeling plants

L-studio/cpfg is a plant modeling software system designed for Windows 95/98/NT platforms. Its key components are the L-system-based plant simulator cpfg and the modeling environment called L-studio. We overview version 1.0 of this system from the user's perspective.

Analysis of steady-state heat transfer through mid-crustal vertical cracks with upward throughflow in hydrothermal systems

We conduct a theoretical analysis of steady-state heat transfer problems through mid-crustal vertical cracks with upward throughflow in hydrothermal systems. In particular, we derive analytical solutions for both the far field and near field of the system. In order to investigate the contribution of the forced advection to the total temperature of the system, two concepts, namely the critical Peclet number and the critical permeability of the system, have been presented and discussed in this paper. The analytical solution for the far field of the system indicates that if the pore-fluid pressure gradient in the crust is lithostatic, the critical permeability of the system can be used to determine whether or not the contribution of the forced advection to the total temperature of the system is negligible. Otherwise, the critical Peclet number should be used. For a crust of moderate thickness, the critical permeability is of the order of magnitude of 10(-20) m(2), under which heat conduction is the overwhelming mechanism to transfer heat energy, even though the pore-fluid pressure gradient in the crust is lithostatic. Furthermore, the lower bound analytical solution for the near field of the system demonstrates that the permeable vertical cracks in the middle crust can efficiently transfer heat energy from the lower crust to the upper crust of the Earth. Copyright (C) 2002 John Wiley Sons, Ltd.

Effect of material anisotropy on the onset of convective flow in three-dimensional fluid-saturated faults

In order to investigate the effect of material anisotropy on convective instability of three-dimensional fluid-saturated faults, an exact analytical solution for the critical Rayleigh number of three-dimensional convective flow has been obtained. Using this critical Rayleigh number, effects of different permeability ratios and thermal conductivity ratios on convective instability of a vertically oriented three-dimensional fault have been examined in detail. It has been recognized that (1) if the fault material is isotropic in the horizontal direction, the horizontal to vertical permeability ratio has a significant effect on the critical Rayleigh number of the three-dimensional fault system, but the horizontal to vertical thermal conductivity ratio has little influence on the convective instability of the system, and (2) if the fault material is isotropic in the fault plane, the thermal conductivity ratio of the fault normal to plane has a considerable effect on the critical Rayleigh number of the three-dimensional fault system, but the effect of the permeability ratio of the fault normal to plane on the critical Rayleigh number of three-dimensional convective flow is negligible.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Systemic Model for Diagnosis of the Micro, Small and Medium Enterprises from Two Cities from the Countryside of the State of So Paulo in Brazil

This article presents a proposal of a systemic model composed for the micro and small companies (MSE) of the region of Ribeiro Preto and the agents which influenced their environment. The proposed model was based on Stafford Beer`s (Diagnosing the system for organizations. Chichester, Wiley, 1985) systemic methodologies VSM (Viable System Model) and on Werner Ulrich`s (1983) CSH (Critical Systems Heuristics). The VSM is a model for the diagnosis of the structure of an organization and of its flows of information through the application of the cybernetics concepts (Narvarte, In El Modelo del Sistema Viable-MSV: experiencias de su aplicacin en Chile. Proyecto Cerebro Colectivo del IAS, Santiago, 2001). On the other hand, CSH focus on the context of the social group applied to the systemic vision as a counterpoint to the organizational management view considered by the VSM. MSE of Ribeiro Preto and Sertozinho had been analyzed as organizations inserted in systems that relate and integrate with other systems concerning the public administration, entities of representation and promotion agencies. The research questions: which are the bonds of interaction among the subsystems in this process and who are the agents involved? The systemic approach not only diagnosed a social group, formed by MSE of Ribeiro Preto and Sertozinho, public authorities and support entities, but could also delineate answers that aimed the clarification of obscure questions generating financial assistance to the formularization of efficient actions for the development of this system.

Thoracic vibrations in stingless bees (Melipona seminigra): resonances of the thorax influence vibrations associated with flight but not those associated with sound production

Bees generate thoracic vibrations with their indirect flight muscles in various behavioural contexts. The main frequency component of non-flight vibrations, during which the wings are usually folded over the abdomen, is higher than that of thoracic vibrations that drive the wing movements for flight. So far, this has been concluded from an increase in natural frequency of the oscillating system in association with the wing adduction. In the present study, we measured the thoracic oscillations in stingless bees during stationary flight and during two types of non-flight behaviour, annoyance buzzing and forager communication, using laser vibrometry. As expected, the flight vibrations met all tested assumptions for resonant oscillations: slow build-up and decay of amplitude; increased frequency following reduction of the inertial load; and decreased frequency following an increase of the mass of the oscillating system. Resonances, however, do not play a significant role in the generation of non-flight vibrations. The strong decrease in main frequency at the end of the pulses indicates that these were driven at a frequency higher than the natural frequency of the system. Despite significant differences regarding the main frequency components and their oscillation amplitudes, the mechanism of generation is apparently similar in annoyance buzzing and forager vibrations. Both types of non-flight vibration induced oscillations of the wings and the legs in a similar way. Since these body parts transform thoracic oscillations into airborne sounds and substrate vibrations, annoyance buzzing can also be used to study mechanisms of signal generation and transmission potentially relevant in forager communication under controlled conditions.

A history of agricultural production at Neolithic Catalhoyuk East, Turkey

A history of agricultural production is proposed for Neolithic Catalhoyuk East, central Turkey, using archaeobotanical, environmental, population and settlement studies. In the aceramic early phase of site occupation, intensive strategies developed as changes in population and environment caused stress on food supplies produced within a limited territory. Food exchange may have been part of the social means by which Catalhoyuk and nearby contemporary settlements amalgamated into the single site of the main occupation phase. Population change, inherited territories and continuing environmental impact led to the development of an extensive system of agriculture using widely dispersed dry soils, with an intensive regime applied to nearby alluvial soils. Social tensions caused by the evolution of this system contributed to the fissioning of the site by the Chalcolithic.

Double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the double diffusion-driven convective instability of three-dimensional fluid-saturated geological fault zones when they are heated uniformly from below. The fault zone is assumed to be more permeable than its surrounding rocks. In particular, we have derived exact analytical solutions to the total critical Rayleigh numbers of the double diffusion-driven convective flow. Using the corresponding total critical Rayleigh numbers, the double diffusion-driven convective instability of a fluid-saturated three-dimensional geological fault zone system has been investigated. The related theoretical analysis demonstrates that: (1) The relative higher concentration of the chemical species at the top of the three-dimensional geological fault zone system can destabilize the convective flow of the system, while the relative lower concentration of the chemical species at the top of the three-dimensional geological fault zone system can stabilize the convective flow of the system. (2) The double diffusion-driven convective flow modes of the three-dimensional geological fault zone system are very close each other and therefore, the system may have the similar chance to pick up different double diffusion-driven convective flow modes, especially in the case of the fault thickness to height ratio approaching 0. (3) The significant influence of the chemical species diffusion on the convective instability of the three-dimensional geological fault zone system implies that the seawater intrusion into the surface of the Earth is a potential mechanism to trigger the convective flow in the shallow three-dimensional geological fault zone system.

Numerical simulation of double-diffusion driven convective flow and rock alteration in three-dimensional fluid-saturated geological fault zones

Numerical methods are used to simulate the double-diffusion driven convective pore-fluid flow and rock alteration in three-dimensional fluid-saturated geological fault zones. The double diffusion is caused by a combination of both the positive upward temperature gradient and the positive downward salinity concentration gradient within a three-dimensional fluid-saturated geological fault zone, which is assumed to be more permeable than its surrounding rocks. In order to ensure the physical meaningfulness of the obtained numerical solutions, the numerical method used in this study is validated by a benchmark problem, for which the analytical solution to the critical Rayleigh number of the system is available. The theoretical value of the critical Rayleigh number of a three-dimensional fluid-saturated geological fault zone system can be used to judge whether or not the double-diffusion driven convective pore-fluid flow can take place within the system. After the possibility of triggering the double-diffusion driven convective pore-fluid flow is theoretically validated for the numerical model of a three-dimensional fluid-saturated geological fault zone system, the corresponding numerical solutions for the convective flow and temperature are directly coupled with a geochemical system. Through the numerical simulation of the coupled system between the convective fluid flow, heat transfer, mass transport and chemical reactions, we have investigated the effect of the double-diffusion driven convective pore-fluid flow on the rock alteration, which is the direct consequence of mineral redistribution due to its dissolution, transportation and precipitation, within the three-dimensional fluid-saturated geological fault zone system. (c) 2005 Elsevier B.V. All rights reserved.

Analysis of dynamic models for structural insight and model reduction .2. A multi-stage compressor shutdown case-study

In this second paper, the three structural measures which have been developed are used in the modelling of a three stage centrifugal synthesis gas compressor. The goal of this case study is to determine the essential mathematical structure which must be incorporated into the compressor model to accurately model the shutdown of this system. A simple, accurate and functional model of the system is created via three structural measures. It was found that the model can be correctly reduced into its basic modes and that the order of the differential system can be reduced from 51(st) to 20(th). Of the 31 differential equational 21 reduce to algebraic relations, 8 become constants and 2 can be deleted thereby increasing the algebraic set from 70 to 91 equations. An interpretation is also obtained as to which physical phenomena are dominating the dynamics of the compressor add whether the compressor will enter surge during the shutdown. Comparisons of the reduced model performance against the full model are given, showing the accuracy and applicability of the approach. Copyright (C) 1996 Elsevier Science Ltd

Prediction of impact forces in a vibratory ball mill using an inverse technique

To understand the dynamic mechanisms of the mechanical milling process in a vibratory mill, it is necessary to determine the characteristics of the impact forces associated with the collision events. However, it is difficult to directly measure the impact force in an operating mill. This paper describes an inverse technique for the prediction of impact forces from acceleration measurements on a vibratory ball mill. The characteristics of the vibratory mill have been investigated by the modal testing technique, and its system modes have been identified. In the modelling of the system vibration response to the impact forces, two modal equations have been used to describe the modal responses. The superposition of the modal responses gives rise to the total response of the system. A method based on an optimisation approach has been developed to predict the impact forces by minimising the difference between the measured acceleration of the vibratory ball mill and the predicted acceleration from the solution of the modal equations. The predicted and measured impact forces are in good agreement. Copyright (C) 1996 Elsevier Science Ltd.