Visualization of Flow Field Around a Vibrating Pipeline Within an Equilibrium Scour Hole

Summary

The goal of the protocol is to enable visualization of the detailed flow fields and determination of the near-boundary shear and normal stresses within an equilibrium scour hole induced by a vibrating pipeline.

Abstract

An experimental method is presented in this paper to facilitate visualization of the detailed flow fields and determination of the near-boundary shear and normal stresses within an equilibrium scour hole induced by a vibrating pipeline. This method involves the implementation of a pipeline vibration system in a straight flume, a time-resolved particle image velocimetry (PIV) system for pipeline displacement tracking and flow fields measurements. The displacement time-series of the vibrating pipeline are obtained by using the cross-correlation algorithms. The steps for processing raw particle laden images obtained by using the time-resolved PIV are described. The detailed instantaneous flow fields around the vibrating pipeline at different vibrating phases are calculated by using a multiple-time-interval cross-correlation algorithm to avoid displacement bias error in the flow regions with a large velocity gradient. By applying the wavelet transform technique, the captured images that have the same vibrating phase are accurately cataloged before the phase-averaged velocity fields are obtained. The key advantages of the flow measurement technique described in this paper are that it has a very high temporal and spatial resolution and can be simultaneously used to obtain the pipeline dynamics, flow fields, and near-boundary flow stresses. By using this technique, more in-depth studies of the 2-dimensional flow field in a complex environment, such as that around a vibrating pipeline, can be conducted to better understand the associated sophisticated scour mechanism.

Introduction

Subsea pipelines are widely used in offshore environments for the purpose of fluid or hydro-carbon products conveyance. When a pipeline is placed on an erodible seabed, a scour hole around the pipeline is likely to form because of the waves, currents or dynamic motions of the pipeline itself (forced-vibration or vortex-induced-vibration)^1,2. To improve the understanding of the scour mechanism around a subsea pipeline, measurements of the turbulent flow fields and estimations of the bed shear and normal stresses within the pipeline-fluid-seabed interaction region are essential in addition to measurements of the scour hole dimension^{1,2,3,4,5,6,7}. In an environment where the bed shear and normal stresses are extremely difficult to be determined because the flow field is unsteady and the bottom boundary is rough, measured instantaneous near-boundary stresses (at approximately 2 mm above the boundary) could be used as their surrogate^8,9. In the past few decades, scour around a vibrating pipeline has been studied and published without quantitatively presenting the values of the sophisticated flow fields around the pipeline within the scour hole^{3,4,5,10,11,12,13,14,15,16,17,18}. Therefore, the goal of this method paper is to provide a novel experimental protocol for visualizing the detailed flow fields and to determine the near-boundary shear and normal stresses within an equilibrium scour hole induced by a forced vibrating pipeline. It should be noted that the pipeline-fluid-seabed interaction process in this study is in a quiescent water environment rather than those with unidirectional currents and waves.

This experimental method consists of two important components, namely, (1) simulation of pipeline (forced) vibrations; and (2) measurements of the flow fields around the pipeline. In the first component, the vibrating pipeline was simulated in an experimental flume by using a vibrating system, which has a servo motor, two connecting springs, and pipeline supporting frames. Different vibration frequencies and amplitudes can be simulated by adjusting the motor speed and location of the connecting springs. In the second component, the time-resolved particle image velocimetry (PIV) and wavelet transform techniques were adopted to obtain high temporal and spatial resolution flow field data at different pipeline vibration phases. The time-resolved PIV system consists of a continuous wave laser, a high-speed camera, seeding particles, and cross-correlation algorithms. Although PIV techniques have been widely used in obtaining steady turbulent flow fields^{19,20,21,22,23,24,25}, applications in complex unsteady flow field conditions, such as cases of pipeline-fluids-seabed interaction, are relatively limited^8,9,26,27. The reason probably is because traditional single-time-interval cross-correlation algorithm of PIV techniques is unable to accurately capture the flow features in unsteady flow fields where a relatively high velocity gradient is present^9,20. The method described in this paper can solve this problem by using the multiple-time-interval cross-correlation algorithm^9,28.

Protocol

1. Laboratory safety check

1. Review the safety rules relating to the use of the laser and flume system.
2. Ensure that the safety training requirements of the laboratory have been met.
   NOTE: In this experiment, a set of 5W air-cooling continuous wave laser with a wavelength of 532 nm and a glass-sided straight flume (Figure 1) with dimensions of 11 m length, 0.6 m width, and 0.6 m depth are used. The basic safety recommendations for these two apparatuses are as follows:
   1. Check the potential reflection surfaces in the laser line-of-sight prior to testing; wear safety goggles when operating the laser device.
   2. Avoid having eyes at the level of the laser beam during the experiments and be careful of reflected laser lights when handling the optical elements or reflective tools.
   3. Ensure that the water hose does not fall off and that there is always no water overflowing from the flume.

2. Flume and seabed model setup

1. Prepare the erodible seabed model located in the middle of the flume.
   NOTE: The sediment material used in this study was a uniformly distributed medium sand with a median grain size d₅₀ = 0.45 mm, relative submerged particle density Δ = 1.65 and geometric standard deviation σ_g = 1.30.
2. Compact and level the seabed using a sand leveler.
3. Slowly fill the flume with a water hose and make sure that a flat seabed surface is intact during the filling process; stop filling when the water level has reached a depth of 0.4 m above the seabed.
4. Clear the flume top platform and glass for setting up the pipeline model and PIV system.

3. Pipeline model and vibration system setup

1. Use a prefabricated pipeline model in the form of an acrylic cylinder with a diameter of 35 mm and a length of 0.56 m.
2. Mount the pipeline model on an aluminum supporting frame, which, in turn, is connected by two springs to a moveable pole on another fixed frame that is locked on the top rails of the flume, as illustrated in Figure 2. Fix the supporting frame inside the fixed frame by using four bearings to ensure that the supporting frame could freely vibrate only along the vertical direction (Figure 2).
3. Use a connecting rod to tie the moveable pole to a servo motor mounted on the top of the fixed frame. In this study, the weight of the assembled vibrating system, including the pipeline model and the aluminum frames, is 1.445 kg, which has an equivalent mass ratio (m*) of 2.682; a natural frequency (f_N) of 0.82 Hz; and damping ratio (ζ) of 0.124.
4. Adjust the moveable pole and the supporting frame to obtain a certain gap ratio between the pipeline and seabed.  In this study, G/D =1, where G is the vertical distance between the bottom of the pipeline and initial seabed surface; and D is the pipeline diameter.
5. Turn on the servo motor to induce a forced vibration on the pipeline; adjust the supporting frames and four bearings to ensure that the pipeline vibration is along the vertical direction. Turn off the servo motor when adjustments of the supporting frames have completed.
6. Compact and level the seabed again before running the experiment if the seabed model is disturbed in 3.5.

4. PIV setup

1. Place the laser device on the top of the flume and install the laser sheet forming optics.
2. Turn on the laser device and adjust the laser sheet forming optics so that an illuminated flat sheet inside the field-of-interest is formed.
   NOTE: In this study, the illuminated green laser sheet is 1.5 mm thick, parallel to the flume glass walls and is cast downward into the water along the centerline of the flume. The field-of-interest of this study refers to the interaction region of the pipeline-fluid-seabed and is confined to the right half side of the pipeline. The shadow of the pipeline will be seen on the left half side of the pipeline.
3. Set up the high-speed camera.
   NOTE: For this study, a high-speed camera with 12-gigabyte-memory storage and a maximum resolution of 2.3 Mpx (1920×1200) is used (e.g., Phantom Miro LAB 320). The detailed operation procedures are as follows:
   1. Mount the lens with appropriate focal length on the high-speed camera. Screw the high-speed camera onto a height-adjustable tripod; adjust the camera to the level of the observation region with its axis perpendicular to the illuminated laser sheet.
      NOTE: This study uses a 60 mm prime lens at its maximum aperture of f/2.8.
   2. Connect the camera to the computer by using an Ethernet cable and turn on the camera control software (e.g., Phantom PCC 2.6); turn on the camera and connect it to the computer in the camera control software interface.
   3. Adjust the tripod to ensure that the field-of-view of the camera covers the pipeline-fluid-seabed interaction region; level the camera using the built-in bubble level on the tripod; tune the focus ring on the lens to ensure that the laser sheet is clear on the focal plane.

5. Experimental setup optimization and calibration

1. Add PIV seeding particles to the test section of the flume.
   NOTE: The seeding particles used in this study were aluminum powders with a diameter of 10 μm and a specific density of 2.7.
2. Enhance the light intensity of the laser sheet if necessary.
3. Verify the focus of the camera by observing the illuminated seeding particles on the laser sheet through a live camera view on the computer; fine-tune the focus ring, if necessary, to ensure that the seeding particles are sharp and in focus.
4. Place a calibration ruler inside the field-of-view on the plane of the laser sheet and capture one calibration image.
   NOTE: The adopted resolution of the image in this study was 1600 × 1200 pixels.
5. Select a proper sampling rate for data collection.
   NOTE: The chosen sampling rate should ensure that the seeding particle displacement within a pair of images is less than 50% of the maximum interrogation window length. In this study, the maximum interrogation window size is 32 × 32 pixels and the adopted sampling rate is 200 frames per second.
6. Turn off the laser and camera when steps 5.1-5.5 are completed.

6. Running the experiment and data collection

1. Place a transparent acrylic plate (20 mm thick) below the laser source and on the water surface, to suppress water surface fluctuations, and ensure tranquil optical access for the laser light.
2. Turn on the servo motor to induce forced vibrations on the pipeline model.
   NOTE: In this study, the induced frequency of the servo motor is f₀ = 0.3 Hz.
3. Keep the vibration system running for (t =)1440 min to obtain a quasi-equilibrium scour hole beneath the vibrating pipeline.
4. Turn on the laser and adjust the output power to the optimized intensity. Turn on the camera and camera control software and apply the calibrated settings to the camera. Turn off the background lights in the laboratory.
5. Start recording the seeding particle-laden flow field image with the sampling rate selected in 5.6 by clicking the Capture Bottom in the camera software control software.
   NOTE: For every single recording in this study, the camera storage allows 1,000 images to be captured.
6. Once the data collection is completed, review the recorded image quality and check if the seeding particle density per interrogation window (32×32 pixels) is greater than 8. Save the recorded file if satisfied, otherwise, the seeding density is increased by slowly injecting seeding solutions in the observation region, and repeat steps 6.3-6.5.
7. Repeat steps 6.3-6.5 to collect more data sets.
   NOTE: For this study, more than 20,000 images were taken to ensure that enough raw data are obtained for calculating the flow velocities, vorticities, turbulence, and near-boundary stresses.
8. Turn off the laser device, camera, and server motor when all the data collections are completed; turn on the background lights in the laboratory.

7. Data processing

1. Open the software; click the File folder button on the toolbar and load the calibration image taken in step 5.4.
   NOTE: Use the data processing program for pipeline displacement tracking and flow field calculations software (e.g., PISIOU).
2. Click the Scale setup button on the toolbar; measure a known distance on the calibration image to calculate the scale of the image.
   NOTE: The calculated image scale was 0.1694 mm/pixel.
3. Click the Origin button on the toolbar; set the origin of the coordinates on each image.
4. Extract the displacement time-series of the vibrating pipeline from the recorded images.
   1. Load the raw images taken in step 6. Then click on the Parameter panel, enter the number of data files and the sample rate.
   2. Apply the Low pass filter in the Image filter menu.
      NOTE: This operation will allow the edge of the pipeline (target to be tracked) to be readily recognized on the processed images (see Figure 3a).
   3. In the toolbar, click the PTV module. Then click Tracing point button, select the center point of the pipeline. Go to PTV tools, adjust Gamma, Light Gate and Median Filter to single out the pipeline outline in the image. Click the Object tracking button on the toolbar; select the target region (i.e. the pipeline) on the processed image and track the displacement of the vibrating pipeline from consecutive processed images; record the displacement time-series, η(t), of the vibrating pipeline for subsequent flow field data processes (see Figure 4).
   4. Export and save the pipeline displacement time-series data for further calculations.
5. Determine instantaneous velocity fields from the recorded images.
   1. Go to PTV tools, click Default button to resume the raw image for subsequent PIV analysis. Deactivate the PTV module by clicking the PTV module. Open the Parameter panel on the toolbar; specify the velocity vector calculation parameter.
      NOTE: In this study, a multi-pass iteration process is adopted as the interrogation windows, which started from 32 × 32 pixels, then passed with 16 × 16 pixels, and ended with 8 × 8 pixels; all passes use a 50% overlap between adjacent sub-windows.
   2. Apply the Laplacian filter function in the Image filter menu to the raw images to highlight the seeding particles and filter out undesirable scattering light (see Figure 3c).
   3. Click the Boundary button on the toolbar, set the geometric mask on the images to exclude the seabed region for further calculation. Click Boundary Save button to save the boundary data.
   4. Click the Run button on the toolbar to calculate the instantaneous velocity fields for different vibrating phases using the cross-correlation method.
      NOTE: In this study, a multi-time interval algorithm is adopted to reduce the bias error due to high velocity gradient in the flow field (see Figure 5). The adopted multiple-time intervals for cross-correlation calculations are Δt, 3Δt, 9Δt and 21Δt (Δt = 5 ms). The satisfactory correlation criterion is greater than 70%.
   5. Export and save the instantaneous velocity fields data for further analysis.
6. Determine the phase-averaged velocity fields from the calculated instantaneous velocity fields with the algorithm as described in Newland 1994^29,30  and Hsieh 2008²⁸.
   NOTE: The calculation procedures for this step is described as follows:
   1. Apply the wavelet transform function to the displacement time-series, η(t), of the vibrating pipeline to obtain the instantaneous phase for each instantaneous velocity field. The wavelet transform function is defined as:
                         (1)
      where W_Ψ is wavelet coefficient; α and β are the scale and translation parameters, respectively; the function Ψ is the Morlet function and is calculated as ; the superscript “*” denotes the complex conjugate. The instantaneous phases, Φ, of the vibrating pipeline that correspond to the different pipeline displacements can be calculated from:
                       (2)
   2. Average the instantaneous velocity fields with the same phase to obtain the phase-averaged velocity fields.
   3. Determine the flow vorticity, ω₂, in the calculated phase-averaged velocity fields from:
                  (3)
      where and are phase-averaged velocities along x and y directions.
7. Load the calculated phase-averaged velocity and vorticity data in Tecplot software for visualization.
8. Determine the near-boundary shear and normal stresses from the calculated instantaneous velocity fields with the algorithm as described in Hsieh et al. 2016 ⁹. The calculation procedures for this step is described as follows:
   1. Extract the near-boundary velocity data (0-5 mm above the seabed) from the calculated phase-averaged velocity flow fields.
   2. Compute the near-boundary shear stresses, t_s and normal stresses, t_n, along the scour profile (approximately 2 mm above the scour hole boundary) for different phases within one vibrating cycle. Note: The calculation equations are as follows:
      ,               (4)
      where, μ = dynamic viscosity of the fluid (herein taken as 1×10^-3 Pa∙s); u_p = near-boundary velocity parallel to the bed; u_n = near-boundary velocity perpendicular to the bed; _n = normal distance from the bed.
9. Load the calculated near-boundary shear and normal stresses data in a software (e.g., Tecplot) for visualization.

Results

An example of the comparison between the raw image and processed image of the pipeline displacements tracking and instantaneous velocity calculation is shown in Figure 3. As shown in Figure 3b, the seeding particles and noise in the raw image are filtered out and the shining pipeline edge is retained to obtain the displacement time series. As shown in Figures 3c, light scatters/reflections around the seeding particles, pipeline ...

Discussion

The protocol presented in this paper describes a method for visualization of the two-dimensional flow fields and determination of the near-boundary flow stress fields around a forced vibrating pipeline in an equilibrium scour hole by using the PIV techniques. Since the designed pipeline motion is one-dimensional along the y direction, preparing and adjusting the pipeline model and vibration system to fulfill this objective are critical prerequisites for a successful outcome. Any undesirable motions of the pipeli...

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Camera control software	Vision Research	Phantom PCC 2.6	Camera control, image data acquisition and processing
Camera lens	Nikon Chiyoda		Nikor  60mm, f=2.8 prime lens
Continuous wave laser	Beijing Laserwave optoelectronics technology co. ltd.		PIV Laser source; Nd:YAG laser, 532 nm; air-cooling
High-speed camera	Vision Research	Phantom Miro LAB 320	Image data recording
Laser sheet forming optics	Thorlabs Inc		Transform the point laser to a thin laser sheet
Pipeline model	ZONCEPZ SOLUTIONS		Acrylic cylinder with a diameter of 35 mm
Pipeline vibration system	ZONCEPZ SOLUTIONS		Consists of a sever motor, two connecting springs and pipeline supporting frames.
PIV calcuation software	AXESEA Engineering Technology Limited Co.	PISIOU	Image data processing for obtaining flow fields and pipeline displacements
PIV seeding materials	Shimakyu		Aluminum powder with a diameter of 10um
Recirculating flume	SZU ENGINEERING PTE LTD		Glass-sided, 11 m long, 0.6 m wide, and 0.6 m deep
Tri-pod	MANFROTTO	SKU MT190GOC4US 410	Camara supporting