Pulse sequence design for MR velocity mapping of complex flow: Notes on the necessity of low echo times

Lowering of the echo time (TE) has been proposed as a way to reduce effects of phase dispersion in MR velocity mapping, because a low TE reduces sensitivity to higher-order motion terms while first-order velocity sensitivity is maintained. Methods of lowering TE involves the use of extreme gradient ramp times and gradient strengths as well as reduction of the duration of transmit/receive windows, the latter method causing decrements in image resolution. When reducing higher-order sensitivity, however, it is not the overall TE that is the critical parameter, but rather the time pattern of the gradients used in the experiment. Hence, changes in TE without subsequent variations in gradient pattern would, according to theory, not affect quantitative measurements of complex flow and vice versa. In this study, we experimentally demonstrate this relation and utilize the experience to create a sequence robust towards complex flow without sacrifices in image resolution. Our experimental observations show that variations in TE alone while maintaining the time course of the velocity-encoding gradient does not significantly affect measurements of through-plane average complex flow in the studied velocity range. A parameter that cannot be measured as accurately if TE is increased is the peak flow. A phase mapping sequence with prolonged TE from 3 ms to 5 ms but with short duration of the velocity-encoding (section-selective) gradient and improved in-plane resolution was demonstrated in vivo.     

Application of linear optimization techniques to MRI phase contrast blood flow measurements

The goal of this study was to use linear optimization techniques as a systematic method of cine phase contrast pulse sequence design and to apply this technique to the measurement of blood flow in vivo. The optimized waveforms were validated in a constant flow phantom with average velocities ranging from 5 to 50 cm/s. The same optimized sequence was also run in a segmented k-space variation with five phase encoding lines per segment. The magnetic resonance (MR) derived velocity measurements were accurate over the entire range of velocities tested (p < .05) in both cases. The same optimized pulse sequence was applied to the measurement of flow in the main pulmonary artery of five normal volunteers and compared with stroke volumes and cardiac outputs calculated from right ventricular volume measurements. These measurements showed a mean difference between the MR phase contrast calculated stroke volume and the volumetric stroke volume measurement of 9.8 ± 11.6%. The mean difference between the calculated phase contrast cardiac output and the volumetric cardiac output was 4.4 ± 10%. These results imply that optimization techniques are an efficient method for designing cine phase contrast pulse sequences.     

In-plane flow velocity quantification along the phase encoding axis in MRI

In-plane flow quantification in MRI offers the potential for assessing vessel patency, and both volume flow rate and flow velocity. These techniques will have definite future impact on MR angiography. The method used in this paper employs motion artifact suppression technique (MAST) gradients to refocus spins travelling along any of the three imaging axes while encoding the velocity component along the phase encoding axis. This method has several advantages over in-plane flow quantification along the read axis. Primarily, flow voids due to complete spin dephasing can be eliminated (or reduced), wider velocity limits can be measured, and gradients can be designed which are sensitive to only velocity along the phase axis with no additional effect from higher order derivatives, or motion along the read axis. Flow phantom studies, carried out on 19 mm inside diameter glass tubes, have produced accurate results for flow rates ranging from 0.6 gallons per minute (GPM) to 2.5 GPM, corresponding to a mean velocity range from 13.2 cm/sec to 55.3 cm/sec. Reynolds numbers varied from 2,700 to 11,500. Errors were less than or equal to 8% over the range of flow rates studied.     

A rotating phantom for the study of flow effects in MR imaging

A common type of phantom used for the study of flow effects in MR imaging is the tube phantom, where a liquid passes through a set of tubes placed in the main magnetic field of an MR scanner. Among the disadvantages with this type of phantom are that a distribution of velocities is present in each tube, and that quantifications of flow effects using tube phantoms may be very time-consuming. In this work, we describe the design and the properties of a rotating wheel flow phantom used for quantification of the effects of flow through the imaging plane as well as in the imaging plane. The proposed phantom is constructed as a rotating gel-filled wheel, surrounded by static volumes filled with the same gel, and the evaluation of the information from rotating and static parts is made with a specially designed computer program. The phantom can be used as a plug flow phantom covering simultaneously an interchangeable velocity interval, which at present has the range −52 mm/s, +52 mm/s. It is shown that the phantom gives adequate information on the dependence of pixel content on first-order motion in MR modulus and phase images. Among the fields of application are rapid calibration of MR imaging units for flow determination using phase information, as well as testing of pulse sequence characteristics and verification of theoretical predictions concerning the flow dependence in MR images.     

Determining uterine blood flow in pregnancy with magnetic resonance imaging

Objective

The purpose of this study is to determine the feasibility of measuring total uterine blood flow in pregnancy using magnetic resonance imaging (MRI) technique.

Methods

Uterine blood flow was determined in pregnant women in whom MRI was being carried out to assess a fetal anomaly. A two-dimensional time-of-flight magnetic resonance (MR) angiogram sequence was performed. Scout images and a peripherally gated phase contrast MR sequence were planned to study simultaneous blood flow in the uterine and ovarian arteries.

Results

The MR pelvic angiogram sequence was completed in 13 women. The uterine arteries were visualized and their cross-sectional area determined. The complexity of the pelvic blood supply prevented the calculation of blood flow velocity and, thus, total uterine blood flow.

Conclusion

The measurement of total uterine blood flow during pregnancy was not possible using our MR technique. The ovarian vessels were not consistently visualized. Doppler ultrasonography remains the best modality by which to estimate total uterine blood flow in pregnancy.     

Direct measurement of quantum phase gradients in superfluid 4He flow

We report a new kind of experiment in which we generate a known superfluid velocity in a straight tube and directly determine the phase difference across the tube's ends using a superfluid matter wave interferometer. By so doing, we quantitatively verify the relation between the superfluid velocity and the phase gradient of the condensate macroscopic wave function. Within the systematic error of the measurement (approximately 10%) we find nu(s)=Planck's/m4InvertedDeltaphi.     

On the accuracy of EPI-based phase contrast velocimetry

For over a decade, echo-planar imaging (EPI) has been used in both the medical and applied sciences to capture velocity fields of fluid flows. However, previous studies have not rigorously confirmed the accuracy of the measurements or sought to understand the limitations of the technique. In this study, a bipolar gradient was added to a flow-compensated EPI pulse sequence to obtain rapid phase contrast images of steady and unsteady flows through two step stenoses. For steady Re = 100 and 258 flows, accuracy was measured through systematic comparisons with CFD simulations, mass flow rate measurements, and spin echo phase contrast images. On average, the EPI image data exhibited velocity errors of 5 to 10 percent, while mass was conserved to within 5.6 percent at each axial position. Compared to spin-echo phase contrast images, the EPI images have 50 percent lower signal-to-noise ratio, larger local velocity errors, and similar mass conservation characteristics. An unsteady flow was then examined by starting a pump and allowing it to reach a steady Re = 100 flow. Accuracy in this case was measured by the consistency between mass flow rate measurements at different axial positions. Images taken at 0.3 s intervals captured the velocity field evolution and showed that 50 to 100 percent errors occur when the flow changes on a time scale faster than the image acquisition time.     

Noninvasive quantification of fluid mechanical energy losses in the total cavopulmonary connection with magnetic resonance phase velocity mapping

A major determinant of the success of surgical vascular modifications, such as the total cavopulmonary connection (TCPC), is the energetic efficiency that is assessed by calculating the mechanical energy loss of blood flow through the new connection. Currently, however, to determine the energy loss, invasive pressure measurements are necessary. Therefore, this study evaluated the feasibility of the viscous dissipation (VD) method, which has the potential to provide the energy loss without the need for invasive pressure measurements. Two experimental phantoms, a U-shaped tube and a glass TCPC, were scanned in a magnetic resonance (MR) imaging scanner and the images were used to construct computational models of both geometries. MR phase velocity mapping (PVM) acquisitions of all three spatial components of the fluid velocity were made in both phantoms and the VD was calculated. VD results from MR PVM experiments were compared with VD results from computational fluid dynamics (CFD) simulations on the image-based computational models. The results showed an overall agreement between MR PVM and CFD. There was a similar ascending tendency in the VD values as the image spatial resolution increased. The most accurate computations of the energy loss were achieved for a CFD grid density that was too high for MR to achieve under current MR system capabilities (in-plane pixel size of less than 0.4 mm). Nevertheless, the agreement between the MR PVM and the CFD VD results under the same resolution settings suggests that the VD method implemented with a clinical imaging modality such as MR has good potential to quantify the energy loss in vascular geometries such as the TCPC.     

Three-dimensional phase contrast MR cerebral venography with zero filling interpolation in the slice encoding direction.

The purpose of this study was to evaluate the magnetic resonance (MR) cerebral venography findings of a three-dimensional phase contrast MR sequence with zero filling interpolation of the data in the slice encoding direction. Fifty volunteers were enrolled in the study. Images were obtained on a 1.5 MR imaging system with acquisition time of 12 min. MIP images were reconstructed throughout the entire imaging volume. A grading scale system was used to assess dural venous sinuses, major deep veins, cortical, and cortical eponymic veins. Inferior group of dural venous sinuses, inferior sagittal sinus, and cortical eponymic veins were poorly demonstrated. Score of the superior sagittal sinus, the straight sinus, the confluence of the superior sinus group, the right transverse and sigmoid sinuses, the internal veins, and the vein of Galen was excellent. The score of the left transverse and sigmoid sinuses was good. In conclusion, when using zero filling interpolation of the data in a three-dimensional phase contrast MR cerebral venography sequence, the superior group of dural venous sinuses and main major deep veins are demonstrated with good conspicuity.     

A method for MR quantification of flow velocities in blood and CSF using interleaved gradient-echo pulse sequences

The aim of this study was to establish a rapid method for in vivo quantification of a large range of flow velocities using phase information. A basic gradient-echo sequence was constructed, in which flow was encoded along the slice selection direction by variation of the amplitude of a bipolar gradient without changes in sequence timings. The influence of field inhomogeneities and eddy currents was studied in a 1.5 T scanner. From the basic sequence, interleaved sequences for calibration and in vivo flow determination were constructed, and flow information was obtained by pairwise subtraction of velocity-encoded from velocity non-encoded phase images. Calibration was performed in a nongated mode using flow phantoms, and the results were compared with theoretically calculated encoding efficiencies. In vivo flow was studied in healthy volunteers in three different areas using cardiac gating; central blood flow in the great thoracic vessels, peripheral blood flow in the popliteal vessels, and flow of cerebrospinal fluid (CSF) in the cerebral aqueduct. The results show good agreement with results obtained with other techniques. The proposed method for flow determination was shown to be rapid and flexible, and we thus conclude that it seems well suited for routine clinical MR examinations.     

Quantitative analysis of PC MRI velocity maps: pulsatile flow in cylindrical vessels.

The accuracy of MR phase contrast (PC) velocity mapping, and the subsequent derivation of wall shear stress (WSS) values, has been quantitatively assessed. Using a retrospectively gated PC gradient-echo technique, the temporal-spatial velocity fields were measured for pulsatile flow in a rigid cylindrical vessel. The experimental data were compared with values derived from the Womersley solution of the Navier-Stokes equations. For a sinusoidal waveform, the overall root-mean-square (rms) difference between the measured and analytical velocities corresponded to 13% of the peak fluid velocity. The WSS derived from the data displayed a 14% rms difference with the analytical model. As an example of a more complicated flow, a triangular saw-tooth waveform was deconstructed into its Fourier components. Velocity maps and the WSS were calculated by the superposition of the individual solutions, weighted by the Fourier series coefficient, for each harmonic. The velocity and experimentally derived WSS agreed with the analytical results (4% and 12% rms difference, respectively). Evaluation of the analytical models allowed an estimate of the inherent accuracy in the measurement of velocity maps and WSS values.     

Advantages of using multiple-echo image combination and asymmetric triangular phase masking in magnetic resonance venography at 3 T

The present work explores the possibility of localizing veins with magnetic resonance venography using susceptibility weighted imaging. It also seeks new approaches, directed by the spatial specificity of activated brain regions, that have sufficient precision for practical use in functional MRI studies. A 3D flow compensated multiple gradient echo sequence, featuring optimized T2* weighting within a reasonable time of acquisition (11 min) and a small voxel size (0.5x0.5x1 mm3), was used to acquire MR images at 3 T. Post-processing consisted of homodyne filtering, linear phase scaling and magnitude masking prior to minimum intensity projection (mIP). The multiple echo approach provided a satisfactory (48+/-7%) increase in signal-to-noise ratio with respect to conventional methods. Specific features of the blood oxygenation level-dependent phase effect were simulated and used for designing and exploring different phase masking methods in relation to vessel morphology and MRI voxel geometry. As with simulations, the best results were obtained with an asymmetric triangular phase masking, featuring an improved venographic contrast without any increase in the full-width at half-maximum. The multiple echo approach provided satisfactory vessel localization capacity by using asymmetric triangular phase masking and a 4-mm-thick mIP. The venographic contrast obtained enabled the detection of vessels with diameter down to approximately 500 microm, suggesting the applicability of the proposed method as an additional technique in fMRI studies.     

MRI of intracranial sinovenous thrombosis: the role of phase imaging

Conventional spin-echo magnetic resonance (MR) imaging of venous thrombosis is complicated by the variable appearance produced by the stage of blood clot degradation and velocity of blood flow. Phase MR imaging is a simple method based primarily on whether protons are stationary or moving. A case of superior sagittal sinus thrombosis demonstrates the utility of phase imaging.     

Three-dimensional phase contrast velocity mapping acquisition improves wall shear stress estimation in vivo

The purpose of this study was to evaluate and apply high-resolution three-dimensional phase contrast mapping for estimation of wall shear stress in vivo. A silicon pipe of 8 mm diameter with a 8.3 ml/s steady flow and the entrance of the carotid bifurcation in 10 young healthy volunteers aged 23.6 +/- 3.1 years was studied. Very high resolution three-dimensional and two-dimensional phase contrast mapping sequences with spatial resolutions of 0.31 x 0.31 x 1.5 mm(3) and 0.31 x 0.31 x 3 mm(3), respectively, were compared in vivo and in vitro. Wall shear stress was calculated using multi-sectored, three-dimensional paraboloid fitting. In comparison to the two-dimensional measurements, the three-dimensional method with only half the slice thickness gave higher signal-to-noise ratio and velocity-to-noise ratios both in vivo and in vitro. Wall shear stress derived from the three-dimensional velocity measurements did not differ from the two-dimensional velocity measurements either in vitro or in vivo. Mean wall shear stress was lowest and oscillatory shear index was highest at the outer wall, towards the carotid bifurcation. Three-dimensional velocity mapping increases resolution and image quality and allows estimation of wall shear stress patterns circumferentially and longitudinally in human arteries.     

Aortoiliac imaging by projective phase sensitive MR angiography: effects of triggering and timing of data acquisition on image quality

To assess the ability of projective phase sensitive magnetic resonance (MR) angiography to visualize the aortoiliac vascular segment, and to determine the effects of triggering and timing of data acquisition om image quality, we studied 18 healthy volunteers, mean age 33.3 +/- 11 years, by color Doppler imaging and by MR angiography. MR angiography was performed at 1.5 T using a flow-adjustable gradient-echo (FLAG) sequence operated in both ECG-triggered and non-triggered acquisition modes. The images were graded in a blinded fashion by two independent observers. The data were analyzed using Pearson's chi-square analysis. Eighteen triggered time-resolved and 17 non-triggered, time-averaged MR angiograms consisting of 252 and 17 angiographic images, (AI) respectively, were analyzed. In the triggered mode 69 (27.4%) AI and in the non-triggered mode 2 (11.8%) AI were diagnostic. At least one triggered diagnostic AI was obtained in each subject. The image grades were not statistically different between observers (kappa = 0.6686). In the triggered mode diagnostic images were acquired within +/- 90 msec of the peak systolic flow velocity determined by Doppler. The proportion of diagnostic images in the triggered mode was highest (73.3%) within a 30-msec interval before the peak flow. In healthy subjects the aortoiliac segment is reliably visualized by FLAG MR angiography. The optimum results are achieved using the triggered acquisition mode and timing acquisition to the initial 180 msec of the abdominal aortic systolic flow pulse.     

油-气-水三相流超声传感器持气率测量*

为探索油-气-水三相流持气率测量难题,该文开展了脉冲透射式超声传感器持气率测量动态实验研究。首先,利用超声传感器与光纤传感器组合,测取了油-气-水三相流中段塞流、混状流、泡状流的响应信号;其次,提取了超声脉冲信号的最大值序列来反映不同流型时超声传感器响应特性,同时,借助双头光纤传感器与相关测速法,计算得到了流体中气泡弦长序列;最后,结合流型与泡径信息,利用超声传感器测量了不同流型下持气率,并分析了不同流型持气率预测的误差来源,为其他油-气-水三相流持气率测量传感器设计提供了借鉴。     

Accurate velocity mapping with FAcE

Spin dislocation between the slice selection, phase encoding, and frequency encoding is a source of image distortions. Two strategies can be pursued to improve the appearance of moving spins in an image. Either the sequence is made equally sensitive to velocity-dependent dislocation artifacts for all spatial directions or the sensitivity is reduced with a shorter echo time. The first approach increases the dislocation for the phase-encoding direction and is therefore not useful if velocity maps with minimal distortion are the goal. FAcE (FID acquired echoes) is a sequence with separate sampling of the left and right k-space half-planes that allows for very short echo times. It was applied for velocity mapping of flow in the slice select direction. Special attention was paid to a compact design of the velocity-encoding select gradient to achieve short echo times even with high velocity sensitivity. Artifacts introduced by in-plane motion were studied for FAcE and conventional gradient-echo sequences, both in phantom experiments and simulation. FAcE allows for very short echo times with inherent motion compensation of the frequency-encoding gradient. Thus, both motion-related dislocation artifacts and signal voids due to coherence loss in regions with irregular flow are minimal.     

Functional Magnetic Resonance Imaging in Intact Plants—Quantitative Observation of Flow in Plant Vessels

Quantitative magnetic resonance (MR) images of flow velocities in intact corn plants were acquired using magnetization-prepared MR microscopy. A phase contrast flow imaging technique was used to quantitate water flow velocities and total volume flow rates in small xylem vessels. The simultaneous measurement of the transpiration of the whole plant was achieved by using a closed climate chamber within the MR magnet. The total volume flow rate and the transpiration values were in close correlation. Functional magnetic resonance imaging in intact plants was performed by light stimulation of the transpiration inside of the magnet. The change in the flow velocities in the xylem vessels of single vascular bundles was in correlation with the changes in the transpiration. Significant differences were observed between the xylem vessels in different vascular bundles. Furthermore, flow velocity measurements were performed on excised plant stems and visualized by the uptake of the MR contrast agent, gadolinium-diethylenetriamine pentaacetic acid (Gd-DTPA). A comparison between the phase contrast flow imaging and the contrast media uptake showed to be in good agreement with each other.     

Cornu's spiral in the Fresnel regime studied using ultrasound: a phase study

A simple experimental technique for measuring the phase and amplitude of diffracting ultrasound wave [A. Hitachi and M. Takata, Am. J. Phys. 78, 678 (2010)] has been applied to diffracting objects with straight edges as a demonstration of the Cornu spiral. Babinet's principle is studied observing the ultrasound field behind a slit and a complementary strip obstacle and is verified directly by comparing vectors (phasors) in the complex plane. The phase of the diffracted wave observed in the geometrical shadow of the straight screen has the form of a cylindrical wave originating at the edge of the straight screen as the boundary diffraction wave proposed by Young. In addition, the incident wave has a phase delay of π/4 behind the wave passing through on the center line of the slit, the plane of symmetry, has been observed as predicted by Huygens-Fresnel diffraction theory.     

Synchronized EPI phase contrast velocimetry in a mixing reactor

Notwithstanding its widespread use in cardiovascular and functional MRI studies, Echo Planar Imaging (EPI) has only recently been subjected to systematic validation studies. Most velocity measurement studies employing such ultrafast MRI methods involve the use of phantoms characterized by rigid or deformable solid motion. The current implementation involves a rotating phantom (angular velocity up to 10.5 rpm) with a superimposed swirling liquid flow (with axial velocities ranging between 0.145 and 0.27 cm/s) of water doped with copper sulfate. The standard implementation of single-shot EPI with phase contrast velocity encoding allows the complete mapping of the Eulerian velocity field in slices perpendicular to the rotation axis following a subtractive procedure requiring the synchronized acquisition of each velocity component on each selected transverse slice during two revolutions of the rotor. The image acquisition time is 100 ms (per velocity component) at each 64 x 64 slice. In addition to acquiring full-field velocity data for future direct comparisons with other techniques, EPI is employed here for the first time to reconstruct the three-dimensional flow field between the blades of a partitioned pipe mixer.