Flow Instability Matlab Tutorial Bangalore, January 2010 Matthew Juniper, Dan Henningson, Peter Schmid 1
|
|
- Ophelia Grant
- 5 years ago
- Views:
Transcription
1 Flow Instability Matlab Tutorial Bangalore, January Matthew Juniper, Dan Henningson, Peter Schmid Simple differentiation matrices Imagine a vector that represents the values of a function y at N equispaced points between < x <. We can create a simple matrix that, at each point, approximates d/dx of this vector:. prog.m function [D] = prog (N) % Create a nd order central difference matrix, D, that differentiates the % vector on which it operates. The vector represents the values of y at N % equispaced points between - < x <. % % Operate this matrix on a sin wave and compare the result with the % analytical solution. 3 N = 7 y(x) y (x) D y(x) To do For small N, e.g. N = 7, run with the command line D = prog (N). Examine matrix D to understand why it differentiates the vector on which it acts. Vary N and run with the command line D = prog (N); to examine its behaviour.. prog.m Repeat this for a th order central difference matrix, D. This tutorial is based on one given at JNCASR by Peter Schmid in January 9
2 3 N = 7 y(x) y (x) D y(x) To do Same as for D = prog (N). Examine matrix D to understand why it is more accurate than that in prog. Vary N and compare with the results of prog..3 prog 5.m Repeat this for a Chebyshev differentiation matrix 3 N = 7 y(x) y (x) D y(x) To do For small N, e.g. N = 7, run with the command line D = prog 5(N). Examine matrix D and note the differences with those from prog and prog. Vary N and compare with the previous results. Eigenvalue problems. prog.m We will now solve a simple eigenvalue problem: d u dx + λu = with homogenous boundary conditions (u = at x = ±). This is the solution for a wave on a string that is held at both ends. The first eigenvalue corresponds to the fundamental
3 mode. Higher eigenvalues correspond to higher harmonics. % Calculate the Chebyshev matrix and the nd order differential operator matrix. [D,x] = cheb(n); D = D^; % Apply the homogenous boundary conditions D = D(:N,:N); % Find the eigenvalues and corresponding eigenvectors [V,Lam] = eig(-d); eig =.*/pi eig =.*/pi eig 3 = 9.*/pi eig =.5*/pi eig 5 = */pi eig = */pi eig 7 = */pi eig 8 =.8335*/pi eig 9 = 8.799*/pi To do Look at how the second order difference matrix is created. Look at how the homogenous boundary conditions are applied. How many eigenvalues are there? What happens as you reduce N? Why? How would you change the boundary conditions to enforce a zero gradient at x = ± (harder)?. prog 7.m We will now solve another eigenvalue problem of the Sturm-Liouville type: ( d x dy ) + x dy dx dx dx = λy with boundary conditions y = at x = and x =. 3
4 To do Watch the eigenvalues and eigenfunctions change as you increase N..3 prog 9.m We will now solve the Orr-Sommerfeld eigenvalue problem for planar Poiseuille flow, assuming that the perturbations are two-dimensional. (Later on we will allow for perturbations to vary in the direction into the page, which introduces the Squire modes.) This problem is adapted from p.m in Spectral Methods in Matlab by Lloyd N. Trefethen (SIAM ). The Orr-Sommerfeld equation is: ( d R ) ( u dx u d d dx + u iu i( x ) ( u d ) ) dx u u = λ dx u with boundary conditions: du = u = at x = ± dx We express this as a generalized eigenvalue problem: Aû = λbû We plot the eigenvalues in the complex plane. In this formulation, the eigenvalue with the highest real part is the most unstable / least stable. In this formulation, the axial wavenumber has been taken to be. The figures below are for Chebyshev points and are created by typing prog 9().
5 eigenvalues for N =.9.8 Highlighted eigenvector.7. imag(λ).. amplitude real(λ).5.5 x To do Click on eigenvalues to plot the corresponding eigenvector. Investigate the different branches of the eigenvalue plot. What happens as you increase N? Change the Reynolds number, R, by editing the program. 3 The Orr-Sommerfeld-Squire system This section is a practical demonstration of parts of. of Peter Schmid s 7 Annual Review paper, Schmid (7). We will start by considering planar Poiseuille flow at Reynolds number R. Unlike in the previous section, we will allow the perturbation to vary in the direction into the page (spanwise). We consider waves in the streamwise direction with wavenumber α and waves in the spanwise direction with wavenumber β. 3. prog.m We will look at planar Poiseuille flow first because it can be compared with the previous section. prog.m calculates the Orr-Sommerfeld-Squire matrix, which is matrix A in equation (.a) of Schmid (7). It plots out the eigenvalues of this matrix and the Orr-Sommerfeld eigenvector and Squire eigenvector of the eigenvalue that is clicked on by the user. e.g. type [OSSQ] = prog (,.,.,577);. function [OSSQ] = prog (N,alp,beta,Re) % Calculate the Orr-Sommerfeld-Squire matrix for planar Poiseuille flow. % % N : Number of Chebyshev points % alp : streamwise wavenumber % beta : spanwise wavenumber % Re : Reynolds number 5
6 eigenvalue number =. Orr Sommerfeld eigenfunction.3. amplitude amplitude.5.5 x Squire eigenfunction.5.5 x To do Compare the eigenvalue plot with that from prog 9.m at the same value of N. Click on eigenvalues to plot the corresponding eigenvectors. Identify the Orr-Sommerfeld eigenvalues and the Squire eigenvalues. Change the Reynolds number and N. prog.m returns the matrix OSSQ. Create your own routine to plot the eigenvalues of OSSQ. Watch out for the scale on your graph. 3. PoiMatrix.m Some of the Orr-Sommerfeld-Squire eigenvalues have little influence on the behaviour of the system and can safely be discarded. The function PoiMatrix.m cuts out the modes outside a defined window. It then calls function OSSQMatrix, which computes the matrix, Q. This matrix has the same eigenvalues as the reduced OSSQ matrix and is used to calculate transient growth. This is on the Maximum Amplification slide of Dan Henningson s notes. This can be run with the command line [Q,invF] = PoiMatrix(,.,.,577);, for instance. In the rest of this tutorial, we will used this reduced matrix, Q. To do Plot the eigenvalues of Q with your own routine. Compare the eigenvalue plot with that of prog 9.m at the same value of N. (Optional: plot out some of the corresponding eigenfunctions). Normal and non-normal matrices This section looks at two simple matrices, one of which is normal (i.e. has orthogonal eigenvectors) and the other of which is non-normal: A n = 5 i.3 + i A nn = 5 i.3 + i
7 . The numerical range The numerical range is defined on page 3 of Schmid (7). The maximum protrusion of the numerical range into the unstable half-plane determines the maximum energy growth at t = + i.e. it is a measure of the growth rate at very small times. Eigenvalues, on the other hand, give a measure of the growth or decay at large times. The program AnPlot.m (A normal plot) plots the eigenvalues and numerical range of matrix A n. The numerical range is a set of straight lines joining the eigenvalues. It therefore does not extend further into the right half plane than the right-most eigenvalue. This means that the growth rate at small times is the same as that at large times. The program AnnPlot.m (A non-normal plot) plots the eigenvalues and numerical range of matrix A nn. The numerical range is loop around the eigenvalues and extends significantly into the right half plane. This means that the growth rate at small times will be higher than that given by the right-most eigenvalue alone. a normal matrix unstable region eigenvalues numerical range a non normal matrix unstable region eigenvalues numerical range λ i λ i 8 λ r 8 λ r To do range. Make up your own non-normal matrices and plot their eigenvalues and numerical. Pseudospectra, a.k.a. the resolvent norm Pseudospectra and the resolvent norm are defined four different but equivalent ways in Trefethen and Embree () and on slide Pseudospectra, resolvents and sensitivity of Dan Henningson s notes. When a random matrix E with norm ɛ is added to a matrix L, the eigenvalues of L move slightly. The maximum amount that they will move is given by the pseudospectra of L. They will not move outside the pseudospectra contour level at height ɛ. The amount that the pseudospectra extend into the right half plane gives an indication of the maximum transient growth possible with matrix L. (For a clearer definition, see the first three chapters of Trefethen and Embree or the video of Nick Trefethen s talk at the first AIM Meeting, available on the AIM website.) 7
8 The program RAnPlot.m plots the pseudospectra of A n and RAnnPlot.m plots the pseudospectra of A nn. The contours show log of the resolvent norm at that value of σ. pseudospectra of a normal matrix pseudospectra of a non normal matrix 8 8 The programs RAnPlot shotgun.m and RAnnPlot shotgun.m plot the pseudospectra of A n and A nn followed by the eigenvalues of matrices that have been perturbed by an amount ɛ. In the figures below, ɛ = log ().. pseudospectra of a normal matrix pseudospectra of a non normal matrix To do Vary ɛ and re-run RAnPlot shotgun.m and RAnnPlot shotgun.m. Repeat for your own non-normal matrices..3 Transient growth If a matrix, A, is non-normal, it can exhibit transient growth. The maximum possible transient growth at time T is called G(T ) and is the maximum singular value of the matrix exponential of A T. Transient growth can occur whether or not the matrix has unstable eigenvalues. The function TransientGrowth simple.m calculates G(T ) for the range specified by the user. Try TransientGrowth simple(a,[ 5],). 8
9 To do Calculate the transient growth for a non-normal matrix with stable eigenvalues. Change the non-normality of the matrix. See what happens to the transient growth. 5 Analysis of the Orr-Sommerfeld-Squire matrix You now have all the tools you need to investigate the Orr-Sommerfeld-Squire matrix with Poiseuille flow as the base flow. Try combining section with section 3. To get you started, the program start.m creates this figure: Eigenvalues, pseudospectra and numerical range for Poiseuille flow The Couette directory contains functions that generate the Orr-Sommerfeld-Squire matrix for Couette flow. To generate the matrix, use the command [Q,invF] = CouMatrix(N,alp,beta,Re). The Pipe directory contains functions that generate the axisymmetric Orr-Sommerfeld-Squire matrix for round pipe flow. If you have time, try creating the Orr-Sommerfeld-Squire matrix for your own base flows. 9
SG2221 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 informationTransition to turbulence in plane Poiseuille flow
Proceedings of the 55th Israel Annual Conference on Aerospace Sciences, Tel-Aviv & Haifa, Israel, February 25-26, 2015 ThL2T5.1 Transition to turbulence in plane Poiseuille flow F. Roizner, M. Karp and
More informationOptimal Control of Plane Poiseuille Flow
Optimal Control of Plane Poiseuille Flow Workshop on Flow Control, Poitiers 11-14th Oct 04 J. Mckernan, J.F.Whidborne, G.Papadakis Cranfield University, U.K. King s College, London, U.K. Optimal Control
More informationLinear Instability of asymmetric Poiseuille flows
Report no. 7/ Linear Instability of asymmetric Poiseuille flows Dick Kachuma & Ian J. Sobey We compute solutions for the Orr-Sommerfeld equations for the case of an asymmetric Poiseuille-like parallel
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 informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow Fulvio Martinelli, Maurizio Quadrio, John McKernan and James F. Whidborne Abstract Subcritical flows may experience large
More informationStability Analysis of a Plane Poisseuille Flow. Course Project: Water Waves and Hydrodynamic Stability 2016 Matthias Steinhausen
Stability Analysis of a Plane Poisseuille Flow Course Project: Water Waves and Hydrodynamic Stability 2016 Matthias Steinhausen y Plane Poiseuille Flow (PPF): Characteristics y 1 0.5 Velocity profile for
More information25 Self-Assessment. 26 Eigenvalue Stability for a Linear ODE. Chapter 7 Eigenvalue Stability
Chapter 7 Eigenvalue Stability As we have seen, while numerical methods can be convergent, they can still exhibit instabilities as the number of timesteps n increases for finite t. For example, when applying
More informationÉCOLE DE PHYSIQUE DES HOUCHES. Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at
ÉCOLE DE PHYSIQUE DES HOUCHES STABILITY OF FLUID FLOWS Carlo COSSU Institut de Mécanique des Fluides Toulouse CNRS Université de Toulouse, also at Département de Mécanique, École Polytechinique http://www.imft.fr/carlo-cossu/
More informationLinear stability of channel entrance flow
Linear stability of channel entrance flow Damien Biau DIAM, University of Genova, Via Montallegro 1, 16145 Genova, Italy Abstract The spatial stability analysis of two dimensional, steady channel flow
More informationTransient growth on boundary layer streaks. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Transient growth on boundary layer streaks Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Primary instability: TS vs TG Secondary instability of the streaks Movie of streaks
More informationStability of Shear Flow
Stability of Shear Flow notes by Zhan Wang and Sam Potter Revised by FW WHOI GFD Lecture 3 June, 011 A look at energy stability, valid for all amplitudes, and linear stability for shear flows. 1 Nonlinear
More informationStochastic excitation of streaky boundary layers. Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden
Stochastic excitation of streaky boundary layers Jérôme Hœpffner Luca Brandt, Dan Henningson Department of Mechanics, KTH, Sweden Boundary layer excited by free-stream turbulence Fully turbulent inflow
More informationLINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP
MAGNETOHYDRODYNAMICS Vol. 48 (2012), No. 1, pp. 147 155 LINEAR STABILITY ANALYSIS FOR THE HARTMANN FLOW WITH INTERFACIAL SLIP N. Vetcha, S. Smolentsev, M. Abdou Fusion Science and Technology Center, UCLA,
More informationInitial-value problem for shear flows
Initial-value problem for shear flows 1 Mathematical Framework The early transient and long asymptotic behaviour is studied using the initialvalue problem formulation for two typical shear flows, the plane
More informationNon-modal disturbances growth in a viscous mixing layer flow
Non-modal disturbances growth in a viscous miing layer flow Helena Vitoshkin and Aleander Yu. Gelfgat School of Mechanical Engineering, Faculty of Engineering, Tel-Aviv University, Ramat Aviv, Tel-Aviv
More informationGeneral introduction to Hydrodynamic Instabilities
KTH ROYAL INSTITUTE OF TECHNOLOGY General introduction to Hydrodynamic Instabilities L. Brandt & J.-Ch. Loiseau KTH Mechanics, November 2015 Luca Brandt Professor at KTH Mechanics Email: luca@mech.kth.se
More informationAnalysis of fluid systems: stability, receptivity, sensitivity
Analysis of fluid systems: stability, receptivity, sensitivity Peter J. Schmid Department of Mathematics Imperial College London London, SW7 AZ, United Kingdom peter.schmid@imperial.ac.uk Luca Brandt Linné
More informationNon-modal stability analysis of stratified two-phase channel flows
Non-modal stability analysis of stratified two-phase channel flows I. Barmak a), A. Yu. Gelfgat, A. Ullmann, and N. Brauner School of Mechanical Engineering, Tel Aviv University, Tel Aviv 69978, Israel
More informationOrdinary Differential Equations II
Ordinary Differential Equations II February 23 2017 Separation of variables Wave eq. (PDE) 2 u t (t, x) = 2 u 2 c2 (t, x), x2 c > 0 constant. Describes small vibrations in a homogeneous string. u(t, x)
More informationTransient growth for streak-streamwise vortex interactions
Physics Letters A 358 006) 431 437 www.elsevier.com/locate/pla Transient growth for streak-streamwise vortex interactions Lina Kim, Jeff Moehlis Department of Mechanical Engineering, University of California,
More informationPhysics 250 Green s functions for ordinary differential equations
Physics 25 Green s functions for ordinary differential equations Peter Young November 25, 27 Homogeneous Equations We have already discussed second order linear homogeneous differential equations, which
More informationModal and non-modal stability of particle-laden channel flow. Abstract
Modal and non-modal stability of particle-laden channel flow Joy Klinkenberg,, 2 H. C. de Lange, and Luca Brandt 2 TU/e, Mechanical Engineering, 56 MB, Eindhoven, The Netherlands 2 Linné Flow Centre, KTH
More informationNonmodal stability analysis
Nonmodal stability analysis Peter Schmid Laboratoire d Hydrodynamique (LadHyX) CNRS Ecole Polytechnique, Palaiseau Motivation: bring together expertise in combustion physics, fluid dynamics, stability
More informationLecture IX. Definition 1 A non-singular Sturm 1 -Liouville 2 problem consists of a second order linear differential equation of the form.
Lecture IX Abstract When solving PDEs it is often necessary to represent the solution in terms of a series of orthogonal functions. One way to obtain an orthogonal family of functions is by solving a particular
More informationMATH 552 Spectral Methods Spring Homework Set 5 - SOLUTIONS
MATH 55 Spectral Methods Spring 9 Homework Set 5 - SOLUTIONS. Suppose you are given an n n linear system Ax = f where the matrix A is tridiagonal b c a b c. A =.........,. a n b n c n a n b n with x =
More informationEnergy Transfer Analysis of Turbulent Plane Couette Flow
Energy Transfer Analysis of Turbulent Plane Couette Flow Satish C. Reddy! and Petros J. Ioannou2 1 Department of Mathematics, Oregon State University, Corvallis, OR 97331 USA, reddy@math.orst.edu 2 Department
More informationADVANCED ENGINEERING MATHEMATICS MATLAB
ADVANCED ENGINEERING MATHEMATICS WITH MATLAB THIRD EDITION Dean G. Duffy Contents Dedication Contents Acknowledgments Author Introduction List of Definitions Chapter 1: Complex Variables 1.1 Complex Numbers
More informationStability of Plane Couette Flow and Pipe Poiseuille Flow PER-OLOV ÅSÉN
Stability of Plane Couette Flow and Pipe Poiseuille Flow PER-OLOV ÅSÉN Doctoral Thesis Stockholm, Sweden, 2007 TRITA-CSC-A 2007:7 ISSN 1653-5723 ISRN KTH/CSC/A--07/07--SE ISBN 978-91-7178-651-7 CSC Skolan
More informationRésonance et contrôle en cavité ouverte
Résonance et contrôle en cavité ouverte Jérôme Hœpffner KTH, Sweden Avec Espen Åkervik, Uwe Ehrenstein, Dan Henningson Outline The flow case Investigation tools resonance Reduced dynamic model for feedback
More informationModeling flow statistics using convex optimization
Modeling flow statistics using convex optimization Abstract A method is proposed to estimate the covariance of disturbances to a stable linear system when its state covariance is known and a dynamic model
More informationMünsteranian Torturials on Nonlinear Science. edited by Uwe Thiele, Oliver Kamps, Svetlana Gurevich. Continuation
Münsteranian Torturials on Nonlinear Science edited by Uwe Thiele, Oliver Kamps, Svetlana Gurevich Continuation SH and SHPER: Steady states of the Swift-Hohenberg equation Uwe Thiele with the support of
More informationFeedback control of transient energy growth in subcritical plane Poiseuille flow
Feedback control of transient energy growth in subcritical plane Poiseuille flow F. Martinelli, M. Quadrio, J. McKernan, J.F. Whidborne LadHyX, École Polytechnique; Politecnico di Milano; King s College,
More informationFlow Transition in Plane Couette Flow
Flow Transition in Plane Couette Flow Hua-Shu Dou 1,, Boo Cheong Khoo, and Khoon Seng Yeo 1 Temasek Laboratories, National University of Singapore, Singapore 11960 Fluid Mechanics Division, Department
More informationTHE INITIAL-VALUE PROBLEM FOR VISCOUS CHANNEL FLOWS 1. W.O. Criminale. University ofwashington. Seattle, Washington T.L.
THE INITIAL-VALUE PROBLEM FOR VISCOUS CHANNEL FLOWS 1 W.O. Criminale Department of Applied Mathematics University ofwashington Seattle, Washington 98195 T.L. Jackson Institute for Computer Applications
More informationTransient growth and minimal defects: Two possible initial paths of transition to turbulence in plane shear flows
PHYSICS OF FLUIDS VOLUME 16, NUMBER 10 OCTOBER 2004 Transient growth and minimal defects: Two possible initial paths of transition to turbulence in plane shear flows Damien Biau a) and Alessandro Bottaro
More informationChapter 1, Section 1.2, Example 9 (page 13) and Exercise 29 (page 15). Use the Uniqueness Tool. Select the option ẋ = x
Use of Tools from Interactive Differential Equations with the texts Fundamentals of Differential Equations, 5th edition and Fundamentals of Differential Equations and Boundary Value Problems, 3rd edition
More informationAND NONLINEAR SCIENCE SERIES. Partial Differential. Equations with MATLAB. Matthew P. Coleman. CRC Press J Taylor & Francis Croup
CHAPMAN & HALL/CRC APPLIED MATHEMATICS AND NONLINEAR SCIENCE SERIES An Introduction to Partial Differential Equations with MATLAB Second Edition Matthew P Coleman Fairfield University Connecticut, USA»C)
More informationRobust control of resistive wall modes using pseudospectra
Robust control of resistive wall modes using pseudospectra M. Sempf, P. Merkel, E. Strumberger, C. Tichmann, and S. Günter Max-Planck-Institut für Plasmaphysik, EURATOM Association, Garching, Germany GOTiT
More informationConvective instability and transient growth in flow over a backwardfacing
Convective instability and transient growth in flow over a backwardfacing step Journal: Manuscript ID: mss type: Date Submitted by the Author: Complete List of Authors: Keyword: Journal of Fluid Mechanics
More information1 Singular Value Decomposition and Principal Component
Singular Value Decomposition and Principal Component Analysis In these lectures we discuss the SVD and the PCA, two of the most widely used tools in machine learning. Principal Component Analysis (PCA)
More informationIntroduction to Mathematical Physics
Introduction to Mathematical Physics Methods and Concepts Second Edition Chun Wa Wong Department of Physics and Astronomy University of California Los Angeles OXFORD UNIVERSITY PRESS Contents 1 Vectors
More information4 Shear flow instabilities
4 Shear flow instabilities Here we consider the linear stability of a uni-directional base flow in a channel u B (y) u B = 0 in the region y [y,y 2. (35) 0 y = y 2 y z x u (y) B y = y In Sec. 4. we derive
More informationReduced-Order Models for Feedback Control of Transient Energy Growth
9th AIAA Flow Control Conference, 5 9 June, Atlanta, Georgia. -9 Reduced-Order Models for Feedback Control of Transient Energy Growth Aniketh Kalur and Maziar S. Hemati Aerospace Engineering and Mechanics,
More informationTransient development of perturbations in a stratified turbulent shear flow
Institute of Thermal Turbomachinery Prof. Dr.-Ing. H.-J. Bauer Transient development of perturbations in a stratified turbulent shear flow Diploma Thesis of Carlos Yáñez Vico Supervisor: Dr.-Ing. Manuel
More informationDynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario
Dynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario NDAMS Workshop @ YITP 1 st November 2011 Meheboob Alam and Priyanka Shukla Engineering Mechanics Unit
More informationLinear stability of plane Poiseuille flow over a generalized Stokes layer
Linear stability of plane Poiseuille flow over a generalized Stokes layer Maurizio Quadrio 1,, Fulvio Martinelli and Peter J Schmid 1 Dip. Ing. Aerospaziale, Politecnico di Milano, Campus Bovisa, I-0156
More information1 Sylvester equations
1 Sylvester equations Notes for 2016-11-02 The Sylvester equation (or the special case of the Lyapunov equation) is a matrix equation of the form AX + XB = C where A R m m, B R n n, B R m n, are known,
More informationAccurate estimation for the closed-loop robust control using global modes
1/21 G. Tissot, GDR CDD 216. 1/21 Accurate estimation for the closed-loop robust control using global modes application to the Ginzburg-Landau equation Gilles Tissot, Jean-Pierre Raymond Institut de Mathématiques
More informationOptimal Control of Plane Poiseuille Flow
Optimal Control of Plane Poiseuille Flow Dr G.Papadakis Division of Engineering, King s College London george.papadakis@kcl.ac.uk AIM Workshop on Flow Control King s College London, 17/18 Sept 29 AIM workshop
More informationConvective instability and transient growth in flow over a backward-facing step
J. Fluid Mech. (28), vol. 63, pp. 271 34. c 28 Cambridge University Press doi:1.117/s221128119 Printed in the United Kingdom 271 Convective instability and transient growth in flow over a backward-facing
More informationNumerically computing zeros of the Evans Function
Numerically computing zeros of the Evans Function Rebekah Coggin Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty Advisor: Todd Kapitula Department of Mathematics
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 informationGENERALIZED STABILITY OF THE TWO-LAYER MODEL
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
More informationProblem set 3: Solutions Math 207B, Winter Suppose that u(x) is a non-zero solution of the eigenvalue problem. (u ) 2 dx, u 2 dx.
Problem set 3: Solutions Math 27B, Winter 216 1. Suppose that u(x) is a non-zero solution of the eigenvalue problem u = λu < x < 1, u() =, u(1) =. Show that λ = (u ) 2 dx u2 dx. Deduce that every eigenvalue
More informationEnergy amplification in channel flows of viscoelastic fluids
J. Fluid Mech. (2008), vol. 601, pp. 407 424. c 2008 Cambridge University Press doi:10.1017/s0022112008000633 Printed in the United Kingdom 407 Energy amplification in channel flows of viscoelastic fluids
More informationSpecial Functions of Mathematical Physics
Arnold F. Nikiforov Vasilii B. Uvarov Special Functions of Mathematical Physics A Unified Introduction with Applications Translated from the Russian by Ralph P. Boas 1988 Birkhäuser Basel Boston Table
More information(Linear equations) Applied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB (Linear equations) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots
More informationSystems of Linear Differential Equations Chapter 7
Systems of Linear Differential Equations Chapter 7 Doreen De Leon Department of Mathematics, California State University, Fresno June 22, 25 Motivating Examples: Applications of Systems of First Order
More informationStability of circular Poiseuille Couette flow to axisymmetric disturbances
J. Fluid Mech. 2004), vol. 500, pp. 169 210. c 2004 Cambridge University Press DOI: 10.1017/S0022112003007158 Printed in the United Kingdom 169 Stability of circular Poiseuille Couette flow to axisymmetric
More informationDegree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m.
Degree Master of Science in Mathematical Modelling and Scientific Computing Mathematical Methods I Thursday, 12th January 2012, 9:30 a.m.- 11:30 a.m. Candidates should submit answers to a maximum of four
More informationThe issue of the relevance of modal theory to rectangular duct instability is discussed
1 Supplemental Appendix 4: On global vs. local eigenmodes Rectangular duct vs. plane Poiseuille flow The issue of the relevance of modal theory to rectangular duct instability is discussed in the main
More informationWILEY. Differential Equations with MATLAB (Third Edition) Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg
Differential Equations with MATLAB (Third Edition) Updated for MATLAB 2011b (7.13), Simulink 7.8, and Symbolic Math Toolbox 5.7 Brian R. Hunt Ronald L. Lipsman John E. Osborn Jonathan M. Rosenberg All
More informationMAE 210C FLUID MECHANICS III SPRING 2015 Rayleigh s equation for axisymmetric flows Forthestabilityanalysisofparallelaxisymmetricflowswithbasicvelocity profileū = (U x,u r,u θ ) = (U(r),0,0) weintroducemodalperturbationsoftheform(v
More informationMath 46, Applied Math (Spring 2009): Final
Math 46, Applied Math (Spring 2009): Final 3 hours, 80 points total, 9 questions worth varying numbers of points 1. [8 points] Find an approximate solution to the following initial-value problem which
More informationReduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD
3rd AIAA Flow Control Conference, 5 8 June 26, San Francisco Reduced-Order Modeling of Channel Flow Using Traveling POD and Balanced POD M. Ilak and C. W. Rowley Dept. of Mechanical and Aerospace Engineering,
More informationChapter #4 EEE8086-EEE8115. Robust and Adaptive Control Systems
Chapter #4 Robust and Adaptive Control Systems Nonlinear Dynamics.... Linear Combination.... Equilibrium points... 3 3. Linearisation... 5 4. Limit cycles... 3 5. Bifurcations... 4 6. Stability... 6 7.
More informationApplied Linear Algebra
Applied Linear Algebra Peter J. Olver School of Mathematics University of Minnesota Minneapolis, MN 55455 olver@math.umn.edu http://www.math.umn.edu/ olver Chehrzad Shakiban Department of Mathematics University
More informationSection 1.1 Algorithms. Key terms: Algorithm definition. Example from Trapezoidal Rule Outline of corresponding algorithm Absolute Error
Section 1.1 Algorithms Key terms: Algorithm definition Example from Trapezoidal Rule Outline of corresponding algorithm Absolute Error Approximating square roots Iterative method Diagram of a general iterative
More informationLaboratory handout 5 Mode shapes and resonance
laboratory handouts, me 34 82 Laboratory handout 5 Mode shapes and resonance In this handout, material and assignments marked as optional can be skipped when preparing for the lab, but may provide a useful
More informationAkademisk avhandling som med tillstνand av Kungl Tekniska Högskolan framlägges till offentlig granskning för avläggande av teknisk doktorsexamen freda
Stability results for viscous shock waves and plane Couette flow Mattias Liefvendahl Stockholm 2001 Doctoral Dissertation Royal Institute of Technology Department of Numerical Analysis and Computer Science
More informationOn an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University
On an Eigenvalue Problem Involving Legendre Functions by Nicholas J. Rose North Carolina State University njrose@math.ncsu.edu 1. INTRODUCTION. The classical eigenvalue problem for the Legendre Polynomials
More informationHopf bifurcations to quasi-periodic solutions for the two-dimensional plane Poiseuille flow
Hopf bifurcations to quasi-periodic solutions for the two-dimensional plane Poiseuille flow Pablo S Casas Departamento de Matemática Aplicada I, Universidad Politécnica de Cataluña, Diagonal, 647 08028
More informationThe need for de-aliasing in a Chebyshev pseudo-spectral method
The need for de-aliasing in a Chebyshev pseudo-spectral method Markus Uhlmann Potsdam Institut für Klimafolgenforschung, D-442 Potsdam uhlmann@pik-potsdam.de (March 2) Abstract In the present report, we
More informationModal Decomposition Methods on Aerodynamic Flows
498 Modal Decomposition Methods on Aerodynamic Flows Mohammad Niaz Murshed University of Michigan Department of Aerospace Engineering February 5, 2018 Abstract- This paper discusses the significance of
More information4 Power Series Solutions: Frobenius Method
4 Power Series Solutions: Frobenius Method Now the ODE adventure takes us to series solutions for ODEs, a technique A & W, that is often viable, valuable and informative. These can be readily applied Sec.
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationMATH 333: Partial Differential Equations
MATH 333: Partial Differential Equations Problem Set 9, Final version Due Date: Tues., Nov. 29, 2011 Relevant sources: Farlow s book: Lessons 9, 37 39 MacCluer s book: Chapter 3 44 Show that the Poisson
More informationNon-Self-Adjoint Operators and Pseudospectra
Non-Self-Adjoint Operators and Pseudospectra E.B. Davies King s College London 1 The Relevance of Pseudospectra The theory of pseudospectra is a way of imposing some order on the chaos of non-self-adjoint
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More information2 Introduction to perturbation methods
2 Introduction to perturbation methods 2.1 What are perturbation methods? Perturbation methods are methods which rely on there being a dimensionless parameter in the problem that is relatively small: ε
More information1 A complete Fourier series solution
Math 128 Notes 13 In this last set of notes I will try to tie up some loose ends. 1 A complete Fourier series solution First here is an example of the full solution of a pde by Fourier series. Consider
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 informationv(x, 0) = g(x) where g(x) = f(x) U(x). The solution is where b n = 2 g(x) sin(nπx) dx. (c) As t, we have v(x, t) 0 and u(x, t) U(x).
Problem set 4: Solutions Math 27B, Winter216 1. The following nonhomogeneous IBVP describes heat flow in a rod whose ends are held at temperatures u, u 1 : u t = u xx < x < 1, t > u(, t) = u, u(1, t) =
More informationMATH 2203 Exam 3 Version 2 Solutions Instructions mathematical correctness clarity of presentation complete sentences
MATH 2203 Exam 3 (Version 2) Solutions March 6, 2015 S. F. Ellermeyer Name Instructions. Your work on this exam will be graded according to two criteria: mathematical correctness and clarity of presentation.
More informationHigh Speed Couette Poiseuille Flow Stability in Reverse Flow Conditions
High Speed Couette Poiseuille Flow Stability in Conditions ALĐ ASLAN EBRĐNÇ Engine & Power Train System Engineering Ford Otosan - Kocaeli TURKEY SERKAN ÖZGEN Aeronatical and Aeronatics Engineering METU
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationTURBULENT BOUNDARY LAYERS: CURRENT RESEARCH AND FUTURE DIRECTIONS
TURBULENT BOUNDARY LAYERS: CURRENT RESEARCH AND FUTURE DIRECTIONS Beverley J. McKeon Graduate Aerospace Laboratories California Institute of Technology http://mckeon.caltech.edu * AFOSR #s FA9550-08-1-0049,
More informationA fast and well-conditioned spectral method: The US method
A fast and well-conditioned spectral method: The US method Alex Townsend University of Oxford Leslie Fox Prize, 24th of June 23 Partially based on: S. Olver & T., A fast and well-conditioned spectral method,
More informationPerturbation Energy Production in Pipe Flow over a Range of Reynolds Numbers using Resolvent Analysis
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition - 8 January 29, Orlando, Florida AIAA 29-3 Perturbation Energy Production in Pipe Flow over a Range of Reynolds
More informationModal and non-modal stability analysis of a channel flow seeded with light particles
Modal and non-modal stability analysis of a channel flow seeded with light particles By Joy Klinkenberg 1,2, H.C. de Lange 1 and Luca Brandt 2 1 TU/e, Mechanical Engineering, 5600 MB, Eindhoven, The Netherlands
More informationForced Mechanical Vibrations
Forced Mechanical Vibrations Today we use methods for solving nonhomogeneous second order linear differential equations to study the behavior of mechanical systems.. Forcing: Transient and Steady State
More informationRevision of the thermal chip simulation, J.E. Akin, Rice University
Revision of the thermal chip simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 3: Wednesday, Jan 9
Problem du jour Week 3: Wednesday, Jan 9 1. As a function of matrix dimension, what is the asymptotic complexity of computing a determinant using the Laplace expansion (cofactor expansion) that you probably
More informationProblem Set 8 Mar 5, 2004 Due Mar 10, 2004 ACM 95b/100b 3pm at Firestone 303 E. Sterl Phinney (2 pts) Include grading section number
Problem Set 8 Mar 5, 24 Due Mar 1, 24 ACM 95b/1b 3pm at Firestone 33 E. Sterl Phinney (2 pts) Include grading section number Useful Readings: For Green s functions, see class notes and refs on PS7 (esp
More informationNotes on Eigenvalues, Singular Values and QR
Notes on Eigenvalues, Singular Values and QR Michael Overton, Numerical Computing, Spring 2017 March 30, 2017 1 Eigenvalues Everyone who has studied linear algebra knows the definition: given a square
More informationMathematical Modeling using Partial Differential Equations (PDE s)
Mathematical Modeling using Partial Differential Equations (PDE s) 145. Physical Models: heat conduction, vibration. 146. Mathematical Models: why build them. The solution to the mathematical model will
More information18.06 Problem Set 9 - Solutions Due Wednesday, 21 November 2007 at 4 pm in
8.6 Problem Set 9 - Solutions Due Wednesday, 2 November 27 at 4 pm in 2-6. Problem : (5) When A = SΛS is a real-symmetric (or Hermitian) matrix, its eigenvectors can be chosen orthonormal and hence S =
More informationDifferential Equations
Differential Equations Theory, Technique, and Practice George F. Simmons and Steven G. Krantz Higher Education Boston Burr Ridge, IL Dubuque, IA Madison, Wl New York San Francisco St. Louis Bangkok Bogota
More information