Computation of the linear elastic properties of random porous materials with a wide variety of microstructure

A finite-element method is used to study the elastic properties of random three-dimensional porous materials with highly interconnected pores. We show that Young's modulus, E, is practically independent of Poisson's ratio of the solid phase, nu(s), over the entire solid fraction range, and Poisson's ratio, nu, becomes independent of nu(s) as the percolation threshold is approached. We represent this behaviour of nu in a flow diagram. This interesting but approximate behaviour is very similar to the exactly known behaviour in two-dimensional porous materials. In addition, the behaviour of nu versus nu(s) appears to imply that information in the dilute porosity limit can affect behaviour in the percolation threshold limit. We summarize the finite-element results in terms of simple structure-property relations, instead of tables of data, to make it easier to apply the computational results. Without using accurate numerical computations, one is limited to various effective medium theories and rigorous approximations like bounds and expansions. The accuracy of these equations is unknown for general porous media. To verify a particular theory it is important to check that it predicts both isotropic elastic moduli, i.e. prediction of Young's modulus alone is necessary but not sufficient. The subtleties of Poisson's ratio behaviour actually provide a very effective method for showing differences between the theories and demonstrating their ranges of validity. We find that for moderate- to high-porosity materials, none of the analytical theories is accurate and, at present, numerical techniques must be relied upon.

The adult ventral nerve cord as a phylogenetic character in brachyceran Diptera

Insect ganglia are often composed of fused segmental units or neuromeres. We estimated the evolution of the ventral nerve cord (VNC) in higher Diptera by comparing the patterns of neuromere fusion among 33 families of the Brachycera. Variation within families is uncommon, and VNC architecture does not appear to be influenced by body shape. The outgroup pattern, seen in lower Diptera, is fusion of neuromeres belonging to thoracic segments 1 and 2 (T1 and T2), and fusion of neuromeres derived from T3 and abdominal segment 1 (A1). In the abdomen, neuromeres A7-10 are fused into the terminal abdominal ganglion (TAG). Increased neuromere fusion is a feature of the Brachycera. No brachyceran shows less fusion than the outgroups. We established six pattern elements; (1) fusion of T1 and T2, (2) fusion of T3 and A1, (3) fusion of the T1/T2 andT3/A1 ganglia, (4) increase in the number of neuromeres comprising the TAG, (5) anteriorward fusion of abdominal neuromeres, and (6) the complete fusion of thoracic and abdominal neuromeres into a synganglion. States 1 and 2 are present in the outgroup lower Diptera, and state 3 in the Xylophagomorpha, Stratiomyomorpha, Tabanomorpha and Cyclorrhapha. State 4 is a feature of all Eremoneura. State 5 is present in Cyclorrhapha only, and state 6, fusion into a synganglion, has evolved at least 4 times in the Eremoneura. Synapomorphies are provided for the Cyclorrhapha and Muscoidea, and a grouping of three basal brachyceran infraorders Xylophagomorpha, Stratiomyomorpha and Tabanomorpha. The patterns of fusion suggest that VNC architecture has evolved irreversibly, in accordance with Dollo's law.

Universal quantum computation and simulation using any entangling Hamiltonian and local unitaries

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? We provide an efficient algorithm to simulate any desired two-body Hamiltonian evolution using any fixed two-body entangling n-qubit Hamiltonian and local unitary operations. It follows that universal quantum computation can be performed using any entangling interaction and local unitary operations.

Practical scheme for quantum computation with any two-qubit entangling gate

Which gates are universal for quantum computation? Although it is well known that certain gates on two-level quantum systems (qubits), such as the controlled-NOT, are universal when assisted by arbitrary one-qubit gates, it has only recently become clear precisely what class of two-qubit gates is universal in this sense. We present an elementary proof that any entangling two-qubit gate is universal for quantum computation, when assisted by one-qubit gates. A proof of this result for systems of arbitrary finite dimension has been provided by Brylinski and Brylinski; however, their proof relies on a long argument using advanced mathematics. In contrast, our proof provides a simple constructive procedure which is close to optimal and experimentally practical.

Preparing encoded states in an oscillator

Recently a scheme has been proposed for constructing quantum error-correcting codes that embed a finite-dimensional code space in the infinite-dimensional Hilbert space of a system described by continuous quantum variables. One of the difficult steps in this scheme is the preparation of the encoded states. We show how these states can be generated by coupling a continuous quantum variable to a single qubit. An ion trap quantum computer provides a natural setting for a continuous system coupled to a qubit. We discuss how encoded states may be generated in an ion trap.

Simulation of insect movement with respect to plant architecture and morphogenesis

Developments in computer and three dimensional (3D) digitiser technologies have made it possible to keep track of the broad range of data required to simulate an insect moving around or over the highly heterogeneous habitat of a plant's surface. Properties of plant parts vary within a complex canopy architecture, and insect damage can induce further changes that affect an animal's movements, development and likelihood of survival. Models of plant architectural development based on Lindenmayer systems (L-systems) serve as dynamic platforms for simulation of insect movement, providing ail explicit model of the developing 3D structure of a plant as well as allowing physiological processes associated with plant growth and responses to damage to be described and Simulated. Simple examples of the use of the L-system formalism to model insect movement, operating Lit different spatial scales-from insects foraging on an individual plant to insects flying around plants in a field-are presented. Such models can be used to explore questions about the consequences of changes in environmental architecture and configuration on host finding, exploitation and its population consequences. In effect this model is a 'virtual ecosystem' laboratory to address local as well as landscape-level questions pertinent to plant-insect interactions, taking plant architecture into account. (C) 2002 Elsevier Science B.V. All rights reserved.

A methodology for building energy modelling and calibration in warm climates

In the last 7 years, a method has been developed to analyse building energy performance using computer simulation, in Brazil. The method combines analysis of building design plans and documentation, walk-through visits, electric and thermal measurements and the use of an energy simulation tool (DOE-2.1E code), The method was used to model more than 15 office buildings (more than 200 000 m(2)), located between 12.5degrees and 27.5degrees South latitude. The paper describes the basic methodology, with data for one building and presents additional results for other six cases. (C) 2002 Elsevier Science Ltd. All rights reserved.

Rapid field calculations for the effect of ferromagnetic material in MRI magnet design

Magnetic resonance imaging (MRI) magnets have very stringent constraints on the homogeneity of the static magnetic field that they generate over desired imaging regions. The magnet system also preferably generates very little stray field external to its structure, so that ease of siting and safety are assured. This work concentrates on deriving, means of rapidly computing the effect of 'cold' and 'warm' ferromagnetic material in or around the superconducting magnet system, so as to facilitate the automated design of hybrid material MR magnets. A complete scheme for the direct calculation of the spherical harmonics of the magnetic field generated by a circular ring of ferromagnetic material is derived under the conditions of arbitrary external magnetizing fields. The magnetic field produced by the superconducting coils in the system is computed using previously developed methods. The final, hybrid algorithm is fast enough for use in large-scale optimization methods. The resultant fields from a practical example of a 4 T, clinical MRI magnet containing both superconducting coils and magnetic material are presented.

An augmented Lanczos algorithm for the efficient computation of a dot-product of a function of a large sparse symmetric matrix

The Lanczos algorithm is appreciated in many situations due to its speed. and economy of storage. However, the advantage that the Lanczos basis vectors need not be kept is lost when the algorithm is used to compute the action of a matrix function on a vector. Either the basis vectors need to be kept, or the Lanczos process needs to be applied twice. In this study we describe an augmented Lanczos algorithm to compute a dot product relative to a function of a large sparse symmetric matrix, without keeping the basis vectors.

Stratigraphy, facies architecture and tectonic implications of the Upper Devonian to Lower Carboniferous Campwyn Volcanics of the northern New England Fold Belt

Upper Devonian to Lower Carboniferous strata of the Campwyn Volcanics of east central Queensland preserve a substantial sequence of first-cycle volcaniclastic sedimentary and coeval volcanic rocks that record prolonged volcanic activity along the northern New England Fold Belt. The style and scale of volcanism varied with time, producing an Upper Devonian sequence of mafic volcano-sedimentary rocks overlain by a rhyolitic ignimbrite-dominated sequence that passes upward into a Lower Carboniferous limestone-bearing sedimentary sequence. We define two facies associations for the Campwyn Volcanics. A lower facies association is dominated by mafic volcanic-derived sedimentary breccias with subordinate primary mafic volcanic rocks comprising predominantly hyaloclastite and peperite. Sedimentary breccias record episodic and high energy, subaqueous depositional events with clastic material sourced from a mafic lava-dominated terrain. Some breccias contain a high proportion of attenuated dense, glassy mafic juvenile clasts, suggesting a syn-eruptive origin. The lower facies association coarsens upwards from a lithic sand-dominated sequence through a thick interval of pebble- to boulder-grade polymict volcaniclastic breccias, culminating in facies that demonstrate subaerial exposure. The silicic upper facies association marks a significant change in eruptive style, magma composition and the nature of eruptive sources, as well as the widespread development of subaerial depositional conditions. Crystal-rich, high-grade, low- to high-silica rhyolite ignimbrites dominate the base of this facies association. Biostratigraphic age controls indicate that the ignimbrite-bearing sequences are Famennian to lower-mid Tournaisian in age. The ignimbrites represent extra-caldera facies with individual units up to 40 m thick and mostly lacking coarse lithic breccias. Thick deposits of pyroclastic material interbedded with fine-grained siliceous sandstone and mudstone (locally radiolarian-bearing) were deposited from pyroclastic flows that crossed palaeoshorelines or represent syn-eruptive, resedimented pyroclastic material. Some block-bearing lithic-pumice-crystal breccias may also reflect more proximal subaqueous silicic explosive eruptions. Crystal-lithic sandstones interbedded with, and overlying the ignimbrites, contain abundant detrital volcanic quartz and feldspar derived from the pyroclastic deposits. Limestone is common in the upper part of the upper facies association, and several beds are oolitic (cf. Rockhampton Group of the Yarrol terrane). Overall, the upper facies association fines upward and is transgressive, recording a return to shallow-marine conditions. Palaeocurrent data from all stratigraphic levels in the Campwyn Volcanics indicate that the regional sediment-dispersal direction was to the northwest, and opposed to the generally accepted notion of easterly sediment dispersal from a volcanic arc source. The silicic upper facies association correlates in age and lithology to Early Carboniferous silicic volcanism in the Drummond (Cycle 1) and Burdekin Basins, Connors Arch, and in the Yarrol terranes of eastern Queensland. The widespread development of silicic volcanism in the Early Carboniferous indicates that silicic (rift-related) magmatism was not restricted to the Drummond Basin, but was part of a more substantial silicic igneous province.

Heritability of adult body height: A comparative study of twin cohorts in eight countries

A major component of variation in body height is due to genetic differences, but environmental factors have a substantial contributory effect. In this study we aimed to analyse whether the genetic architecture of body height varies between affluent western societies. We analysed twin data from eight countries comprising 30,111 complete twin pairs by using the univariate genetic model of the Mx statistical package. Body height and zygosity were self-reported in seven populations and measured directly in one population. We found that there was substantial variation in mean body height between countries; body height was least in Italy (177 cm in men and 163 cm in women) and greatest in the Netherlands (184 cm and 171 cm, respectively). In men there was no corresponding variation in heritability of body height, heritability estimates ranging from 0.87 to 0.93 in populations under an additive genes/unique environment (AE) model. Among women the heritability estimates were generally lower than among men with greater variation between countries, ranging from 0.68 to 0.84 when an additive genes/shared environment/unique environment (ACE) model was used. In four populations where an AE model fit equally well or better, heritability ranged from 0.89 to 0.93. This difference between the sexes was mainly due to the effect of the shared environmental component of variance, which appears to be more important among women than among men in our study populations. Our results indicate that, in general, there are only minor differences in the genetic architecture of height between affluent Caucasian populations, especially among men.

Numerical study of hydrogenic effective mass theory for an impurity P donor in Si in the presence of an electric field and interfaces

In this paper we examine the effects of varying several experimental parameters in the Kane quantum computer architecture: A-gate voltage, the qubit depth below the silicon oxide barrier, and the back gate depth to explore how these variables affect the electron density of the donor electron. In particular, we calculate the resonance frequency of the donor nuclei as a function of these parameters. To do this we calculated the donor electron wave function variationally using an effective-mass Hamiltonian approach, using a basis of deformed hydrogenic orbitals. This approach was then extended to include the electric-field Hamiltonian and the silicon host geometry. We found that the phosphorous donor electron wave function was very sensitive to all the experimental variables studied in our work, and thus to optimize the operation of these devices it is necessary to control all parameters varied in this paper.

Genetic variation in foliar carbon isotope composition in relation to tree growth and foliar nitrogen concentration in clones of the F-1 hybrid between slash pine and Caribbean pine

The objectives of this study were: (1) to quantify the genetic variation in foliar carbon isotope composition (delta(13)C) of 122 clones of ca. 4-year-old F-1 hybrids between slash pine (Pinus elliottii Engelm var. elliottii) and Caribbean pine (Pinus caribaea var. hondurensis Barr.,et Golf.) grown at two field experimental sites with different water and nitrogen availability in southeast Queensland, Australia, in relation to tree growth and foliar nitrogen concentration (N-mass); and (2) to assess the potential of using delta(13)C measurements, in the foliage materials collected from the clone hedges at nursery and the 4-year-old tree canopies in the field, as an indirect index of tree water use efficiency for selecting elite F-1 hybrid pine clones with improved tree growth. There were significant differences in foliar delta(13)C between the nursery hedges and the 4-year-old tree canopies in the field, between the summer and winter seasons, between the two experimental sites, and between the upper outer and lower outer canopy positions sampled. This indicates that delta(13)C measurements in the foliage materials are significantly influenced by the sampling techniques and environmental conditions. Significant differences in foliar delta(13)C, at the upper outer canopy in both field experiments in summer and winter, were detected between the clones, and between the female parents of the clones. Clone means of tree height at age ca. 3 years were positively related to those of the upper outer canopy delta(13)C at both experimental sites in winter, but only for the wetter site in summer. There were positive, linear relationships between clone means of canopy delta(13)C and those of canopy N-mass, indicating that canopy photosynthetic capacity might be an important factor regulating the clonal variation in canopy delta(13)C. Significant correlations were found between clone means of canopy delta(13)C at both experimental sites in summer and winter, and between those at the upper outer and lower outer canopy positions. Mean clone delta(13)C for the nursery hedges was only positively related to mean clone stem diameter at 1.3 m height at age 3 years on the wetter site. The clone by site interaction for foliar delta(13)C at the upper outer canopy was significant only in summer. Overall, the relatively high genetic variance components for foliar delta(13)C and significant, positive correlations between clone means of foliar delta(13)C and tree growth have highlighted the potential of using foliar delta(13)C measurements for assisting in selection of the elite F-1 hybrid pine clones with improved tree growth. (C) 2002 Elsevier Science B.V. All rights reserved.

A yield architecture framework to explain adaptation of pearl millet to environmental stress

Functional knowledge of the physiological basis of crop adaptation to stress is a prerequisite for exploiting specific adaptation to stress environments in breeding programs. This paper presents an analysis of yield components for pearl millet, to explain the specific adaptation of local landraces to stress environments in Rajasthan, India. Six genotypes, ranging from high-tillering traditional landraces to low-tillering open-pollinated modern cultivars, were grown in 20 experiments, covering a range of nonstress and drought stress patterns. In each experiment, yield components (particle number, grain number, 100 grain mass) were measured separately for main shoots, basal tillers, and nodal tillers. Under optimum conditions, landraces had a significantly lower grain yield than the cultivars, but no significant differences were observed at yield levels around 1 ton ha(-1). This genotype x environment interaction for grain yield was due to a difference in yield strategy, where landraces aimed at minimising the risk of a crop failure under stress conditions, and modem cultivars aimed at maximising yield potential under optimum conditions. A key aspect of the adaptation of landraces was the small size of the main shoot panicle, as it minimised (1) the loss of productive tillers during stem elongation; (2) the delay in anthesis if mid-season drought occurs; and (3) the reduction in panicle productivity of the basal tillers under stress. In addition, a low investment in structural panicle weight, relative to vegetative crop growth rate, promoted the production of nodal tillers, providing a mechanism to compensate for reduced basal tiller productivity if stress occurred around anthesis. A low maximum 100 grain mass also ensured individual grain mass was little affected by environmental conditions. The strategy of the high-tillering landraces carries a yield penalty under optimum conditions, but is expected to minimise the risk of a crop failure, particularly if mid-season drought stress occurs. The yield architecture of low-tillering varieties, by contrast, will be suited to end-of-season drought stress, provided anthesis is early. Application of the above adaptation mechanisms into a breeding program could enable the identification of plant types that match the prevalent stress patterns in the target environments. (C) 2003 E.J. van Oosterom. Published by Elsevier Science B.V. All rights reserved.

Contrasts in the basins of attraction of structurally identical iterative root finding methods

A numerical comparison is performed between three methods of third order with the same structure, namely BSC, Halley’s and Euler–Chebyshev’s methods. As the behavior of an iterative method applied to a nonlinear equation can be highly sensitive to the starting points, the numerical comparison is carried out, allowing for complex starting points and for complex roots, on the basins of attraction in the complex plane. Several examples of algebraic and transcendental equations are presented.