<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ENG</journal-id><journal-title-group><journal-title>Engineering</journal-title></journal-title-group><issn pub-type="epub">1947-3931</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/eng.2015.76028</article-id><article-id pub-id-type="publisher-id">ENG-57298</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Pressure Surge Dependence on Valve Operations in a Pipeline Loading System
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>arinaadaa</surname><given-names>Thaddeus Lebele-Alawa</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Felix</surname><given-names>Ezekiel Oparadike</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Mechanical Engineering, Rivers State University of Science and Technology, Port Harcourt, 
Nigeria</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>lebele-alawa.thaddeus@ust.edu.ng(ATL)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>17</day><month>06</month><year>2015</year></pub-date><volume>07</volume><issue>06</issue><fpage>322</fpage><lpage>330</lpage><history><date date-type="received"><day>19</day>	<month>May</month>	<year>2015</year></date><date date-type="rev-recd"><day>accepted</day>	<month>20</month>	<year>June</year>	</date><date date-type="accepted"><day>23</day>	<month>June</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  This paper discusses the influence of valve operations on pressure surge in a pipeline. The valve is a protective type which remains open in a pipeline loading system during normal operation but shut down the system when there is an emergency such as storm. The data for the study were obtained from measurements at Agbada 1 flow station as well as log sheets. Also, calculations were made using existing and derived formulas to obtain the values of Crude Oil and Pipe parameters that could not be measured directly or derived from data or log sheets. Surge analysis was carried out on the pipeline system to ascertain changes in pressure and flow rates along the pipeline following valve shut down at any time using developed pressure and flow equations. The results of the simulation analysis showed remarkable changes in the fluid pressure and flow rates along the pipe on shut down at any time. The pipeline recorded the highest pressure of 37.4 bar against initial pressure of 25 bar at length 6000 m in 1.5 second valve closure. There is also remarkable pressure drop along the pipe capable of reducing the crude oil pressure below its vapour pressure. The flow is turbulent even before valve operation with Reynolds number as high as 57024.53. The model equations compute changes in pressure and flow rates at different points in a pipeline installed with emergency-relief coupling valve. This enables point of extreme and low pressure to be detected accurately in a pipeline which guides the engineer while positioning surge suppression devices which cushion the effects of pressure surge in any pipeline.
 
</p></abstract><kwd-group><kwd>Pressure Surge</kwd><kwd> Valve</kwd><kwd> Pipeline</kwd><kwd> Crude Oil</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>This paper addresses the issue of pressure surge arising from transmission of fluid (crude oil) in pipe following the instantaneous closure of an emergency relief coupling valve (ERV) installed in a pipeline system. The pressure changes at different points on the pipe arising from the shutdown of the ERV when transmitting a given quantity of crude oil in a given pipe diameter over given pipe length is to be determined with the associated closure time t given. This valve serves as a protection valve and it is always open during normal operations of the system. It is installed to close when there is an event such as storm that can damage pipe machines which leads to outbreak of fire or crude oil leakage. When this valve closes instantaneously, that is at t = 0, the initial head rises to infinity. However, in practice, it is not possible to close this valve instantaneously as it always takes some time.</p><p>Usually, most pipelines have installed surge suppression devices. In this case, surge suppression, to an extent is taken care of. What if these suppression devices fail or there is an installation failure? It means the whole system would have problem if the ERV eventually shut down. But this shut down to safeguard the system, leads to a fluctuating pressure in the system that threaten the system make up. This pressure changes is computed using developed model. Also of interest is the hydraulic force developed in the bend as a result of the ERV closure at time t.</p><p>Stone [<xref ref-type="bibr" rid="scirp.57298-ref1">1</xref>] carried out the analysis of avoiding pressure surge in pipeline system. According to the research, there is one prime requirement to avoid pressure surge damage to piping and pipeline. This has to do with bearing in mind that there is likely to be surge in a pipeline system either by measurement or by engineering analysis. The research concluded that damages arising from surge events are; catastrophic failure of the pipeline system or equipment and fatigue failure of the pipeline supports, instrumentation, equipment and components. Hence, pipe work can be designed to withstand the damaging effects of pressure surges. This becomes necessary where conventional means of mitigating surge pressures cannot be employed such as when handling radioactive, highly corrosive or lethal fluids, where no fluid is allowed to escape. The research concluded that increase in pipe wall thickness, flange rating and pipe supports can be designed to prevent catastrophic failure. In order to prevent an increase in surge fatigue damage, devices such as variable speed drives for pumps and slow closing valves should be considered. Also rerouting of pipelines can avoid a profile that is conducive to column separation. For a particular pipeline, it may be possible to use a thermoplastic or GRP material rather than a ferrous pipe material. This applies to low head pipelines found in the mining water and waste water industry where high temperatures do not occur. Therefore, the reduced modulus, results in a reduced celerity or wave speed. The modulus is the prime variable in establishing the celerity.</p><p>Pressure surges in fluid systems often lead to pipe damages [<xref ref-type="bibr" rid="scirp.57298-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.57298-ref3">3</xref>] . They are caused by an abrupt change in flow speed and can have amplitudes higher than the usual pressure conditions in normal operation. Preventing pressure surges is an important consideration in the design and reliable operation of pipes. The analysis concluded that changes in pipe system such as an increase in through put quantities or a change of operating mode can bring about the danger of pressure surges that may excessively stress the system. The research analyzed entire pipe system with regard to pressure surges. On the basis of its findings, the suggestions for solutions include adjustment of the valve closing times, additional air vessels, and use of bypass valves, ventilation or bleeding valves, Non-return valves and pump fly wheel mass.</p><p>Cla-val [<xref ref-type="bibr" rid="scirp.57298-ref4">4</xref>] carried out research on complex pipelines with undulating terrain where column separation in the pipe may occur. Hence, the research stated that detail pipeline profile as well as air valve size and placement when pump stops due to power failure causes down surge wave which travels down the pipe. Low pressure wave is reflected at the pipe end boundary back to the pump with a resulting high pressure spike generated at the pump. This high pressure spike damages the piping. A test was carried out using a standard relief valve. The standard relief valve could not open fast enough to dissipate the excess pressure spike but on testing with cal-val 52-03 surge anticipator valve, it opened on the initial down surge pressure and stays open to dissipate the returning high pressure wave.</p><p>Jung, et al. [<xref ref-type="bibr" rid="scirp.57298-ref5">5</xref>] showed in their research that surge modeling is important to safeguard against breaches in distribution system integrity. Several cases were considered in this research. Hence the surge analysis results of the pressure-sensitive demand model were compared with those of pressure-insensitive (constant) demand formulations.</p><p>This paper reviews the engineering analysis behind pressure surge propagation in a pipeline system. The analysis is carried out on a 254 mm pipe of length 6 kilometers from the emergency relief-coupling valve (ERV) installation point. The pipeline carries crude oil whose flow is aided by pumps from a flow station. The pipe terminates at a loading point. The system is made up of two valves namely the normal-shut off valve (NSV) and the emergency relief coupling valve (ERV). Of interest to this project is the closure time of the ERV at time, t and the pressure changes which accompanies this closure. The pressure change in the system results to two basic anomalies namely surge pressure and cavitations or column separation. These two anomalies threaten the operation of pipeline system especially pipe machines such as pumps, fittings and also leads to spillage.</p></sec><sec id="s2"><title>2. Materials and Methods</title><p>Data were collected from direct measurements, log sheets for ten months, and relevant charts and manuals. Some of the pertinent parameters were: crude oil density, pipeline length, pipe diameter, flow rate, flow pressure, pipe roughness size, fittings, viscosity of crude oil, Young’s Modulus of the pipe, emergency-relief coupling valve installed closure or shutdown time.</p><p>In laminar flow, frictional forces are due to viscous drag and are independent of pipe conditions but in turbulent flow, viscous shear predominates. Total friction is affected by pipe surface roughness. The pipeline under investigation undergoes turbulent flow. Therefore, friction and head losses in the pipe are calculated using the following equations [<xref ref-type="bibr" rid="scirp.57298-ref6">6</xref>] :</p><disp-formula id="scirp.57298-formula248"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x5.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula249"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x6.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula250"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x7.png"  xlink:type="simple"/></disp-formula><p>The following pressure wave equations for pipe flow were used to evaluate pressure changes due to gradual valve closure by the emergency-relief coupling valve [<xref ref-type="bibr" rid="scirp.57298-ref7">7</xref>] :</p><p>Retardation of crude oil (fluid)</p><disp-formula id="scirp.57298-formula251"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x8.png"  xlink:type="simple"/></disp-formula><p>Force available for producing retardation = mass &#215; retardation</p><disp-formula id="scirp.57298-formula252"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x9.png"  xlink:type="simple"/></disp-formula><p>Force due to pressure wave</p><disp-formula id="scirp.57298-formula253"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x10.png"  xlink:type="simple"/></disp-formula><p>Therefore,</p><disp-formula id="scirp.57298-formula254"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x11.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula255"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x12.png"  xlink:type="simple"/></disp-formula><p>But</p><disp-formula id="scirp.57298-formula256"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x13.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula257"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x14.png"  xlink:type="simple"/></disp-formula><p>Therefore,</p><disp-formula id="scirp.57298-formula258"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x15.png"  xlink:type="simple"/></disp-formula><p>Therefore, Equation (8) becomes</p><disp-formula id="scirp.57298-formula259"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x16.png"  xlink:type="simple"/></disp-formula><p>Equation (12) is used to model the pressure change at any point in the pipeline arising from shut down of pipeline by emergency-relief coupling valve (ERV).</p><p>The Reynolds number of the fluid following valve closure at time t is deduced using equation below;</p><disp-formula id="scirp.57298-formula260"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x17.png"  xlink:type="simple"/></disp-formula><p>The wave speed of the system is calculated using equation below [<xref ref-type="bibr" rid="scirp.57298-ref6">6</xref>] :</p><disp-formula id="scirp.57298-formula261"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x18.png"  xlink:type="simple"/></disp-formula><p>Unbalanced hydraulic forces are normally created by pressure surge phenomenon and are particularly dangerous. The following equations are used to calculate the magnitude of hydraulic forces propagated [<xref ref-type="bibr" rid="scirp.57298-ref7">7</xref>] :</p><p>Force generated in the 45 degrees elbow or bend;</p><disp-formula id="scirp.57298-formula262"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x19.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula263"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x20.png"  xlink:type="simple"/></disp-formula><p>Force generated in 90 degrees elbow or bend;</p><disp-formula id="scirp.57298-formula264"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x21.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula265"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x22.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula266"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x23.png"  xlink:type="simple"/></disp-formula><p>When the valve (ERV) is shut at time, t the crude oil moving with velocity, v, at a temperature, T, is brought to rest adiabatically, hence reaching a final temperature, T<sub>O</sub>, when at rest. The flow equation becomes [<xref ref-type="bibr" rid="scirp.57298-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.57298-ref9">9</xref>] :</p><disp-formula id="scirp.57298-formula267"><graphic  xlink:href="http://html.scirp.org/file/3-8102391x24.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.57298-formula268"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x25.png"  xlink:type="simple"/></disp-formula><p>At the valve, the stagnation pressure, P<sub>o</sub>, would be:</p><disp-formula id="scirp.57298-formula269"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x26.png"  xlink:type="simple"/></disp-formula><p>Also, the crude oil would have enthalpy:</p><disp-formula id="scirp.57298-formula270"><label>(22)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/3-8102391x27.png"  xlink:type="simple"/></disp-formula></sec><sec id="s3"><title>3. Results and Discussions</title><p>The results are presented in Tables 1-12.</p><p>From the results shown above, the flow velocity is 1.12 m/s but upon shut down of the system in 1.5 seconds, the velocity of flow is 53.7 m/s which is as a result of the high pressure of 37.23 bar attained as shown in <xref ref-type="table" rid="table4">Table 4</xref>. But as the valve closure rate decreased, there is a decrease in the pressure rise, the final flow velocity, V<sub>1</sub> also decreased. The remarkable increase in final velocity V<sub>1</sub> in the pipe at different closure rate can result to several changes in fluid parameters. From the result, the pipe recorded the least final velocity V<sub>1</sub> of 25.1 m/s when the</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Fluid data</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Fluid density, ρ (Kg/m<sup>3</sup>)</th><th align="center" valign="middle" >Flow rate Q (m<sup>3</sup>/s)</th><th align="center" valign="middle" >Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Viscosity of fluid &#181; (Kg/ms)</th></tr></thead><tr><td align="center" valign="middle" >847</td><td align="center" valign="middle" >0.0557</td><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4.15 &#215; 10<sup>−3</sup></td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Pipe data</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Pipe length L (m)</th><th align="center" valign="middle" >Pipe diameter d (m)</th><th align="center" valign="middle" >Wall roughness size e (m)</th><th align="center" valign="middle" >k-factor 90˚ elbow</th><th align="center" valign="middle" >k-factor 45˚ elbow</th></tr></thead><tr><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >0.254</td><td align="center" valign="middle" >0.0002</td><td align="center" valign="middle" >0.9</td><td align="center" valign="middle" >0.4</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Result of Reynolds number re of the various fluids in varied valve closure time (s)</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="2"  >Bonny Light</th><th align="center" valign="middle"  colspan="2"  >Bonny medium</th><th align="center" valign="middle"  colspan="2"  >Water</th></tr></thead><tr><td align="center" valign="middle" >Re</td><td align="center" valign="middle" >t (s)</td><td align="center" valign="middle" >Re</td><td align="center" valign="middle" >t (s)</td><td align="center" valign="middle" >Re</td><td align="center" valign="middle" >t (s)</td></tr><tr><td align="center" valign="middle" >2.1 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >1.5</td><td align="center" valign="middle" >3.5 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >1.5</td><td align="center" valign="middle" >9.0 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >1.5</td></tr><tr><td align="center" valign="middle" >1.2 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >2.5</td><td align="center" valign="middle" >2.5 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >2.5</td><td align="center" valign="middle" >5.3 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >2.5</td></tr><tr><td align="center" valign="middle" >1.0 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >1.5 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >3</td><td align="center" valign="middle" >4.5 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >3</td></tr><tr><td align="center" valign="middle" >8.0 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >9.3 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >4</td><td align="center" valign="middle" >3.3 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >6.4 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >7.2 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >5</td><td align="center" valign="middle" >2.7 &#215; 10<sup>8</sup></td><td align="center" valign="middle" >5</td></tr></tbody></table></table-wrap><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Pressure and flow rate changes along the pipe in 1.5 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >744757.644</td><td align="center" valign="middle" >1755214</td><td align="center" valign="middle" >2.93E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >1489515.3</td><td align="center" valign="middle" >1010484.7</td><td align="center" valign="middle" >1.42E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >2234272.9</td><td align="center" valign="middle" >265726.1</td><td align="center" valign="middle" >3.37E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >2979030.5</td><td align="center" valign="middle" >−479031.2</td><td align="center" valign="middle" >−5.65E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >3723788.22</td><td align="center" valign="middle" >−1223788.2</td><td align="center" valign="middle" >−1.37E+13</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Pressure and Flow rate changes along the pipe in 2.5 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >68445.4</td><td align="center" valign="middle" >2431554.1</td><td align="center" valign="middle" >2.11E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >1236893.2</td><td align="center" valign="middle" >1263106.8</td><td align="center" valign="middle" >1.31E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >1816345.3</td><td align="center" valign="middle" >683654.7</td><td align="center" valign="middle" >2.18E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >2039443.2</td><td align="center" valign="middle" >460556.8</td><td align="center" valign="middle" >3.16E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >2234053.1</td><td align="center" valign="middle" >265946.9</td><td align="center" valign="middle" >1.13E+13</td></tr></tbody></table></table-wrap><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Pressure and Flow rate changes along the pipe in 1.5 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >744757.644</td><td align="center" valign="middle" >1755214</td><td align="center" valign="middle" >2.93E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >1489515.3</td><td align="center" valign="middle" >1010484.7</td><td align="center" valign="middle" >1.42E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >2234272.9</td><td align="center" valign="middle" >265726.1</td><td align="center" valign="middle" >3.37E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >2979030.5</td><td align="center" valign="middle" >−479031.2</td><td align="center" valign="middle" >−5.65E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >3723788.22</td><td align="center" valign="middle" >−1223788.2</td><td align="center" valign="middle" >−1.37E+13</td></tr></tbody></table></table-wrap><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> Pressure and Flow rate changes along the pipe in 2.5 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >68445.4</td><td align="center" valign="middle" >2431554.1</td><td align="center" valign="middle" >2.11E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >1236893.2</td><td align="center" valign="middle" >1263106.8</td><td align="center" valign="middle" >1.31E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >1816345.3</td><td align="center" valign="middle" >683654.7</td><td align="center" valign="middle" >2.18E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >2039443.2</td><td align="center" valign="middle" >460556.8</td><td align="center" valign="middle" >3.16E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >2234053.1</td><td align="center" valign="middle" >265946.9</td><td align="center" valign="middle" >1.13E+13</td></tr></tbody></table></table-wrap><table-wrap id="table8" ><label><xref ref-type="table" rid="table8">Table 8</xref></label><caption><title> Pressure and Flow rate changes along the pipe in 3 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >522289.8</td><td align="center" valign="middle" >2447710.2</td><td align="center" valign="middle" >2.41E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >1144993.5</td><td align="center" valign="middle" >1355006.5</td><td align="center" valign="middle" >3.66E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >1511922.4</td><td align="center" valign="middle" >988077.6</td><td align="center" valign="middle" >2.00E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >1623344.9</td><td align="center" valign="middle" >876655.1</td><td align="center" valign="middle" >2.16E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >1861234.3</td><td align="center" valign="middle" >688765.7</td><td align="center" valign="middle" >1.33E+13</td></tr></tbody></table></table-wrap><table-wrap id="table9" ><label><xref ref-type="table" rid="table9">Table 9</xref></label><caption><title> Pressure and Flow rate changes along the pipe in 4 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >40138.7</td><td align="center" valign="middle" >2459861.3</td><td align="center" valign="middle" >1.11E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >924778.1</td><td align="center" valign="middle" >1575221.9</td><td align="center" valign="middle" >2.66E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >1023496.3</td><td align="center" valign="middle" >1496503.7</td><td align="center" valign="middle" >1.56E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >1122496.4</td><td align="center" valign="middle" >1377503.6</td><td align="center" valign="middle" >2.37E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >1396445.3</td><td align="center" valign="middle" >1103554.7</td><td align="center" valign="middle" >3.49E+11</td></tr></tbody></table></table-wrap><table-wrap id="table10" ><label><xref ref-type="table" rid="table1">Table 1</xref>0</label><caption><title> Pressure and Flow rate changes along the pipe in 5 s valve shut down time</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Initial Pressure P (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Length of Pipe (m)</th><th align="center" valign="middle" >Pressure Change P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Pressure Difference P-P<sub>L</sub> (N/m<sup>2</sup>)</th><th align="center" valign="middle" >Flow Rate Change Q<sub>L</sub> (m<sup>3</sup>/s)</th></tr></thead><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >1200</td><td align="center" valign="middle" >41289.4</td><td align="center" valign="middle" >2458710.6</td><td align="center" valign="middle" >2.0E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >2400</td><td align="center" valign="middle" >712594.6</td><td align="center" valign="middle" >1787405.4</td><td align="center" valign="middle" >2.11E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >3600</td><td align="center" valign="middle" >946438.7</td><td align="center" valign="middle" >1553561.3</td><td align="center" valign="middle" >3.08E+12</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >4800</td><td align="center" valign="middle" >1247891.3</td><td align="center" valign="middle" >1252108.7</td><td align="center" valign="middle" >3.44E+13</td></tr><tr><td align="center" valign="middle" >25 &#215; 10<sup>5</sup></td><td align="center" valign="middle" >6000</td><td align="center" valign="middle" >1355663.2</td><td align="center" valign="middle" >1144336.8</td><td align="center" valign="middle" >1.0E+13</td></tr></tbody></table></table-wrap><table-wrap id="table11" ><label><xref ref-type="table" rid="table1">Table 1</xref>1</label><caption><title> Hydraulic force F (KN) result at varied valve closure time t (s)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Valve closure time, t (s)</th><th align="center" valign="middle" >Force, F (45 deg) (KN)</th><th align="center" valign="middle" >Force, F (90 deg (KN)</th><th align="center" valign="middle" >V<sub>1</sub> (m/s)</th></tr></thead><tr><td align="center" valign="middle" >1.5</td><td align="center" valign="middle" >18291.4</td><td align="center" valign="middle" >34651</td><td align="center" valign="middle" >53.7</td></tr><tr><td align="center" valign="middle" >2.5</td><td align="center" valign="middle" >6544.1</td><td align="center" valign="middle" >12552</td><td align="center" valign="middle" >25.1</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >4276.4</td><td align="center" valign="middle" >8552.7</td><td align="center" valign="middle" >38.9</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >2509.9</td><td align="center" valign="middle" >4670.9</td><td align="center" valign="middle" >51.2</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >1758.8</td><td align="center" valign="middle" >2906.5</td><td align="center" valign="middle" >57.2</td></tr></tbody></table></table-wrap><table-wrap id="table12" ><label><xref ref-type="table" rid="table1">Table 1</xref>2</label><caption><title> Stagnation values at varied closure time (Bonny light)</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Valve closure time, t (s)</th><th align="center" valign="middle" >V (m/s)</th><th align="center" valign="middle" >T<sub>o</sub> (k)</th><th align="center" valign="middle" >P<sub>O</sub> (m/s)</th><th align="center" valign="middle" >H<sub>O</sub> (KJ/kg)</th></tr></thead><tr><td align="center" valign="middle" >1.5</td><td align="center" valign="middle" >53.7</td><td align="center" valign="middle" >2174.4</td><td align="center" valign="middle" >12.2</td><td align="center" valign="middle" >1515.9</td></tr><tr><td align="center" valign="middle" >2.5</td><td align="center" valign="middle" >25.1</td><td align="center" valign="middle" >834.6</td><td align="center" valign="middle" >2.7</td><td align="center" valign="middle" >389.11</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >38.9</td><td align="center" valign="middle" >1359.6</td><td align="center" valign="middle" >6.4</td><td align="center" valign="middle" >830.7</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >51.0</td><td align="center" valign="middle" >2006.4</td><td align="center" valign="middle" >11.0</td><td align="center" valign="middle" >1374.6</td></tr><tr><td align="center" valign="middle" >5</td><td align="center" valign="middle" >57.2</td><td align="center" valign="middle" >2405.2</td><td align="center" valign="middle" >13.0</td><td align="center" valign="middle" >1710.0</td></tr></tbody></table></table-wrap><p>valve closed in 2.5 seconds as shown in <xref ref-type="table" rid="table1">Table 1</xref>1. This velocity should be adjudged to the best valve closure rate suitable to maintain forty one thousand barrels per day capacity (0.0557 m/s) of the system.</p><p>But on reduction of the flow rate (volume) and the valve maintains 2.5 seconds closure time, the velocity dropped considerably which indicates that the pipeline operates at high flow rate initially. The modeling shows that for any volume flow rate, Q m<sup>3</sup>/s there is a valve closure time suitable for the pipeline in which pressure changes can be tolerated even at surge suppression device failure. The capacity of the pipeline under investigation is 41,000 barrels/day (0.0557 m<sup>3</sup>/s) operating at 25 bar initial pressure and installed emergency relief- coupling valve closure time of 1.5 seconds, upon the shutdown of the system by the valve in response to an emergency, the pressure attained at 6000 m length of pipe is 37.23 bar as against 25 bar operating pressure despite the pressure loss in the flow. This pressure attained on closure at 1.5 seconds is very high for the system and leads the system to surge. Consider the result in <xref ref-type="table" rid="table1">Table 1</xref>1, the flow velocity is 1.12 m/s but upon shut down of the system in 1.5 seconds, the velocity of flow is 53.7 m/s which is as a result of the high pressure of 37.23 bar attained as shown in <xref ref-type="table" rid="table4">Table 4</xref>. But as the valve closure rate decreased, there is a decrease in the pressure rise, the final flow velocity, V<sub>1</sub> also decreased. The remarkable increase in final velocity V<sub>1</sub> in the pipe at different closure rate can result to several changes in fluid parameters. From the result, the pipe recorded the least final velocity V<sub>1</sub> of 25.1 m/s when the valve closed in 2.5 seconds as shown in <xref ref-type="table" rid="table1">Table 1</xref>1. This velocity should be adjudged to the best valve closure rate suitable to maintain forty one thousand barrels per day capacity (0.0557 m/s) of the system. At 4.0 seconds, it experienced a gradual fall in the hydraulic force. It can be concluded that hydraulic force propagation decreases with an increase in valve closure time. The same trend followed in the 45˚ elbow, though the highest hydraulic force is 18,000 KN when the valve closed in 1.1 seconds whereas the 90˚ elbow has the highest hydraulic force value of 37,000 KN when the valve closed in 1.5 seconds. It can also be deduced that hydraulic force propagation depends on the type of elbow used in the pipeline. Therefore, the research does not advice against the use of 90˚ elbows since its usage depends on the slope of the soil. Though it is suggested that instead of using more 90˚ elbow in pipeline transmission, effort should be made to use more 45˚ elbow.</p><p>Consider what happens if the inlet flow rate (volume) was increased to 0.06109 m<sup>3</sup>/s and the valve set to close in 2.5 seconds. The maximum pressure rise was 24.2 bar and a decrease in flow rate of 0.04752 m<sup>3</sup>/s, the maximum pressure rise was 18.8 bar. It shows that at an increased flow rate, at time t, the pressure change would be increasing capable of surge pressure. But at a decreased flow rate at time t, the pressure change would be decreasing capable of falling below vapor pressure of the crude oil. Hence, modeling the closure time t at a given flow rate is quite necessary for satisfactory result. As stated earlier, there is always a closure time limit for the valve over which pressure changes can be tolerated owing to level of pressure losses in the system. This is particularly shown in <xref ref-type="table" rid="table4">Table 4</xref> where a flow rate of 0.0557 m<sup>3</sup>/s of crude oil was maintained at different valve closure time. Hence, the analysis shows that the limit between 2.5 seconds to 3.0 seconds valve closure rate can be tolerated. This fact has not been considered previously owing to over-dependence on surge suppression devices which on its own has shown high failure rate. The failure of these devices has led to adverse environmental anomalies like spillage, excessive cost of maintenance of pipeline and pipe machines.</p></sec><sec id="s4"><title>4. Conclusion</title><p>Results of the analysis identified valve operation as a major cause of surge propagation. Surge is a potentially damaging flow instability that limits the flow operation of pipe machines. The rate at which the emergency relief-coupling valve (ERV) closes is significant in pressure surge propagation. At faster closure rate, the pressure of the system rises such that can destroy the pipeline mechanically, while at slower valve closure rate, there is high tendency of vacuum separation in the system due to a remarkable fall in the pressure of the system. Hence, there is always a limit of valve closure time suitable for any pipeline at any given flow rate. Hydraulic forces propagation decreases with increasing emergency relief-coupling valve (ERV) closure time. Leakages from fittings (bends) are mostly due to the propagation of hydraulic forces arising from valve closure. Pressure change in a pipeline installed with emergency relief-coupling valve increases with, fluid density, flow rate (volume) pipe length and decreases with pipe diameter and valve closure rate upon shut down. The fluid Reynolds number decreases with increasing emergency relief-coupling valve closure rate. As a result of the movement of the elbow, created by the shutdown of ERV at time t, the joint will become weaker and in extreme cases, they may come apart, hence enable air to be trapped into the system. An emergency relief-coupling valve (ERV) should be installed in a system to guard against storm; hence proper design of its operation should be considered before installation.</p></sec><sec id="s5"><title>Notations</title><p>A Area of pipe, m<sup>2</sup></p><p>C Shock wave speed, m/s</p><p>C<sub>1</sub> Pipe distensibility (restrain factor)</p><p>C<sub>p</sub><sub> </sub> Specific heat of fluid at constant pressure, kJ/kgK</p><p>d Diameter of pipe, m</p><p>E<sub>y</sub> Young’s Modulus for the pipe material, N/m<sup>2</sup></p><p>ERV Emergency relief coupling valve</p><p>g Acceleration due to gravity, m/s<sup>2</sup></p><p>H<sub>f</sub> Head loss due to friction, m</p><p>K Fluid bulk Modulus, N/m<sup>2</sup></p><p>L Length of pipe, m</p><p>ṁ Mass flow rate, Kg/s</p><p>NSV Normal shut-off valve</p><p>P Pipe pressure, N/m<sup>2 </sup></p><p>P<sub>o</sub> Stagnation pressure, N/m<sup>2</sup></p><p>Re Reynolds number</p><p>t Period (time), s</p><p>t<sub>h</sub> Pipe thickness, m</p><p>T Temperature, K</p><p>T<sub>o</sub> Stagnation temperature, K</p><p>V Flow velocity, m/s</p><p>e Pipe Roughness Size</p><p>F Force, N</p><p>F<sub>x</sub> Force in direction of x, N</p><p>F<sub>y</sub> Force in direction of y, N</p><p>f Pipe Friction Factor</p><p>H Total Head Loss m</p><p>k Fitting Loss Coefficient</p><p>K<sub>1 </sub>Bulk Modulus of Fluid<sup> </sup></p><p>Q Volumetric Flow Rate, m<sup>3</sup>/s</p><p>Q<sub>L</sub><sub> </sub>Change in Flow Rate, m<sup>3</sup>/s</p><p>R<sub>L</sub><sub> </sub>Pressure ratio</p><p>t Time, s</p><p>t<sub>n</sub> Pipe thickness, m</p></sec><sec id="s6"><title>Greek Letters</title><p>ρ Density of the Fluid, kg/m<sup>3</sup></p><p>&#181; Dynamic Viscosity of Fluid, kg/m/s</p><p>τ Wall Shear Stress, N/m<sup>2 </sup></p></sec><sec id="s7"><title>Suffixes</title><p>1. Initial state</p><p>2. Final state</p></sec></body><back><ref-list><title>References</title><ref id="scirp.57298-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Stone, G.D. (2013) Avoiding Pressure Surge Damage in Pipeline Systems. SE Asia, Feb Publishers, Australia.</mixed-citation></ref><ref id="scirp.57298-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Voith Industrial Services (2014) Pressure Surges: Analysis and Prevention. www.voith.com</mixed-citation></ref><ref id="scirp.57298-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Lebele-Alawa, B.T. and Oparadike, F.E. (2011) Analysis of the Effects of Valve Propagated Pressure Surge on Pipe Flow. Engineering, 3, 1098-1101. http://dx.doi.org/10.4236/eng.2011.311137</mixed-citation></ref><ref id="scirp.57298-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Cla-Val, A. (2013) Surge Analysis and Control. www.cal-val.com</mixed-citation></ref><ref id="scirp.57298-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Jung, B.S., Boulos, P.F. and Don, J.W. (2009) Effect of Pressure-Sensitive Demand on Surge Analysis. Journal of American Water Works Association, (AWWA), 101, 103-111.</mixed-citation></ref><ref id="scirp.57298-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Sunrise Systems Customer Education Training Manual (2008) Modeling Pipeline and Utility Systems Using Pipenet Vision. www.sunrise-sys.com</mixed-citation></ref><ref id="scirp.57298-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Rajput, R.K. (2004) A Textbook of Fluid Mechanics and Hydraulic Machines. S. Chand and Company Ltd., New Delhi.</mixed-citation></ref><ref id="scirp.57298-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Eastop, T.D. and McConkey, A. (1993) Applied Thermodynamics for Engineering Technologists. Prentice Hall, Singapore.</mixed-citation></ref><ref id="scirp.57298-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Lebele-Alawa, B.T. (2007) Effects of Compressor-Blade Profile Change on Gas-Turbine Performance. Ph.D. Thesis, Rivers State University of Science and Technology, Port Harcourt.</mixed-citation></ref></ref-list></back></article>