Relief Valve Modeling (English Units)

Relief Valve Modeling (Metric Units)

Summary

A system pumping crude oil over a long distance experiences a sudden valve closure. The pump is not able to stop before the valve is completely closed. A relief valve will be considered for surge mitigation.

Topics Covered

  • Modeling relief valve profiles during pressure surges

  • Examining the effect of valve rate limits on a system

  • Observing differences between relief valve profiles

Required Knowledge

This example assumes the user has already worked through the Beginner: Valve Closure example, or has a level of knowledge consistent with that topic. You can also watch the AFT Impulse Quick Start Video (English Units) on the AFT website, as it covers the majority of the topics discussed in the Valve Closure example.

Model Files

This example uses the following files, which are 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 Fluid Properties Group

  1. Open Analysis Setup from the toolbar or from the Analysis menu.

  2. Open the Fluid panel then define the fluid:

    1. Fluid Library = User Specified Fluid

    2. Name = Crude Oil

    3. Density = 0.85 S.G. water

    4. Dynamic Viscosity = 4 centipose

    5. Bulk Modulus = 200000 psia

    6. Vapor Pressure = 10 psia

Step 3. 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.

We will first evaluate the system without a relief valve to see if the maximum pressure of 350 psi is exceeded.

Figure 1: Workspace with no relief valve added

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 Model tab

    1. Pipe Material = Steel - ANSI

    2. Size = 24 inch

    3. Type = STD (schedule 20)

    4. Friction Model Data Set = Standard

    5. Length = Use table below

Pipe Length (feet)
1 20
2 10000
3 20
4 20

Junction Properties

To define the pump, a pre-built library will be used. We will need to first connect the library as follows

  1. From the Library menu select Library Manager.

  2. At the bottom, click Add Existing Library.

  3. Browse to the Relief Valve Modeling.dat library file in the Impulse Examples folder and click OK. The library should now appear in the Library Browser under the name Relief Valve Modeling. Make sure the box is checked next to this library.

  4. Click OK to exit the Library Manager.

Now proceed with defining the junctions.

  1. J1 Reservoir

    1. Tank Model = Infinite Reservoir

    2. Liquid Surface Elevation = 10 feet

    3. Liquid Surface Pressure = 10 atm

    4. Pipe Depth = 10 feet

  2. J2 Pump

    1. Library Jct = US - Crude Oil Pump

    2. Inlet Elevation = 0 feet

Note: The Model Combined Check Valve and Pump Option is used in this example to simplify the model. For higher accuracy it is generally recommended to model a check valve at the pump discharge using a Check Valve junction instead. See the Check Valve Modeling example for more information.

Figure 2: Pump Model tab with pump information from library

Figure 3: Pump Properties Transient Data

  1. J3 Branch

    1. Elevation = 0 feet

  2. J4 Valve

    1. Inlet Elevation = 0 feet

    2. Loss Model tab

      1. Valve Data Source = User Specified

      2. Loss Model = Cv

      3. Cv = 10000

    3. Transient tab

      1. Transient Special Condition = None

      2. Initiation of Transient = Time

      3. Transient Data = Absolute Values

Time (seconds) Cv
0 10000
1 2000
2 0
30 0

Figure 4: Valve Properties Transient Data

  1. J5 Assigned Pressure

    1. Elevation = 0 feet

    2. Pressure = 12 atm

    3. Pressure Specification = Stagnation

Ø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 4. Define the Pipe Sectioning and Output Group

ØOpen Analysis Setup and open the Sectioning panel. When the panel is first opened it will automatically search for the best option for one to five sections in the controlling pipe. The results will be displayed in the table at the top. Select the row to use one section in the controlling pipe. 

Step 5. Define the Transient Control Group

ØOpen the Simulation Mode/Duration panel in the Transient Control group. Enter 30 seconds for the Stop Time.

All groups should now be complete and the model is ready to run. If all groups in Analysis Setup have a green checkmark then Click OK and proceed. Otherwise, enter the missing information.

Step 6. 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 Output button at the bottom of the Solution Progress window.

Step 7. Examine the Output

The maximum pressure requirement for the line is 350 psia. A profile graph should be created to review the maximum pressure along the line as follows:

  1. Go to the Workspace and select all of the pipes

  2. Right-click the Workspace and select Generate Graph for Pipes > Pressure Static

  3. Click Add Graph to List and name it Static Pressure Profile

The maximum/minimum pressure profile graph can be seen in Figure 5. The pressure goes above the maximum pressure requirement for the majority of the line. To decrease the maximum pressure we will add a relief valve to the line near to the closing valve, which is the source of the pressure surge.

Figure 5: Maximum/minimum pressure profile along the pipeline

Step 8. Create Child Scenarios

We will create two scenarios to compare the model with and without a relief valve.

ØIn the Scenario Manager right-click the Base Scenario and click Create Child. Name the scenario No Relief Valve. Repeat this process to create a second scenario, Passive Relief Valve.

Now we will add the relief valve to the Passive Relief Valve scenario. Make sure the relief valve scenario is active in the Workspace by double-clicking the scenario name, then add a relief valve line as shown in Figure 6.

Figure 6: Workspace with relief valve and vessel added to the line

The crude oil will need to be contained in a closed vessel at the exit of the relief valve. For this example we will assume that the transient is fast, thus the pressure in the closed vessel will not change during the transient. The surge tank junction will be used to represent the closed, pressurized vessel. If the pressure in the closed vessel were to change significantly during the transient, then the surge tank would not be a good representation of this.

Define the additional junctions and pipes as follows.

Pipe Properties

  1. Pipe Model tab

    1. Pipe Material = Steel - ANSI

    2. Size = 12 inch

    3. Type = STD

    4. Length = 20 feet

Junction Properties

  1. J6 Relief Valve

    1. Inlet Elevation = 0 feet

    2. Valve Model tab

      1. Valve Model = Internal

      2. Profiles = Passive

      3. Pressure Type = Pressure

      4. Constant Backpressure (Balanced) = Checked

      5. Overpressure = 325 psia

      6. Set Pressure = 300 psia

      7. Blowdown Pressure = 250 psia

    3. Loss Model tab

      1. Loss Model = Cv (Variable)

      2. Variable Data = Linear based on setpoints

      3. Full Open Cv = 500

Note: The Ignore Relief Valve option on the Optional tab can be used to keep the relief valve closed. This is useful for running cases with no relief valve in the system, or where the relief valve fails to open without deleting the relief valve junction.

Figure 7: Relief Valve Model Properties

  1. J7 Surge Tank

    1. Cross-Sectional Area = Constant

    2. Value = 1000 feet2

    3. Initial State of Liquid Height = Liquid Height Known (Possible Inflow/Outflow)

    4. Surface Pressure = 11 atm

    5. Initial Liquid Height = 0 feet

    6. Junction Elevation = 0 feet

Figure 8: Surge Tank Properties

Step 9. Define the Pipe Sectioning and Output Group

New pipes were added to the model, so the model needs to be sectioned again.

ØOpen Analysis Setup and open the Sectioning panel. When the panel is first opened it will automatically search for the best option for one to five sections in the controlling pipe. The results will be displayed in the table at the top. Select the row to use one section in the controlling pipe.

Step 10. Run the Model

Click Run Model from the toolbar or from the Analysis menu. When complete, click the Graph Results button at the bottom of the Solution Progress window.

Step 11. Examine the Output

Reload the Static Pressure Profile graph item by double-clicking the graph item in the Graph List Manager. The pressure profile should now appear as shown in Figure 9 below.

With the relief valve the maximum pressure reached in the pipeline is about 319 psia, which is below the required maximum pressure of 350 psia.

Now create a Transient Jct graph for the relief valve, and select Cv as the parameter to be graphed (Figure 10). Save this graph as Relief Valve Cv in the Graph List Manager.

Taking a closer look at the graph it can be seen that the valve opens over about 0.2 seconds, and closes over about 1.7 seconds(time from maximum Cv to zero). This is by design - the relief valve should open quickly to reduce pressure, but close more gradually to avoid creating a pressure wave. When defining the relief valve for this model we did not take into account the actual closing time that will occur at the relief valve. We will add rate limits to the model to better account for the closing time at the valve.

Figure 9: Maximum/minimum pressure profile along the pipeline with a passive relief valve

Figure 10: Relief valve Cv vs time

Step 12. Add Rate Limits to the Relief Valve

A relief valve's rate limit represents the rate at which the valve can open and close. If no rate limit is defined for the relief valve, it is possible for the valve to open/close nearly instantaneously. An instant opening for a relief valve is idealized, as the high pressures can be instantly alleviated in the pipeline once the valve setpoint is reached. Instant closure of the relief valve would represent the worst case scenario, as the rapid closure would produce the maximum pressure rise at the relief valve. In reality it is not physically possible for the relief valve to close in one time step. In addition, relief valve manufacturer's often design the valves to close slowly to avoid surge events that could occur due to a sudden closure. Adding rate limits allows for a more realistic representation of how fast the valve position can change over time.

Add rate limits to the relief valve as follows:

  1. In the Workspace window, under Scenario Manager, create a new child scenario under the Passive Relief Valve scenario named With Rate Limits, and load it as the active scenario.

  2. Open the properties window for Relief Valve J6.

  3. On the Loss Model tab check the boxes next to both Apply Closing Limit and Apply Opening Limit (Figure 11)

  4. Add a Max Closing Cv change of 150 and a Max Opening Cv change of 1500

  5. Click OK

  6. Run the Model, click on Graph Results

  7. Open your saved graph Relief Valve Cv under Graph List Manager

  8. In the Graph Parameters section, click the Multi-scenario button, then select the With Rate Limits and Passive Relief Valve scenarios.

  9. Click Generate to create the graph as shown in Figure 12.

  10. Repeat steps 7 and 8 above for the Static Pressure Profile graph, as shown in Figure 13.

In Figure 12 it can be seen that though the relief valve still opens in about 0.2 seconds, the relief valve now closes at a slower rate of about 3.1 seconds.

From Figure 13 it can be seen that the slower opening speed results in a higher maximum pressure in the pipeline since the relief valve takes longer to fully open. The longer closing time practically eliminates surge waves caused by the closure of the relief valve, which causes a slightly lower minimum pressure along the pipeline.

Figure 11: Relief Valve Loss Model tab with opening and closing rate limits defined

Figure 12: Relief valve Cv profile with and without rate limits

Figure 13: Comparison of maximum and minimum pressures along the pipeline when the passive relief valve is defined with and without rate limits

Step 13. Rupture Disk

We are now going to compare the passive valve above to a rupture disk relief valve. A rupture disk will burst open when the set pressure is met, and will remain open after that. The rupture disk will still relieve to a closed tank to keep the crude oil contained.

Update the model as follows:

  1. In the Workspace window, under Scenario Manager, right click on the Passive Relief Valve scenario and choose Clone Without Children. Name the new scenario Rupture Disk.

  2. Double click on the Relief Valve J6.

  3. In the Profiles list select Rupture Disk (Figure 14)

  4. Set Pressure (dP) = 150 psid

Note: The rupture disk will open on a pressure difference, unlike the passive relief valve which was defined using a fixed setpoint at which the relief valve would open. This is because the passive relief valve modeled previously had a mechanism to maintain a constant back-pressure on the valve, which is not used with the rupture disk. The pressure difference setpoint for the rupture disk was found by subtracting the outlet tank pressure from the desired opening pressure.

  1. In the Loss Model tab, Cv (Constant) = 300

  2. Click OK to close the window.

Figure 14: Inline Exit Rupture Disk Properties

The scenario is now complete and ready to be run.

ØRun the Model and select Graph Results.

Load the Relief Valve Cv graph from the Graph List Manager, as is shown in Figure 15. Unlike the passive relief valve type, the rupture disk opens instantly, and remains open, as can be seen in the graph.

Figure 15: Relief valve Cv over time for the rupture disk relief valve

Next we want to compare the pressure results for each of the relief valve types.

  1. Reload the Static Pressure Profile graph.

  2. In the Graph Parameters area in the Quick Access Panel click the Multi-scenario dropdown and choose Select Scenarios

  3. Select both the Rupture Disk and With Rate Limits scenarios, then click OK

  4. Click Generate to incorporate the additional scenario (Figure 16)

The rupture disk generally maintains a lower maximum pressure in the pipeline, and overall reaches a lower minimum pressure. This makes sense considering that the rupture disk is allowed to open instantaneously (over one time step), thus the pressure in the pipeline is lowered faster for that scenario. In reality the rupture disk will not open instantaneously, but will still likely open faster than a passive relief valve type.

A lower minimum pressure occurs with the rupture disk because the rupture disk never closes.

Figure 16: Comparison of maximum and minimum pressures with the passive and rupture disk relief valve types

Summary

Relief Valves provide pressure relief to systems during pressure surges. Impulse can model relief valves as an internal valve model, exit valve model, or inline exit valve model. An internal valve model connects two pipes acting as a gate providing relief from elsewhere in the system when necessary. An exit valve essentially models a spray nozzle at the end of a pipeline. An inline exit valve provides relief to the atmosphere along the length of some piping. These valve models can be profiled in a number of ways. In this example we worked with graphing results with a passive profile and rupture disk profile.