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.

NMR shim coil design utilising a rapid spherical harmonic calculation method

A rapid spherical harmonic calculation method is used for the design of Nuclear Magnetic Resonance shim coils. The aim is to design each shim such that it generates a field described purely by a single spherical harmonic. By applying simulated annealing techniques, coil arrangements are produced through the optimal positioning of current-carrying circular arc conductors of rectangular cross-section. This involves minimizing the undesirable harmonies in relation to a target harmonic. The design method is flexible enough to be applied for the production of coil arrangements that generate fields consisting significantly of either zonal or tesseral harmonics. Results are presented for several coil designs which generate tesseral harmonics of degree one.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Fuel stagnation temperature effects on mixing with supersonic combustion flows

An experimental study of the effect of fuel stagnation temperature on mixing in a supersonic hydrogen-air flame is described, The combustor consisted of a constant-area rectangular duct with a centrally located fuel-injection strut that spanned the width. A high-enthalpy stream of air was supplied by a free-piston shock tunnel, and heated hydrogen fuel, supplied by a gun-tunnel, was injected into the freestream as a coflowing planar jet. The freestream total enthalpies were 5.6, 6.5, and 9 MJ/kg, and fuel stagnation temperatures were 300, 450, and 700 K, Raising the fuel stagnation temperature increased the fuel velocity to be near that of the airstream and resulted in a decrease in the mixing rate, Even as the fuel and air velocities became equal, significant mixing still occurred because of a large difference in density, Increasing the freestream enthalpy reduced the difference between the initial air temperature and the adiabatic flame temperature, which in turn reduced the heat addition, and subsequently, the amount of pressure rise in the duct.

Robust algorithms for solving stochastic partial differential equations

A robust semi-implicit central partial difference algorithm for the numerical solution of coupled stochastic parabolic partial differential equations (PDEs) is described. This can be used for calculating correlation functions of systems of interacting stochastic fields. Such field equations can arise in the description of Hamiltonian and open systems in the physics of nonlinear processes, and may include multiplicative noise sources. The algorithm can be used for studying the properties of nonlinear quantum or classical field theories. The general approach is outlined and applied to a specific example, namely the quantum statistical fluctuations of ultra-short optical pulses in chi((2)) parametric waveguides. This example uses a non-diagonal coherent state representation, and correctly predicts the sub-shot noise level spectral fluctuations observed in homodyne detection measurements. It is expected that the methods used wilt be applicable for higher-order correlation functions and other physical problems as well. A stochastic differencing technique for reducing sampling errors is also introduced. This involves solving nonlinear stochastic parabolic PDEs in combination with a reference process, which uses the Wigner representation in the example presented here. A computer implementation on MIMD parallel architectures is discussed. (C) 1997 Academic Press.

Beach water table fluctuations due to wave run-up: Capillarity effects

High-frequency beach water table fluctuations due to wave run-up and rundown have been observed in the field [Waddell, 1976]. Such fluctuations affect the infiltration/exfiltration process across the beach face and the interstitial oxygenation process in the beach ecosystem. Accurate representation of high-frequency water table fluctuations is of importance in the modeling of (1) the interaction between seawater and groundwater, more important, the effects on swash sediment transport and (2) the biological activities in the beach ecosystem. Capillarity effects provide a mechanism for high-frequency water table fluctuations. Previous modeling approaches adopted the assumption of saturated flow only and failed to predict the propagation of high-frequency fluctuations in the aquifer. In this paper we develop a modified kinematic boundary condition (kbc) for the water table which incorporates capillarity effects. The application of this kbc in a boundary element model enables the simulation of high-frequency water table fluctuations due to wave run-up. Numerical tests were carried out for a rectangular domain with small-amplitude oscillations; the behavior of water table responses was found to be similar to that predicted by an analytical solution based on the one-dimensional Boussinesq equation. The model was also applied to simulate the water table response to wave run-up on a doping beach. The results showed similar features of water table fluctuations observed in the field. In particular, these fluctuations are standing wave-like with the amplitude becoming increasingly damped inland. We conclude that the modified kbc presented here is a reasonable approximation of capillarity effects on beach water table fluctuations. However, further model validation is necessary before the model can confidently be used to simulate high-frequency water table fluctuations due to wave run-up.

Rapid computation of static fields produced by thick circular solenoids

A straightforward method is proposed for computing the magnetic field produced by a circular coil that contains a large number of turns wound onto a solenoid of rectangular cross section. The coil is thus approximated by a circular ring containing a continuous constant current density, which is very close to the real situation when sire of rectangular cross section is used. All that is required is to evaluate two functions, which are defined as integrals of periodic quantities; this is done accurately and efficiently using trapezoidal-rule quadrature. The solution can be obtained so rapidly that this procedure is ideally suited for use in stochastic optimization, An example is given, in which this approach is combined with a simulated annealing routine to optimize shielded profile coils for NMR.

The structures of Langmuir-Blodgett films of fatty acids and their salts

Recent advances in several experimental techniques have enabled detailed structural information to be obtained for floating (Langmuir) monolayers and Langmuir-Blodgett films. These techniques are described briefly and their application to the study of films of fatty acids and their salts is discussed. Floating monolayers on aqueous subphases have been shown to possess a complex polymorphism with phases whose structures may be compared to those of smectic mesophases. However, only those phases that exist at high surface pressures are normally used in Langmuir-Blodgett (LB) deposition. In single LB monolayers of fatty acids and fatty acid salts the acyl chains are in the all-cans conformation with their long axes normal to the substrate. The in-plane molecular packing is hexagonal with long-range bond orientational order and short-range positional order: known as the hexatic-B structure. This structure is found irrespective of the phase of the parent floating monolayer. The structures of multilayer LB films are similar to the structures of their bulk crystals, consisting of stacked bilayer lamellae. Each lamella is formed from two monolayers of fatty acid molecules or ions arranged head to head and held together by hydrogen bonding between pairs of acids or ionic bonding through the divalent cations. With acids the acyl chains are tilted with respect to the substrate normal and have a monoclinic structure, whereas the salts with divalent cations may have the chains normal to the substrate or tilted. The in-plane structures are usually centred rectangular with the chains in the trans conformation and packed in a herringbone pattern, Multilayer films of the acids show only a single-step order-disorder transition at the malting point, This temperature tends to rise as the number of layers increases. Complex changes occur when multilayer films of the salts are heated. Disorder of the chains begins at low temperatures but the arrangement of the head groups does not alter until the melting temperature is reached, Slow heating to a temperature just below the melting temperature gives, with some salts, a radical change in phase. The lamellar structure disappears and a new phase consisting of cylindrical rods lying parallel to the substrate surface and stacked in a hexagonal pattern is formed, In each rod the cations are aligned along the central axis surrounded by the disordered acyl chains. (C) 2001 Elsevier Science B,V. All rights reserved.

Effect of frother type and concentration on the water recovery and entrainment recovery relationship

A study was undertaken to determine the effects of five industrial frothers on the relationship between entrainment recovery and water recovery. The experiments were conducted using a Batequip rectangular 60-L steel flotation cell that was run in parallel with the first cell in a copper/silver prefloat rougher circuit. The concentration of the frothers ranged from 2.5 to 8.4 g/t. The experiments were repeated at three different froth depths. A power function was found to fit this relationship with no discontinuity in the relationship caused by the changes in frother type and concentration and froth depth.

X-ray scattering studies of mixed langmuir monolayers and Langmuir-Blodgett films of a noncentrosymmetric porphyrin with cadmium arachidate

The structures of mixed Langmuir (floating) monolayers and Langmuir-Blodgett (LB) films of a phenanthroline-porphyrin with cadmium arachidate (PhenPor + CdAr) have been investigated by synchrotron X-ray grazing incidence diffraction (GIXD) and specular X-ray reflectivity (SXR). GIXD measurements of the floating monolayers showed only one peak, arising from the CdAr domains in the films, at a scattering angle of 21.5 degrees. This is consistent with a hexagonal structure (alpha = 4.77 Angstrom). The correlation length in these domains is 250 Angstrom. GMD measurements of the LB films, however, show two sets of diffraction features: one arises from CdAr domains with a rectangular in-plane structure (alpha = 7.44 Angstrom and b = 4.90 Angstrom) and a correlation length of 85 Angstrom; the other is from porphyrin domains with an oblique in-plane structure (alpha (p) 15.2 Angstrom, b(p) = 8.86 Angstrom, and gamma (p) = 80 degrees) and a correlation length of 105 Angstrom. These dimensions are consistent with the surface pressure-area isotherm measurements and indicate that the two components are immiscible. The thickness of the bilayer is 57 Angstrom, and there is no correlation between the bilayers. Introduction of a trigger compound does not alter the structure of the films but slightly increases the bilayer thickness. The SXR measurements of the floating monolayers also support the suggested immiscibility of the two components in the films.

How cracks in SiOx-coated polyester films affect gas permeation

In this paper theoretical models have been established that can account for the gas transmission through nanocomposite laminates, consisting of an oxide layer of finite permeability containing defects, on a polymer sheet of finite thickness. The defect shapes can either be in the form of long cracks or rectangular holes. The models offer a choice of exact numerical calculations or fast and intuitive analytical approximations. The experimental measurements of oxygen permeation through four different SiOx/poly (ethylene terephthalate) samples that were strained to produce distributions or cracks showed good agreement when compared with predicted results from the approximate analytic model. As a consequence of this observation, a key practical conclusion is that, because of the logarithmic dependence of transmission on the width of a crack, for a given strain it is better to have a small number of large cracks rather than a large number of small cracks. (C) 2001 Elsevier Science B.V. All rights reserved.

Exploring through cover - the integrated interpretation of high resolution aeromagnetic, airborne electromagnetic and ground gravity data from the Grant's Patch area, Eastern Goldfields Province, Archaean Yigarn Craton Part B; Gravity inversion as a

The cost and risk associated with mineral exploration in Australia increases significantly as companies move into deeper regolith-covered terrain. The ability to map the bedrock and the depth of weathering within an area has the potential to decrease this risk and increase the effectiveness of exploration programs. This paper is the second in a trilogy concerning the Grant's Patch area of the Eastern Goldfields. The recent development of the VPmg potential field inversion program in conjunction with the acquisition of high-resolution gravity data over an area with extensive drilling provided an opportunity to evaluate three-dimensional gravity inversion as a bedrock and regolith mapping tool. An apparent density model of the study area was constructed, with the ground represented as adjoining 200 m by 200 m vertical rectangular prisms. During inversion VPmg incrementally adjusted the density of each prism until the free-air gravity response of the model replicated the observed data. For the Grant's Patch study area, this image of the apparent density values proved easier to interpret than the Bouguer gravity image. A regolith layer was introduced into the model and realistic fresh-rock densities assigned to each basement prism according to its interpreted lithology. With the basement and regolith densities fixed, the VPmg inversion algorithm adjusted the depth to fresh basement until the misfit between the calculated and observed gravity response was minimised. The resulting geometry of the bedrock/regolith contact largely replicated the base of weathering indicated by drilling with predicted depth of weathering values from gravity inversion typically within 15% of those logged during RAB and RC drilling.

Calorimetric demonstration of the potential of molecular crowding to emulate the effect of an allosteric activator on pyruvate kinase kinetics

A method based on isothermal calorimetry is described for the direct kinetic assay of pyruvate kinase. In agreement with earlier findings based on the standard coupled assay system for this enzyme in the presence of a fixed ADP concentration, the essentially rectangular hyperbolic dependence of initial velocity upon phosphoenolpyruvate concentration is rendered sigmoidal by the allosteric inhibitor phenylalanine. This effect of phenylalanine can be countered by including a high concentration of a space- filling osmolyte such as proline in the reaction mixtures. This investigation thus affords a dramatic example that illustrates the need to consider potential consequences of thermodynamic nonideality on the kinetics of enzyme reactions in crowded molecular environments such as the cell cytoplasm.

Skin-friction measurements in high-enthalpy hypersonic boundary layers

Skin-friction measurements are reported for high-enthalpy and high-Mach-number laminar, transitional and turbulent boundary layers. The measurements were performed in a free-piston shock tunnel with air-flow Mach number, stagnation enthalpy and Reynolds numbers in the ranges of 4.4-6.7, 3-13 MJ kg(-1) and 0.16 x 10(6)-21 x 10(6), respectively. Wall temperatures were near 300 K and this resulted in ratios of wall enthalpy to flow-stagnation enthalpy in the range of 0.1-0.02. The experiments were performed using rectangular ducts. The measurements were accomplished using a new skin-friction gauge that was developed for impulse facility testing. The gauge was an acceleration compensated piezoelectric transducer and had a lowest natural frequency near 40 kHz. Turbulent skin-friction levels were measured to within a typical uncertainty of +/-7%. The systematic uncertainty in measured skin-friction coefficient was high for the tested laminar conditions; however, to within experimental uncertainty, the skin-friction and heat-transfer measurements were in agreement with the laminar theory of van Driest (1952). For predicting turbulent skin-friction coefficient, it was established that, for the range of Mach numbers and Reynolds numbers of the experiments, with cold walls and boundary layers approaching the turbulent equilibrium state, the Spalding & Chi (1964) method was the most suitable of the theories tested. It was also established that if the heat transfer rate to the wall is to be predicted, then the Spalding & Chi (1964) method should be used in conjunction with a Reynolds analogy factor near unity. If more accurate results are required, then an experimentally observed relationship between the Reynolds analogy factor and the skin-friction coefficient may be applied.