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.

Numerical modelling of double diffusion driven reactive flow transport in deformable fluid-saturated porous media with particular consideration of temperature-dependent chemical reaction rates

Numerical methods ave used to solve double diffusion driven reactive flow transport problems in deformable fluid-saturated porous media. in particular, thp temperature dependent reaction rate in the non-equilibrium chemical reactions is considered. A general numerical solution method, which is a combination of the finite difference method in FLAG and the finite element method in FIDAP, to solve the fully coupled problem involving material deformation, pore-fluid flow, heat transfer and species transport/chemical reactions in deformable fluid-saturated porous media has been developed The coupled problem is divided into two subproblems which are solved interactively until the convergence requirement is met. Owing to the approximate nature of the numerical method, if is essential to justify the numerical solutions through some kind of theoretical analysis. This has been highlighted in this paper The related numerical results, which are justified by the theoretical analysis, have demonstrated that the proposed solution method is useful for and applicable to a wide range of fully coupled problems in the field of science and engineering.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

An assumption-driven case-specific model editor

A case sensitive intelligent model editor has been developed for constructing consistent lumped dynamic process models and for simplifying them using modelling assumptions. The approach is based on a systematic assumption-driven modelling procedure and on the syntax and semantics of process,models and the simplifying assumptions.

Functional significance of the beta-hairpin in the insecticidal neurotoxin omega-atracotoxin-Hv1a

omega -Atracotoxin-Hv1a is an insect-specific neurotoxin whose phylogenetic specificity derives from its ability to antagonize insect, but not vertebrate, voltage-gated calcium channels. In order to help understand its mechanism of action and to enhance its utility as a lead compound for insecticide development, we used a combination of protein engineering and site-directed mutagenesis to probe the toxin for key functional regions. First, we constructed a Hairpinless mutant in which the C-terminal beta -hairpin, which is highly conserved in this family of neurotoxins, was excised without affecting the fold of the residual disulfide-rich core of the toxin. The Hairpinless mutant was devoid of insecticidal activity, indicating the functional importance of the hairpin. We subsequently developed a highly efficient system for production of recombinant toxin and then probed the hairpin for key functional residues using alanine-scanning mutagenesis followed by a second round of mutagenesis based on initial hits from the alanine scan. This revealed that two spatially proximal residues, Asn(27) and Arg(35), form a contiguous molecular surface that is essential for toxin activity. We propose that this surface of the beta -hairpin is a key site for interaction of the toxin with insect calcium channels.

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.

Characterization of optically driven fluid stress fields with optical tweezers

We present a controlled stress microviscometer with applications to complex fluids. It generates and measures microscopic fluid velocity fields, based on dual beam optical tweezers. This allows an investigation of bulk viscous properties and local inhomogeneities at the probe particle surface. The accuracy of the method is demonstrated in water. In a complex fluid model (hyaluronic acid), we observe a strong deviation of the flow field from classical behavior. Knowledge of the deviation together with an optical torque measurement is used to determine the bulk viscosity. Furthermore, we model the observed deviation and derive microscopic parameters.

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.

A new algorithm and refined bounds for extended gcd computation

Extended gcd computation is interesting itself. It also plays a fundamental role in other calculations. We present a new algorithm for solving the extended gcd problem. This algorithm has a particularly simple description and is practical. It also provides refined bounds on the size of the multipliers obtained.

Resonance fluorescence and Autler-Townes spectra of a two-level atom driven by two fields of equal frequencies

We study the effects of driving a two-level atom by two intense field modes that have equal frequencies but are otherwise distinguishable; the intensity of one mode is also assumed to be greater than that of the other. We calculate first the dressed states of the system, and then its resonance fluorescence and Autler-Townes absorption spectra. We find that the energy spectrum of the doubly dressed atom consists of a ladder of doublet continua. These continua manifest themselves in the fluorescence spectrum, where they produce continua at the positions of the Mellow sideband frequencies omega(L)+/-2 Omega of the strong field, and in the Autler-Townes absorption spectrum, which becomes a two-continuum doublet.

Resonance fluorescence spectrum of a two-level atom driven by a bichromatic field in a squeezed vacuum

The steady-state resonance fluorescence spectrum of a two-level atom driven by a bichromatic field in a broadband squeezed vacuum is studied. When the carrier frequency of the squeezed vacuum is tuned to the frequency of the central spectral line, anomalous spectral features, such as hole burning and dispersive profiles, can occur at the central line. We show that these features appear for wider, and experimentally more convenient, ranges of the parameters than in the case of monochromatic excitation. ?he absence of a coherent spectral component at the central line makes any experimental attempt to observe these features much easier. We also discuss the general features of the spectrum. When the carrier frequency of the squeezed vacuum is tuned to the first odd or even sidebands, the spectrum is asymmetric and only the sidebands an sensitive to phase. For appropriate choices of the phase the linewidths or only the odd or even sidebands can be reduced. A dressed-stale interpretation is provided.

Transport in bidisperse adsorbents: Significance of the macroscopic adsorbate flux

Transport in bidisperse adsorbents is investigated here, while incorporating a two-dimensional model for adsorbate diffusion in the microparticles. The latter treatment permits consideration of the macropore concentration variation around the microparticle surface, and thereby predicts an adsorbate through-flux on the macroscopic coordinate. Such a through-flux has earlier been postulated in the literature, but with unrealistic mechanistic justification. The new model therefore resolves the existing ambiguity in this regard, and covers the entire spectrum of behaviour between microparticle and macropore diffusion control. Computational results show that if the macroscopic adsorbate flux, ignored in the conventional analysis, has a significant contribution to the total flux under macropore control conditions then it is always important even when the microparticle diffusion resistance is not negligible. The effect of various parameters such as relative microparticle size and isotherm heterogeneity on the uptake is also studied and discussed. (C) 1997 Elsevier Science Ltd.

Effect of the Bloch-Siegert shift in a strongly driven transition: Asymmetric Autler-Townes profile

Under the conditions of the rotating wave approximation (RWA), a transition strongly driven by a resonant oscillating field displays the well known symmetric Autler-Townes doublet. However, if the counter-rotating component, neglected in the RWA, is taken into account, the Bloch-Siegert shift gives rise to an Autler-Townes doublet of unequal intensity even in the case of a resonant driving field. This effect is investigated theoretically in a V-shaped three-level double-resonance configuration and the results are presented in this paper. An interesting observation is that the level of asymmetry not only depends on the driving-field intensity but also on the characteristics of the driven system including relaxation rates and equilibrium population distributions.

Alternatives for parallel Krylov subspace basis computation

Numerical methods related to Krylov subspaces are widely used in large sparse numerical linear algebra. Vectors in these subspaces are manipulated via their representation onto orthonormal bases. Nowadays, on serial computers, the method of Arnoldi is considered as a reliable technique for constructing such bases. However, although easily parallelizable, this technique is not as scalable as expected for communications. In this work we examine alternative methods aimed at overcoming this drawback. Since they retrieve upon completion the same information as Arnoldi's algorithm does, they enable us to design a wide family of stable and scalable Krylov approximation methods for various parallel environments. We present timing results obtained from their implementation on two distributed-memory multiprocessor supercomputers: the Intel Paragon and the IBM Scalable POWERparallel SP2. (C) 1997 by John Wiley & Sons, Ltd.