Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Generalized quasi mode theory of macroscopic canonical quantization in cavity quantum electrodynamics and quantum optics I. Theory

The quasi mode theory of macroscopic quantization in quantum optics and cavity QED developed by Dalton, Barnett and Knight is generalized. This generalization allows for cases in which two or more quasi permittivities, along with their associated mode functions, are needed to describe the classical optics device. It brings problems such as reflection and refraction at a dielectric boundary, the linear coupler, and the coupling of two optical cavities within the scope of the theory. For the most part, the results that are obtained here are simple generalizations of those obtained in previous work. However the coupling constants, which are of great importance in applications of the theory, are shown to contain significant additional terms which cannot be 'guessed' from the simpler forms. The expressions for the coupling constants suggest that the critical factor in determining the strength of coupling between a pair of quasi modes is their degree of spatial overlap. In an accompanying paper a fully quantum theoretic derivation of the laws of reflection and refraction at a boundary is given as an illustration of the generalized theory. The quasi mode picture of this process involves the annihilation of a photon travelling in the incident region quasi mode, and the subsequent creation of a photon in either the incident region or transmitted region quasi modes.

Generalized quasi mode theory of macroscopic canonical quantization in cavity quantum electrodynamics and quantum optics II. Application to reflection and refraction at a dielectric interface

The generalization of the quasi mode theory of macroscopic quantization in quantum optics and cavity QED presented in the previous paper, is applied to provide a fully quantum theoretic derivation of the laws of reflection and refraction at a boundary. The quasi mode picture of this process involves the annihilation of a photon travelling in the incident region quasi mode, and the subsequent creation of a photon in either the incident region or transmitted region quasi modes. The derivation of the laws of reflection and refraction is achieved through the dual application of the quasi mode theory and a quantum scattering theory based on the Heisenberg picture. Formal expressions from scattering theory are given for the reflection and transmission coefficients. The behaviour of the intensity for a localized one photon wave packet coming in at time minus infinity from the incident direction is examined and it is shown that at time plus infinity, the light intensity is only significant where the classical laws of reflection and refraction predict. The occurrence of both refraction and reflection is dependent upon the quasi mode theory coupling constants between incident and transmitted region quasi modes being nonzero, and it is seen that the contributions to such coupling constants come from the overlap of the mode functions in the boundary layer region, as might be expected from a microscopic theory.

Complex natural resonances of conducting planar objects buried in a dielectric half-space

A scheme is presented to incorporate a mixed potential integral equation (MPIE) using Michalski's formulation C with the method of moments (MoM) for analyzing the scattering of a plane wave from conducting planar objects buried in a dielectric half-space. The robust complex image method with a two-level approximation is used for the calculation of the Green's functions for the half-space. To further speed up the computation, an interpolation technique for filling the matrix is employed. While the induced current distributions on the object's surface are obtained in the frequency domain, the corresponding time domain responses are calculated via the inverse fast Fourier transform (FFT), The complex natural resonances of targets are then extracted from the late time response using the generalized pencil-of-function (GPOF) method. We investigate the pole trajectories as we vary the distance between strips and the depth and orientation of single, buried strips, The variation from the pole position of a single strip in a homogeneous dielectric medium was only a few percent for most of these parameter variations.

Do generalized visual training programmes for sport really work? An experimental investigation

We assessed the effectiveness of two generalized visual training programmes in enhancing visual and motor performance for racquet sports. Forty young participants were assigned equally to groups undertaking visual training using Revien and Gabor's Sports Vision programme (Group 1), visual training using Revien's Eyerobics (Group 2), a placebo condition involving reading (Group 3) and a control condition involving physical practice only (Group 4). Measures of basic visual function and of sport-specific motor performance were obtained from all participants before and immediately after a 4-week training period. Significant pre- to post-training differences were evident on some of the measures; however, these were not group-dependent. Contrary to the claims made by proponents of generalized visual training, we found no evidence that the visual training programmes led to improvements in either vision or motor performance above and beyond those resulting simply from test familiarity.

Two-link approximation schemes for loss networks with linear structure and trunk reservation

Loss networks have long been used to model various types of telecommunication network, including circuit-switched networks. Such networks often use admission controls, such as trunk reservation, to optimize revenue or stabilize the behaviour of the network. Unfortunately, an exact analysis of such networks is not usually possible, and reduced-load approximations such as the Erlang Fixed Point (EFP) approximation have been widely used. The performance of these approximations is typically very good for networks without controls, under several regimes. There is evidence, however, that in networks with controls, these approximations will in general perform less well. We propose an extension to the EFP approximation that gives marked improvement for a simple ring-shaped network with trunk reservation. It is based on the idea of considering pairs of links together, thus making greater allowance for dependencies between neighbouring links than does the EFP approximation, which only considers links in isolation.

Field quantization, photons and non-Hermitean modes

Field quantization in unstable optical systems is treated by expanding the vector potential in terms of non-Hermitean (Fox-Li) modes. We define non-Hermitean modes and their adjoints in both the cavity and external regions and make use of the important bi-orthogonality relationships that exist within each mode set. We employ a standard canonical quantization procedure involving the introduction of generalized coordinates and momenta for the electromagnetic (EM) field. Three-dimensional systems are treated, making use of the paraxial and monochromaticity approximations for the cavity non-Hermitean modes. We show that the quantum EM field is equivalent to a set of quantum harmonic oscillators (QHOs), associated with either the cavity or the external region non-Hermitean modes, and thus confirming the validity of the photon model in unstable optical systems. Unlike in the conventional (Hermitean mode) case, the annihilation and creation operators we define for each QHO are not Hermitean adjoints. It is shown that the quantum Hamiltonian for the EM field is the sum of non-commuting cavity and external region contributions, each of which can be expressed as a sum of independent QHO Hamiltonians for each non-Hermitean mode, except that the external field Hamiltonian also includes a coupling term responsible for external non-Hermitean mode photon exchange processes. The non-commutativity of certain cavity and external region annihilation and creation operators is associated with cavity energy gain and loss processes, and may be described in terms of surface integrals involving cavity and external region non-Hermitean mode functions on the cavity-external region boundary. Using the essential states approach and the rotating wave approximation, our results are applied to the spontaneous decay of a two-level atom inside an unstable cavity. We find that atomic transitions leading to cavity non-Hermitean mode photon absorption are associated with a different coupling constant to that for transitions leading to photon emission, a feature consequent on the use of non-Hermitean mode functions. We show that under certain conditions the spontaneous decay rate is enhanced by the Petermann factor.

On the fast approximation of Green's functions in MPIE formulations for planar layered media

The numerical implementation of the complex image approach for the Green's function of a mixed-potential integralequation formulation is examined and is found to be limited to low values of k(0) rho (in this context k(0) rho = 2 pirho/ lambda(0), where rho is the distance between the source and the field points of the Green's function and lambda(0) is the free space wavelength). This is a clear limitation for problems of large dimension or high frequency where this limit is easily exceeded. This paper examines the various strategies and proposes a hybrid method whereby most of the above problems can be avoided. An efficient integral method that is valid for large k(0) rho is combined with the complex image method in order to take advantage of the relative merits of both schemes. It is found that a wide overlapping region exists between the two techniques allowing a very efficient and consistent approach for accurately calculating the Green's functions. In this paper, the method developed for the computation of the Green's function is used for planar structures containing both lossless and lossy media.

Micro-mechanical approximation to the stress field around an inclined fiber of FRCC

Matrix spalling or crushing is one of the important mechanisms of fiber-matrix interaction of fiber reinforced cementitious composites (FRCC). The fiber pullout mechanisms have been extensively studied for an aligned fiber but matrix failure is rarely investigated since it is thought not to be a major affect. However, for an inclined fiber, the matrix failure should not be neglected. Due to the complex process of matrix spalling, experimental investigation and analytical study of this mechanism are rarely found in literature. In this paper, it is assumed that the load transfer is concentrated within the short length of the inclined fiber from the exit point towards anchored end and follows the exponential law. The Mindlin formulation is employed to calculate the 3D stress field. The simulation gives much information about this field. The 3D approximation of the stress state around an inclined fiber helps to qualitatively understand the mechanism of matrix failure. Finally, a spalling criterion is proposed by which matrix spalling occurs only when the stress in a certain volume, rather than the stress at a small point, exceeds the material strength. This implies some local stress redistribution after first yield. The stress redistribution results in more energy input and higher load bearing capacity of the matrix. In accordance with this hypothesis, the evolution of matrix spalling is demonstrated. The accurate prediction of matrix spalling needs the careful determination of the parameters in this model. This is the work of further study. (C) 2002 Elsevier Science Ltd. All rights reserved.

Modelling activated sludge flocculation using population balances

Activated sludge flocculation was modelled using population balances. The model followed the dynamics of activated sludge flocculation providing a good approximation of the change in mean floe size with time. Increasing the average velocity gradient decreased the final floe size. The breakage rate coefficient and collision efficiency also varied with the average velocity gradient. A power law relationship was found for the increase in breakage rate coefficient with increasing average velocity gradient. Further investigation will be conducted to determine the relationship between the collision efficiency and particle size to provide a better approximation of dynamic changes in the floe size distribution during flocculation. (C) 2002 Published by Elsevier Science B.V.

Testing the link between the latitudinal gradient in species richness and rates of molecular evolution

Numerous hypotheses have been proposed to explain latitudinal gradients in species richness, but all are subject to ongoing debate. Here we examine Rohde's (1978, 1992) hypothesis, which proposes that climatic conditions at low latitudes lead to elevated rates of speciation. This hypothesis predicts that rates of molecular evolution should increase towards lower latitudes, but this prediction has never been tested. We discuss potential links between rates of molecular evolution and latitudinal diversity gradients, and present the first test of latitudinal variation in rates of molecular evolution. Using 45 phylogenetically independent, latitudinally separated pairs of bird species and higher taxa, we compare rates of evolution of two mitochondrial genes and DNA-DNA hybridization distances. We find no support for an effect of latitude on rate of molecular evolution. This result casts doubt on the generality of a key component of Rohde's hypothesis linking climate and speciation.

Non-compact generalized variational inequalities for quasi-monotone and hemi-continuous operators with applications

Some results are obtained for non-compact cases in topological vector spaces for the existence problem of solutions for some set-valued variational inequalities with quasi-monotone and lower hemi-continuous operators, and with quasi-semi-monotone and upper hemi-continuous operators. Some applications are given in non-reflexive Banach spaces for these existence problems of solutions and for perturbation problems for these set-valued variational inequalities with quasi-monotone and quasi-semi-monotone operators.