Vectrino+ noise levels

Up to Velocimeters - single point

Vectrino+ noise levels

Posted by Stig Grafsrønningen at June 02. 2010

Hi

I'm using the Vectrino+ with the standard head. I'm trying correct the measurement data to compute Reynolds stresses. The method I'm using assumes that, and only works, if the noise level is the same for all velocity components. However, my measurements show that the noise levels aren't the same at all. The Z1 and Z2 components have a noise level that is one order of magnitude lower than the X and Y components.

Why are the noise levels different? (see attached figure) Is there a post processing of the Z1 and Z2 components?

 

Or:

According to the manual "USER GUIDE" October 2004, rev. C. p. 18

"This means that the Vectrino is more sensitive to the Z-velocity (the component

parallel to the transmit beam) than it is to the X- or Y-velocity. Consequently, the

Z-velocity component yields a lower measurement uncertainty."

 

Is it beacuse of the geometric configuration that the noise level for the Z-velocities are lower? In that case, what I understand is that I cannot use a correction method which assumes an equal noise level for all velocity components.

 

Regards,

 

Stig Grafsrønningen

Attachments

Re: Vectrino+ noise levels

Posted by P.J. Rusello at June 02. 2010

Hi Stig,

Regardless of manufacturer, all acoustic velocimeters equipped a standard head like the Vectrino you are using have lower noise levels for the Z component of velocity because of the head geometry. This is because the Z component is aligned closely with the bistatic axis of the transmitter-receiver pair and is more directly measured than the X or Y components.

The Vectrino also differs from other velocimeters by having a separate channel for each receiver rather than multiplexed sampling where the receivers would share electronics. This is what allows it to sample up to 4x faster than other velocimeters at 200 Hz.

This means each receiver will have different (but very similar) noise levels. Noise levels for each receiver are reported in the header for a data file. Typical levels are 10-12 dB. Check the section called "Velocity Header" for this information.

Finally, the noise for the velocity components can be regarded as white noise (equal energy at all frequencies) and uncorrelated with the noise from another velocity component. Reynolds stress calculations (or any other covariance quantity) should be fairly robust to noise contamination after averaging because of this.

P.J.

Re: Vectrino+ noise levels

Posted by Stig Grafsrønningen at June 04. 2010

Hi

Thanks for the clarification! I guess what I meant is that I need the noise level for each beam (and not velocity) to be the same. A transformation from XYZ to beam coordinates and a FFT should show this.

The Reynols stress calculations should come out fairly easy when taking advantage of the fact that the noise is uncorrelated. I computed the normal stress in the Z-direction (using both signals for Z) for a laminar flow, which should be 0. The variance for these signals were about 1E-6. The computed normal stress was about 1E-18, which translates to machine precision.

However, I've encountered another problem that I hope you could help me with. I have the feeling that I'm doing something wrong when transforming the velocities from XYZ to beam velocities using the transformation matrix. What I've done so far is:

beam_velocities = inv(T)*U(X,Y,Z1,Z2)-->

U[beam1,beam2,beam3,beam4]=inv[1.9949         0   -1.9963         0; 0    1.9741         0   -1.9666;0.5190         0    0.5134         0; 0    0.5022         0    0.5315]*[U,V,W1,W2].  where T is the 4x4 Transformation matrix.

By doing this I get out velocitis that look good,but which velocity is for which beam? I hope you can clarify this for me.

I also tried to compute the vector for each beam in order to find out where the data belonged. However, I must be doing something wrong. It doesn't ad up. By taking:

inv(T)*[1 0 0 0 ; 0 1 0 0; 0 0 1; 0 0 1 0] one migh(?t!??) get vectors pointing from the sampling volume perpendicular on to the beams. Nowing which beam is in the x and y direction, this should give me enough information to find the out which velocity that belongs to which beam. The result from this for my transformation matrix was:  [0.2492         0    0.9690;      0    0.2609    0.9655; -0.2519         0    0.9683;      0   -0.2466    0.9692 ] --> the length of each of these vectors is 1.0 which is a good sign.

Taking the first vector: [0.2492         0    0.9690]  This vector should point in the positive X direction and in the Z direction. (i.e. I know which beam it is)

However, this is were things doesn't ad up. The angle between the Z axis and this vector should be 30° (according to the manual), this is not through because theta=acos(0.9690)= 14.3°. Which is approx. half of the 30° it should be.

Got any tips? I have really no idea.

 

Regards,

Stig Grafsrønningen

Re: Vectrino+ noise levels

Posted by P.J. Rusello at June 04. 2010
Hi Stig,
I don't quite have time to work through your math, but the answer about which beam is which is pretty easy based on the rules of matrix multiplication. For the transformation matrix T, each column represents a beam velocity and each row a velocity component (X,Y,Z1,Z2), so when you multiply T * V_beam where T is 4x4 and V_beam is 4 x number of velocity samples, you get each row of the resultant matrix representing X,Y,Z1,Z2 in that order going down a column. The reverse is also true, so inv(T) * V_xyz results in a matrix 4 x number of samples with each column having rows representing Beams 1,2,3 and 4 in that order.
 
Physically, Beam 1is marked with a red band. I forget which way the numbering goes, but it's pretty easy to figure out. Just place the head in water and cover each beam while looking at the probe check display and you'll see the amplitudes for the appropriate beam drop.
 
Finally, remember that the beam velocities are measured on the bistatic axis (essentially the angle bisector between the transmitter and receiver) which for the Vectrino is 15º (i.e. half of 30º). So, without having checked your math I think you've got things right but were incorrect on the angle you expected.
 
Hope that clarifies things, but feel free to ask any further questions. I haven't had my coffee yet so the first paragraph might be a bit hard to follow.
 
P.J.
 
Powered by Ploneboard
Document Actions
Log in


Forgot your password?
New user?