Finite difference time domain (FDTD) method for modeling the effect of switched gradients on the human body in MRI

Numerical modeling of the eddy currents induced in the human body by the pulsed field gradients in MRI presents a difficult computational problem. It requires an efficient and accurate computational method for high spatial resolution analyses with a relatively low input frequency. In this article, a new technique is described which allows the finite difference time domain (FDTD) method to be efficiently applied over a very large frequency range, including low frequencies. This is not the case in conventional FDTD-based methods. A method of implementing streamline gradients in FDTD is presented, as well as comparative analyses which show that the correct source injection in the FDTD simulation plays a crucial rule in obtaining accurate solutions. In particular, making use of the derivative of the input source waveform is shown to provide distinct benefits in accuracy over direct source injection. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and the source injection method has been verified against examples with analytical solutions. Results are presented showing the spatial distribution of gradient-induced electric fields and eddy currents in a complete body model.

Evolution of additive and nonadditive genetic variance in development time along a cline in Drosophila serrata

Latitudinal clines provide natural systems that may allow the effect of natural selection on the genetic variance to be determined. Ten clinal populations of Drosophila serrata collected from the eastern coast of Australia were used to examine clinal patterns in the trait mean and genetic variance of the life-history trait egg-to-adult development time. Development time significantly lengthened from tropical areas to temperate areas. The additive genetic variance for development time in each population was not associated with latitude but was associated with the population mean development time. Additive genetic variance tended to be larger in populations with more extreme development times and appeared to be consistent with allele frequency change. In contrast, the nonadditive genetic variance was not associated with the population mean but was associated with latitude. Levels of nonadditive genetic variance were greatest in the region of the cline where the gradient in the change in mean was greatest, consistent with Barton's (1999) conjecture that the generation of linkage disequilibrium may become an important component of the genetic variance in systems with a spatially varying optimum.

Lateral epicondylalgia: a musculoskeletal physiotherapy perspective

Tennis elbow or lateral epicondylalgia (LE) is a challenging musculoskeletal condition to treat. This is largely due to the lack of research-based evidence of the clinical efficacy of the myriad of treatment approaches espoused in the literature. In view of this, successful rehabilitation of LE is based on choosing treatments that address the physical impairments found during clinical examination. The primary physical impairment in LE is a deficit in grip strength predominately due to pain and its consequences on motor function. Hence the mainstay of successful management of this condition is therapeutic exercise, providing it is not pain provocative. Adjunctive procedures such as manipulative therapy and sports taping techniques have recently been shown to provide substantial initial pain relief. Early relief of pain in the rehabilitation program helps accelerate recovery and most importantly motivates the client to persist with the therapeutic exercise program. The manipulative therapy and taping treatments presented in this masterclass warrant consideration in the clinical best practice management of LE, and serve as a model for other similar musculoskeletal conditions. (C) 2003 Elsevier Science Ltd. All rights reserved.

A new method for analysing the equilibrium and time-dependent behaviour of Markovian models

Many large-scale stochastic systems, such as telecommunications networks, can be modelled using a continuous-time Markov chain. However, it is frequently the case that a satisfactory analysis of their time-dependent, or even equilibrium, behaviour is impossible. In this paper, we propose a new method of analyzing Markovian models, whereby the existing transition structure is replaced by a more amenable one. Using rates of transition given by the equilibrium expected rates of the corresponding transitions of the original chain, we are able to approximate its behaviour. We present two formulations of the idea of expected rates. The first provides a method for analysing time-dependent behaviour, while the second provides a highly accurate means of analysing equilibrium behaviour. We shall illustrate our approach with reference to a variety of models, giving particular attention to queueing and loss networks. (C) 2003 Elsevier Ltd. All rights reserved.

Life cycle evolution in the Digenea: a new perspective from phylogeny

We use a new molecular phylogeny, developed from small and large subunit ribosomal RNA genes, to explore evolution of the digenean life cycle. Our approach is to map character states on the phylogeny and then use parsimony to infer how the character evolved. We conclude that, plesiomorphically, digenean miracidia hatched from eggs and penetrated gastropod first intermediate hosts externally. Fork-tailed cercariae were produced in rediae and emerged from the snail to be eaten directly by the teleost definitive host. These plesiomorphic characters are seen in extant Bivesiculidae. We infer that external encystment and the use of second intermediate hosts are derived from this behaviour and that second intermediate hosts have been adopted repeatedly. Tetrapod definitive hosts have also been adopted repeatedly. The new phylogeny proposes a basal dichotomy between 'Diplostomida' (Diplostomoidea, Schistosomatoidea and Brachylaimoidea) and 'Plagiorchiida' (all other digeneans). There is no evidence for coevolution between these clades and groups of gastropods. The most primitive life cycles are seen in basal Plagiorchiida. Basal Diplostomida have three-host life cycles and are associated with tetrapods. The blood flukes (Schistosomatoidea) are inferred to have derived their two-host life cycles by abbreviating three-host cycles. Diplostomida have no adult stages in fishes except by life cycle abbreviation. We present and test a radical hypothesis that the blood-fluke cycle is plesiomorphic within the Diplostomida.

The variety of variables in automated real-time refinement

The refinement calculus is a well-established theory for deriving program code from specifications. Recent research has extended the theory to handle timing requirements, as well as functional ones, and we have developed an interactive programming tool based on these extensions. Through a number of case studies completed using the tool, this paper explains how the tool helps the programmer by supporting the many forms of variables needed in the theory. These include simple state variables as in the untimed calculus, trace variables that model the evolution of properties over time, auxiliary variables that exist only to support formal reasoning, subroutine parameters, and variables shared between parallel processes.

A high definition, finite difference time domain method

A high definition, finite difference time domain (HD-FDTD) method is presented in this paper. This new method allows the FDTD method to be efficiently applied over a very large frequency range including low frequencies, which are problematic for conventional FDTD methods. In the method, no alterations to the properties of either the source or the transmission media are required. The method is essentially frequency independent and has been verified against analytical solutions within the frequency range 50 Hz-1 GHz. As an example of the lower frequency range, the method has been applied to the problem of induced eddy currents in the human body resulting from the pulsed magnetic field gradients of an MRI system. The new method only requires approximately 0.3% of the source period to obtain an accurate solution. (C) 2003 Elsevier Science Inc. All rights reserved.

An international perspective on sunken state vessels as underwater cultural heritage

On 2 November 2001, the General Assembly of the United Nations Scientific, Economic and Cultural Organisation (UNESCO) adopted the convention on the Protection of the Underwater Cultural Heritage. Among the many complex issues addressed in the convention is the legal status of sunken state-owned vessels, including warships. Prior to the adoption of this convention, no conventional or customary international law existed with regards to the question of abandonment of state-owned vessels or the application of the principle of sovereign immunity to sunken state vessels. While difficulties between coastal states and maritime and former colonial powers resulted in a regime that does not comprehensively address the issues, the convention does provide some guidance in this regard and may provide a basis for further development.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.