Simulation of Multi-batch Driven Pipelines
|
|
- Dana Walton
- 6 years ago
- Views:
Transcription
1 Simulation of Multi-batch Driven Pipelines DRAGO MATKO, SAŠO BLAŽIČ University of Ljubljana Faculty of Electrical Engineering Tržaška 25, SI-000 Ljubljana SLOVENIA GERHARD GEIGER University of Applied Sciences Gelsenkirchen FB Elektrotechnik Neidenburger Str. 0, D Gelsenkirchen GERMANY Abstract: The problem of modelling and simulating pipelines that are used for transporting different fluids is addressed in the paper. The problem is solved by including fluid density in the model beside pressure and velocity of the medium. First, the system of nonlinear partial differential equations is derived. Then, the obtained model is linearised and transformed into the transfer function form with three inputs and three outputs. Admittance form of model description is presented in the paper. Since the transfer function is transcendent, it cannot be simulated using classical tools. Rational transfer function approximation of the model were used and validated on the real industrial pipeline. It was also compared to the model that does not take the changes in fluid density into account. The latter model cannot cope with batch changes whereas the proposed one can. Key Words: Pipeline, Fluid dynamics simulation, Transcendent transfer function, Lumped parameter model Introduction Real time transient model (RTTM) based leak monitoring systems require a sophisticated mathematical model of the flow in pipelines. The so called water hammer equations are relatively simple mathematical models assuming isentropic flow; they are obtained using the principles of mass and momentum conservation []. However, in the case when different fluids are transported through the same pipeline, the above model is not adequate. The water hammer equation can easily be extended, and this enables simplified description of multi-product-flows with multiple products or batches being transported at the same time in one pipeline. Up to now, there is no analytical solution for this nonlinear, partial differential equation system available. Instead, numerical solution techniques like the method of characteristics can be used [2]. Another possibility to solve the problem is to use linearisation and Laplace transformation techniques in order to get a frequency domain description [3]. This leads to a simplified pipeline model with lumped parameters. We hereby get some advantages: the classical system theory for Multi-Input Multi-Output (MIMO) systems can be used, e.g. for controller design and system identification. The resulting algorithms are less time-consuming and hence better suited for critical real time applications. Additionally, the analysis of fluid transients caused by leaks is much easier. In Section 2 the nonlinear model of a pipeline is given that takes into account multiple fluids being transported. It is linearised, and in Section 3 a simplified model with lumped parameters is referred to. The obtained models is compared to the one derived in [3] that assumes constant density of the fluid. The results are validated on the real industrial pipeline in Section 4. At the end some conclusions are given. 2 Mathematical model The nonlinear model of a pipeline that takes into account multiple fluids being transportedis as follows [4] dρ dt + ρ v = 0 () dv dt + p λv v + g sinα + ρ 2D = 0 (2) dp dt a2dρ dt = 0 (3) where D is the diameter of the pipeline, ρ(x) and v(x) the density and the velocity of the fluid respectively, λ(v) the friction coefficient, a the (isentropic) speed of sound in the fluid and p(x) the pressure. The quantity g sinα is the x-component of the standard gravity vector g. By replacing total derivatives with partial ones, multiplying Eq. (2) with ρ, and by inserting Eq. () ISSN: Page 53 ISBN:
2 v[m/s] ρ[kg/m 3 ] p[bar].5 2 x x[m] x[m] x[m] Figure : Steady state solution of the pipeline: pressure p, velocity v, and density ρ into Eq. (3) we get ρ t + v ρ + ρ v = 0 (4) ρ v v + ρv t + p ρλv v + ρg sinα + 2D = 0 (5) p t + v p + a2 ρ v = 0 (6) In order to simplify this system of partial differential equations, first the steady state will be evaluated. It is obtained by setting to zero all partial derivatives with respect to time ( t := 0) which yields the following set of ordinary differential equations (with respect to x): v d ρ + ρd v = 0 (7) ρ v d v + d p ρλ v v + ρg sin α + 2D = 0 (8) v d p d v + a2 ρ = 0 (9) This set of equations cannot be solved analytically. It was solved by MATLAB - SIMULINK for the pipeline used for the verification of the model (see Section 4). The results for p, v and ρ are shown in Fig.. The convective terms ρv v p and v in Eqs. (5) and (6), respectively, are small compared to the derivatives with respect to time ρ v p t and t, respectively, and will be neglected. Next, Eqs. (5) and (6) will be linearised around steady-state solution. In order to do this, new variables will be introduced: ṽ(x, t) = v(x, t) v(x) (0) p(x, t) = p(x, t) p(x) () ρ(x, t) = ρ(x, t) ρ(x) (2) After linearising Eqs. (5) and (6), and taking into account Eqs. (8) and (9), the following is obtained: ρ ṽ t + p ρλ v + ρg sinα + D ṽ + λ v v 2D ρ = 0 (3) p ṽ d v + a2 ρ + a2 ρ t = 0 (4) The last term in Eq. (4) can be neglected due to the fact that d v is very small (see Fig. ). Using notations that are common in the analysis of electrical transmission lines: and by denoting we get T L = ρ (5) R = ρλ v D (6) C = a 2 ρ (7) = g sinα + λ v v 2D L ṽ + Rṽ + T ρ = p t C p t = ṽ (8) (9) (20) Since the fluid in the pipeline is almost incompressible, the dynamics of ρ in Eq. (4) are relatively slow, compared to those of v and p in Eqs. (5) and (6). As a consequence, the velocity profile along the pipeline is nearly constant within each time instant. Due to a special reason that will be explained later, the derivative v will not be neglected in Eq. (4). Rather, it will be approximated with a very small constant v x. Linearising (4) and considering (2) then yields: + v(x) + ρ(x, t) t ρ(x, t) + ṽ(x, t) d ρ(x) + + ρ(x, t) v x = 0 (2) Note that d ρ is very small (see Fig. ). It could be neglected in Eq. (2), but it will not be at this point since its influence is similar to the influence of d v (see Eq. 7). The first step in obtaining analytical solution of the system is to solve Eq. (2) by applying Laplace transformation on it. It is assumed that it is permissible to interchange the order of differentiation with respect to x and the taking of the Laplace transform. ISSN: Page 54 ISBN:
3 The consequence is that the first order ordinary differential equation is obtained [5] s ρ(x, s) ρ(x,0) + ṽ(x, s) d ρ(x) + d ρ(x, s) + v(x) + v x ρ(x, s) = 0 (22) It is assumed here that the system rests at time t = 0, i.e. ρ(x,0) = 0, p(x,0) = 0, ṽ(x,0) = 0. Equation (22) can then be transformed to d ρ(x, s) + s + v x v(x) ρ(x, s) = d ρ(x) ṽ(x, s) v(x) (23) The solution of the homogenous part of Eq. (23) would be with ρ(x, s) = ρ 0 (s)e mx (24) ρ 0 (s) = ρ(0, s) (25) m = s + v x v(x) (26) if m (or v) did not depend on x. Since the latter assumption is violated and Eq. (23) is not homogenous, the candidate for the solution is ρ(x, s) = ρ (x, s)e mx (27) Inserting Eq. (27) into Eq. (23), and taking into account yields dm = s + v x d v v 2 = m v d v dρ (x, s) e mx + ρ (x, s)e mx ( m x dm ) + mρ (x, s)e mx = d ρ(x) ṽ(x, s) v(x) dρ (x, s) = ρ (x, s)m x v d v emxd ρ(x) (28) ṽ(x, s) v(x) (29) Since d v are very small, it follows from Eq. (29) that ρ can be regarded as independent of x (at least on the interval of interest x [0, L p ]). Therefore, the solution (24) will be treated in the paper. Similarly, dependence of physical parameters of the pipeline (L, R, C, T, and m) on x shall be neglected in the rest of the paper since they are only functions of v and ρ. and d ρ In the following, only the deviation model will be considered. To simplify the notation, the tildes will be omitted in the equations. The varaibles p, v, and ρ will stand for the deviations of the respective variables from the stationary values. The next step in the derivation of a simple model of the pipeline is the analytical solution of linear Eqs. (9) and (20). Performing the Laplace transformation on them yields dp(x, s) (Ls + R)V (x, s) = Tρ(x, s) (30) dv (x, s) Cs P(x, s) = (3) where stationary initial conditions v and p were assumed, i.e. v(x,0) = 0, p(x,0) = 0. By differentiating Eqs. (30) and (3) with respect to x dv (x, s) (Ls + R) dp(x, s) Cs = d2 P(x, s) dρ(x, s) 2 T (32) = d2 V (x, s) 2 (33) the following equations are obtained using Eqns. (30) and (3) (Ls + R) Cs P(x, s) = d2 P(x, s) dρ(x, s) 2 + T (34) (Ls + R) Cs V (x, s) = d2 V (x, s) 2 CsTρ (35) Taking into account Eq.(24) we obtain: (Ls + R) Cs P(x, s) = d2 P(x, s) 2 (Ls + R) Cs V (x, s) = d2 V (x, s) 2 Tρ 0 (s)me mx CsTρ 0 (s)e mx which are known as wave equations. Their solutions are P(x, s) = C (s)e nx + C 2 (s)e nx Tρ 0(s)m n 2 m 2 e mx (36) V (x, s) = C 3 (s)e nx + C 4 (s)e nx CsTρ 0(s) n 2 m 2 e mx where (37) n 2 = (Ls + R) Cs (38) ISSN: Page 55 ISBN:
4 The four expressions, C (s), C 2 (s), C 3 (s), and C 4 (s) are not completely independent. By introducing (37) into (30) and by differentiating (36) with respect to x, we obtain (Ls + R)[C 3 (s)e nx + C 4 (s)e nx CsTρ 0 (s) n 2 m 2 e mx ] = nc (s)e nx nc 2 (s)e nx Tρ 0 (s)m 2 n 2 m 2 e mx Tρ 0 (s)e mx (39) From Eq. (39), the following relations between C (s), C 3 (s), and C 2 (s), C 4 (s), respectively, are obtained C (s) C 3 (s) = Ls + R Ls + R = = (40) n Cs C 2 (s) C 4 (s) = Ls + R = n Ls + R Cs = (4) since [ ] (Ls + R)Cs Tρ 0 (s)e mx n 2 m 2 m2 n 2 m 2 0 (42) The term is called the characteristic impedance. The complete solution is obtained applying the boundary conditions. First, using the boundary conditions for x = 0 P(0, s) = P 0 (s) = C (s) + C 2 (s) Tρ 0(s)m n 2 m 2 = (C 3 (s) C 4 (s)) Tρ 0(s)m n 2 m 2 (43) V (0, s) = V 0 (s) = C 3 (s) + C 4 (s) CsTρ 0(s) n 2 m 2 (44) the coefficients C 3 (s) and C 4 (s) are obtained in the following form C 3 (s) = 2 V 0(s) + P 0 (s) + Tρ 0(s) 2 2 n m (45) C 4 (s) = 2 V 0(s) P 0 (s) + Tρ 0(s) 2 2 n + m (46) Next, the boundary conditions for x = L p are used in Eqs. (36) and (37) P(L p, s) = P L (s) = [C 3 (s)e nlp (47) C 4 (s)e nlp ] Tρ 0(s)m n 2 m 2 e mlp V (L p, s) = V L (s) = C 3 (s)e nlp + (48) + C 4 (s)e nlp CsTρ 0(s) n 2 m 2 e mlp By introducing (45, 46) into (47, 48) the inverse chain representation of the pipeline is obtained P L (s) = V 0 (s)sinh(nl p ) + P 0 (s)cosh(nl p ) + + Tmρ [ ] 0(s) n 2 m 2 cosh(nl p ) n m sinh(nl p) e mlp (49) V L (s) = V 0 (s)cosh(nl p ) P 0 (s)sinh(nl p ) + + Tmρ [ ] 0(s) n Z k (n 2 m 2 ) m cosh(nl p) sinh(nl p ) n m e mlp (50) Expressions (49) and (50) can be simplified by the evaluation of n m using Eqs. (26), (38), (5), (6), and (7): ) ( ρs + n ρλ v (Ls + R)Cs m = (s + v x )/ v = D a 2 ρ s (s + v x )/ v = v s 2 + λ v D s a (s + v x ) 2 (5) Note that n m since the sound speed is much bigger than the fluid speed (a v). Consequently, the terms with n m in square brackets in Eqs. (49) and (50) can be neglected. By taking into account m n 2 m 2 = ( ( ) m n 2 m m) v = s + v x and (52) n Z k = Cs (53) Eqs. (49) and (50) take the following form P L (s) = V 0 (s)sinh(nl p ) + P 0 (s)cosh(nl p ) T vρ ( ) 0(s) cosh(nl p ) e mlp (54) s + v x V L (s) = V 0 (s)cosh(nl p ) P 0 (s)sinh(nl p ) T vρ 0(s) Z k (s + v x ) ( sinh(nl p)) (55) The constant v x is very small but positive (what can be seen from Fig. ). Consequently, s+ v x is a stable transfer function. In the time scope of interest this transfer function is equivalent to an integrator, and will be replaced by one in the equations of the model. This is the reason why v x was not neglected ISSN: Page 56 ISBN:
5 in the early phase of the model derivation. Similarly, the term e mlp in Eq. (54) becomes e mlp = e s+ vx L v p e Lp v s = e τs (56) and can be interpreted as a pure delay system (τ is the time needed for the fluid to reach the outlet from the inlet of the pipeline transport delay). Apart from Eqs. (54) and (55), an additional equation is needed for the complete description of the system. This equation defines the density at the pipeline outlet and can be obtained by setting x to L p in Eq. (24): ρ(l p, s) = ρ L (s) = ρ 0 (s)e mlp = ρ 0 (s)e τs (57) Finally, we arrive to the noncausal (the transfer function matrix goes to as s ) representation of the pipeline: P L (s) V L (s) = [A P 0 (s) A 2 ] V 0 (s) (58) ρ L (s) ρ 0 (s) where A = A 2 = cosh(nl p ) sinh(nl p ) sinh(nl p ) cosh(nl p ) 0 0 Tv 0 s (cosh(nl p ) e τs ) Tv 0 Z k s sinh(nl p) e τs The linearised model of the pipeline (58) can be written in different forms which differ from each other with respect to the model inputs (independent quantities) and outputs (dependent quantities), however here only the admittance representation (inputs P 0, P L, ρ 0, outputs V 0, V L, ρ L ) will be given V 0 P 0 V L = [A 3 A 4 ] P L (59) ρ L ρ 0 where Z K coth(nl p ) sinh(nl p) A 3 = sinh(nl p) coth(nl p ) (60) 0 0 ( Tv 0 Z k s coth(nl p ) sinh(nl e τs) ( p) A 4 = Tv 0 Z k s sinh(nl coth(nl p) p)e τs) (6) e τs It should be noted that the form (59) represents causal models, whilst Eqn. (58) represent a noncausal model. Therefore, the latter cannot be realised by means of simulation. This completes the derivation of the transfer functions of the linearised pipeline. The resulting transfer functions are transcendent. In the next section, their approximations by rational transfer functions will be given. 3 Simplified pipeline model with lumped parameters In this section, the rational transfer functions of the pipeline will be refered to since such transfer functions are much easier to simulate, e.g. when designing real time transient model based leak monitoring systems. There are two different transcendent functions (e τs is the third one, but this one will be left as it is, since it can be simulated directly): coth(nl p ), and sinh(nl p).their approximations with rational transfer functions (obtained by expanding the transcendent transfer functions into a Taylor series) can be found in [3]. 4 Validation of the simplified model The models were validated on a real pipeline with the following data: length of the pipeline L p = 9854 m, velocity of sound a = 059 m/s, friction coefficient λ = 0.058, gravity constant g = 9.8 m/s 2, diameter D = m, and inclination α = rad. Simulations of the plant were performed on the lumped parameter model in admittance form. The most significant change in the fluid density ρ occurs during a batch change, so the data was recorded in the operation phase where the batch was changed twice. Due to operational reasons the pipeline must be stopped before and after a batch change. The stationary operation between 7000 s and 4000 s was chosen as the operating point. Figure 2 shows the time courses of the (measured) density ρ at the inlet (dash-dot line) and outlet (measured density solid line; model response dotted line) of the pipeline, respectively. The batch changes (000 s to 5300 s and 4800 s to 9000 s) can be seen clearly. The model response corresponds perfectly to the measured data for the second batch change, where the current operating conditions meet the chosen operating point. However, during the first change some missagreement can be noticed between the measured and the simulated density at the pipeline outlet. The difference is expectable since the system is not situated in the operating point where the linearised model was ISSN: Page 57 ISBN:
6 th WSEAS International Conference on FLUID MECHANICS (FLUIDS'08) Acapulco, Mexico, January 25-27, Measured data 850 ρ[kg/m 3 ] without ρ with ρ 650 v inlet [m/s] Time [s] Figure 2: The time courses of the density ρ measured ρ at the inlet (dash-dot line), measured ρ at the outlet (solid line) and simulated ρ (dotted line) obtained. Next, the velocity part of the lumped parameters model will be validated and compared to the model presented in [3]. The latter was obtained based on assumption that the density of the fluid being transported is constant all the time. The model proposed here takes the fluid changes into account. The comparison will allow us to estimate the benefit of the extended model. Simulation results of the both models will be compared to the real plant data. In Figure 3 three time courses of the fluid velocity at the inlet of the pipeline are shown: the measured one and responses of two models (with and without the consideration of the density ρ). It can be seen that the proposed model which takes into consideration the density ρ can cope with the changeable operating conditions whereas the model from [3] cannot. Note that the latter model performs very well between 7000 s and 4000 s where current density is approximately equal to the constant density of the simple model. But after the batch change the simulation results are not very satisfactory. In Figure 4 a detail of the Figure 3 is depicted, where the transient phase can be seen in detail. Good coincidence between the measured data and the proposed model response can be established. 5 Conclusion The results explained in the previous sections show much better results with multibatch pipelines as previous methods References: [] V. L. Streeter, E. B. Wylie. Fluid Transients. McGraw-Hill, New York, 978. x Simulated response with ρ without ρ Time [s] Figure 3: The time courses of the fluid velocity at the inlet of the pipeline measured velocity, model response without and with the consideration of the density ρ v inlet [m/s] Measured data without ρ Time [s] with ρ Figure 4: A detail of Figure 3 [2] A. H. Shapiro. The Dynamics and Thermodynamics of Compressible Fluid Flow, volume. The Ronald Press Company, 953. [3] D. Matko, G. Geiger, W. Gregoritza. Pipeline simulation techniques. Mathematics and Computers in Simulation, 52:2 230, [4] Blažič S., D.Matko, I. Škrjanc G.Geiger. Mathematical modelling of pipelines transporting different fluids. Math. comput. model. dyn. syst., :95 09, [5] D. H. Himmelblau, K. B. Bischoff. Process Analysis and Simulation, Deterministic Systems. John Wiley & Sons, New York, 968. x 0 4 x 0 4 ISSN: Page 58 ISBN:
Thermodynamics ENGR360-MEP112 LECTURE 7
Thermodynamics ENGR360-MEP11 LECTURE 7 Thermodynamics ENGR360/MEP11 Objectives: 1. Conservation of mass principle.. Conservation of energy principle applied to control volumes (first law of thermodynamics).
More informationCHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS NOOR ALIZA AHMAD
CHAPTER 3 BASIC EQUATIONS IN FLUID MECHANICS 1 INTRODUCTION Flow often referred as an ideal fluid. We presume that such a fluid has no viscosity. However, this is an idealized situation that does not exist.
More informationDevelopment of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More information2.25 Advanced Fluid Mechanics
MIT Department of Mechanical Engineering.5 Advanced Fluid Mechanics Problem 4.05 This problem is from Advanced Fluid Mechanics Problems by A.H. Shapiro and A.A. Sonin Consider the frictionless, steady
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationTheoretical Models of Chemical Processes
Theoretical Models of Chemical Processes Dr. M. A. A. Shoukat Choudhury 1 Rationale for Dynamic Models 1. Improve understanding of the process 2. Train Plant operating personnel 3. Develop control strategy
More informationWhere does Bernoulli's Equation come from?
Where does Bernoulli's Equation come from? Introduction By now, you have seen the following equation many times, using it to solve simple fluid problems. P ρ + v + gz = constant (along a streamline) This
More informationTransient Phenomena in Liquid/Gas Flow in Pipelines
Proceedings of the International Conference on Heat Transfer and Fluid Flow Prague, Czech Republic, August 11-12, 214 Paper No. 71 Transient Phenomena in Liquid/Gas Flow in Pipelines Zohra Ouchiha, S.
More informationInvestigation of a nonlinear dynamic hydraulic system model through the energy analysis approach
Journal of Mechanical Science and Technology 3 (009) 973~979 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-9x DOI.07/s6-009-081- Investigation of a nonlinear dynamic hydraulic
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationChapter Four fluid flow mass, energy, Bernoulli and momentum
4-1Conservation of Mass Principle Consider a control volume of arbitrary shape, as shown in Fig (4-1). Figure (4-1): the differential control volume and differential control volume (Total mass entering
More informationIntroduction to Physical Acoustics
Introduction to Physical Acoustics Class webpage CMSC 828D: Algorithms and systems for capture and playback of spatial audio. www.umiacs.umd.edu/~ramani/cmsc828d_audio Send me a test email message with
More informationFLOWS IN GAS PIPELINES
FLOWS IN GAS PIPELINES 1. Introduction This problem was brought to the 1989Mathematics in Industry Study Group by Mr R. Calvert from The Pipeline Authority (TPA). TPA owns and operates a natural gas pipeline
More informationWebster s horn model on Bernoulli flow
Webster s horn model on Bernoulli flow Aalto University, Dept. Mathematics and Systems Analysis January 5th, 2018 Incompressible, steady Bernoulli principle Consider a straight tube Ω R 3 havin circular
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More information( ρ Adx) = ρ AV t (1) ρ A V ρ + dx A + dxv + dx x x x Neglecting higher order derivatives and using the definition of material derivative, Eq. (1) can
16 th Australasian Fluid Mechanics Conference Crown Plaza, Gold Coast, Australia -7 December 007 Numerical Simulation of Pig Motion through Gas Pipelines S.M. Hosseinalipour 1, A. Zarif Khalili 1, and
More informationTotal energy in volume
General Heat Transfer Equations (Set #3) ChE 1B Fundamental Energy Postulate time rate of change of internal +kinetic energy = rate of heat transfer + surface work transfer (viscous & other deformations)
More informationBERNOULLI EQUATION. The motion of a fluid is usually extremely complex.
BERNOULLI EQUATION The motion of a fluid is usually extremely complex. The study of a fluid at rest, or in relative equilibrium, was simplified by the absence of shear stress, but when a fluid flows over
More informationFundamentals of Acoustics
Fundamentals of Acoustics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy of Sciences
More informationChE 6303 Advanced Process Control
ChE 6303 Advanced Process Control Teacher: Dr. M. A. A. Shoukat Choudhury, Email: shoukat@buet.ac.bd Syllabus: 1. SISO control systems: Review of the concepts of process dynamics and control, process models,
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationDevelopment of Dynamic Models. Chapter 2. Illustrative Example: A Blending Process
Development of Dynamic Models Illustrative Example: A Blending Process An unsteady-state mass balance for the blending system: rate of accumulation rate of rate of = of mass in the tank mass in mass out
More informationNavier-Stokes Equation: Principle of Conservation of Momentum
Navier-tokes Equation: Principle of Conservation of Momentum R. hankar ubramanian Department of Chemical and Biomolecular Engineering Clarkson University Newton formulated the principle of conservation
More informationIsentropic Efficiency in Engineering Thermodynamics
June 21, 2010 Isentropic Efficiency in Engineering Thermodynamics Introduction This article is a summary of selected parts of chapters 4, 5 and 6 in the textbook by Moran and Shapiro (2008. The intent
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationDetailed Derivation of Fanno Flow Relationships
Detailed Derivation of Fanno Flow Relationships Matthew MacLean, Ph.D. Version. October 03 v. fixed a sign error on Eq. 5 ! My motivation for writing this note results from preparing course notes for the
More informationApplications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin
Applications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin Faculty of Mechanical Engineering and Naval Architecture, Zagreb 10000, Croatia, gregor.cvijetic@gmail.com The Harmonic
More informationPh.D. Qualifying Exam in Fluid Mechanics
Student ID Department of Mechanical Engineering Michigan State University East Lansing, Michigan Ph.D. Qualifying Exam in Fluid Mechanics Closed book and Notes, Some basic equations are provided on an
More informationNodalization. The student should be able to develop, with justification, a node-link diagram given a thermalhydraulic system.
Nodalization 3-1 Chapter 3 Nodalization 3.1 Introduction 3.1.1 Chapter content This chapter focusses on establishing a rationale for, and the setting up of, the geometric representation of thermalhydraulic
More informationObjectives. Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation
Objectives Conservation of mass principle: Mass Equation The Bernoulli equation Conservation of energy principle: Energy equation Conservation of Mass Conservation of Mass Mass, like energy, is a conserved
More informationModeling of Hydraulic Control Valves
Modeling of Hydraulic Control Valves PETR CHALUPA, JAKUB NOVAK, VLADIMIR BOBAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nam. T.G. Masaryka 5555, 76 1
More informationMathematical Modeling of Chemical Processes. Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University
Mathematical Modeling of Chemical Processes Aisha Osman Mohamed Ahmed Department of Chemical Engineering Faculty of Engineering, Red Sea University Chapter Objectives End of this chapter, you should be
More informationConservation of Mass. Computational Fluid Dynamics. The Equations Governing Fluid Motion
http://www.nd.edu/~gtryggva/cfd-course/ http://www.nd.edu/~gtryggva/cfd-course/ Computational Fluid Dynamics Lecture 4 January 30, 2017 The Equations Governing Fluid Motion Grétar Tryggvason Outline Derivation
More information5 ENERGY EQUATION OF FLUID MOTION
5 ENERGY EQUATION OF FLUID MOTION 5.1 Introduction In order to develop the equations that describe a flow, it is assumed that fluids are subject to certain fundamental laws of physics. The pertinent laws
More informationIJSRD - International Journal for Scientific Research & Development Vol. 1, Issue 8, 2013 ISSN (online):
IJSRD - International Journal for Scientific Research & Development Vol. 1, Issue 8, 2013 ISSN (online): 2321-0613 Development and theoretical analysis of mathematical expressions for change of entropy
More informationA SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations
A SHORT INTRODUCTION TO TWO-PHASE FLOWS Two-phase flows balance equations Hervé Lemonnier DM2S/STMF/LIEFT, CEA/Grenoble, 38054 Grenoble Cedex 9 Ph. +33(0)4 38 78 45 40 herve.lemonnier@cea.fr, herve.lemonnier.sci.free.fr/tpf/tpf.htm
More informationLast name: First name: Student ID: Discussion: You solution procedure should be legible and complete for full credit (use scratch paper as needed).
University of California, Berkeley Mechanical Engineering ME 106, Fluid Mechanics ODK/Midterm 2, Fall 2015 Last name: First name: Student ID: Discussion: Notes: You solution procedure should be legible
More informationNumerical Simulation of Pressure Surge with the Method of Characteristics
Numerical Simulation of Pressure Surge with the Method of Characteristics R. Fiereder 02.04.2009 Saint Petersburg, Russia Content Motivation Governing Equations Continuity Equation Momentum Equation Method
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationThermodynamics I Spring 1432/1433H (2011/2012H) Saturday, Wednesday 8:00am - 10:00am & Monday 8:00am - 9:00am MEP 261 Class ZA
Thermodynamics I Spring 1432/1433H (2011/2012H) Saturday, Wednesday 8:00am - 10:00am & Monday 8:00am - 9:00am MEP 261 Class ZA Dr. Walid A. Aissa Associate Professor, Mech. Engg. Dept. Faculty of Engineering
More informationCentrifugal pumps - characteristics
University of Ljubljana Faculty of mechanical engineering Askerceva 6 1000 Ljubljana, Slovenija telefon: 01 477 1 00 faks: 01 51 85 67 www.fs.uni-lj.si e-mail: dekanat@fs.uni-lj.si Laboratory for Heat
More informationA new solution of Euler s equation of motion with explicit expression of helicity
A new solution of Euler s equation of motion with explicit expression of helicity TSUTOMU KAMBE Former Professor (Physics), University of Tokyo Home: Higashi-yama 2--3, Meguro-ku, Tokyo, 53-0043 JAPAN
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 43 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Treatment of Boundary Conditions These slides are partially based on the recommended textbook: Culbert
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Essential physics for game developers Introduction The primary issues Let s move virtual objects Kinematics: description
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More information34.3. Resisted Motion. Introduction. Prerequisites. Learning Outcomes
Resisted Motion 34.3 Introduction This Section returns to the simple models of projectiles considered in Section 34.1. It explores the magnitude of air resistance effects and the effects of including simple
More informationAerodynamics. Basic Aerodynamics. Continuity equation (mass conserved) Some thermodynamics. Energy equation (energy conserved)
Flow with no friction (inviscid) Aerodynamics Basic Aerodynamics Continuity equation (mass conserved) Flow with friction (viscous) Momentum equation (F = ma) 1. Euler s equation 2. Bernoulli s equation
More informationMass of fluid leaving per unit time
5 ENERGY EQUATION OF FLUID MOTION 5.1 Eulerian Approach & Control Volume In order to develop the equations that describe a flow, it is assumed that fluids are subject to certain fundamental laws of physics.
More informationIn this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor
Lect- 3 In this lecture... Centrifugal compressors Thermodynamics of centrifugal compressors Components of a centrifugal compressor Centrifugal compressors Centrifugal compressors were used in the first
More informationPressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers
Pressure Losses for Fluid Flow Through Abrupt Area Contraction in Compact Heat Exchangers Undergraduate Research Spring 004 By Bryan J. Johnson Under Direction of Rehnberg Professor of Ch.E. Bruce A. Finlayson
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationCHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW
CHAPTER 4 BOUNDARY LAYER FLOW APPLICATION TO EXTERNAL FLOW 4.1 Introduction Boundary layer concept (Prandtl 1904): Eliminate selected terms in the governing equations Two key questions (1) What are the
More information43. A person sits on a freely spinning lab stool that has no friction in its axle. When this person extends her arms,
43. A person sits on a freely spinning lab stool that has no friction in its axle. When this person extends her arms, A) her moment of inertia increases and her rotational kinetic energy remains the same.
More information( ) ( = ) = ( ) ( ) ( )
( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note
More informationVarious lecture notes for
Various lecture notes for 18311. R. R. Rosales (MIT, Math. Dept., 2-337) April 12, 2013 Abstract Notes, both complete and/or incomplete, for MIT s 18.311 (Principles of Applied Mathematics). These notes
More informationDetermination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure
Strojniški vestnik - Journal of Mechanical Engineering 55(009)4, 7-4 Paper received: 7.10.008 UDC 61.64 Paper accepted: 0.04.009 Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering
More informationApplied Gas Dynamics Flow With Friction and Heat Transfer
Applied Gas Dynamics Flow With Friction and Heat Transfer Ethirajan Rathakrishnan Applied Gas Dynamics, John Wiley & Sons (Asia) Pte Ltd c 2010 Ethirajan Rathakrishnan 1 / 121 Introduction So far, we have
More informationn v molecules will pass per unit time through the area from left to
3 iscosity and Heat Conduction in Gas Dynamics Equations of One-Dimensional Gas Flow The dissipative processes - viscosity (internal friction) and heat conduction - are connected with existence of molecular
More informationModels of ocean circulation are all based on the equations of motion.
Equations of motion Models of ocean circulation are all based on the equations of motion. Only in simple cases the equations of motion can be solved analytically, usually they must be solved numerically.
More informationWaves in plasma. Denis Gialis
Waves in plasma Denis Gialis This is a short introduction on waves in a non-relativistic plasma. We will consider a plasma of electrons and protons which is fully ionized, nonrelativistic and homogeneous.
More informationHeriot-Watt University. A pressure-based estimate of synthetic jet velocity Persoons, Tim; O'Donovan, Tadhg. Heriot-Watt University.
Heriot-Watt University Heriot-Watt University Research Gateway A pressure-based estimate of synthetic jet velocity Persoons, Tim; O'Donovan, Tadhg Published in: Physics of Fluids DOI: 10.1063/1.2823560
More informationDynamic Characteristics of Double-Pipe Heat Exchangers
Dynamic Characteristics of Double-Pipe Heat Exchangers WILLIAM C. COHEN AND ERNEST F. JOHNSON Princeton University, Princeton, N. J. The performance of automatically controlled process plants depends on
More informationAE/ME 339. K. M. Isaac. 9/22/2005 Topic 6 FluidFlowEquations_Introduction. Computational Fluid Dynamics (AE/ME 339) MAEEM Dept.
AE/ME 339 Computational Fluid Dynamics (CFD) 1...in the phrase computational fluid dynamics the word computational is simply an adjective to fluid dynamics.... -John D. Anderson 2 1 Equations of Fluid
More informationIn this section, mathematical description of the motion of fluid elements moving in a flow field is
Jun. 05, 015 Chapter 6. Differential Analysis of Fluid Flow 6.1 Fluid Element Kinematics In this section, mathematical description of the motion of fluid elements moving in a flow field is given. A small
More informationChapter 5: Mass, Bernoulli, and Energy Equations
Chapter 5: Mass, Bernoulli, and Energy Equations Introduction This chapter deals with 3 equations commonly used in fluid mechanics The mass equation is an expression of the conservation of mass principle.
More informationDepartment of Energy Fundamentals Handbook. THERMODYNAMICS, HEAT TRANSFER, AND FLUID FLOW, Module 3 Fluid Flow
Department of Energy Fundamentals Handbook THERMODYNAMICS, HEAT TRANSFER, AND FLUID FLOW, Module 3 REFERENCES REFERENCES Streeter, Victor L., Fluid Mechanics, 5th Edition, McGraw-Hill, New York, ISBN 07-062191-9.
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
More informationMODEL BASED PIPELINE MONITORING WITH LEAK DETECTION
MODEL BASED PIPELINE MONITORING WITH LEAK DETECTION Espen Hauge Ole Morten Aamo, John-Morten Godhavn Department of Engineering Cybernetics, Norwegian University of Science and Technology, Trondheim, Norway
More informationImpact of a Jet. Experiment 4. Purpose. Apparatus. Theory. Symmetric Jet
Experiment 4 Impact of a Jet Purpose The purpose of this experiment is to demonstrate and verify the integral momentum equation. The force generated by a jet of water deflected by an impact surface is
More informationMOMENTUM PRINCIPLE. Review: Last time, we derived the Reynolds Transport Theorem: Chapter 6. where B is any extensive property (proportional to mass),
Chapter 6 MOMENTUM PRINCIPLE Review: Last time, we derived the Reynolds Transport Theorem: where B is any extensive property (proportional to mass), and b is the corresponding intensive property (B / m
More informationPhysics 451/551 Theoretical Mechanics. G. A. Krafft Old Dominion University Jefferson Lab Lecture 18
Physics 451/551 Theoretical Mechanics G. A. Krafft Old Dominion University Jefferson Lab Lecture 18 Theoretical Mechanics Fall 018 Properties of Sound Sound Waves Requires medium for propagation Mainly
More informationRational derivation of the Boussinesq approximation
Rational derivation of the Boussinesq approximation Kiyoshi Maruyama Department of Earth and Ocean Sciences, National Defense Academy, Yokosuka, Kanagawa 239-8686, Japan February 22, 2019 Abstract This
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 5
.9 Numerical Fluid Mechanics Fall 011 Lecture 5 REVIEW Lecture 4 Roots of nonlinear equations: Open Methods Fixed-point Iteration (General method or Picard Iteration), with examples Iteration rule: x g(
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationA First Course on Kinetics and Reaction Engineering Example 23.1
Example 23.1 parameter. Problem Purpose This problem illustrates the transient analysis of a CSTR following a change in an operating Problem Statement Recall the isothermal 4430 cm 3 steady-state chemostat
More information1.4 Unit Step & Unit Impulse Functions
1.4 Unit Step & Unit Impulse Functions 1.4.1 The Discrete-Time Unit Impulse and Unit-Step Sequences Unit Impulse Function: δ n = ቊ 0, 1, n 0 n = 0 Figure 1.28: Discrete-time Unit Impulse (sample) 1 [n]
More informationMobile Robot Control on a Reference Path
Proceedings of the 3th Mediterranean Conference on Control and Automation Limassol, Cyprus, June 7-9, 5 WeM5-6 Mobile Robot Control on a Reference Path Gregor Klančar, Drago Matko, Sašo Blažič Abstract
More informationLecture 3: 1. Lecture 3.
Lecture 3: 1 Lecture 3. Lecture 3: 2 Plan for today Summary of the key points of the last lecture. Review of vector and tensor products : the dot product (or inner product ) and the cross product (or vector
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More information0.2. CONSERVATION LAW FOR FLUID 9
0.2. CONSERVATION LAW FOR FLUID 9 Consider x-component of Eq. (26), we have D(ρu) + ρu( v) dv t = ρg x dv t S pi ds, (27) where ρg x is the x-component of the bodily force, and the surface integral is
More informationSimulation of low pressure water hammer
IOP Conference Series: arth and nvironmental Science Simulation of low pressure water hammer To cite this article: D Himr and V Haán 2010 IOP Conf. Ser.: arth nviron. Sci. 12 012087 View the article online
More informationIntroduction to Turbomachinery
1. Coordinate System Introduction to Turbomachinery Since there are stationary and rotating blades in turbomachines, they tend to form a cylindrical form, represented in three directions; 1. Axial 2. Radial
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More information¹ NTNU - Norwegian University of Science and Technology ² Polytec R&D Institute ³ Gassco AS
Accuracy of 1D natural gas transmission models Jan Fredrik Helgaker ¹ T. Ytrehus ¹, A. Oosterkamp ², L. I. Langelandsvik ³, W. Postvoll ³ GERG Academic Network Event, 14th June 2012 ¹ NTNU - Norwegian
More informationModeling and Analysis of Dynamic Systems
Modeling and Analysis of Dynamic Systems Dr. Guillaume Ducard Fall 2017 Institute for Dynamic Systems and Control ETH Zurich, Switzerland G. Ducard c 1 / 21 Outline 1 Lecture 4: Modeling Tools for Mechanical
More informationPIPE FLOWS: LECTURE /04/2017. Yesterday, for the example problem Δp = f(v, ρ, μ, L, D) We came up with the non dimensional relation
/04/07 ECTURE 4 PIPE FOWS: Yesterday, for the example problem Δp = f(v, ρ, μ,, ) We came up with the non dimensional relation f (, ) 3 V or, p f(, ) You can plot π versus π with π 3 as a parameter. Or,
More information1 One-dimensional analysis
One-dimensional analysis. Introduction The simplest models for gas liquid flow systems are ones for which the velocity is uniform over a cross-section and unidirectional. This includes flows in a long
More informationOPAC102. The Acoustic Wave Equation
OPAC102 The Acoustic Wave Equation Acoustic waves in fluid Acoustic waves constitute one kind of pressure fluctuation that can exist in a compressible fluid. The restoring forces responsible for propagating
More information7.5 Elastic Buckling Columns and Buckling
7.5 Elastic Buckling The initial theory of the buckling of columns was worked out by Euler in 1757, a nice example of a theory preceding the application, the application mainly being for the later invented
More informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationHydraulics Prof. Dr. Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati
Hydraulics Prof. Dr. Arup Kumar Sarma Department of Civil Engineering Indian Institute of Technology, Guwahati Module No. # 08 Pipe Flow Lecture No. # 05 Water Hammer and Surge Tank Energy cannot be consumed
More informationDesign of Partial Enclosures. D. W. Herrin, Ph.D., P.E. University of Kentucky Department of Mechanical Engineering
D. W. Herrin, Ph.D., P.E. Department of Mechanical Engineering Reference 1. Ver, I. L., and Beranek, L. L. (2005). Control Engineering: Principles and Applications. John Wiley and Sons. 2. Sharp, B. H.
More informationChapter 1. Introduction to Nonlinear Space Plasma Physics
Chapter 1. Introduction to Nonlinear Space Plasma Physics The goal of this course, Nonlinear Space Plasma Physics, is to explore the formation, evolution, propagation, and characteristics of the large
More informationAnalytical and Numerical Investigation of Transient Gas Blow down
Analytical and Numerical Investigation of Transient Gas Blow down Hameed Balassim Mahood, Farhan Lafta Rasheed & Abdulsattar K. Abbas Ministry of Sciences and Technology, Baghdad, Iraq Received: August
More information4.1 LAWS OF MECHANICS - Review
4.1 LAWS OF MECHANICS - Review Ch4 9 SYSTEM System: Moving Fluid Definitions: System is defined as an arbitrary quantity of mass of fixed identity. Surrounding is everything external to this system. Boundary
More informationHomework #4 Solution. μ 1. μ 2
Homework #4 Solution 4.20 in Middleman We have two viscous liquids that are immiscible (e.g. water and oil), layered between two solid surfaces, where the top boundary is translating: y = B y = kb y =
More informationM E 320 Professor John M. Cimbala Lecture 10. The Reynolds Transport Theorem (RTT) (Section 4-6)
M E 320 Professor John M. Cimbala Lecture 10 Today, we will: Discuss the Reynolds Transport Theorem (RTT) Show how the RTT applies to the conservation laws Begin Chapter 5 Conservation Laws D. The Reynolds
More information