PD Pulsation Study - PFA (Metric Units)

PD Pulsation Study - PFA (English Units)

Summary

A system with positive displacement pumps has a known pulsation issue in a section of piping with a closed valve. This issue is caused by the system’s natural frequency matching up with certain PD pump frequencies, resulting in severe pressure oscillations at the closed valve. Use AFT Impulse and the PFA (Pulsation Frequency Analysis) module to discover the problematic pump frequencies and what the pressure response at the closed valve is at these specific pump frequencies.

Note: This example can only be run if you have a license of the PFA module

Topics Covered

  • Ringing systems to find natural acoustic frequencies

  • Finding system excitation frequencies

  • Determining 'worst case' PD pump speeds at excitation frequencies

  • Calculating pressure response of system at above PD speeds

Required Knowledge

This example assumes the user has already worked through the Walk-Through Examples section, and has a level of knowledge consistent with the topics covered there. If this is not the case, please review the Walk-Through Examples, beginning with the Beginner: Valve Closure example. You can also watch the AFT Impulse Quick Start Video (Metric Units) on the AFT website, as it covers the majority of the topics discussed in the Valve Closure example.

This example is based on the paper by James M. Blanding and Trey Walters.

Model File

This example uses the following file, which is installed in the Examples folder as part of the AFT Impulse installation:

Step 1. Start AFT Impulse

From the Start Menu choose the AFT Impulse 9 folder and select AFT Impulse 9.

To ensure that your results are the same as those presented in this documentation, this example should be run using all default AFT Impulse settings, unless you are specifically instructed to do otherwise.

Step 2. Define the Modules Panel

First the PFA add-on modules needs to be activated for use in the model. To do this, open Analysis Setup from the toolbar or Analysis menu. On the Modules panel, check the box next to Activate PFA. After the PFA module is activated, the Use option should being selected by default.

Step 3. Define the Fluid Properties Group

Open the Fluid panel then define the fluid:

  1. Fluid Library = User Specified Fluid

  2. Name = PFA Fluid

  3. Density = 0.92 gram/cm3

  4. Dynamic Viscosity = 0.95 centipoise

  5. Bulk Modulus = 9653 bar

Step 4. Define the Pipes and Junctions Group

At this point, the first two groups are completed in Analysis Setup. The next undefined group is the Pipes and Junctions group. To define this group, the model needs to be assembled with all pipes and junctions fully defined. Click OK to save and exit Analysis Setup then assemble the model as shown in the figure below.

Figure 1: Layout of PD Plunger Pulsation Study example

The system is in place but now we need to enter the input data for the pipes and junctions. Double-click each pipe and junction and enter the following data in the properties window.

Pipe Properties

  1. Pipe Material = (User Specified)

  2. Inner Diameter = Use table below

  3. Length = Use table below

  4. Wavespeed = User Specified Wavespeed

    1. Wavespeed = 915 meters/sec

  5. Friction Model Data Set = User Specified

  6. Selection = Absolute Roughness

  7. Value = 0.00152 cm

Pipe Inner Diameter (cm) Length (meters)
1 30.5 1
2 30.5 1
3 5.5 12
4 5.5 17
5 8.3 10
6 4.3 6
7 8.3 15
8 4.3 1
9 8.3 9

Junction Properties

Note: For the sake of simplicity, assigned flow junctions J1 and J2 are being used to represent the PD pumps. The Pulsation Setup window will be used later to place the necessary transient data in the assigned flow junction for the pump being analyzed.

  1. Assigned Flows J1 & J2

    1. Elevation = 0 meters

    2. Volumetric Flow Rate = 8 m3/hr

  2. Area Changes J3 & J4

    1. Inlet Elevation = 0 meters

    2. Type = Abrupt Transition (Cylindrical)

  3. Tee/Wye J5 & J6

    1. Elevation = 0 meters

    2. Loss Model = Simple (no loss)

  4. Dead End J7

    1. Elevation = 0 meters

  5. Area Changes J8 & J9

    1. Inlet Elevation = 0 meters

    2. Type = Conical Transition

    3. Angle (μ) = 90 degrees

  6. Spray Nozzle J10

    1. Elevation = 0 meters

    2. Loss Model = Cd Spray (Discharge Coefficient)

    3. Geometry = Spray Nozzle

    4. Exit Properties = Pressure

    5. Exit Pressure = 0 barG

    6. Cd (Discharge Coefficient) = 1

    7. Discharge Flow Area = 0.148 cm2

ØTurn on the Show Object Status from the View menu to verify if all data is entered. If so, the Pipes and Junctions group in Analysis Setup will have a check mark. If not, the uncompleted pipes or junctions will have their number shown in red. If this happens, go back to the uncompleted pipes or junctions and enter the missing data.

Step 5. Define the Pulsation Setup Group

Defining the initial pulse or ring is an important step in determining the natural frequencies that different stations in the model will respond to. This is the numerical equivalent of hitting the system with a hammer to see how it responds at different frequencies. This first step does not represent a realistic operating condition, but rather reveals what frequencies of real world operation may present operational problems.

Open Analysis Setup to see the Pulsation Setup group below the Transient Control group. Use the Pulsation Setup group to specify input regarding the initial pulse and other information about the pulsation source. There are two panels in the Pulsation Setup group: the Pulse Setup and PD Pump Settings panels. The PD Pump Settings panel only requires input if the source of pulsation is a positive displacement pump. Because the source of pulsation in our model is a positive displacement pump, we will need to complete the input on both of these tabs. Enter the data below to define the Pulse Setup panel.

  1. Applied At = J2 (Assigned Flow)

  2. Start Time = 0 seconds

  3. Magnitude = Automatic (Twice Steady-State Flow)

  4. Peak = 16.00 m3/hr

  5. Cutoff Frequency = PD Pump Characteristics

  6. Filter Type = Chebyshev

  7. Order = 10

  8. Ripples Allowed = 1 dB

  9. Pulsation Source = Positive Displacement (PD) Pump

  10. Number of Heads = 3

  11. Pump Speed = 250 RPM

  12. Speed Buffer = 15 Percent (%)

Figure 2: Pulsation Setup window with Input

Now that all required input has been entered in the Pulse Setup panel, enter the following information on the PD Pump Settings panel (see Figure 3).

  1. Minimum Pump Speed = 125 RPM

  2. Bore Diameter = 7.6 cm

  3. Stroke = 5.2 cm

  4. Rod Length = 25.4 cm

  5. Clearance = 70 Percent (%)

  6. Suction Pressure = 4.5 bar

  7. Discharge Pressure = 415 bar

Figure 3: PD Pump Settings tab with Input

Step 6. Define the Pipe Sectioning and Output Group

Open the Sectioning panel in Analysis Setup. While using the AFT Impulse PFA module, the pulsation setup needs to be completed prior to sectioning the pipes in the model. To section the model’s pipes for the simulation, complete the following:

  1. Select the significant pipes that should consider wavespeed variance. In this example we will consider the variance in all modeled pipes, except for Pipe P8. You can see that this pipe has been excluded in Figure 4. Pipe P8 is excluded because it represents a flow meter and not an actual pipe. Units such as flow meters and dampeners may be better modeled with a pipe than with a junction if they have significant length and associated acoustic interaction aspects in the physical system. Ensure that the box next to each pipe is selected under Consider Variances, except for Pipe P8.

  2. Click Update Results and select one section in the controlling pipe by clicking on one section in the table under Search Results. This window is displayed in Figure 4.

  3. Select Show Pulsation Graphs located at the bottom left of the Sectioning panel. A graph of the pulse will appear (see Figure 5). You can also view graphs for the FFT of the Pulse, the LPF of the Pulse FFT, and the FFT of the LPF in the Graph Parameter graph that appears. The LPF of Pulse is shown in Figure 6. After you are done viewing these graphs, click Close to return to the Section Pipes window.

  4. Click OK to save and close Analysis Setup.

  5. Click OK for the pop-up window notifying you that the forcing function has been placed on J2 (the junction that was selected in the Pulsation Setup window) as transient data.

  6. Open J2 to verify that the numerical ring has been applied to the transient (see Figure 7 for the transient input). This transient is a brief jump in flow rate at junction J2.

Figure 4: Section Pipes window

Figure 5: Graph of Pulse in this example

Figure 6: Graph of LPF of Pulse

Figure 7: Assigned Flow Properties with Transient from Pulsation Setup. The Show Graph button under the Transient Data table can be selected to view the graph of this data

Step 7. Run the Model

Click Run Model from the toolbar or from the Analysis menu. This will open the Solution Progress window. This window allows you to watch the progress of the Steady-State and Transient Solvers. When complete, click the Graph Results button at the bottom of the Solution Progress window.

Step 8. Examine the Results

  1. Find excitation frequencies to study

    1. On the Graph Control tab in the Quick Access Panel, click the Frequency tab to begin generating an Excitation Frequency Analysis graph

    2. We need to determine which frequencies excite the pump speeds within the pump speed range we have specified. For the purposes of this example, we will graph the Excitation Frequency graph for 10 different pipe stations. Note that it is the engineer's responsibility to evaluate all pipe stations that could be excited by various frequencies and to perform the necessary analysis.

    3. On the Graph Results tab, select the Frequency tab and select the inlet and outlet stations for the following pipes: 1, 2, 3, 5, and 6 (see Figure 8).

Figure 8: Graph Parameters selected for Excitation Frequency Analysis graph

    1. Click Generate. The resulting graph (shown in Figure 9) represents the frequencies that are excited by the pulse defined in the Pulsation Setup window.

Figure 9: Excitation Frequency Analysis graph

    1. From the graph in Figure 9, we see that there are several frequencies that produce a large pressure magnitude. Points on the plot with a large response magnitude are easily identified by finding local maxima on the graph. These maxima are the acoustic natural frequencies that are potentially the most damaging when excited by the identified pump speeds.

    2. Note that the graph in Figure 9 can quickly be changed to display the magnitude on a logarithmic scale in order to allow for easier viewing of the magnitude behavior changes with a changing frequency. To view the Excitation Frequency Analysis graph with the magnitude displayed on a logarithmic scale, generate the frequency graph as discussed above. Right-click on either axis and select the box next to Logarithmic. See Figure 10 for an image of the graph in Figure 9 with the amplitude changed to be displayed on a logarithmic scale.

Figure 10: Excitation Frequency Analysis Graph with Magnitude on Logarithmic Scale

  1. Evaluate excitation frequencies

    1. The identified frequencies can now be evaluated to determine pump speeds that would cause them to be excited, making them potentially problematic. Follow the steps below to select and evaluate the excitation frequencies.

    2. Left-click near the frequency of about 33.24 Hz and drag the mouse over the peak on the graph near this frequency. This will highlight and flag the local maximum. Figure 11 shows the frequency graph after the dialog flag appears.

Figure 11: Excitation Frequency Analysis with Dialog Flag

    1. Right-click on the dialog flag that appears on the peak of this frequency as shown in Figure 12.

Figure 12: Right Click on the Dialog flag to evaluate the frequency

    1. Select Evaluate Excitation Frequency to add this frequency to the Pump RPM Evaluation panel, which is labeled with Excitation Frequencies and Pump Speeds (RPM) (see Figure 13)

Figure 13: Pump RPM Evaluation panel

    1. Repeat steps a through c for the local maxima at the frequencies of 28.98 Hz and 6.89 Hz. Figure 14 shows the Pump RPM Evaluation Panel with the excitation pump speeds displayed.

Figure 14: Pump RPM Evaluation Panel with Excitation Frequencies and their Harmonic Multiples Displayed

    1. When choosing low or high frequencies on the graph, it is possible to select a frequency which is excited in the system but does not have a pump speed within the minimum and maximum pump speeds for the analysis. The minimum pump speed is defined directly in the PD Pump Settings in the Pulsation Setup window, while the maximum is calculated using the Pump Speed and Speed Buffer defined under Pulsation Source in Pulsation Setup. For example, try selecting the largest maxima, which occurs at 6.1 Hz. The message in Figure 15 appears.

Figure 15: Message Indicating No Excitation Pump Speeds Between the Minimum and Maximum Pump Speeds Exist

  1. Determine the pressure response of the system at excited frequencies

    1. In the Pump RPM Evaluation panel, corresponding pump speeds can be selected for study. At these speeds, the pressure response is important to characterize. By default, all speeds in the Pump RPM Evaluation panel are pre-selected and must be clicked to be turned off. Scenarios can then be created that will use the information from the PD Pump Settings to create a child Pump Speed scenario that will have the Flow vs. Time profile at the selected speed (RPM) that will excite the evaluated frequencies.

    2. We will create child scenarios for all six pump speeds in the Pump RPM Evaluation panel. To do this, ensure that all pump speeds are selected (not italicized).

    3. Click Create Scenarios Based On Selected Pump Speeds…

    4. A dialog box will appear stating that six scenarios will be created based on the selected speeds in the Pump RPM Evaluation panel. Click Yes to continue creating the child scenarios.

    5. Load the 133 RPM (33.2 Hz) scenario.

    6. The scenario will now have different transient information in J2, modeling the pump performance at 133 RPM. Go to the Workspace and open the Properties window for J2. On the Transient tab click Show Graph under the Transient Data table to verify the Pump Flow Forcing Function has been applied, as shown in Figure 16.

Figure 16: New Assigned Flow transient representing the PD Pump operating at 133 RPM

Run the 133 RPM (33.2 Hz) scenario to see the pressure response. Often, the most noticeable pressure response will occur in stagnant branched lines. For the purposes of this example, we will evaluate the pressure responses through pipe P6. However, engineers must ensure that all pipes that could experience significant pressure responses are analyzed.

  1. To create a pressure profile graph in pipe P6, go to the Graph Results tab and select Profile for the graph type

  2. Check the box for Pipe P6

  3. Select Pressure Stagnation for the Parameter

  4. Use bar for the units

  5. Click Generate

Figure 17 shows the profile plot for Pipe P6 at a pump speed of 133 RPM .

Figure 17: Max/Min Pressure Profile of Pipe P6 for 133 RPM

The pressure vs. time can be displayed for any pipe station. Plot the pressure vs. time for pipe P6 by following the directions below:

  1. In the Quick Access Panel, select Transient Pipe for the graph type

  2. Expand the pipe station menu for Pipe P6

  3. Add the outlet and inlet of Pipe P6 to the Graph These Pipes/Stations list

  4. Select Pressure Stagnation for the Parameter

  5. Use bar for the units

  6. Set the Time Frame to User Specified with a Start Time of 0 and a Stop Time of 2 seconds.

  7. Click Generate (see Figure 18).

Figure 18: Pressure Oscillations at Inlet and Outlet of Pipe P6 for 133 RPM

After looking at these graphs, go to the Output tab, which should open to the Warnings tab. The PFA module will check the Peak-to-Peak pressures and compare them to the permissible values in the API-674 standard. If the limits from the standard are violated in the model, a warning will appear. In this case, the the Peak-to-Peak pressures exceed the allowable limits. See Figure 19. Additionally, the vapor pressure margin is checked. This is the margin between the minimum pressure and 10% over the vapor pressure, per the API-674 standard. For a full summary, click on the Pulsation Summary tab, as is shown in Figure 20.

Figure 19: Warnings section the Output window with information about API-674 being exceeded

Figure 20: Pulsation Summary in the Output window

For the purposes of this example, we will repeat these steps for the pump speed of 221.6 RPM because this pump speed causes the largest pressure oscillations. We will show these graphs here. Though not shown, multi-scenario graphing was used to determine the pump speed that caused the maximum oscillations. When performing an analysis on your system, all pump speeds that excite the system should be evaluated.

Figure 21 shows the Max/Min Pressure Profile for P6 for the pump speed of 221.6 RPM .

Figure 21: Max/Min Pressure Profile of Pipe P6 for 221.6 RPM

Notice the large difference between the maximum and minimum pressures at the end of the pipe, located next to the Dead End junction.

Figure 22 shows the Max/Min Pressure Profile for P6 for the pump speed of 221.6 RPM.

Figure 22: Pressure Oscillations at Inlet and Outlet of Pipe P6 for 221.6 RPM

The large pressure oscillation in Pipe P6 can be seen in Figure 22. It is important to note that the inlet of Pipe P6 does not exhibit as large of a response as the outlet does. This oscillation, if allowed to continue, can lead to mechanical and support fatigue and thus failure.

Return to the Output tab of the Graph results window for the scenario with the pump speed of 221.6 RPM, and make note of the message in red at the bottom right of the Graph Results screen that states that Warnings Exist. This scenario also violates API-674 pulsation and pressure limits similar to the 133 RPM scenario.

PFA Analysis Summary

AFT Impulse and the PFA (Pulsation Frequency Analysis) module were used to discover the problematic pump frequencies occurring in a system with a known pulsation issue. The software was used to identify these frequencies, and after the frequencies were identified, a corresponding pump operating speed was examined to find the pressure response of the system. Using graphing features, the pressure response was plotted and studied. The steps followed lead to the successful pulsation analysis of the system.