Optimal Control and Estimation:

Size: px
Start display at page:

Download "Optimal Control and Estimation:"

Transcription

1 Optimal Control and Estimation: An ICFD Retrospective N.K. Nichols University of Reading ICFD 25th Anniversary

2 With thanks to: CEGB / National Power SSEB BP Exploration British Gas Institute of Hydrology (CEH) HR Wallingford DRA RAE Bedford Met Office for support for 26 PhD students and 7 post-docs to work on this research

3 Outline 1. Introduction 2. Tidal Energy 3. Flight Control 4. Numerical Weather Prediction 5. Conclusions / Future

4 1. Introduction

5

6

7 Control problem: Find the input or control function that produces a given response from the system State estimation problem: Given measured output from system, find the state of the system at a specified time.

8 2. Tidal Energy

9

10 Control is discontinuous!

11 Results No control Energy = With control Energy =

12 Severn Barrage Model

13 Fluid conservation Momentum conservation Boundary conditions Periodic conditions Control constraints

14

15

16 3. Automatic Flight Control Systems

17

18

19

20

21

22 2

23

24 Bifurcation Diagrams Jet Fighter No control Robust decoupling control

25 4. Numerical Weather Prediction Met Office Website

26 Significant Properties: Very large number of unknowns ( ) Few observations ( ) System nonlinear unstable/chaotic Multi-scale dynamics Stochastic errors in model and observations

27 State Estimation Aim: Find the best estimate (analysis) of the true state of a system, consistent with both observations and the system dynamics given: Numerical prediction model Observations of the system (over time) Background state (prior) Estimates of the errors

28 State Estimation Aim: Find the best estimate (analysis) of the true state of a system, consistent with both observations and the system dynamics: State Numerical prediction model Observations of the system (over time) Background state (prior) Estimates of the errors

29 Best Unbiased Estimate min J ( x 0 ) = min 1 ( x 0 x b )T B 1 ( x 0 x b ) 2 n + ( H i [x i ] y io )T R i 1 ( H i [x i ] y io ) i= 0 subject to where x b - Background state (prior estimate) y i - Observations H i - Observation operator B - Background error covariance matrix R i - Observation error covariance matrix

30 Algorithm Uses Approximate Gauss-Newton method: Solves sequence of linearized least squares problems by inner gradient iteration procedure Finds gradients by adjoint integration Truncates inner loop iterations Uses approximate linear system models Theoretical convergence results obtained by reference to GaussNewton method (SIOPT).

31 New Research Find approximate linear system models using optimal reduced order modeling techniques from control theory to improve computational efficiency Test feasibility of approach: compare solutions using Low resolution linear model Optimal reduced order model for a simple shallow water flow model.

32 Norm of analysis error for 1-D SWE model Log Error Low Res Model of order = 200 Norm = vs Reduced Model of order = 200 Norm = Low Res Model of order = 200 Norm = vs Reduced Model of order = 80 Norm = Component of state Red (dotted) = Low Resolution Model Green (dashed) = Reduced Order Model

33 5. Conclusions / Future Many exciting open problems exist with valuable applications New research areas include: - PDE constrained optimization - Stochastic PDE control and estimation (data assimilation) - Model Order Reduction - Uncertainty in prediction

34 References A.S. Lawless, N.K. Nichols, C. Boess and A. Bunse-Gerstner, Using model reduction methods within incremental four-dimensional variational data assimilation, Monthly Weather Review, 136, 2008, S. Gratton, A.S. Lawless and N.K. Nichols, Approximate GaussNewton methods for nonlinear least squares problems, SIAM Journal on Optimization, 18, 2007, D.M. Littleboy and N.K. Nichols, Modal Coupling in Linear Control Systems Using Robust Eigenstructure Assignment, Linear Algebra and Applications, , 1998, L.P. Gibson, N.K. Nichols and D.M. Littleboy, Bifurcation analysis of eigenstructure assignment control in a simple nonlinear aircraft model, J. of Guidance, Control and Dynamics, Vol. 21, 1998, J. Kautsky, N.K. Nichols and P. Van Dooren, Robust pole assignment in linear state feedback, International J. of Control, 41, 1985, N.R.C. Birkett, N.K. Nichols and D.A.C. Nicol, Dynamic models for optimal control of tidal power schemes, ICE Tidal Power Symposium, Thomas Telford, London, 1987, N.R.C. Birkett and N.K. Nichols, Optimal control problems in tidal power generation, Industrial Numerical Analysis - Case Histories, (eds. C. Elliott and S. McKee), Clarendon Press, Oxford, 1986, N.R.C. Birkett, B.M. Count, N.K. Nichols and D.A.C. Nicol, Optimal control problems in tidal power generation - full estuary model, Applied Optimization Techniques in Energy Problems, (ed. Hj. Wacker), B.G. Teubner, Stuttgart, 1985, N.R.C. Birkett, B.M. Count and N.K. Nichols, Optimal control problems in tidal power, Water Power and Dam Construction, Jan. 1984,

Optimal state estimation for numerical weather prediction using reduced order models

Optimal state estimation for numerical weather prediction using reduced order models Optimal state estimation for numerical weather prediction using reduced order models A.S. Lawless 1, C. Boess 1, N.K. Nichols 1 and A. Bunse-Gerstner 2 1 School of Mathematical and Physical Sciences, University

More information

A comparison of the model error and state formulations of weak-constraint 4D-Var

A comparison of the model error and state formulations of weak-constraint 4D-Var A comparison of the model error and state formulations of weak-constraint 4D-Var A. El-Said, A.S. Lawless and N.K. Nichols School of Mathematical and Physical Sciences University of Reading Work supported

More information

Using Model Reduction Methods within Incremental Four-Dimensional Variational Data Assimilation

Using Model Reduction Methods within Incremental Four-Dimensional Variational Data Assimilation APRIL 2008 L A W L E S S E T A L. 1511 Using Model Reduction Methods within Incremental Four-Dimensional Variational Data Assimilation A. S. LAWLESS AND N. K. NICHOLS Department of Mathematics, University

More information

Mathematical Concepts of Data Assimilation

Mathematical Concepts of Data Assimilation Mathematical Concepts of Data Assimilation N.K. Nichols 1 Introduction Environmental systems can be realistically described by mathematical and numerical models of the system dynamics. These models can

More information

Variational assimilation Practical considerations. Amos S. Lawless

Variational assimilation Practical considerations. Amos S. Lawless Variational assimilation Practical considerations Amos S. Lawless a.s.lawless@reading.ac.uk 4D-Var problem ] [ ] [ 2 2 min i i i n i i T i i i b T b h h J y R y B,,, n i i i i f subject to Minimization

More information

Department of Mathematics and Statistics

Department of Mathematics and Statistics School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-01-17 15 August 01 Data assimilation with correlated observation errors: analysis accuracy with approximate

More information

The University of Reading

The University of Reading The University of Reading Radial Velocity Assimilation and Experiments with a Simple Shallow Water Model S.J. Rennie 2 and S.L. Dance 1,2 NUMERICAL ANALYSIS REPORT 1/2008 1 Department of Mathematics 2

More information

New Applications and Challenges In Data Assimilation

New Applications and Challenges In Data Assimilation New Applications and Challenges In Data Assimilation Met Office Nancy Nichols University of Reading 1. Observation Part Errors 1. Applications Coupled Ocean-Atmosphere Ensemble covariances for coupled

More information

DATA ASSIMILATION FOR FLOOD FORECASTING

DATA 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 information

Introduction to System Identification and Adaptive Control

Introduction to System Identification and Adaptive Control Introduction to System Identification and Adaptive Control A. Khaki Sedigh Control Systems Group Faculty of Electrical and Computer Engineering K. N. Toosi University of Technology May 2009 Introduction

More information

OPTIMAL ESTIMATION of DYNAMIC SYSTEMS

OPTIMAL ESTIMATION of DYNAMIC SYSTEMS CHAPMAN & HALL/CRC APPLIED MATHEMATICS -. AND NONLINEAR SCIENCE SERIES OPTIMAL ESTIMATION of DYNAMIC SYSTEMS John L Crassidis and John L. Junkins CHAPMAN & HALL/CRC A CRC Press Company Boca Raton London

More information

Uncertainty quantification for Wavefield Reconstruction Inversion

Uncertainty quantification for Wavefield Reconstruction Inversion Uncertainty quantification for Wavefield Reconstruction Inversion Zhilong Fang *, Chia Ying Lee, Curt Da Silva *, Felix J. Herrmann *, and Rachel Kuske * Seismic Laboratory for Imaging and Modeling (SLIM),

More information

Numerical Optimization of Partial Differential Equations

Numerical Optimization of Partial Differential Equations Numerical Optimization of Partial Differential Equations Bartosz Protas Department of Mathematics & Statistics McMaster University, Hamilton, Ontario, Canada URL: http://www.math.mcmaster.ca/bprotas Rencontres

More information

Active Stabilization of Unstable System Under Bounded Control with Application to Active Flutter Suppression Problem

Active Stabilization of Unstable System Under Bounded Control with Application to Active Flutter Suppression Problem Active Stabilization of Unstable System Under Bounded Control with Application to Active Flutter Suppression Problem Student: Supervisor: Second Reader: Konstantin Vikhorev Prof. Mikhail Goman Dr. Junsheng

More information

c 2007 Society for Industrial and Applied Mathematics

c 2007 Society for Industrial and Applied Mathematics SIAM J. OPTIM. Vol. 18, No. 1, pp. 106 13 c 007 Society for Industrial and Applied Mathematics APPROXIMATE GAUSS NEWTON METHODS FOR NONLINEAR LEAST SQUARES PROBLEMS S. GRATTON, A. S. LAWLESS, AND N. K.

More information

OPTIMAL CONTROL AND ESTIMATION

OPTIMAL 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 information

4. DATA ASSIMILATION FUNDAMENTALS

4. DATA ASSIMILATION FUNDAMENTALS 4. DATA ASSIMILATION FUNDAMENTALS... [the atmosphere] "is a chaotic system in which errors introduced into the system can grow with time... As a consequence, data assimilation is a struggle between chaotic

More information

Variational data assimilation

Variational data assimilation Background and methods NCEO, Dept. of Meteorology, Univ. of Reading 710 March 2018, Univ. of Reading Bayes' Theorem Bayes' Theorem p(x y) = posterior distribution = p(x) p(y x) p(y) prior distribution

More information

The Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance

The Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July

More information

Forecasting and data assimilation

Forecasting and data assimilation Supported by the National Science Foundation DMS Forecasting and data assimilation Outline Numerical models Kalman Filter Ensembles Douglas Nychka, Thomas Bengtsson, Chris Snyder Geophysical Statistics

More information

Inverse Problems and Optimal Design in Electricity and Magnetism

Inverse Problems and Optimal Design in Electricity and Magnetism Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini

More information

How to Validate Stochastic Finite Element Models from Uncertain Experimental Modal Data Yves Govers

How to Validate Stochastic Finite Element Models from Uncertain Experimental Modal Data Yves Govers How to Validate Stochastic Finite Element Models from Uncertain Experimental Modal Data Yves Govers Slide 1 Outline/ Motivation Validation of Finite Element Models on basis of modal data (eigenfrequencies

More information

A Comparative Study on Automatic Flight Control for small UAV

A Comparative Study on Automatic Flight Control for small UAV Proceedings of the 5 th International Conference of Control, Dynamic Systems, and Robotics (CDSR'18) Niagara Falls, Canada June 7 9, 18 Paper No. 13 DOI: 1.11159/cdsr18.13 A Comparative Study on Automatic

More information

Background and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience

Background and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Background and observation error covariances Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading

More information

ADJONT-BASED ANALYSIS OF OBSERVATION IMPACT ON TROPICAL CYCLONE INTENSITY FORECASTS

ADJONT-BASED ANALYSIS OF OBSERVATION IMPACT ON TROPICAL CYCLONE INTENSITY FORECASTS 7A.3 ADJONT-BASED ANALYSIS OF OBSERVATION IMPACT ON TROPICAL CYCLONE INTENSITY FORECASTS Brett T. Hoover* and Chris S. Velden Cooperative Institute for Meteorological Satellite Studies, Space Science and

More information

HADLEY CELL EXPANSION IN TODAY S CLIMATE AND PALEOCLIMATES

HADLEY CELL EXPANSION IN TODAY S CLIMATE AND PALEOCLIMATES HADLEY CELL EXPANSION IN TODAY S CLIMATE AND PALEOCLIMATES Bill Langford University Professor Emeritus Department of Mathematics and Statistics University of Guelph, Canada Presented to the BioM&S Symposium

More information

Operational Perspectives on Hydrologic Model Data Assimilation

Operational Perspectives on Hydrologic Model Data Assimilation Operational Perspectives on Hydrologic Model Data Assimilation Rob Hartman Hydrologist in Charge NOAA / National Weather Service California-Nevada River Forecast Center Sacramento, CA USA Outline Operational

More information

Feedback Control of Aerodynamic Flows

Feedback Control of Aerodynamic Flows 44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-843 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Feedback Control of Aerodynamic

More information

The Lifted Newton Method and Its Use in Optimization

The Lifted Newton Method and Its Use in Optimization The Lifted Newton Method and Its Use in Optimization Moritz Diehl Optimization in Engineering Center (OPTEC), K.U. Leuven, Belgium joint work with Jan Albersmeyer (U. Heidelberg) ENSIACET, Toulouse, February

More information

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev

SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,

More information

1 Number Systems and Errors 1

1 Number Systems and Errors 1 Contents 1 Number Systems and Errors 1 1.1 Introduction................................ 1 1.2 Number Representation and Base of Numbers............. 1 1.2.1 Normalized Floating-point Representation...........

More information

Introduction to Model Order Reduction

Introduction to Model Order Reduction Introduction to Model Order Reduction Lecture 1: Introduction and overview Henrik Sandberg Kin Cheong Sou Automatic Control Lab, KTH ACCESS Specialized Course Graduate level Ht 2010, period 1 1 Overview

More information

Optimal Control of Weakly Coupled Systems and Applications

Optimal Control of Weakly Coupled Systems and Applications Optimal Control of Weakly Coupled Systems and Applications HIGH ACCURACY TECHNIQUES Z. Gajić, M-T. Lim, D. Skatarić W-C. Su, and V. Kecman Taylor & Francis (CRC Press, Dekker) 2008 Preface This book is

More information

A NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION

A NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION A NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION ALIREZA ESNA ASHARI, RAMINE NIKOUKHAH, AND STEPHEN L. CAMPBELL Abstract. The problem of maximizing a quadratic function subject to an ellipsoidal

More information

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII

CONTROL SYSTEMS, ROBOTICS AND AUTOMATION CONTENTS VOLUME VII CONTENTS VOLUME VII Control of Linear Multivariable Systems 1 Katsuhisa Furuta,Tokyo Denki University, School of Science and Engineering, Ishizaka, Hatoyama, Saitama, Japan 1. Linear Multivariable Systems

More information

DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION

DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial

More information

Adaptive Data Assimilation and Multi-Model Fusion

Adaptive Data Assimilation and Multi-Model Fusion Adaptive Data Assimilation and Multi-Model Fusion Pierre F.J. Lermusiaux, Oleg G. Logoutov and Patrick J. Haley Jr. Mechanical Engineering and Ocean Science and Engineering, MIT We thank: Allan R. Robinson

More information

A Computational Framework for Quantifying and Optimizing the Performance of Observational Networks in 4D-Var Data Assimilation

A Computational Framework for Quantifying and Optimizing the Performance of Observational Networks in 4D-Var Data Assimilation A Computational Framework for Quantifying and Optimizing the Performance of Observational Networks in 4D-Var Data Assimilation Alexandru Cioaca Computational Science Laboratory (CSL) Department of Computer

More information

R. Balan. Splaiul Independentei 313, Bucharest, ROMANIA D. Aur

R. Balan. Splaiul Independentei 313, Bucharest, ROMANIA D. Aur An On-line Robust Stabilizer R. Balan University "Politehnica" of Bucharest, Department of Automatic Control and Computers, Splaiul Independentei 313, 77206 Bucharest, ROMANIA radu@karla.indinf.pub.ro

More information

Observer design for rotating shafts excited by unbalances

Observer design for rotating shafts excited by unbalances Observer design for rotating shafts excited by unbalances R. S. Schittenhelm, Z. Wang, S. Rinderknecht Institute for Mechatronic Systems in Mechanical Engineering, Technische Universität Darmstadt, Germany

More information

Weak Constraints 4D-Var

Weak Constraints 4D-Var Weak Constraints 4D-Var Yannick Trémolet ECMWF Training Course - Data Assimilation May 1, 2012 Yannick Trémolet Weak Constraints 4D-Var May 1, 2012 1 / 30 Outline 1 Introduction 2 The Maximum Likelihood

More information

A Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization

A Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization A Family of Preconditioned Iteratively Regularized Methods For Nonlinear Minimization Alexandra Smirnova Rosemary A Renaut March 27, 2008 Abstract The preconditioned iteratively regularized Gauss-Newton

More information

Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method

Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method Mathematical Problems in Engineering Volume 1, Article ID 693453, 1 pages doi:11155/1/693453 Research Article A Note on the Solutions of the Van der Pol and Duffing Equations Using a Linearisation Method

More information

Parallel Algorithms for Four-Dimensional Variational Data Assimilation

Parallel Algorithms for Four-Dimensional Variational Data Assimilation Parallel Algorithms for Four-Dimensional Variational Data Assimilation Mie Fisher ECMWF October 24, 2011 Mie Fisher (ECMWF) Parallel 4D-Var October 24, 2011 1 / 37 Brief Introduction to 4D-Var Four-Dimensional

More information

The ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations

The ECMWF Hybrid 4D-Var and Ensemble of Data Assimilations The Hybrid 4D-Var and Ensemble of Data Assimilations Lars Isaksen, Massimo Bonavita and Elias Holm Data Assimilation Section lars.isaksen@ecmwf.int Acknowledgements to: Mike Fisher and Marta Janiskova

More information

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi 4th EnKF Workshop, April 2010 Relationship

More information

Intro to Linear & Nonlinear Optimization

Intro to Linear & Nonlinear Optimization ECE 174 Intro to Linear & Nonlinear Optimization i Ken Kreutz-Delgado ECE Department, UCSD Contact Information Course Website Accessible from http://dsp.ucsd.edu/~kreutz Instructor Ken Kreutz-Delgado kreutz@ece.ucsd.eduucsd

More information

NUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group

NUMERICAL METHODS. lor CHEMICAL ENGINEERS. Using Excel', VBA, and MATLAB* VICTOR J. LAW. CRC Press. Taylor & Francis Group NUMERICAL METHODS lor CHEMICAL ENGINEERS Using Excel', VBA, and MATLAB* VICTOR J. LAW CRC Press Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Croup,

More information

Introduction to Applied Linear Algebra with MATLAB

Introduction to Applied Linear Algebra with MATLAB Sigam Series in Applied Mathematics Volume 7 Rizwan Butt Introduction to Applied Linear Algebra with MATLAB Heldermann Verlag Contents Number Systems and Errors 1 1.1 Introduction 1 1.2 Number Representation

More information

SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS

SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS BIT 0006-3835/00/4004-0726 $15.00 2000, Vol. 40, No. 4, pp. 726 734 c Swets & Zeitlinger SYMMETRIC PROJECTION METHODS FOR DIFFERENTIAL EQUATIONS ON MANIFOLDS E. HAIRER Section de mathématiques, Université

More information

On the convergence of the iterative solution of the likelihood equations

On the convergence of the iterative solution of the likelihood equations On the convergence of the iterative solution of the likelihood equations R. Moddemeijer University of Groningen, Department of Computing Science, P.O. Box 800, NL-9700 AV Groningen, The Netherlands, e-mail:

More information

Model-based Fault Diagnosis Techniques Design Schemes, Algorithms, and Tools

Model-based Fault Diagnosis Techniques Design Schemes, Algorithms, and Tools Steven X. Ding Model-based Fault Diagnosis Techniques Design Schemes, Algorithms, and Tools Springer Notation XIX Part I Introduction, basic concepts and preliminaries 1 Introduction 3 1.1 Basic concepts

More information

Data Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality

Data Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality Data Assimilation for Weather Forecasting: Reducing the Curse of Dimensionality Philippe Toint (with S. Gratton, S. Gürol, M. Rincon-Camacho, E. Simon and J. Tshimanga) University of Namur, Belgium Leverhulme

More information

Is My CFD Mesh Adequate? A Quantitative Answer

Is My CFD Mesh Adequate? A Quantitative Answer Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011

More information

Regional Hazardous Weather Advisory Centres (RHWACs)

Regional Hazardous Weather Advisory Centres (RHWACs) Regional Hazardous Weather Advisory Centres (RHWACs) The following outlines the criteria for the selection of RHWACs based on operational and functional requirements 1. Basic Principles The RHWAC must:

More information

A Sobolev trust-region method for numerical solution of the Ginz

A Sobolev trust-region method for numerical solution of the Ginz A Sobolev trust-region method for numerical solution of the Ginzburg-Landau equations Robert J. Renka Parimah Kazemi Department of Computer Science & Engineering University of North Texas June 6, 2012

More information

Research and Development of Advanced Radar Data Quality Control and Assimilation for Nowcasting and Forecasting Severe Storms

Research and Development of Advanced Radar Data Quality Control and Assimilation for Nowcasting and Forecasting Severe Storms DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Research and Development of Advanced Radar Data Quality Control and Assimilation for Nowcasting and Forecasting Severe

More information

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF

Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi ECODYC10, Dresden 28 January 2010 Relationship

More information

SIMULATION OF THE STABILITY LOSS OF THE VON MISES TRUSS IN AN UNSYMMETRICAL STRESS STATE

SIMULATION OF THE STABILITY LOSS OF THE VON MISES TRUSS IN AN UNSYMMETRICAL STRESS STATE Engineering MECHANICS, Vol. 14, 2007 1 SIMULATION OF THE STABILITY LOSS OF THE VON MISES TRUSS IN AN UNSYMMETRICAL STRESS STATE Petr Frantík* A hypothesis regarding the dynamical process of stability loss

More information

Active Flutter Control using an Adjoint Method

Active Flutter Control using an Adjoint Method 44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an

More information

Digital Control Engineering Analysis and Design

Digital Control Engineering Analysis and Design Digital Control Engineering Analysis and Design M. Sami Fadali Antonio Visioli AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO Academic Press is

More information

Dynamic 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. 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 information

Computational challenges in Numerical Weather Prediction

Computational challenges in Numerical Weather Prediction Computational challenges in Numerical Weather Prediction Mike Cullen Oxford 15 September 2008 Contents This presentation covers the following areas Historical background Current challenges Why does it

More information

Suboptimal Open-loop Control Using POD. Stefan Volkwein

Suboptimal Open-loop Control Using POD. Stefan Volkwein Institute for Mathematics and Scientific Computing University of Graz, Austria PhD program in Mathematics for Technology Catania, May 22, 2007 Motivation Optimal control of evolution problems: min J(y,

More information

nonrobust estimation The n measurement vectors taken together give the vector X R N. The unknown parameter vector is P R M.

nonrobust estimation The n measurement vectors taken together give the vector X R N. The unknown parameter vector is P R M. Introduction to nonlinear LS estimation R. I. Hartley and A. Zisserman: Multiple View Geometry in Computer Vision. Cambridge University Press, 2ed., 2004. After Chapter 5 and Appendix 6. We will use x

More information

The Inversion Problem: solving parameters inversion and assimilation problems

The Inversion Problem: solving parameters inversion and assimilation problems The Inversion Problem: solving parameters inversion and assimilation problems UE Numerical Methods Workshop Romain Brossier romain.brossier@univ-grenoble-alpes.fr ISTerre, Univ. Grenoble Alpes Master 08/09/2016

More information

Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit

Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Convergence of Square Root Ensemble Kalman Filters in the Large Ensemble Limit Evan Kwiatkowski, Jan Mandel University of Colorado Denver December 11, 2014 OUTLINE 2 Data Assimilation Bayesian Estimation

More information

Conditioning of the Weak-Constraint Variational Data Assimilation Problem for Numerical Weather Prediction. Adam El-Said

Conditioning of the Weak-Constraint Variational Data Assimilation Problem for Numerical Weather Prediction. Adam El-Said THE UNIVERSITY OF READING DEPARTMENT OF MATHEMATICS AND STATISTICS Conditioning of the Weak-Constraint Variational Data Assimilation Problem for Numerical Weather Prediction Adam El-Said Thesis submitted

More information

Iterative Feedback Tuning

Iterative Feedback Tuning Iterative Feedback Tuning Michel Gevers CESAME - UCL Louvain-la-Neuve Belgium Collaboration : H. Hjalmarsson, S. Gunnarsson, O. Lequin, E. Bosmans, L. Triest, M. Mossberg Outline Problem formulation Iterative

More information

Numerical Methods for Large-Scale Nonlinear Equations

Numerical Methods for Large-Scale Nonlinear Equations Slide 1 Numerical Methods for Large-Scale Nonlinear Equations Homer Walker MA 512 April 28, 2005 Inexact Newton and Newton Krylov Methods a. Newton-iterative and inexact Newton methods. Slide 2 i. Formulation

More information

Doppler radial wind spatially correlated observation error: operational implementation and initial results

Doppler radial wind spatially correlated observation error: operational implementation and initial results Doppler radial wind spatially correlated observation error: operational implementation and initial results D. Simonin, J. Waller, G. Kelly, S. Ballard,, S. Dance, N. Nichols (Met Office, University of

More information

Fixed Order H Controller for Quarter Car Active Suspension System

Fixed Order H Controller for Quarter Car Active Suspension System Fixed Order H Controller for Quarter Car Active Suspension System B. Erol, A. Delibaşı Abstract This paper presents an LMI based fixed-order controller design for quarter car active suspension system in

More information

7.1 Sampling Error The Need for Sampling Distributions

7.1 Sampling Error The Need for Sampling Distributions 7.1 Sampling Error The Need for Sampling Distributions Tom Lewis Fall Term 2009 Tom Lewis () 7.1 Sampling Error The Need for Sampling Distributions Fall Term 2009 1 / 5 Outline 1 Tom Lewis () 7.1 Sampling

More information

Non-linear least squares

Non-linear least squares Non-linear least squares Concept of non-linear least squares We have extensively studied linear least squares or linear regression. We see that there is a unique regression line that can be determined

More information

Stochastic Analogues to Deterministic Optimizers

Stochastic Analogues to Deterministic Optimizers Stochastic Analogues to Deterministic Optimizers ISMP 2018 Bordeaux, France Vivak Patel Presented by: Mihai Anitescu July 6, 2018 1 Apology I apologize for not being here to give this talk myself. I injured

More information

Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra

Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs

More information

Review on Aircraft Gain Scheduling

Review on Aircraft Gain Scheduling Review on Aircraft Gain Scheduling Z. Y. Kung * and I. F. Nusyirwan a Department of Aeronautical Engineering, Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310 Skudai, Johor, Malaysia.

More information

Riccati difference equations to non linear extended Kalman filter constraints

Riccati difference equations to non linear extended Kalman filter constraints International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R

More information

Ocean data assimilation for reanalysis

Ocean data assimilation for reanalysis Ocean data assimilation for reanalysis Matt Martin. ERA-CLIM2 Symposium, University of Bern, 14 th December 2017. Contents Introduction. On-going developments to improve ocean data assimilation for reanalysis.

More information

Eigenstructure Assignment for Helicopter Hover Control

Eigenstructure Assignment for Helicopter Hover Control Proceedings of the 17th World Congress The International Federation of Automatic Control Eigenstructure Assignment for Helicopter Hover Control Andrew Pomfret Stuart Griffin Tim Clarke Department of Electronics,

More information

Effects of a convective GWD parameterization in the global forecast system of the Met Office Unified Model in Korea

Effects of a convective GWD parameterization in the global forecast system of the Met Office Unified Model in Korea Effects of a convective GWD parameterization in the global forecast system of the Met Office Unified Model in Korea Young-Ha Kim 1, Hye-Yeong Chun 1, and Dong-Joon Kim 2 1 Yonsei University, Seoul, Korea

More information

Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model

Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Newton-Krylov-Schwarz Method for a Spherical Shallow Water Model Chao Yang 1 and Xiao-Chuan Cai 2 1 Institute of Software, Chinese Academy of Sciences, Beijing 100190, P. R. China, yang@mail.rdcps.ac.cn

More information

Basic Aspects of Discretization

Basic Aspects of Discretization Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)

More information

Linear System Theory. Wonhee Kim Lecture 1. March 7, 2018

Linear System Theory. Wonhee Kim Lecture 1. March 7, 2018 Linear System Theory Wonhee Kim Lecture 1 March 7, 2018 1 / 22 Overview Course Information Prerequisites Course Outline What is Control Engineering? Examples of Control Systems Structure of Control Systems

More information

Uncertainty quantification for inverse problems with a weak wave-equation constraint

Uncertainty quantification for inverse problems with a weak wave-equation constraint Uncertainty quantification for inverse problems with a weak wave-equation constraint Zhilong Fang*, Curt Da Silva*, Rachel Kuske** and Felix J. Herrmann* *Seismic Laboratory for Imaging and Modeling (SLIM),

More information

An Optimal Control Approach to Sensor / Actuator Placement for Optimal Control of High Performance Buildings

An Optimal Control Approach to Sensor / Actuator Placement for Optimal Control of High Performance Buildings Purdue University Purdue e-pubs International High Performance Buildings Conference School of Mechanical Engineering 01 An Optimal Control Approach to Sensor / Actuator Placement for Optimal Control of

More information

Revision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation

Revision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation Revision of TR-9-25: A Hybrid Variational/Ensemble ilter Approach to Data Assimilation Adrian Sandu 1 and Haiyan Cheng 1 Computational Science Laboratory Department of Computer Science Virginia Polytechnic

More information

Scientific Data Computing: Lecture 3

Scientific Data Computing: Lecture 3 Scientific Data Computing: Lecture 3 Benson Muite benson.muite@ut.ee 23 April 2018 Outline Monday 10-12, Liivi 2-207 Monday 12-14, Liivi 2-205 Topics Introduction, statistical methods and their applications

More information

Convergence of the Ensemble Kalman Filter in Hilbert Space

Convergence of the Ensemble Kalman Filter in Hilbert Space Convergence of the Ensemble Kalman Filter in Hilbert Space Jan Mandel Center for Computational Mathematics Department of Mathematical and Statistical Sciences University of Colorado Denver Parts based

More information

Correlations of control variables in variational data assimilation

Correlations of control variables in variational data assimilation Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 37: 62 63, April 2 A Correlations of control variables in variational data assimilation D. Katz a,a.s.lawless* a,n.k.nichols

More information

6.252 NONLINEAR PROGRAMMING LECTURE 10 ALTERNATIVES TO GRADIENT PROJECTION LECTURE OUTLINE. Three Alternatives/Remedies for Gradient Projection

6.252 NONLINEAR PROGRAMMING LECTURE 10 ALTERNATIVES TO GRADIENT PROJECTION LECTURE OUTLINE. Three Alternatives/Remedies for Gradient Projection 6.252 NONLINEAR PROGRAMMING LECTURE 10 ALTERNATIVES TO GRADIENT PROJECTION LECTURE OUTLINE Three Alternatives/Remedies for Gradient Projection Two-Metric Projection Methods Manifold Suboptimization Methods

More information

Intro to Linear & Nonlinear Optimization

Intro to Linear & Nonlinear Optimization ECE 174 Intro to Linear & Nonlinear Optimization Ken Kreutz-Delgado ECE Department UCSD Version 10.5.2017 Contact Information Fall 2017 Course Website Accessible from http://dsp.ucsd.edu/~kreutz/; Piazza:

More information

SCOTT DAWSON. stdawson/ 1200 E. California Blvd, MC Pasadena, CA (609)

SCOTT DAWSON.   stdawson/ 1200 E. California Blvd, MC Pasadena, CA (609) RESEARCH INTERESTS SCOTT DAWSON http://www.caltech.edu/ stdawson/ 1200 E. California Blvd, MC 105-50 Pasadena, CA 91125 +1 (609) 356-8267 stdawson@caltech.edu Modeling, optimization and control in fluid

More information

Impact of GPS and TMI Precipitable Water Data on Mesoscale Numerical Weather Prediction Model Forecasts

Impact of GPS and TMI Precipitable Water Data on Mesoscale Numerical Weather Prediction Model Forecasts Journal of the Meteorological Society of Japan, Vol. 82, No. 1B, pp. 453--457, 2004 453 Impact of GPS and TMI Precipitable Water Data on Mesoscale Numerical Weather Prediction Model Forecasts Ko KOIZUMI

More information

Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake

Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme

More information

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers

A Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned

More information

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY

A MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,

More information

Satellites, Weather and Climate Module??: Polar Vortex

Satellites, Weather and Climate Module??: Polar Vortex Satellites, Weather and Climate Module??: Polar Vortex SWAC Jan 2014 AKA Circumpolar Vortex Science or Hype? Will there be one this year? Today s objectives Pre and Post exams What is the Polar Vortex

More information

Some ideas for Ensemble Kalman Filter

Some ideas for Ensemble Kalman Filter Some ideas for Ensemble Kalman Filter Former students and Eugenia Kalnay UMCP Acknowledgements: UMD Chaos-Weather Group: Brian Hunt, Istvan Szunyogh, Ed Ott and Jim Yorke, Kayo Ide, and students Former

More information

Proc. 9th IFAC/IFORS/IMACS/IFIP/ Symposium on Large Scale Systems: Theory and Applications (LSS 2001), 2001, pp

Proc. 9th IFAC/IFORS/IMACS/IFIP/ Symposium on Large Scale Systems: Theory and Applications (LSS 2001), 2001, pp INVARIANT POLYHEDRA AND CONTROL OF LARGE SCALE SYSTEMS Jean-Claude Hennet LAAS - CNRS 7, Avenue du Colonel Roche, 31077 Toulouse Cedex 4, FRANCE PHONE: (+33) 5 61 33 63 13, FAX: (+33) 5 61 33 69 36, e-mail:

More information