A hybrid, inverse approach to the design of magnetic resonance imaging magnets

This paper describes a hybrid numerical method of an inverse approach to the design of compact magnetic resonance imaging magnets. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first, kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. The emphasis of this work is on the optimal design of short MRI magnets. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric MRI magnets as well as asymmetric magnets. The results highlight that the method can be used to obtain a compact MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1 m in length, significantly shorter than current designs. Viable asymmetric magnet designs, in which the edge of the homogeneous region is very close to one end of the magnet system are also presented. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. (C) 2000 American Association of Physicists in Medicine. [S0094-2405(00)00303-5].

Effect of density on spatial distribution, pit formation and pit diameter of Myrmeleon acer Walker, (Neuroptera : Myrmeleontidae): patterns and processes

The effect of increasing population density on the formation of pits, their size and spatial distribution, and on levels of mortality was examined in the antlion Myrmeleon acer Walker. Antlions were kept at densities ranging from 0.4 to 12.8 individuals per 100 cm(2). The distribution of pits was regular or uniform across all densities, but antlions constructed proportionally fewer and smaller pits as density increased. Mortality through cannibalism was very low and only occurred at densities greater than five individuals per 100 cm(2). Antlions in artificially crowded situations frequently relocated their pits and when more space became available, individuals became more dispersed with time. Redistribution of this species results from active avoidance of other antlions and sand throwing associated with pit construction and maintenance, rather than any attempt to optimise prey capture per se.

Comparative analysis of dendritic cell density and total number in commonly transplanted organs: morphometric estimation in normal mice

Dendritic cells (DC) are considered to be the major cell type responsible for induction of primary immune responses. While they have been shown to play a critical role in eliciting allosensitization via the direct pathway, there is evidence that maturational and/or activational heterogeneity between DC in different donor organs may be crucial to allograft outcome. Despite such an important perceived role for DC, no accurate estimates of their number in commonly transplanted organs have been reported. Therefore, leukocytes and DC were visualized and enumerated in cryostat sections of normal mouse (C57BL/10, B10.BR, C3H) liver, heart, kidney and pancreas by immunohistochemistry (CD45 and MHC class II staining, respectively). Total immunopositive cell number and MHC class II+ cell density (C57BL/10 mice only) were estimated using established morphometric techniques - the fractionator and disector principles, respectively. Liver contained considerably more leukocytes (similar to 5-20 x 10(6)) and DC (similar to 1-3 x 10(6)) than the other organs examined (pancreas: similar to 0.6 x 10(6) and similar to 0.35 x 10(6): heart: similar to 0.8 x 10(6) and similar to 0.4 x 10(6); kidney similar to 1.2 x 10(6) and 0.65 x 10(6), respectively). In liver, DC comprised a lower proportion of all leukocytes (similar to 15-25%) than in the other parenchymal organs examined (similar to 40-60%). Comparatively, DC density in C57BL/10 mice was heart > kidney > pancreas much greater than liver (similar to 6.6 x 10(6), 5 x 10(6), 4.5 x 10(6) and 1.1 x 10(6) cells/cm(3), respectively). When compared to previously published data on allograft survival, the results indicate that the absolute number of MHC class II+ DC present in a donor organ is a poor predictor of graft outcome. Survival of solid organ allografts is more closely related to the density of the donor DC network within the graft. (C) 2000 Elsevier Science B.V. All rights reserved.

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.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Race and sex effects on the association between muscle strength, soft tissue, and bone mineral density in healthy elders: The Health, Aging, and Body Composition Study

Two factors generally reported to influence bone density are body composition and muscle strength. However, it is unclear if these relationships are consistent across race and sex, especially in older persons. If differences do exist by race and/or sex, then strategies to maintain bone mass or minimize bone loss in older adults may need to be modified accordingly. Therefore, we examined the independent effects of bone mineral-free lean mass (LM), fat mass (FM), and muscle strength on regional and whole body bone mineral density (BMD) in a cohort of 2619 well-functioning older adults participating in the Health, Aging, and Body Composition (Health ABC) Study with complete measures. Participants included 738 white women, 599 black women, 827 white men, and 455 black men aged 70-79 years. BMD (g/cm(2)) of the femoral neck, whole body, upper and lower limb, and whole body and upper limb bone mineral-free LM and FM was assessed by dual-energy X-ray absorptiometry (DXA). Handgrip strength and knee extensor torque were determined by dynamometry. In analyses stratified by race and sex and adjusted for a number of confounders, LM was a significant (p < 0.001) determinant of BMD, except in white women for the lower limb and whole body. In women, FM also was an independent contributor to BMD at the femoral neck, and both PM and muscle strength contributed to limb BMD. The following were the respective Beta-weights (regression coefficients for standardized data, Std beta) and percent difference in BMD per unit (7.5 kg) LM: femoral neck, 0.202-0.386 and 4.7-6.9 %; lower limb,.0.209-0.357 and 2.9-3.5%; whole body, 0.239-0.484 and 3.0-4.7 %; and upper limb (unit = 0.5 kg), 0.231-0.407 and 3.1-3.4%. Adjusting for bone size (bone mineral apparent density [BMAD]) or body size BMD/height) diminished the importance of LM, and the contributory effect of FM became more pronounced. These results indicate that LM and FM were associated with bone mineral depending on the bone site and bone index used. Where differences did occur, they were primarily by sex not race. To preserve BMD, maintaining or increasing LM in the elderly would appear to be an appropriate strategy, regardless of race or sex.

A numerical study of flexural buckling of foliated rock slopes

The occurrence of foliated rock masses is common in mining environment. Methods employing continuum approximation in describing the deformation of such rock masses possess a clear advantage over methods where each rock layer and each inter-layer interface (joint) is explicitly modelled. In devising such a continuum model it is imperative that moment (couple) stresses and internal rotations associated with the bending of the rock layers be properly incorporated in the model formulation. Such an approach will lead to a Cosserat-type theory. In the present model, the behaviour of the intact rock layer is assumed to be linearly elastic and the joints are assumed to be elastic-perfectly plastic. Condition of slip at the interfaces are determined by a Mohr-Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformations. The model is incorporated into the finite element program AFENA and validated against an analytical solution of elementary buckling problems of a layered medium under gravity loading. A design chart suitable for assessing the stability of slopes in foliated rock masses against flexural buckling failure has been developed. The design chart is easy to use and provides a quick estimate of critical loading factors for slopes in foliated rock masses. It is shown that the model based on Euler's buckling theory as proposed by Cavers (Rock Mechanics and Rock Engineering 1981; 14:87-104) substantially overestimates the critical heights for a vertical slope and underestimates the same for sub-vertical slopes. Copyright (C) 2001 John Wiley & Sons, Ltd.

Finite element modelling of heat transfer through permeable cracks in hydrothermal systems with upward throughflow

We use the finite element method to model the heat transfer phenomenon through permeable cracks in hydrothermal systems with upward throughflow. Since the finite element method is an approximate numerical method, the method must be validated before it is used to soh,e any new, kind of problem. However, the analytical solution, which can be used to validate the finite element method and other numerical methods, is rather limited in the literature, especially, for the problem considered here. Keeping this in mind, we have derived analytical solutions for the temperature distribution along the vertical axis of a crack in a fluid-saturated porous layer. After the finite element method is validated by comparing the numerical solution with the analytical solution for the same benchmark problem, it is used to investigate the pore-fluid flow and heat transfer in layered hydrothermal systems with vertical permeable cracks. The related analytical and numerical results have demonstrated that vertical cracks are effective and efficient members to transfer heat energy from the bottom section to the top section in hydrothermal systems with upward throughflow.

Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Aspects of the biology of the yabby Trypea australiensis (Dana) (Decapoda : Thalassinidea) and the potential of burrow counts as an indirect measure of population density

The thalassinidean shrimp Trypea australiensis (the yabby) commonly occurs on intertidal sandflats and subtidal regions of sheltered embayments and estuaries along the east coast of Australia and is harvested commercially and recreationally for use as bait by anglers. The potential for counts of burrow openings to provide a reliable indirect estimate of the abundance of yabbies was examined on intertidal sandflats on North Stradbroke Island (Queensland, Australia). The relationship between the number of burrow openings and the abundance of yabbies was generally poor and also varied significantly through time, casting doubt on previous estimates of abundance for this species based on unvalidated hole counts. Spatial and temporal variation in population density, the size at maturity and the reproductive period of the yabby were also assessed. Except for an initial peak in abundance as a result of recruitment, the density of yabbies was constant throughout the study but considerably less than that estimated from a previous study in the same area. Ovigerous females were recorded at 3 mm carapace length (CL) which is smaller than previously recorded for this species and thalassinideans in general. Small ovigerous females were found throughout the study, including the summer months, which is unusual for thalassinideans in the intertidal zone. It was hypothesised that in the intertidal zone, small female yabbies may be able to balance the metabolic demands of reproduction and respiration at higher temperatures than can larger females allowing them to reproduce in the warmer months.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

A randomized trial comparing hormone replacement therapy (HRT) and HRT plus calcitriol in the treatment of postmenopausal osteoporosis with vertebral fractures: Benefit of the combination on total body and hip density

We report a prospective, randomized, multi-center, open-label 2-year trial of 81 postmenopausal women aged 53-79 years with at least one minimal-trauma vertebral fracture (VF) and low (T-score below 2) lumbar bone mineral density (BMD). Group HRT received piperazine estrone sulfate (PES) 0.625 - 1.25 mg/d +/- medroxyprogesterone acetate (MPA) 2.5 - 5 mg/d,- group HRT/D received HRT plus calcitriol 0.25 mug bd. All with a baseline dietary calcium (Ca) of < I g/d received Ca carbonate 0.6 g nocte. Final data were on 66 - 70 patients. On HRT/D, significant (P < 0.001) BNID increases from baseline by DXA were at total body - head, trochanter, Ward's, total hip, inter-trochanter and femoral shaft (% group mean Delta 4.2, 6.1, 9.3. 3.7. 3.3 and 3.3%, respectively). On HRT, at these significant Deltas were restricted to the trochanter and sites. si Wards. Significant advantages of HRT/D over HRT were in BMD of total body (- head), total hip and trochanter (all P = 0.01). The differences in mean Delta at these sites were 1.3, 2.6 and 3.9%. At the following, both groups Improved significantly -lumbar spine (AP and lateral), forearm shaft and ultradistal tibia/fibula. The weightbearing, site - specific benefits of the combination associated with significant suppression of parathyroid hormone-suggest a beneficial effect on cortical bone. Suppression of bone turnover was significantly greater on HRT/D (serum osteocalcin P = 0.024 and urinary hydroxyproline/creatinine ratio P = 0.035). There was no significant difference in the number of patients who developed fresh VFs during the trial (HRT 8/36, 22%; HRT/D 4/34, 12% - intention to treat); likewise in the number who developed incident nonvertebral fractures. This Is the first study comparing the 2 treatments in a fracture population. The results indicate a significant benefit of calcitriol combined with HRT on total body BMD and on BNID at the hip, the major site of osteoporotic fracture.

A 2D numerical model for simulating the physics of fault systems

Simulations provide a powerful means to help gain the understanding of crustal fault system physics required to progress towards the goal of earthquake forecasting. Cellular Automata are efficient enough to probe system dynamics but their simplifications render interpretations questionable. In contrast, sophisticated elasto-dynamic models yield more convincing results but are too computationally demanding to explore phase space. To help bridge this gap, we develop a simple 2D elastodynamic model of parallel fault systems. The model is discretised onto a triangular lattice and faults are specified as split nodes along horizontal rows in the lattice. A simple numerical approach is presented for calculating the forces at medium and split nodes such that general nonlinear frictional constitutive relations can be modeled along faults. Single and multi-fault simulation examples are presented using a nonlinear frictional relation that is slip and slip-rate dependent in order to illustrate the model.

Effects of settler size and density on early post-settlement survival of Ciona intestinalis in the field

Effects of variation in larval quality on post-metamorphic performance in marine invertebrates are increasingly apparent. Recently, it has been shown that variation in offspring size can also strongly affect post-settlement survival, but variation in environmental conditions can mediate this effect. The quality of habitat into which marine invertebrate larvae settle can vary markedly, and 1 influence on quality is the number of conspecifics present. We tested the effects of settler size and settler density on early (1 wk after settlement) post-settlement survival in the field for the solitary ascidian Ciona intestinalis. Larger settlers survived better than smaller settlers, within and among groups of siblings. Increases in the density of settlers decreased survival, but the density-dependent effects were much stronger for smaller settlers. We suggest that larger settlers are better able to cope with intra-specific competition because they have greater energetic reserves or a greater capacity to feed than smaller settlers.