Abstract
Aiming at turbulent measurements in opaque suspensions, a simplistic methodology for measuring the turbulent stresses with phase-contrast magnetic resonance velocimetry is described. The method relies on flow-compensated and flow-encoding protocols with the flow encoding gradient normal to the slice. The experimental data is compared with direct numerical simulations (DNS), both directly but also, more importantly, after spatial averaging of the DNS data that resembles the measurement and data treatment of the experimental data. The results show that the most important MRI data (streamwise velocity, streamwise variance and Reynolds shear stress) is reliable up to at least \({\bar{r}} = 0.75\) without any correction, paving the way for dearly needed turbulence and stress measurements in opaque suspensions.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Magnetic resonance velocimetry (MRV) has been gaining popularity in engineering design due to its versatility in measuring laminar and turbulent quantities in complex geometries and opaque suspensions. In many engineering applications, measuring turbulent quantities enables design modifications to enhance production and reduce losses. This is crucial in, for example, the pulp and paper industry where not only do fibre suspensions alter turbulent quantities, but the flowfield in many devices used in the paper making process are complex in nature (MacKenzie et al. 2014).
Methods such as ultrasonic doppler velocimetry (UDV) and electrical impedance tomography (EIT) have conventionally been used for measuring the flowfield of opaque suspensions, whereas particle image velocimetry (PIV) and laser doppler velocimetry (LDV) are more commonly used for measuring the flowfield of transparent fluids. An advantage of MRV over other experimental techniques is that it can measure the full 3D velocity field in a few minutes for both transparent and opaque fluids (Caprihan and Fukushima 1990). The limitation of MRV is that in many cases computational fluid dynamics (CFD) simulations or PIV data are needed to provide details about the turbulent velocity fluctuations and Reynolds stresses (Elkins et al. 2009; Iaccarino and Elkins 2006). If MRV is to become a reliable method for measuring the flowfield of opaque suspensions at high Reynolds numbers, where conventional methods are unsuitable, it is important to extend the capacity towards measuring turbulent quantities. In doing so, it should be shown that MRV can produce rms values and Reynolds shear stress that agree with direct numerical simulations (DNS) in a well-defined geometry. This paper details our efforts in utilising phase-contrast MRV to measure the variance and covariance components of the Reynolds stress tensor for water in pipe flow. It will be shown that MRV provides reliable and accurate data in a large portion of the pipe. Works on turbulent flow in fibre suspensions will follow.
Measurements of diffusion with magnetic resonance has a long history (Stejskal and Tanner 1965) and this ability can be used to quantify turbulence through turbulent diffusion. There have been numerous investigations focusing on deriving a method to measure the mean and rms component of velocity using MRV (Dyverfeldt et al. 2006; Gao and Gore 1991; Gatenby and Gore 1994; Li et al. 1994; Newling et al. 2004; Kuethe 1989). Here, the focus will be on methods where the rms value of the velocity is found by assuming turbulent diffusion to be much greater than the molecular self-diffusion of the moving spins. In doing so, turbulence can be modelled as the loss of the net magnetisation through intravoxel phase dispersion. Spins within a voxel which precess at a different phase will be incoherent to each other, reducing the magnitude of the phase-encoded image.
Depending on method, the rms value of velocity is calculated by modelling the phase-variance for different time intervals, or by assuming a phase correlation to the first gradient magnetic moment (see Eq. 3). Both methods are required to assume a spin velocity distribution of a known form. The works of Seymour and Callaghan (1996), Scheven et al. (2005), Caprihan and Seymour (2000) and Jullien and Lemonnier (2012) shows promise towards measuring turbulent diffusivity where the velocity distribution follows a non-Gaussian behaviour, examples of which include the inner layer in fully turbulent pipe flow. By comparing the evolution of the normalised rms widths, Scheven et al. (2005) found that the second and third moments of the perturbed probability distribution for molecular displacements could be found in terms of the measured moments from the original distribution. Protocols aiming at direct measurements of the velocity distribution has been demonstrated to obtain good agreement with DNS results (Jullien and Lemonnier 2012).
The difficulty inherent in applying this type of analysis in fully turbulent flows surrounds the inability to control the displacement field in regions where the velocity distribution is not Gaussian, largely a result of inadequate spatial resolution near the wall. Similar to Seymour and Callaghan (1996); Seymour et al. (2000), an effective diffusivity can be found by limiting the dispersion coefficient to low q-value data, namely where the attenuation function is logarithmically linear with increasing q-values (i.e. Gaussian approximation). In the work of Gao and Gore (1991), a velocity autocorrelation for the phase-variance and a Gaussian phase distribution was used to derive an expression for signal loss due to phase dispersion.
The more recent work of Dyverfeldt et al. (2006) is similar as that of Gao and Gore (1991), however, they assume that a phase-contrast MRI signal in the presence of a first gradient magnetic moment is governed by the spin velocity distribution, magnetic field inhomogeneities and phase accumulation from the velocity field. With the use of two phase-contrast MRI scans, the effect of field inhomogeneities can be eliminated, and the rms of the velocity can be found. Both works use similar approaches towards measuring an rms value of velocity, the difference coming from the assumption that the correction for gradient amplitude duration, separation time, and rise time cancel when two successive scans, which differ in first gradient moment, are compared.
In this work, we adopt the method of Dyverfeldt et al. (2006), a general form for measuring the rms component of velocity with phase-contrast MRV. In the implementation of 2D phase-contrast MRV available to us, the flow encoding vector is limited to only being defined normal to the slice orientation. When measuring quantities such as the Reynolds stresses, multiple encoding directions would be required to isolate the effects of spanwise perturbations in the flow. Without simultaneously encoding the streamwise and spanwise directions (Elkins et al. 2009) it is difficult to acquire quantitative information about the Reynolds stresses due to coordinate system transforms and slice orientation. We discuss here the methodology required to measure the variance and covariance components of the Reynolds stress tensor in a straight pipe with a conventional 2D phase-contrast MRV sequence. We discuss the reliability of this method by comparing our results to DNS data for a range of Reynolds numbers.
2 Theory and methods
2.1 Theory
Magnetic resonance velocimetry (Caprihan and Fukushima 1990; Gladden 1994) works from the basis that hydrogen nuclei are dipoles. In the presence of an external magnetic field, spin packets of hydrogen nuclei are represented by their net magnetisation vector, \({\mathbf {M}},\) which is oriented along the direction of the magnetic field. When the spins are relaxed, hydrogen nuclei will align in the direction of the external magnetic field while precessing at the Larmor frequency, \(\omega _0 = \gamma |{\mathbf {B}}|/2\pi ,\) in the direction of the external magnetic field. Here, \({\mathbf {B}}\) is the magnetic field strength, and \(\gamma\) is the gyromagnetic ratio. The orientation of the spin magnetisation vector, \({\mathbf {M}},\) can be altered by applying radio-frequency (RF) pulses that match the Larmor frequency of the precessing spins. This effectively tilts the spins away from their original alignment, which upon turning off the RF-pulse, will induce relaxation of the spins towards the external magnetic field direction. This relaxation process is described by the Bloch (1946) equation
where T1 and T2 are the spin–lattice and spin–spin relaxation times, respectively. If, for example, the external magnetic field direction is along \({\hat{k}},\) then T1 describes how long it will take for the longitudinal component of \({\mathbf {M}}\) to recover to \(63\%\) of its original value (\(M_0\)). T2 is the time it takes for the transverse component of \({\mathbf {M}}\) to lose \(63\%\) of its original value. As the magnetisation returns to its original alignment, a time varying signal is induced in the receiver coils. This signal is detected in the MR system and is the foundation for magnetic resonance imaging (MRI) and MRV. Note that a short T2 relaxation time, and a long T1 relaxation time would result in a reduction in the signal intensity at shorter times.
To generate 2D images, spatially varying magnetic field gradients are used to encode the phase and frequency of the moving spins. The spins will accumulate phase according to
where \(\phi _0\) is the phase accumulation due to magnetic field inhomogeneities, \({\mathbf {G}}\) is the magnetic field gradient strength, and \({\mathbf {r}}\) is the position. The Larmor frequency of the spins will change according to \(\omega = \gamma (|{\mathbf {B}}|+{\mathbf {G}}\cdot {\mathbf {r}})/2\pi .\) If spins are moving at a velocity parallel to a magnetic field gradient, they will acquire a phase shift which is proportional to their velocity \({\mathbf {u}}.\) The total phase accumulation of a moving spin can then be written as
where, \(\mathbf {M_0}\) and \(\mathbf {M_1}\) are the zeroth and first order gradient magnetic moments, respectively. Conventionally, bipolar flow encoding gradients are used to isolate the contribution of velocity on the total phase accumulation of moving spins. Since the accumulation of phase from static spins is independent of gradient polarity, the velocity of the moving spins can be calculated from the net phase shift between the bipolar gradients. The velocity can then be directly calculated by
where,
2.2 Turbulence MRV measurements
2.2.1 Velocity variance in the direction of the slice
In MRV, intravoxel phase-dispersion within a spin packet reduces the amplitude of the acquired signal from the MR system, which in turn can be modelled as an rms value in velocity. In this work, the method of Dyverfeldt et al. (2006) was adopted, which first starts with considering the signal \(S({\mathbf {k_v}})\) within a voxel for given flow and encoding parameters as
where \(s({\mathbf {u}})\) is the spin velocity distribution, defined in this work to be Gaussian, namely as
where \(\sigma\) is the standard deviation. Equation 6 shows that for phase-contrast MRV, the signal acquired from the MR system within a voxel is governed by magnetic field inhomogeneities and intravoxel phase accumulation. Previous works (Dyverfeldt et al. 2006; Gao and Gore 1991) found that if a velocity distribution is assumed, such as the Gaussian distribution shown in Eq. 7, then \(S({\mathbf {k_v}}) = g(u_{\rm{{rms}}}).\) Substituting Eq. 7 into Eq. 6 results in
By combining the results of two scans with different first gradient moments, the effect of magnetic field inhomogeneities can be eliminated. Since the variant term, \(\sigma ^2,\) is located in the real part of the exponent, the magnitude portion of \(S({\mathbf {k_v}})\) gives
The flow encoding vector can be defined to measure the velocity in any direction on a 3D cartesian coordinate system. If the encoding direction is defined as: \({\mathbf {e}} = x{\hat{i}} + y{\hat{j}} + z{\hat{k}},\) where \(\hat{k}\) is along the streamwise direction z, and \(\hat{j}\) and \(\hat{i}\) are along the spanwise directions y and x, respectively, then Fig. 1 shows the resulting slice orientation for flow encoding vectors \({\mathbf {e_3}} = 0{\hat{i}} + 0{\hat{j}} + 1{\hat{k}}\) and \({\mathbf {e_9}} = \sqrt{2}/2{\hat{i}} + 0{\hat{j}} + \sqrt{2}/2{\hat{k}}\) (remember that our flow encoding gradient can only be normal to the slice).
The tilting of the orientation of the slice to flow encode a straight pipe along a direction in the \({\hat{i}}{\hat{k}}\) or \({\hat{j}}{\hat{k}}\)-plane results in an elliptic image as shown in Fig. 1. With extrinsic matrix rotations the coordinate system defined with \(\mathbf {e_3}\) (i.e. normal to the streamwise direction) can be derived from \(\mathbf {e_9},\) and the Reynolds stress tensor can thus be computed using the measured velocity variances from multiple encoding directions.
2.2.2 Reynolds stresses
The method adopted here for determining the Reynolds stresses was originally used for hot-wire anemometry or single probe LDV measurements (e.g. Inasawa et al. 2003), and has been extended in this work to account for the use of multiple flow encoding directions. To compare the results obtained from MRV to DNS, a relationship between the variance of the measured velocity, flow encoded in a Cartesian system, to the Reynolds stress tensor in cylindrical coordinates can be achieved by decomposing the measured velocity variance to the averaged variance and covariance components of the stress tensor. Table 1 summarises the flow encoding vector directions and the resulting variance and covariance components found from selectively decomposing lines from the 2D images.
To exemplify Table 1, the instantaneous velocity decomposition for a flow encoding vector \(\mathbf {e}\) throughout the cross-section of the pipe is detailed. The tangential position within the pipe is defined as
where y and x define the position of a voxel on the Cartesian coordinate system. The instantaneous measured velocity for a flow encoding vector \(\mathbf {e}\) can then be shown to be
where \(u_z,\) \(u_\theta ,\) and \(u_r\) are the fluctuating components of velocity, and \(e_{\hat{i}}\) and \(e_{\hat{k}}\) are the components of the encoding vector as described in Table 1. A relationship between the variance of the measured velocity, \(\langle u_{M}^2 \rangle ,\) and the Reynolds stress components can be found by expanding Eq. 11, squaring the result and averaging to get
To calculate the variance and covariance components in Eq. 12 from the measured MRV data \(\langle u_{M}^2 \rangle ,\) the set of linear equations, defined by the flow encoding vectors in the \({\hat{i}}{\hat{k}}\)-plane and \(\hat{k}\)-plane (i.e. \(\langle u_z^2 \rangle\)) shown in Table 1, are solved using a least-squares fit. If a line of data along \(\Phi _o = 0\) is selected, the linear system to be solved can be expressed as
where,
If a line of data is selected along \(\Phi _o = \pi ,\) the terms \(\langle u_r^2 \rangle\) and \(\langle u_zu_r \rangle\) in Eq. 13 would be replaced with \(\langle u_\theta ^2 \rangle\) and \(\langle u_zu_\theta \rangle ,\) respectively. In this work, the covariant turbulent stresses \(\langle u_z u_\theta \rangle\) and \(\langle u_r u_\theta \rangle\) were not investigated.
2.3 Extracting measured data from DNS
In order for the data obtained from Eq. 13 to be accurate, the artefacts induced by spatial averaging in combination with tilting of the slice must be limited. Data from two published DNS studies (Ahn et al. 2015; Khoury et al. 2013) is used to evaluate these effects. First, a flowfield of the measured quantities in Eq. 11 for multiple flow encoding vectors is created. An example of the velocity variance, \(\langle u_M^2 \rangle ^{1/2}\) is shown in Fig. 2. As a comparison, the velocity variance measured directly from MRV with an identical flow encoding vector is shown. It can be seen that qualitatively DNS and MRI agree well, where both exhibit the same asymmetric variance over the cross-section of the pipe. The asymmetry is a result of slice obliqueness and averaging along the slice thickness direction.
Using elemental matrix rotations, the 2D field of \(\langle u_{M}^2 \rangle\) from DNS was rotated such that the xy-plane is normal to a defined flow encoding vector. An illustration of the rotation is shown in Fig. 3. To investigate the effect of slice thickness and averaging, the variance of the measured quantity \(\langle u_{M}^2 \rangle\) in the DNS was averaged normal to the plane, similar to the operation of the MRI. This can be expressed as
where \(x_c^\prime\) and \(y_c^\prime\) are the coordinates along the centre of the slice, and \(h(2\alpha )\) is the number of positions along \(z^\prime\) that were averaged. For all numerical simulations, \(\Delta z^\prime\) was set to \(\Delta y^\prime .\)
The DNS data was interpolated to match the MRI grid in the \(x^\prime y^\prime\)-plane by averaging according to
where FOV is the MRI’s field-of-view, and \({\rm{IPR}}_y\) is the y-component of the MRI’s in-plan matrix resolution (see Sect. 2.4). The spatially matched DNS data can then be rotated back into the xy-plane by
where z is zero \(\forall\) (x, y). The variance and covariance components are finally calculated using Eq. 13. Note that a linear interpolation along \(\Phi _o\) was used to match the spatial resolution of the rotated data to that acquired with \(e_{\hat{k}} = 1.\)
2.4 MR system
All of the experiments in this work were performed on a 1T MR system (Aspect Imaging, Shoham, Israel) with a 30 G/cm peak gradient strength. The MR system composes of 3 gradient coils and 1 radio frequency coil. Data was collected using a variety of flow encoding directions where the field-of-view (FOV) ranged from 38 to 56 mm, depending on the orientation of the slice. All measurements were made with water with no contrast agent. The T1 and T2 relaxation times are assumed to be approximately 4 and 2 s, respectively. The in-plane matrix resolution (IPR) was chosen to be \(128 \times 128\) pixels. The PC-MRV sequence, shown in Fig. 4, was set with a slice thickness of 5 mm or \(D_{\rm{pipe}}/7\) and a repetition time (TR) defined as either 10 ms for \(Re_\tau \le 321\) or 6.1 ms for \(Re_\tau \ge 597.\) The echo time (TE) for all Reynolds numbers was 3 ms and the flip angle was 30\(^\circ .\) Note that in this work, \(Re_\tau = \rho D_{\rm{pipe}} u_\tau /\mu ,\) where \(u_\tau\) is the friction velocity, found experimentally via the pressure gradient relation \(\rho u_\tau ^2 = \tau _w = D_{\rm{pipe}}|{\rm{d}}P/{\rm{d}}z|/4.\) Details on the measurements can be found in Sect. 2.5. The density and viscosity of water, \(\rho\) and \(\mu\) respectively, where used for all calculations.
To increase the signal-to-noise ratio, five signal acquisitions, or combined images, were prescribed for each experiment where, depending on Reynolds number (see Table 2), between 25 and 120 experiments were repeated. With a TR time of 10 ms the total scan time for 25 experiments at one bipolar gradient amplitude was 6.2 min. The velocity encode (VENC) for measuring the mean component of flow was set such that the maximum velocity in the phase-unwrapped image was always less than the VENC. The VENC defined for the turbulent stress measurements was chosen to be \(15\%\) of the velocity encode used to measure the mean component of flow at the same \(Re_\tau .\) That particular value was chosen because the DNS data showed that no variant or covariant rms components exceeded \(14\%\) of the centreline velocity, therefore, the resolution in phase variance was increased by choosing a VENC to be \(15\%\) of what was used to measure the mean flow.
2.5 Flow loop
A centrifugal pump (Flygt model no. 040-6716260) provided flow, ranging from 0.11 to 2.5 l/s. The volumetric flow rate was measured using a Krohne aquaflux 070 flowmeter and the pressure gradient was measured using a Fuji electric FCX-AII differential pressure transducer. Both the flow meter and the differential pressure transducer were connected to a NI-9205 DAQ device for data logging. Flow exits the pump through approximately 10 m of flexible 5 cm I.D. tubing, before splitting into two 5 cm I.D. PVC pipes. The split flow enters a single 5 cm I.D. PVC with a 90\(^\circ\) bend before constricting to the 3.5 cm I.D. acrylic test section. The MRV is located more than 140 pipe diameters from the constriction.
Figure 5 shows the friction factor data for a range of Reynolds numbers measured in the test rig. It can be seen that a transition from laminar flow occurs near \(Re \approx 2000,\) where the data agrees well to the Colebrook friction factor line (Colebrook and White 1936), with an assumed roughness height of 1.6 \(\upmu\)m, from \(Re = 4000\) to 10,000. Note that the Reynolds number was experimentally calculated as \(Re = 4\rho Q_{\rm{flow}}/\pi \mu D_{\rm{pipe}},\) where \(Q_{\rm{flow}}\) is the volumetric flow rate measured from the Krohne aquaflux flowmeter. The data suggests that the experimental rig agrees with the typical range for transition from laminar to turbulent flow. The position of the differential pressure transducer was found to have no effect on the measured friction factor, indicating the flow is fully developed for all Reynolds numbers within the MRI test section.
3 Results
The results presented here are broken into two sections: evaluation of the errors inherent in the MRI measurements by utilisation of the DNS data, and results from the MRI. The evaluation of the error was made for \(Re_\tau = 180,\) 360, 550, 1000, and 3000, whereas eight Reynolds numbers were measured using MRV. A summary of the measurements can be found in Table 2.
3.1 Evaluation of errors in the measurements
The error evaluation using DNS data is used to determine the optimal combination of encoding vectors to minimize the effect of spatial averaging along the rotated streamwise direction, \(z^\prime ,\) and secondly, the effect of slice thickness on the deviation from the un-averaged \(\langle u_zu_r \rangle\) and \(\langle u_ru_r \rangle\) components. Nine pairs of encoding vectors were considered, ranging from
to
in increments of \(\pi /36\) (see Table 1). The covariance and variance components, \(\langle u_zu_r \rangle\) and \(\langle u_ru_r \rangle ,\) were calculated using Eq. 13 for multiple combinations of encoding pairs. Figure 6 shows the L2-norm of \(\langle u_zu_r \rangle\) and \(\langle u_ru_r \rangle ,\) normalised to the streamwise centreline velocity \(U_\infty ,\) for every combination of flow encoding vector pairs. This can be expressed as, for example,
In Eq. 15 the difference in the covariance component \(\langle u_z u_r \rangle\) is defined as,
where \(\langle u_z u_r \rangle _{\rm{measured}}\) is obtained by averaging DNS data and is a function of N. Note that \(\bar{r}\) in Eq. 15 is the non-dimensional radius and N is the number of encoding pairs ranging from 1 to 9. The L2-norm here is used as a quantitative measure of the deviation from the original DNS data due to spatial averaging of the MRV for each combination of encoding vectors. It can be seen from Fig. 6 that the L2-norm for the value of \(\langle u_zu_r \rangle\) is always less than \(3\%\) of \(U_\infty ,\) regardless of the number of encoding pairs considered in Eq. 13. The L2-norm in the rms value of \(\langle u_ru_r \rangle ,\) however, was found to deviate up to \(43\%\) of \(U_{\infty }\) if only \(\pm {\mathbf {e}}_9 = \pm \sin \frac{\pi }{4}{\hat{i}} + 0{\hat{j}} + \cos \frac{\pi }{4}{\hat{k}}\) was used. This is a result of the gradient in the rms values along \(z^\prime ,\) which increases with decreasing size of the outer layer (i.e. \((1-r)^+>50)\) or \(Re_\tau .\) Note that \((1-r)^+\) is the distance from the wall measured in viscous lengths, calculated only with DNS as \((1-r)^+ = u_\tau (R-r)/\nu .\)
In Fig. 7 the effect of slice thickness is shown in the diagnostic plot for \(Re_\tau = 180\) and 1000. The result of only spatially condensing the data to match the resolution of the MRI (refer to Sect. 2.3) is illustrated with \(2\alpha = 0,\) where it can be seen that \(| \Delta \langle u_zu_r \rangle ^{1/2}|\) and \(| \Delta \langle u_ru_r\rangle ^{1/2}|\) for \(Re_\tau = 1000\) is greater than with \(Re_\tau = 180\) near \(U/U_\infty = 0.65.\) This correlates well to \((1-r)^+ = 50\) or \({\bar{r}} = 0.95\) for \(Re_\tau = 1000,\) which effectively means that the MRI will likely be unable to measure beyond the peak in \(\langle u_zu_z \rangle\) for high Reynolds number flows purely due to spatial resolution. Increasing the slice thickness to \(2\alpha = D/7\) and 2D/7, it can be seen that the effect of averaging along \(z^\prime\) resulted in \(\Delta \langle u_zu_r \rangle ^{1/2}\) values always below \(1\%\) of \(U_\infty\) for either \(Re_\tau\) equal to 180 or 1000 in the outer layer.
The most profound effect of the slice thickness is on the radial variance \(\langle u_ru_r \rangle\) for \(Re_\tau = 180,\) where \(| \Delta \langle u_ru_r\rangle ^{1/2}|\) varied between 2 and \(8\%\) of \(U_\infty\) within \(U/U_\infty = 0.8\) and 0.6, respectively. As mentioned, this is a result of the outer layer starting at \({\bar{r}} = 0.72\) or \(U/U_\infty = 0.83.\) The radial variance as well diverged from zero near \((1-r)^+ = 50\) or \(U/U_\infty \approx 0.65\) for \(Re_\tau = 1000.\) The error in both \(\langle u_zu_r \rangle\) and \(\langle u_ru_r \rangle\) near \(U/U_\infty = 1\) was found to exist for all Reynolds numbers, and is a result of the steep gradient in all rms components within \(U/U_\infty \in (0.92,0.96).\)
3.2 MRV results
The experimental rig used for this work is primarily used for analysing multiphase cellulose fibre suspension flows, and as a result does not contain a settling chamber to remove swirl from pipe bends and pumping, or oscillations from connections in the piping system. This is because fibre suspensions tend to flocculate and clog settling chambers. Upon measuring the mean tangential velocity profile it was found that the flow of water exhibited relatively high magnitudes of swirl. Figure 8 shows the magnitude of the swirl measured from the MRI in a format similar to the diagnostic plot. A swirl number was defined as, \(S = U_{\theta ,w}/U_\infty ,\) where \(U_{\theta ,w}\) is the maximum swirl component measured nearest to the wall. Conventionally, a swirl number is used for rotating pipe flows, however, it has been adopted in this work to account for the effects of the swirl on the mean streamwise velocity component. Of the Reynolds numbers investigated in this work, it was found that the swirl number ranged from 0.1 to 0.7. This indicates that near the wall there is shear layer from the swirling component, leading to an increase in the friction velocity at lower bulk velocities. To account for the swirl when comparing the rms values measured from MRV to DNS, a scaling factor was introduced on the centreline velocity measured from the MRI. This scaling is justified with the work of Örlü and Alfredsson (2008), who studied the mean and rms of the streamwise component in a swirling jet emanating from a fully developed axially rotating pipe flow. It was found that a rotating pipe of swirl number 0.5 reduced the centreline velocity by almost \(12\%\) at the pipe exit, whereas the rms component was unaffected.
We adopt the outer layer scaling law proposed by Oberlack (1999) for a rotating pipe flow, which can be expressed as
where \(\chi\) is a function of the velocity ratio and \(\varphi\) is a constant. Since the form of \(\chi\) is unknown, and the swirl cannot be controlled for a given \(Re_\tau ,\) a first order empirical swirl-correction for \(U_\infty\) was used. Assuming that \(\varphi = 1\) and \(\chi\) is a non-uniform expansion, chosen to be
with limits,
Substituting \(\chi\) and \(\varepsilon\) into Eq. 16 and rearranging yields
For simplicity, the corrected centreline velocity will be defined as
Figure 9 compares the MRV data in the diagnostic plot, with and without scaling in \(U_\infty ,\) to DNS. It can be seen that without the scaling, the normalised rms values deviate further from the DNS data for Reynolds numbers which display strong swirl (e.g. \(Re_\tau = 160\)). Applying the swirl correction to the rms value improves the agreement to DNS for \(Re_\tau \ge 321,\) however, the correlation is still limited to either the resolution of the MRI, or the onset of the viscous wall region, \((1-r)^+ < 50.\) It must be noted that we make no attempt to model the swirl in either the viscous wall region, or the outer layer. It is assumed that the effect of the swirling component on the distribution of the mean flow is negligible in the outer layer, and that the swirl correction fails for \(Re_\tau = 160\) and 2283 as a result of \({\rm{d}}U_\theta /{\rm{d}}r\) being highly non-linear in the outer layer.
To establish the effect of encoding angles on the calculated variance and covariance components from Eq. 13, the L2-norm was calculated for all possible combinations of encoding pairs ranging from \({\mathbf {e}}_4\) to \({\mathbf {e}}_9.\) Figure 10 presents the minimum L2-norm for each Reynolds number considered in this work, compared against DNS data at an \(Re_\tau\) closest to that measured with MRI. It can be seen that, depending on encoding angles, the L2-norm of \(\langle u_zu_r \rangle\) was below 0.06 for all Reynolds numbers. The encoding angles which minimised the L2-norm for three Reynolds numbers, \(Re_\tau = 321,\) 1115, and 1827 are highlighted. It was found that \({\pm }45^\circ\) or \({\pm }\pi /4\) was a common encoding pair between the Reynolds numbers shown in Fig. 10. It appears that the encoding vector \({\mathbf {e}}_9\) gives leading order accuracy for the covariance component, where higher order accuracy is achieved by acquiring data at encoding angles which approach the direct measurement of the streamwise variance (i.e. \({\mathbf {e}}_3\)).
The L2-norm of the radial variance was found to range between 0.18 and 0.46 for \(Re_\tau = 160\)–2283. It was found that out of all possible combinations, the minimum L2-norm was acquired with only \({\mathbf {e}}_9.\) This is not surprising considering that an encoding angle of 90\(^\circ\) or \(\pi /2\) would directly measure the radial variance, however, only at the centre of the pipe. Based upon the corrected DNS, measuring the radial variance will be a challenge as steep encoding angles result in increased error but approach a direct measurement of the radial variance as \(e_{\hat{i}} \rightarrow 1.\) A reduction in slice thickness would undoubtedly reduce the averaging effects, but decreases the signal-to-noise ratio. It also must be noted that the effects of the swirl on the radial variance are not quantified in this work, which may well be contributing to the divergence from DNS. Limiting the swirl and optimising the MR measurement is required before a quantitative analysis on \(\langle u_ru_r \rangle\) is possible.
Figure 11 presents the DNS and MRV (see Table 2 for marker summary) covariance rms component in the diagnostic format for all the Reynolds numbers considered in this work. The data shown was calculated using the encoding angles which minimised the L2-norm, the values of which can be seen in Fig. 10. The MRV data measured at \(Re_\tau = 321\) and 597 (top figure) was found to agree to DNS results at \(Re_\tau = 360\) (green) and 550 (red), respectively. The middle plot in Fig. 11 shows good agreement for \(Re_\tau = 922\) and 1115 when compared to DNS acquired at 1000. The normalised covariance measured from MRI at \(Re_\tau = 1600\) is slightly discontinuous when compared to DNS near \(U/U_{\infty } = 0.95,\) a possible result of the correction factor failing to account for the swirl. Similarly, Fig. 11 (bottom) indicates that the swirl correction factor \(\xi\) failed for \(Re_\tau = 180,\) 1827, and 2283. The covariance rms measured at \(Re_\tau = 2283\) was significantly lower than the DNS data both with and without the swirl correction factor, an analogous result to the streamwise variance shown in Fig. 9.
Figure 12 presents the radial distribution of the normalised turbulent stress for the MRV cases with minimal swirl. The turbulent stress and position were defined as
where \(u_\tau\) is the friction velocity determined experimentally from the axial pressure gradient measured along the test section. The MRV results are in excelled agreement to the corrected DNS where the discrepancies are largely a result of the uncertainty in wall position. To overcome this uncertainty, it is suggested to diagnose turbulent stress data when using titled slices in the format shown in Figs. 9 and 11.
4 Conclusions
A method for measuring the turbulent stresses using phase-contrast magnetic resonance velocimetry with the flow encoding gradient normal to the slice has been detailed. Data from DNS was used to quantify the effects of spatial averaging from the MRI, where it was found that slice thickness and encoding angle had a much more profound effect on the radial variance rather than the covariance component. Results from MRV were presented in a diagnostic format for the streamwise variance and a covariance component of the turbulent stress tensor. It was found that the streamwise variance agreed well to DNS once a correction factor was applied to account for the swirl within the test section. The swirl is a deficiency in the flow loop since it is optimised to handle suspensions, not to create a perfect single phase flow.
The covariance rms \(\langle u_zu_r \rangle ^{1/2}\) calculated from the MRV data was in good agreement to DNS for Reynolds numbers with minimal swirl without any corrections. The discrepancies between DNS and MRV measurements are likely a result of the additional shear caused by the swirl. Quantifying the angular and axial momentum along the test section would likely rectify swirl effects, however, is outside the scope of this work.
The use of tilted slices was found to be effective in measuring the covariance component \(\langle u_zu_r \rangle\) and the streamwise variance of the turbulent stress tensor for \(r/R < 0.75\) with phase-contrast MRV. Excellent agreement with DNS was obtained in this region.
A particularly interesting application for this method will be measurements in opaque suspensions, where turbulent stress data is crucial for modelling but hard to obtain with optical or intrusive methods.
References
Ahn J, Lee J, Lee J, Kang J, Sung H (2015) Direct numerical simulation of a 30R long turbulent pipe flow at Re\(_\tau\) = 3008. Phys Fluids 27(6):065110–065115
Bloch F (1946) Nuclear induction. Phys Rev 70(7–8):460–474
Caprihan A, Fukushima E (1990) Flow measurements by NMR. Phys Rep 198(4):195–235
Caprihan A, Seymour J (2000) Correlation time and diffusion coefficient imaging: application to a granular flow system. J Magn Reson 144:96
Colebrook C, White C (1939) Experiments with fluid friction in roughened pipes. J Inst Civ Eng I1:133–156
Dyverfeldt P, Sigfridsson A, Kvitting J, Ebbers T (2006) Quantification of intravoxel velocity standard deviation and turbulence intensity by generalizing phase-contrast MRI. Magn Reson Med 56(4):850–858
Elkins C, Alley M, Saetran L, Eaton J (2009) Three-dimensional magnetic resonance velocimetry measurements of turbulence quantities in complex flow. Exp Fluids 46:285–296
Gao J, Gore J (1991) Turbulent flow effects on NMR imaging: measurement of turbulent intensity. Med Phys 18(5):1045–1051
Gatenby J, Gore J (1994) Mapping of turbulent intensity by magnetic resonance imaging. J Magn Reson A 121:193–200
Gladden L (1994) Nuclear magnetic resonance in chemical engineering: principles and applications. Chem Eng Sci 49(20):3339–3408
Iaccarino G, Elkins C (2006) Towards rapid analysis of turbulent flows in complex internal passages. Flow Turbul Combust 77:27–39
Inasawa A, Lundell F, Matsubara M, Kohama Y, Alfredsson P (2003) Velocity statistics and flow structures observed in bypass transition using stereo PTV. Exp Fluids 34:242–252
Jullien P, Lemonnier H (2012) On the validation of magnetic resonance velocimetry in single-phase turbulent pipe flows. J Magn Res 216:101
Khoury GE, Schaltter P, Noorani A, Fischer P, Brethouwer G, Johansson A (2013) Direct numerical simulation of turbulent pipe flow at moderately high Reynolds number. Flow Turbul Combust 91:474–495
Kuethe D (1989) Measuring distributions of diffusivity in turbulent fluids with magnetic-resonance imaging. Phys Rev A 40(8):4542
Li T, Seymour J, Powell R, McCarthy K, Ödberg L, McCarthy M (1994) Turbulent pipe flow studied by time-averaged NMR imaging: measurements of velocity profile and turbulent intensity. Magn Reson Imaging 12(6):923–934
MacKenzie J, Martinez D, Olson J (2014) A quantitative analysis of turbulent drag reduction in a hydrocyclone. Can J Chem Eng 92(8):1432–1443
Newling B, Poirier C, Zhi Y, Rioux J, Coristine A, Roach D, Balcom B (2004) Velocity imaging of highly turbulent gas flow. Phys Rev Lett 93(15):154503
Oberlack M (1999) Similarity in non-rotating and rotating turbulent pipe flows. J Fluid Mech 379:1–22
Örlü R, Alfredsson P (2008) An experimental study of the near-field mixing characteristics of a swirling jet. Flow Turbul Combust 80:51–90
Scheven U, Verganelakis D, Harris R, Johns M, Gladden L (2005) Quantitative nuclear magnetic resonance measurements of preasymptotic dispersion in flow through porous media. Phys Fluids 17:117107
Seymour J, Caprihan A, Altobelli S, Fukushima E (2000) Pulsed gradient spin echo nuclear magnetic resonance imaging of diffusion in granular flow. Phys Rev Lett 84(2):266
Seymour J, Callaghan P (1996) “Flow-diffraction” structural characterization and measurement of hydrodynamic dispersion in porous media by PGSE NMR. J Magn Reson A 122:90
Stejskal E, Tanner J (1965) Spin diffusion measurements: spin echoes in the presence of a time-dependent field gradient. J Chem Phys 42:288
Acknowledgements
The Bo Rydin Foundation for Scientific Research is thanked for funding of the project. The Nils and Dorthi Troëdsson Foundation for Scientific Research supports AS as adjunct professor at KTH.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
MacKenzie, J., Söderberg, D., Swerin, A. et al. Turbulent stress measurements with phase-contrast magnetic resonance through tilted slices. Exp Fluids 58, 51 (2017). https://doi.org/10.1007/s00348-017-2328-8
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00348-017-2328-8