The design of planar gradient coils. Part I: A winding path correction method

In this work a new approach for designing planar gradient coils is outlined for the use in an existing MRI apparatus. A technique that allows for gradient field corrections inside the diameter-sensitive volume is deliberated. These corrections are brought about by making changes to the wire paths that constitute the coil windings, and hence, is called the path correction method. The existing well-known target held method is used to gauge the performance of a typical gradient coil. The gradient coil design methodology is demonstrated for planar openable gradient coils that can be inserted into an existing MRI apparatus. The path corrected gradient coil is compared to the coil obtained using the target field method. It is shown that using a wire path correction with optimized variables, winding patterns that can deliver high magnetic gradient field strengths and large imaging regions can be obtained.

Tidal influence on behaviour of a coastal aquifer adjacent to a low-relief estuary

The tidal influence on groundwater hydrodynamics, salt-water intrusion and submarine groundwater discharge from coastal/estuarine aquifers is poorly quantified for systems with a mildly sloping beach, in contrast to the case where a vertical beach face is assumed. We investigated the effect of beach slope for a coastal aquifer adjacent to a low-relief estuary, where industrial waste was emplaced over the aquifer. The waste was suspected to discharge leachate towards the estuary. Field observations at various locations showed that tidally induced groundwater head fluctuations were skewed temporally. Frequency analysis suggested that the fluctuation amplitudes decreased exponentially and the phase-tags increased Linearly for the primary tidal signals as they propagated inland. Salinisation zones were observed in the bottom part of the estuary and near the beach surface. Flow and transport processes in a cross-section perpendicular to the estuary were simulated using SEAWAT-2000, which is capable of depicting density-dependent flow and multi-species transport. The simulations showed that the modelled water table fluctuations were in good agreement with the monitored data. Further simulations were conducted to gain insight into the effects of beach slope. In particular the limiting case of a vertical beach face was considered. The simulations showed that density difference and tidal forcing drive a more complex hydrodynamic pattern for the mildly sloping beach than the vertical beach, as well as a profound asymmetry in tidally induced water table fluctuations and enhanced salt-water intrusion. The simulation results also indicated that contaminant transport from the aquifer to the estuary was affected by the tide, where for the mildly sloping beach, the tide tended to intensify the vertical mass exchange in the vicinity of the shorelines, (c) 2005 Elsevier B.V. All rights reserved.

Numerical modelling of tide-induced beach water table fluctuations

Field studies have shown that the elevation of the beach groundwater table varies with the tide and such variations affect significantly beach erosion or accretion. In this paper, we present a BEM (Boundary Element Method) model for simulating the tidal fluctuation of the beach groundwater table. The model solves the two-dimensional flow equation subject to free and moving boundary conditions, including the seepage dynamics at the beach face. The simulated seepage faces were found to agree with the predictions of a simple model (Turner, 1993). The advantage of the present model is, however, that it can be used with little modification to simulate more complicated cases, e.g., surface recharge from rainfall and drainage in the aquifer may be included (the latter is related to beach dewatering technique). The model also simulated well the field data of Nielsen (1990). In particular, the model replicated three distinct features of local water table fluctuations: steep rising phase versus flat falling phase, amplitude attenuation and phase lagging.

Extraction of electrical characteristics from pixels of multifrequency EIT images

Computer modelling has shown that electrical characteristics of individual pixels may be extracted from within multiple-frequency electrical impedance tomography (MFEIT) images formed using a reference data set obtained from a purely resistive, homogeneous medium. In some applications it is desirable to extract the electrical characteristics of individual pixels from images where a purely resistive, homogeneous reference data set is not available. One such application of the technique of MFEIT is to allow the acquisition of in vivo images using reference data sets obtained from a non-homogeneous medium with a reactive component. However, the reactive component of the reference data set introduces difficulties with the extraction of the true electrical characteristics from the image pixels. This study was a preliminary investigation of a technique to extract electrical parameters from multifrequency images when the reference data set has a reactive component. Unlike the situation in which a homogenous, resistive data set is available, it is not possible to obtain the impedance and phase information directly from the image pixel values of the MFEIT images data set, as the phase of the reactive reference is not known. The method reported here to extract the electrical characteristics (the Cole-Cole plot) initially assumes that this phase angle is zero. With this assumption, an impedance spectrum can be directly extracted from the image set. To obtain the true Cole-Cole plot a correction must be applied to account for the inherent rotation of the extracted impedance spectrum about the origin, which is a result of the assumption. This work shows that the angle of rotation associated with the reactive component of the reference data set may be determined using a priori knowledge of the distribution of frequencies of the Cole-Cole plot. Using this angle of rotation, the true Cole-Cole plot can be obtained from the impedance spectrum extracted from the MFEIT image data set. The method was investigated using simulated data, both with and without noise, and also for image data obtained in vitro. The in vitro studies involved 32 logarithmically spaced frequencies from 4 kHz up to 1 MHz and demonstrated that differences between the true characteristics and those of the impedance spectrum were reduced significantly after application of the correction technique. The differences between the extracted parameters and the true values prior to correction were in the range from 16% to 70%. Following application of the correction technique the differences were reduced to less than 5%. The parameters obtained from the Cole-Cole plot may be useful as a characterization of the nature and health of the imaged tissues.

Structure optimization combining soft-core interaction functions, the diffusion equation method, and molecular dynamics

Smoothing the potential energy surface for structure optimization is a general and commonly applied strategy. We propose a combination of soft-core potential energy functions and a variation of the diffusion equation method to smooth potential energy surfaces, which is applicable to complex systems such as protein structures; The performance of the method was demonstrated by comparison with simulated annealing using the refinement of the undecapeptide Cyclosporin A as a test case. Simulations were repeated many times using different initial conditions and structures since the methods are heuristic and results are only meaningful in a statistical sense.

A two-dimensional fuzzy random model of soil pore structure

A new conceptual model for soil pore-solid structure is formalized. Soil pore-solid structure is proposed to comprise spatially abutting elements each with a value which is its membership to the fuzzy set ''pore,'' termed porosity. These values have a range between zero (all solid) and unity (all pore). Images are used to represent structures in which the elements are pixels and the value of each is a porosity. Two-dimensional random fields are generated by allocating each pixel a porosity by independently sampling a statistical distribution. These random fields are reorganized into other pore-solid structural types by selecting parent points which have a specified local region of influence. Pixels of larger or smaller porosity are aggregated about the parent points and within the region of interest by controlled swapping of pixels in the image. This creates local regions of homogeneity within the random field. This is similar to the process known as simulated annealing. The resulting structures are characterized using one-and two-dimensional variograms and functions describing their connectivity. A variety of examples of structures created by the model is presented and compared. Extension to three dimensions presents no theoretical difficulties and is currently under development.

Rapid computation of static fields produced by thick circular solenoids

A straightforward method is proposed for computing the magnetic field produced by a circular coil that contains a large number of turns wound onto a solenoid of rectangular cross section. The coil is thus approximated by a circular ring containing a continuous constant current density, which is very close to the real situation when sire of rectangular cross section is used. All that is required is to evaluate two functions, which are defined as integrals of periodic quantities; this is done accurately and efficiently using trapezoidal-rule quadrature. The solution can be obtained so rapidly that this procedure is ideally suited for use in stochastic optimization, An example is given, in which this approach is combined with a simulated annealing routine to optimize shielded profile coils for NMR.

Robust parameter settings of evolutionary algorithms for the optimisation of agricultural systems models

Numerical optimisation methods are being more commonly applied to agricultural systems models, to identify the most profitable management strategies. The available optimisation algorithms are reviewed and compared, with literature and our studies identifying evolutionary algorithms (including genetic algorithms) as superior in this regard to simulated annealing, tabu search, hill-climbing, and direct-search methods. Results of a complex beef property optimisation, using a real-value genetic algorithm, are presented. The relative contributions of the range of operational options and parameters of this method are discussed, and general recommendations listed to assist practitioners applying evolutionary algorithms to the solution of agricultural systems. (C) 2001 Elsevier Science Ltd. All rights reserved.

The impact of an expert's gender on jurors' decisions

This study uses a simulated civil trial to examine the effect of a male expert's testimony in a male-dominated industry as compared to a female expert's testimony in a traditionally female-dominated industry. ... As noted by Cooper et al., research on persuasion has reliably demonstrated that, under conditions of message complexity, people rely on heuristic cues rather than the content of the message when judging its validity. ... Similarly, Swenson, Nash, and Roos determined that a female expert witness in a child custody dispute was perceived as possessing greater expertise than a male expert, although this difference was only marginally significant. Findings from an unpublished dissertation, which investigated the influence of expert gender in a case involving child sexual abuse, also found some support, in terms of whether or not jurors reached a verdict in a specified period of time or remained hung, for the hypothesis that a female expert would be more influential than her male counterpart. ... Within each of these trial domains (construction, women's clothing), the second experimental variable was manipulated by varying the gender of the plaintiff's expert witness, with half of the participants receiving testimony from a female expert (Dr. Elizabeth Pinder) and half of the participants receiving testimony from a male expert (Dr. Michael Pinder).

Three-dimensional structure of RTD-1, a cyclic antimicrobial defensin from rhesus macaque leukocytes

Most mammalian defensins are cationic peptides of 29-42 amino acids long, stabilized by three disulfide bonds. However, recently Tang et al. (1999, Science 286, 498-502) reported the isolation of a new defensin type found in the leukocytes of rhesus macaques. In contrast to all the other defensins found so far, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue [Tang et al. (1999) Science 286, 498-502]. To elucidate the three-dimensional structure of RTD-1 and its open chain analogue, both peptides were synthesized using solid-phase peptide synthesis and tert-butyloxycarbonyl chemistry. The structures of both peptides in aqueous solution were determined from two-dimensional H-1 NMR data recorded at 500 and 750 MHz. Structural constraints consisting of interproton distances and dihedral angles were used as input for simulated-annealing calculations and water refinement with the program CNS. RTD-1 and its open chain analogue oRTD-1 adopt very similar structures in water. Both comprise an extended beta -hairpin structure with turns at one or both ends. The turns are well defined within themselves and seem to be flexible with respect to the extended regions of the molecules. Although the two strands of the beta -sheet are connected by three disulfide bonds, this region displays a degree of flexibility. The structural similarity of RTD-1 and its open chain analogue oRTD-1, as well as their comparable degree of flexibility, support the theory that the additional charges at the termini of the open chain analogue rather than overall differences in structure or flexibility are the cause for oRTD-1's lower antimicrobial activity. In contrast to numerous other antimicrobial peptides, RTD-1 does not display any amphiphilic character, even though surface models of RTD-1 exhibit a certain clustering of positive charges. Some amide protons of RTD-1 that should be solvent-exposed in monomeric beta -sheet structures show low-temperature coefficients, suggesting the possible presence of weak intermolecular hydrogen bonds.

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.

Denitrification, leaching and immobilisation of N-15-labelled nitrate in winter under windrowed harvesting residues in hoop pine plantations of 1-3 years old in subtropical Australia

A field study was carried out to investigate the impacts of windrowed harvesting residues on denitrification, immobilisation and leaching of N-15-labelled nitrate applied at 20 kg N ha(-1) to microplots in second-rotation hoop pine (Araucaria cunninghamii) plantations of 1-3 years old in southeast Queensland, Australia. The PVC microplots were 235 mm in diameter and 150 mm. long, and driven into the 100 mm soil. There were three replications of such microplots for each of the six treatments which were areas just under and between 1-, 2- and 3-year-old windrows of harvesting residues. Based on gaseous N losses estimated by the difference between the recoveries of bromide (Br) applied at 100 kg Br ha(-1) and N-15-labelled nitrate, denitrification was highest (23% based on N-15 loss) in the areas just under the 1-year-old windrows 25 days after a simulated 75 mm rainfall and following several natural rainfall events. There was no significant difference in N-15 losses (14-17%) among the other treatments. The N-15 immobilisation rate was highest for microplots in the areas between the 1-year-old windrows and generally higher for microplots in the areas just under the windrows (30-39%) than that (26-30%) between the windrows. Direct measurement of N-15 gas emissions (N-15(2) + (N2O)-N-15) confirmed that the highest denitrification rate occurred in the microplots under the 1-year-old windrows although the gaseous N-15 loss calculated by gas emission was only about one-quarter that estimated by the N-15 mass balance method. A significant, positive linear relationship (P < 0.05) existed between the gaseous N-15 losses measured by the two methods used. The research indicates that considerable mineral N could be lost via denitrification during the critical inter-rotation period and early phase of the second rotation. However, the impacts of windrowed harvesting residues on N losses via denitrification might only last for a period of about 2 years. Published by Elsevier Science B.V.

Determination of breeding strategies for beef cattle bull breeders selling bulls into two competing markets on a non-linear price grid

Computer simulation was used to suggest potential selection strategies for beef cattle breeders with different mixes of clients between two potential markets. The traditional market paid on the basis of carcass weight (CWT), while a new market considered marbling grade in addition to CWT as a basis for payment. Both markets instituted discounts for CWT in excess of 340 kg and light carcasses below 300 kg. Herds were simulated for each price category on the carcass weight grid for the new market. This enabled the establishment of phenotypic relationships among the traits examined [CWT, percent intramuscular fat (IMF), carcass value in the traditional market, carcass value in the new market, and the expected proportion of progeny in elite price cells in the new market pricing grid]. The appropriateness of breeding goals was assessed on the basis of client satisfaction. Satisfaction was determined by the equitable distribution of available stock between markets combined with the assessment of the utility of the animal within the market to which it was assigned. The best goal for breeders with predominantly traditional clients was a CWT in excess of 330 kg, while that for breeders with predominantly new market clients was a CWT of between 310 and 329 kg and with a marbling grade of AAA in the Ontario carcass pricing system. For breeders who wished to satisfy both new and traditional clients, the optimal CWT was 310-329 kg and the optimal marbling grade was AA-AAA. This combination resulted in satisfaction levels of greater than 75% among clients, regardless of the distribution of the clients between the traditional and new marketplaces.

The stochastic design of force-minimized compact magnets for high-field magnetic resonance imaging applications

New designs for force-minimized compact high-field clinical MRI magnets are described. The design method is a modified simulated annealing (SA) procedure which includes Maxwell forces in the error function to be minimized. This permits an automated force reduction in the magnet designs while controlling the overall dimensions of the system. As SA optimization requires many iterations to achieve a final design, it is important that each iteration in the procedure is rapid. We have therefore developed a rapid force calculation algorithm. Novel designs for short 3- and 4-T clinical MRI systems are presented in which force reduction has been invoked. The final designs provide large homogeneous regions and reduced stray fields in remarkable short magnets. A shielded 4-T design that is approximately 30% shorter than current designs is presented. This novel magnet generates a full 50-cm diameter homogeneous region.

The plausibility of long-wavelength stress correlation or stress magnitude as a mechanism for precursory seismicity: Results from two simple elastic models

Observations of accelerating seismic activity prior to large earthquakes in natural fault systems have raised hopes for intermediate-term eartquake forecasting. If this phenomena does exist, then what causes it to occur? Recent theoretical work suggests that the accelerating seismic release sequence is a symptom of increasing long-wavelength stress correlation in the fault region. A more traditional explanation, based on Reid's elastic rebound theory, argues that an accelerating sequence of seismic energy release could be a consequence of increasing stress in a fault system whose stress moment release is dominated by large events. Both of these theories are examined using two discrete models of seismicity: a Burridge-Knopoff block-slider model and an elastic continuum based model. Both models display an accelerating release of seismic energy prior to large simulated earthquakes. In both models there is a correlation between the rate of seismic energy release with the total root-mean-squared stress and the level of long-wavelength stress correlation. Furthermore, both models exhibit a systematic increase in the number of large events at high stress and high long-wavelength stress correlation levels. These results suggest that either explanation is plausible for the accelerating moment release in the models examined. A statistical model based on the Burridge-Knopoff block-slider is constructed which indicates that stress alone is sufficient to produce accelerating release of seismic energy with time prior to a large earthquake.