Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws
|
|
- Tamsyn Hannah Ellis
- 5 years ago
- Views:
Transcription
1 Stability and boundary stabilization of physical networks represented by 1-D hyperbolic balance laws G. Bastin
2 Georges Bastin Jean Michel Coron Progress in Nonlinear Differential Equations and Their Applications Subseries in Control 88 Stability and Boundary Stabilization of 1-D Hyperbolic Systems Outline Hyperbolic 2x2 systems of balance laws Steady-state and characteristic form Networks of balance laws Boundary conditions Exponential stability Lyapunov stability analysis Real-life application : hydraulic networks Non uniform systems
3 Open channels, St Venant equations (1871) th + x (hv) = tv + x (gh v2 )=gs Cv 2 /h h = water level, v = water velocity, g= gravity, S = canal slope, C = friction coe cient Road traffic, Aw-Rascle equations (2) t + x ( v) = t(v + p( )) + v x v + p( )) = (V ( ) v)/ = tra c density, v = tra c velocity, p( ) = tra c pressure, V ( ) = preferential velocity, = constant time constant
4 Telegrapher equations, Heaviside t (Li)+@ x v t (Cv)+@ x i = Ri, Gv, i =currentintensity,v = voltage, L = line inductance, C = capacitance, R = line resistance, G = dielectric admittance (per unit length). Optical fibers : Raman amplifiers t S + x S = t P p@ x P = p s S + p P s SP, psp, S(t, x) = transmitted signal power, P (t, x) = pump laser beam power, s and p = propagation velocities, s and p = attenuation coe cients, s and p = amplification gains
5 Euler isentropic equations t x ( v) t ( v)+@ x ( v 2 + P ( )) = C v v = gas density, v = gas velocity, P =pressure,c = friction coe cient Fluid flow in elastic tubes (e.g. blood t A x (AV t (AV )+@ x AV 2 + applea 2 = CV A = tube cross-section, V = fluid velocity,, apple,c = constant coe cients
6 Hyperbolic 2x2 systems of balance laws Space x [,L] Time t [, + ) State Y (t, x) y 1 (t, x) y 2 (t, x) ty + x f(y )=g(y ) ty + A(Y ) x Y = g(y ) A(Y ) has 2 distinct real eigenvalues
7 Uniform steady state ty + A(Y ) x Y = g(y ) A steady-state is a constant solution Y (t, x) Ȳ which satisfies the equation g(ȳ )= and (obviously) the state equation tȳ + A(Ȳ ) xȳ = g(ȳ ) Road traffic = density v = velocity Open channels h = water depth v = velocity Steady state v = V ( ) Steady state r gs v = C h (Toricelli formula)
8 Characteristic form Hyperbolic system : ty + A(Y ) x Y = g(y ) Change of coordinates : (Y )= 1 2 ( Riemann coordinates t + c1 ( ) c 2 ( ) with c 1 ( ) = c 2 ( ) eigenvalues of A(Y x = h( ) ( Characteristic velocities ) The change of coordinates (Y ) is defined up to a constant. It can therefore be selected such that (Ȳ )= h() =.
9 Physical netwoks of 2x2 hyperbolic systems (e.g. hydraulic networks (irrigation, waterways) or road tra c networks) directed graph n edges one 2x2 hyperbolic system of balance laws attached to each arc L x arrow = orientation of the x t i n+i + ci ( ) c n+i ( ) x n+i = h i n+i (i =1,...,n)
10 Characteristic t + C( )@ x = h( ) (characteristic velocities) C( ) = diag(c 1 ( ),...,c 2n ( )) c i ( ) > c i ( ) < i =1,...,n i = n +1,...,2n Notations + =( 1, 2,..., n ) Boundary conditions + (t, ) (t, L) =( n+1,..., 2n ) = G + (t, L) (t, ) x 2 [,L] L
11 Boundary conditions = Physical constraints + (t, ) (t, L) = G + (t, L) (t, ) F (Y (t, ),Y(t, L)) = Road traffic 1 = density v = velocity 3 Open channels h = water depth v = velocity Flow conservation at a junction 3(t, )v 3 (t, ) = 1(t, L)v 1 (t, L)+ 2 (t, L)v 2 (t, L) Modelling of hydraulic gates h 2 (t, )v 2 (t, ) = h 1 (t, L) u 3/2
12 Boundary conditions = boundary feedback control Road traffic = density v = velocity q = v = flux Open channels h = water depth v = velocity 1 3 h 2 Feedback implementation of ramp metering Feedback control of water depth in navigable rivers G function of the control tuning parameters : How to design the control laws to make the boundary conditions stabilizing?
13 t 2 [, +1) x 2 t + C( )@ x = h( ) + (t, ) + (t, L) = G (t, L) (t, ) (,x)= (x) Conditions on C, h and G such that =is exponentially stable?
14 t 2 [, +1) x 2 t + C( )@ x = h( ) + (t, ) + (t, L) = G (t, L) (t, ) (,x)= (x) Conditions on C, h and G such that =is exponentially stable? In a linear space of functions from [,L] into R 2n Definition of exponential stability for the norm k k 9 ",, such that solutions are defined 8 t 2 [, +1) and k k 6 " )k (t,.)k 6 e t k k (+ usual compatibility conditions adapted to )
15 Lyapunov stability : we start with the linear case System Boundary condition t + x = B (t, ) = K (t, L) = diag{ i > } Lyapunov function V = Z L T (t, x)p (t, x)e µx dx µ> P = diag {p i > } V = Z L (t, x) µp [B T P + PB] (t, x)e µx dx T (t, L)[P e µl K T P K] (t, L)
16 System t + x = B (t, ) = K (t, L) (,x)= (x) V = V = Z L Z L Lyapunov function T (t, x)p (t, x)e µx dx (t, x) µp [B T P + PB] (t, x)e µx dx T (t, L)[P e µl K T P K] (t, L) If 9 P (diag. pos.) and µ> such that: 1) µp [B T P + PB] 2) P e µl K T P K ( : positive definite) Then 9 > s.t. V 6 V ) exponential stability for L 2 -norm (B-C chapter 5)
17 Let us now consider the nonlinear case System Boundary t + C( )@ x = h( ) c i ( ) > (t, ) = G( (t, L)) Notations: = C(), B = h (), K = G () (linearization) If 9 P (diag. pos.) and µ> such that: 1) µp [B T P + PB] 2) P e µl K T P K then the steady-state is exponentially stable for the H 2 -norm (B-C chapter 6)
18 Let us now consider the nonlinear case System Boundary t + C( )@ x = h( ) c i ( ) > (t, ) = G( (t, L)) Notations: = C(), B = h (), K = G () (linearization) If 9 P (diag. pos.) and µ> such that: 1) µp [B T P + PB] 2) P e µl K T P K then the steady-state is exponentially stable for the H 2 -norm Lyapunov function V = Z L ( T P + t P t + tt P tt )e µx dx (B-C chapter 6)
19 Special case: a more explicit stability t + C( )@ x = h( ) (t, ) = G( (t, L)) = C(), B = h (), K = G () Theorem If B =,orifkbksu ciently small, then is exponentially stable for H 2 -norm if 2 (K) < 1 2 (K) =min(k K 1 k, positive diagonal) (boundary damping) (t, ) = K (t, L) (k k : 2-norm) (B-C chapters 4 and 6)
20 Special case: Li Ta Tsien t + C( )@ x = h( ) (t, ) = G( (t, L)) = C(), B = h (), K = G () Theorem If B =,orifkbk su ciently small, then is exponentially stable for C 1 -norm if ( K ) < 1 ( K ) = spectral radius of the matrix h i K ij (boundary damping) (t, ) = K (t, L) (Li Ta Tsien 1994 ) (B-C chapter 6)
21 H 2 /C 1 exponential t + C( )@ x = h( ) (t, ) = G( (t, L)) = C(), B = h (), K = G () For every K 2 M m,m 2 (K) 6 ( K ) Example: for K = 1 1, (K) = p 2 < 2= ( K ) (B-C chapter 4 + Appendix) For the C 1 -norms, the condition ( K ) < 1 is su cient for the stability the condition 2 (K) < 1isnot su cient (Coron and Nguyen 215)
22 H 2 /C 1 exponential t + C( )@ x = h( ) (t, ) = G( (t, L)) = C(), B = h (), K = G () For every K 2 M m,m Example: for K = 2 (K) 6 ( K ) 1 1, (K) = p 2 < 2= ( K ) (B-C chapter 4 + Appendix) For the C 1 -norms, the condition ( K ) < 1 is su cient for the stability the condition 2 (K) < 1isnot su cient There are boundary conditions that are su ciently damping for the H 2 -norm but not for the C 1 -norm
23 Another special t + C( )@ x = h( ) (t, ) = G( (t, L)) = C(), B = h (), K = G () If 9 P (diag. pos.) such that: 1) B T P + PB 4 (e.g. dissipative friction term) 2) 2 (K) =(k K 1 k < 1 with 2 = P (boundary damping)
24 The stability conditions hold also, with appropriate modifications, for systems with positive and negative characteristic velocities c i ( ), as we shall see in the example.
25 Example: Boundary control for a channel with multiple pools Pool i H i V i u i H i+1 V i+1 x x Pool i +1 u i+1 Meuse river (Belgium)
26 Real-life application : level control in navigable rivers. Meuse river (Belgium).
27 Model : Saint-Venant equations Pool i H i V i u i H i+1 V i+1 x x Pool i +1 u i+1 t Hi V i + x 1 2 V i 2 H i V i + gh i Boundary conditions = g[s i C i V 2 i H 1 slope i ] friction, i =1...,n. (width = W ) 1) Conservation of flows H i (t, L)V i (t, L) =H i+1 (t, )V i+1 (t, ) i =1,...,n 1 H i (t, L)V i (t, L) =k G H i (t, L) u i (t) 3 2) Gate models i =1,...,n 3) Input flow WH 1 (t, )V 1 (t, ) = Q (t) = disturbance input in navigable rivers or control input in irrigation networks
28 Riemann coordinates i =(V i V i )+(H i H i ) Boundary control design g H i n+i =(V i V i ) (H i H i ) (H i,v i ) = steady-state g H i i =1,..., n. Control objective n+i(t, L) = k i i (t, L) i =1,...,n k i = control tuning parameters Gate model H i (t, L)V i (t, L) =k G H i (t, L) u i (t) 3 i =1,...,n u i (t) =H i (t, L) Control law H i H i (t, L) 1 k i g Si H (H i (t, L) H k G 1+k i ) + i i H i C = level set points 2/3 i =1,...,n
29 Linearized Saint-Venant equations in Riemann coordinates i t i n+i x i x i n+i = i i i i i n+i n+i = V i gh i < < i = V i + gh i < n+i < i < i = gs i 1 V i 1 2 gh i < i = gs i 1 V i gh i
30 Lyapunov t i n+i x i n+i Subcritical flow x i n+i = < n+i < i < i < i i i i i i n+i V = n L i=1 p i 2 i e µx + p n+i 2 n+ie µx dx, p i,p n+i,µ> If parameters p i selected such that p i i = p n+i i p i+1 = p i and µ su ciently small (dissipative friction) control tuning parameters k i such that k 2 i i i n+i i < 1 (boundary damping) steady-state exponentially stable (B-C chapter 5)
31 Experimental result (Dinant) meters above sea level change of set point water level low flow rate high flow rate gate levels = control action 16 to 23 october 212
32 The nonuniform case : example Now we consider a pool of a prismatic horizontal open channel with a rectangular cross section and a unit width. H(x, t) = water depth V (x, t) = water velocity C = friction coe cient th + x (HV )= V 2 tv + x 2 + gh + gc V 2 H = Steady-state H (x),v (x) for a constant flow rate Q = H (x)v (x) dv dx = gc Q (V (x)) 5 gq (V (x)) 3 Nonuniform steady state H (x) V (x) L x
33 th + x (HV )= V 2 tv + x 2 + gh + gc V 2 H = Linearization about the steady-state h(t, x) :=H(t, x) H (x), v(t, x) :=V (t, x) V (x) Linearized model in physical coordinates h t v t + V H g V h x v x + V x gc V 2 H H x V x +2gC V H h = v
34 h t v t + V H g V h x v x + V x gc V 2 H H x V x +2gC V H h v = Lyapunov function V = Z L gh 2 + H v 2 dx = Z L h g 1 A H h v dx dv dt = Z L Y T N(x)Y dx Y T M(x)Y L Y = h v g 2 CV 3 H (gh V 2 ) N(x) = gcv 2 H 1 gcv 2 H C 2gCV 3 (gh V 2 ) +4gCV A M(x) = gv gh gh H V
35 h t v t + V H g V h x v x + V x gc V 2 H H x V x +2gC V H h v = Lyapunov function V = Subcritical flow (i.e. fluvial): gh V 2 > Z L gh 2 + H v 2 dx = Z L h g 1 A H h v dx dv dt = Z L Y T N(x)Y dx Y T M(x)Y L Y = h v g 2 CV 3 H (gh V 2 ) N(x) = gcv 2 H 1 gcv 2 H C 2gCV 3 (gh V 2 ) +4gCV A M(x) = gv gh gh H V positive definite
36 h t v t + V H g V h x v x + V x gc V 2 H H x V x +2gC V H h v = Lyapunov function V = Subcritical flow (i.e. fluvial): gh V 2 > Z L gh 2 + H v 2 dx = Z L h g 1 A H h v dx dv dt = Z L Y T N(x)Y dx Y T M(x)Y L boundary damping conditions? Y = h v g 2 CV 3 H (gh V 2 ) N(x) = gcv 2 H 1 gcv 2 H C 2gCV 3 (gh V 2 ) +4gCV A M(x) = gv gh gh H V positive definite
37 disturbance Q (t) H(t, x) V (t, x) control Q L (t) L x control Q L (t) =Q (t)+k P (H(t, L) H sp ) feedforward feedback
38 disturbance Q (t) H(t, x) V (t, x) control Q L (t) L x measured measured set point control Q L (t) =Q (t)+k P (H(t, L) H sp ) control tuning
39 disturbance Q (t) H(t, x) V (t, x) control Q L (t) L x control Q L (t) =Q (t)+k P (H(t, L) H sp ) Y T M(x)Y L > =) > boundary damping V apple () gh H () (gh () V 2 ())h 2 (L) V 2 (L) (t, )+ H (2k P V (L)) + V (L) (L) H (L) k2 P h 2 (t, L) =) exponential stability if k P is su ciently large...
40 back : closed loop front : open loop Thank you Meuse river
41
ON LYAPUNOV STABILITY OF LINEARISED SAINT-VENANT EQUATIONS FOR A SLOPING CHANNEL. Georges Bastin. Jean-Michel Coron. Brigitte d Andréa-Novel
NETWORKS AND HETEROGENEOUS MEDIA doi:.3934/nhm.9.4.77 c American Institute of Mathematical Sciences Volume 4, Number, June 9 pp. 77 87 ON LYAPUNOV STABILITY OF LINEARISED SAINT-VENANT EQUATIONS FOR A SLOPING
More informationCoupling conditions for transport problems on networks governed by conservation laws
Coupling conditions for transport problems on networks governed by conservation laws Michael Herty IPAM, LA, April 2009 (RWTH 2009) Transport Eq s on Networks 1 / 41 Outline of the Talk Scope: Boundary
More informationExponential stability of PI control for Saint-Venant equations
Exponential stability of PI control for Saint-Venant equations Georges Bastin and Jean-Michel Coron February 28, 218 Abstract We consider open channels represented by Saint-Venant equations that are monitored
More informationStability and output regulation for a cascaded network of 2 2 hyperbolic systems with PI control
Stability and output regulation for a cascaded network of 2 2 hyperbolic systems with PI control Ngoc-Tu TRINH, Vincent ANDRIEU and Cheng-Zhong XU Laboratory LAGEP, Batiment CPE, University of Claude Bernard
More informationRobust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations
1 Robust boundary control of systems of conservation laws: theoretical study with numerical and experimental validations LAAS - CNRS Toulouse (France) Workshop on Irrigation Channels and Related Problems
More informationStability and Boundary Stabilization of 1-D Hyperbolic Systems. Georges Bastin and Jean-Michel Coron
Stability and Boundary Stabilization of 1-D Hyperbolic Systems Georges Bastin and Jean-Michel Coron April 18, 216 2 c G. Bastin and J-M. Coron, 216. Contents Preface 9 1 Hyperbolic Systems of Balance Laws
More informationBoundary Feedback Control and Lyapunov Stability Analysis for Physical Networks of 2 2 Hyperbolic Balance Laws
Boundary Feedback Control and Lyapunov Stability Analysis for Physical Networks of 2 2 Hyperbolic Balance Laws Georges Bastin, Jean-Michel Coron, Brigitte D Andréa-Novel To cite this version: Georges Bastin,
More informationExponential Decay: From Semi-Global to Global
Exponential Decay: From Semi-Global to Global Martin Gugat Benasque 2013 Martin Gugat (FAU) From semi-global to global. 1 / 13 How to get solutions that are global in time? 1 Sometimes, the semiglobal
More informationSome results on the stabilization and on the controllability of hyperbolic systems
Some results on the stabilization and on the controllability of hyperbolic systems Jean-Michel Coron http://www.ann.jussieu.fr/~coron/ Irrigation Channels and Related Problems Maiori (Salerno, Italy) October
More informationBoundary Control of Open Channels With Numerical and Experimental Validations Valérie Dos Santos and Christophe Prieur
1252 IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 16, NO. 6, NOVEMBER 2008 Boundary Control of Open Channels With Numerical and Experimental Validations Valérie Dos Santos and Christophe Prieur
More informationStabilité des systèmes hyperboliques avec des commutations
Stabilité des systèmes hyperboliques avec des commutations Christophe PRIEUR CNRS, Gipsa-lab, Grenoble, France Séminaire EDP, Université de Versailles 9 octobre 2014 1/34 C. Prieur UVSQ, Octobre 2014 2/34
More informationStabilité des systèmes hyperboliques avec des commutations
Stabilité des systèmes hyperboliques avec des commutations Christophe PRIEUR CNRS, Gipsa-lab, Grenoble, France Contrôles, Problèmes inverses et Applications Clermont-Ferrand Septembre 2014 1/34 C. Prieur
More informationBlock 3 Open channel flow
Numerical Hydraulics Block 3 Open channel flow Markus Holzner Contents of the course Block 1 The equations Block Computation of pressure surges Block 3 Open channel flow (flow in rivers) Block 4 Numerical
More informationIndex. Index. More information. in this web service Cambridge University Press
A-type elements, 4 7, 18, 31, 168, 198, 202, 219, 220, 222, 225 A-type variables. See Across variable ac current, 172, 251 ac induction motor, 251 Acceleration rotational, 30 translational, 16 Accumulator,
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source of preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationIndex. Chromatography, 3 Condition of C 0 compatibility, 8, 9 Condition of C 1 compatibility, 8, 9, 33, 36, 38, 39, 41, 43, 50, 53, 56, 58, 59, 62
References 1. Bressan, A., Čanić, S., Garavello, M., Herty, M., Piccoli, B.: Flows on networks: recent results and perspectives. EMS Surv. Math. Sci. 1, 47 111 (2014) 2. Coron, J.-M., Wang, Z.: Controllability
More informationNumerical Hydraulics
ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow
More informationMath 308 Exam I Practice Problems
Math 308 Exam I Practice Problems This review should not be used as your sole source for preparation for the exam. You should also re-work all examples given in lecture and all suggested homework problems..
More informationOPEN CHANNEL FLOW. Computer Applications. Numerical Methods and. Roland Jeppson. CRC Press UNIVERSITATSB'BUOTHEK TECHNISCHE. INFORMATlONSBiBUOTHEK
OPEN CHANNEL FLOW Numerical Methods and Computer Applications Roland Jeppson TECHNISCHE INFORMATlONSBiBUOTHEK UNIVERSITATSB'BUOTHEK HANNOVER Si. i. CRC Press Taylor &.Francis Group Boca Raton London New
More informationLyapunov functions for switched linear hyperbolic systems
Lyapunov functions for switched linear hyperbolic systems Christophe PRIEUR Antoine GIRARD and Emmanuel WITRANT UJF-Grenoble 1/CNRS, Grenoble Image Parole Signal Automatique (GIPSA-lab), UMR 5216, B.P.
More informationSolutions of Math 53 Midterm Exam I
Solutions of Math 53 Midterm Exam I Problem 1: (1) [8 points] Draw a direction field for the given differential equation y 0 = t + y. (2) [8 points] Based on the direction field, determine the behavior
More informationElementary Differential Equations
Elementary Differential Equations George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 310 George Voutsadakis (LSSU) Differential Equations January 2014 1 /
More informationPhys101 First Major-061 Zero Version Coordinator: Abdelmonem Monday, October 30, 2006 Page: 1
Coordinator: Abdelmonem Monday, October 30, 006 Page: 1 Q1. An aluminum cylinder of density.70 g/cm 3, a radius of.30 cm, and a height of 1.40 m has the mass of: A) 6.8 kg B) 45.1 kg C) 13.8 kg D) 8.50
More informationInfinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems
Infinite-dimensional nonlinear predictive controller design for open-channel hydraulic systems D. Georges, Control Systems Dept - Gipsa-lab, Grenoble INP Workshop on Irrigation Channels and Related Problems,
More information1.060 Engineering Mechanics II Spring Problem Set 8
1.060 Engineering Mechanics II Spring 2006 Due on Monday, May 1st Problem Set 8 Important note: Please start a new sheet of paper for each problem in the problem set. Write the names of the group members
More informationVARIED FLOW IN OPEN CHANNELS
Chapter 15 Open Channels vs. Closed Conduits VARIED FLOW IN OPEN CHANNELS Fluid Mechanics, Spring Term 2011 In a closed conduit there can be a pressure gradient that drives the flow. An open channel has
More informationClosed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationFirst order differential equations
First order differential equations Samy Tindel Purdue University Differential equations and linear algebra - MA 262 Taken from Differential equations and linear algebra by Goode and Annin Samy T. First
More informationClosed duct flows are full of fluid, have no free surface within, and are driven by a pressure gradient along the duct axis.
OPEN CHANNEL FLOW Open channel flow is a flow of liquid, basically water in a conduit with a free surface. The open channel flows are driven by gravity alone, and the pressure gradient at the atmospheric
More informationBasic Theory of Differential Equations
page 104 104 CHAPTER 1 First-Order Differential Equations 16. The following initial-value problem arises in the analysis of a cable suspended between two fixed points y = 1 a 1 + (y ) 2, y(0) = a, y (0)
More informationLyapunov-based Stability Analysis for a MIMO Counter-Propagating Raman Amplifier. Daniel Beauchamp
Lyapunov-based Stability Analysis for a MIMO Counter-Propagating Raman Amplifier by Daniel Beauchamp A thesis submitted in conformity with the requirements for the degree of Master of Applied Science Graduate
More information6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities
6.207/14.15: Networks Lectures 4, 5 & 6: Linear Dynamics, Markov Chains, Centralities 1 Outline Outline Dynamical systems. Linear and Non-linear. Convergence. Linear algebra and Lyapunov functions. Markov
More information2.004 Dynamics and Control II Spring 2008
MIT OpenCourseWare http://ocwmitedu 00 Dynamics and Control II Spring 00 For information about citing these materials or our Terms of Use, visit: http://ocwmitedu/terms Massachusetts Institute of Technology
More informationThe Impulse-Momentum Principle
Chapter 6 /60 The Impulse-Momentum Principle F F Chapter 6 The Impulse-Momentum Principle /60 Contents 6.0 Introduction 6. The Linear Impulse-Momentum Equation 6. Pipe Flow Applications 6.3 Open Channel
More informationdt 2 The Order of a differential equation is the order of the highest derivative that occurs in the equation. Example The differential equation
Lecture 18 : Direction Fields and Euler s Method A Differential Equation is an equation relating an unknown function and one or more of its derivatives. Examples Population growth : dp dp = kp, or = kp
More informationLecture 12: Transcritical flow over an obstacle
Lecture 12: Transcritical flow over an obstacle Lecturer: Roger Grimshaw. Write-up: Erinna Chen June 22, 2009 1 Introduction The flow of a fluid over an obstacle is a classical and fundamental problem
More informationFirst Order Differential Equations Lecture 3
First Order Differential Equations Lecture 3 Dibyajyoti Deb 3.1. Outline of Lecture Differences Between Linear and Nonlinear Equations Exact Equations and Integrating Factors 3.. Differences between Linear
More informationLinear Hyperbolic Systems
Linear Hyperbolic Systems Professor Dr E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8, 2014 1 / 56 We study some basic
More informationLecture Note for Open Channel Hydraulics
Chapter -one Introduction to Open Channel Hydraulics 1.1 Definitions Simply stated, Open channel flow is a flow of liquid in a conduit with free space. Open channel flow is particularly applied to understand
More informationApplication of the Muskingum-Cunge method for dam break flood routing F. Macchione Dipartimento di Difesa del Suolo, Universita delta Calabria,
Application of the Muskingum-Cunge method for dam break flood routing F. Macchione Dipartimento di Difesa del Suolo, Universita delta Calabria, Abstract This paper deals with the application of the Muskingum-Cunge
More informationBasic Procedures for Common Problems
Basic Procedures for Common Problems ECHE 550, Fall 2002 Steady State Multivariable Modeling and Control 1 Determine what variables are available to manipulate (inputs, u) and what variables are available
More informationNotes on RKDG methods for shallow-water equations in canal networks
Notes on RKDG methods for shallow-water equations in canal networks Maya Briani, Benedetto Piccoli, Jing-Mei Qiu December 9, 015 Abstract PDE models for network flows are used in a number of different
More informationChapter 2 Notes, Kohler & Johnson 2e
Contents 2 First Order Differential Equations 2 2.1 First Order Equations - Existence and Uniqueness Theorems......... 2 2.2 Linear First Order Differential Equations.................... 5 2.2.1 First
More informationPassivity-based Control of Euler-Lagrange Systems
Romeo Ortega, Antonio Loria, Per Johan Nicklasson and Hebertt Sira-Ramfrez Passivity-based Control of Euler-Lagrange Systems Mechanical, Electrical and Electromechanical Applications Springer Contents
More informationMath Applied Differential Equations
Math 256 - Applied Differential Equations Notes Existence and Uniqueness The following theorem gives sufficient conditions for the existence and uniqueness of a solution to the IVP for first order nonlinear
More informationLaboratory studies on colliding gravity currents. Qiang ZHONG. Environmental Fluid Dynamics Group University of Notre Dame October
Laboratory studies on colliding gravity currents Qiang ZHONG Environmental Fluid Dynamics Group University of Notre Dame October 8 25 Outlines Introduction 2 Experiments 3 Currently Result 4 Conclusion
More informationControl system design for concrete irrigation channels
Control system design for concrete irrigation channels Timm Strecker 1, Michael Cantoni 2 and Ole Morten Aamo 1 Abstract Concrete channels find use at the periphery of irrigation networks, for expansion
More informationDifferential Equation (DE): An equation relating an unknown function and one or more of its derivatives.
Lexicon Differential Equation (DE): An equation relating an unknown function and one or more of its derivatives. Ordinary Differential Equation (ODE): A differential equation that contains only ordinary
More informationMath 2a Prac Lectures on Differential Equations
Math 2a Prac Lectures on Differential Equations Prof. Dinakar Ramakrishnan 272 Sloan, 253-37 Caltech Office Hours: Fridays 4 5 PM Based on notes taken in class by Stephanie Laga, with a few added comments
More informationP = 2Rθ. The previous Manning formulas are used to predict V o and Q for uniform flow when the above expressions are substituted for A, P, and R h.
Uniform Flow in a Partly Full, Circular Pipe Fig. 10.6 shows a partly full, circular pipe with uniform flow. Since frictional resistance increases with wetted perimeter, but volume flow rate increases
More information20D - Homework Assignment 4
Brian Bowers (TA for Hui Sun) MATH 0D Homework Assignment November, 03 0D - Homework Assignment First, I will give a brief overview of how to use variation of parameters. () Ensure that the differential
More informationSection 11.1 What is a Differential Equation?
1 Section 11.1 What is a Differential Equation? Example 1 Suppose a ball is dropped from the top of a building of height 50 meters. Let h(t) denote the height of the ball after t seconds, then it is known
More informationProblems in VLSI design
Problems in VLSI design wire and transistor sizing signal delay in RC circuits transistor and wire sizing Elmore delay minimization via GP dominant time constant minimization via SDP placement problems
More informationA Multiphysics Strategy for Free Surface Flows
A Multiphysics Strategy for Free Surface Flows Edie Miglio, Simona Perotto, and Fausto Saleri MOX, Modeling and Scientific Computing, Department of Mathematics, Politecnico of Milano, via Bonardi 9, I-133
More informationRutgers University Department of Physics & Astronomy. 01:750:271 Honors Physics I Fall Lecture 7 JJ II. Home Page. Title Page.
Rutgers University Department of Physics & Astronomy 01:750:271 Honors Physics Fall 2015 Lecture 7 Page 1 of 31 Midterm 1: Monday October 10th Motion in one, two and three dimensions Forces and Motion
More informationAppendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2)
Appendix A: Exercise Problems on Classical Feedback Control Theory (Chaps. 1 and 2) For all calculations in this book, you can use the MathCad software or any other mathematical software that you are familiar
More informationNumerical simulation of shallow water equations and comparison with experimental data.
Numerical simulation of shallow water equations and comparison with experimental data. Mikolaj Szydlarski, Chiara Simeoni To cite this version: Mikolaj Szydlarski, Chiara Simeoni. Numerical simulation
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationSinusoidal Steady State Analysis (AC Analysis) Part II
Sinusoidal Steady State Analysis (AC Analysis) Part II Amin Electronics and Electrical Communications Engineering Department (EECE) Cairo University elc.n102.eng@gmail.com http://scholar.cu.edu.eg/refky/
More informationQ. 1 Q. 25 carry one mark each.
GATE 5 SET- ELECTRONICS AND COMMUNICATION ENGINEERING - EC Q. Q. 5 carry one mark each. Q. The bilateral Laplace transform of a function is if a t b f() t = otherwise (A) a b s (B) s e ( a b) s (C) e as
More informationChapter 6 The Impulse-Momentum Principle
Chapter 6 The Impulse-Momentum Principle 6. The Linear Impulse-Momentum Equation 6. Pipe Flow Applications 6.3 Open Channel Flow Applications 6.4 The Angular Impulse-Momentum Principle Objectives: - Develop
More informationYour exam contains 5 problems. The entire exam is worth 70 points. Your exam should contain 6 pages; please make sure you have a complete exam.
MATH 124 (PEZZOLI) WINTER 2017 MIDTERM #2 NAME TA:. Section: Instructions: Your exam contains 5 problems. The entire exam is worth 70 points. Your exam should contain 6 pages; please make sure you have
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Linearization and Characteristic Relations 1 / 31 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
More informationHW2 Solutions. MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22]
HW2 Solutions MATH 20D Fall 2013 Prof: Sun Hui TA: Zezhou Zhang (David) October 14, 2013 Checklist: Section 2.6: 1, 3, 6, 8, 10, 15, [20, 22] Section 3.1: 1, 2, 3, 9, 16, 18, 20, 23 Section 3.2: 1, 2,
More informationModule 4: One-Dimensional Kinematics
4.1 Introduction Module 4: One-Dimensional Kinematics Kinematics is the mathematical description of motion. The term is derived from the Greek word kinema, meaning movement. In order to quantify motion,
More informationYell if you have any questions
Class 36: Outline Hour 1: Concept Review / Overview PRS Questions Possible Exam Questions Hour : Sample Exam Yell if you have any questions P36-1 Before Starting All of your grades should now be posted
More informationSummer School on Mathematical Control Theory. Control of biological processes
united nations educational, scientific and cultural organization the i international centre for theoretical physics international atomic energy agency SMR1 327/25 Summer School on Mathematical Control
More informationDifferential Equations & Separation of Variables
Differential Equations & Separation of Variables SUGGESTED REFERENCE MATERIAL: As you work through the problems listed below, you should reference Chapter 8. of the recommended textbook (or the equivalent
More informationChapter 7 Applications of Integration
Chapter 7 Applications of Integration 7.1 Area of a Region Between Two Curves 7.2 Volume: The Disk Method 7.3 Volume: The Shell Method 7.4 Arc Length and Surfaces of Revolution 7.5 Work 7.6 Moments, Centers
More informationScattering Parameters
Berkeley Scattering Parameters Prof. Ali M. Niknejad U.C. Berkeley Copyright c 2016 by Ali M. Niknejad September 7, 2017 1 / 57 Scattering Parameters 2 / 57 Scattering Matrix Voltages and currents are
More informationA system is defined as a combination of components (elements) that act together to perform a certain objective. System dynamics deal with:
Chapter 1 Introduction to System Dynamics A. Bazoune 1.1 INTRODUCTION A system is defined as a combination of components (elements) that act together to perform a certain objective. System dynamics deal
More informationDynamic Systems. Modeling and Analysis. Hung V. Vu. Ramin S. Esfandiari. THE McGRAW-HILL COMPANIES, INC. California State University, Long Beach
Dynamic Systems Modeling and Analysis Hung V. Vu California State University, Long Beach Ramin S. Esfandiari California State University, Long Beach THE McGRAW-HILL COMPANIES, INC. New York St. Louis San
More informationMATH 32A: MIDTERM 1 REVIEW. 1. Vectors. v v = 1 22
MATH 3A: MIDTERM 1 REVIEW JOE HUGHES 1. Let v = 3,, 3. a. Find e v. Solution: v = 9 + 4 + 9 =, so 1. Vectors e v = 1 v v = 1 3,, 3 b. Find the vectors parallel to v which lie on the sphere of radius two
More informationDifferential Equations Spring 2007 Assignments
Differential Equations Spring 2007 Assignments Homework 1, due 1/10/7 Read the first two chapters of the book up to the end of section 2.4. Prepare for the first quiz on Friday 10th January (material up
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 42 Flows with a Free Surface Part II Good morning. I welcome you to this session
More informationAP Physics C Mechanics Objectives
AP Physics C Mechanics Objectives I. KINEMATICS A. Motion in One Dimension 1. The relationships among position, velocity and acceleration a. Given a graph of position vs. time, identify or sketch a graph
More informationFluvial Dynamics. M. I. Bursik ublearns.buffalo.edu October 26, Home Page. Title Page. Contents. Page 1 of 18. Go Back. Full Screen. Close.
Page 1 of 18 Fluvial Dynamics M. I. Bursik ublearns.buffalo.edu October 26, 2008 1. Fluvial Dynamics We want to understand a little of the basic physics of water flow and particle transport, as so much
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationMATH 4B Differential Equations, Fall 2016 Final Exam Study Guide
MATH 4B Differential Equations, Fall 2016 Final Exam Study Guide GENERAL INFORMATION AND FINAL EXAM RULES The exam will have a duration of 3 hours. No extra time will be given. Failing to submit your solutions
More informationFOUNDATION STUDIES EXAMINATIONS April PHYSICS First Paper February Program 2007
FOUNDATION STUDIES EXAMINATIONS April 2007 HYSICS First aper February rogram 2007 Time allowed hour for writing 0 minutes for reading This paper consists of 3 questions printed on 5 pages. LEASE CHECK
More informationPUTNAM PROBLEMS DIFFERENTIAL EQUATIONS. First Order Equations. p(x)dx)) = q(x) exp(
PUTNAM PROBLEMS DIFFERENTIAL EQUATIONS First Order Equations 1. Linear y + p(x)y = q(x) Muliply through by the integrating factor exp( p(x)) to obtain (y exp( p(x))) = q(x) exp( p(x)). 2. Separation of
More informationChapter 4 DYNAMICS OF FLUID FLOW
Faculty Of Engineering at Shobra nd Year Civil - 016 Chapter 4 DYNAMICS OF FLUID FLOW 4-1 Types of Energy 4- Euler s Equation 4-3 Bernoulli s Equation 4-4 Total Energy Line (TEL) and Hydraulic Grade Line
More information1 Motion of a single particle - Linear momentum, work and energy principle
1 Motion of a single particle - Linear momentum, work and energy principle 1.1 In-class problem A block of mass m slides down a frictionless incline (see Fig.). The block is released at height h above
More informationBeaver Creek Corridor Design and Analysis. By: Alex Previte
Beaver Creek Corridor Design and Analysis By: Alex Previte Overview Introduction Key concepts Model Development Design Accuracy Conclusion Refresh v = Beaver Creek Site = Wittenberg Introduction Low head
More informationTwo-Dimensional Simulation of Subcritical Flow at a Combining Junction: Luxury or Necessity?
Two-Dimensional Simulation of Subcritical Flow at a Combining Junction: Luxury or Necessity? Rabih Ghostine 1 ; Robert Mose ; José Vazquez 3 ; Abdellah Ghenaim 4 ; and Caroline Grégoire 5 Abstract: Classically,
More informationLab 7: Nonuniform Flow and Open Channel Transitions
CE 3620: Water Resources Engineering Spring 2015 Lab 7: Nonuniform Flow and Open Channel Transitions BACKGROUND An open channel transition may be defined as a change either in the direction, slope, or
More informationFeedback Control Systems
ME Homework #0 Feedback Control Systems Last Updated November 06 Text problem 67 (Revised Chapter 6 Homework Problems- attached) 65 Chapter 6 Homework Problems 65 Transient Response of a Second Order Model
More information= Q:An Qn% icx=zv. A, with Bn = T- n. Modelling of irrigation channel dynamics for controller design
Modelling of irrigation channel dynamics for controller design Jean-Pierre BAUME Jacques SAU Pierre-Olivier Malaterre Cemagref, Division Irrigation, 361 rue J-F. ISTIL Bat. 201, RC, UniversitC Claude Cemagref,
More informationIntegrator Windup
3.5.2. Integrator Windup 3.5.2.1. Definition So far we have mainly been concerned with linear behaviour, as is often the case with analysis and design of control systems. There is, however, one nonlinear
More informationMath 4B Notes. Written by Victoria Kala SH 6432u Office Hours: T 12:45 1:45pm Last updated 7/24/2016
Math 4B Notes Written by Victoria Kala vtkala@math.ucsb.edu SH 6432u Office Hours: T 2:45 :45pm Last updated 7/24/206 Classification of Differential Equations The order of a differential equation is the
More information28.2 Classification of Jumps
28.2 Classification of Jumps As mentioned earlier, the supercritical flow Froude number influences the characteristics of the hydraulic jump. Bradley and Peterka, after extensive experimental investigations,
More informationExamination paper for TMA4195 Mathematical Modeling
Department of Mathematical Sciences Examination paper for TMA4195 Mathematical Modeling Academic contact during examination: Harald Hanche-Olsen Phone: 73 59 35 25 Examination date: December 14, 2016 Examination
More informationWe will assume straight channels with simple geometries (prismatic channels) and steady state flow (in time).
56 Review Drag & Lift Laminar vs Turbulent Boundary Layer Turbulent boundary layers stay attached to bodies longer Narrower wake! Lower pressure drag! 8. Open-Channel Flow Pipe/duct flow closed, full,
More informationECE2262 Electric Circuit
ECE2262 Electric Circuit Chapter 7: FIRST AND SECOND-ORDER RL AND RC CIRCUITS Response to First-Order RL and RC Circuits Response to Second-Order RL and RC Circuits 1 2 7.1. Introduction 3 4 In dc steady
More information11.1 Mass Density. Fluids are materials that can flow, and they include both gases and liquids. The mass density of a liquid or gas is an
Chapter 11 Fluids 11.1 Mass Density Fluids are materials that can flow, and they include both gases and liquids. The mass density of a liquid or gas is an important factor that determines its behavior
More informationDifferential equations, comprehensive exam topics and sample questions
Differential equations, comprehensive exam topics and sample questions Topics covered ODE s: Chapters -5, 7, from Elementary Differential Equations by Edwards and Penney, 6th edition.. Exact solutions
More informationBoundary feedback control in networks of open channels
Available online at www.sciencedirect.com Automatica 39 (3) 1365 1376 www.elsevier.com/locate/automatica Boundary feedback control in networks of open channels J. de Halleux a, C. Prieur b, J.-M. Coron
More informationCHEMICAL ENGINEERING
GATE - CHEMICAL ENGINEERING SAMPLE THEORY BERNOULLI'S EQUATION HEAT TRANSFER: RADIATION CONTROLLER MODES For IIT-JAM, JNU, GATE, NET, NIMCET and Other Entrance Exams -C-8, Sheela Chowdhary Road, Talwandi,
More informationLecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30
289 Upcoming labs: Lecture 12 Lab 20: Internal model control (finish up) Lab 22: Force or Torque control experiments [Integrative] (2-3 sessions) Final Exam on 12/21/2015 (Monday)10:30-12:30 Today: Recap
More information