Assimilation de Données, Observabilité et Quantification d Incertitudes appliqués à l hydraulique à surface libre
|
|
- Mitchell Hampton
- 5 years ago
- Views:
Transcription
1 Assimilation de Données, Observabilité et Quantification d Incertitudes appliqués à l hydraulique à surface libre Pierre-Olivier Malaterre, Igor Gejadze, Hind Oubanas UMR G-EAU, Irstea, Montpellier, France Pour mieux affirmer ses missions, le Cemagref devient Irstea Colloque National d Assimilation 2 décembre 206
2 Plan Systèmes étudiés : canaux d'irrigation, fleuves Questions : problèmes inverses & incertitudes Modèles : Saint-Venant D-.5D non linéaire (SIC) Modèle Linéaire Tangent, Adjoint Filtre de Kalman, pour un canal d irrigation (X, Q p, C d ) FK, Evolution des Erreurs, Observabilité 4D-Var, données SWOT sur la Garonne 4D-Var, quantification d'incertitudes Conclusions
3 Irrigation canal Water for irrigation, industries and domestic uses Automatic gates SCADA systems
4 Control objectives Distribute water to users (offtaes) Minimize water losses Satisfy constraints (levels, discharges, etc) By mean of control actions (manual or automatic) on cross devices (gates, weirs, pumps) Automatic control = type of inverse problem
5 Inverse problems for controlled Irrigation Canals Known: Bathymetry, boundary conditions (up, mid, down) Observations: Water levels close to cross devices Unnowns / Uncertainties (active & passive): Offtae withdrawal or inflow (Q p ) Cross device discharge coefficients (C d ) Hydraulic state (Q, Z) (e.g. state feedbac ctrl) Friction coefficients Objective: water loss or theft, fault detection, model update for operation & maintenance, etc.
6 Surface Water and Ocean Topography (SWOT) mission (202) Scientific requirements Observable river width > 00 m Height accuracy Slope accuracy Width accuracy Data collection 0 cm over area > m².7 cm/m over area > m 5% of the evaluated river 90% of all ocean/continents within the orbit during 90% of the operational time. August 22, 206 hind.oubanas@irstea.fr 6
7 Inverse problems for rivers observed from space Known: River length, mean slope, mean annual discharge!? boundary condition (down)!? Observations: Water levels locally, or globally (SWOT swaths), local slope, width Unnowns / Uncertainties (active & passive): Boundary conditions (Q up ) Bathymetry, Friction (strong effect) Hydraulic state (Q, Z) outside observations windows (time and space) Tributaries inflows (Q p ) Objective: flood control, resource availability, navigation, water balance
8 Example (5-pool irrigation canal) y y 2 y 3 y4 y 5 y 6 y 7 y 8 y 9 y 0 Q p Q p2 Q p3 Q p4 Q p5 5 pools, 5 unnown withdrawals Q p, 0 water level measurements Y
9 SIC² : Simulation and Integration for Controls of Canals - D/.5D hydraulic model (Multiple beds + Storage areas) - Flow dynamics of rivers - Irrigation canals, drainage networ, etc. - Based on the full Saint-Venant equations - TLM & Adjoint (Tapenade) Gejadze & Malaterre 206 August 22, 206 hind.oubanas@irstea.fr 9
10 SIC² : Simulation and Integration for Controls of Canals Regular sections : t A + x Q = Q L t Q + x Q 2 /A + ga x Z = gas f + C Q L v, S f = - D/.5D hydraulic model (Multiple beds + Storage areas) - Flow dynamics of rivers - Irrigation canals, drainage networ, etc. - Based on the full Saint-Venant equations - TLM & Adjoint (Tapenade) Q 2 K 2 S A 2, t [0, T] R4/3 Q(x, t): Discharge Z(x, t): Water level A x, t : Wetted cross sectional area v x, t = Q/A: Mean velocity R: Hydraulic radius K S : Stricler coefficient C : Lateral discharge coefficient Gejadze & Malaterre 206 August 22, 206 hind.oubanas@irstea.fr 0
11 SIC² : Simulation and Integration for Controls of Canals Regular sections : t A + x Q = Q L t Q + x Q 2 /A + ga x Z = gas f + C Q L v, S f = Singular sections : Q S3,i Q S3,i+ = 0, Q S3,i = F(Z S3,i, Z S3,i+, C d S 3,i ) - D/.5D hydraulic model (Multiple beds + Storage areas) - Flow dynamics of rivers - Irrigation canals, drainage networ, etc. - Based on the full Saint-Venant equations - TLM & Adjoint (Tapenade) Q 2 K 2 S A 2, t [0, T] R4/3 Q(x, t): Discharge Z(x, t): Water level A x, t : Wetted cross sectional area v x, t = Q/A: Mean velocity R: Hydraulic radius K S : Stricler coefficient C : Lateral discharge coefficient C d : Cross device coefficient S i, : th cross section of the reach i Gejadze & Malaterre 206 August 22, 206 hind.oubanas@irstea.fr
12 SIC² : Simulation and Integration for Controls of Canals Regular sections : t A + x Q = Q L t Q + x Q 2 /A + ga x Z = gas f + C Q L v, S f = Singular sections : Q S3,i Q S3,i+ = 0, Q S3,i = F(Z S3,i, Z S3,i+, C d S 3,i ) Boundary conditions: - D/.5D hydraulic model (Multiple beds + Storage areas) - Flow dynamics of rivers - Irrigation canals, drainage networ, etc. - Based on the full Saint-Venant equations - TLM & Adjoint (Tapenade) Q 2 K 2 S A 2, t [0, T] R4/3 Upstream nodes : Q t or Z(t) Downstream nodes : Q t = f(z t, p rc ) Q(x, t): Discharge Z(x, t): Water level A x, t : Wetted cross sectional area v x, t = Q/A: Mean velocity R: Hydraulic radius K S : Stricler coefficient C : Lateral discharge coefficient C d : Cross device coefficient S i, : th cross section of the reach i p rc : rating curve parameters Gejadze & Malaterre 206 August 22, 206 hind.oubanas@irstea.fr 2
13 LQG (Linear Quadratic Gaussian) Perturbations (w) Q p Canal Sorties contrôlées (z) Z av Commandes (u) LQG Mesures (y)
14 Design du contrôleur LQG Minimisation d'un critère J : Contraintes dynamiques : J y y Q y y u R u T y T N 2 0 ( ( ) *( ))..( ( ) *( )) ( ).. ( ) ) (. ) ( ) (. ) (. ) (. ) (. ) ( C x y u B u B B u A x x c c p p u : vecteur des ouvertures des 5 vannes u p : vecteur des débits aux 5 prises u c : vecteur des coefficients de débit aux 5 ouvrages en travers
15 Solution optimale Solution optimale (retour d état) : u( ) K ( ). x( ) LQ? K LQ matrice gain du contrôleur LQ, solution d une équation de Riccati
16 Observateur d'état et de perturbation Cx y y y L u u y y L u u y y L u B u B B u A x x c c c p p p c c p p ˆ ˆ ˆ).( ˆ ˆ ˆ).( ˆ ˆ ˆ).( ˆ. ˆ.. ˆ. ˆ L, L p, L c : par placement de pole (A-LC) pour avoir convergence de l'erreur de reconstruction de l état x vers 0 => Observateur de Luenberger (théorie de l'observabilité, théorie de la détectabilité, principe de séparation) L, L p, L c : Minimisation de l'erreur (variance) de reconstruction => Filtre de Kalman
17 Kalman Filter equations Q A AP P BU AX X t ~ ~ ~ ˆ ~ ˆ t t P C K I P CX Y K X X R C CP C P K ~ ~ ˆ ˆ ˆ ~ ~ ~ c p c p c p c p v u u x C y w u B u u x B B A u u x Initial linear model + noise with given covariance matrices Q and R Augmented state for the reconstruction of x, u p (Q p ) and u c (C d ) No observability requirement!!!
18 Uncertainty reduction (on the 5 Q p and 5 C d ) y y 2 y 3 y4 y 5 y 6 y 7 y 8 y 9 y 0 Cd Cd 2 Q p Q p2 Cd 3 Cd 4 Cd 5 Q p3 Q p4 Q p5 Is it possible to identify Q p and C d from limited water level measurements Y?
19 Scenario tested with a Kalman Filter Time step for the simulations is 300 seconds. T end = 48h The 5 offtaes are operated at time 6 h (-0.00 m3/s) and again at time 20 h (+0.00 m3/s) returning to the initial discharge values of m3/s All cross regulators eep the same gate positions (open loop), but the discharge coefficient of their gates is changed from 0.82 to 0.66 at time 30 mn, then to 0.72 at time 6 h and returning to initial value 0.82 at 24 h The initial state X 0 is taen using random variables For the matrix Q, we choose σ Q = 0. m 3 /s and σ Z = 0. m for the normal hydraulic states and σ QC = for the augmented states (Q p and C d ). For the matrix R, we choose σ R = 0.0 m There are 0 water level measurements (y, y 2,, y 0 ), or less!
20 State reconstruction Twin exp X (Q & Z) 4 2 X Q 0 nx= Time (h) X Z Time (h)
21 State reconstruction Twin exp - C d Cd Reconstruction of discharges coefficients Cd at cross devices Cd Cd2 Cd3 Cd4 Cd5 Good and fast convergence Despite wrong initial state X 0 (No additional noise) Time (h)
22 State reconstruction Twin exp - Q p Qp Reconstruction of discharges Qp at offtaes Qp Qp2 Qp3 Qp4 Qp5 Good and fast convergence Despite wrong initial state X 0 (No additional noise) Time (h)
23 No convergence towards the correct values (case with 9 measurements only: y 8 removed) Reconstruction of discharges Qp at offtaes wrong! Qp Qp2 Qp3 Qp4 Qp Reconstruction of discharges coefficients Cd at cross devices wrong! Cd Cd2 Cd3 Cd4 Cd5 0 0 Qp Cd wrong! Time (h) Time (h) Timing is good but not the values Very bad for Q p3 and Q p4 (the 2 influenced a lot by y 8 ), surprisingly good for C d (except C d4 for the same reasons)
24 Variances of Pb on X Z Variances of Pb on Qp Increase of the local variance (close to y 8 ) during time iterations: we now that the values are not good State index Variance Measurement points State index Variances of Pb on X Q Variances of Pb on X Z State index Variances of Pb on Cd Variances of Pb on Qp Variance State index Variances of Pb on X Q Variances of Pb on Cd State index State index
25 Convergence study o a a o t a t a a t a o t a a t a o t a a t a o t a a t b a a t b b a t a a K C K A C K A K C K X X A C K X X A X A X A C C X A C K X A X A B U C X A C C B U C X A C K X A X A B U C X A C X C K X A B U X A X C Y K B U X A X X C Y K X X X ) ( 0 lim lim 0 0 et si? a a a a E K C A A A C K A A C K A E E E E A C K A E We want that the estimation error converges toward 0 Convergence condition : If for > n this term is Schur =0 or
26 With the 0 measurements: Convergence of the reconstruction error -> 0 Max eig. (A-KCA) = < Then the error between the estimated state X a and the true state X t is converging towards 0 (within the linear assumptions) Poles of the Kalman Filter during iterations /T /T /T 0.9/T 0.8/T 0.8/T 0.7/T 0.6/T Poles 0.5/T of the 0.4/T Kalman Filter. Max = /T /T /T /T 0./T 0.7/T 0.3/T /T 0.4/T 0.5/T
27 Non convergence if only 9 measurements or less (e.g.: remove y 8 ) Max eigenvalue = Poles of the Kalman Filter during iterations /T /T /T 0.9/T 0.8/T 0.8/T 0.7/T 0.6/T 0.5/T Poles 0.4/T of the Kalman Filter. Max = /T /T /T /T 0./T 0.7/T 0.3/T /T 0.4/T 0.5/T
28 Questions to solve Does it depend only on A, C? Is-it possible to choose or analyze (Q,R) so that the resulting gain matrix K, leads to a Schur (A-KCA) matrix?
29 Luenberger Observer Kalman Filter Observability Observability C CA 2 CA O n CA O full ran = n K such that p, K such that poles(a-kc)=p A-KC Schur Convergence of estimation? Condition on A, C??? (Q,R) K such that? (Q,R) K such that A-KCA Schur Convergence de of l estimation
30 Detectability Detectability non observables states of O are stables K such that A-KC Schur Convergence of estimation Rm: Same property for A-KCA
31 Convergence Sufficient condition for convergence : Let (A, C) pair be detectable Let Q matrix factorizable in Q=ΓΓ T, such that pair (A, Γ) be stabilizable Let R matrix be definite positive Then the Riccati algebraic equation for the Kalman Gain has and only positive solution and the corresponding matrix (A-KC) is Schur Rm: depends on A, C and Q Def: (A,B) controllable if ran(b AB A 2 B A n- B) = n Def: stabilizable if non controllable poles are stable
32 Influence de Q sur l'observabilité Non Oui "Optimal"
33 Conclusions (cadre linéaire FK). Les 0 mesures de niveau permettent d identifier les 5 Q p et 5 C d 2. Le problème inverse est bien posé, et l erreur de reconstruction de l état analysé tend bien vers 0 3. Dès qu on enlève mesure on perd ces résultats : l erreur ne tends pas vers 0, certaines composantes du vecteur de contrôle ne sont pas bien reconstituées 4. Un test simple a priori (eig value) ou a posteriori (variance) permet de vérifier cette propriété qu on souhaite 5. L observabilité (A,C) serait suffisante, mais n est pas nécessaire. Et difficile à vérifier en grande dimension (nx 0 2, 0 3 ) 6. La détectabilité (A,C) + qq autres propriétés sont suffisantes 7. Les covariances des erreurs de modèle Q sont importantes dans ces conditions 8. Le cadre linéaire précédent propose des outils puissants mais a des limites (linéarité, dimension, sensibilité au bruit, etc.) -> 4D-Var pour la suite
34 Surface Water and Ocean Topography (SWOT) mission Scientific requirements Observable river width > 00 m Height accuracy Slope accuracy Width accuracy Data collection 0 cm over area > m².7 cm/m over area > m 5% of the evaluated river 90% of all ocean/continents within the orbit during 90% of the operational time. August 22, 206 hind.oubanas@irstea.fr 34
35 Data assimilation method Variational approach Observations: Y = Y t T + ε o Y, ε o ~N 0, O, R = E ε o ε o Bacground: U b = U t T + ε b U, ε b ~N 0, B, B = E ε b ε b Control vector: U U Classical cost function: Tihonov cost function: Iterative regularization: J U = 2 R 2 G U Y J U, α = 2 R 2 G U Y B 2 U U b 2 + α2 2 B 2 U U b, α > 0, Tihonov regularization parameter J W = 2 R 2 G U b + B 2W, U 0 Y, J W ~χ 2 (M) 2 U = U b + B 2W, U U : Control vector U b U : Bacground Y Y : Observation vector M : Observation space dimension G U Y : Nonlinear mapping operator B : Bacground covariance matrix R: Observation covariance matrix August 22,
36 Data assimilation method Variational approach Limited-memory Broyden Fletcher Goldfarb Shanno (L-BFGS) method : W i+ = W i + β i H i J W i W, W 0 = 0 Gradient of the cost function Adjoint Model J W i W = B 2 G U b + B 2W, U 0 R G U b + B 2W, U 0 Y Automatic differentiation TAPENADE (INRIA) (Gejadze & Malaterre 206) August 22, 206 hind.oubanas@irstea.fr 36
37 Experimental framewor Study area / period Garonne River France Downstream reach : 50 m Mean width : 70m Mean slope : 28cm/m T c ~24h August 22, 206 hind.oubanas@irstea.fr 37
38 Experiment () : Estimation of Q, given K S and Z b Estimation of upstream discharge Q assuming nown the bed level Z b and the friction coefficient K S, investigating the influence of the SWOT temporal frequency. Identical twin experiments framewor Assimilation of water surface elevation Z observations Observations error σ = 0 cm Observations time period: from day to 5 days Observation spatial sampling: each 0 m The first guess on discharge is taen as the mean annual value Bathymetry and friction assumed nown Sequential version (DA sub-window: 75-day period) rrmse(q) = T T 0 Q estimate t Q true t 2 dt /2 August 22, 206 hind.oubanas@irstea.fr 38
39 Experiment () : Estimation of Q, given K S and Z b Discharge hydrograph at Tonneins from 0/0/200 to 3/05/200 (a) -day, (b) 2-day, (c) 4-day, (d) 5-day obervation period (a) (b) (c) (d) Q rrmse 2.% 9.5% 2.9% 8.2% August 22, 206 hind.oubanas@irstea.fr 39
40 Experiment () : Estimation of Q, given K S and Z b Nyquist sampling theorem The sampling frequency should be at least twice the highest frequency contained in the signal. Discharge hydrograph at Tonneins from 0/0/200 to 3/05/200 (a) -day, (b) 2-day, (c) 4-day, (d) 5-day obervation period (a) (b) (c) (d) Q rrmse 2.% 9.5% 2.9% 8.2% August 22, 206 hind.oubanas@irstea.fr 40
41 Experiment (2) : Simultaneous estimation of Q and K S, given exact Z b Estimation of upstream discharge Q under uncertainty in the friction coefficient K S, assuming nown the bed level Z b First guess on the friction coefficient is taen as a 20% error of the mean value Bathymetry assumed nown Sequential version (DA sub-window: 75-day period) August 22, 206 hind.oubanas@irstea.fr 4
42 Experiment (2) : Simultaneous estimation of Q and K S, given exact Z b A-I A-II Q rrmse 2.9% 2.6% K S rrmse 20.4% 3.4% A-I : Estimation of Q solely using the first guess on K S. A-II : Estimation of Q and K S. August 22, 206 hind.oubanas@irstea.fr 42
43 Experiment (2) : Simultaneous estimation of Q and K S, given exact Z b A-I A-II Q rrmse 2.9% 2.6% K S rrmse 20.4% 3.4% A-I : Estimation of Q solely using the first guess on K S. A-II : Estimation of Q and K S. Experiment (3) : Simultaneous estimation of Q and Z b, given exact K S Estimation of upstream discharge Q under uncertainty in the bed level Z b, assuming nown the friction coefficient K S First guess on Bed level is derived from the perturbed steady flow WSE Friction coefficient and cross-sections shape assumed nown Sequential version (DA sub-window: 75-day period) August 22, 206 hind.oubanas@irstea.fr 43
44 Experiment (2) : Simultaneous estimation of Q and K S, given exact Z b A-I A-II Q rrmse 2.9% 2.6% K S rrmse 20.4% 3.4% A-I : Estimation of Q solely using the first guess on K S. A-II : Estimation of Q and K S. Experiment (3) : Simultaneous estimation of Q and Z b, given exact K S B-I B-II Q rrmse 50% 3.8% Z b rrmse 5.7% 4.9% B-I : Estimation of Q solely using the first guess on Z b. B-II : Estimation of Q and Z b. August 22, 206 hind.oubanas@irstea.fr 44
45 Experiment (4) : Simultaneous estimation of Q, K S and Z b Estimation of upstream discharge Q, under uncertainty in the bed level Z b and the friction coefficient K S Identical twin experiments framewor Assimilation of water surface elevation Z observations Observations error σ = 0 cm Observations time period: day (up to 5 days) Observation spatial sampling: each 0 m The first guess on discharge is taen as the mean annual value First guess on the friction coefficient is taen as a 20% error of the mean value First guess on Bed level is derived from the perturbed steady flow WSE Cross-sections shape assumed nown Sequential version (DA sub-window: 75-day period) August 22, 206 hind.oubanas@irstea.fr 45
46 Experiment (4) : Simultaneous estimation of Q, K S and Z b C-I: Estimation of Q solely using the first guess on K S and Z b C-II: Estimation of Q, K S and Z b C-III: Estimation of Q and Z b, using the first guess on K S rrmse C-I C-II C-III Q 40.5% 7.% 5.% K S 3% 24.4% 3% Z b 5.7% 4.7% 4.5% November 2, 206 PhD Irstea/CLS 46
47 Experiment (4) : Simultaneous estimation of Q, K S and Z b C-I: Estimation of Q solely using the first guess on K S and Z b C-II: Estimation of Q, K S and Z b C-III: Estimation of Q and Z b, using the first guess on K S Equifinality issue!! rrmse C-I C-II C-III Q 40.5% 7.% 5. K S 3% 24.4% 3% Z b 5.7% 4.7% 4.5% November 2, 206 PhD Irstea/CLS 47
48 Control set design: notations
49 Control set design: variational DA
50 Control set design: goal-function covariance
51 Control set design: partial control case
52 Control set design: partial control covariance
53 Control set design: implementation
54 Control set design: algorithm
55 Control set design: test configuration
56 Control set design: NA results, case A
57 Control set design: NA results, case B
58 Conclusion 2 (cadre variationnel). Le concept de choix du vecteur de contrôle a été présenté. Ce concept est utile pour les modèles ayant des incertitudes nombreuses, complexes parmi leurs entrées. En particulier, les erreurs de modèle peuvent être incluses dans ces entrées du modèle. 2. La méthode présentée permet de quantifier la performance obtenue pour toute combinaison d entrées actives (sous-partie du vecteur des entrées), permettant de faire apparaitre les sous-vecteurs suffisants. Le choix entre ces options de sousvecteurs suffisants peut ensuite être fait selon des considérations de solvabilité et de robustesse. 3. La méthode est une généralisation de l approche variationnelle classique de Quantification ou Réduction des Incertitudes. Elle n utilise pas de formalisme matriciel et est donc adaptée aux modèles de grande dimension. 4. La méthode a été appliquée dans le domaine de l hydraulique à surface libre et a prouvé son intérêt. Les résultats obtenus illustrent et justifient toutes les conclusions que nous avions tirées de nos tests d AD avec SIC (et son adjoint). 5. La même approche peut être étendue au cas ou les entrées (vecteur de contrôle) sont divisées entre une part active, passive et de nuisance. Des développements peuvent aussi être proposés pour obtenir une performance globale plutôt que locale.
59 Thans! Questions? Thans for your attention
Outils de Recherche Opérationnelle en Génie MTH Astuce de modélisation en Programmation Linéaire
Outils de Recherche Opérationnelle en Génie MTH 8414 Astuce de modélisation en Programmation Linéaire Résumé Les problèmes ne se présentent pas toujours sous une forme qui soit naturellement linéaire.
More informationApprentissage automatique Machine à vecteurs de support - motivation
Apprentissage automatique Machine à vecteurs de support - motivation RÉGRESSION À NOYAU régression à noyau Algorithme de régression à noyau entraînement : prédiction : a = (K + λi N ) 1 t. y(x) =k(x) T
More informationApprentissage automatique Méthodes à noyaux - motivation
Apprentissage automatique Méthodes à noyaux - motivation MODÉLISATION NON-LINÉAIRE prédicteur non-linéaire On a vu plusieurs algorithmes qui produisent des modèles linéaires (régression ou classification)
More informationA MULTIVARIABLE APPROACH FOR THE COMMAND OF CANAL DE PROVENCE AIX NORD WATER SUPPLY SUBSYSTEM
A MULTIVARIABLE APPROACH FOR THE COMMAND OF CANAL DE PROVENCE AIX NORD WATER SUPPLY SUBSYSTEM Yann Viala 1 Pierre-Olivier Malaterre 2 Jean-Luc Deltour 1 Franck Sanfilippo 1 Jacques Sau 3 ABSTRACT The Canal
More informationDETERMINING HIGH VOLTAGE CABLE CONDUCTOR TEMPERATURES. Guy Van der Veken. Euromold, Belgium. INVESTIGATIONS. INTRODUCTION.
DETERMINING HIGH VOLTAGE CABLE CONDUCTOR TEMPERATURES. Guy Van der Veken. Euromold, Belgium. INTRODUCTION. INVESTIGATIONS. Type tests on MV cable accessories are described in CENELEC HD68 and HD69 documents.
More informationThe epsilon method: analysis of seepage beneath an impervious dam with sheet pile on a layered soil
The epsilon method: analysis of seepage beneath an impervious dam with sheet pile on a layered soil Zheng-yi Feng and Jonathan T.H. Wu 59 Abstract: An approximate solution method, referred to as the epsilon
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 informationAssimilation des Observations et Traitement des Incertitudes en Météorologie
Assimilation des Observations et Traitement des Incertitudes en Météorologie Olivier Talagrand Laboratoire de Météorologie Dynamique, Paris 4èmes Rencontres Météo/MathAppli Météo-France, Toulouse, 25 Mars
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationApprentissage automatique Classification linéaire - fonction discriminante
Apprentissage automatique Classification linéaire - fonction discriminante TYPES D APPRENTISSAGE apprentissage supervisé, classification, régression L apprentissage supervisé est lorsqu on a une cible
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 informationThe use of L-moments for regionalizing flow records in the Rio Uruguai basin: a case study
Regionalization in Ifylwltm (Proceedings of the Ljubljana Symposium, April 1990). IAHS Publ. no. 191, 1990. The use of L-moments for regionalizing flow records in the Rio Uruguai basin: a case study ROBM
More informationPoisson s ratio effect of slope stability calculations
Poisson s ratio effect of slope stability calculations Murray Fredlund, & Robert Thode SoilVision Systems Ltd., Saskatoon, SK, Canada ABSTRACT This paper presents the results of a study on the effect of
More informationOptimisation par réduction d incertitudes : application à la recherche d idéotypes
: application à la recherche d idéotypes Victor Picheny 1, D. Da Silva et E. Costes 2 Rencontres du réseau Mexico, Toulouse 23 mai 2014 1. INRA MIAT 2. INRA GAP Plan de l exposé 1 Introduction : recherche
More informationContaminant Isolation By Cutoff Walls: Reconsideration Of Mass Fluxes
GeoHalifax29/GéoHalifax29 Contaminant Isolation By Cutoff Walls: Reconsideration Of Mass Fluxes Christopher. Neville S.S. Papadopulos & ssociates, Inc., Waterloo, Ontario BSTRCT Cutoff alls are used frequently
More informationK. FUJITA INTRODUCTION. Dr., Managing Director of Hazama-Gumi, Ltd. K. UEDA. Deputy Director, Institute of Technology, Hazama-Gumi, Ltd. M.
A METHOD TO PREDICT THE LOAD-DISPLACEMENT RELATIONSHIP OF GROUND ANCHORS Modèle pour calculer la relation charge-déplacement des ancrages dans les sols by K. FUJITA Dr., Managing Director of Hazama-Gumi,
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationMACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS
MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS M.L. Kavvas and L.Liang UCD J.Amorocho Hydraulics Laboratory University of California, Davis, CA 95616, USA Uncertainties
More informationKato s inequality when u is a measure. L inégalité de Kato lorsque u est une mesure
Kato s inequality when u is a measure L inégalité de Kato lorsque u est une mesure Haïm Brezis a,b, Augusto C. Ponce a,b, a Laboratoire Jacques-Louis Lions, Université Pierre et Marie Curie, BC 187, 4
More informationBest linear unbiased prediction when error vector is correlated with other random vectors in the model
Best linear unbiased prediction when error vector is correlated with other random vectors in the model L.R. Schaeffer, C.R. Henderson To cite this version: L.R. Schaeffer, C.R. Henderson. Best linear unbiased
More informationValidation of an adjoint method for compressible channel flow sensitivities
9 ème Congrès Français de Mécanique Marseille, -8 août 9 Validation of an method for compressible channel flow sensitivities L. MORE-GABARRO, P. CAHALIFAUD, C. AIRIAU Université de oulouse ; CNRS, INP,
More information= 0 otherwise. Eu(n) = 0 and Eu(n)u(m) = δ n m
A-AE 567 Final Homework Spring 212 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You
More informationThe measurement and description of rill erosion
The hydrology of areas of low precipitation L'hydrologie des régions à faibles précipitations (Proceedings of the Canberra Symposium, December 1979; Actes du Colloque de Canberra, décembre 1979): IAHS-AISH
More informationContent. Content. Introduction. T. Chateau. Computer Vision. Introduction. Outil projectif permettant l acquisition d une scène 3D sur un plan 2D
Content Modèle de caméra T Chateau Lamea/Gravir/ComSee, Blaie Pacal Univerit Computer Viion 2 Content La projection perpective Changement de repère objet/caméra Changement de repère caméra/image Changement
More informationx(n + 1) = Ax(n) and y(n) = Cx(n) + 2v(n) and C = x(0) = ξ 1 ξ 2 Ex(0)x(0) = I
A-AE 567 Final Homework Spring 213 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You
More informationOPTIMAL CONTROL AND ESTIMATION
OPTIMAL CONTROL AND ESTIMATION Robert F. Stengel Department of Mechanical and Aerospace Engineering Princeton University, Princeton, New Jersey DOVER PUBLICATIONS, INC. New York CONTENTS 1. INTRODUCTION
More informationLecture 10: River Channels
GEOG415 Lecture 10: River Channels 10-1 Importance of channel characteristics Prediction of flow was the sole purpose of hydrology, and still is a very important aspect of hydrology. - Water balance gives
More informationControl Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation
Lecture 9: State Feedback and s [IFAC PB Ch 9] State Feedback s Disturbance Estimation & Integral Action Control Design Many factors to consider, for example: Attenuation of load disturbances Reduction
More information( ) 2 ( kg) ( 9.80 m/s 2
Chapitre 1 Charges et champs électriques [9 au 1 mai] DEVOIR : 1.78, 1.84, 1.56, 1.90, 1.71 1.1. Charge électrique et structure de la matière À lire rapidement. Concepts déjà familiers. 1.. Conducteurs,
More informationSediment yield and availability for two reservoir drainage basins in central Luzon, Philippines
Sediment Budgets (Proceedings of the Porto Alegre Symposium, December 1988). IAHS Publ. no. 174, 1988. Sediment yield and availability for two reservoir drainage basins in central Luzon, Philippines SUE
More informationOptimal synthesis of sensor networks
Université de Liège Faculté des Sciences Appliquées Département de Chimie Appliquée Laboratoire d Analyse et de Synthèse des Systèmes Chimiques Optimal synthesis of sensor networks Carine Gerkens Thèse
More information4DEnVar. Four-Dimensional Ensemble-Variational Data Assimilation. Colloque National sur l'assimilation de données
Four-Dimensional Ensemble-Variational Data Assimilation 4DEnVar Colloque National sur l'assimilation de données Andrew Lorenc, Toulouse France. 1-3 décembre 2014 Crown copyright Met Office 4DEnVar: Topics
More informationApport altimétrie pour modélisation hydrologique des grands fleuves
Apport altimétrie pour modélisation hydrologique des grands fleuves C. Emery, V. Pedinotti, S. Biancamaria, A. Boone, A. Paris, S. Calmant, S. Ricci, P.-A. Garambois, B. Decharme, Workshop Niger, 14-15
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationarxiv:cs/ v1 [cs.dm] 21 Apr 2005
arxiv:cs/0504090v1 [cs.dm] 21 Apr 2005 Abstract Discrete Morse Theory for free chain complexes Théorie de Morse pour des complexes de chaines libres Dmitry N. Kozlov Eidgenössische Technische Hochschule,
More informationMesurer des déplacements ET des contraintes par correlation mécanique d images
Mesurer des déplacements ET des contraintes par correlation mécanique d images J. RÉTHORÉ a, A. LEYGUE a, M. CORET a, L. STAINIER a, E. VERRON a a. Institut de Recherche en Génie Civil et Mécanique (GeM)
More informationSTUDY of DUALITY IN BOND GRAPH MODELS. APPLICATION TO CONTROL LAW SYNTHESIS
STUDY of DUALITY IN BOND GRAPH MODELS. APPLICATION TO CONTROL LAW SYNTHESIS Stefan Lichiardopol To cite this version: Stefan Lichiardopol. STUDY of DUALITY IN BOND GRAPH MODELS. APPLICATION TO CON- TROL
More informationFrom Newton s Fluxions to Virtual Microscopes
From Newton s Fluxions to Virtual Microscopes Jacques Bair and Valerie Henry J.Bair@ulg.ac.be V.Henry@ulg.ac.be Key words : Fluxions ; tangent lines ; hyperreal numbers ; dual numbers ; virtual microscopes
More informationExercise sheet n Compute the eigenvalues and the eigenvectors of the following matrices. C =
L2 - UE MAT334 Exercise sheet n 7 Eigenvalues and eigenvectors 1. Compute the eigenvalues and the eigenvectors of the following matrices. 1 1 1 2 3 4 4 1 4 B = 1 1 1 1 1 1 1 1 1 C = Which of the previous
More informationON 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 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 informationAERT 2013 [CA'NTI 19] ALGORITHMES DE COMMANDE NUMÉRIQUE OPTIMALE DES TURBINES ÉOLIENNES
AER 2013 [CA'NI 19] ALGORIHMES DE COMMANDE NUMÉRIQUE OPIMALE DES URBINES ÉOLIENNES Eng. Raluca MAEESCU Dr.Eng Andreea PINEA Prof.Dr.Eng. Nikolai CHRISOV Prof.Dr.Eng. Dan SEFANOIU Eng. Raluca MAEESCU CONEN
More informationDATA ASSIMILATION FOR FLOOD FORECASTING
DATA ASSIMILATION FOR FLOOD FORECASTING Arnold Heemin Delft University of Technology 09/16/14 1 Data assimilation is the incorporation of measurement into a numerical model to improve the model results
More informationOn variational data assimilation for 1D and 2D fluvial hydraulics
On variational data assimilation for D and D fluvial hydraulics I. Gejadze, M. Honnorat, FX Le Dimet, and J. Monnier LJK - MOISE project-team, Grenoble, France. Contact: Jerome.Monnier@imag.fr Civil Engineering
More informationADAPTIVE SLIDING-MODE CONTROL: SIMPLEX-TYPE METHOD. Ta-Tau Chen
ADAPTIVE SLIDING-MODE CONTROL: SIMPLEX-TYPE METHOD Ta-Tau Chen Department of Electrical Engineering Kun Shan University, Tainan, Taiwan, R.O.C. E-mail: cdd0221@yahoo.com.tw ICETI 2014 AA1004_SCI No. 15-CSME-17,
More informationRelative Merits of 4D-Var and Ensemble Kalman Filter
Relative Merits of 4D-Var and Ensemble Kalman Filter Andrew Lorenc Met Office, Exeter International summer school on Atmospheric and Oceanic Sciences (ISSAOS) "Atmospheric Data Assimilation". August 29
More informationOn a multivariate implementation of the Gibbs sampler
Note On a multivariate implementation of the Gibbs sampler LA García-Cortés, D Sorensen* National Institute of Animal Science, Research Center Foulum, PB 39, DK-8830 Tjele, Denmark (Received 2 August 1995;
More informationTHÈSE. présentée et soutenue publiquement le 04/10/2013 pour l obtention du. Doctorat de l Université de Lorraine. spécialité automatique par
Ecole doctorale IAEM Lorraine Département de formation doctorale en automatique UFR ESSTIN Identification de Systèmes Dynamiques Hybrides : géométrie, parcimonie, et non-linéarités Hybrid Dynamical System
More informationANNALES. FLORENT BALACHEFF, ERAN MAKOVER, HUGO PARLIER Systole growth for finite area hyperbolic surfaces
ANNALES DE LA FACULTÉ DES SCIENCES Mathématiques FLORENT BALACHEFF, ERAN MAKOVER, HUGO PARLIER Systole growth for finite area hyperbolic surfaces Tome XXIII, n o 1 (2014), p. 175-180.
More informationOn the direct kinematics of planar parallel manipulators: special architectures and number of solutions
On the direct kinematics of planar parallel manipulators: special architectures and number of solutions by Clément M. Gosselin and Jean-Pierre Merlet Département de Génie Mécanique Université Laval Ste-Foy,
More informationDeveloping rules of thumb for groundwater modelling in large open pit mine design
Developing rules of thumb for groundwater modelling in large open pit mine design Jim Hazzard, Branko Damjanac, Christine Detournay & Loren Lorig Itasca Consulting Group, Minneapolis, MN, USA ABSTRACT
More informationSelection on selected records
Selection on selected records B. GOFFINET I.N.R.A., Laboratoire de Biometrie, Centre de Recherches de Toulouse, chemin de Borde-Rouge, F 31320 Castanet- Tolosan Summary. The problem of selecting individuals
More informationAnalysis of the quality of suspended sediment data
Sediment Budgets (Proceedings of the Porto Alegre Symposium, December 1988). IAHS Publ. no. 174, 1988. Analysis of the quality of suspended sediment data F. R. SEMMELMANN & A. E. LANNA Institute of Hydraulic
More informationA Stochastic Approach For The Range Evaluation
A Stochastic Approach For The Range Evaluation Andrei Banciu To cite this version: Andrei Banciu. A Stochastic Approach For The Range Evaluation. Signal and Image processing. Université Rennes 1, 2012.
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
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 informationNonlinear error dynamics for cycled data assimilation methods
Nonlinear error dynamics for cycled data assimilation methods A J F Moodey 1, A S Lawless 1,2, P J van Leeuwen 2, R W E Potthast 1,3 1 Department of Mathematics and Statistics, University of Reading, UK.
More informationContrôle multi-objectifs d ordre réduit
Contrôle multi-objectifs d ordre réduit Christian Fischer To cite this version: Christian Fischer. Contrôle multi-objectifs d ordre réduit. Autre [cs.oh]. École Nationale Supérieure des Mines de Paris,
More informationSteady State Kalman Filter
Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:
More informationVariational Data Assimilation Current Status
Variational Data Assimilation Current Status Eĺıas Valur Hólm with contributions from Mike Fisher and Yannick Trémolet ECMWF NORDITA Solar and stellar dynamos and cycles October 2009 Eĺıas Valur Hólm (ECMWF)
More informationA DIFFERENT APPROACH TO MULTIPLE CORRESPONDENCE ANALYSIS (MCA) THAN THAT OF SPECIFIC MCA. Odysseas E. MOSCHIDIS 1
Math. Sci. hum / Mathematics and Social Sciences 47 e année, n 86, 009), p. 77-88) A DIFFERENT APPROACH TO MULTIPLE CORRESPONDENCE ANALYSIS MCA) THAN THAT OF SPECIFIC MCA Odysseas E. MOSCHIDIS RÉSUMÉ Un
More informationModélisation & simulation de la génération de champs magnetiques par des écoulements de métaux liquides. Wietze Herreman
Modélisation & simulation de la génération de champs magnetiques par des écoulements de métaux liquides Wietze Herreman 19ième Colloque Alain Bouyssy!!"#$%&'()*+(,#-*.#*+( )/01+"2(!!!!!!! Origine des champs
More informationFeed-Forward Control of Open Channel Flow Using Differential Flatness
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY 1 Feed-Forward Control of Open Channel Flow Using Differential Flatness Tarek S. Rabbani, Florent Di Meglio, Xavier Litrico, and Alexandre M. Bayen Abstract
More information(Extended) Kalman Filter
(Extended) Kalman Filter Brian Hunt 7 June 2013 Goals of Data Assimilation (DA) Estimate the state of a system based on both current and all past observations of the system, using a model for the system
More informationAgeostrophic instabilities of a front in a stratified rotating fluid
8 ème Congrès Français de Mécanique Grenoble, 27-3 août 27 Ageostrophic instabilities of a front in a stratified rotating fluid J. Gula, R. Plougonven & V. Zeitlin Laboratoire de Météorologie Dynamique
More informationTOLERABLE HAZARD RATE FOR FUNCTION WITH INDEPENDENT SAFETY BARRIER ACTING AS FAILURE DETECTION AND NEGATION MECHANISM
e Congrès de maîtrise des risques et de sûreté de fonctionnement - Saint-Malo -3 octobre 6 TOLERABLE HAZARD RATE FOR FUNCTION WITH INDEPENDENT SAFETY BARRIER ACTING AS FAILURE DETECTION AND NEGATION MECHANISM
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationFeed-Forward Control of Open Channel Flow Using Differential Flatness
IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY, VOL. 18, NO. 1, JANUARY 2010 213 Feed-Forward Control of Open Channel Flow Using Differential Flatness Tarek S. Rabbani, Florent Di Meglio, Xavier Litrico,
More informationSediment storage requirements for reservoirs
Challenges in African Hydrology and Water Resources (Proceedings of the Harare Symposium, July 1984). IAHS Publ. no. 144. Sediment storage requirements for reservoirs INTRODUCTION T, C, KABELL The Hydrological
More informationEstimation of the transmissmty of the Santiago aquifer, Chile, using different geostatistical methods
Groundwater Management: Quantity and Quality (Proceedings of the Benidorm Symposium, October 1989). IAHS Publ. no. 188,1989. Estimation of the transmissmty of the Santiago aquifer, Chile, using different
More informationCRITICAL SHEAR STRESSES FOR EROSION AND DEPOSITION OF FINE SUSPENDED SEDIMENTS IN THE FRASER RIVER DOE FRAP
CRITICAL SHEAR STRESSES FOR EROSION AND DEPOSITION OF FINE SUSPENDED SEDIMENTS IN THE FRASER RIVER DOE FRAP 1994-13 Prepared for : Environment Canada Environmental Conservation Branch Science Division
More informationSome elements for improving interpretation of concrete electrical resistivity
Some elements for improving interpretation of concrete electrical resistivity Jean-François LATASTE 1, Stéphane LAURENS 2, Andrzej MOCZKO 3 1 Université Bordeaux 1, Talence, France, jf.lataste@ghymac.u-bordeaux1.fr
More informationNonlinear Observers. Jaime A. Moreno. Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México
Nonlinear Observers Jaime A. Moreno JMorenoP@ii.unam.mx Eléctrica y Computación Instituto de Ingeniería Universidad Nacional Autónoma de México XVI Congreso Latinoamericano de Control Automático October
More informationBasis Function Selection Criterion for Modal Monitoring of Non Stationary Systems ABSTRACT RÉSUMÉ
Basis Function Selection Criterion for Modal Monitoring of Non Stationary Systems Li W. 1, Vu V. H. 1, Liu Z. 1, Thomas M. 1 and Hazel B. 2 Zhaoheng.Liu@etsmtl.ca, Marc.Thomas@etsmtl.ca 1 Dynamo laboratory,
More informationApplication of computational fluid dynamics to spray drying
Application of computational fluid dynamics to spray drying Simon Lo To cite this version: Simon Lo. Application of computational fluid dynamics to spray drying. Le Lait, INRA Editions, 2005, 85 (4-5),
More informationUnconditionally stable scheme for Riccati equation
ESAIM: Proceedings, Vol. 8, 2, 39-52 Contrôle des systèmes gouvernés par des équations aux dérivées partielles http://www.emath.fr/proc/vol.8/ Unconditionally stable scheme for Riccati equation François
More informationGate. Dam. Bank. z(1) u(2) pool. z(2) u(3) u(1) z(4) z(3) u(4) z(5) u(5) w(1) (2) (3) w(4) w(5)
` controller design for a high-order 5-pool irrigation canal system Pierre-Olivier Malaterre, Mustafa Khammash 2 Keywords: control of irrigation canals, civil structures, ` control, robustness, integrator,
More informationRunup and uncertainty quantification: sensitivity analysis via ANOVA decomposition
Runup and uncertainty quantification: sensitivity analysis via ANOVA decomposition Mario Ricchiuto, Pietro Marco Congedo, Argiris I. Delis To cite this version: Mario Ricchiuto, Pietro Marco Congedo, Argiris
More informationAssimilation de données dans Polyphemus
Assimilation de données dans Polyphemus Action ADOQA Lin Wu et Vivien Mallet Lin.Wu@inria.fr Projet CLIME (INRIA / ENPC) Assimilation de données dans Polyphemus p. 1 Introduction System structure of Polyphemus
More informationLecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004
MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical
More informationLa question posée (en français, avec des mots justes ; pour un calcul, l'objectif doit être clairement écrit formellement)
Exercise : You have to make one ton of mayonnaise sauce using 95 % oil, 2.5 % egg yolk, 2.5 % vinegar. What is the minimum energy that you have to spend? Calculation for mayonnaise Hervé 4th October 2013
More informationFlatness-based control of an irrigation canal using SCADA
Flatness-based control of an irrigation canal using SCADA Tarek Rabbani, Simon Munier, David Dorchies, Pierre-Olivier Malaterre, Alexandre Bayen and Xavier Litrico With a population of more than six billion
More informationThèse de Doctorat de L'Université Paris-Saclay. L'Université Paris-Sud. Inria Saclay Ile-de-France
NNT : 2016SACLS459 Thèse de Doctorat de L'Université Paris-Saclay préparée à L'Université Paris-Sud au sein de Inria Saclay Ile-de-France ÉCOLE DOCTORALE N 580 Sciences et technologies de l'information
More informationème Congrès annuel, Section technique, ATPPC th Annual Meeting, PAPTAC
2000-86 ème Congrès annuel, Section technique, ATPPC 2000-86th Annual Meeting, PAPTAC Use of components of formation for predicting print quality and physical properties of newsprint Jean-Philippe Bernié,
More informationSnow mapping and hydrological forecasting by airborne Y-ray spectrometry in northern Sweden
Hydrological Applications of Remote Sensing and Remote Data Transmission (Proceedings of the Hamburg Symposium, August 1983). IAHS Publ. no. 145. Snow mapping and hydrological forecasting by airborne Y-ray
More informationTransition from transient Theis wells to steady Thiem wells
Mythological Sciences Journal des Sciences Hydrologiques, 43(6) December 1998 859 Transition from transient Theis wells to steady Thiem wells WILLEM J. ZAADNOORDIJK* /WACO, Consultants for Water and Environment,
More informationCRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN
Criteria for the choice of flood routing methods in natural... CRITERIA FOR THE CHOICE OF FLOOD ROUTING METHODS IN NATURAL CHANNELS WITH OVERBANK FLOWS Roger Moussa 1 Abstract: The classification of river
More information39.1 Gradually Varied Unsteady Flow
39.1 Gradually Varied Unsteady Flow Gradually varied unsteady low occurs when the low variables such as the low depth and velocity do not change rapidly in time and space. Such lows are very common in
More informationOn a multilinear character sum of Burgess
On a multilinear character sum of Burgess J. Bourgain IAS M. Chang UCR January 20, 2010 Abstract Let p be a sufficiently large prime and (L i ) 1 i n a nondegenerate system of linear forms in n variables
More informationAIR SYSTEMS MODELING AND CONTROL FOR TURBOCHARGED ENGINES
AIR SYSTEMS MODELING AND CONTROL FOR TURBOCHARGED ENGINES Philippe Moulin To cite this version: Philippe Moulin. AIR SYSTEMS MODELING AND CONTROL FOR TURBOCHARGED EN- GINES. Automatic. École Nationale
More informationLUC CHARROIS. E. COS M E 2, M. D U M O N T 1, M. L A FAY S S E 1, S. M O R I N 1, Q. L I B O I S 2, G. P I C A R D 2 a n d L. 1,2
4 T H WOR KS HOP R E MOT E S E N S I N G A N D MODE L IN G OF S UR FACES P RO PE RT IES Towards the assimilation of MODIS reflectances into the detailled snowpack model SURFEX/ISBA- Crocus LUC CHARROIS
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationEstimation of monthly river runoff data on the basis of satellite imagery
Hydrological Applications of Remote Sensing and Remote Data Transmission (Proceedings of the Hamburg Symposium, August 1983). IAHS Publ. no. 145. Estimation of monthly river runoff data on the basis of
More informationEvaluation of Muskingum - Cunge model. irrigation advance phase in Shahid
ICID 21 st International Congress on Irrigation and Drainage, 15-23 ICID 21 st October Congress, 211, Tehran, Tehran, October Iran 211 R.56.5/Poster/8 Evaluation of Muskingum - Cunge model with Sirmod
More informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
More informationEcoulements turbulents des eaux peu profondes Turbulent flows in shallow water
Ecoulements turbulents des eaux peu profondes Turbulent flows in shallow water SERGEY GAVRILYUK a, HENRI GOUIN b a. Université d Aix-Marseille & C.N.R.S. U.M.R. 6595, IUSTI, Project SMASH, 5 rue E. Fermi,
More informationA STRATEGY FOR FREE-VARIABLE TABLEAUX
A STRATEGY FOR FREE-VARIABLE TABLEAUX FOR VARIANTS OF QUANTIFIED MODAL LOGICS Virginie Thion Université de Paris-Sud, L.R.I. E-mail : Virginie.Thion@lri.fr May 2002 Abstract In previous papers [CMC00,
More informationIt s a Small World After All Calculus without s and s
It s a Small World After All Calculus without s and s Dan Sloughter Department of Mathematics Furman University November 18, 2004 Smallworld p1/39 L Hôpital s axiom Guillaume François Antoine Marquis de
More informationA Singularity-Free Method for the Time-Dependent Schrödinger Equation for Nonlinear Molecules
A Singularity-Free Method for the Time-Dependent Schrödinger Equation for Nonlinear Molecules A.D.Bandrauk a,, Hui-Zhong Lu a a Labo. de Chimie Théorique, Faculté des Sciences, Université de Sherbrooke,
More information