Water Hammer

Download as pdf or txt
Download as pdf or txt
You are on page 1of 49

Analysis for Water Hammer considering the effect of Fluid Structure Interaction in Straight Pipes

A THESIS SUBMITTED IN THE PARTIAL FULFILLMENTS OF THE REQUIREMENTS FOR THE DEGREE OF BACHELOR OF TECHNOLOGY In MECHANICAL ENGINEERING

By

SUYASH MISHRA Roll number 108ME050

Department of Mechanical Engineering National Institute of Technology, Rourkela May 2012

Analysis for Water Hammer considering the effect of Fluid Structure Interaction in Straight Pipes

A THESIS SUBMITTED IN THE PARTIAL FULFILLMENTS OF THE REQUIREMENTS FOR THE DEGREE OF BACHELOR OF TECHNOLOGY In MECHANICAL ENGINEERING

By

SUYASH MISHRA Roll number 108ME050

Under the Guidance of


Prof. J. SRINIVAS

Department of Mechanical Engineering National Institute of Technology, Rourkela May 2012

National Institute of Technology Rourkela

CERTIFICATE
This is to certify that the work in this thesis entitled Analysis for Water Hammer considering the effect of Fluid Structure Interaction in Straight Pipes by Suyash Mishra (108ME050), has been carried out under my supervision in partial fulfillment of the requirements for the degree of Bachelor of Technology in Mechanical Engineering during session 2011- 2012 in the Department of Mechanical Engineering, National Institute of Technology, Rourkela. To the best of my knowledge, this work has not been submitted to any other University/Institute for the award of any degree or diploma.

Date 10th May, 2012 Place- Rourkela Prof. J. Srinivas (Supervisor) Associate Professor Dept. of Mechanical Engineering National Institute of Technology Rourkela

Mechanical Engineering 08-12

ACKNOWLEDGEMENT

I would like to express my deep sense of gratitude and respect to my supervisor Prof. J.Srinivas for his excellent guidance, valuable suggestions and endless support. I am also thankful to Prof. S. Chakraborty, Department of Mathematics, NIT Rourkela for his constant support and encouragement. Last, but not the least I extend my sincere thanks to other faculty members of the Department of Mechanical Engineering, NIT Rourkela and my dear friends for their valuable advice in every stage for successful completion of this project report.

SUYASH MISHRA Roll no-108ME050 Mechanical Engineering Department National Institute of Technology, Rourkela

Mechanical Engineering 08-12

II

CONTENTS
Certificate Acknowledgement List of Figures and Tables Abstract List of Figures and Table Nomenclature 1. INTRODUCTION 1.1 1.2 1.3 Problem statement Literature Review Scope and Objectives of the work 1 4 7 I II IV V VI VII

2. MATHEMATICAL MODELLING 2.1 2.2 Classical water hammer theory Modification in classical water hammer equation 2.2.1 2.2.2 2.3 Governing Equation of fluid Governing Equation of solid(Pipe) 8 9 9 10

Initial and Boundary Condition 2.3.1 2.3.2 Initial Condition Boundary Condition 11 11

2.4

Numerical Solution 2.4.1 2.4.2 2.4.3 2.4.4 2.4.5 Finite element Method Displacement and Velocity Potential Formulation Displacement, Velocity Potential and Pressure Formulation Interpolation Function Matrix Element formulation in terms of shape function 11 12 14 15 16
III

Mechanical Engineering 08-12

2.4.6

Solution

17

3. RESULTS and DISCUSSION 3.1 Results for Displacement, Velocity Potential and Pressure Formulation 3.1.1 3.1.2 3.1.3 3.1.4 3.1.5 3.2 Displacement Vs Time Velocity Vs Time Pressure Vs Time Force on solid(pipe) Vs Time Force on Fluid Vs Time 19 20 21 23 24

Results and observation for Displacement and Velocity Potential 3.2.1 3.2.2 3.2.3 3.2.4 Displacement Vs Time Velocity Vs Time Force on Solid( Pipe) Vs Time Force on Liquid Vs Time 25 25 26 26 28 32

3.3 3.4

Comparison of individual effect of coupling Comparison with u- and u--P

4. CONCLUSIONS 4.1 4.2 Summary Future scope of work 33 34 35 36

REFRENCES APPENDIX

Mechanical Engineering 08-12

IV

ABSTRACT
Water hammer is a serious problem in flow pipelines of nuclear power plant as fluid structure interaction significantly affects the response of pipe components and performance of such plants as a whole. This phenomenon is basically observed when there is sudden increase in flow velocity, abrupt change in cross section and if the flow is forced fully to stop. The purpose of this study is to use a finite element approach for studying such a highly nonlinear problem in order to investigate the effect of fluid structure interaction (FSI) in pipe lines. In this study, the transient water hammer problem is attempted by considering the effect of all three Poisson, friction and junction couplings. The continuity and momentum equations are written in terms of pressure and velocity of the flow. Along with this the pipe-structural equations in axial direction is considered. Compared to conventional water hammer relations, present work accounts the Poisson coupling as an important term in the fluid dynamic equations. The complete interaction of fluid and pipe comes with boundary condition, which is taken as gradual valve closure. There are various solution approaches such as method of characteristics (MOC), finite element-MOC mixed method as well as finite element approach alone. In present work, finite element approach is employed to mesh the fluid and solid regions. In FEM method two kind of formulation has been done to account for symmetric matrix formulation and shape function used in determining the mass stiffness and damping matrix element. This kind of formulation uses velocity potential () and hydrostatic pressure (p) as unknowns in fluid domain and axial displacement (u) as unknown in solid domain. A MATLAB code has been written to calculate result for Finite element formulation by solving the nonlinear second-order differential equations using Newmark-beta scheme.

Mechanical Engineering 08-12

LIST OF FIGURE
Page no. Fig.1 Symmetric diagram for interaction and excitation in liquid pipe system Fig.2.4.2.1 Pipe element with two degree of freedom at each node Fig2.4.3 Pipe element with 3 degree of freedom at three nodes Fig.3.1.1 Displacement (m) Vs Time (sec) for u--P method for 3 nodes Fig.3.1.2 Velocity (m/s) Vs Time (sec) for u--P method for 3 nodes Fig.3.1.3 Pressure Variation between classical water hammer and Fem Approach Fig.3.1.4 Pressure (N/m^2) Vs Time (sec) for u--P method for 3 nodes Fig.3.1.5 Force on structure (N) Vs Time (sec) for u--P method for 3 nodes Fig.3.1.6 Force on fluid (N) Vs Time (sec) for u--P method for 3 nodes Fig3.2.1 Displacement (m) Vs Time (sec) for u- method Fig.3.2.2 Velocity (m/sec) Vs Time (sec) for u- method Fig.3.2.3 Force on structure (N) Vs Time (sec) for u- method Fig.3.2.4Force on fluid (N) Vs Time (sec) for u- method 3 14 15 19 20 21 21 22 23 24 24 25 25

Fig.3.3.1 Variation of Displacement (m) Vs Time (sec) for u- method with consideration of all coupling 28 Fig.3.3.2 Variation of Velocity (m/s) Vs Time (sec) for u- method with consideration of all coupling 28 Fig.3.3.3 Variation of Force on structure (m) Vs Time (sec) for u- method with consideration of all coupling 29 Fig.3.3.4 Variation of Force on fluid (N) Vs Time (sec) for u- method with consideration of all coupling 29 Fig.3.4.1 Comparison of Displacement (m) Vs Time (sec) for u- and u--P method with consideration of all coupling 30 Table2.6.1 Specification of piping system
Mechanical Engineering 08-12
17

VI

NOMENCLATURE
E- Young modulus (N/m^2) Kf -Bulk Modulus (N/m^2) D-Diameter (m) e- Thickness (m) f- Friction factor L-Length (m) C0-Initial wave Velocity (m/sec^2) Ms- Mass of pipe (Kg/m) Mf-Mass of fluid (Kg/m) - Density (Kg/m^3) FsForce due to fluid on structure Ff-Force due to structure on fluid Fv-Force due to Valve Cv- Valve closure coefficient f- Density of fluid(Kg/m3) s- Density of solid(pipe) (Kg/m3) - Axial stress in pipe (N/m2) - Poisson ratio - Strain Tensor - Velocity Potential - Pressure Ac-Area of cross section (m2) As- Area of Solid(Pipe) (m2) Af- Area of Pipe (m2)
Mechanical Engineering 08-12 VII

Chapter 1 INTRODUCTION
1.1 Proble m State ment Study of water hammer phenomenon has become very important not only in terms of safety reason but also from the economics point of view as any mishap can result to low performance of plant. Various attempt like analytical modeling, computational modeling and experimental modeling has been made to get a greater insight of the problem but still it is most challenging problem to deal with because of its highly non linear nature and coupling between two domain of physics solid and fluid. Water Hammer is phenomenon of sudden rise of pressure due to closing of valve or due to any other similar cause and subsequent propagating of wave along the pipe causing noise and excitation. The rise of pressure may ever burst the pipe. The magnitude of pressure depends upon speed of valve closure, velocity of flow, length of pipe, elasticity of pipe and flowing fluid. A classical water hammer approach has been established which describe the phenomenon in more general way and predicts correctly extreme pressure and velocity of the wave but this theory could not predict or measure damping and sudden change in wave traversing along the pipe. The reason behind above failure of theory is non-consideration of different coupling mechanism like friction, Poisson and junction coupling which governs the exact nature of water hammer phenomenon.

Mechanical Engineering 08-12

The waves exert acoustic pressure on the solid and the solid and fluid are coupled because of this force exerted by fluid flow on solid. This fluid flow results in deformation of structure and deformation of structure cause the change in fluid flow. This coupled deformation of fluid and solid is term as fluid structure interaction. All the three coupling play a major role in above mechanism. Water hammer phenomenon is caused by number of phenomenon mentioned below: 1. Pressure wave in single phase flow 2. Water column separation 3. Flow into voided section 4. Acceleration of a water slug in steam lines 5. Steam condensation induced water hammer 6. Two phase transient 7. Valve Dynamics/ Instabilities In all above major reason for water hammer in plant setup we will confine our self in dealing with valve dynamics / instabilities problem for single phase flow which is also termed as Delft Hydraulics Benchmark Problem. The FSI has got great attention in recent years because of safety issue, reliability of plant set up, environmental concern in pipe delivery system and plant performance. The Fluid Structure interaction process basically deals with transfer of momentum and forces to pipe system and fluid contained in it in an unsteady. The excitation process may be caused by sudden change in flow and pressure or by some mechanical action namely sudden closure of valve. This

Mechanical Engineering 08-12

interaction result in vibration of pipe and abrupt increase in velocity and pressure of fluid. The resulting load is transferred to pipe supporting system. There are three type of coupling governing fluid structure interaction in pipe system 1. Friction Coupling 2. Poisson Coupling 3. Junction coupling Friction coupling is due to transient liquid shear stresses resisting relatively axial motion between fluid and pipe wall. These stress act at interface between fluid and pipe. The effect of this coupling as compared to other two is less significant. Poisson coupling is associated to hoop or circumferential stress deflection produced by liquid pressure transient which is transferred to axial stress in the ration of Poisson ratio. The axial stress and axial strain deflection traverse along pipe with speed of sound approximately. Junction coupling is most important coupling among three and its impact in most significant. It results because of reaction set up by unbalanced pressure forces and by change in fluid momentum at different position in pipe s uch as valve, bend etc. Sources of excitation are not concentrated only on liquid part but also they are equally contributed from solid side as in case of valve motion the change in structure moment results in liquid transient and cause vibration. Fig.1 shows the bi-directional effect of these coupling on fluid and pipe system.
Friction, Poisson and Junction coupling

Valve Action

Fluid

Pipe

Pipe Rupture

Fig.1 Symmetric diagram for interaction and excitation in liquid pipe system

Mechanical Engineering 08-12

There may be a deviation from practical data due to several other factors. The highly non linear nature of problem as it is coupling to physical world make the solution more tough and it get worst if we consider highly turbulent flow. The Navier-strokes equation for fluid and Euler equation for pipe has to be solved simultaneously with approximate boundary condition. In this project we have use a finite element method to solve this problem showing the effect of each coupling individually. 1.2 Lite rature Review Fluid structure interaction is one of the major studies in water hammer analysis : Tijsseling [1] carried out a very detailed review of transient phenomena in liquid- filled pipe systems. He dealt with water hammer, cavitations, structural dynamics and fluid-structure interaction (FSI). His main focus was on the history of FSI research in time domain. Wiggert and Tijsseling [2] attempted to succinctly summarize the essential mechanisms that cause FSI, and present relevant data that describe the phenomenon. In addition, the various numerical and analytical methods that have been developed to successfully predict FSI has been described. Heinsbroek [3] compared two different ways, MOC and finite element method for solving structural equation. Thos study showed that the FSI in pipeline system can adequately be investigated by application of MOC and FEM for hydraulic and structure of piping system, respectively. Tijsseling [4] illustrated the MOC approach for obtaining exact solution of FSI four equation model and it is generally referred in most of the latest papers for validating the results with benchmark examples.

Mechanical Engineering 08-12

Wiggert et al. [5] used a one dimensional wave formulation in both the liquid reaches and the piping structure resulting in five wave components and fourteen variables. MOC approach is introduced in this paper. Murat [6] studied finite element analyses with scheme of Finite difference method and method of characteristics to simulate aspects of water hammer. Le et al. [7] has independently solved the main FSI equation analytically. In both studies only junction coupling is considered but only for an unrestrained valve. In addition to the time domain analysis many researchers have studies the theoretical and experimental aspects in frequency domain. Ahmadi and Keramat [8] studied the analysis of water hammer with fluidstructure interaction (FSI) to investigate junction coupling effects. Junction coupling effects were studied in various types of discrete points, such as pumps, valves and branches. The emphasis was placed on an unrestrained pump and branch in the system and the associated relations were derived for modelling them. Proposed relations were considered as boundary conditions for the numerical modelling which was implemented using the finite element method for the structural equations and the method of characteristics for the hydraulic equations. Mahmoodi et al.[9] determined the position of loss of coolant accident in nuclear power plants using the transient vibration signal from a pipe rupture. A finite element formulation is implemented to include the effect of fluid-structure interaction. The coupled equations of fluid motion and pipe displacement are solved. Kermat et al. [10] investigated water hammer with FSI and viscoelasticity concerns pipes made of plastic, where FSI effects are more significant than in steel pipes, simply because they are more flexible (although with thicker walls, their modulus of elasticity is lower and their

Mechanical Engineering 08-12

Poissons ratio is higher). This makes FSI analysis necessary to reliably estimate the ultimate pipe stresses, elbow displacements and anchor forces, especially for designs with flexibly supported pipes. It was found herein that in the early moments of the transient event FSI is significant. Olson and Bathe [11] have formulated finite element procedure by establishing direct symmetric model based on fluid velocity potential. Sandberg [12] has proposed a symmetric finite element model for coupled acoustic vibration between fluid and structure and impact of eigen value shifting in complex system. Lee and Park [13] presented a methodology to transform coupled pipe dynamics equation into linear form about steady state value of fluid pressure and velocity. Spectral element model is used and compared with finite element model. Xiaodong [14] in his book described the fundamental procedure for solving FSI effect by use of virtual work method and establishing formulation of displacement-Velocity potential and displacement- velocity potential- Pressure method. Elghariani [15] in his thesis has discussed formulation of equations considering effect of FSI and solving it with MOC with friction coupling as a major source of investigation. Other works [16-17] in literature mentioned the finite element formulation for the fully coupled dynamic equations of motion to include the effect of fluidstructure interaction (FSI) and applied to a pipeline system used in nuclear reactors. 1.3 Objective and scope of work The main objective of the project is show the impact of each coupling individually in water hammer analysis for single phase fluid in simple reservoir-pipeline and valve system. Two finite element formulations are used to predict fluid structure interaction effect in water hammer. The results of two formulations have been analyzed to check the significant advantage of one over the other.

Mechanical Engineering 08-12

Chapter-2 MATHEMATICAL MODELLING


This chapter deals with mathematical modeling of physics present in Water hammer analysis. The formulation will give details of classical water hammer and change in equation if effect of FSI is considered. 2.1 Classical Water Hamme r Theory We consider following assumption for the study of classical water hammer theory. 1. One dimensional model is considered with average cross sectional velocity and pressure. 2. Friction is considered to be constant in our study. 3. The pipe is assumed to be fully filled in due course of study and any possibility of column separation or cavity formation is avoided. 4. There is no fluid other than water in pipe. 5. The density and other structural property as well as thermodynamic property remain constant. 6. The pipe is assumed to be straight, thin walled, linear elastic and of circular cross- section containing a weakly compressible fluid. 7. Additional effect of damping, friction and gravity are ignored in our study. As inertia forces are neglected in radial direction for both liquid and pipe wall therefore hoop stress is linearly related to pressure: (1)

Mechanical Engineering 08-12

where

-Hoop Stress, p- pressure, R- radius of pipe and e- thickness of pipe.

Classical water hammer includes only equation of momentum and equation of continuity: (2)

(3)

(4)

where above equation representing propagation of elastic wave through fluid in pipe. Speed. 2.2 Modification for classical wate r hammer equation

-Wave

This section describes the formulation of pipe dynamics equation represent the water hammer phenomenon. 2.2.1 Governing equation For Fluid Equation for conservation of momentum and conservation of mass for fluid considering effect of FSI as per Wiggert et al. [2] are: (5)

(6)

(7)

Mechanical Engineering 08-12

The friction effect is not considered in our study as its role in water hammer effect is very negligible as compared to other coupling. If equation 5 and 6 are differentiated with respect to axial distance and time respectively two wave equation are observed e ither in terms of pressure or velocity. These are elliptical in nature and are among best suited for FEM method. As valve closure is a function of velocity so equation of wave will be represented in terms of velocity: (8) 2.2.2 Governing equation for Structural part: The axial motion of pipe is described again in pipe dynamics equation expressed by Wiggert [2]
(9)

This is also wave equation elliptical in nature without consideration of friction and pressure head due to gravity. The 2D beam element with 2 degree of freedom at each node is modeled using pipe dynamics equation. 2.3Inital and Boundary condition: Appropriate initial and boundary condition are given to solve the benchmark problem of valve closure. 2.3.1 Initial condition 1. V(x, t=0, 0) = Vo 2. u (x, t= 0, 0) = 0 3. P(x, t= 0, 0) = (

Mechanical Engineering 08-12

There first and second derivative of above variable is taken as 0. 2.3.2 Boundary condition For boundary condition as we have to simulate Poisson and Junction coupling some equation expressing physical process in boundary phenomenon should be mentioned. As junction coupling is basically interaction between fluid and solid as in case of valve we have to express valve phenomenon with some equation. For a simple reservoir-pipe- valve system with structural velocity zero and pressure head constant the boundary condition for valve is

(10)

(11)

Ccalibration =1.475 2.4 Numerical Solution In this part important issues regarding problem formulation according to solution method will be discussed. Two methodologies have been discussed in subsequent section for analyzing water hammer effect considering effect of FSI is u- and u--P formulation. It is known from basic knowledge of FEM that this method is not good for solving fluid equation as it is incapable of travelling discontinuities but in our case it is of less important as in our reservoir-pipe and valve system chances of random change in pressure and velocity are negligible. The general equation stating finite element method which considers damping, stiffness and mass element is

Mechanical Engineering 08-12

10

(12) M- Mass Matrix, C- Damping Matrix, K- Stiffness Matrix. 2.4.1 Finite element method A number of Finite element formulations have been proposed to model fluid for FSI analysis. Among those formulation Displacement Velocity potential and Displacement- Velocity Potential-Pressure formulation are major breakthrough and are presented here. The comparison of results between them has been done in order to understand the difference between pros and cons of the formulation. These are the linear model and ignore material as well as geometric non linear ties. FSI interface motions are considered to be insignificant. 2.4.2 Displacement and Velocity potential formulation Based on Hamilton principle virtual work statement of problem is mentioned below: (12)

(13)

(14)

The total virtual force acting on structure due to flow of fluid and the resultant forces is given by (15) Where over dot means time derivative and prime means spatial derivative. E signifies to young modulus, L is length of pipe, As area of pipe and T1 and T2 are the axial forces acting on pipe at L=0 and L=L.
Mechanical Engineering 08-12 11

For one dimensional solid part the principle the Variational Indicator in defined as (16) Where = strain tensor. Cs = material stress strain tensor, = Surface traction vector, region. For compressible, inviscid fluid which experiences only small displacement and One dimensional fluid flow equation is: .

=density, u= displacement vector, = Body force vector, I = FSI

= FSI interface traction vector,

(17)

= velocity potential, k=fluid bulk modulus;

= fluid density, n = unit outward vector

Kinematic and dynamic matching conditions along the FSI interface for linking solid and fluid part: (18)

(19)

As per Hamilton principle all variation vanishes at t1 and t2. We solve equation15 and 16 by putting 17 and 18 (20)

(21)

Mechanical Engineering 08-12

12

Above two matrix are expressed as discretized form and in the form of matrix:

u1, 1

u2,2

Fig.2.4.2.1 Pipe element with two degree of freedom at each node

(22)

Solution for displacement Velocity potential formulation: [d]- displacement field, [F]= force vector,[di]and [do]=initial condition ,

[d]=

, [F]=

, [di]=

, [do]=

(23)

2.4.3 Displacement, velocity Potential and pressure formulation In this kind of formulation we need one hydrostatic pressure unknown P to replace one nodal velocity potential in last section. The discretized form of equation obtained from the similar procedure followed in above section in form of matrix is u2, 2, p2 u1, 1, p1 u1,1, p1

Fig2.4.3 Pipe element with 3 degree of freedom at three nodes

Mechanical Engineering 08-12

13

(24)

Solution for displacement velocity potential and pressure formulation: [d]- Displacement field, [F]= force vector,[di]and [do]=initial condition

(25)

The above two formulation have been used for study of fluid structure interaction. We have assumed an in viscid, ir-rotational compressible fluid with small motion in contact with elastic solids (pipe). In , we use velocity potential as state variable for fluid and displacement

for solid. In second type of formulation we have introduce pressure because this method offers symmetry formulation with reduced fluid degree of freedom, it can be easy incorporated in displacement based program. The coupling node in structural and fluid domain has same global co-ordinates; hence interface element has zero mass. The mass, stiffness and damping matrix has been established with force matrix representing basic finite element formulation. In next section we will calculate the matrix element. 2.4.4 Inte rpolation or Shape function As we are solving two case when pipe has two nodes with 2 degree of freedom at each node in first formulation discussed in 2.4.2 and in 2.4.3 we consider a pipe with three node with 3 degree

Mechanical Engineering 08-12

14

of freedom at each node. The finite element model is modeled by using displacement and fluid field. (26)

Where respectively. , and

represent the shape function for structural and fluid part are the displacement vectors. The selection o f

shape function depends upon degree of freedom. Isoparametric function are used for all degree of freedom related to fluid, structure and interface. In u- formulation we have 2 degree of freedom per node therefore we choose a linear shape function and this shape function can be used for structure and fluid for making calculation easy. 0xL (27)

and For second case of u--P formulation as degree of freedom per node are 3. To approximate this behavior with fluid finite element we have used quadratic shape function as pressure field is a solution to Helmholtz equation which is characterized by sinusoida l variation in space. 0xL (28) ,

Mechanical Engineering 08-12

15

2.4.5 Matrix element formulation in terms of shape function Using these shape function and expressing for mass, stiffness and coupling matrix have been derived from the above discussion in 2.4.2 and 2.4.3. They are defined as 1. Mass Matrix Element:

and

(29)

2. Stiffness matrix element , and , , (30)

3. Damping(coupling) matrix element

and (31)

4.

(a) Force Matrix for u- formulation

and

(32)

(b) Force matrix for u--P formulation

and

(33)

Mechanical Engineering 08-12

16

2.4.6 Solution Procedure As the simulation of water hammer including FSI effect is modeled for one dimension mathematical model we choose the numerical values of various parameters used by Serial no. 1 2 3 4 5 6 7 8 9 10 11 12 Parameter Length(m) Outer diameter(m) Thickness(m) Youngs Modulus Bulk Modulus Poissons Ratio Gravitational acc.(m/sec2) Fluid velocity(m/sec) Pressure behind valve(N/m2 ) Density of pipe(Kg/m3 ) Density of fluid(Kg/m3 ) Wave velocity(a*)(m/sec) Table2.1 Specification of piping system A MATLAB program was written following New mark scheme with alpha=0.5 and gamma=0.1677 to solve for all degree of freedom u, p and and obtain their variation with time. To examine the validation of proposed model the results were compared with Ahmadi and Keramat [8].There is no friction and valve closure time in 4 sec. Numerical Value 20 0.813 0.008 2.1*1011 N/m2 2.1*109 N/m2 0.3 9.8 1 0 7900 1000 1020

Mechanical Engineering 08-12

17

Chapter-3 RESULT AND DISCUSSION


3.1 Result and analysis for Displacement, Velocity Potential and Pressure Formulation The variation of various parameters like displacement amplitude, velocity, pressure and forces with time has been discussed 3.1.1 Displacement variation with time

Fig.3.1.1 shows the variation of displacement at different nodes with time. It can be observed from the graph that at node 3 which is actually place for junction coupling suffers maximum displacement accompanied by fluctuation as valve is closed gradually with time. We can also observe that as we move from valve end to reservoir end (from node 3 to node 1) the displacement amplitude decreases. The results are in good accordance with experimental result.[8]

0.008 0.007 0.006 0.005 0.004 0.003 0.002 0.001 0 -0.001 0 -0.002 -0.003

Displacement(m)

Displacement(m) Node 1
Displacement (m) Node 2 Displacement(m) Node 3

Time(sec)

Fig.3.1.1: Displacement (m) Vs Time (sec) fo r u--P method for 3 nodes

Mechanical Engineering 08-12

18

3.1.2 Velocity variation with Time Discussion The variation of velocity with time is very well expressed from Fig.3.1.2 as it is in great similarity with experimental result. The variation of velocity with length of pipe (different node) suggests that as the valve will be closed the velocity will increase randomly from valve end to reservoir end. The quantitative value gives an idea of large hydraulic force acting on structure.

40 30 20

Velocity(m/sec)

10
0 0 -10 -20 1 2 3 4

Velocity(m/sec) at node 1

Velocity(m/sec) at 5 node2
Velocity(m/sec) at node 3

-30
-40 -50

Time(sec)

Fig .3.1.2 Velocity (m/s) Vs Time (sec) for u--P method for 3 nodes

Mechanical Engineering 08-12

19

3.1.3Pressure variation with time Fig.3.1.3 shows the comparison between pressure amplitude between classical water hammer solution and FEM approach. There is been deviation in pressure amplitude measured from both method. The deviation is because of FSI effect considered in our approach. Classical water hammer result gives pressure amplitude of 300000 N/m2 where as our simulation result predicts peak pressure of 1000000 N/m2 . This result is corresponding to FSI effect caused by various coupling. (34) It also shows the random increase in pressure amplitude which is attributed to junction coupling effect, as the valve is closed gradually the pressure increases and velocity decreases(from previous result). It is also suggested from above graph that pressure amplitude falls with the length of pipe from valve to reservoir end. Peaks in the graph suggest the possible reason for structure deformation as the change in pipe amplitude is not gradual. The experimental results are in great accordance with it.

Mechanical Engineering 08-12

20

1000000 800000 600000 Pressure(N/m^2) 400000 200000 0

Pressure classical water hammer

Pressure( N/m^2)at node 1

Pressure( N/m^2)at node 2 Pressure( N/m^2)at node 3

-200000 -400000
-600000 -800000

Time(s)

Fig.3.1.3 Pressure (N/ m^2) Vs Time (sec) for u--p method for 3 nodes

3.1.4 Force on Solid (pipe) variation with time at each pipe node: Fig.3.1.4 presents a quantitative idea of force acting on structure (pipe) across the length of pipe. The force is maximum at node 3 which is accordance of theory because of valve closure there will be sudden rise in pressure and therefore force acting on valve due to fluid increases and as the pressure decreases along the length from valve to reservoir the force amplitude decrease. The values for force at different node are helpful in deciding the material of pipe and valve for a particular physical solution and are very important for design consideration.

Mechanical Engineering 08-12

21

Force on solid (N) node 1 Force on solid (N) node 2 1000 800 600 Force on solid (N) node 3

400 200
Force(N)

0
0 -200 -400 -600 -800 -1000 Time(s) 1 2 3 4 5

Fig 3.1.4 Force on structure (N) Vs Time (sec) for u--P method for 3 nodes

3.1.5 Force on fluid vs time variation at all nodes Fig.3.1.5 states that the force acting on fluid due to structure is very less as compared to fluid on structure this is because of highly rigid nature for pipe. For elastic pipes this force will be almost of same order of Force on structure (Pipe). This force gradually increases wave velocity of fluid and this gradual increase in fluid velocity results in gradual increase of pipe displacement.

Mechanical Engineering 08-12

22

Force on fluid(N) at node 1 300 200 Force on fluid(N) at node 2

100
0 Force(N) 0 -100 -200 1 2 3 4

Force on fluid(N) at node 3

-300 -400
-500

Time(sec)

Fig 3.1.5 Force on fluid (N) Vs Time (sec) for u--P method for 3 nodes

3.2 Result and analysis for Displacement and Velocity Potential formulation: Fig.3.2.1-3.2.4 Time history for displacement amplitude, velocity and forces are obtained after solving nonlinear equation; it is found in great similarity with published result [8]. The discussion has already been made in previous section. It is important to consider that forces exerted on solid and fluid are same as compared to previous result, the reason for this it is probably the nature of formulation where we added hydraulic force and force acting on valve to both force acting on structure due to fluid and vice versa.

Mechanical Engineering 08-12

23

0.008 0.006 0.004 0.002

Displacement(m)

0 0 -0.002 0.5 1 1.5 2 2.5 3 3.5 4 4.5

-0.004 -0.006 -0.008 -0.01 -0.012


Time(s)

Fig.3.2.1 Displacement (m) Vs Time (sec) for u- method

12 10

8
6 Velocity(m/sec) 4 2

0
-2 -4 -6 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5

-8

Time(s)

Fig.3.2.2 Velocity (m/sec) Vs Time (sec) for u- method

Mechanical Engineering 08-12

24

300 200

100
Force(N) 0 -100 -200 0 1 2 3 4 5

-300 -400
-500 Time(s)

Fig.3.2.3 Force on structure (N) Vs Time (sec) for u- method

400 300 200

Force(N)

100 0
0 1 2 3 4 5

-100
-200 -300

Time(s)

Fig.3.2.4: Force on flu id (N) Vs Time (sec) for u - method

Mechanical Engineering 08-12

25

3.3 Comparison of individual effect of coupling The comparison is done by taking account of three major cases 1. Considering effect of Poisson coupling 2. Considering effect of Junction coupling 3. Considering effect of both Junction and Poisson coupling Fig.3.3.1-3.3.4 presents a very important prospect of the result as we found that displacement, velocity and force variation is very less as compared to case where both the coupling are considered for formulation. This suggest that in case pipe material has high rigidity value the Poisson coupling is of less importance as its contribution to water hammer is small. If the pipe material is elastic in nature then the effect of Poisson coupling will surely increase and therefore become an important consideration for calculating the water hammer effect in elastic pipes. From Fig.3.3.1-3.3.4 It is found that when we consider only junction coupling without considering the impact of other two coupling the displacement, velocity and force variation do not differ a lot from the case where all coupling have been considered. This suggest that in case of valve closure boundary phenomenon with rigid pipe junction coupling play an important role in water hammer and is reason for pipe and valve rupture. It is also observed that pipe undergo compression at start of valve closure thereby by putting a huge pressure on liquid. It is also evident that both coupling when acted together produced enhanced amplitude of all variables displacement, velocity and forces than individual coupling. This is because when Poisson coupling come into picture it directly relates the pressure with stress of pipe and fluid velocity with change in pipe displacement, therefore both domain(fluid and solid) response to
Mechanical Engineering 08-12 26

change in each other dynamically which is not possible in case of only junction coupling alone. But our factor also depends upon on rigidity modulus of pipe material if it is high the role of Poisson coupling will be small and if its is small then its role be large. In our case Poisson coupling contributes to small deformation or change in Displacement amplitude.

0
-0.0005

0.5

1.5

2.5

-0.001
Displacement(m) -0.0015 -0.002 -0.0025

-0.003 Displacement(P)(m)
-0.0035 -0.004 Displacement(J)(m) Displacement(m) Time(s)

-0.0045

Fig.3.3.1 Variation of Displacement (m) Vs Time (sec) for u- method with consideration of all coupling

Mechanical Engineering 08-12

27

Velocity(P) (m/sec)

4
3.5 Velocity(m/sec)

Velocity(j)(m/sec) Velocity(m/sec)

3
2.5

2
1.5

1
0.5 0 0 0.5 1 1.5 Time(sec) 2 2.5

Fig 3.3.2: Variation of Velocity (m/s) Vs Time (sec) for u- method with consideration of all coupling

Force(S(p)(N)

20 0 -20 -40
Force (N) 0 0.5 1 1.5 2 2.5

Force(S)(j)(N) Force on Solid(N)

-60 -80
-100 -120 -140 -160

Time(s)

Fig.3.3.3 Variation of Force on structure (m) Vs Time (sec) for u- method with consideration of all coupling

Mechanical Engineering 08-12

28

Variation of Force on fluid due to different coupling


140 120

100
Force(N) 80 60 40

Force(L)(p) (N) Force(L)(j)(N) force on liquid( N)

20
0

0.5

1
Time(s)

1.5

2.5

Fig.3.3.4 Variation of Force on fluid (N) Vs Time (sec) fo r u- method with consideration of all coupling

3.4

Validating use of u--p formulation as compared to u-.

From Fig.3.4.1It is found that displacement amplitude variation with time calculated from both u- and u--P formulation is very close to each other. This numerical result validates the theoretical prediction of Olson and Sandberg [9, 10] in which they described that both the formulation results in almost same solution. It is also evident from our study that using u- -p formulation gives more accurate result as we have introduced one more degree of freedom called pressure and calculates its variation at FSI interface.

Mechanical Engineering 08-12

29

0.008
Displacementu-(m) 0.006 0.004 Displacement(m)u--p Node 1

0.002 Displacement(m)
0 0 -0.002 1 2 3 4 5

-0.004
-0.006 -0.008

-0.01
-0.012

Time(s)

Fig .3.4.1: Co mparison of Displacement (m) Vs Time (sec) for u- and u--P method with consideration of all coupling

Mechanical Engineering 08-12

30

Chapter 4 CONCLUSIONS
4.1 Summary of Result This report attempts to get a coupled solution for water hammer problem in straight pipe with valve closure at an end and reservoir at other end. Time histories for displacement, velocity, pressure and forces were compared along the length of pipe. We also tried to show effect of Junction and Poisson coupling as compared to classical water hammer. Some of the important conclusions drawn out of study are given below. [1] The variation of Displacement suggests that it is maximum near valve and slowly reduces towards reservoir. [2] The pressure peaks are found to be as large as 1000000 N/m2 , this data reveals the destructive nature of this phenomenon as the pressure has raised 10 times more than atmospheric pressure which is great concern for safe design. [3] The force acting on structure is more than the forces acting on fluid. The values obtained from our study provide a great idea to set up designer about valve and pipe material selection so that it can bear such instant high loads. [4] For non FSI simulation the predicted value are as much as 45% less for velocity as compared to simulation carried out by considering the effect of FSI which suggests the importance of our study and its use for safe design.

Mechanical Engineering 08-12

31

[5] The Junction coupling will play a major role in water hammer than poison coupling. In our result Junction coupling is contributing to water hammer effect by more than 70% [6] The u- and u--p formulation results do not vary much as suggested by OLSON, this theoretical prediction has been validated .The result also shows that u--p formulation predicted the variation of displacement and pressure in great accordance to practical data. [7] This kind of formulation for predicting water hammer is validated as compared to other method like MOC, mixed MOC-FEM and other FEM techniques. 4.2 Future Scope of Work The method discussed over here provides a great deal of understanding of such difficult non linear formulation by FEM. It is advised to work with interpolation function and follow nonlinear integration scheme to predict result more accurately. It is also suggested to take this simulation to next higher level and simulate it for two dimensions or with mix phase of fluid to give more practical touch to problem. The variety of pipe with different stiffness and rigidity can be tested and analyzed so that suitable selection of pipe and valve can be made.

Mechanical Engineering 08-12

32

REFERENCES
1. Tijsseling, A.S., Fluid-structure interaction in liquid- filled pipe systems: a review. J. Fluids and Struct. 10, 109146, 1996. 2. Wiggert, D.C , Tijsseling A.S. , Fluid transients and fluid-structure interaction in flexible liquid- filled piping, ASME Journal of Appl. Mech. Rev, Vol. 54, pp 455481,2001. 3. Heinsbroek A., Fluid Structure interaction in non-rigid pipeline system, Journal of Nuclear Engineering Design, Vol. 172, 123-135, 1996. 4. A.S. Tijsseling, Exact solution of linear hyperbolic four-equation system in axial liquidpipe vibration, J. Fluids and Structures, 18, 179196, 2003. 5. Wiggert, D.C, Hatfield F.J., Struckenbruck S., Analysis of liquid and structural transients in piping by method of characteristics, ASME Journal of Fluid Engineering, Vol. 109, pp 161-165, 1987 6. Murat,E., Simulation of passive Water Hammer in Pipes, Maters Thesis MIT, 35-36, 1996. 7. Li, Q.S., Yang, K., Zhang, L., Zhang, N., Frequency domain analysis of fluid structure interaction in liquid- filled pipe systems by transfer matrix method. Int. J. Mech. Sci. Vol.44, pp 20672087, 2002. 8. Ahmadi,A. and Keramatn, A., Investigation of fluidstructure interaction with various types of junction coupling, Journal of Fluids and Structures, Vol.26, pp.11231141,2010. pp

Mechanical Engineering 08-12

33

9. Mahmoodi.R, Shahriari.M, Zolfaghari.A, Minuchehr, An advanced method of determination of loss of coolant accident in nuclear power plant, Journal of Nuclear Engineering and Design, Vol. 241, pp 2013-2019, 2011. 10. Keramat A, Tijsseling A.S., Hou Q., Ahmadi Ahmad, Fluidstructure interaction with pipe-wall viscoelasticity during water hammer, Journal of Fluids and Structures, Vol. 28, pp 434-455, 2012. 11. Olson and Bathe, Analysis of fluid structure interactions. A direct symmetric coupled formulation based on the fluid velocity potential, Journal of computers and structure. Vol. (21), pp 21-32, 1985 12. Sandberg.G, A symmetric finite element formulation of acoustic fluid structure interaction analysis, Journal of sound and vibration, Vol. 123(3), pp 507-515. 1988. 13. Lee.U and Park.J, Spectral element modeling and analysis of a pipeline conveying internal unsteady fluid, Journal of Fluids and Structure, Vol.22, pp 273-292. 2006 14. Xiaodong,W. Fundamental of Fluid Structure Interactions: Analytical and computational approach. Elsevier Publication, Oxford, UK 313-328. 2008. 15. Elghariani,S.K., Transient analysis of fluidstructure interaction in straight pipe Master Thesis Universiti Teknologi Malaysia, pp 17-34, 2007. 16. Jayaraj, K., Ganesan, N., Padmanabham, C., Model reduction for parametric instability analysis of shells conveying fluid Journal of Sound and Vibration, Vol.262, pp 633 649, 2003. 17. Sreejith, B., Jayaraj,K., Ganesan, N., Padmanabhan, C., Chellapandi P. and Selvaraj,P. Finite element analysis of fluidstructure interaction in pipeline systems, Nuclear Engineering and Design, Vol. 227, pp.313322, 2004.

Mechanical Engineering 08-12

34

18. Ghidaoui, M.S., Zhao, M., McInnis, D.A. and Axworthy, D.A., A review of water hammer theory and practice, Applied Mechanics Reviews, Vol.58, pp.49-76, 2005.

Mechanical Engineering 08-12

35

APPENDIX
NUMERICAL VALUE FOR MATRIX and MATLAB PROGRAM
(a) Numerical values for mass, stiffness and damping matrices have been calculated and presented in matrix format The formulation of matrix for both u- .

The formulation for matrix for and u- -P are:

Mechanical Engineering 08-12

36

(b) MATLAB Program Matlab Code of Integration of N*N matrix using a Newmark Scheme: The code predict time dependent values for displacement, velocity, pressure and forces. The graph is obtained with the help of excel. This code provide flexibility to user to N*N dimensional matrix depending upon nature of problem. alpha=0.5 and gamma =0.1677 b.1 For u- formulation (2*2 matrix)
function [u,ud,udd] = Water Hammer( ) format long n=input('no. of time increments required n='); dt=input('time interval d t ='); t(1)=0; for i=1:n t(i+1)=t(i)+dt; Cv=1.78*(i*dt)^0.239; end

Mechanical Engineering 08-12

37

disp('input the stiffness and mass matrices'); k=input('stiffness matrix='); m=input('mass matrix='); c=input('damping matrix c='); disp('stiffness matrix k='); disp(k); disp('mass matrix m='); disp(m); disp('damping matrix c='); disp(c); u=[0;0;1;1]; ud=[0;0;0;0]; udd=[0;0;0;0]; alpha=input('integration constantl alpha ='); gamma=input('integration constant gamma='); a0=1/(alpha*dt^2); a1=gamma/(alpha*dt); a2=1/(alpha*dt); a3=(1/(2*alpha))-1; a4=(gamma/alpha)-1; a5=((gamma/alpha)-2)*(dt/2); a6=dt*(1-gamma); a7=gamma*dt; keff=k+a0*m+a1*c; for i=1:n f=[-3.34-35.76*Cv; 5.31+25.26*Cv; -2-35.76*Cv; 3+25.26*Cv]; disp('force matrix f='); disp(f); s = m*(a0*u+a2*ud+a3*udd)+c*(a1*u+a4*ud+a5*udd); feff=f+s; disp('feffective ::') disp(feff) temp1=u; u=inv(keff)*feff; temp=udd; udd=a0*(u-temp1)-a2*(ud)-a3*(udd); ud=ud+a6*temp+a7*udd; disp('upadated value of u:') disp(u) % plot(t(i)+dt,u(i)); disp('upadated value of ud:') disp(ud) disp('upadated value of udd:') disp(udd) end end Input:

stiffness matrix=[3950000000 -3950000000 0 0;-3950000000 3950000000 0 0;0 0 6666.66 3333.33;0 0 3333.33 6666.66]; mass matrix=[6393.160 3196.630 0 0;3196.630 6393.260 0 0; 0 0 3.16*10^-3 1.58*10^-3; 0 0 1.58*10^-3 3.16*10^-3]; damping matrix c=[30 -30 0 0;-30 30 0 0 ;0 0 30 -30;0 0 -30 30];

Mechanical Engineering 08-12

38

b.2 For u--P formulation


function [u,ud,udd] = Water Hammer( ) format long n=input('no. of time increments required n='); dt=input('time interval d t ='); t(1)=0; for i=1:n t(i+1)=t(i)+dt; Cv=1.78*(i*dt)^0.239; end disp('input the stiffness and mass matrices'); k=input('stiffness matrix='); m=input('mass matrix='); c=input('damping matrix c='); disp('stiffness matrix k='); disp(k); disp('mass matrix m='); disp(m); disp('damping matrix c='); disp(c); u=[0;0;0;1;1;1;0;0;0]; ud=[0;0;0;0;0;0;0;0;0]; udd=[0;0;0;0;0;0;0;0;0]; alpha=input('integration constantl alpha ='); gamma=input('integration constant gamma='); a0=1/(alpha*dt^2); a1=gamma/(alpha*dt);a2=1/(alpha*dt); a3=(1/(2*alpha))-1;a4=(gamma/alpha)-1;a5=((gamma/alpha)-2)*(dt/2); a6=dt*(1-gamma);a7=gamma*dt; keff=k+a0*m+a1*c; for i=1:n f=[94.7 ;104.76; -199.95;190; 210; -400; 1.404*Cv; 51.63*Cv;-64.10*Cv]; disp('force matrix f='); disp(f); s = m*(a0*u+a2*ud+a3*udd)+c*(a1*u+a4*ud+a5*udd); feff=f+s; disp('feffective ::') disp(feff) temp1=u; u=inv(keff)*feff; temp=udd; udd=a0*(u-temp1)-a2*(ud)-a3*(udd); ud=ud+a6*temp+a7*udd; disp('upadated value of u:') disp(u) % plot(t(i)+dt,u(i)); disp('upadated value of ud:') disp(ud) disp('upadated value of udd:') disp(udd) end end

Mechanical Engineering 08-12

39

Start

Flowchart

input

exit

no

yes

increment value of t and CV

calculate const ao-a7, calculate keff,

exit

no i=i+1

calculate t(i+1)=t(i)+dt

display f, f eff,u,udd, ud

stop Mechanical Engineering 08-12 40

You might also like