Exact Solution of the Linear, 4D-Discrete, Implicit, Semi-Lagrangian Dynamic Equations with Orographic Forcing
|
|
- Anastasia Fowler
- 6 years ago
- Views:
Transcription
1 1 Exact Solution of the Linear, 4D-Discrete, Implicit, Semi-Lagrangian Dynamic Equations with Orographic Forcing Marko Zirk, Rein Rõõm Tartu University, Estonia
2 1 1 Introduction Exact solution of completely discrete linear equations has multiple outputs: Discovery of the exact structure of the solution of a non-stationary numerical model Investigation of specific details of flows for optional (non-homogenous) wind and temperature stratification Investigation of the influence of vertical discretization (variable z) on the quality of numerical models Investigation of the time step t size impact on the numerical stability and precision of numerical schemes Testing of adiabatic cores of NWP models
3 2 2 Domain of integration Geometry is plane: p s (x, y) p = const. Orography (taken into account in boundary conditions but not in geometry) h is presented by 2D Agnesi profile h(x, y) = 1 + (x/a (h height, a x, a y -half-widths): x ) 2 + (y/a y ) 2. The reference state of the atmosphere is given by wind-speed U(p) and temperature T(p). The Lagrangian time-derivative is approximated as d dt = t + U(p) x.
4 3 3 Continuous linear equations Linear form of White equations (White 1989, Rm 21). Notation is standard, except of ϕ which is the nonhydrostatic geopotential fluctuation. dω dt = R p H 2 T p2 H 2 ϕ p, du dt = ϕ f u, dt dt = θω p, (1) u + ω p =, dp s dt + da dt = ω p. Reference state parameters: H(p) = RT(p), p θ(p) = R c p T(p) p T p. Vertical boundary conditions ϕ p = gh(p ) p s p. ω p= =.
5 4 4 Vertical grid and difference operators Vertical grid is the Arakawa C-grid. Fields ϕ k, u k, pressure difference p k and non-dimensional height difference ζ k are located on full levels. Pressure p k+1/2, T k+1/2, H k+1/2, θ k+1/2, ω k+1/2 and T k+1/2 belong to the half-levels. Vertical discretization is originated by p k+1/2. Non-dimensional height ζ k+1/2 : p k+1/2 = p klev+1/2 e ζ k+1/2, ζ k+1/2 = ln(p klev+1/2 /p k+1/2 ). Full-level height and pressure: height differences ζ k = ζ k 1/2 + ζ k+1/2 2, p k = p klev+1/2 e ζ k. ζ k = ζ k 1/2 ζ k+1/2, ζ k+1/2 = ζ k ζ k+1.
6 5 The vertical difference derivatives are ( ) ψ = ψ ( ) k ψ, p p k ζ k p k k+1/2 = ψ k+1/2 p k+1/2 ζ k+1/2. 5 Horizontal grid and difference operators ψ n ijk = ψ(x i, y j, p k, t n ), ψ n 1 i l k,jk = ψ(x i U k t, y j, p k, t n 1 ), l k = U k t x, non-dimensional displacement along x-axes during t. δ t ψijk n = ψn ij,k ψn 1 i l k,j,k t δ i ψ = 1 x δ j ψ = 1 y, ψ n ijk = ψn ijk + ψn 1 i l k,jk 2 ( ) ψi+1/2 ψ i 1/2, δi+1/2 ψ = 1 x (ψ i+1 ψ i ), ( ) ψj+1/2 ψ j 1/2, δj+1/2 ψ = 1 y (ψ j+1 ψ j )
7 6 ψ x i = 1 2 ψ y j = 1 2 ( ψi 1/2 + ψ i+1/2 ), ψ x i+1/2 = 1 2 (ψ i + ψ i+1 ), ( ψj 1/2 + ψ j+1/2 ), ψ y j+1/2 = 1 2 (ψ j + ψ j+1 ) 6 Semi-Implicit, Semi-Lagrangian Equations Equations (1) become in SISL presentation ( p δ t ωijk+1/2 n = R T H )k+1/2 2 n ijk+1/2 ( p H 2 ζ ) k+1/2 ϕ n ijk+1/2, δ t u n i+1/2,jk = (δ xϕ) n i+1/2,jk + fvxyn i+1/2,jk, δ t vij+1/2,k n = (δ yϕ) n ij+1/2,k fuxyn ij+1/2,k, ( ) δ t T n θ ijk+1/2 = ω n ijk+1/2 p, (2) k+1/2 ( u) n ijk + ωn ijk (p ζ) k =,
8 7 δ t (p s) n ij,klev+1/2 + δ t(a) ij,klev+1/2 = ω n ijklev+1/2. 7 Embedding into 3D Fourier space In following, the discrete Fourier transformation is employed in the form f i = 2N 1 q= ˆf q e iη qi, kus ˆf q = 1 2N 2N 1 i= f q e iη qi ja η q = π q N. Orography in the 2D Fourier basis is a ij a(x i, y j ) = qr â qr e i(ηx q i+ηy r j), Using short notation A general presentation holds ψ n ijk = { } T ijk+1/2 n, ωn ijk+1/2, un i+1/2jk, vn ij+1/2k, ϕn ijk, (p s) n ij ψ n ijk = qrs ˆψ s qrke i(ηx q i+ηy r j ds qrk n), (3)
9 8 where d s qrk are the non-dimensional frequencies. 8 Normal mode equations: (3) (2) : iν s qrk+1/2ˆωs qr,k+1/2 = R ( p H 2 )k+1/2 ( ) p ˆT qrk+1/2 s H 2 ζ k+1/2 ˆϕ s qrk+1/2, iν s qrkû s qrk = iµ x q ˆϕ s qrk +F qrˆv s qrk, iν s qrkˆv s qrk = iµ y r ˆϕ s qrk F qrû s qrk, (4) iν s qrk+1/2 ˆT s qrk+1/2 = θ k+1/2 p k+1/2 ˆω s qrk+1/2, i[(µx q) û s qrk+(µ y r) ˆv s qrk]+ ˆωs qrk (p ζ) k =. F qr = fa x q(a y r), a x q = 1 + eiηx q 2, a y r = 1 + eiηy r 2 ν s qrk = 2 t tan[1 2 (ηx q l k d s qrk)]., µ x q = eiηx q 1 i x, µy r = eiη y r 1. i y Stationary case s = : d qrk = ν s qrk = ν qk = 2 t tan(1 2 ηx q l k ). Non-stationary case s : d s qrk = η x q l k 2 d s qr ν s qrk = ν s qr = 2 t tan(ds qr).
10 9 9 Wave equation for vertical velocity ˆω s qrk+1/2 From normal mode equations, the wave equation follows for ˆω s qrk+1/2 (Lˆω s ) qrk+1/2 + λ s qrk+1/2ˆωs qrk+1/2 =, (5) where λ s qrk+1/2 = H2 k+1/2 µ2 qr (ν s qrk+1/2 )2 N 2 k+1/2 F qr 2 (ν s qrk+1/2 )2. (Lˆω s ) qrk+1/2 = L + qrk+1/2 ˆωs qrk+1 L qrk+1/2 ˆωs qrk, L + qrk+1/2 = νs qrk+1 ν s qrk+1/2 e ζ k+1/2 ζ k+1/2 ζ k+1, L qrk+1/2 = νs qrk ν s qrk+1/2 e ζ k/2 ζ k+1/2 ζ k, ν s qrk = (νs qrk )2 F qr 2 ν s qrk.
11 1 1 Stationary solution 1.1 Main recurrence formula Solutions of Eq. (5) are sought in the form k χ k+1/2 = χ 1/2 c j, χ 1/2 = 1. (8) j= From (5), a nonlinear, two-point recurrence follows for growth factors c k L + k+1/2 (c k+1 1) + L k+1/2 (1/c k 1) + λ k+1/2 =. (9) The initial value c(), required to start this simple straightforward formula, can be established, assuming the homogeneous layer in the top in the free buoyancy-wave radiating state. The solution then presents ω qrk+1/2 = Cχ k+1/2, where C is specified from boundary condition on the surface.
12 Special case of homogeneous stratification Homogeneous atmosphere : ζ k = ζ k+1/2 = ζ, Homogeneous atmosphere : λ qrk = λ qr = H 2 µ 2 qr Consequently c k = c, χ k+1/2 = χ 1/2 c k, L ± k+1/2 = e ζ/2 ( ζ) 2, (ν qr) 2 N 2 F qr 2 (ν qr) 2. where c satisfies equation e ζ/2 (c 1) + e ζ/2 (c 1 1) + ( ζ) 2 λ =. The solution of this equation is: ( trapped wave, Q 1 : c = e ζ/2 Q + ) Q 2 1, where free wave, Q < 1 : c = e ( ζ/2 i ξ), Q Q qr = cosh( ζ/2) 1 2 ( ζ)2 λ qr. 1 Q 2 ξ = Arctan, Q
13 12 11 Non-stationary free waves Assuming ν s qrk+1/2 = νs qr, the Eq. (5) simplifies, as the coefficients L ± become independent of the frequency (and hence, of horizontal wave-numbers q, r): L + k+1/2 = p k+1/2/p k+1 ζ k+1/2 ζ k+1 = e ζk+1/2 ζ k+1/2 ζ k+1, L k+1/2 = p k+1/2/p k ζ k+1/2 ζ k = e ζ k/2 ζ k+1/2 ζ k. Quest quantities are the frequencies ν s qr in λ and the corresponding eigenfunctions. The solution presents as a combination of functions (8), found from recurrence (9). The wanted frequencies are selected by the homogeneous bottom boundary condition.
14 13 12 Examples: homogeneous stratification 3 V z : u = 5 m/s, T = 265 K, h =.2 km, a x = 2 km, zlev = 3, xlev = 124, z =.1 km, x =.4 km X, km V z : u = 1 m/s, T = 265 K, h =.2 km, a x = 2 km, zlev = 3, xlev = 124, z =.1 km, x =.4 km X, km
15 14 13 Examples: BW refraction on tropopause 3 3 Vz: D(Vz)=.1m/s,U= 5m/s,gam=4.7K/km,h=1m,a= 3km,dx=.55km,Mlev=3,dz=1m T U γ=4.7 K/km T, K 5 1 U, m/s X, km 3 3 Vz: D(Vz)=.1m/s,U= 5m/s,gam=9.5K/km,h=1m,a= 3km,dx=.5km,Mlev=3,dz=1m T U γ=9.5 K/km T, K U, m/s X, km
16 Vz: D(Vz)=.1m/s,U= 5-2m/s,gam=9.5K/km,h=1m,a= 3km,dx=.5km,Mlev=3,dz=1m T U γ=9.5 K/km T, K U, m/s X, km 3 3 Vz: D(Vz)=.1m/s,U= 1-5m/s,gam=9.5K/km,h=1m,a= 3km,dx=.5km,Mlev=3,dz=1m T U γ=8.8 K/km T, K U, m/s X, km
17 16 14 Examples: Hyperbolic reference wind profile Reference temperature and wind ω: u = m/s, γ =.65 K/m, a x = 2 km, h =.2 km, xlev = 496, zlev = 3, x =.4 km, z = 1 m 2 T U γ=6.5k/km p, hpa p,hpa T, K 2 U, m/s x, km ω: u = m/s, γ =.65 K/m, a x = 2 km, h =.2 km, xlev = 496, zlev = 3, x =.4 km, z = 1 m 2 p, hpa x, km
18 17 15 Example: Testing NH HIRLAM HIRLAM, V z : D(V z )=.5m/s,a x =3km, h=1m, MLEV=1,dx=.55km,dt=3s,6 steps Reference temperature and wind T γ=8. K/km U p, hpa γ=4.5 K/km 8 1 p,hpa T, K 2 U, m/s X, km V z : D(V z )=.5m/s; U,T-HIRLAM,h=1m,a x = 3km,dx=.55km,MLEV=2,dz=1m Reference temperature and wind T γ=8. K/km U p, hpa γ=4.5 K/km 8 1 p,hpa T, K 2 U, m/s X, km
19 18 16 References White, A.A.,1989: An extended version of nonhydrostatic, pressure coordinate model. QJRMS115, Robert, A., T. L. Yee, H. Richie, 1985: A semi-lagrangian and semiimplicit integration scheme for multi-level atmospheric models. Mon. Weather Rev., 113, Rõõm, R., 21: Nonhydrostatic adiabatic kernel for HIRLAM. Part I: Fundametals of nonhydrostatic dynamics in pressure-related coordinates. HIRLAM Technical Report, 48, 26p.
20 19 Vz: D(Vz)=.1m/s,U,T-HIRLAM,h=1m,a= 3km,dx=.55km,Mlev=22,dz=1m X, km
21 2 Linear, Vz: D(Vz)=.1m/s,U,T-HIRLAM,h=1m,a= 3km,dx=.55km,Mlev=22,dz=1m X, km
22 21 HIRLAM, V z : D(V z )=.1m/s,a x =3km, h=1m, MLEV=1,dx=.55km,dt=3s,6 steps X, km
Linear model for investigation of nonlinear NWP model accuracy. Marko Zirk, University of Tartu
Linear model for investigation of nonlinear NWP model accuracy Marko Zirk, University of Tartu Introduction A method for finding numerical solution of non-hydrostatic linear equations of atmospheric dynamics
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 informationOverview of the Numerics of the ECMWF. Atmospheric Forecast Model
Overview of the Numerics of the Atmospheric Forecast Model M. Hortal Seminar 6 Sept 2004 Slide 1 Characteristics of the model Hydrostatic shallow-atmosphere approimation Pressure-based hybrid vertical
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 informationModeling Interactions of Soliton Trains. Effects of External Potentials. Part II
Modeling Interactions of Soliton Trains. Effects of External Potentials. Part II Michail Todorov Department of Applied Mathematics and Computer Science Technical University of Sofia, Bulgaria Work done
More informationThe Finite Element Method
The Finite Element Method 3D Problems Heat Transfer and Elasticity Read: Chapter 14 CONTENTS Finite element models of 3-D Heat Transfer Finite element model of 3-D Elasticity Typical 3-D Finite Elements
More informationA Uniform PV Framework for Balanced Dynamics
A Uniform PV Framework for Balanced Dynamics vertical structure of the troposphere surface quasigeostrophic models Dave Muraki Simon Fraser University Greg Hakim University of Washington Chris Snyder NCAR
More informationSemi-implicit methods, nonlinear balance, and regularized equations
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 1 6 (7 Published online 9 January 7 in Wiley InterScience (www.interscience.wiley.com.1 Semi-implicit methods, nonlinear balance, and regularized equations
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationFigure 1. adiabatically. The change in pressure experienced by the parcel is. dp = -ρ o gξ
6. Internal waves Consider a continuously stratified fluid with ρ o (z) the vertical density profile. z p' ξ p ρ ρ ο (z) Figure 1. Figure by MIT OpenCourseWare. At a point P raise a parcel of water by
More informationDefn: A linear recurrence relation has constant coefficients if the a i s are constant. Ch 14: 1, 4, 5, 10, 13, 18, 22, 24, 25 and
7.4: linear homogeneous recurrence relation: Defn: A recurrence relation is linear if h n = a 1 (n)h n 1 + a (n)h n +... + a k (n)h n k + b(n) A recurrence relation has order k if a k 0 Ex: Derangement
More informationIntegrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows
Integrodifferential Hyperbolic Equations and its Application for 2-D Rotational Fluid Flows Alexander Chesnokov Lavrentyev Institute of Hydrodynamics Novosibirsk, Russia chesnokov@hydro.nsc.ru July 14,
More informationExact Solutions of the Einstein Equations
Notes from phz 6607, Special and General Relativity University of Florida, Fall 2004, Detweiler Exact Solutions of the Einstein Equations These notes are not a substitute in any manner for class lectures.
More informationQuasi-geostrophic system
Quasi-eostrophic system (or, why we love elliptic equations for QGPV) Charney s QG the motion of lare-scale atmospheric disturbances is overned by Laws of conservation of potential temperature and potential
More informationRome - May 12th Université Paris-Diderot - Laboratoire Jacques-Louis Lions. Mean field games equations with quadratic
Université Paris-Diderot - Laboratoire Jacques-Louis Lions Rome - May 12th 2011 Hamiltonian MFG Hamiltonian on the domain [0, T ] Ω, Ω standing for (0, 1) d : (HJB) (K) t u + σ2 2 u + 1 2 u 2 = f (x, m)
More informationPAPER 71 COSMOLOGY. Attempt THREE questions There are seven questions in total The questions carry equal weight
MATHEMATICAL TRIPOS Part III Friday 31 May 00 9 to 1 PAPER 71 COSMOLOGY Attempt THREE questions There are seven questions in total The questions carry equal weight You may make free use of the information
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 informationOn explicit integration of two non-holonomic problems
On explicit integration of two non-holonomic problems Alexey V. Borisov 1 1 Institute of Computer Sciences, Izhevsk, Russia Generalized Chaplygin systems Equations of motion of the generalized Chaplygin
More informationA Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations
A Nodal High-Order Discontinuous Galerkin Dynamical Core for Climate Simulations Institute for Mathematics Applied to the Geosciences (IMAGe) National Center for Atmospheric Research (NCAR) Boulder CO
More informationACOUSTIC BLACK HOLES. MASSIMILIANO RINALDI Université de Genève
ACOUSTIC BLACK HOLES MASSIMILIANO RINALDI Université de Genève OUTLINE Prelude: GR vs QM Hawking Radiation: a primer Acoustic Black Holes Hawking Radiation in Acoustic Black Holes Acoustic Black Holes
More information(TRAVELLING) 1D WAVES. 1. Transversal & Longitudinal Waves
(TRAVELLING) 1D WAVES 1. Transversal & Longitudinal Waves Objectives After studying this chapter you should be able to: Derive 1D wave equation for transversal and longitudinal Relate propagation speed
More informationle LPTMS en Bretagne... photo extraite du site
le LPTMS en Bretagne... 1 photo extraite du site http://www.chateau-du-val.com le LPTMS en Bretagne... 1 2 Quantum signature of analog Hawking radiation in momentum space Nicolas Pavloff LPTMS, CNRS, Univ.
More informationChapter 5. Fundamentals of Atmospheric Modeling
Overhead Slides for Chapter 5 of Fundamentals of Atmospheric Modeling by Mark Z. Jacobson Department of Civil & Environmental Engineering Stanford University Stanford, CA 94305-4020 January 30, 2002 Altitude
More informationWater wave radiation by a submerged rough disc
radiation by a Leandro Farina Institute of Mathematics & Post-Graduation Program in Geosciences UFRGS Porto Alegre, Brazil BCAM, Bilbao, 14 June 2012. Outline 1 2 3 TLP Mars platform. Projected for suporting
More informationModule-3: Kinematics
Module-3: Kinematics Lecture-20: Material and Spatial Time Derivatives The velocity acceleration and velocity gradient are important quantities of kinematics. Here we discuss the description of these kinematic
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 informationBoundary Layers: Stratified Fluids
Boundary Layers: Stratified Fluids Lecture 3 by Jeroen Hazewinkel continued from lecture 2 Using w I = /(2σS) 2 T I, the interior of the cylinder is be described by 2σS 2 T I + 1/2 4 2 This result can
More informationImpact of Turbulence on the Intensity of Hurricanes in Numerical Models* Richard Rotunno NCAR
Impact of Turbulence on the Intensity of Hurricanes in Numerical Models* Richard Rotunno NCAR *Based on: Bryan, G. H., and R. Rotunno, 2009: The maximum intensity of tropical cyclones in axisymmetric numerical
More informationσ 2 + π = 0 while σ satisfies a cubic equation λf 2, σ 3 +f + β = 0 the second derivatives of the potential are = λ(σ 2 f 2 )δ ij, π i π j
PHY 396 K. Solutions for problem set #4. Problem 1a: The linear sigma model has scalar potential V σ, π = λ 8 σ + π f βσ. S.1 Any local minimum of this potential satisfies and V = λ π V σ = λ σ + π f =
More informationMassive Spinors and ds/cft Correspondence
Massive Spinors and ds/cft Correspondence Farhang Loran arxiv:hep-th/00135v3 16 Jun 00 Department of Physics, Isfahan University of Technology IUT) Isfahan, Iran, Institute for Studies in Theoretical Physics
More information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationWeather Forecasting Models in Met Éireann. Eoin Whelan UCD Seminar 3 rd April 2012
Weather Forecasting Models in Met Éireann Eoin Whelan UCD Seminar 3 rd April 2012 Overview Background HIRLAM Models Local Implementation Verification Development work Background Met Éireann Dept of the
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
More informationPreliminary Examination - Day 2 May 16, 2014
UNL - Department of Physics and Astronomy Preliminary Examination - Day May 6, 04 This test covers the topics of Thermodynamics and Statistical Mechanics (Topic ) and Mechanics (Topic ) Each topic has
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationChange of Measure formula and the Hellinger Distance of two Lévy Processes
Change of Measure formula and the Hellinger Distance of two Lévy Processes Erika Hausenblas University of Salzburg, Austria Change of Measure formula and the Hellinger Distance of two Lévy Processes p.1
More informationEquations of linear stellar oscillations
Chapter 4 Equations of linear stellar oscillations In the present chapter the equations governing small oscillations around a spherical equilibrium state are derived. The general equations were presented
More informationUniversity of Illinois at Chicago Department of Physics
University of Illinois at Chicago Department of Physics Electromagnetism Qualifying Examination January 4, 2017 9.00 am - 12.00 pm Full credit can be achieved from completely correct answers to 4 questions.
More informationEvolution of semiclassical Wigner function (the higher dimensio
Evolution of semiclassical Wigner function (the higher dimensional case) Workshop on Fast Computations in Phase Space, WPI-Vienna, November 2008 Dept. Appl. Math., Univ. Crete & IACM-FORTH 1 2 3 4 5 6
More informationStable Semi-Discrete Schemes for the 2D Incompressible Euler Equations
Stable Semi-Discrete Schemes for the D Incompressible Euler Equations Doron Levy Department of Mathematics Stanford University dlevy@math.stanford.edu http://math.stanford.edu/ dlevy Maryland, May 004
More informationIntroduction to Path Integrals
Introduction to Path Integrals Consider ordinary quantum mechanics of a single particle in one space dimension. Let s work in the coordinate space and study the evolution kernel Ut B, x B ; T A, x A )
More informationPAPER 310 COSMOLOGY. Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight.
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 9:00 am to 12:00 pm PAPER 310 COSMOLOGY Attempt no more than THREE questions. There are FOUR questions in total. The questions carry equal weight. STATIONERY
More informationWHY BLACK HOLES PHYSICS?
WHY BLACK HOLES PHYSICS? Nicolò Petri 13/10/2015 Nicolò Petri 13/10/2015 1 / 13 General motivations I Find a microscopic description of gravity, compatibile with the Standard Model (SM) and whose low-energy
More informationGlobal Mandelbrot transformation for Maslov tunnel pseudodifferential equations
Український математичний вiсник Том 5 2008), 1, 46 58 Global Mandelbrot transformation for Maslov tunnel pseudodifferential equations Vladimir G. Danilov Presented by E. Ya. Khruslov) Abstract. We construct
More informationToday s Lecture: Atmosphere finish primitive equations, mostly thermodynamics
Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics Reference Peixoto and Oort, Sec. 3.1, 3.2, 3.4, 3.5 (but skip the discussion of oceans until next week); Ch. 10 Thermodynamic
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 information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation / 59 Agenda Agenda
More informationChapter 2. The continuous equations
Chapter. The continuous equations Fig. 1.: Schematic of a forecast with slowly varying weather-related variations and superimposed high frequency Lamb waves. Note that even though the forecast of the slow
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationGAUSSIAN PROCESSES AND THE LOCAL TIMES OF SYMMETRIC LÉVY PROCESSES
GAUSSIAN PROCESSES AND THE LOCAL TIMES OF SYMMETRIC LÉVY PROCESSES MICHAEL B. MARCUS AND JAY ROSEN CITY COLLEGE AND COLLEGE OF STATEN ISLAND Abstract. We give a relatively simple proof of the necessary
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationHopping transport in disordered solids
Hopping transport in disordered solids Dominique Spehner Institut Fourier, Grenoble, France Workshop on Quantum Transport, Lexington, March 17, 2006 p. 1 Outline of the talk Hopping transport Models for
More informationMeteorology 6150 Cloud System Modeling
Meteorology 6150 Cloud System Modeling Steve Krueger Spring 2009 1 Fundamental Equations 1.1 The Basic Equations 1.1.1 Equation of motion The movement of air in the atmosphere is governed by Newton s Second
More informationChemistry 532 Practice Final Exam Fall 2012 Solutions
Chemistry 53 Practice Final Exam Fall Solutions x e ax dx π a 3/ ; π sin 3 xdx 4 3 π cos nx dx π; sin θ cos θ + K x n e ax dx n! a n+ ; r r r r ˆL h r ˆL z h i φ ˆL x i hsin φ + cot θ cos φ θ φ ) ˆLy i
More informationA Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models
A Strategy for the Development of Coupled Ocean- Atmosphere Discontinuous Galerkin Models Frank Giraldo Department of Applied Math, Naval Postgraduate School, Monterey CA 93943 Collaborators: Jim Kelly
More informationAnalytical Mechanics: Elastic Deformation
Analytical Mechanics: Elastic Deformation Shinichi Hirai Dept. Robotics, Ritsumeikan Univ. Shinichi Hirai (Dept. Robotics, Ritsumeikan Univ.) Analytical Mechanics: Elastic Deformation 1 / 60 Agenda Agenda
More informationGFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability
GFD 2012 Lecture 1: Dynamics of Coherent Structures and their Impact on Transport and Predictability Jeffrey B. Weiss; notes by Duncan Hewitt and Pedram Hassanzadeh June 18, 2012 1 Introduction 1.1 What
More informationModelling the atmosphere. Hennie Kelder University of Technology Eindhoven
Modelling the atmosphere Hennie Kelder University of Technology Eindhoven Content Basics of the atmosphere Atmospheric dynamics Large scale circulation Planetary waves Brewer-Dobson circulation Some Key
More informationLecture 11: Internal solitary waves in the ocean
Lecture 11: Internal solitary waves in the ocean Lecturer: Roger Grimshaw. Write-up: Yiping Ma. June 19, 2009 1 Introduction In Lecture 6, we sketched a derivation of the KdV equation applicable to internal
More informationThe Evolution of Large-Amplitude Internal Gravity Wavepackets
The Evolution of Large-Amplitude Internal Gravity Wavepackets Sutherland, Bruce R. and Brown, Geoffrey L. University of Alberta Environmental and Industrial Fluid Dynamics Laboratory Edmonton, Alberta,
More informationAppendix A Notational Conventions
Appendix A Notational Conventions Throughout the book we use Einstein s implicit summation convention: repeated indices in an expression are automatically summed over. We work in natural units where the
More informationNOTES ON NEWLANDER-NIRENBERG THEOREM XU WANG
NOTES ON NEWLANDER-NIRENBERG THEOREM XU WANG Abstract. In this short note we shall recall the classical Newler-Nirenberg theorem its vector bundle version. We shall also recall an L 2 -Hörmer-proof given
More informationMagnetic wells in dimension three
Magnetic wells in dimension three Yuri A. Kordyukov joint with Bernard Helffer & Nicolas Raymond & San Vũ Ngọc Magnetic Fields and Semiclassical Analysis Rennes, May 21, 2015 Yuri A. Kordyukov (Ufa) Magnetic
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationFINAL EXAM MATH303 Theory of Ordinary Differential Equations. Spring dx dt = x + 3y dy dt = x y.
FINAL EXAM MATH0 Theory of Ordinary Differential Equations There are 5 problems on 2 pages. Spring 2009. 25 points Consider the linear plane autonomous system x + y x y. Find a fundamental matrix of the
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
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 informationMath 4263 Homework Set 1
Homework Set 1 1. Solve the following PDE/BVP 2. Solve the following PDE/BVP 2u t + 3u x = 0 u (x, 0) = sin (x) u x + e x u y = 0 u (0, y) = y 2 3. (a) Find the curves γ : t (x (t), y (t)) such that that
More informationEstimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations.
Estimation of State Noise for the Ensemble Kalman filter algorithm for 2D shallow water equations. May 6, 2009 Motivation Constitutive Equations EnKF algorithm Some results Method Navier Stokes equations
More informationTheory of Cosmological Perturbations
Theory of Cosmological Perturbations Part III CMB anisotropy 1. Photon propagation equation Definitions Lorentz-invariant distribution function: fp µ, x µ ) Lorentz-invariant volume element on momentum
More informationIdentification Methods for Structural Systems. Prof. Dr. Eleni Chatzi Lecture March, 2016
Prof. Dr. Eleni Chatzi Lecture 4-09. March, 2016 Fundamentals Overview Multiple DOF Systems State-space Formulation Eigenvalue Analysis The Mode Superposition Method The effect of Damping on Structural
More informationHorizontal reduction of pressure to mean sea level
Horizontal reduction of pressure to mean sea level Henrik Feddersen, DMI email: hf@dmi.dk October 4 Abstract. A horizontal reduction of pressure to mean sea level used in DMI-HIRLAM is documented. This
More informationChapter 1. Governing Equations of GFD. 1.1 Mass continuity
Chapter 1 Governing Equations of GFD The fluid dynamical governing equations consist of an equation for mass continuity, one for the momentum budget, and one or more additional equations to account for
More informationGauge-Higgs Unification and the LHC
Gauge-Higgs Unification and the LHC If the Higgs boson is 124 or 126 or? GeV with SM couplings, Explain SM Higgs. with non-sm couplings, is not seen at LHC, Higgs is stable. Higgs does not exist. 2 If
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationApplied Mathematics Masters Examination Fall 2016, August 18, 1 4 pm.
Applied Mathematics Masters Examination Fall 16, August 18, 1 4 pm. Each of the fifteen numbered questions is worth points. All questions will be graded, but your score for the examination will be the
More informationSolvability of the Incompressible Navier-Stokes Equations on a Moving Domain with Surface Tension. 1. Introduction
Solvability of the Incompressible Navier-Stokes Equations on a Moving Domain with Surface Tension Hantaek Bae Courant Institute of Mathematical Sciences, New York University 251 Mercer Street, New York,
More informationIMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS
Approved for public release; distribution is unlimited IMPACTS OF SIGMA COORDINATES ON THE EULER AND NAVIER-STOKES EQUATIONS USING CONTINUOUS/DISCONTINUOUS GALERKIN METHODS Sean L. Gibbons Captain, United
More informationUsing hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling
Using hybrid isentropic-sigma vertical coordinates in nonhydrostatic atmospheric modeling Michael D. Toy NCAR Earth System Laboratory/ Advanced Study Program Boulder, Colorado, USA 9th International SRNWP
More informationAccelerator Physics Homework #3 P470 (Problems: 1-5)
Accelerator Physics Homework #3 P470 (Problems: -5). Particle motion in the presence of magnetic field errors is (Sect. II.2) y + K(s)y = B Bρ, where y stands for either x or z. Here B = B z for x motion,
More informationCandidates must show on each answer book the type of calculator used. Log Tables, Statistical Tables and Graph Paper are available on request.
UNIVERSITY OF EAST ANGLIA School of Mathematics Spring Semester Examination 2004 FLUID DYNAMICS Time allowed: 3 hours Attempt Question 1 and FOUR other questions. Candidates must show on each answer book
More informationPhysics 506 Winter 2004
Physics 506 Winter 004 G. Raithel January 6, 004 Disclaimer: The purpose of these notes is to provide you with a general list of topics that were covered in class. The notes are not a substitute for reading
More informationNumerical solution of the Schrödinger equation on unbounded domains
Numerical solution of the Schrödinger equation on unbounded domains Maike Schulte Institute for Numerical and Applied Mathematics Westfälische Wilhelms-Universität Münster Cetraro, September 6 OUTLINE
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 informationNumerics and Control of PDEs Lecture 7. IFCAM IISc Bangalore. Feedback stabilization of a 1D nonlinear model
1/3 Numerics and Control of PDEs Lecture 7 IFCAM IISc Bangalore July August, 13 Feedback stabilization of a 1D nonlinear model Mythily R., Praveen C., Jean-Pierre R. /3 Plan of lecture 7 1. The nonlinear
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationA2 Assignment lambda Cover Sheet. Ready. Done BP. Question. Aa C4 Integration 1 1. C4 Integration 3
A Assignment lambda Cover Sheet Name: Question Done BP Ready Topic Comment Drill Mock Exam Aa C4 Integration sin x+ x+ c 4 Ab C4 Integration e x + c Ac C4 Integration ln x 5 + c Ba C Show root change of
More informationNONLINEAR CONTINUUM FORMULATIONS CONTENTS
NONLINEAR CONTINUUM FORMULATIONS CONTENTS Introduction to nonlinear continuum mechanics Descriptions of motion Measures of stresses and strains Updated and Total Lagrangian formulations Continuum shell
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 information9 The conservation theorems: Lecture 23
9 The conservation theorems: Lecture 23 9.1 Energy Conservation (a) For energy to be conserved we expect that the total energy density (energy per volume ) u tot to obey a conservation law t u tot + i
More informationSimulating orographic precipitation: Sensitivity to physics parameterizations and model numerics
Simulating orographic precipitation: Sensitivity to physics parameterizations and model numerics 2nd COPS-Meeting, 27 June 2005 Günther Zängl Overview A highly idealized test of numerical model errors
More informationSound Propagation in the Nocturnal Boundary Layer. Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert
Sound Propagation in the Nocturnal Boundary Layer Roger Waxler Carrick Talmadge Xiao Di Kenneth Gilbert The Propagation of Sound Outdoors (over flat ground) The atmosphere is a gas under the influence
More information1 Continuous-time chains, finite state space
Université Paris Diderot 208 Markov chains Exercises 3 Continuous-time chains, finite state space Exercise Consider a continuous-time taking values in {, 2, 3}, with generator 2 2. 2 2 0. Draw the diagramm
More informationActive Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: Hydrodynamics of SP Hard Rods
Active Matter Lectures for the 2011 ICTP School on Mathematics and Physics of Soft and Biological Matter Lecture 3: of SP Hard Rods M. Cristina Marchetti Syracuse University Baskaran & MCM, PRE 77 (2008);
More informationTrapped ghost wormholes and regular black holes. The stability problem
Trapped ghost wormholes and regular black holes. The stability problem Kirill Bronnikov in collab. with Sergei Bolokhov, Arislan Makhmudov, Milena Skvortsova (VNIIMS, Moscow; RUDN University, Moscow; MEPhI,
More informationThe KPP minimal speed within large drift in two dimensions
The KPP minimal speed within large drift in two dimensions Mohammad El Smaily Joint work with Stéphane Kirsch University of British Columbia & Pacific Institute for the Mathematical Sciences Banff, March-2010
More informationBackground Error Covariance Modelling
Background Error Covariance Modelling Mike Fisher Slide 1 Outline Diagnosing the Statistics of Background Error using Ensembles of Analyses Modelling the Statistics in Spectral Space - Relaxing constraints
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More information