A four-stage index 2 Diagonally Implicit Runge-Kutta method

High index Differential Algebraic Equations (DAEs) force standard numerical methods to lower order. Implicit Runge-Kutta methods such as RADAU5 handle high index problems but their fully implicit structure creates significant overhead costs for large problems. Singly Diagonally Implicit Runge-Kutta (SDIRK) methods offer lower costs for integration. This paper derives a four-stage, index 2 Explicit Singly Diagonally Implicit Runge-Kutta (ESDIRK) method. By introducing an explicit first stage, the method achieves second order stage calculations. After deriving and solving appropriate order conditions., numerical examples are used to test the proposed method using fixed and variable step size implementations. (C) 2001 IMACS. Published by Elsevier Science B.V. All rights reserved.

Cross-reactivity of Sydney funnel-web spider antivenom: neutralization of the in vitro toxicity of other Australian funnel-web (Atrax and Hadronyche) spider venoms

Australian funnel-web spiders are recognized as one of the most venomous spiders to humans world-wide. Funnel-web spider antivenom (FWS AV) reverses clinical effects of envenomation from the bite of Atrax robustus and a small number of related Hadronyche species. This study assessed the in vitro efficacy of FWS AV in neutralization of the effects of funnel-web spider venoms, collected from various locations along the eastern seaboard of Australia, in an isolated chick biventer cervicis nerve-muscle preparation. Venoms were separated by SDS-PAGE electrophoresis to compare protein composition and transblotted for Western blotting and incubation with FWS AV. SDS-PAGE of venoms revealed similar low and high molecular weight protein bands. Western blotting with FWS AV showed similar antivenom binding with protein bands in all the venoms tested. Male funnel-web spider venoms (7/7) and female venoms (5110) produced muscle contracture and fasciculation when applied to the nerve-muscle preparation. Venom effects were reversed by subsequent application of FWS AV or prevented by pretreatment of the preparation with antivenom. FWS AV appears to reverse the in vitro toxicity of a number of funnel-web spider venoms from the eastern seaboard of Australia. FWS AV should be effective in the treatment of envenomation from most, if not all, species of Australian funnel-web spiders. (C) 2001 Elsevier Science Ltd. All rights reserved.

Application of a simplified method based on regular regime approach to determine the effective moisture diffusivity of mixture of low molecular weight sugars and maltodextrin during desorption

The binary diffusivities of water in low molecular weight sugars; fructose, sucrose and a high molecular weight carbohydrate; maltodextrin (DE 11) and the effective diffusivities of water in mixtures of these sugars (sucrose, glucose, fructose) and maltodextrin (DE 11) were determined using a simplified procedure based on the Regular Regime Approach. The effective diffusivity of these mixtures exhibited both the concentration and molecular weight dependence. Surface stickiness was observed in all samples during desorption, with fructose exhibiting the highest and maltodextrin the lowest. (C) 2002 Elsevier Science Ltd. All rights reserved.

Comparison of an ELISA and an LC/MS/MS method for measuring tacrolimus concentrations and making dosage decisions in transplant recipients

This study compared an enzyme-linked immunosorbent assay (ELISA) to a liquid chromatography-tandem mass spectrometry (LC/MS/MS) technique for measurement of tacrolimus concentrations in adult kidney and liver transplant recipients, and investigated how assay choice influenced pharmacokinetic parameter estimates and drug dosage decisions. Tacrolimus concentrations measured by both ELISA and LC/MS/MS from 29 kidney (n = 98 samples) and 27 liver (n = 97 samples) transplant recipients were used to evaluate the performance of these methods in the clinical setting. Tacrolimus concentrations measured by the two techniques were compared via regression analysis. Population pharmacokinetic models were developed independently using ELISA and LC/MS/MS data from 76 kidney recipients. Derived kinetic parameters were used to formulate typical dosing regimens for concentration targeting. Dosage recommendations for the two assays were compared. The relation between LC/MS/MS and ELISA measurements was best described by the regression equation ELISA = 1.02 . (LC/MS/MS) + 0.14 in kidney recipients, and ELISA = 1.12 . (LC/MS/MS) - 0.87 in liver recipients. ELISA displayed less accuracy than LC/MS/MS at lower tacrolimus concentrations. Population pharmacokinetic models based on ELISA and LC/MS/MS data were similar with residual random errors of 4.1 ng/mL and 3.7 ng/mL, respectively. Assay choice gave rise to dosage prediction differences ranging from 0% to 30%. ELISA measurements of tacrolimus are not automatically interchangeable with LC/MS/MS values. Assay differences were greatest in adult liver recipients, probably reflecting periods of liver dysfunction and impaired biliary secretion of metabolites. While the majority of data collected in this study suggested assay differences in adult kidney recipients were minimal, findings of ELISA dosage underpredictions of up to 25% in the long term must be investigated further.

Linguistic evaluation method for the environmental assessment of infrastructure projects

A soft linguistic evaluation method is proposed for the environmental assessment of physical infrastructure projects based on fuzzy relations. Infrastructure projects are characterized in terms of linguistic expressions of 'performance' with respect to factors or impacts and the 'importance' of those factors/impacts. A simple example is developed to illustrate the method in the context of three road infrastructure projects assessed against five factors/impacts. In addition, a means to include hard or crisp factors is presented and illustrated with respect to a sixth factor.

Beneficial effects of air inclusions on the performance of ethylene vinyl acetate (EVA) mouthguard material

Objective: To investigate the impact characteristics of an ethylene vinyl acetate (EVA) mouthguard material with regulated air inclusions, which included various air cell volumes and wall thickness between air cells. In particular, the aim was to identify the magnitude and direction of forces within the impacts. Method: EVA mouthguard material, A mm thick and with and without air inclusions, was impacted with a constant force impact pendulum with an energy of 4.4 J and a velocity of 3 m/s. Transmitted forces through the EVA material were measured using an accelerometer, which also allowed the determination of force direction and magnitude within the impacts. Results: Statistically significant reductions in the transmitted forces were observed with all the air inclusion materials when compared with EVA without air inclusions. Maximum transmitted force through one air inclusion material was reduced by 32%. Force rebound was eliminated in one material, and reduced second force impulses were observed in all the air inclusion materials. Conclusion: The regulated air inclusions improved the impact characteristics of the EVA mouthguard material, the material most commonly used in mouthguards world wide.

Evaluating Australia's National Medicines Policy using geographical mapping

Background: There has been a proliferation of quality use of medicines activities in Australia since the 1990s. However, knowledge of the nature and extent of these activities was lacking. A mechanism was required to map the activities to enable their coordination. Aims: To develop a geographical mapping facility as an evaluative tool to assist the planning and implementation of Australia's policy on the quality use of medicines. Methods: A web-based database incorporating geographical mapping software was developed. Quality use of medicines projects implemented across the country was identified from project listings funded by the Quality Use of Medicines Evaluation Program, the National Health and Medical Research Council, Mental Health Strategy, Rural Health Support, Education and Training Program, the Healthy Seniors Initiative, the General Practice Evaluation Program and the Drug Utilisation Evaluation Network. In addition, projects were identified through direct mail to persons working in the field. Results: The Quality Use of Medicines Mapping Project (QUMMP) was developed, providing a Web-based database that can be continuously updated. This database showed the distribution of quality use of medicines activities by: (i) geographical region, (ii) project type, (iii) target group, (iv) stakeholder involvement, (v) funding body and (vi) evaluation method. At September 2001, the database included 901 projects. Sixty-two per cent of projects had been conducted in Australian capital cities, where approximately 63% of the population reside. Distribution of projects varied between States. In Western Australia and Queensland, 36 and 73 projects had been conducted, respectively, representing approximately two projects per 100 000 people. By comparison, in South Australia and Tasmania approximately seven projects per 100 000 people were recorded, with six per 100 000 people in Victoria and three per 100 000 people in New South Wales. Rural and remote areas of the country had more limited project activity. Conclusions: The mapping of projects by geographical location enabled easy identification of high and low activity areas. Analysis of the types of projects undertaken in each region enabled identification of target groups that had not been involved or services that had not yet been developed. This served as a powerful tool for policy planning and implementation and will be used to support the continued implementation of Australia's policy on the quality use of medicines.

In vivo DNA electrotransfer

The use of electrotransfer for DNA delivery to prokaryotic cells, and eukaryotic cells in vitro, has been well known and widely used for many years. However, it is only recently that electric fields have been used to enhance DNA transfer to animal cells in vivo, and this is known as DNA electrotransfer or in vivo DNA electroporation. Some of the advantages of this method of somatic cell gene transfer are that it is a simple method that can be used to transfer almost any DNA construct to animal cells and tissues in vivo; multiple constructs can be co-transfected; it is equally applicable to dividing and nondividing cells; the DNA of interest does not need to be subeloned into a specific viral transfer vector and there is no need for the production of high titre viral stocks; and, as no viral genes are expressed there is less chance of an adverse immunologic reaction to vector sequences. The ease with which efficient in vivo gene transfer can be achieved with in vivo DNA electrotransfer is now allowing genetic analysis to be applied to a number of classic animal model systems where transgenic and embryonic stem cell techniques are not well developed, but for which a wealth of detailed descriptive embryological information is available, or surgical manipulation is much more feasible. As well as exciting applications in developmental biology, in vivo DNA electrotransfer is also being used to transfer genes to skeletal muscle and drive expression of therapeutically active proteins, and to examine exogenous gene and protein function in normal adult cells situated within the complex environment of a tissue and organ system in vivo. Thus, in effect providing the in vivo equivalent of the in vitro transient transfection assay. As the widespread use of in vivo electroporation has really only just begun, it is likely that the future will hold many more applications for this technology in basic research, biotechnology and clinical research areas.

Assessing clinical reasoning: a method to monitor its development in a PBL curriculum

The aim of this study was to develop and trial a method to monitor the evolution of clinical reasoning in a PBL curriculum that is suitable for use in a large medical school. Termed Clinical Reasoning Problems (CRPs), it is based on the notion that clinical reasoning is dependent on the identification and correct interpretation of certain critical clinical features. Each problem consists of a clinical scenario comprising presentation, history and physical examination. Based on this information, subjects are asked to nominate the two most likely diagnoses and to list the clinical features that they considered in formulating their diagnoses, indicating whether these features supported or opposed the nominated diagnoses. Students at different levels of medical training completed a set of 10 CRPs as well as the Diagnostic Thinking Inventory, a self-reporting questionnaire designed to assess reasoning style. Responses were scored against those of a reference group of general practitioners. Results indicate that the CRPs are an easily administered, reliable and valid assessment of clinical reasoning, able to successfully monitor its development throughout medical training. Consequently, they can be employed to assess clinical reasoning skill in individual students and to evaluate the success of undergraduate medical schools in providing effective tuition in clinical reasoning.

Application of the maximum entropy method to the (2 + 1)D four-fermion model

We investigate spectral functions extracted using the maximum entropy method from correlators measured in lattice simulations of the (2+1)-dimensional four-fermion model. This model is particularly interesting because it has both a chirally broken phase with a rich spectrum of mesonic bound states and a symmetric phase where there are only resonances. In the broken phase we study the elementary fermion, pion, sigma, and massive pseudoscalar meson; our results confirm the Goldstone nature of the π and permit an estimate of the meson binding energy. We have, however, seen no signal of σ→ππ decay as the chiral limit is approached. In the symmetric phase we observe a resonance of nonzero width in qualitative agreement with analytic expectations; in addition the ultraviolet behavior of the spectral functions is consistent with the large nonperturbative anomalous dimension for fermion composite operators expected in this model.

Conditional simulation of random fields by successive residuals

This paper presents a new approach to the LU decomposition method for the simulation of stationary and ergodic random fields. The approach overcomes the size limitations of LU and is suitable for any size simulation. The proposed approach can facilitate fast updating of generated realizations with new data, when appropriate, without repeating the full simulation process. Based on a novel column partitioning of the L matrix, expressed in terms of successive conditional covariance matrices, the approach presented here demonstrates that LU simulation is equivalent to the successive solution of kriging residual estimates plus random terms. Consequently, it can be used for the LU decomposition of matrices of any size. The simulation approach is termed conditional simulation by successive residuals as at each step, a small set (group) of random variables is simulated with a LU decomposition of a matrix of updated conditional covariance of residuals. The simulated group is then used to estimate residuals without the need to solve large systems of equations.

Tracing the fate of N-15-enriched feed in an intensive shrimp system

The fate of N-15-nitrogen-enriched formulated feed fed to shrimp was traced through the food web in shallow, outdoor tank systems (1000 1) stocked with shrimp. Triplicate tanks containing shrimp water with and without sediment were used to identify the role of the natural biota in the water column and sediment in processing dietary nitrogen (N). A preliminary experiment demonstrated that N-15-nitrogen-enriched feed products could be detected in the food web. Based on this, a 15-day experiment was conducted. The ammonium (NH4+) pool in the water column became rapidly enriched (within one day) with N-15-nitrogen after shrimp were fed N-15-enriched feed. By day 15, 6% of the added N-15-nitrogen was in this fraction in the 'sediment' tanks compared with 0.4% in the 'no sediment' tanks. The particulate fraction in the water column, principally autotrophic nanoflagellates, accounted for 4-5% of the N-15-nitrogen fed to shrimp after one day. This increased to 16% in the 'no sediment' treatment, and decreased to 2% in the 'sediment' treatment by day 15. It appears that dietary N was more accessible to the phytoplankton community in the absence of sediment. The difference is possibly because a proportion of the dietary N was buried in the sediment in the 'sediment' treatment, making it unavailable to the phytoplankton. Alternatively, the dietary N was retained in the NH4+ pool in the water column since phytoplankton growth, and hence, N utilization was lower in the 'sediment' treatment. The lower growth of phytoplankton in the 'sediment' treatment appeared to be related to higher turbidity, and hence, lower light availability for growth. The percentage N-15-nitrogen detected in the sediment was only 6% despite the high capacity for sedimentation of the large biomass of plankton detritus and shrimp waste. This suggests rapid remineralization of organic waste by the microbial community in the sediment resulting in diffusion of inorganic N sources into the water column. It is likely that most of the dietary N will ultimately be removed from the tank system by water discharges. Our study showed that N-15-nitrogen derived from aquaculture feed can be processed by the microbial community in outdoor aquaculture systems and provides a method for determining the effect of dietary N on ecosystems. However, a significant amount of the dietary N was not retained by the natural biota and is likely to be present in the soluble organic fraction. (C) 2002 Elsevier Science B.V. All rights reserved.

MR image-based measurement of rates of change in volumes of brain structures. Part 1: Method and validation

A detailed analysis procedure is described for evaluating rates of volumetric change in brain structures based on structural magnetic resonance (MR) images. In this procedure, a series of image processing tools have been employed to address the problems encountered in measuring rates of change based on structural MR images. These tools include an algorithm for intensity non-uniforniity correction, a robust algorithm for three-dimensional image registration with sub-voxel precision and an algorithm for brain tissue segmentation. However, a unique feature in the procedure is the use of a fractional volume model that has been developed to provide a quantitative measure for the partial volume effect. With this model, the fractional constituent tissue volumes are evaluated for voxels at the tissue boundary that manifest partial volume effect, thus allowing tissue boundaries be defined at a sub-voxel level and in an automated fashion. Validation studies are presented on key algorithms including segmentation and registration. An overall assessment of the method is provided through the evaluation of the rates of brain atrophy in a group of normal elderly subjects for which the rate of brain atrophy due to normal aging is predictably small. An application of the method is given in Part 11 where the rates of brain atrophy in various brain regions are studied in relation to normal aging and Alzheimer's disease. (C) 2002 Elsevier Science Inc. All rights reserved.

A design method for superconducting MRI magnets with ferromagnetic material

The emphasis of this work is on the optimal design of MRI magnets with both superconducting coils and ferromagnetic rings. The work is directed to the automated design of MRI magnet systems containing superconducting wire and both `cold' and `warm' iron. Details of the optimization procedure are given and the results show combined superconducting and iron material MRI magnets with excellent field characteristics. Strong, homogeneous central magnetic fields are produced with little stray or external field leakage. The field calculations are performed using a semi-analytical method for both current coil and iron material sources. Design examples for symmetric, open and asymmetric clinical MRI magnets containing both superconducting coils and ferromagnetic material are presented.

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.