Analysis Methods in Atmospheric and Oceanic Science
|
|
- Erik West
- 5 years ago
- Views:
Transcription
1 Analysis Methods in Atmospheric and Oceanic Science AOSC 65 Partial Differential Equations Week 1, Day 3 1 Nov 014 1
2 Partial Differential Equation An equation that defines the relationship between an unknown function (or functions) of several independent variables and their partial derivatives with respect to those variables
3 3 Partial Differential Equations u (,) xt = K (,) xt (1)
4 4 Partial Differential Equations u (,) xt = K (,) xt (1) What is u?
5 5 Partial Differential Equations What is u? u (,) xt = K (,) xt Heat: measure of total energy that depends on mass Temperature: measure of molecular motion that is independent of mass (1) See
6 6 Partial Differential Equations What is u? What is K? u (,) xt = K (,) xt Heat: measure of total energy that depends on mass Temperature: measure of molecular motion that is independent of mass (1) See
7 7 Partial Differential Equations u (,) xt = K (,) xt (1) Alternate forms of this equation: cρ ( xyzt,,, ) = kxyz (,, ) ( xyzt,,, ) + kxyz (,, ) ( xyzt,,, ) y y + kxyz (,, ) ( xyzt,,, ) () z z
8 8 Partial Differential Equations u (,) xt = K (,) xt (1) Alternate forms of this equation: cρ ( xyzt,,, ) = kxyz (,, ) ( xyzt,,, ) + kxyz (,, ) ( xyzt,,, ) y y + kxyz (,, ) ( xyzt,,, ) () z z c is specific heat ρ is density and k is thermal conductivity where
9 9 Partial Differential Equations u (,) xt = K (,) xt (1) Alternate forms of this equation: cρ ( xyzt,,, ) = kxyz (,, ) ( xyzt,,, ) + kxyz (,, ) ( xyzt,,, ) y y + kxyz (,, ) ( xyzt,,, ) () z z Equation can be re-written (and solved in a straight forward manner) if c, ρ, and k are constant wrt spatial dimensions: cρ u u u ( xyzt,,, ) = ( xyzt,,, ) + ( xyzt,,, ) + ( xyzt,,, ) k y z (3)
10 10 Partial Differential Equations Equation can be re-written (and solved in a straight forward manner) if c, ρ, and k are constant wrt spatial dimensions: c u u u u ρ ( xyzt,,, ) = ( xyzt,,, ) + ( xyzt,,, ) + ( xyzt,,, ) k y z (3) How might a physicist write this same equation?
11 Partial Differential Equations Equation can be re-written (and solved in a straight forward manner) if c, ρ, and k are constant wrt spatial dimensions: c u u u u ρ ( xyzt,,, ) = ( xyzt,,, ) + ( xyzt,,, ) + ( xyzt,,, ) k y z (3) How might a physicist write this same equation? T ρ Cp = Q + ( k T) where T is temperature, T is the gradient of temperature is the divergence operator ρ is the density C is the heat capacity and Q k (4) 11 p is the heat added to the system is the thermal conductivity
12 1 Partial Differential Equations u (,) xt = K (,) xt Another application of this same equation in AOSC that does not involve Temperature or Heat: (1) Carbon Carbon (z, t) = K ( zt, ) (5) z where Carbon is the Dissolved Inorganic Carbon content of the ocean z is the depth of the ocean K is the eddy diffusion coefficient
13 13 Partial Differential Equations A classic example of PDEs from oceanography: For the barotropic, tropical ocean, the linear inviscid governing equations (Gill, 198) are: βyv v + βyu v + y v y = g = g = 1 H + βy η η y η v + y (1(6a) a), (1(6b) b), (6c) (1c ), + βv = 0 (6d) (1d )
14 14 Partial Differential Equations A classic example of PDEs from atmospheric sciences: For the barotropic atmosphere the governing equations (Anderson and McCreary, 1985) are: U βyv V + βyu U V + y P = ru P = rv y = r ( (7a) a), ((7b) b), P ((7c) c), c
15 15 Partial Differential Equations Coupled ocean/atmosphere system: = + = + + =, 1,, 0 t H y v x u y g yu t v H x g yv t u x η η β ρ τ η β + = + = + =,,, c Q c P r y V x U rv y P yu t V ru x P yv t U β β The ocean The atmosphere
16 16 Partial Differential Equations One more example of a classic PDE from AOSC: θ Potential Vorticity: PV = g ( ζ θ + f ) p (8a) where: PV is potential vorticity f is Coriolis Parameter, Ω sin(latitude), where Ω is angular speed of Earth s rotation (s 1 ) g is gravitational acceleration (m s ) p is pressure (mbar, or force/area kg m s / m ) θ is potential temperature (K) ζ θ is relative isentropic vorticity (s 1 ) ζ θ v = evaluated on an isentropic surface y (8b) Units: K m s 1 kg 1. The quantity 10 6 K m s 1 kg 1 is often called MKS PV units
17 17 u (,) xt = K (,) xt Typically in atmospheric and oceanic sciences, PDEs are solved numerically rather than analytically: although in many cases analytic manipulation is key to understanding the solutions (e.g., wave propagation) or even being able to proceed with a numerical solution Three aspects of Numerical Solutions to PDEs that must be considered: 1. Initial Conditions and Boundary Conditions. Numerical Stability 3. Numerical Diffusion Note: material presented in this part of the lecture is based on Chapter 1 (Numerical Solutions to PDEs) of Numerical Analysis, 8 th Edition, by Richard L. Burden and J. Douglass Faires (1)
18 18 u (,) xt = K (,) xt Typically in atmospheric and oceanic sciences, PDEs are solved numerically rather than analytically: although in many cases analytic manipulation is key to understanding the solutions (e.g., wave propagation) or even being able to proceed with a numerical solution Three aspects of Numerical Solutions to PDEs that must be considered: 1. Initial Conditions and Boundary Conditions. Numerical Stability 3. Numerical Diffusion We ll consider all three aspects today: we re ambitious (1)
19 19 u (,) xt = K (,) xt 1. Initial Conditions (IC) and Boundary Conditions (BC) For the heat diffusion equation, must typically specify distribution of temperature (or heat) at all locations for t = 0 ux (,0) = f( x) Also, must specify either value of temperature (or heat) at boundary: 0 (1) (IC) u(0, t) = u and u( L, t) = u where the rod has length L (BCa) L or flux of temperature (or heat) at boundary: (0, t) = 0 and ( Lt, ) = 0 if the rod is insulated on both ends (BCb)
20 0 u (,) xt = K (,) xt 1. Initial Conditions (IC) and Boundary Conditions (BC) For the heat diffusion equation, must typically specify distribution of temperature (or heat) at all locations for t = 0 ux (,0) = f( x) Also, must specify either value of temperature (or heat) at boundary: 0 (1) (IC) u(0, t) = u and u( L, t) = u where the rod has length L (BCa) L or flux of temperature (or heat) at boundary: (0, t) = 0 and ( Lt, ) = 0 if the rod is insulated on both ends (BCb) Let s examine how our solution depends on IC and BC!
21 1 u (,) xt = K (,) xt (1) 1. Initial Conditions (IC) and Boundary Conditions (BC) For the heat diffusion equation, must typically specify distribution of temperature (or heat) at the all locations for t = 0 ux (,0) = f ( x ) (IC) Also, must specify either value of temperature (or heat) at boundary: u (0, t ) = u and u ( L, t ) = u where the rod has length L (BCa) 0 L or flux of temperature (or heat) at boundary: u u (0, t ) = 0 and ( Lt, ) = 0 if the rod is insulated (BCb) x x on both ends Let s examine how our solution depends on IC and BC! Can show that for BCa Can show that for BCb ul u lim uxt (, ) = u0 + x t L lim uxt (, ) = Constant t 0
22 1. Initial Conditions (IC) and Boundary Conditions (BC) For the heat diffusion equation, must typically specify distribution of temperature (or heat) at the all locations for t = 0 ux (,0) = f ( x ) (IC) Also, must specify either value of temperature (or heat) at boundary: u (0, t ) = u and u ( L, t ) = u where the rod has length L (BCa) 0 u (,) xt = K (,) xt (1) L or flux of temperature (or heat) at boundary: u u (0, t ) = 0 and ( Lt, ) = 0 if the rod is insulated (BCb) x x on both ends Let s examine how our solution depends on IC and BC! Can show that for BCa ul u0 lim uxt (, ) = u0 + x t L Can show that for BCb lim uxt (, ) = Constant What happened to the IC?!? t
23 3 u (,) xt = K (,) xt (1) u(0, t) = ult (, ) = 0, t> 0 and ux (,0) = f( x), 0 x L (IC) & (BCa). Numerical Stability Select a time-step size k ux ( i, tj) ux ( i, tj 1) k u ( xi, tj) + + k
24 4 u (,) xt = K (,) xt (1) u(0, t) = ult (, ) = 0, t> 0 and ux (,0) = f( x), 0 x L (IC) & (BCa). Numerical Stability Select a time-step size k ux ( i, tj) ux ( i, tj 1) k u ( xi, tj) + + k Select a spatial-step size h 4 u ux ( i + ht, j) ux ( i, tj) + ux ( i ht, j) h u ( x, ) i tj + h 1 4
25 5 u (,) xt = K (,) xt (1) u(0, t) = ult (, ) = 0, t> 0 and ux (,0) = f( x), 0 x L (IC) & (BCa). Numerical Stability Select a time-step size k ux ( i, tj) ux ( i, tj 1) k u ( xi, tj) + + k Select a spatial-step size h Then: 4 u ux ( i + ht, j) ux ( i, tj) + ux ( i ht, j) h u ( x, ) i tj + h 1 4 u u u u + u K k h i, j i, j 1 i + 1, j i, j i 1, j = 0
26 u (,) xt = K (,) xt (1) u(0, t) = ult (, ) = 0, t> 0 and ux (,0) = f( x), 0 x L (IC) & (BCa). Numerical Stability The relation: u u u u + u K k h i, j i, j 1 i + 1, j i, j i 1, j combined with the IC & BC imply: = 0 (1+ λ) λ 0 0 u u 1, j 1, j 1 λ (1 + λ) λ u, j u, j 1 λ (1 + λ) λ = λ λ 0 λ (1 + λ) um 1, j um 1, j 1 where λ = Kk ( / h) 6
27 u (,) xt = K (,) xt (1) u(0, t) = ult (, ) = 0, t> 0 and ux (,0) = f( x), 0 x L (IC) & (BCa). Numerical Stability The relation: u u u u + u K k h i, j i, j 1 i + 1, j i, j i 1, j combined with the IC & BC imply: 7 = 0 (1+ λ) λ 0 0 u u 1, j 1, j 1 λ (1 + λ) λ u, j u, j 1 λ (1 + λ) λ = λ λ 0 λ (1 + λ) um 1, j um 1, j 1 Another matrix equation, finally! where λ = Kk ( / h)
28 . Numerical Stability has everything to do with the properties of the matrix! It can be shown that, due to the properties of the eigenvalues of this matrix, it s inverse must exist: therefore, can always find values of u at time j based on knowledge of values at time j 1 (1+ λ) λ 0 0 u u 1, j 1, j 1 λ (1 + λ) λ u, j u, j 1 λ (1 + λ) λ = λ λ 0 λ (1 + λ) um 1, j um 1, j 1 where λ = Kk ( / h) 8
29 . Numerical Stability has everything to do with the properties of the matrix! It can be shown that, due to the properties of the eigenvalues of this matrix, it s inverse must exist: therefore, can always find values of u at time j based on knowledge of values at time j 1 This solution backward difference or implicit soln is said to be unconditionally stable Interestingly, had we defined: (1+ λ) λ 0 0 u1, j u1, j 1 we λ would (1 have λ) arrived λ at a different matrix that can only be inverted + when k, h, and K satisfy a specific relation. u This, j solution, u, j 1called the λ (1 + λ) λ = forward difference or explicit soln is said to be conditionally stable λ λ um 1, j um 1, j 1 0 λ (1 + λ) where λ = Kk ( / h) 9
30 30. Numerical Stability has everything to do with the properties of the matrix! It can be shown that, due to the properties of the eigenvalues of this matrix, it s inverse must exist: therefore, can always find values of u at time j based on knowledge of values at time j 1 This solution backward difference or implicit soln is said to be unconditionally stable Interestingly, had we defined: u ux ( i, tj+ 1) ux ( i, tj) k u ( xi, tj) + k we would have arrived at a different matrix that can only be inverted when k, h, and K satisfy a specific relation. This solution, called the forward difference or explicit soln is said to be conditionally stable
31 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Note: material presented in this part of the lecture is based on an excellent presentation by Anne Douglass (NASA GSFC) that was given at the University of Toronto, available on the web at: 31
32 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Physical solution 3
33 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Physical solution Upwind material mixes uniformly in neighboring grid box. Mathematical solution 33
34 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Physical solution Critical aspect: what happens at next time step! Upwind material mixes uniformly in neighboring grid box. Mathematical solution 34
35 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Physical solution Next time step: some of the material may appear here, which disagrees with physical solution (i.e., non-physical!) Upwind material mixes uniformly in neighboring grid box. Mathematical solution 35
36 3. Numerical Diffusion Consider the continuity equation: Chemical = v Chemical + Production Loss Let s examine the advective term & consider transport of a square wave to the right with velocity v: v t = 0 Move the material 0.1 grid box to the right: Upwind material mixes uniformly in neighboring grid box. NUMERICAL DIFFUSION Physical solution Next time step: some of the material will appear here, which disagrees with physical solution (i.e., is non-physical) Mathematical solution 36
37 37 Prather, M. J., Numerical advection by conservation of second-order moments, JGR, 91, , 1986 developed an algorithm that almost exactly accounts for the advective transport of a constituent in grid box of appropriate size of 3D climate model Advantages: Non-diffusive, mass-conserving, numerically stable slopes slopes first order moments (overshoot and undershoot)! 1 start with a step function upwind (very diffusive) second order moments with limiter (no undershoot) f 0 x second order moments (small overshoot and undershoot) Disadvantage: Computational requirements (need nine 3D arrays for each constituent) Take home message: if conducting 3D calculations of the distribution of a constituent that involves the solution of a PDE, will often have to consider (and account for) effects of numerical diffusion (there are numerous schemes to consider)
38 Prather, M. J., Numerical advection by conservation of second-order moments, JGR, 91, , 1986 developed an algorithm that almost exactly accounts for the advective transport of a constituent in grid box of appropriate size of 3D climate model Advantages: Non-diffusive, mass-conserving, numerically stable slopes slopes first order moments (overshoot and undershoot)! 1 start with a step function upwind (very diffusive) second order moments with limiter (no undershoot) f 0 x second order moments (small overshoot and undershoot) Disadvantage: Computational requirements (need nine 3D arrays for each constituent) Next Mon: review last two assignments, discuss project logistics, & quiz for those who have opted out of HW#1 Rest of Class: focus on projects 38
t tendency advection convergence twisting baroclinicity
RELATIVE VORTICITY EQUATION Newton s law in a rotating frame in z-coordinate (frictionless): U + U U = 2Ω U Φ α p U + U U 2 + ( U) U = 2Ω U Φ α p Applying to both sides, and noting ω U and using identities
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationLecture 7. Heat Conduction. BENG 221: Mathematical Methods in Bioengineering. References
BENG 221: Mathematical Methods in Bioengineering Lecture 7 Heat Conduction References Haberman APDE, Ch. 1. http://en.wikipedia.org/wiki/heat_equation Lecture 5 BENG 221 M. Intaglietta Heat conduction
More informationPHYS 432 Physics of Fluids: Instabilities
PHYS 432 Physics of Fluids: Instabilities 1. Internal gravity waves Background state being perturbed: A stratified fluid in hydrostatic balance. It can be constant density like the ocean or compressible
More informationChapter 3. Finite Difference Methods for Hyperbolic Equations Introduction Linear convection 1-D wave equation
Chapter 3. Finite Difference Methods for Hyperbolic Equations 3.1. Introduction Most hyperbolic problems involve the transport of fluid properties. In the equations of motion, the term describing the transport
More informationChapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction
Chapter 2. Quasi-Geostrophic Theory: Formulation (review) 2.1 Introduction For most of the course we will be concerned with instabilities that an be analyzed by the quasi-geostrophic equations. These are
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationPhysics 8 Monday, October 28, 2013
Physics 8 Monday, October 28, 2013 Turn in HW8 today. I ll make them less difficult in the future! Rotation is a hard topic. And these were hard problems. HW9 (due Friday) is 7 conceptual + 8 calculation
More informationIntroduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes Heat equation The simple parabolic PDE with the initial values u t = K 2 u 2 x u(0, x) = u 0 (x) and some boundary conditions
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationNumerical Hydraulics
ETHZ, Fall 017 Numerical Hydraulics Assignment 4 Numerical solution of 1D solute transport using Matlab http://www.bafg.de/ http://warholian.com Numerical Hydraulics Assignment 4 ETH 017 1 Introduction
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationPDE and Boundary-Value Problems Winter Term 2014/2015
PDE and Boundary-Value Problems Winter Term 2014/2015 Lecture 13 Saarland University 5. Januar 2015 c Daria Apushkinskaya (UdS) PDE and BVP lecture 13 5. Januar 2015 1 / 35 Purpose of Lesson To interpretate
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More information2.2. Methods for Obtaining FD Expressions. There are several methods, and we will look at a few:
.. Methods for Obtaining FD Expressions There are several methods, and we will look at a few: ) Taylor series expansion the most common, but purely mathematical. ) Polynomial fitting or interpolation the
More informationEAS372 Open Book Final Exam 11 April, 2013
EAS372 Open Book Final Exam 11 April, 2013 Professor: J.D. Wilson Time available: 2 hours Value: 30% Please check the Terminology, Equations and Data section before beginning your responses. Answer all
More information= vorticity dilution + tilting horizontal vortices + microscopic solenoid
4.4 Vorticity Eq 4.4.1 Cartesian Coordinates Because ζ = ˆk V, gives D(ζ + f) x minus [v momentum eq. in Cartesian Coordinates] y [u momentum eq. in Cartesian Coordinates] = vorticity dilution + tilting
More informationFDM for parabolic equations
FDM for parabolic equations Consider the heat equation where Well-posed problem Existence & Uniqueness Mass & Energy decreasing FDM for parabolic equations CNFD Crank-Nicolson + 2 nd order finite difference
More informationTutorial 2. Introduction to numerical schemes
236861 Numerical Geometry of Images Tutorial 2 Introduction to numerical schemes c 2012 Classifying PDEs Looking at the PDE Au xx + 2Bu xy + Cu yy + Du x + Eu y + Fu +.. = 0, and its discriminant, B 2
More informationFundamentals Physics
Fundamentals Physics And Differential Equations 1 Dynamics Dynamics of a material point Ideal case, but often sufficient Dynamics of a solid Including rotation, torques 2 Position, Velocity, Acceleration
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationPhysics 106a, Caltech 13 November, Lecture 13: Action, Hamilton-Jacobi Theory. Action-Angle Variables
Physics 06a, Caltech 3 November, 08 Lecture 3: Action, Hamilton-Jacobi Theory Starred sections are advanced topics for interest and future reference. The unstarred material will not be tested on the final
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Implicit Schemes for the Model Problem The Crank-Nicolson scheme and θ-scheme
More informationQuasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States
Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States See CalTech 2005 paper on course web site Free troposphere assumed to have moist adiabatic lapse rate (s* does
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationThe Shallow Water Equations
If you have not already done so, you are strongly encouraged to read the companion file on the non-divergent barotropic vorticity equation, before proceeding to this shallow water case. We do not repeat
More informationAnswers to Problem Set Number 04 for MIT (Spring 2008)
Answers to Problem Set Number 04 for 18.311 MIT (Spring 008) Rodolfo R. Rosales (MIT, Math. Dept., room -337, Cambridge, MA 0139). March 17, 008. Course TA: Timothy Nguyen, MIT, Dept. of Mathematics, Cambridge,
More informationMA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20
MA 201: Partial Differential Equations D Alembert s Solution Lecture - 7 MA 201 (2016), PDE 1 / 20 MA 201 (2016), PDE 2 / 20 Vibrating string and the wave equation Consider a stretched string of length
More informationModel equations for planetary and synoptic scale atmospheric motions associated with different background stratification
Model equations for planetary and synoptic scale atmospheric motions associated with different background stratification Stamen Dolaptchiev & Rupert Klein Potsdam Institute for Climate Impact Research
More informationME Computational Fluid Mechanics Lecture 5
ME - 733 Computational Fluid Mechanics Lecture 5 Dr./ Ahmed Nagib Elmekawy Dec. 20, 2018 Elliptic PDEs: Finite Difference Formulation Using central difference formulation, the so called five-point formula
More informationParameterizing large-scale dynamics using the weak temperature gradient approximation
Parameterizing large-scale dynamics using the weak temperature gradient approximation Adam Sobel Columbia University NCAR IMAGe Workshop, Nov. 3 2005 In the tropics, our picture of the dynamics should
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
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 informationBurgers equation - a first look at fluid mechanics and non-linear partial differential equations
Burgers equation - a first look at fluid mechanics and non-linear partial differential equations In this assignment you will solve Burgers equation, which is useo model for example gas dynamics anraffic
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Ordinary Differential Equations Week 12, Day 1 1 Differential Equations are central to Atmospheric and Ocean Sciences They provide quantitative
More informationImplicit Scheme for the Heat Equation
Implicit Scheme for the Heat Equation Implicit scheme for the one-dimensional heat equation Once again we consider the one-dimensional heat equation where we seek a u(x, t) satisfying u t = νu xx + f(x,
More informationA Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations
A Space-Time Expansion Discontinuous Galerkin Scheme with Local Time-Stepping for the Ideal and Viscous MHD Equations Ch. Altmann, G. Gassner, F. Lörcher, C.-D. Munz Numerical Flow Models for Controlled
More informationAdvection / Hyperbolic PDEs. PHY 604: Computational Methods in Physics and Astrophysics II
Advection / Hyperbolic PDEs Notes In addition to the slides and code examples, my notes on PDEs with the finite-volume method are up online: https://github.com/open-astrophysics-bookshelf/numerical_exercises
More informationFUNDAMENTALS OF FINITE DIFFERENCE METHODS
FUNDAMENTALS OF FINITE DIFFERENCE METHODS By Deep Gupta 3 rd Year undergraduate, Mechanical Engg. Deptt., IIT Bombay Supervised by: Prof. Gautam Biswas, IIT Kanpur Acknowledgements It has been a pleasure
More informationLast time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler:
Lecture 7 18.086 Last time: Diffusion - Numerical scheme (FD) Heat equation is dissipative, so why not try Forward Euler: U j,n+1 t U j,n = U j+1,n 2U j,n + U j 1,n x 2 Expected accuracy: O(Δt) in time,
More informationContents of lecture 2b. Lectures 2a & 2b. Physical vs. computational coordinates [2] Physical vs. computational coordinates [1]
Contents of lecture b Lectures a & b P. A. Zegeling Mathematical Institute Utrecht University The Netherlands Parameter-free non-singular moving grids in D: Theory & properties Application to resistive
More informationBy convention, C > 0 for counterclockwise flow, hence the contour must be counterclockwise.
Chapter 4 4.1 The Circulation Theorem Circulation is a measure of rotation. It is calculated for a closed contour by taking the line integral of the velocity component tangent to the contour evaluated
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
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 informationEAS372 Open Book Final Exam 11 April, 2013
EAS372 Open Book Final Exam 11 April, 2013 Professor: J.D. Wilson Time available: 2 hours Value: 30% Please check the Terminology, Equations and Data section before beginning your responses. Answer all
More informationOcean General Circulation Modeling
Lecture 8 Ocean General Circulation Modeling 9. The equations of motion: Navier-Stokes The governing equations for a real fluid are the Navier-Stokes equations (conservation of linear momentum and mass
More informationLecture 3: Hamilton-Jacobi-Bellman Equations. Distributional Macroeconomics. Benjamin Moll. Part II of ECON Harvard University, Spring
Lecture 3: Hamilton-Jacobi-Bellman Equations Distributional Macroeconomics Part II of ECON 2149 Benjamin Moll Harvard University, Spring 2018 1 Outline 1. Hamilton-Jacobi-Bellman equations in deterministic
More informationToday s Lecture (Lecture 5): General circulation of the atmosphere
Climate Dynamics (Summer Semester 2017) J. Mülmenstädt Today s Lecture (Lecture 5): General circulation of the atmosphere Reference Hartmann, Global Physical Climatology (1994), Ch. 2, 3, 6 Peixoto and
More informationLecture 3 Numerical Solutions to the Transport Equation
Lecture 3 Numerical Solutions to the Transport Equation Introduction I There are numerous methods for solving the transport problem numerically. First we must recognize that we need to solve two problems:
More informationAdvection in two dimensions
Lecture 0 Advection in two dimensions 6. Stability of multiple terms (in multiple dimensions) When we analyzed the stability of time-stepping methods we tended to consider either a single damping term
More informationTransformed Eulerian Mean
Chapter 15 Transformed Eulerian Mean In the last few lectures we introduced some fundamental ideas on 1) the properties of turbulent flows in rotating stratified environments, like the ocean and the atmosphere,
More informationEffective Depth of Ekman Layer.
5.5: Ekman Pumping Effective Depth of Ekman Layer. 2 Effective Depth of Ekman Layer. Defining γ = f/2k, we derived the solution u = u g (1 e γz cos γz) v = u g e γz sin γz corresponding to the Ekman spiral.
More informationPhysical Oceanography, MSCI 3001 Oceanographic Processes, MSCI Dr. Katrin Meissner Week 5.
Physical Oceanography, MSCI 3001 Oceanographic Processes, MSCI 5004 Dr. Katrin Meissner k.meissner@unsw.e.au Week 5 Ocean Dynamics Transport of Volume, Heat & Salt Flux: Amount of heat, salt or volume
More informationChapter 10: Partial Differential Equations
1.1: Introduction Chapter 1: Partial Differential Equations Definition: A differential equations whose dependent variable varies with respect to more than one independent variable is called a partial differential
More informationLecture Notes on Numerical Schemes for Flow and Transport Problems
Lecture Notes on Numerical Schemes for Flow and Transport Problems by Sri Redeki Pudaprasetya sr pudap@math.itb.ac.id Department of Mathematics Faculty of Mathematics and Natural Sciences Bandung Institute
More informationIterative Methods for Ax=b
1 FUNDAMENTALS 1 Iterative Methods for Ax=b 1 Fundamentals consider the solution of the set of simultaneous equations Ax = b where A is a square matrix, n n and b is a right hand vector. We write the iterative
More informationGravity Waves. Lecture 5: Waves in Atmosphere. Waves in the Atmosphere and Oceans. Internal Gravity (Buoyancy) Waves 2/9/2017
Lecture 5: Waves in Atmosphere Perturbation Method Properties of Wave Shallow Water Model Gravity Waves Rossby Waves Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature
More informationApproximations of diffusions. Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine
Lecture 3b Approximations of diffusions Mathématiques appliquées (MATH0504-1) B. Dewals, Ch. Geuzaine V1.1 04/10/2018 1 Learning objectives Become aware of the existence of stability conditions for the
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationThe Advanced Research WRF (ARW) Dynamics Solver
Dynamics: Introduction The Advanced Research WRF (ARW) Dynamics Solver 1. What is a dynamics solver? 2. Variables and coordinates 3. Equations 4. Time integration scheme 5. Grid staggering 6. Advection
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationBaroclinic wave. Atmospheric Dynamics: lecture 14 18/12/15. Topics. Chapter 9: Baroclinic waves and cyclogenesis. What is a baroclinic wave?
Atmospheric Dynamics: lecture 14 (http://www.staff.science.uu.nl/~delde102/) Topics Chapter 9: Baroclinic waves and cyclogenesis What is a baroclinic wave? Quasi-geostrophic equations Omega equation Original
More informationIntroduction to Partial Differential Equations
Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat
More informationDiscretization of Convection Diffusion type equation
Discretization of Convection Diffusion type equation 10 th Indo German Winter Academy 2011 By, Rajesh Sridhar, Indian Institute of Technology Madras Guides: Prof. Vivek V. Buwa Prof. Suman Chakraborty
More informationNeeds work : define boundary conditions and fluxes before, change slides Useful definitions and conservation equations
Needs work : define boundary conditions and fluxes before, change slides 1-2-3 Useful definitions and conservation equations Turbulent Kinetic energy The fluxes are crucial to define our boundary conditions,
More informationFinite Difference Method for PDE. Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36
Finite Difference Method for PDE Y V S S Sanyasiraju Professor, Department of Mathematics IIT Madras, Chennai 36 1 Classification of the Partial Differential Equations Consider a scalar second order partial
More informationFinite Difference Methods for
CE 601: Numerical Methods Lecture 33 Finite Difference Methods for PDEs Course Coordinator: Course Coordinator: Dr. Suresh A. Kartha, Associate Professor, Department of Civil Engineering, IIT Guwahati.
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationGoverning Equations and Scaling in the Tropics
Governing Equations and Scaling in the Tropics M 1 ( ) e R ε er Tropical v Midlatitude Meteorology Why is the general circulation and synoptic weather systems in the tropics different to the those in the
More information2.29 Numerical Fluid Mechanics Spring 2015 Lecture 13
REVIEW Lecture 12: Spring 2015 Lecture 13 Grid-Refinement and Error estimation Estimation of the order of convergence and of the discretization error Richardson s extrapolation and Iterative improvements
More informationQ ( q(m, t 0 ) n) S t.
THE HEAT EQUATION The main equations that we will be dealing with are the heat equation, the wave equation, and the potential equation. We use simple physical principles to show how these equations are
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 information1/3/2011. This course discusses the physical laws that govern atmosphere/ocean motions.
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More information22.615, MHD Theory of Fusion Systems Prof. Freidberg Lecture 14: Formulation of the Stability Problem
.65, MHD Theory of Fusion Systems Prof. Freidberg Lecture 4: Formulation of the Stability Problem Hierarchy of Formulations of the MHD Stability Problem for Arbitrary 3-D Systems. Linearized equations
More informationA stable treatment of conservative thermodynamic variables for semi-implicit semi-lagrangian dynamical cores
A stable treatment of conservative thermodynamic variables for semi-implicit semi-lagrangian dynamical cores Kevin Viner Naval Research Laboratory, Monterey, CA September 26, 2012 Kevin Viner (NRL) PDE
More informationSAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany
WAVES IN THE OCEANS Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany Keywords: Wind waves, dispersion, internal waves, inertial oscillations, inertial waves,
More information3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1
Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University A Model Problem and Its Difference Approximations 1-D Initial Boundary Value
More informationBasics on Numerical Methods for Hyperbolic Equations
Basics on Numerical Methods for Hyperbolic Equations Professor Dr. E F Toro Laboratory of Applied Mathematics University of Trento, Italy eleuterio.toro@unitn.it http://www.ing.unitn.it/toro October 8,
More informationWave mo(on in the tropics. Shayne McGregor Lecturer (Climate), EAE, Monash Associate Inves(gator (ARC CSS)
Wave mo(on in the tropics Shayne McGregor Lecturer (Climate), EAE, Monash Associate Inves(gator (ARC CSS) Equatorially trapped waves The linear SWM Equatorial Kelvin waves Other equatorially trapped waves
More information2.3. Quantitative Properties of Finite Difference Schemes. Reading: Tannehill et al. Sections and
3 Quantitative Properties of Finite Difference Schemes 31 Consistency, Convergence and Stability of FD schemes Reading: Tannehill et al Sections 333 and 334 Three important properties of FD schemes: Consistency
More information+ ω = 0, (1) (b) In geometric height coordinates in the rotating frame of the Earth, momentum balance for an inviscid fluid is given by
Problem Sheet 1: Due Thurs 3rd Feb 1. Primitive equations in different coordinate systems (a) Using Lagrangian considerations and starting from an infinitesimal mass element in cartesian coordinates (x,y,z)
More informationLecture 1: Introduction to Linear and Non-Linear Waves
Lecture 1: Introduction to Linear and Non-Linear Waves Lecturer: Harvey Segur. Write-up: Michael Bates June 15, 2009 1 Introduction to Water Waves 1.1 Motivation and Basic Properties There are many types
More informationSpotlight on Laplace s Equation
16 Spotlight on Laplace s Equation Reference: Sections 1.1,1.2, and 1.5. Laplace s equation is the undriven, linear, second-order PDE 2 u = (1) We defined diffusivity on page 587. where 2 is the Laplacian
More informationCache Oblivious Stencil Computations
Cache Oblivious Stencil Computations S. HUNOLD J. L. TRÄFF F. VERSACI Lectures on High Performance Computing 13 April 2015 F. Versaci (TU Wien) Cache Oblivious Stencil Computations 13 April 2015 1 / 19
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Introduction to Hyperbolic Equations The Hyperbolic Equations n-d 1st Order Linear
More informationAtmospheric Dynamics: lecture 11
Atmospheric Dynamics: lecture 11 (http://www.staff.science.uu.nl/~delde102/) Chapter 9 Baroclinic waves and cyclogenesis What is a baroclinic wave? Quasi-geostrophic equations Omega equation Original articles:
More informationOcean dynamics: the wind-driven circulation
Ocean dynamics: the wind-driven circulation Weston Anderson March 13, 2017 Contents 1 Introduction 1 2 The wind driven circulation (Ekman Transport) 3 3 Sverdrup flow 5 4 Western boundary currents (western
More information7 Hyperbolic Differential Equations
Numerical Analysis of Differential Equations 243 7 Hyperbolic Differential Equations While parabolic equations model diffusion processes, hyperbolic equations model wave propagation and transport phenomena.
More informationMA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation
MA 201, Mathematics III, July-November 2016, Partial Differential Equations: 1D wave equation (contd.) and 1D heat conduction equation Lecture 12 Lecture 12 MA 201, PDE (2016) 1 / 24 Formal Solution of
More information2 Transport of heat, momentum and potential vorticity
Transport of heat, momentum and potential vorticity. Conventional mean momentum equation We ll write the inviscid equation of onal motion (we ll here be using log-pressure coordinates with = H ln p=p,
More informationMechanics Physics 151
Mechanics Physics 151 Lecture 4 Continuous Systems and Fields (Chapter 13) What We Did Last Time Built Lagrangian formalism for continuous system Lagrangian L Lagrange s equation = L dxdydz Derived simple
More information7 The General Circulation
7 The General Circulation 7.1 The axisymmetric state At the beginning of the class, we discussed the nonlinear, inviscid, axisymmetric theory of the meridional structure of the atmosphere. The important
More informationIn this section of notes, we look at the calculation of forces and torques for a manipulator in two settings:
Introduction Up to this point we have considered only the kinematics of a manipulator. That is, only the specification of motion without regard to the forces and torques required to cause motion In this
More informationSIO 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes
NAME: SIO 210 Introduction to Physical Oceanography Mid-term examination November 3, 2014; 1 hour 20 minutes Closed book; one sheet of your own notes is allowed. A calculator is allowed. (100 total points.)
More informationSolution of Differential Equation by Finite Difference Method
NUMERICAL ANALYSIS University of Babylon/ College of Engineering/ Mechanical Engineering Dep. Lecturer : Dr. Rafel Hekmat Class : 3 rd B.Sc Solution of Differential Equation by Finite Difference Method
More informationThe Euler Equation of Gas-Dynamics
The Euler Equation of Gas-Dynamics A. Mignone October 24, 217 In this lecture we study some properties of the Euler equations of gasdynamics, + (u) = ( ) u + u u + p = a p + u p + γp u = where, p and u
More informationA semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate
A semi-implicit non-hydrostatic covariant dynamical kernel using spectral representation in the horizontal and a height based vertical coordinate Juan Simarro and Mariano Hortal AEMET Agencia Estatal de
More information