The displacement and strain field of three-dimensional rheologic model of earthquake preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Nine classes of integrable boundary conditions for the eight-state supersymmetric fermion model

Nine classes of integrable boundary conditions for the eight-state supersymmetric model of strongly correlated fermions are presented. The boundary systems are solved by using the coordinate Bethe ansatz method and the Bethe ansatz equations for all nine cases are given.

The Bazhanov-Stroganov model from 3D approach

We apply a three-dimensional approach to describe a new parametrization of the L-operators for the two-dimensional Bazhanov-Stroganov (BS) integrable spin model related to the chiral Potts model. This parametrization is based on the solution of the associated classical discrete integrable system. Using a three-dimensional vertex satisfying a modified tetrahedron equation, we construct an operator which generalizes the BS quantum intertwining matrix S. This operator describes the isospectral deformations of the integrable BS model.

A new two-parameter integrable model of strongly correlated fermions with quantum superalgebra symmetry

A new two-parameter integrable model with quantum superalgebra U-q[gl(3/1)] symmetry is proposed, which is an eight-state fermions model with correlated single-particle and pair hoppings as well as uncorrelated triple-particle hopping. The model is solved and the Bethe ansatz equations are obtained.

Separation and hedging results with state-contingent production

A state-contingent model of production under uncertainty is developed and compared with more traditional models of production under uncertainty. Producer behaviour with both production and price risk, in the presence and in the absence of futures and forward markets, is analysed in this state-contingent framework. Conditions for the optimal hedge to be positive or negative are derived. We also show that, under plausible conditions, a risk-averse producer facing price uncertainty and the ability to hedge price risk will never willingly adopt a nonstochastic technology. New separation results, which hold in the presence of both price and production risk, are then developed. These separation results generalize Townsend's spanning results by reducing the number of necessary forward markets by one.

Ring flexibility within tricyclic antidepressant drugs

The internal flexibility of the central seven-membered ring of a series of tricyclic antidepressant drugs (TCAs), imipramine {l}, amitriptyline {2}, doxepin {3}, and dothiepin {4}, has been investigated by H-1 and C-13 nuclear magnetic (NMR) techniques. Two dynamic processes were examined: ring inversion and bridge flexing. H-1 NMR lineshape analysis was used to obtain ring inversion barriers for 2-4. These studies yielded energy barriers of 14.3, 16.7, and 15.7 +/- 0.6 kcal/mol for the hydrochloride salts of doxepin, dothiepin, and amitriptyline, respectively. The barriers for the corresponding free bases were lower by 0.6 kcal/mol on average. (CT1)-C-13 relaxation measurements were used to determine the degree of bridge flexing associated with the central seven-membered ring for all four compounds. By fitting the T-1 data to a two-state jump model, lifetimes and amplitudes of rapid bridge flexing motions were determined. The results show that imipramine has the fastest rate of bridge flexing, followed by amitriptyline, doxepin, and dothiepin. The pharmacological profiles of the TCAs are complex and they interact with many receptor sites, resulting in numerous side effects and a general lack of understanding of their precise mode of action in different anxiety-related disorders. They all have similar three-dimensional structures, which makes it difficult to rationalize their differing relative potency in different assays/clinical settings. However, the clear finding here that there are significantly different degrees of internal mobility suggests that molecular dynamics should be an additional factor considered when trying to understand the mode of action of this clinically important family of molecules. (C) 2001 Wiley-Liss, Inc. and the American Pharmaceutical Association J Pharm Sci 90:713-721, 2001.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Why Scotland needs more than just a new migration policy

The Scottish Executive has adopted a policy to combat Scotland's declining population by encouraging inward migration. Using a multi-state population model this paper presents nine long-term population scenarios for Scotland using three net international migration levels and three fertility paths. The results show inward migration can slow population decline but makes little difference to population ageing. Without a higher fertility rate Scotland's population will become demographically unsustainable. Our simulations show that a higher fertility rate substantially reduces the future ageing.

Estimating trends and seasonality in coronary heart disease

We present two methods of estimating the trend, seasonality and noise in time series of coronary heart disease events. In contrast to previous work we use a non-linear trend, allow multiple seasonal components, and carefully examine the residuals from the fitted model. We show the importance of estimating these three aspects of the observed data to aid insight of the underlying process, although our major focus is on the seasonal components. For one method we allow the seasonal effects to vary over time and show how this helps the understanding of the association between coronary heart disease and varying temperature patterns. Copyright (C) 2004 John Wiley Sons, Ltd.

Dependence of transient dynamics in a class-C laser upon variation of inversion with time

The transient statistics of a gain-switched coherently pumped class-C laser displays a linear correlation between the first passage time and subsequent peak intensity. Measurements are reported showing a positive or negative sign of this linear correlation, controlled through the switching time and the laser detuning. Further measurements of the small-signal laser gain combined with calculations involving a three-level laser model indicate that this sign fundamentally depends upon the way the laser inversion varies during the gain switching, despite the added dynamics of the laser polarization in the class-C laser. [S1050-2947(97)07112-6].

Applying linear time-varying constraints to econometric models: With an application to demand systems

When linear equality constraints are invariant through time they can be incorporated into estimation by restricted least squares. If, however, the constraints are time-varying, this standard methodology cannot be applied. In this paper we show how to incorporate linear time-varying constraints into the estimation of econometric models. The method involves the augmentation of the observation equation of a state-space model prior to estimation by the Kalman filter. Numerical optimisation routines are used for the estimation. A simple example drawn from demand analysis is used to illustrate the method and its application.

A framework for extreme-event control in wastewater treatment

In this paper an approach to extreme event control in wastewater treatment plant operation by use of automatic supervisory control is discussed. The framework presented is based on the fact that different operational conditions manifest themselves as clusters in a multivariate measurement space. These clusters are identified and linked to specific and corresponding events by use of principal component analysis and fuzzy c-means clustering. A reduced system model is assigned to each type of extreme event and used to calculate appropriate local controller set points. In earlier work we have shown that this approach is applicable to wastewater treatment control using look-up tables to determine current set points. In this work we focus on the automatic determination of appropriate set points by use of steady state and dynamic predictions. The performance of a relatively simple steady-state supervisory controller is compared with that of a model predictive supervisory controller. Also, a look-up table approach is included in the comparison, as it provides a simple and robust alternative to the steady-state and model predictive controllers, The methodology is illustrated in a simulation study.

Growth and compaction behaviour of copper concentrate granules in a rotating drum

In order to understand the growth and compaction behaviour of chalcopyrite (copper concentrate), batch granulation tests were carried out using a rotating drum. The granule growth exhibited induction-type behaviour, as defined by Iveson and Litster [AIChE J. 44 (1998) 15 10]. There were two consecutive stages during granulation: the induction stage, during which the granules are gradually being compacted and little or no growth occurs, and the rapid growth stage, which starts when the granules have become surface wet and are rapidly growing. In agreement with earlier findings. an increased amount of binder liquid shortened the induction time. The compaction behaviour was also investigated. A displaced volume method was adopted to determine the porosity of the granules. It was shown that this technique had a limitation as it was unable to detect the reduction of the volumes of the granule pores after the granules had become surface wet. Due to this, some of the measurements were not suited for fitting a three-parameter empirical model. Attempts were made to determine whether the rapid growth stage started with the pore saturation exceeding a certain critical value, but due to the scatter in the porosity measurements and the fact that some of the measurements could not be used, it was not possible to determine a critical pore saturation, However, the porosity measurements clearly demonstrated that the porosity of the granules decreased during the induction stage of an experiment and that when rapid growth occurred, the granules had a pore saturation was around 0.85. This value was slightly lower than unity, which is most likely due to trapped air bubbles. (C) 2002 Published by Elsevier Science B.V.

Actively shielded multi-layer gradient coil designs with improved cooling properties

In standard cylindrical gradient coils consisting of a single layer of wires, a limiting factor in achieving very large magnetic field gradients is the rapid increase in coil resistance with efficiency. This is a particular problem in small-bore scanners, such as those used for MR microscopy. By adopting a multi-layer design in which the coil wires are allowed to spread out into multiple layers wound at increasing radii, a more favourable scaling of resistance with efficiency is achieved, thus allowing the design of more powerful gradient coils with acceptable resistance values. Previously this approach has been applied to the design of unshielded, longitudinal, and transverse gradient coils. Here, the multi-layer approach has been extended to allow the design of actively shielded multi-layer gradient coils, and also to produce coils exhibiting enhanced cooling characteristics. An iterative approach to modelling the steady-state temperature distribution within the coil has also been developed. Results indicate that a good level of screening can be achieved in multi-layer coils, that small versions of such coils can yield higher efficiencies at fixed resistance than conventional two-layer (primary and screen) coils, and that performance improves as the number of layers of increases. Simulations show that by optimising multi-layer coils for cooling it is possible to achieve significantly higher gradient strengths at a fixed maximum operating temperature. A four-layer coil of 8 mm inner diameter has been constructed and used to test the steady-state temperature model. (C) 2003 Elsevier Inc. All rights reserved.

Population pharmacokinetics and association between A77 1726 plasma concentrations and disease activity measures following administration of leflunomide to people with rheumatoid arthritis

Aims To investigate the concentration-effect relationship and pharmacokinetics of leflunomide in patients with rheumatoid arthritis (RA). Methods Data were collected from 23 RA patients on leflunomide therapy (as sole disease modifying antirheumatic drug (DMARD)) for at least 3 months. Main measures were A77 1726 (active metabolite of leflunomide) plasma concentrations and disease activity measures including pain, duration/intensity of morning stiffness, and SF-36 survey. A population estimate was sought for apparent clearance (CL/F ) and volume of distribution was fixed (0.155 l kg(-1)). Factors screened for influence on CL/F were weight, age, gender and estimated creatinine clearance. Results Significantly higher A77 1726 concentrations were seen in patients with less swollen joints and with higher SF-36 mental summary scores than in those with measures indicating more active disease (P < 0.05); concentration-effect trends were seen with five other disease activity measures. Statistical analysis of all disease activity measures showed that mean A77 1726 concentrations in groups with greater control of disease activity were significantly higher than those in whom the measures indicated less desirable control (P < 0.05). There was large between subject variability in the dose-concentration relationship. A steady-state infusion model best described the pharmacokinetic data. Inclusion of age as a covariate decreased interindividual variability (P < 0.01), but this would not be clinically important in terms of dosage changes. Final parameter estimate (% CV interindividual variability) for CL/F was 0.0184 l h(-1) (50%) (95% CI 0.0146, 0.0222). Residual (unexplained) variability (% CV) was 8.5%. Conclusions This study of leflunomide in patients using the drug clinically indicated a concentration-effect relationship. From our data, a plasma A77 1726 concentration of 50 mg l(-1) is more likely to indicate someone with less active disease than is a concentration around 30 mg l(-1). The marked variability in pharmacokinetics suggests a place for individualized dosing of leflunomide in RA therapy.