Biomedical subjects

Stuart Crozier

Publications and source records attributed to Stuart Crozier.

At least 19 recordsLinked to original sources

PubMed2006

The role of quadratus lumborum asymmetry in the occurrence of lesions in the lumbar vertebrae of cricket fast bowlers.

In cricket fast bowlers an increased incidence of stress fractures or lesions in the L4 pars interarticularis is observed, which shows a strong statistical correlation with the presence of hypertrophy in the contralateral Quadratus Lumborum (QL) muscle. This study aims to find a physical explanation for this correlation. A mathematical model was used to estimate the forces and moments on the L3 and L4 vertebrae in six postures attained during fast bowling. These forces and moments were used in finite element models to estimate the stresses in the pars interarticularis. Two scenarios were examined per posture: symmetric QL muscles, and right QL muscle volume 30% enlarged. Influence of muscle activation was also investigated. QL asymmetry only correlates with significant stress increases when stress levels are relatively low. When stress levels are high, due to extreme posture or muscle activation, asymmetry only causes small stress changes, suggesting that asymmetry is not the cause of stress fractures in the pars. There are even indications that asymmetry might help to reduce stresses, but more detailed knowledge of the size and activation of the lumbar muscles is needed to confirm this.

Athletic Injuries↗

PubMed2006

Simulation of Brugada syndrome using cellular and three-dimensional whole-heart modeling approaches.

Brugada syndrome (BS) is a genetic disease identified by an abnormal electrocardiogram (ECG) (mainly abnormal ECGs associated with right bundle branch block and ST-elevation in right precordial leads). BS can lead to increased risk of sudden cardiac death. Experimental studies on human ventricular myocardium with BS have been limited due to difficulties in obtaining data. Thus, the use of computer simulation is an important alternative. Most previous BS simulations were based on animal heart cell models. However, due to species differences, the use of human heart cell models, especially a model with three-dimensional whole-heart anatomical structure, is needed. In this study, we developed a model of the human ventricular action potential (AP) based on refining the ten Tusscher et al (2004 Am. J. Physiol. Heart Circ. Physiol. 286 H1573-89) model to incorporate newly available experimental data of some major ionic currents of human ventricular myocytes. These modified channels include the L-type calcium current (I(CaL)), fast sodium current (I(Na)), transient outward potassium current (I(to)), rapidly and slowly delayed rectifier potassium currents (I(Kr) and I(Ks)) and inward rectifier potassium current (I(Ki)). Transmural heterogeneity of APs for epicardial, endocardial and mid-myocardial (M) cells was simulated by varying the maximum conductance of I(Ks) and I(to). The modified AP models were then used to simulate the effects of BS on cellular AP and body surface potentials using a three-dimensional dynamic heart-torso model. Our main findings are as follows. (1) BS has little effect on the AP of endocardial or mid-myocardial cells, but has a large impact on the AP of epicardial cells. (2) A likely region of BS with abnormal cell AP is near the right ventricular outflow track, and the resulting ST-segment elevation is located in the median precordium area. These simulation results are consistent with experimental findings reported in the literature. The model can reproduce a variety of electrophysiological behaviors and provides a good basis for understanding the genesis of abnormal ECG under the condition of BS disease.

Action Potentials↗

PubMed2006

Effect of cardiac motion on body surface electrocardiographic potentials: an MRI-based simulation study.

This paper describes an electrical model of cardiac ventricles incorporating real geometry and motion. The heart anatomy and its motion through the cardiac cycle are obtained from segmentations of multiple-slice MRI time sequences; the special conduction system is constructed using an automated mapping procedure from an existing static heart model. The heart model is mounted in an anatomically realistic voxel model of the human body. The cardiac electrical source and surface potentials are determined numerically using both a finite-difference scheme and a boundary-element method with the incorporation of the motion of the heart. The electrocardiograms (ECG) and body surface potential maps are calculated and compared to the static simulation in the resting heart. The simulations demonstrate that introducing motion into the cardiac model modifies the ECG signals, with the most obvious change occurring during the T-wave at peak contraction of the ventricles. Body surface potential maps differ in some local positions during the T-wave, which may be of importance to a number of cardiac models, including those incorporating inverse methods.

Algorithms↗

PubMed2006

Transient temperature rise in a mouse due to low-frequency regional hyperthermia.

A refined nonlinear heat transfer model of a mouse has been developed to simulate the transient temperature rise in a neoplastic tumour and neighbouring tissue during regional hyperthermia using a 150 kHz inductive coil. In this study, we incorporate various bio-energetic enhancements to the heat transfer equation and numerical validations based on experimental findings for the mouse, in terms of nonlinear metabolic heat production, homeothermy, blood perfusion parameters, thermoregulation, psychological and physiological effects. The discretized bio-heat transfer equation has been validated with the commercial software FEMLAB on a canonical multi-sphere object before applying the scheme to the inhomogeneous mouse voxel phantom. The time-dependent numerical results of regional hyperthermia of mouse thigh have been compared with the available experimental temperature results with only a few small disparities. During the first 20 min of local unfocused heating, the temperature in the tumour and the surrounding tissue increased by around 7.5 degrees C. The objective of this preliminary study was to develop a validated electrothermal numerical scheme for inductive hyperthermia of a small mammal with the intention of expanding the model into a complete numerical solution involving ferromagnetic nanoparticles for targeted heating of tumours at low frequencies. In addition, the numerical scheme herein could assist in optimizing and tailoring of focused electromagnetic fields for hyperthermia.

Animals↗

PubMed2005

A theoretical comparison of two optimization methods for radiofrequency drive schemes in high frequency MRI resonators.

In this paper, numerical simulations are used in an attempt to find optimal source profiles for high frequency radiofrequency (RF) volume coils. Biologically loaded, shielded/unshielded circular and elliptical birdcage coils operating at 170 MHz, 300 MHz and 470 MHz are modelled using the FDTD method for both 2D and 3D cases. Taking advantage of the fact that some aspects of the electromagnetic system are linear, two approaches have been proposed for the determination of the drives for individual elements in the RF resonator. The first method is an iterative optimization technique with a kernel for the evaluation of RF fields inside an imaging plane of a human head model using pre-characterized sensitivity profiles of the individual rungs of a resonator; the second method is a regularization-based technique. In the second approach, a sensitivity matrix is explicitly constructed and a regularization procedure is employed to solve the ill-posed problem. Test simulations show that both methods can improve the B(1)-field homogeneity in both focused and non-focused scenarios. While the regularization-based method is more efficient, the first optimization method is more flexible as it can take into account other issues such as controlling SAR or reshaping the resonator structures. It is hoped that these schemes and their extensions will be useful for the determination of multi-element RF drives in a variety of applications.

Computer Simulation↗

PubMed2005

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model.

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy of ECG simulation when it is used in heart-torso model-based body surface potential simulation studies.

Computer Simulation↗

PubMed2005

Focused, eight-element transceive phased array coil for parallel magnetic resonance imaging of the chest--theoretical considerations.

A new transceive system for chest imaging for MRI applications is presented. A focused, eight-element transceive torso phased array coil is designed to investigate transmitting a focused radiofrequency field deep within the torso and to enhance signal homogeneity in the heart region. The system is used in conjunction with the SENSE reconstruction technique to enable focused parallel imaging. A hybrid finite-difference-time-domain/method-of-moments method is used to accurately predict the radiofrequency behavior inside the human torso. The simulation results reported herein demonstrate the feasibility of the design concept, which shows that radiofrequency field focusing with SENSE reconstruction is theoretically achievable.

Algorithms↗

PubMed2005

An inverse methodology for high-frequency RF coil design for MRI with de-emphasized B1 fields.

An inverse methodology for the design of biologically loaded radio-frequency (RF) coils for magnetic resonance imaging applications is described. Free space time-harmonic electromagnetic Green's functions and de-emphasized B1 target fields are used to calculate the current density on the coil cylinder. In theory, with the B1 field de-emphasized in the middle of the RF transverse plane, the calculated current distribution can generate an internal magnetic field that can reduce the central overemphasis effect caused by field/tissue interactions at high frequencies. The current distribution of a head coil operating at 4 T (170 MHz) is calculated using an inverse methodology with de-emphasized B1 target fields. An in-house finite-difference time-domain routine is employed to evaluate B1 field and signal intensity inside a homogenous cylindrical phantom and then a complete human head model. A comparison with a conventional RF birdcage coil is carried out and demonstrates that this method can help in decreasing the normal bright region caused by field/tissue interactions in head images at 170 MHz and higher field strengths.

Brain↗

PubMed2005

The use of unwrapped phase in MR image segmentation: a preliminary study.

This paper considers the problem of tissue classification in 3D MRI. More specifically, a new set of texture features, based on phase information, is used to perform the segmentation of the bones of the knee. The phase information provides a very good discrimination between the bone and the surrounding tissues, but is usually not used due to phase unwrapping problems. We present a method to extract textural information from the phase that does not require phase unwrapping. The textural information extracted from the magnitude and the phase can be combined to perform tissue classification, and used to initialise an active shape model, leading to a more precise segmentation.

Algorithms↗

PubMed2004

Numerical modelling of thermal effects in rats due to high-field magnetic resonance imaging (0.5-1 GHZ).

A finite-difference time-domain (FDTD) thermal model has been developed to compute the temperature elevation in the Sprague Dawley rat due to electromagnetic energy deposition in high-field magnetic resonance imaging (MRI). The field strengths examined ranged from 11.75-23.5 T (corresponding to 1H resonances of 0.5-1 GHz) and an N-stub birdcage resonator was used to both transmit radio-frequency energy and receive the MRI signals. With an in-plane resolution of 1.95 mm, the inhomogeneous rat phantom forms a segmented model of 12 different tissue types, each having its electrical and thermal parameters assigned. The steady-state temperature distribution was calculated using a Pennes 'bioheat' approach. The numerical algorithm used to calculate the induced temperature distribution has been successfully validated against analytical solutions in the form of simplified spherical models with electrical and thermal properties of rat muscle. As well as assisting with the design of MRI experiments and apparatus, the numerical procedures developed in this study could help in future research and design of tumour-treating hyperthermia applicators to be used on rats in vivo.

Animals↗

PubMed2004

A diffusion and T2 relaxation MRI study of the ovine lumbar intervertebral disc under compression in vitro.

The ovine lumbar intervertebral disc is a useful model for the human lumbar disc. We present preliminary estimates of diffusion coefficients and T2 relaxation times in a pilot MRI study of the ovine lumbar intervertebral disc during uniaxial compression in vitro, and identify factors that hamper the ability to accurately monitor the temporal evolution of the effective diffusion tensor at high spatial resolution.

Animals↗

PubMed2004

Electromagnetic fields inside a lossy, multilayered spherical head phantom excited by MRI coils: models and methods.

The precise evaluation of electromagnetic field (EMF) distributions inside biological samples is becoming an increasingly important design requirement for high field MRI systems. In evaluating the induced fields caused by magnetic field gradients and RF transmitter coils, a multilayered dielectric spherical head model is proposed to provide a better understanding of electromagnetic interactions when compared to a traditional homogeneous head phantom. This paper presents Debye potential (DP) and Dyadic Green's function (DGF)-based solutions of the EMFs inside a head-sized, stratified sphere with similar radial conductivity and permittivity profiles as a human head. The DP approach is formulated for the symmetric case in which the source is a circular loop carrying a harmonic-formed current over a wide frequency range. The DGF method is developed for generic cases in which the source may be any kind of RF coil whose current distribution can be evaluated using the method of moments. The calculated EMFs can then be used to deduce MRI imaging parameters. The proposed methods, while not representing the full complexity of a head model, offer advantages in rapid prototyping as the computation times are much lower than a full finite difference time domain calculation using a complex head model. Test examples demonstrate the capability of the proposed models/methods. It is anticipated that this model will be of particular value for high field MRI applications, especially the rapid evaluation of RF resonator (surface and volume coils) and high performance gradient set designs.

Algorithms↗

PubMed2004

A distributed equivalent magnetic current based FDTD method for the calculation of E-fields induced by gradient coils.

This paper evaluates a new, low-frequency finite-difference time-domain method applied to the problem of induced E-fields/eddy currents in the human body resulting from the pulsed magnetic field gradients in MRI. In this algorithm, a distributed equivalent magnetic current is proposed as the electromagnetic source and is obtained by quasistatic calculation of the empty coil's vector potential or measurements therein. This technique circumvents the discretization of complicated gradient coil geometries into a mesh of Yee cells, and thereby enables any type of gradient coil modelling or other complex low frequency sources. The proposed method has been verified against an example with an analytical solution. Results are presented showing the spatial distribution of gradient-induced electric fields in a multi-layered spherical phantom model and a complete body model.

Algorithms↗

PubMed2003

Influence of magnetically-induced E-fields on cardiac electric activity during MRI: A modeling study.

In modern magnetic resonance imaging (MRI), patients are exposed to strong, time-varying gradient magnetic fields that may be able to induce electric fields (E-fields)/currents in tissues approaching the level of physiological significance. In this work we present theoretical investigations into induced E-fields in the thorax, and evaluate their potential influence on cardiac electric activity under the assumption that the sites of maximum E-field correspond to the myocardial stimulation threshold (an abnormal circumstance). Whole-body cylindrical and planar gradient coils were included in the model. The calculations of the induced fields are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic, whole-body model. The potential for cardiac stimulation was evaluated using an electrical model of the heart. Twelve-lead electrocardiogram (ECG) signals were simulated and inspected for arrhythmias caused by the applied fields for both healthy and diseased hearts. The simulations show that the shape of the thorax and the conductive paths significantly influence induced E-fields. In healthy patients, these fields are not sufficient to elicit serious arrhythmias with the use of contemporary gradient sets. However, raising the strength and number of repeated switching episodes of gradients, as is certainly possible in local chest gradient sets, could expose patients to increased risk. For patients with cardiac disease, the risk factors are elevated. By the use of this model, the sensitivity of cardiac pathologies, such as abnormal conductive pathways, to the induced fields generated by an MRI sequence can be investigated.

Action Potentials↗

PubMed2003

Calculation of electric fields induced by body and head motion in high-field MRI.

In modern magnetic resonance imaging (MRI), patients are exposed to strong, nonuniform static magnetic fields outside the central imaging region, in which the movement of the body may be able to induce electric currents in tissues which could be possibly harmful. This paper presents theoretical investigations into the spatial distribution of induced electric fields and currents in the patient when moving into the MRI scanner and also for head motion at various positions in the magnet. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic, full-body, male model. 3D field profiles from an actively shielded 4T magnet system are used and the body model projected through the field profile with a range of velocities. The simulation shows that it possible to induce electric fields/currents near the level of physiological significance under some circumstances and provides insight into the spatial characteristics of the induced fields. The results are extrapolated to very high field strengths and tabulated data shows the expected induced currents and fields with both movement velocity and field strength.

Computer Simulation↗

PubMed2003

On the induced electric field gradients in the human body for magnetic stimulation by gradient coils in MRI.

Prior theoretical studies indicate that the negative spatial derivative of the electric field induced by magnetic stimulation may be one of the main factors contributing to depolarization of the nerve fiber. This paper studies this parameter for peripheral nerve stimulation (PNS) induced by time-varying gradient fields during MRI scans. The numerical calculations are based on an efficient, quasi-static, finite-difference scheme and an anatomically realistic human, full-body model. Whole-body cylindrical and planar gradient sets in MRI systems and various input signals have been explored. The spatial distributions of the induced electric field and their gradients are calculated and attempts are made to correlate these areas with reported experimental stimulation data. The induced electrical field pattern is similar for both the planar coils and cylindrical coils. This study provides some insight into the spatial characteristics of the induced field gradients for PNS in MRI, which may be used to further evaluate the sites where magnetic stimulation is likely to occur and to optimize gradient coil design.

Adult↗

PubMed2002

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.

Anatomy, Cross-Sectional↗

PubMed2002

A method for the design of MRI radiofrequency coils based on triangular and pulse basis functions.

This paper presents a numerical technique for the design of an RF coil for asymmetric magnetic resonance imaging (MRI) systems. The formulation is based on an inverse approach where the cylindrical surface currents are expressed in terms of a combination of sub-domain basis functions: triangular and pulse functions. With the homogeneous transverse magnetic field specified in a spherical region, a functional method is applied to obtain the unknown current coefficients. The current distribution is then transformed to a conductor pattern by use of a stream function technique. Preliminary MR images acquired using a prototype RF coil are presented and validate the design method.

Electromagnetic Fields↗