Measurement uncertainty analysis of field-programmable gate-array-based, real-time signal processing for ultrasound flow imaging

Research in magnetohydrodynamics (MHD) aims to understand the complex interactions of electrically conductive fluids and magnetic fields. A promising approach for investigating complex instationary flow phenomena are lab-scale experiments with low-melting alloys. They require a noninvasive flow instrumentation for opaque liquids with a high spatiotemporal resolution, a low velocity uncertainty and a long measurement duration. Ultrasound Doppler velocimetry can achieve multiplane, multicomponential flow imaging with multiple linear ultrasound arrays. However the average raw data output amounts to 1.2GBs−1 at a frame rate of 33Hz in a typical configuration for 200 transducers. This usually prevents long-duration measurements when offline signal processing is used. In this paper, we propose an online signal-processing chain for pulsed-wave Doppler velocimetry that is tailored to the specific requirements of flow imaging for lab-scale experiments. The trade-off between measurement uncertainty and computational complexity is evaluated for different algorithmic variants in relation to the Cramér–Rao bound. By utilizing selected approximations and parameter choices, a prepossessing could be efficiently implemented on a field-programmable gate array (FPGA), enabling a typical reduction of the data bandwidth of 6.5 : 1 and online flow visualization. We validated the performance of the signal processing on a test rig, yielding a velocity standard deviation that is a factor of 3 above the theoretical limit despite a low computational complexity. Potential applications for this signal processing include multihour flow measurements during a crystal-growth process and closed-loop velocity feedback for model experiments.


Introduction
Many important industrial processes, such as continuous steel casting and photovoltaic wafer production, involve metal or semiconductor melt flows. The quality of the product and the energy efficiency of the process strongly depends on the flow behavior of the liquid (Müller and Friedrich, 2010;Gardin et al., 1995;Yasuda et al., 2007). A noncontact way of influencing the flow of electrically conductive melts is the application of magnetic fields that introduce Lorentz forces to the fluid. Investigating the interaction of a magnetic field with the flow pattern and optimizing the spatiotemporal structure of the magnetic field for different applications are subjects of ongoing research in magnetohydrodynamics (MHD). Besides numerical simulations, low-temperature, model-scale experiments are important tools for MHD investigations (Eckert et al., 2007b). They often require advanced flow instrumentation for visualizing complex and instationary flows in opaque liquids. A typical set of requirements for MHD research are as follows: -Noninvasiveness -the influence of the instrumentation to the flow should be negligible (Eckert et al., 2007a).
-Flow imaging capability -the fluid's velocity should be visualized in multiple planes (2D) with two or three ve-locity components (2c or 3c) in order to adequately represent complex flow patterns.
-Spatial resolution -the relevant flow structures have to be resolved, typically in the range of 10 mm (Timmel et al., 2011).
-Long measurement duration -flow phenomena on different timescales should be adequately captured; for instance, rapid spontaneous changes of the flow regime in a rotating flow (Galindo et al., 2017) or in multihour model experiments of the semiconductor crystallization process (Thieme et al., 2017).
-Capability of near-wall measurements -in typical MHD experiments, the metal melt is contained in a vessel. The vicinity of the wall is especially important because the Lorentz force is often concentrated in this region. Contrary to, for instance, medical applications, the walls can be seen as completely stationary in most cases.
-Online capability -conducting long-running MHD experiments requires the ability to examine the data during the duration of the measurement. Some model experiments in the semiconductor crystallization process even benefit from an active control of parameters, like magnetic field intensity and temperature gradient, based on the feedback from online velocity data to stabilize the flow (Thieme et al., 2017).
A measurement system for flow mapping of opaque liquids, namely the ultrasound array Doppler velocimeter (UADV; Nauber et al., 2013a, b), was presented in previous publications. It extends the pulsed-wave Doppler principle (Takeda, 1986;Baker, 1970) by employing multiple linear sensor arrays to achieve multiplane, two-componential flow imaging. The sensors are designed to achieve a lateral resolution of ≈ 3 mm in Galinstan (GaInSn). A combination of spatial-and time-division multiplexing allows one to parallelize the scanning process for a planar velocity map; hence increasing the temporal resolution compared to a strict sequential scan. However, online processing of the data for 200 transducer elements simultaneously on 32 channels at a temporal resolution typically of 33 Hz overburdens PC-based hardware with 1.2 GBs −1 . Therefore, only discontinuous offline measurements could be performed with a limited duration of a few seconds. This severely impedes the usability of the UADV in the context of MHD experiments and restricts the investigations into stationary or periodic flows.
Although several investigations on the measurement uncertainty of Doppler velocity estimation methods for laserbased instrumentation (Fischer et al., 2010), for flow-rate measurements in a pipe (Furuichi, 2013), and for blood-flow measurements in the human body (Lovstakken et al., 2007) have been performed, no comprehensive measurement uncertainty budget in the context of instrumenting an MHD experiment has been presented to the knowledge of the authors.
This paper provides a signal-processing chain that is tailored to the specific requirements of MHD model experiments and shows a real-time implementation using a fieldprogrammable gate array (FPGA). It enables the UADV system to perform long-duration measurements with high frame rates and online flow visualization. Furthermore, we evaluate the measurement uncertainty of the whole UADV system in the context of MHD experiments and present an uncertainty budget according to the methodology proposed by the "Guide to the expression of uncertainty in measurement" (GUM; JCGM, 2008) for a typical configuration.

Measurement principle
In pulsed-wave ultrasound Doppler velocimetry (PW-UDV), short bursts are emitted periodically with a pulse repetition frequency f PR (Baker, 1970). The emission times t s = n b /f PR span the so-called slow-time axis t s , with n b = 0. . .N EPP being the bursts number. The emitted bursts usually consist of N periods periods of a sinusoidal wave, with the frequency f 0 . As the bursts travel through the fluid, scattering particles reflect a fraction of the signal back to the ultrasound transceiver. The received echo signal z(t f , t s ) is acquired, starting from the emission time along the fast-time axis t f . Figure 1 depicts an example of the echo signal for a single moving scattering particle.
The movement of a scattering particle leads to a phase shift of the echo signal between multiple burst emissions (Kasai et al., 1985). The mean phase shift per time unit, expressed as mean frequency f d , is related to the velocity v for a given speed of sound c by the following: with f tx denoting the mean frequency of the received signal burst and c v. The mean phase shift per time unit f d can be interpreted as a Doppler frequency shift f d (Kasai et al., 1985); hence the name Doppler velocimetry. The time since the burst emission t f corresponds to the distance d between the scattering particle and the transducer according to the following equation: This allows a spatially resolved flow measurement along the axis of the transducer, given that the scattering particles follow the motion of the fluid with negligible slip. The axial resolution can be estimated with the following (Jensen, 1996): An example of a single particle moving away from the transducer is given.
The lateral resolution x is given by the width of the ultrasound beam, which is a result of the transducer geometry, the frequency f 0 , and the speed of sound c in the fluid. The temporal resolution t of the velocity measurement is determined through the following:

Ultrasound array Doppler velocimeter
The ultrasound array Doppler velocimeter (UADV) is a modular research platform developed at the Laboratory of Measurement and Sensor System Technique (MST) for flow imaging in opaque liquids with PW-UDV. It is flexible and especially well suited for instrumenting a wide range of experiments in the field of MHD. The hardware of the UADV consists of individually configurable modules driving 25 ultrasound transducers each. It can be scaled to support up to 200 transducers in various configurations; for instance, in four linear arrays which can be individually parameterized regarding ultrasound frequency, pulse shape and length, and pulse-repetition frequency (Nauber et al., 2016;Büttner et al., 2013). A module of the UADV consists of an arbitrary function generator and a power amplifier for generating parameterizable burst signals which are routed through a programmable switching matrix and a transmit/receive switch to the transducers. The received echo signals are amplified with a para-metric gain and routed to the digitization unit. A single microcontroller-driven control unit provides the overall synchronization and the communication with the host PC. Using a combined spatial-and time-division multiplexing scheme, an ultrasound transducer array can scan a measurement plane at higher rates than a strict sequential scan. The UADV supports four independent digitization channels per module. The detailed description of the measurement system is given in Nauber et al. (2016).

Overview
The signal processing for PW-UDV can generally be classified into wideband and narrowband techniques; a comprehensive comparison is given by Torp et al. (1993). While parts of the signal processing in the radio frequency (RF) band can be realized in analog circuitry (Shung, 2015), fully digital implementations have found widespread use in the last decades because of the availability of fast digitizers and the increased flexibility and robustness of such approaches. To simultaneously handle a high number (e.g., 32) of channels through a fully digital signal-processing chain, very large data bandwidths have to be processed. This can be achieved by utilizing the parallel-processing capability of a field-programmable gate array (FPGA). Especially narrowband algorithms are very suitable for FPGA-based implementations, due to their low computational complexity (Alam and Parker, 2003;Loupas et al., 1995a). Therefore, this paper focuses on investigating the most common narrowband method, the velocity estimator by Kasai et al. (1985) and the extensions proposed by Loupas et al. (1995b).
A typical narrowband signal-processing chain is shown in Fig. 2. In this fully digital realization, the slow time t s is sampled for each burst n b at t s = n b 1 f PR and the fast time t f is sampled with a frequency f s as follows: The signals are then bandpass filtered to reduce noise contributions outside of the bandwidth of the transmitted ultrasound signal. A quadrature demodulation is performed, consisting of a Hilbert transform and a subsequent down sampling. Static echoes are removed through a clutter reduction filter (CRF) and the velocities are estimated by an autocorrelation.

Quadrature demodulation
In order to meet the assumptions of the narrowband signal processing and to reduce the influence of noise, a bandpass filtering is performed as follows: with the filter coefficients c i . In order to maximize the SNR for signals with additive white Gaussian noise, a matched filter is used (Turin, 1960) as follows: with the transmitted signal s tx with N tx samples. The result of the quadrature demodulation is a complex signal h unfilt (k, n b ) in the baseband, which can be sampled at a lower rate (reduction by a factor of n sub ) than the raw signal, as follows: with N b = 0, 1, . . . N epp , k = 0, 1, . . . K and the Hilbert transform signalẑ (k, n b ) (with 90 • phase shift with respect to z ).

Clutter-reduction filtering
A common problem of ultrasound Doppler flow measurements is distinguishing between static echoes originating from the walls (the so-called clutter) and echoes originating from scattering particles. Multiple reflections from the transmitted burst inside the wall superimpose the signal from scatter particles in the vicinity of the wall. For this problem, a multitude of signal-processing methods were proposed, most of them based on digital filters (finite impulse response (FIR) or infinite impulse response (IIR) filters) with various initialization techniques (Lee et al., 2009). With these methods, the clutter is distinguished from the particle echoes by a velocity close to zero, respectively, by a Doppler frequency shift close to zero. Because filtering will influence the spectrum of the signal, a bias may be introduced to the subsequent velocity estimation, depending on the frequency cutoff. For typical MHD experimental setups, the wall can be assumed to be completely stationary (in contrast to, e.g., medical applications where clutter is often constituted by slowly moving tissue; cf. Jensen, 1996); therefore, a steep cutoff at a frequency of zero is desirable. The simplest and computationally most efficient approach is to filter the constant component of the demodulated IQ signal by subtracting its mean value, which is the equivalent of applying a very narrowband, high-pass filter (Thomas and Hall, 1994;Jensen, 1996;Torp, 1997;Bjaerum et al., 2002) as follows: As the filter is noncausal, all N EPP samples have to be acquired before the result can be computed.

One-dimensional autocorrelation algorithm
A widely used approach for velocity estimation is the autocorrelation method proposed by Kasai et al. (1985), which operates solely in the domain of IQ-demodulated echo signals and therefore can be implemented very efficiently (Alam and Parker, 2003). It uses the properties of the signals' discrete autocorrelation function as follows: where its values at a lag of 1 relate to the center of mass of the signal's power density spectrum through the Wiener-Khinchin theorem. As shown by Kasai et al. (1985) and Jensen (1996), the mean Doppler shift f d can be approximated through evaluating the autocorrelation function at a lag of n b = 1 slow-time samples as follows: This autocorrelation computation can be expressed solely by repeatedly multiplying accumulate operations and therefore can be implemented very efficiently. Kasai's method approximates the center frequency f d of the received signal with the frequency of the emitted signal as follows: Being based on a phase estimation, the Kasai algorithm is inherently limited in the maximum measurable velocity. Given the 2π -phase ambiguity in Eq. (11), the measurable velocity range resulting from Eq. (1) is (Jensen, 1996) as follows:

Two-dimensional autocorrelation algorithm
An extension of Kasai's autocorrelation method is proposed by Loupas et al. to improve its performance in the following two regards (Loupas et al., 1995b): 1. The assumption of an unchanged center frequency of an ultrasound burst throughout emission, propagation inside the fluid and reception is discarded. This allows one to account for the effect of frequency-dependent attenuation, which is present in most relevant fluids. By explicitly estimating the center frequency of the received signal, a systematic velocity error stemming from the relationship in Eq. (1) v ∝ 1/f tx is avoided.
2. An information loss occurs if only a narrow-band part of a broadband echo signal is processed. Hence a better estimation of the velocity is achieved by including a larger part of the signal spectrum.
Both aspects are addressed by increasing the dimensionality of Kasai's autocorrelation; instead of just correlating along the slow-time axis, a 2D autocorrelation along the slow-and fast-time axis is performed. An autocorrelation with a lag of one fast-time sample yields the estimate of the center frequency as follows: Furthermore, the estimation of the frequencies f tx and f d can be performed using M samples per gate, as follows: The extension of the Kasai autocorrelation algorithm potentially improves the estimation performance while still preserving a low computational complexity.

Online-capable, FPGA-based signal-processing implementation
In order to provide online capability, the signal processing depicted in Fig. 2 has been realized on an FPGA (NI PXIe-7965R; National Instruments, Austin, Texas, USA). The FPGA communicates with a host PC through a peripheral component interconnect express (PCIe) bus and has the ability to stream data through direct memory access into the main memory of the PC. The amplified echo signals (U PP, max = 1 V) are digitized through an A/D converter module (NI-5752; National Instruments, Austin, Texas, USA) for n ch = 32 channels at an externally provided sampling rate 32 MHz > f s > 50 MHz with a quantization of 12 bit. The raw data rate r ADC at this stage is as follows: r ADC = n ch · n sampbytes · f frame · N EPP · N sw · N gates · K. (17) Data are processed as signed 16 bit integer (n sampbytes = 2 B) and, for a typical configuration as listed in Table 1, the data rate amounts to 1.2 GBs −1 .
This data bandwidth is hardly suitable for continuous streaming to a storage device over a long duration (> 1 h) with common PC hardware. Therefore, raw data are only briefly retrieved for debugging purposes or for low framerate measurements and are otherwise not transferred to the host.
The signal-processing steps that perform an IQ demodulation (bandpass filtering, Hilbert transform and down sampling) are significantly reduced in their computational complexity by fixing the ratio of the sampling frequency f s to the ultrasound center frequency f 0 at f tx /f s = 1/4. The matched filter can be realized for a sinusoidal transmit signal at f tx with N periods periods, assuming f tx ≈ f rx with only trivial filter coefficients c i , as follows: This allows one to implement the filtering without multiplication operations, only negations and additions are needed. To provide a low computational complexity approximation of the Hilbert transform for a narrowband case, a fixed time delay can be employed (Kantz et al., 2012) as follows: where f tx /f s = 1/4. The signal processing up to this point contains just the summation, negation and storage primitives and therefore can be implemented on an FPGA with modest resources. The data rate r IQ at this stage for a typical configuration is given by the following: Through the data reduction of 6.5 : 1, the data rate at this stage is r IQ = 185 MBs −1 for a typical configuration, as listed in Table 1. A continuous data streaming to a storage device can be sustained for a long duration at this rate.

Theoretical limit of measurement uncertainty
In order to characterize the performance of a signalprocessing algorithm, it is not only helpful to have relative data compared to other algorithms but also to relate it to a fundamental limit of attainable precision. This absolute limit of uncertainty can be provided by means of the estimation theory using the Cramér-Rao bound (CRB; Radhakrishna Rao, 1945;Cramér, 1946). Given a suitable signal model, the CRB represents the lowest possible variance for estimating a parameter from the signal with an unbiased estimator. In the following, a simple signal model for a discrete time idealized ultrasound echo is described and a derivation of the CRB for velocity estimation is given. A simple approximation of the ultrasound echo signal realizations x(k, n b , θ , σ 2 n ) consists of a sinusoidal signal s(k, n b , θ ) superimposed with additive white Gaussian noise n(k, n b , σ 2 n ) sparsely and periodically sampled in the fast-(k) and slow-time (n b ) axis as follows: with and A being the amplitude of the scattering particles' echo, ϕ 0 a constant phase, and n(k, n b , σ 2 n ) Gaussian white noise, with a variance σ 2 n and zero mean. The unknown quantities are as follows: The CRB provides the lower boundary for the variance of an estimatorθ i according to the inequality, as follows: with I(θ ) being the Fisher information matrix, as follows: Kay (1993) provided a formula for the case when the probability density function p(x, θ ) of the signal model x, Eq. (22), is a Gaussian joint probability function as follows: The differentiation of s(k, n b , θ ) with respect to the unknown quantities is performed analytically, while the matrix inversion was performed numerically using MATLAB (The MathWorks, Inc., Natick, Massachusetts, USA). The resulting CRB for the velocity uncertainty as a function of the signal-to-noise ratio (SNR) is given in Fig. 4c-d. It has a slope of −20 dB/decade, which is consistent with the CRB of other Doppler-based signal-processing problems (Fischer et al., 2010;Chan et al., 2012;Demirli and Saniie, 2001).

UADV measurements on a reference experiment
For an experimental characterization of the measurement performance of the UADV system, a test rig based on the linear translation of a single scattering object is used (Fig. 3). It consists of a linear stage (41.121.102E; OWIS GmbH, Staufen, Germany) that is mounted over a glass tank with the dimensions of 212 × 81 × 135 mm 3 . It moves a scattering object (glass fiber with a spherical tip, and diameter of 0.6 mm, mounted in a hollow needle) with a constant velocity through water (ϑ = 20 • C; c = 1480 ms −1 ). The ultrasound sensor array is mounted on the front wall of the tank and therefore insonates through an 8 mm glass wall and a water-based ultrasound couplant.
In order to trace back the measurement results of the UADV to the definitions of the respective units in the SI system, a simultaneous measurement of the relative position and velocity was done with a vibrometer (OFV-503; Polytech, Waldbronn, Germany; displacement decoder DD-900 and velocity decoder VD-09). A retroreflective tape (3M Scotchlite) was attached to the shaft of the scattering object's mount. For a velocity set point of 10 mms −1 , a standard deviation of the velocity σ v, ref, rel = 0.178 % was determined for the linear stage-vibrometer combination (for the same averaging time as the UADV system).
A total of 130 measurement cycles were conducted, consisting of a constant translation away from the front wall of the tank with a velocity set-point v ref = 10 mms −1 and the respective backward motion. Of the continuously obtained UADV measurements, only those that originate from two defined positions near to and far from the wall during the movement away from the ultrasound transducer (Richter Sensor and Transducer Technology, Germany) are selected in the postprocessing. The clutter-to-signal ratio (CSR) is CSR 1 = −7.3 dB and CSR 2 = −19.0 dB, respectively. To ensure a common time base for vibrometer and UADV measurements, the trigger signal of the UADV is acquired simultaneously with the velocity and position signals. To test the performance under different SNR conditions, white Gaussian noise was added to the raw digitized signals to achieve SNR = −6, −3, . . . , 12 dB. Four algorithmic variants were compared as follows: -(DEF) -the 1D Kasai velocity estimator without clutter filtering, as described in Sect. 3.4 -(CRF) -the 1D Kasai velocity estimator with a clutter filtering according to Sect. 3.3

Number of repetitions
) -the 2D velocity estimator as described in Sect. 3.5 with clutter filtering but without an estimation of f tx -(CRF 2D RF) -the 2D velocity estimator as with clutter filtering including the estimation of f tx .
The parameterization of the experiment and of the algorithms is listed in Table 2. Figure 4 shows the relative systematic deviation from the reference velocity and the relative velocity standard deviation of the tested algorithms. For the low-CSR case (far from the wall) at SNR = 12 dB, it can be seen that a slight negative bias of (DEF) is turned into a positive bias through clutter filtering (CRF) and (CRF 2D). This is compensated by the frequency estimation of (CRF 2D RF), which shows the lowest deviations of all variants for SNR ≥ 3 dB. For the high-CSR case (near to the wall), the variant (DEF) without clutter filter has a relative deviation of v/v ref < −42 %. Through clutter filtering this strong negative bias is turned into a positive bias, which increases with lower SNR for (CRF) and (CRF 2D). The RF estimation of (CRF 2D RF) gives the lowest systematic bias for SNR ≥ 0 dB. The relative standard deviations of all variants of the Kasai's algorithm do not reach the CRB for the given signal model, which is consistent with the findings of Chan et al. (2012). The lowest standard deviations are consistently provided by the variants (CRF 2D) and (CRF 2D RF), which come as close as a factor of 3 to the CRB by using more samples per gate than (DEF) and (CRF).    (Fig. 4) for an SNR of 5dB; uniformly distributed in the interval v rel = ±0.08% Value of the speed of sound of the fluid c GaInSn = 2740 ms −1 (given by Morley et al., 2008, without a measurement uncertainty) Type A estimation based on Proffit and Carome (1962); σ c, rel = 0.03 % 0.03 % Unknown systematic variations of the speed of sound in the fluid due to temperature changes Type B estimation for the sound-speed temperature coefficient of liquid gallium (Proffit and Carome, 1962;Popel et al., 2005); Influence of the spatial resolution from the finite width of the sound field Type B estimation for a beam width x = 3 mm and typical velocity gradients of MHD experiments estimated from numerical simulation; dv dx = 0.16 mm −1 · v max (Galindo et al., 2017); v = 1/2 · dv dx · b v rel = ±24 % 13.9 % Total uncertainty i σ 2 v, rel, i = 13.9 % For the given experimental data, the algorithm variant (CRF 2D RF) provides a suitable trade-off between systematic and standard deviation and computational complexity.

Measurement uncertainty budget of the UADV in liquid metal
A measurement uncertainty budget according to the GUM (JCGM, 2008) is used to assess the contributions of measurement uncertainty for the UADV system. Based on Eq.
(1), the measurand v is derived from the quantities f d , c and f rx . Furthermore, the direct influence of the spatial averaging over the flow within the ultrasound beam width is considered. In Table 3, the uncertainties contribution of these quantities are given for a typical MHD experiment.
For the uncertainties of f d and f rx , the results of Sect. 5.2 are transferred from the reference experiments in water to typical measurement conditions in low-melting liquid metals. The maximum relative systematic deviation and standard deviation for both investigated CSR and a typical SNR of SNR = 5 dB are used to calculate the equivalent uncertainty of the velocity. The influence of an uncertainty in the fluid's speed of sound, c, is estimated by the uncertainty of the measurement of this quantity, in the literature and the temperature dependence, assuming a temperature gradient of ϑ = ±10 K. The uncertainty arising from spatial averaging through the ultrasound beam characteristics is calculated by assuming a lateral averaging of x = 3 mm and velocity gradients of numerical simulations of typical MHD experiments (Galindo et al., 2017).
It can be seen that the biggest contribution to the velocity uncertainty of the UADV measurement system for typical MHD settings with σ v, rel = 13.9 % stems from the spatial averaging over lateral resolution given by the ultrasound beam width of the unfocused transducers. This provides the most promising starting point for further improvements regarding the measurement uncertainty of the UADV system. Furthermore, it justifies the approximations taken for computationally efficiently implementing the signal processing, even though lower uncertainty algorithms exist that approach the CRB (Chan et al., 2012) because signal processing is not the limiting factor in the measurement uncertainty budget.

Example of liquid metal flow imaging
To demonstrate the capabilities of the ultrasound array Doppler velocimeter (UADV) with the proposed signal processing, it is applied to a simple MHD experiment. A cubic vessel with the dimensions of 67×67×67 mm 3 is filled with GaInSn and a 25-element linear transducer array (Richter Sensor and Transducer Technology, Germany) is attached to insonify the central horizontal plane (cf. Fig. 5a). With the application of a horizontally counterclockwise rotating magnet field, a counterclockwise central vortex forms. The UADV measures the velocity component along the axis of the transducers (d axis) with the parameterization given in Table 2 and with f PR = 200 Hz. The resulting planar flow image, using signal-processing variant (CRF 2D RF), is shown in Fig. 5b and c.

Conclusions
Experimental research in the field of MHD can benefit from online, noninvasive flow imaging for investigating fundamental phenomena, such as flow instabilities and optimizing industrial processes. We describe an online-capable signal processing for pulsed-wave Doppler velocimetry that is tailored to the specific requirements of lab-scale model experiments. It is based on a 2D autocorrelator, which allows for a reduction of systematic and stochastic errors through explicitly estimating the RF and utilizing multiple samples per gate. We optimized the signal processing for low computational complexity and implemented substantial parts on an FPGA. A typical reduction of the data bandwidth of 6.5 : 1 enables continuous data streaming to PC hardware.
We evaluated the performance of the implemented signal processing in a water test rig with a single scattering object and a reference velocity obtained through a laser vibrometer. Two different clutter signal levels emulate a measurement close to and far from a wall. A velocity standard deviation of σ v, rel = 0.4 % was found, which is about 3 times the fundamental limit of the uncertainty, the CRB, for velocity estimation. The systematic deviation is v rel = ±0.08 %.
We investigated the measurement uncertainty budget for flow velocity measurements in a typical MHD experimental setup for the low-melting alloy GaInSn. The total measurement uncertainty of σ v, rel = 13.9 % almost solely stems from the effect of spatial averaging over the lateral resolution for flows with high-velocity gradients. This justifies the approximations taken for lowering the computational complexity of the signal processing.
A measurement uncertainty budget of a typical MHD experiment at laboratory scale suggests improvements towards a better lateral resolution. In the context of flow imaging, this can be provided by the focusing and steering of the ultrasound beam using the phased-array principle.
The presented signal processing enables online, multiplane flow visualization with the UADV research platform. A long measurement duration (> 1 h), combined with a high frame rate (> 10 Hz), allows one to investigate complex, instationary flows such as instability phenomena in cubes.
Data availability. Research data are available upon request from the authors.