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.

Derivation of the probability distribution function for the local density of states of a disordered quantum wire via the replica trick and supersymmetry

We consider the statistical properties of the local density of states of a one-dimensional Dirac equation in the presence of various types of disorder with Gaussian white-noise distribution. It is shown how either the replica trick or supersymmetry can be used to calculate exactly all the moments of the local density of states.' Careful attention is paid to how the results change if the local density of states is averaged over atomic length scales. For both the replica trick and supersymmetry the problem is reduced to finding the ground state of a zero-dimensional Hamiltonian which is written solely in terms of a pair of coupled spins which are elements of u(1, 1). This ground state is explicitly found for the particular case of the Dirac equation corresponding to an infinite metallic quantum wire with a single conduction channel. The calculated moments of the local density of states agree with those found previously by Al'tshuler and Prigodin [Sov. Phys. JETP 68 (1989) 198] using a technique based on recursion relations for Feynman diagrams. (C) 2001 Elsevier Science B.V. All rights reserved.

Measurement of qubits

We describe in detail the theory underpinning the measurement of density matrices of a pair of quantum two-level systems (qubits). Our particular emphasis is on qubits realized by the two polarization degrees of freedom of a pair of entangled photons generated in a down-conversion experiment; however, the discussion applies in general, regardless of the actual physical realization. Two techniques are discussed, namely, a tomographic reconstruction (in which the density matrix is linearly related to a set of measured quantities) and a maximum likelihood technique which requires numerical optimization (but has the advantage of producing density matrices that are always non-negative definite). In addition, a detailed error analysis is presented, allowing errors in quantities derived from the density matrix, such as the entropy or entanglement of formation, to be estimated. Examples based on down-conversion experiments are used to illustrate our results.

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.

Influence of seaward boundary condition on contaminant transport in unconfined coastal aquifers

Contaminant transport in coastal aquifers is complicated partly due to the conditions at the seaward boundary including seawater intrusion and tidal variations of sea level. Their inclusion in modelling this system will be computationally expensive. Therefore, it will be instructive to investigate the consequence of simplifying the seaward boundary condition by neglecting the seawater density and tidal variations in numerical predictions of contaminant transport in this zone. This paper presents a comparison of numerical predictions for a simplified seaward boundary condition with experimental results for a corresponding realistic one including a saltwater interface and tidal variations. Different densities for contaminants are considered. The comparison suggests that the neglect of the seawater intrusion and tidal variations does not affect noticeably the overall migration rate of the plume before it reaches the saltwater interface. However, numerical prediction shows that a more dense contaminant travels further seaward and part of the solute mass exits under the sea if the seawater density is not included. This is not consistent with the experimental result, which shows that the contaminant travels upwards to the shoreline along the saltwater interface. Neglect of seawater density, therefore, will result in an underestimation of the exit rate of solute mass around the coastline and fictitious migration paths under the seabed. For a less dense contaminant, neglect of seawater density has little effect on numerical prediction of migration paths. (C) 2001 Elsevier Science B.V. All rights reserved.

'Neighbourhood' size, dispersal and density estimates in the prickly forest skink (Gnypetoscincus queenslandiae) using individual genetic and demographic methods

Dispersal, or the amount of dispersion between an individual's birthplace and that of its offspring, is of great importance in population biology, behavioural ecology and conservation, however, obtaining direct estimates from field data on natural populations can be problematic. The prickly forest skink, Gnypetoscincus queenslandiae, is a rainforest endemic skink from the wet tropics of Australia. Because of its log-dwelling habits and lack of definite nesting sites, a demographic estimate of dispersal distance is difficult to obtain. Neighbourhood size, defined as 4 piD sigma (2) (where D is the population density and sigma (2) the mean axial squared parent-offspring dispersal rate), dispersal and density were estimated directly and indirectly for this species using mark-recapture and microsatellite data, respectively, on lizards captured at a local geographical scale of 3 ha. Mark-recapture data gave a dispersal rate of 843 m(2)/generation (assuming a generation time of 6.5 years), a time-scaled density of 13 635 individuals * generation/km(2) and, hence, a neighbourhood size of 144 individuals. A genetic method based on the multilocus (10 loci) microsatellite genotypes of individuals and their geographical location indicated that there is a significant isolation by distance pattern, and gave a neighbourhood size of 69 individuals, with a 95% confidence interval between 48 and 184. This translates into a dispersal rate of 404 m(2)/generation when using the mark-recapture density estimation, or an estimate of time-scaled population density of 6520 individuals * generation/km(2) when using the mark-recapture dispersal rate estimate. The relationship between the two categories of neighbourhood size, dispersal and density estimates and reasons for any disparities are discussed.

The effect of channel geometry and wall boundary conditions on the formation of extrusion surface instabilities for LLDPE

It is believed that surface instabilities can occur during the extrusion of linear low density polyethylene due to high extensional stresses at the exit of the die. Local crack development can occur at a critical stress level when melt rupture is reached. This high extensional stress results from the rearrangement of the flow at the boundary transition between the wall exit and the free surface. The stress is highest at the extrudate surface and decreases into the bulk of the material. The location of the region where the critical level is reached can determine the amplitude of the extrudate surface distortion, This paper studies the effect of wall slip on the numerically simulated extensional stress level at the die exit and correlates this to the experimentally determined amplitude of the surface instability. The effect of die exit radius and die wall roughness on extrusion surface instabilities is also correlated to the exit stress level in the same way. Whereas full slip may completely suppress the surface instability, a reduction in the exit stress level and instability amplitude is also shown for a rounded die exit and a slight increase in instability is shown to result from a rough die wall. A surface instability map demonstrates how the shear rate for onset of extrusion surface instabilities can be predicted on the basis of melt strength measurements and simulated stress peaks at the exit of the die. (C) 2001 Elsevier Science B.V. All rights reserved.

Maximum Likelihood Estimation of Mixture Densities for Binned and Truncated Multivariate Data

Binning and truncation of data are common in data analysis and machine learning. This paper addresses the problem of fitting mixture densities to multivariate binned and truncated data. The EM approach proposed by McLachlan and Jones (Biometrics, 44: 2, 571-578, 1988) for the univariate case is generalized to multivariate measurements. The multivariate solution requires the evaluation of multidimensional integrals over each bin at each iteration of the EM procedure. Naive implementation of the procedure can lead to computationally inefficient results. To reduce the computational cost a number of straightforward numerical techniques are proposed. Results on simulated data indicate that the proposed methods can achieve significant computational gains with no loss in the accuracy of the final parameter estimates. Furthermore, experimental results suggest that with a sufficient number of bins and data points it is possible to estimate the true underlying density almost as well as if the data were not binned. The paper concludes with a brief description of an application of this approach to diagnosis of iron deficiency anemia, in the context of binned and truncated bivariate measurements of volume and hemoglobin concentration from an individual's red blood cells.

Low-density koala (Phascolarctos cinereus) populations in the mulgalands of south-west Queensland. I. Faecal pellet sampling protocol

Koala (Phascolarctos cinereus) populations in eastern Australia are threatened by land clearing for agricultural and urban development. At the same time, conservation efforts are hindered by a dearth of information about inland populations. Faecal deposits offer a source of information that is readily available and easily collected non-invasively. We detail a faecal pellet sampling protocol that was developed for use in a large rangeland biogeographic region. The method samples trees in belt transects, uses a thorough search at the tree base to quickly identify trees with koala pellets under them, then estimates the abundance of faecal pellets under those trees using 1-m(2) quadrats. There was a strong linear relationship between these estimates and a complete enumeration of pellet abundance under the same trees. We evaluated the accuracy of our method in detecting trees where pellets were present by means of a misclassification index that was weighed more heavily for missed trees that had high numbers of pellets under them. This showed acceptable accuracy in all landforms except riverine, where some trees with large numbers of pellets were missed. Here, accuracy in detecting pellet presence was improved by sampling with quadrats, rather than basal searches. Finally, we developed a method to reliably age pellets and demonstrate how this protocol could be used with the faecal-standing-crop method to derive a regional estimate of absolute koala abundance.

Wolbachia density and virulence attenuation after transfer into a novel host

The factors that control replication rate of the intracellular bacterium Wolbachia pipientis in its insect hosts are unknown and difficult to explore, given the complex interaction of symbiont and host genotypes. Using a strain of Wolbachia that is known to over-replicate and shorten the lifespan of its Drosophila melanogaster host, we have tracked the evolution of replication control in both somatic and reproductive tissues in a novel host/Wolbachia association. After transinfection (the transfer of a Wolbachia strain into a different species) of the over-replicating Wolbachia popcorn strain from D. metanogaster to Drosophila simulans, we demonstrated that initial high densities in the ovaries were in excess of what was required for perfect maternal transmission, and were likely causing reductions in reproductive fitness. Both densities and fitness costs associated with ovary infection rapidly declined in the generations after transinfection. The early death effect in D. simulans attenuated only slightly and was comparable to that induced in D. metanogaster. This study reveals a strong host involvement in Wolbachia replication rates, the independence of density control responses in different tissues, and the strength of natural selection acting on reproductive fitness.

Numerical investigations of flow and energy fields near a thermoacoustic couple

The flow field and the energy transport near thermoacoustic couples are simulated using a 2D full Navier-Stokes solver. The thermoacoustic couple plate is maintained at a constant temperature; plate lengths, which are short and long compared with the particle displacement lengths of the acoustic standing waves, are tested. Also investigated are the effects of plate spacing and the amplitude of the standing wave. Results are examined in the form of energy vectors, particle paths, and overall entropy generation rates. These show that a net heat-pumping effect appears only near the edges of thermoacoustic couple plates, within about a particle displacement distance from the ends. A heat-pumping effect can be seen even on the shortest plates tested when the plate spacing exceeds the thermal penetration depth. It is observed that energy dissipation near the plate increases quadratically as the plate spacing is reduced. The results also indicate that there may be a larger scale vortical motion outside the plates which disappears as the plate spacing is reduced. (C) 2002 Acoustical Society of America.

The relationship of natural killer cell counts, perforin mRNA and CD2 expression to post-exercise natural killer cell activity in humans

The aim of this study was to further investigate the mechanism of suppression of natural killer (NK) cell cytotoxic activity In peripheral blood following strenuous exercise. Blood was collected for analysis of NK cell concentration, cytotoxic activity, CD2 surface expression and perforin gene expression from runners (RUN, n = 6) and resting controls (CONTROL, n = 4) pre-exercise, 0, 1.5, 5, and 24 h following a 60-min treadmill run at 80% of VO2 peak. Natural killer cytotoxic activity, measured using a whole blood chromium release assay, fluctuated minimally in the CONTROL group and increased by 63% and decreased by 43% 0 and 1.5 h post-exercise, respectively, in the RUN group (group x time, P < 0.001). Lytic index (cytotoxic activity per cell) did not change. Perforin mRNA, measured using quantitative real-time polymerase chain reaction (ORT-PCR) decreased from pre- to post-exercise and remained decreased through 24 h, The decrease from pre- to 0 In post-exercise was seen predominately in the RUN group and was inversely correlated r = - 0.95) to pre-exercise perform mRNA. The NK cell surface expression of CD2 (lymphocyte function-associated antigen-2) was determined using fluorescent antibodies and flow cytometry, There was no change in the proportion of NK cells expressing CD2 or CD2 density, We conclude that (1) numerical redistribution accounted for most of the change in NK cytotoxic activity following a strenuous run, (2) decrease in perforin gene expression during the run was inversely related to pre-exercise levels but did not parallel changes in cytotoxic activity, and (3) CD2 surface expression was not affected by exercise.