Multiple-acquisition parallel imaging combined with a transceive array for the amelioration of high-field RF distortion: a modeling study

A new method for ameliorating high-field image distortion caused by radio frequency/tissue interaction is presented and modeled, The proposed method uses, but is not restricted to, a shielded four-element transceive phased array coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both images together, the image distortion can be reduced several-fold. A hybrid finite-difference time-domain/method-of-moments method is used to theoretically demonstrate the method and also to predict the radio frequency behavior inside the human head. in addition, the proposed method is used in conjunction with the GRAPPA reconstruction technique to enable rapid imaging. Simulation results reported herein for IIT (470 MHz) brain imaging applications demonstrate the feasibility of the concept where multiple acquisitions using parallel imaging elements with GRAPPA reconstruction results in improved image quality. (c) 2006 Wiley Periodicals, Inc.

Many-body quantum dynamics of polarization squeezing in optical fibers

We report new experiments that test quantum dynamical predictions of polarization squeezing for ultrashort photonic pulses in a birefringent fiber, including all relevant dissipative effects. This exponentially complex many-body problem is solved by means of a stochastic phase-space method. The squeezing is calculated and compared to experimental data, resulting in excellent quantitative agreement. From the simulations, we identify the physical limits to quantum noise reduction in optical fibers. The research represents a significant experimental test of first-principles time-domain quantum dynamics in a one-dimensional interacting Bose gas coupled to dissipative reservoirs.

High-field magnetic resonance imaging with reduced field/tissue RF artefacts - A modeling study using hybrid MoM/FEM and FDTD technique

Due to complex field/tissue interactions, high-field magnetic resonance (MR) images suffer significant image distortions that result in compromised diagnostic quality. A new method that attempts to remove these distortions is proposed in this paper and is based on the use of transceiver-phased arrays. The proposed system uses, in the examples presented herein, a shielded four-element transceive-phased array head coil and involves performing two separate scans of the same slice with each scan using different excitations during transmission. By optimizing the amplitudes and phases for each scan, antipodal signal profiles can be obtained, and by combining both the images together, the image distortion can be reduced several fold. A combined hybrid method of moments (MoM)/finite element method (FEM) and finite-difference time-domain (FDTD) technique is proposed and used to elucidate the concept of the new method and to accurately evaluate the electromagnetic field (EMF) in a human head model. In addition, the proposed method is used in conjunction with the generalized auto-calibrating partially parallel acquisitions (GRAPPA) reconstruction technique to enable rapid imaging of the two scans. Simulation results reported herein for 11-T (470-MHz) brain imaging applications show that the new method with GRAPPA reconstruction theoretically results in improved image quality and that the proposed combined hybrid MoM/FEM and FDTD technique is. suitable for high-field magnetic resonance imaging (MRI) numerical analysis.

An object-orient designed FDTD simulator: Applications to high field systems

An object-oriented finite-difference time-domain (FDTD) simulator has been developed for electromagnetic study and design applications in Magnetic Resonance Imaging. It is aimed to be a complete FDTD model of an MRI system including all high and low-frequency field generating units and electrical models of the patient. The design method is described and MRI-based numerical examples are presented to illustrate the function of the numerical solver, particular emphasis is placed on high field studies.

Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

Customer information system data pre-processing with feature selection techniques for non-technical losses prediction in an electricity market

Non-technical losses (NTL) identification and prediction are important tasks for many utilities. Data from customer information system (CIS) can be used for NTL analysis. However, in order to accurately and efficiently perform NTL analysis, the original data from CIS need to be pre-processed before any detailed NTL analysis can be carried out. In this paper, we propose a feature selection based method for CIS data pre-processing in order to extract the most relevant information for further analysis such as clustering and classifications. By removing irrelevant and redundant features, feature selection is an essential step in data mining process in finding optimal subset of features to improve the quality of result by giving faster time processing, higher accuracy and simpler results with fewer features. Detailed feature selection analysis is presented in the paper. Both time-domain and load shape data are compared based on the accuracy, consistency and statistical dependencies between features.

Tailoring magnetic field gradient design to magnet cryostat geometry

Eddy currents induced within a magnetic resonance imaging (MRI) cryostat bore during pulsing of gradient coils can be applied constructively together with the gradient currents that generate them, to obtain good quality gradient uniformities within a specified imaging volume over time. This can be achieved by simultaneously optimizing the spatial distribution and temporal pre-emphasis of the gradient coil current, to account for the spatial and temporal variation of the secondary magnetic fields due to the induced eddy currents. This method allows the tailored design of gradient coil/magnet configurations and consequent engineering trade-offs. To compute the transient eddy currents within a realistic cryostat vessel, a low-frequency finite-difference time-domain (FDTD) method using total-field scattered-field (TFSF) scheme has been performed and validated

FDTD study of an UWB radar technique for breast tumor detection

A finite difference time domain (FDTD) method is applied to investigate capabilities of an ultra-wide band (UWB) radar system to detect a breast tumor. The first part of the investigations concerns FDTD simulations of a phantom formed by a plastic container with liquid and a small reflecting target. The second part focuses on a three-dimensional numerical breast model with a small tumor. FDTD simulations are carried out assuming a planar incident wave. Various time snap shots of the electromagnetic field are recorded to learn about the physical phenomenon of reflection and scattering in different layers of the phantom.

Facilitating soil-and plant-water research through the use of TDR

The development of TDR for measurement of soil water content and electrical conductivity has resulted in a large shift in measurement methods for a breadth of soil and hydrological characterization efforts. TDR has also opened new possibilities for soil and plant research. Five examples show how TDR has enhanced our ability to conduct our soil- and plant-water research. (i) Oxygen is necessary for healthy root growth and plant development but quantitative evaluation of the factors controlling oxygen supply in soil depends on knowledge of the soil water content by TDR. With water content information we have modeled successfully some impact of tillage methods on oxygen supply to roots and their growth response. (ii) For field assessment of soil mechanical properties influencing crop growth, water content capability was added to two portable soil strength measuring devices; (a) A TDT (Time Domain Transmittivity)-equipped soil cone penetrometer was used to evaluate seasonal soil strengthwater content relationships. In conventional tillage systems the relationships are dynamic and achieve the more stable no-tillage relationships only relatively late in each growing season; (b) A small TDR transmission line was added to a modified sheargraph that allowed shear strength and water content to be measured simultaneously on the same sample. In addition, the conventional graphing procedure for data acquisition was converted to datalogging using strain gauges. Data acquisition rate was improved by more than a factor of three with improved data quality. (iii) How do drought tolerant plants maintain leaf water content? Non-destructive measurement of TDR water content using a flat serpentine triple wire transmission line replaces more lengthy procedures of measuring relative water content. Two challenges remain: drought-stressed leaves alter salt content, changing electrical conductivity, and drought induced changes in leaf morphology affect TDR measurements. (iv) Remote radar signals are reflected from within the first 2 cm of soil. Appropriate calibration of radar imaging for soil water content can be achieved by a parallel pair of blades separated by 8 cm, reaching 1.7 cm into soil and forming a 20 cm TDR transmission line. The correlation between apparent relative permittivity from TDR and synthetic aperture radar (SAR) backscatter coefficient was 0.57 from an airborne flyover. These five examples highlight the diversity in the application of TDR in soil and plant research.

Mathematical modelling and optimal multivariable control of chemical processes

This work reports the developnent of a mathenatical model and distributed, multi variable computer-control for a pilot plant double-effect climbing-film evaporator. A distributed-parameter model of the plant has been developed and the time-domain model transformed into the Laplace domain. The model has been further transformed into an integral domain conforming to an algebraic ring of polynomials, to eliminate the transcendental terms which arise in the Laplace domain due to the distributed nature of the plant model. This has made possible the application of linear control theories to a set of linear-partial differential equations. The models obtained have well tracked the experimental results of the plant. A distributed-computer network has been interfaced with the plant to implement digital controllers in a hierarchical structure. A modern rnultivariable Wiener-Hopf controller has been applled to the plant model. The application has revealed a limitation condition that the plant matrix should be positive-definite along the infinite frequency axis. A new multi variable control theory has emerged fram this study, which avoids the above limitation. The controller has the structure of the modern Wiener-Hopf controller, but with a unique feature enabling a designer to specify the closed-loop poles in advance and to shape the sensitivity matrix as required. In this way, the method treats directly the interaction problems found in the chemical processes with good tracking and regulation performances. Though the ability of the analytical design methods to determine once and for all whether a given set of specifications can be met is one of its chief advantages over the conventional trial-and-error design procedures. However, one disadvantage that offsets to some degree the enormous advantages is the relatively complicated algebra that must be employed in working out all but the simplest problem. Mathematical algorithms and computer software have been developed to treat some of the mathematical operations defined over the integral domain, such as matrix fraction description, spectral factorization, the Bezout identity, and the general manipulation of polynomial matrices. Hence, the design problems of Wiener-Hopf type of controllers and other similar algebraic design methods can be easily solved.

Cognition, ocular accommodation, and cardiovascular function in emmetropes and late-onset myopes

PURPOSE. To investigate objectively and noninvasively the role of cognitive demand on autonomic control of systemic cardiovascular and ocular accommodative responses in emmetropes and myopes of late-onset. METHODS. Sixteen subjects (10 men, 6 women) aged between 18 and 34 years (mean ± SD: 22.6 ± 4.4 years), eight emmetropes (EMMs; mean spherical equivalent [MSE] refractive error ± SD: 0.05 ± 0.24 D) and eight with late-onset myopia (LOMs; MSE ± SD: -3.66 ± 2.31 D) participated in the study. Subjects viewed stationary numerical digits monocularly within a Badal optical system (at both 0.0 and -3.0 D) while performing a two-alternative, forced-choice paradigm that matched cognitive loading across subjects. Five individually matched cognitive levels of increasing difficulty were used in random order for each subject. Five 20-second, continuous-objective recordings of the accommodative response measured with an open-view infrared autorefractor were obtained for each cognitive level, whereas simultaneous measurement of heart rate was continuously recorded with a finger-mounted piezoelectric pulse transducer for 5 minutes. Fast Fourier transformation of cardiovascular function allowed the relative power of the autonomic components to be assessed in the frequency domain, whereas heart period gave an indication of the time-domain response. RESULTS. Increasing the cognitive demand led to a significant reduction in the accommodative response in all subjects (0.0 D: by -0.35 ± 0.33 D; -3.0 D: by -0.31 ± 0.40 D, P < 0.001). The greater lag of LOMs compared with EMMs was not significant (P = 0.07) at both distance (0.38 ± 0.35 D) and near (0.14 ± 0.42 D). Mean heart period reduced with increasing levels of workload (P < 0.0005). LOMs exhibited a relative elevation in sympathetic system activity compared to EMMs. Within refractive groups, however, accommodative shifts with increasing cognition correlated with parasympathetic activity (r = 0.99, P < 0.001), more than with sympathetic activity (r = 0.62, P > 0.05). CONCLUSIONS. In an equivalent workload paradigm, increasing cognitive demand caused a reduction in accommodative response that was attributable principally to a concurrent reduction in the relative power of the parasympathetic component of the autonomic nervous system (ANS). The disparity in accommodative response between EMMs and LOMs, however, appears to be augmented by changes in the sympathetic nervous component of the systemic ANS. Copyright © Association for Research in Vision and Ophthalmology.

Autonomic correlates of ocular accommodation and cardiovascular function

PURPOSE: To evaluate the hypothesis that objective measures of open- and closed-loop ocular accommodation are related to systemic cardiovascular function, and ipso facto autonomic nervous system activity. METHODS: Sixty subjects (29 male; 31 female) varying in age from 18 to 33 years (average: 20.3 +/- 2.9 years) with a range of refractive errors [mean spherical equivalent (MSE): -7.12 to +1.82 D] participated in the study. Five 20-s continuous objective recordings of the accommodative response, measured with an open-view IR autorefractor (Shin-Nippon SRW-5000), were obtained for a variety of open- and closed-loop accommodative demands while simultaneous continuous measurement of heart rate was recorded with a finger-mounted piezo-electric pulse transducer for 5 min. Fast Fourier Transformation of cardiovascular function allowed the absolute and relative power of the autonomic components to be assessed in the frequency-domain, whereas heart period gave an indication of the time-domain response. RESULTS: Increasing closed-loop accommodative demand led to a concurrent increase in heart rate of approximately 2 beats/min for a 4.0 D increase in accommodative demand. The increase was attributable to a reduction in the absolute (p < 0.05) and normalised (p < 0.001) input of the systemic parasympathetic nervous system, and was unaffected by refractive group. The interaction with refractive group failed to reach significance. CONCLUSIONS: For sustained accommodation effort, the data demonstrate covariation between the oculomotor and cardiovascular systems which implies that a near visual task can significantly influence cardiovascular behaviour. Accommodative effort alone, however, is not a sufficient stimulus to induce autonomic differences between refractive groups. The data suggest that both the oculomotor and cardiovascular systems are predominantly attributable to changes in the systemic parasympathetic nervous system.

Temporal finite elements in the dynamic analysis of structures

This thesis demonstrates that the use of finite elements need not be confined to space alone, but that they may also be used in the time domain, It is shown that finite element methods may be used successfully to obtain the response of systems to applied forces, including, for example, the accelerations in a tall structure subjected to an earthquake shock. It is further demonstrated that at least one of these methods may be considered to be a practical alternative to more usual methods of solution. A detailed investigation of the accuracy and stability of finite element solutions is included, and methods of applications to both single- and multi-degree of freedom systems are described. Solutions using two different temporal finite elements are compared with those obtained by conventional methods, and a comparison of computation times for the different methods is given. The application of finite element methods to distributed systems is described, using both separate discretizations in space and time, and a combined space-time discretization. The inclusion of both viscous and hysteretic damping is shown to add little to the difficulty of the solution. Temporal finite elements are also seen to be of considerable interest when applied to non-linear systems, both when the system parameters are time-dependent and also when they are functions of displacement. Solutions are given for many different examples, and the computer programs used for the finite element methods are included in an Appendix.

Force sensing for measuring human body movement

The research developed in this thesis explores the sensing and inference of human movement in a dynamic way, as opposed to conventional measurement systems, that are only concerned with discrete evaluations of stimuli in sequential time. Typically, conventional approaches are used to infer the dynamic movement of the body; such as vision and motion tracking devices, with either a human diagnosis or complex image processing algorithm to classify the movement. This research is therefore the first of its kind to attempt and provide a movement classifying algorithm through the use of minimal sensing points, with the application for this novel system, to classify human movement during a golf swing. There are two main categories of force sensing. Firstly, array-type systems consisting of many sensing elements, and are the most commonly researched and commercially available. Secondly, reduced force sensing element systems (RFSES) also known as distributive systems have only been recently exploited in the academic world. The fundamental difference between these systems is that array systems handle the data captured from each sensor as unique outputs and suffer the effects of resolution. The effect of resolution, is the error in the load position measurement between sensing elements, as the output is quantized in terms of position. This can be compared to a reduced sensor element system that maximises that data received through the coupling of data from a distribution of sensing points to describe the output in discrete time. Also this can be extended to a coupling of transients in the time domain to describe an activity or dynamic movement. It is the RFSES that is to be examined and exploited in the commercial sector due to its advantages over array-based approaches such as reduced design, computational complexity and cost.

The application of parallel processing methods to vibrational analysis

The trend in modal extraction algorithms is to use all the available frequency response functions data to obtain a global estimate of the natural frequencies, damping ratio and mode shapes. Improvements in transducer and signal processing technology allow the simultaneous measurement of many hundreds of channels of response data. The quantity of data available and the complexity of the extraction algorithms make considerable demands on the available computer power and require a powerful computer or dedicated workstation to perform satisfactorily. An alternative to waiting for faster sequential processors is to implement the algorithm in parallel, for example on a network of Transputers. Parallel architectures are a cost effective means of increasing computational power, and a larger number of response channels would simply require more processors. This thesis considers how two typical modal extraction algorithms, the Rational Fraction Polynomial method and the Ibrahim Time Domain method, may be implemented on a network of transputers. The Rational Fraction Polynomial Method is a well known and robust frequency domain 'curve fitting' algorithm. The Ibrahim Time Domain method is an efficient algorithm that 'curve fits' in the time domain. This thesis reviews the algorithms, considers the problems involved in a parallel implementation, and shows how they were implemented on a real Transputer network.