GENERALIZED STABILITY OF THE TWO-LAYER MODEL
|
|
- Aubrey Carson
- 5 years ago
- Views:
Transcription
1 GENERALIZED STABILITY OF THE TWO-LAYER MODEL The simplest mid-latitude jet model supporting the baroclinic growth mechanism is the two-layer model The equations for the barotropic and baroclinic geostrophic streamfunction perturbations in this model assuming meridional wavenumber zero are: t + U ψ m m x x + β ψ m x + U ψ T T 0 x x t + U ψ T m x x λ ψ T + β ψ T x + U ψ m T + λ ψ x x m 0 Assuming the general temporal dependence: ψ m x, t ˆψ m te ikx ψ T x, t ˆψ T te ikx the resulting dynamical system can be written in standard form: d ˆψm dt ˆψ T βi U k Mik U T ik λ k λ +k U T ik βik U λ +k M ik ˆψm ˆψ T The essential baroclinic dynamics is retained with the simplifying assumption U m β 0 The additional assumptions U T and λ 0 correspond to choosing space scale /λ 0 and time scale /λ 0 U T /λ 0 is taken to have the typical value of 000 km and a typical value for U T is 5 m/s so the time scale is 077 day In order to see the effect of static stability on baroclinic growth the nondimensional Rossby radius /λ is retained as a variable The system under these assumptions is: d ˆψm dt ˆψ T 0 ik ik λ k λ +k 0 ˆψm ˆψ T
2 Writing the above system as dψ Aψ first we solve the eigenproblems for A and A+A A is dt the Hermitian transpose of A: A 0 ik ik λ k λ +k 0 A + A 0 ikλ λ +k ikλ λ +k 0 Recall that the maximum real part of the eigenvalue of A is the growth rate of the perturbation of modal form For a general perturbation to a system dx Ax the instantaneous perturbation dt growth rate can be found by differentiating x x x: d ln x dt xa + A x xx The rhs of the above expression is maximized by the eigenfunction of A + A associated with its necessarily real maximum eigenvalue This eigenvalue is called the numerical range The eigenvalues of A are: λ k σ A ±k λ + k, and the maximum real part of the eigenvalues is shown as a function of k and λ in Fig The eigenvalues of A+A are: σ A+A ±k λ k λ + k
3 5 5 λ k Figure : The modal growth rate as a function of wavenumber k and inverse Rossby radius λ The maximum of these necessarily real eigenvalues is the numerical range and it is shown as a function of k and λ in Fig The eigenvectors of A are: ˆψm ˆψ T ±i λ k λ+k The upper and lower eigenvector streamfunctions are: ˆψ ψm ˆ + ˆ ψ T ˆψ ψˆ M ψ ˆ T Normalizing so that ˆψ : k i +k ˆψ ˆψ ±i k +k
4 5 5 λ k Figure : The maximum instantaneous growth rate the numerical range as a function of wavenumber k and inverse Rossby radius λ From which we see that the phase shift between the upper and lower layer eigenfunction streamfunctions vanishes for k > which is the stable region In the unstable region k < the modes have phase shift θ arctan k +k k +k with negative phase for the growing mode This phase shift between the upper and lower mode streamfunction is shown as a function of k in Fig Turning to the numerical range the eigenfunctions of A + A are: ˆψ ˆψ i ±i The maximum instantaneous growth occurs for upper layer phase π/ relative to the lower
5 θ k Figure : Phase shift between the upper and lower layers for the growing mode as a function of k for λ 5
6 numerical range λ mode growth rate k Figure : Growth rate of the mode compared to the maximum instantaneous growth rate for λ layer; and maximum decay for phase difference π/ The numerical range is greater or equal to the maximum mode growth rate as can be seen in Fig in which these quantities are shown as a function of k for γ The standard solution form Ψτ R τ Ψ0 maps the initial state at t0 to the final state at t τ in terms of the propagator R τ e Aτ written as a matrix exponential expm Aτ in Matlab we can evaluate the e Aτ by noticing that e Aτ Ee Λτ E in which E is the matrix with columns the eigenvectors of A and Λ is the diagonal matrix of the eigenvalues: e Aτ i γ / i γ / e kγ / τ 0 0 e kγ/ τ +iγ/ / iγ / / 6
7 This results in the propagator: e Aτ coshkγ / τ i γ / sinhkγ/ τ iγ / sinhkγ / τ coshkγ / τ Taking our example values k, λ so that γ / the propagator becomes: e Aτ cos τ i sin τ i sin τ cos τ We can find the optimal growth as a function of τ using only eigenanalysis by noticing that R τ is the square root of the first eigenvalue of e A τ e Aτ and the associated eigenvector is the optimal perturbation use a Rayleigh quotient argument to verify this e A τ e Aτ cos τ + sin τ i cos τ sin τ i cos τ sin τ cos τ + sin τ The leading eigenvalue of this matrix is: R τ + sin τ + + sin τ The optimal growth as a function of τ for k and γ is shown in Fig 5 The first maximum occurs at τ max π with growth R τ This maximum is called the global optimal because there is no damping in this problem the upper and lower layers come 7
8 R τ τ Figure 5: Growth of the optimal perturbation for k and λ in and out of phase and this interference results in multiple maxima; there is typically only one global maximum when damping is allowed for The eigenvector at τ max π is: ˆψ ˆψ This corresponds to phase shift of π between the layers with the upper level lagging by twice as much phase as in the case of the numerical range This norm of R τ can be found in general using normr τ in Matlab The norm is the first singular value of the propagator and this growth is attained by the initial perturbation in the form of the right singular vector of the propagator these are found in general using the singular value decomposition command in Matlab: [U,S,V] svdr τ 8
The Singular Value Decomposition
The Singular Value Decomposition Philippe B. Laval KSU Fall 2015 Philippe B. Laval (KSU) SVD Fall 2015 1 / 13 Review of Key Concepts We review some key definitions and results about matrices that will
More informationIn two-dimensional barotropic flow, there is an exact relationship between mass
19. Baroclinic Instability In two-dimensional barotropic flow, there is an exact relationship between mass streamfunction ψ and the conserved quantity, vorticity (η) given by η = 2 ψ.the evolution of the
More informationThe Eady problem of baroclinic instability described in section 19a was shown to
0. The Charney-Stern Theorem The Eady problem of baroclinic instability described in section 19a was shown to be remarkably similar to the Rayleigh instability of barotropic flow described in Chapter 18.
More informationJet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities. Mark Fruman
p. 1/24 Jet Formation in the Equatorial Oceans Through Barotropic and Inertial Instabilities Mark Fruman Bach Lien Hua, Richard Schopp, Marc d Orgeville, Claire Ménesguen LPO IFREMER, Brest, France IAU
More informationCan a Simple Two-Layer Model Capture the Structure of Easterly Waves?
Can a Simple Two-Layer Model Capture the Structure of Easterly Waves? Cheryl L. Lacotta 1 Introduction Most tropical storms in the Atlantic, and even many in the eastern Pacific, are due to disturbances
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS BRIAN F. FARRELL Division of Engineering and Applied Sciences, Harvard University Pierce Hall, 29
More informationOPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS
Stochastics and Dynamics, Vol. 2, No. 3 (22 395 42 c World Scientific Publishing Company OPTIMAL PERTURBATION OF UNCERTAIN SYSTEMS Stoch. Dyn. 22.2:395-42. Downloaded from www.worldscientific.com by HARVARD
More informationChap 3. Linear Algebra
Chap 3. Linear Algebra Outlines 1. Introduction 2. Basis, Representation, and Orthonormalization 3. Linear Algebraic Equations 4. Similarity Transformation 5. Diagonal Form and Jordan Form 6. Functions
More informationAnnounce Statistical Motivation Properties Spectral Theorem Other ODE Theory Spectral Embedding. Eigenproblems I
Eigenproblems I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Eigenproblems I 1 / 33 Announcements Homework 1: Due tonight
More informationMATH 5524 MATRIX THEORY Problem Set 5
MATH 554 MATRIX THEORY Problem Set 5 Posted Tuesday April 07. Due Tuesday 8 April 07. [Late work is due on Wednesday 9 April 07.] Complete any four problems, 5 points each. Recall the definitions of the
More information( ) = 9φ 1, ( ) = 4φ 2.
Chemistry 46 Dr Jean M Standard Homework Problem Set 6 Solutions The Hermitian operator A ˆ is associated with the physical observable A Two of the eigenfunctions of A ˆ are and These eigenfunctions are
More informationComputational Methods CMSC/AMSC/MAPL 460. Eigenvalues and Eigenvectors. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 Eigenvalues and Eigenvectors Ramani Duraiswami, Dept. of Computer Science Eigen Values of a Matrix Recap: A N N matrix A has an eigenvector x (non-zero) with corresponding
More informationAccurate Low-Dimensional Approximation of the Linear Dynamics of Fluid Flow
15 SEPTEMBER 21 FARRELL AND IOANNOU 2771 Accurate Low-Dimensional Approximation of the Linear Dynamics of Fluid Flow BRIAN F. FARRELL Department for Earth and Planetary Sciences, Harvard University, Cambridge,
More information2. Baroclinic Instability and Midlatitude Dynamics
2. Baroclinic Instability and Midlatitude Dynamics Midlatitude Jet Stream Climatology (Atlantic and Pacific) Copyright 26 Emily Shuckburgh, University of Cambridge. Not to be quoted or reproduced without
More informationAnnounce Statistical Motivation ODE Theory Spectral Embedding Properties Spectral Theorem Other. Eigenproblems I
Eigenproblems I CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Doug James (and Justin Solomon) CS 205A: Mathematical Methods Eigenproblems I 1 / 33 Announcements Homework 1: Due tonight
More informationMATH 23 Exam 2 Review Solutions
MATH 23 Exam 2 Review Solutions Problem 1. Use the method of reduction of order to find a second solution of the given differential equation x 2 y (x 0.1875)y = 0, x > 0, y 1 (x) = x 1/4 e 2 x Solution
More informationChapter 13 Instability on non-parallel flow Introduction and formulation
Chapter 13 Instability on non-parallel flow. 13.1 Introduction and formulation We have concentrated our discussion on the instabilities of parallel, zonal flows. There is the largest amount of literature
More informationSpacecraft and Aircraft Dynamics
Spacecraft and Aircraft Dynamics Matthew M. Peet Illinois Institute of Technology Lecture 11: Longitudinal Dynamics Aircraft Dynamics Lecture 11 In this Lecture we will cover: Longitudinal Dynamics: Finding
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationLecture 6. Four postulates of quantum mechanics. The eigenvalue equation. Momentum and energy operators. Dirac delta function. Expectation values
Lecture 6 Four postulates of quantum mechanics The eigenvalue equation Momentum and energy operators Dirac delta function Expectation values Objectives Learn about eigenvalue equations and operators. Learn
More informationj=1 u 1jv 1j. 1/ 2 Lemma 1. An orthogonal set of vectors must be linearly independent.
Lecture Notes: Orthogonal and Symmetric Matrices Yufei Tao Department of Computer Science and Engineering Chinese University of Hong Kong taoyf@cse.cuhk.edu.hk Orthogonal Matrix Definition. Let u = [u
More informationGFD II: Balance Dynamics ATM S 542
GFD II: Balance Dynamics ATM S 542 DARGAN M. W. FRIERSON UNIVERSITY OF WASHINGTON, DEPARTMENT OF ATMOSPHERIC SCIENCES WEEK 8 SLIDES Eady Model Growth rates (imaginary part of frequency) Stable for large
More informationsemi-infinite Eady model
The semi-infinite Eady model Using a non-modal decomposition technique based on the potential vorticity (PV) perspective, the optimal perturbation or singular vector (SV) of the Eady model without upper
More informationESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves
ESCI 343 Atmospheric Dynamics II Lesson 11 - Rossby Waves Reference: An Introduction to Dynamic Meteorology (4 rd edition), J.R. Holton Atmosphere-Ocean Dynamics, A.E. Gill Fundamentals of Atmospheric
More information1. Comparison of stability analysis to previous work
. Comparison of stability analysis to previous work The stability problem (6.4) can be understood in the context of previous work. Benjamin (957) and Yih (963) have studied the stability of fluid flowing
More information9 Rossby Waves. 9.1 Non-divergent barotropic vorticity equation. CSU ATS601 Fall (Holton Chapter 7, Vallis Chapter 5)
9 Rossby Waves (Holton Chapter 7, Vallis Chapter 5) 9.1 Non-divergent barotropic vorticity equation We are now at a point that we can discuss our first fundamental application of the equations of motion:
More informationAFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES
AFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES Victoria Dollar RTG Seminar Research - Spring 2018 April 16, 2018 Victoria Dollar ASU April 16, 2018 1 / 26 Overview Introduction Rossby waves and
More informationChapter 9. Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation
Chapter 9 Geostrophy, Quasi-Geostrophy and the Potential Vorticity Equation 9.1 Geostrophy and scaling. We examined in the last chapter some consequences of the dynamical balances for low frequency, nearly
More informationFrictional Damping of Baroclinic Waves
Frictional Damping of Baroclinic Waves HHH, 26th May 2006 Bob Plant With thanks to Stephen Belcher, Brian Hoskins, Dan Adamson Motivation Control simulation, T+60 Simulation with no boundary layer turbulence,
More informationApproximating Optimal State Estimation
Approximating Optimal State Estimation Brian F. Farrell and Petros J. Ioannou Harvard University Cambridge MA 238, U.S.A farrell@deas.harvard.edu, pji@cc.uoa.gr ABSTRACT Minimizing forecast error requires
More information18.06 Problem Set 8 - Solutions Due Wednesday, 14 November 2007 at 4 pm in
806 Problem Set 8 - Solutions Due Wednesday, 4 November 2007 at 4 pm in 2-06 08 03 Problem : 205+5+5+5 Consider the matrix A 02 07 a Check that A is a positive Markov matrix, and find its steady state
More information1 Distributions (due January 22, 2009)
Distributions (due January 22, 29). The distribution derivative of the locally integrable function ln( x ) is the principal value distribution /x. We know that, φ = lim φ(x) dx. x ɛ x Show that x, φ =
More informationCHMY 564 Homework #3 2Feb17 Due Friday, 10Feb17
564-17 Lec 11 Mon 6feb17 CHMY 564 Homework #3 2Feb17 Due Friday, 10Feb17 1. (a) Find the eigenfunctions and eigenvalues of the x component of linear momentum. (b) Why can a momentum eigenfunction never
More informationLecture 8: Linear Algebra Background
CSE 521: Design and Analysis of Algorithms I Winter 2017 Lecture 8: Linear Algebra Background Lecturer: Shayan Oveis Gharan 2/1/2017 Scribe: Swati Padmanabhan Disclaimer: These notes have not been subjected
More information13-2 Text: 28-30; AB: 1.3.3, 3.2.3, 3.4.2, 3.5, 3.6.2; GvL Eigen2
The QR algorithm The most common method for solving small (dense) eigenvalue problems. The basic algorithm: QR without shifts 1. Until Convergence Do: 2. Compute the QR factorization A = QR 3. Set A :=
More information1. The Polar Decomposition
A PERSONAL INTERVIEW WITH THE SINGULAR VALUE DECOMPOSITION MATAN GAVISH Part. Theory. The Polar Decomposition In what follows, F denotes either R or C. The vector space F n is an inner product space with
More informationBasic Calculus Review
Basic Calculus Review Lorenzo Rosasco ISML Mod. 2 - Machine Learning Vector Spaces Functionals and Operators (Matrices) Vector Space A vector space is a set V with binary operations +: V V V and : R V
More informationMath Theory of Partial Differential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions.
Math 412-501 Theory of Partial ifferential Equations Lecture 3-2: Spectral properties of the Laplacian. Bessel functions. Eigenvalue problem: 2 φ + λφ = 0 in, ( αφ + β φ ) n = 0, where α, β are piecewise
More informationModeling Large-Scale Atmospheric and Oceanic Flows 2
Modeling Large-Scale Atmospheric and Oceanic Flows V. Zeitlin known s of the Laboratoire de Météorologie Dynamique, Univ. P. and M. Curie, Paris Mathematics of the Oceans, Fields Institute, Toronto, 3
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More information1 Directional Derivatives and Differentiability
Wednesday, January 18, 2012 1 Directional Derivatives and Differentiability Let E R N, let f : E R and let x 0 E. Given a direction v R N, let L be the line through x 0 in the direction v, that is, L :=
More informationQuantum Physics Lecture 8
Quantum Physics ecture 8 Steady state Schroedinger Equation (SSSE): eigenvalue & eigenfunction particle in a box re-visited Wavefunctions and energy states normalisation probability density Expectation
More information1 Linearity and Linear Systems
Mathematical Tools for Neuroscience (NEU 34) Princeton University, Spring 26 Jonathan Pillow Lecture 7-8 notes: Linear systems & SVD Linearity and Linear Systems Linear system is a kind of mapping f( x)
More informationdf(x) = h(x) dx Chemistry 4531 Mathematical Preliminaries Spring 2009 I. A Primer on Differential Equations Order of differential equation
Chemistry 4531 Mathematical Preliminaries Spring 009 I. A Primer on Differential Equations Order of differential equation Linearity of differential equation Partial vs. Ordinary Differential Equations
More informationLinear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions
Linear Algebra, part 2 Eigenvalues, eigenvectors and least squares solutions Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Main problem of linear algebra 2: Given
More informationLinear Least Squares. Using SVD Decomposition.
Linear Least Squares. Using SVD Decomposition. Dmitriy Leykekhman Spring 2011 Goals SVD-decomposition. Solving LLS with SVD-decomposition. D. Leykekhman Linear Least Squares 1 SVD Decomposition. For any
More informationNonlinear baroclinic dynamics of surface cyclones crossing a zonal jet
Nonlinear baroclinic dynamics of surface cyclones crossing a zonal jet Jean-Baptiste GILET, Matthieu Plu and Gwendal Rivière CNRM/GAME (Météo-France, CNRS) 3rd THORPEX International Science Symposium Monterey,
More informationEddy PV Fluxes in a One Dimensional Model of Quasi-Geostrophic Turbulence
Eddy PV Fluxes in a One Dimensional Model of Quasi-Geostrophic Turbulence Christos M.Mitas Introduction. Motivation Understanding eddy transport of heat and momentum is crucial to developing closure schemes
More informationExercise 8: Frequency Response of MIMO Systems
Exercise 8: Frequency Response of MIMO Systems 8 Singular Value Decomposition (SVD The Singular Value Decomposition plays a central role in MIMO frequency response analysis Let s recall some concepts from
More informationApplied Mathematics 205. Unit V: Eigenvalue Problems. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit V: Eigenvalue Problems Lecturer: Dr. David Knezevic Unit V: Eigenvalue Problems Chapter V.2: Fundamentals 2 / 31 Eigenvalues and Eigenvectors Eigenvalues and eigenvectors of
More information6 Two-layer shallow water theory.
6 Two-layer shallow water theory. Wewillnowgoontolookatashallowwatersystemthathastwolayersofdifferent density. This is the next level of complexity and a simple starting point for understanding the behaviour
More informationNon-linear Stability of Steady Geophysical Flows
Non-linear Stability of Steady Geophysical Flows Di Qi, and Andrew J. Majda Courant Institute of Mathematical Sciences Fall 2016 Advanced Topics in Applied Math Di Qi, and Andrew J. Majda (CIMS) Non-linear
More informationTutorials in Optimization. Richard Socher
Tutorials in Optimization Richard Socher July 20, 2008 CONTENTS 1 Contents 1 Linear Algebra: Bilinear Form - A Simple Optimization Problem 2 1.1 Definitions........................................ 2 1.2
More informationRadiating Instability of a Meridional Boundary Current
2294 J O U R N A L O F P H Y S I C A L O C E A N O G R A P H Y VOLUME 38 Radiating Instability of a Meridional Boundary Current HRISTINA G. HRISTOVA MIT WHOI Joint Program in Oceanography, Woods Hole,
More informationFiber Optics. Equivalently θ < θ max = cos 1 (n 0 /n 1 ). This is geometrical optics. Needs λ a. Two kinds of fibers:
Waves can be guided not only by conductors, but by dielectrics. Fiber optics cable of silica has nr varying with radius. Simplest: core radius a with n = n 1, surrounded radius b with n = n 0 < n 1. Total
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters, transforms,
More informationChapter 9. Barotropic Instability. 9.1 Linearized governing equations
Chapter 9 Barotropic Instability The ossby wave is the building block of low ossby number geophysical fluid dynamics. In this chapter we learn how ossby waves can interact with each other to produce a
More informationLecture 3: Review of Linear Algebra
ECE 83 Fall 2 Statistical Signal Processing instructor: R Nowak, scribe: R Nowak Lecture 3: Review of Linear Algebra Very often in this course we will represent signals as vectors and operators (eg, filters,
More informationLecture 4: Analysis of MIMO Systems
Lecture 4: Analysis of MIMO Systems Norms The concept of norm will be extremely useful for evaluating signals and systems quantitatively during this course In the following, we will present vector norms
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
More informationMATH 3321 Sample Questions for Exam 3. 3y y, C = Perform the indicated operations, if possible: (a) AC (b) AB (c) B + AC (d) CBA
MATH 33 Sample Questions for Exam 3. Find x and y so that x 4 3 5x 3y + y = 5 5. x = 3/7, y = 49/7. Let A = 3 4, B = 3 5, C = 3 Perform the indicated operations, if possible: a AC b AB c B + AC d CBA AB
More informationDependence of Singular Vector Structure and Evolution on the Choice of Norm
1NOVEMBER 00 KIM AND MORGAN 3099 Dependence of Singular Vector Structure and Evolution on the Choice of Norm HYUN MEE KIM AND MICHAEL C. MORGAN Department of Atmospheric and Oceanic Sciences, University
More information22.2. Applications of Eigenvalues and Eigenvectors. Introduction. Prerequisites. Learning Outcomes
Applications of Eigenvalues and Eigenvectors 22.2 Introduction Many applications of matrices in both engineering and science utilize eigenvalues and, sometimes, eigenvectors. Control theory, vibration
More informationSalmon: Lectures on partial differential equations
6. The wave equation Of the 3 basic equations derived in the previous section, we have already discussed the heat equation, (1) θ t = κθ xx + Q( x,t). In this section we discuss the wave equation, () θ
More information22.3. Repeated Eigenvalues and Symmetric Matrices. Introduction. Prerequisites. Learning Outcomes
Repeated Eigenvalues and Symmetric Matrices. Introduction In this Section we further develop the theory of eigenvalues and eigenvectors in two distinct directions. Firstly we look at matrices where one
More informationPhysics 505 Homework No. 1 Solutions S1-1
Physics 505 Homework No s S- Some Preliminaries Assume A and B are Hermitian operators (a) Show that (AB) B A dx φ ABψ dx (A φ) Bψ dx (B (A φ)) ψ dx (B A φ) ψ End (b) Show that AB [A, B]/2+{A, B}/2 where
More informationSG2221 Wave Motion and Hydrodinamic Stability. MATLAB Project on 2D Poiseuille Flow Alessandro Ceci
SG2221 Wave Motion and Hydrodinamic Stability MATLAB Project on 2D Poiseuille Flow Alessandro Ceci Base Flow 2D steady Incompressible Flow Flow driven by a constant pressure gradient Fully developed flow
More informationReview Sheet 2 Solutions
Review Sheet Solutions 1. If y x 3 x and dx dt 5, find dy dt when x. We have that dy dt 3 x dx dt dx dt 3 x 5 5, and this is equal to 3 5 10 70 when x.. A spherical balloon is being inflated so that its
More informationQuantum Computing Lecture 2. Review of Linear Algebra
Quantum Computing Lecture 2 Review of Linear Algebra Maris Ozols Linear algebra States of a quantum system form a vector space and their transformations are described by linear operators Vector spaces
More informationUnderstanding inertial instability on the f-plane with complete Coriolis force
Understanding inertial instability on the f-plane with complete Coriolis force Abstract Vladimir Zeitlin Laboratory of Dynamical Meteorology, University P. and M. Curie and Ecole Normale Supérieure, Paris,
More informationThe Singular Value Decomposition (SVD) and Principal Component Analysis (PCA)
Chapter 5 The Singular Value Decomposition (SVD) and Principal Component Analysis (PCA) 5.1 Basics of SVD 5.1.1 Review of Key Concepts We review some key definitions and results about matrices that will
More informationr 2. (5) This agrees of course with the traditional multipole way of writing it: for an axial point dipole at the coordinate origin
1. A Digression on Dipoles From the perspective of potential theory, the dipole represents the point source obtained by allowing two point charges of opposite sign to approach each other, in such a way
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer. Lecture 5, January 27, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Fall Semester 2006 Christopher J. Cramer Lecture 5, January 27, 2006 Solved Homework (Homework for grading is also due today) We are told
More informationEE5120 Linear Algebra: Tutorial 7, July-Dec Covers sec 5.3 (only powers of a matrix part), 5.5,5.6 of GS
EE5 Linear Algebra: Tutorial 7, July-Dec 7-8 Covers sec 5. (only powers of a matrix part), 5.5,5. of GS. Prove that the eigenvectors corresponding to different eigenvalues are orthonormal for unitary matrices.
More informationNumerical Methods I Eigenvalue Problems
Numerical Methods I Eigenvalue Problems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 2nd, 2014 A. Donev (Courant Institute) Lecture
More informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
More information11 a 12 a 21 a 11 a 22 a 12 a 21. (C.11) A = The determinant of a product of two matrices is given by AB = A B 1 1 = (C.13) and similarly.
C PROPERTIES OF MATRICES 697 to whether the permutation i 1 i 2 i N is even or odd, respectively Note that I =1 Thus, for a 2 2 matrix, the determinant takes the form A = a 11 a 12 = a a 21 a 11 a 22 a
More informationLecture 2: Computing functions of dense matrices
Lecture 2: Computing functions of dense matrices Paola Boito and Federico Poloni Università di Pisa Pisa - Hokkaido - Roma2 Summer School Pisa, August 27 - September 8, 2018 Introduction In this lecture
More information1 Numerical Approximation to Integrals and Expectation
1 Numerical Approximation to Integrals and Expectation Values Physics 21 Columbia University Copyright A. J. Millis All rights reserved 1.1 Overview These notes cover one approach to the numerical evaluation
More informationPDEs, part 3: Hyperbolic PDEs
PDEs, part 3: Hyperbolic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2011 Hyperbolic equations (Sections 6.4 and 6.5 of Strang). Consider the model problem (the
More informationPhysics 137A Quantum Mechanics Fall 2012 Midterm II - Solutions
Physics 37A Quantum Mechanics Fall 0 Midterm II - Solutions These are the solutions to the exam given to Lecture Problem [5 points] Consider a particle with mass m charge q in a simple harmonic oscillator
More informationComputational Methods CMSC/AMSC/MAPL 460. EigenValue decomposition Singular Value Decomposition. Ramani Duraiswami, Dept. of Computer Science
Computational Methods CMSC/AMSC/MAPL 460 EigenValue decomposition Singular Value Decomposition Ramani Duraiswami, Dept. of Computer Science Hermitian Matrices A square matrix for which A = A H is said
More informationHW #3 Solutions: M552 Spring (c)
HW #3 Solutions: M55 Spring 6. (4.-Trefethen & Bau: parts (a), (c) and (e)) Determine the SVDs of the following matrices (by hand calculation): (a) [ 3 (c) (e) [ ANS: In each case we seek A UΣV. The general
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationAssignment #10: Diagonalization of Symmetric Matrices, Quadratic Forms, Optimization, Singular Value Decomposition. Name:
Assignment #10: Diagonalization of Symmetric Matrices, Quadratic Forms, Optimization, Singular Value Decomposition Due date: Friday, May 4, 2018 (1:35pm) Name: Section Number Assignment #10: Diagonalization
More information18.06 Problem Set 10 - Solutions Due Thursday, 29 November 2007 at 4 pm in
86 Problem Set - Solutions Due Thursday, 29 November 27 at 4 pm in 2-6 Problem : (5=5+5+5) Take any matrix A of the form A = B H CB, where B has full column rank and C is Hermitian and positive-definite
More informationSeparation of Variables in Linear PDE: One-Dimensional Problems
Separation of Variables in Linear PDE: One-Dimensional Problems Now we apply the theory of Hilbert spaces to linear differential equations with partial derivatives (PDE). We start with a particular example,
More informationMath 215 HW #11 Solutions
Math 215 HW #11 Solutions 1 Problem 556 Find the lengths and the inner product of 2 x and y [ 2 + ] Answer: First, x 2 x H x [2 + ] 2 (4 + 16) + 16 36, so x 6 Likewise, so y 6 Finally, x, y x H y [2 +
More informationweak mean flow R. M. Samelson
An effective-β vector for linear planetary waves on a weak mean flow R. M. Samelson College of Oceanic and Atmospheric Sciences 14 COAS Admin Bldg Oregon State University Corvallis, OR 97331-553 USA rsamelson@coas.oregonstate.edu
More informationLinear Algebra using Dirac Notation: Pt. 2
Linear Algebra using Dirac Notation: Pt. 2 PHYS 476Q - Southern Illinois University February 6, 2018 PHYS 476Q - Southern Illinois University Linear Algebra using Dirac Notation: Pt. 2 February 6, 2018
More informationMATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators.
MATH 423 Linear Algebra II Lecture 33: Diagonalization of normal operators. Adjoint operator and adjoint matrix Given a linear operator L on an inner product space V, the adjoint of L is a transformation
More informationLecture 15 Review of Matrix Theory III. Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore
Lecture 15 Review of Matrix Theory III Dr. Radhakant Padhi Asst. Professor Dept. of Aerospace Engineering Indian Institute of Science - Bangalore Matrix An m n matrix is a rectangular or square array of
More information1/27/2010. With this method, all filed variables are separated into. from the basic state: Assumptions 1: : the basic state variables must
Lecture 5: Waves in Atmosphere Perturbation Method With this method, all filed variables are separated into two parts: (a) a basic state part and (b) a deviation from the basic state: Perturbation Method
More informationThe goal of this chapter is to study linear systems of ordinary differential equations: dt,..., dx ) T
1 1 Linear Systems The goal of this chapter is to study linear systems of ordinary differential equations: ẋ = Ax, x(0) = x 0, (1) where x R n, A is an n n matrix and ẋ = dx ( dt = dx1 dt,..., dx ) T n.
More information10 Shallow Water Models
10 Shallow Water Models So far, we have studied the effects due to rotation and stratification in isolation. We then looked at the effects of rotation in a barotropic model, but what about if we add stratification
More informationStrauss PDEs 2e: Section Exercise 4 Page 1 of 6
Strauss PDEs 2e: Section 5.3 - Exercise 4 Page of 6 Exercise 4 Consider the problem u t = ku xx for < x < l, with the boundary conditions u(, t) = U, u x (l, t) =, and the initial condition u(x, ) =, where
More informationDiagnosing Cyclogenesis by Partitioning Energy and Potential Enstrophy in a Linear Quasi-Geostrophic Model
Diagnosing Cyclogenesis by Partitioning Energy and Potential Enstrophy in a Linear Quasi-Geostrophic Model by Daniel Hodyss * and Richard Grotjahn Atmospheric Science Program, Dept. of Land, Air, and Water
More informationThe Generalized Eigenvalue-Eigenvector Problem or How to Simultaneously Diagonalize Two Quadratic Forms
The Generalized Eigenvalue-Eigenvector Problem or How to Simultaneously Diagonalize Two Quadratic orms The following records deal with the generalized eigenvalue problem of two positive definite matrices
More informationLecture 4 Eigenvalue problems
Lecture 4 Eigenvalue problems Weinan E 1,2 and Tiejun Li 2 1 Department of Mathematics, Princeton University, weinan@princeton.edu 2 School of Mathematical Sciences, Peking University, tieli@pku.edu.cn
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More information