The challenges of Linearized Physics (LP) in Data Assimilation (DA)
|
|
- Justin McBride
- 5 years ago
- Views:
Transcription
1 The challenges f Linearized Physics (LP) in Data Assimilatin (DA) Philippe Lpez, ECMWF with thanks t Marta aniskvá (ECMWF). Outline Brief reminder f what DA is abut. Physics in (4D-Var) DA: Why, where and hw? Cnstraints specific t LP. Benefits f LP. Summary and prspects. ECMWF Annual Seminar n Physical Prcesses, Reading, UK, 1-4 September 2015
2 The purpse and ingredients f Data Assimilatin * The purpse f DA is t merge infrmatin cming frm bservatins with a priri infrmatin cming frm a frecast mdel t btain an ptimal 3D representatin f the atmspheric state at a given time (= the analysis ). This 3D analysis then prvides initial cnditins t the numerical frecast mdel. * The main ingredients f DA are: a set f bservatins available ver a perid f typically a few hurs. a previus shrt-range frecast frm the NWP mdel ( backgrund infrmatin). sme statistical descriptin f the errrs f bth bservatins and mdel backgrund. a data assimilatin methd (e.g. variatinal DA: 3D-Var, 4D-Var; EnKF, ). This presentatin will fcus n the use f physical parameterizatins in 4D-Var DA.
3 mdel state analysis time t a 4D-Var All bservatins y are cnsidered at their actual time initial time t 0 y Mdel trajectry frm crrected initial state x a x b Mdel trajectry frm first guess x b Time (UTC) 12-hur assimilatin windw 4D-Var prduces the analysis (x a ) which minimizes the distance t a set f available bservatins (y ) and t sme a priri backgrund infrmatin frm the mdel (x b ), given the respective errrs f bservatins and mdel backgrund.
4 4D-Var Incremental 4D-Var aims at minimizing the fllwing cst functin: δx 0 = 1 2 δx 0 T B δx where: x 0 = x 0 x b 0 (increment; at lwer reslutin). B = backgrund errr cvariance matrix. d i = y i H(M i [x b 0]) (innvatin vectr; at high reslutin). M i = tangent-linear f frecast mdel (t 0 t i ). H = tangent-linear f bservatin peratr. i = time index (4D-Var windw is split in n intervals). = bservatin errr cvariance matrix. R i n i=1 (HM i δx 0 d i ) T (HM i δx 0 d i ) R i 1 δx 0 n 1 T T 1 B δx0 Mi [ ti,t0] H R i HM i δx0 i1 d i Adjint f frecast mdel with simplified linearized physics (simplified: t reduce cmputatinal cst and t avid nnlinear prcesses)
5 In ther wrds, physical parameterizatins are used in 4D-Var DA: * t describe the time evlutin f the mdel state ver the assimilatin windw as accurately as pssible. * t cnvert the mdel state variables (typ. T, wind, humidity, Psurf) int bserved equivalents s that the differences bsmdel can be cmputed (at the time f each bservatin). Side remark: One advantage f the incremental apprach is that the minimizatin f the cst functin can be run at lwer reslutin than the trajectry cmputatins (at ECMWF: 80 km versus 16 km).
6 The chice f physical parameterizatins will affect the 4D-Var results. Example: M: input = mdel state (T,q v ) utput = surface cnvective rainfall rate. acbians f cnvective surface rainfall rate w.r.t. input T and q v M T M q v M T M q v frm Marécal and Mahfuf (2002) Betts-Miller (adjustment scheme) Tiedtke (ECMWF s per mass-flux scheme)
7 Testing the tangent-linear cde The crrectness f the tangent-linear mdel must be assessed by checking that the first-rder Taylr apprximatin is valid: δx M x + λ δx M x lim λ 0 λ Mδx = 1 Example f utput frm a successful tangent-linear test: Tiny perturbatins Machine precisin reached Larger perturbatins Imprvement when perturbatin size decreases
8 Testing the adjint cde The crrectness f the adjint mdel needs t be assessed by checking that it satisfies the mathematical relatinship: δx, δy Mδx, δy = δx, M T δy where M is the tangent-linear mdel and M T is the adjint mdel. Example f utput frm a successful adjint test: <M x, y> = E-01 < x, M T y> = E-01 The difference is times the zer f the machine The adjint test shuld be crrect at the level f machine precisin (typ. 13 t 15 digits fr the entire mdel). Otherwise there must be a bug in the cde!
9 Linearity assumptin Variatinal assimilatin is based n the strng assumptin that the analysis is perfrmed in a (quasi-)linear framewrk. Hwever, physical prcesses are ften nnlinear and their parameterizatins invlve discntinuities r nn-differentiable functins (e.g. switches r threshlds). The trickiest parameterizatins are cnvectin, large-scale clud prcesses and vertical diffusin. Regularizatin needs t be applied: smthing f functins, reductin f sme perturbatins. y riginal tangent in x 0 Precipitatin frmatin rate 0 Dy (nnlinear) new tangent in x 0 Dy (tangent-linear) x 0 Dx (finite size perturbatin) x Clud water amunt
10 An example f spurius TL nise caused by a threshld in the autcnversin frmulatin f the large-scale clud scheme. ~700 hpa znal wind increments [m/s] frm 12h mdel integratin. Nnlinear finite difference: M(x+x) M(x) Thursday 15 March UTC ECMWF Frecast t+12 VT: Friday 16 March UTC Mdel Level 44 **u-velcity 12 Tangent-linear integratin: Mx Thursday 15 March UTC ECMWF Frecast t+12 VT: Friday 16 March UTC Mdel Level **u-velcity frm M. aniskvá -12 with perturbatin reductin in autcnversin -12
11 Linearized physics package used in ECMWF s peratinal 4D-Var (1) Main simplificatins/regularizatins with respect t full nnlinear mdel are highlighted in red. Large-scale cndensatin scheme: [Tmpkins and aniskvá 2004] - based n a unifrm PDF t describe subgrid-scale fluctuatins f ttal water. - melting f snw included. - precipitatin evapratin included. - reductin f clud fractin perturbatin and in autcnversin f clud int rain. Cnvectin scheme: [Lpez and Mreau 2005] - mass-flux apprach [Tiedtke 1989]. - deep cnvectin (CAPE clsure) and shallw cnvectin (q-cnvergence) are treated. - perturbatins f all cnvective quantities are included. - cupling with clud scheme thrugh detrainment f liquid water frm updraught. - sme perturbatins (buyancy, initial updraught vertical velcity) are reduced. Radiatin: TL and AD f lngwave and shrtwave radiatin available [aniskvá et al. 2002] - shrtwave: based n Mrcrette (1991), nly 2 spectral intervals (instead f 6 in nnlinear versin). - lngwave: based n Mrcrette (1989), called every 2 hurs nly.
12 Linearized physics package used in ECMWF s peratinal 4D-Var (2) Vertical diffusin: [aniskvá and Lpez 2013] - mixing in the surface and planetary bundary layers. - based n K-thery and Blackadar mixing length. - exchange cefficients based n: Luis et al. [1982] in stable cnditins, Mnin-Obukhv in unstable cnditins. - mixed-layer parametrizatin and PBL tp entrainment. - Perturbatins f exchange cefficients are smthed (esp. near the surface). Surface scheme: [aniskvá 2014, pers. cmm.] - evlutin f sil, snw and sea-ice temperature. - perturbatin reductin in snw phase change frmulatin. Orgraphic gravity wave drag: [Mahfuf 1999] - subgrid-scale rgraphic effects [Ltt and Miller 1997]. - nly lw-level blcking part is used. Nn-rgraphic gravity wave drag: [Or et al. 2010] - istrpic spectrum f nn-rgraphic gravity waves [Scincca 2003]. - Perturbatins f utput wind tendencies belw 200 hpa reset t zer.
13 Impact f linearized physics n TL apprximatin (1) Znal mean crss-sectin f change in TL errr when TL includes: VDIF + rg. GWD + SURF Relative t adiabatic TL run (50-km resl.; 20 runs; after 12h integr.). Blue = TL errr reductin = Temperature
14 Impact f linearized physics n TL apprximatin (2) Znal mean crss-sectin f change in TL errr when TL includes: VDIF + rg. GWD + SURF + RAD Relative t adiabatic TL run (50-km resl.; 20 runs; after 12h integr.). Blue = TL errr reductin = Temperature
15 Impact f linearized physics n TL apprximatin (3) Znal mean crss-sectin f change in TL errr when TL includes: VDIF + rg. GWD + SURF + RAD + nn-rg GWD + mist physics Relative t adiabatic TL run (50-km resl.; 20 runs; after 12h integr.). Blue = TL errr reductin = Temperature
16 Impact f ECMWF linearized physics n frecast scres Cmparisn f tw T511 L91 4D-Var 3-mnth experiments with & withut full linearized physics: Relative change in frecast anmaly crrelatin. > 0 = aniskvá and Lpez (2013)
17 Influence f integratin length and reslutin n linearity assumptin Cmparisn f pairs f ppsite twin experiments using ECMWF s nnlinear mdel. Time evlutin f crrelatin between M(x+x)M(x) and M(xx)M(x). Inspired frm Walser et al. (2004). P 100 hpa z 10 m The linearity assumptin becmes less valid with integratin time and reslutin, especially clse t the surface and when physics is activated.
18 Summary and prspects (1) Linearized physical parameterizatins have becme essential cmpnents f variatinal data assimilatin systems: Better representatin f the evlutin f the atmspheric state during the minimizatin f the cst functin (via the adjint mdel integratin). Extractin f infrmatin frm bservatins that are strngly affected by physical prcesses (e.g. by cluds r precipitatin). Hwever, there are sme limitatins t the LP apprach: 1) Theretical: The dmain f validity f the linear hypthesis shrinks with increasing reslutin and integratin length. 2) Technical: Linearized mdels require sustained & time-cnsuming attentin: Testing tangent-linear apprximatin and adjint cde. Regularizatins / simplificatins t eliminate any surce f instability. Revisins t ensure gd match with reference nn-linear frecast mdel.
19 Summary and prspects (2) In practice, it all cmes dwn t achieving the best cmprmise between: Realism Cst Linearity Alternative data assimilatin methds exist that d nt require the develpment f linearized cde, but s far nne f them has been able t utperfrm 4D-Var, especially in glbal mdels: Ensemble Kalman Filter (EnKF; still relies n the linearity assumptin), Particle filters (difficult t implement fr high-dimensinal prblems). S what is the future f LP?
20 Frm a small challenge
21 t a much bigger challenge
22 Summary and prspects (3) Eventually, it might becme impractical r even impssible t make LP wrk efficiently at reslutins f a few kilmetres, even if the linearity cnstraint can be relaxed (e.g. by using shrter 4D-Var windw r weak-cnstraint 4D- Var).?
23
24 Simple example f tangent-linear and adjint cdes. Simplified nnlinear cde: Z = a / X 2 + b Y lg(w) Tangent-linear cde: Z = (2 a / X 3 ) X + b lg(w) Y + (b Y / W) W Adjint cde: X* = 0 Y* = 0 W* = 0 X* = X* (2 a / X 3 ) Z* Y* = Y* + b lg(w) Z* W = W* + (b Y / W) Z* Z* = 0
25 Example f nnlinear, tangent-linear and adjint peratrs in the cntext f 4D-Var. Nnlinear peratrs are applied t full fields: mdel initial state x 0 T qv u v P s T q u v P q q s liq ice time t i bservatin equivalent time t 0 = satellite cludy radiances v time t i M H ~ y i Rad Rad Rad ch1 ch2 ch3 M = frecast mdel with physics H = radiative transfer mdel
26 Tangent-linear peratrs are applied t perturbatins: ~ ], [ ch ch ch i ice liq s v i s v Rad Rad Rad q q P v u q T P v u q T t t δy δx H M time t i time t 0 Adjint peratrs are applied t the cst functin gradient: s v i T ice liq s v T ch ch ch P v u q T t t q q P v u q T Rad Rad Rad i / / / / / ], [ / / / / / / / / / / ~ x y M H time t i time t 0
27 Illustratin f discntinuity effect n cst functin shape: Mdel backgrund = {T b, q b }; Observatin = RR bs Simple parametrizatin f rain rate: RR = {q q sat (T)} if q > q sat (T), 0 therwise 1 2 T T T b q q q b [ q qsat( T)] RR RR bs bs 2 q Saturated backgrund q Dry backgrund {T b, q b } min OK N cnvergence! Single minimum f cst functin T Several lcal minima f cst functin T
28 M Frecasts Climate runs N OK? Timing: Yes Pure cding N M M(x+x)M(x) M x M* <M x,y> = <x,m*y> OK? Yes N OK? Yes Debugging and testing (incl. regularizatin) N N OK? OK? Singular Vectrs (EPS) 4D-Var (minim) NL TL AD APPL Stages in the develpment f a new linearized parameterizatin
29 Tangent-linear apprximatin and assciated errr The tangent-linear apprximatin is assessed by cmparing the difference between tw integratins f the full nn-linear mdel: M(xx) M(x) with an integratin f the tangent-linear mdel frm the same initial perturbatin (x = x ana x bg ): Mx The TL errr is then defined as: = < M(xx) M(x) Mx > where < > dentes spatial averaging (e.g. znally r glbally).
30 TOA Impact f linearized physics n TL apprximatin (4) Mean vertical prfile f change in TL errr fr T, U, V and Q when full linearized physics is included in TL cmputatins. Relative t adiabatic TL run (50-km resl.; twenty runs, 12h integ.) < 0 = surface Inclusin f linearized physics leads t better TL apprximatin.
31 Effects f mdel nnlinearities versus effects f reslutin Cmparisn f pairs f ppsite twin experiments using ECMWF s nnlinear mdel. Time evlutin f crrelatin between M(x+x)M(x) and M(xx)M(x). Inspired frm Walser et al. (2004). P 100 hpa z 10 m The effects f nnlinearities are f cmparable magnitude with thse f reslutin, except near the surface when physics is activated (the frmer effects then dminate).
Parametrizations in Data Assimilation
Parametrizatins in Data Assimilatin ECMWF raining Curse 0-4 Mar 04 Philippe Lpez Physical Aspects Sectin, Resear Department, ECMWF (Rm 00) Parametrizatins in Data Assimilatin Intrductin Why are physical
More informationComparison of hybrid ensemble-4dvar with EnKF and 4DVar for regional-scale data assimilation
Cmparisn f hybrid ensemble-4dvar with EnKF and 4DVar fr reginal-scale data assimilatin Jn Pterjy and Fuqing Zhang Department f Meterlgy The Pennsylvania State University Wednesday 18 th December, 2013
More informationApplication of a Coupled EnKF and 4DVar Data Assimilation Method in the Study of Tropical Cyclone Genesis
! Applicatin f a Cupled EnKF and 4DVar Data Assimilatin Methd in the Study f Trpical Cyclne Genesis Ashfrd Reyes, Gregry Jenkins, Jnathan PterJy and Fuqing Zhang Hward University Prgram fr Atmspheric Sciences,
More informationHow do we solve it and what does the solution look like?
Hw d we slve it and what des the slutin l lie? KF/PFs ffer slutins t dynamical systems, nnlinear in general, using predictin and update as data becmes available. Tracing in time r space ffers an ideal
More informationAPPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL
JP2.11 APPLICATION OF THE BRATSETH SCHEME FOR HIGH LATITUDE INTERMITTENT DATA ASSIMILATION USING THE PSU/NCAR MM5 MESOSCALE MODEL Xingang Fan * and Jeffrey S. Tilley University f Alaska Fairbanks, Fairbanks,
More informationOperational Use of the Model Crocus
Operatinal Use f the Mdel Crcus by French Avalanche Frecast Services E.Brun Meterlgie Natinale Centre d'etudes de la Neige BP 44 Dmaine Universitaire 3842 St-Martin d 'eres France ntrductin Since 1971
More informationLeast Squares Optimal Filtering with Multirate Observations
Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical
More informationKinetic Model Completeness
5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins
More informationPattern Recognition 2014 Support Vector Machines
Pattern Recgnitin 2014 Supprt Vectr Machines Ad Feelders Universiteit Utrecht Ad Feelders ( Universiteit Utrecht ) Pattern Recgnitin 1 / 55 Overview 1 Separable Case 2 Kernel Functins 3 Allwing Errrs (Sft
More informationCHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.
MATH 1342 Ch. 24 April 25 and 27, 2013 Page 1 f 5 CHAPTER 24: INFERENCE IN REGRESSION Chapters 4 and 5: Relatinships between tw quantitative variables. Be able t Make a graph (scatterplt) Summarize the
More informationENSC Discrete Time Systems. Project Outline. Semester
ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding
More informationNAME TEMPERATURE AND HUMIDITY. I. Introduction
NAME TEMPERATURE AND HUMIDITY I. Intrductin Temperature is the single mst imprtant factr in determining atmspheric cnditins because it greatly influences: 1. The amunt f water vapr in the air 2. The pssibility
More informationz = Geometric height (m)
13 Z = Geptential height (m) = Lapse rate (6.5 K km -1 ) R = Gas cnstant fr dry air (287 Jkg -1 K) g = Acceleratin f gravity (9.8 ms -2 ) TS = Surface Temperature (K) p = Initial air pressure (Assumptin:
More informationhal , version 1-18 Sep 2007
Ocean Dynamics manuscript N. (will be inserted by the editr) Cmparisn f reduced-rder, sequential and variatinal data assimilatin methds in the trpical Pacific Ocean Céline Rbert 1, Eric Blay 1, Jacques
More informationTangent-linear and adjoint models in data assimilation
Tangent-linear and adjoint models in data assimilation Marta Janisková and Philippe Lopez ECMWF Thanks to: F. Váňa, M.Fielding 2018 Annual Seminar: Earth system assimilation 10-13 September 2018 Tangent-linear
More informationCOMP 551 Applied Machine Learning Lecture 9: Support Vector Machines (cont d)
COMP 551 Applied Machine Learning Lecture 9: Supprt Vectr Machines (cnt d) Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Class web page: www.cs.mcgill.ca/~hvanh2/cmp551 Unless therwise
More informationThe blessing of dimensionality for kernel methods
fr kernel methds Building classifiers in high dimensinal space Pierre Dupnt Pierre.Dupnt@ucluvain.be Classifiers define decisin surfaces in sme feature space where the data is either initially represented
More informationComputational modeling techniques
Cmputatinal mdeling techniques Lecture 11: Mdeling with systems f ODEs In Petre Department f IT, Ab Akademi http://www.users.ab.fi/ipetre/cmpmd/ Mdeling with differential equatins Mdeling strategy Fcus
More informationThe Kullback-Leibler Kernel as a Framework for Discriminant and Localized Representations for Visual Recognition
The Kullback-Leibler Kernel as a Framewrk fr Discriminant and Lcalized Representatins fr Visual Recgnitin Nun Vascncels Purdy H Pedr Mren ECE Department University f Califrnia, San Dieg HP Labs Cambridge
More informationFebruary 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA
February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal
More informationDetermining the Accuracy of Modal Parameter Estimation Methods
Determining the Accuracy f Mdal Parameter Estimatin Methds by Michael Lee Ph.D., P.E. & Mar Richardsn Ph.D. Structural Measurement Systems Milpitas, CA Abstract The mst cmmn type f mdal testing system
More informationModelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA
Mdelling f Clck Behaviur Dn Percival Applied Physics Labratry University f Washingtn Seattle, Washingtn, USA verheads and paper fr talk available at http://faculty.washingtn.edu/dbp/talks.html 1 Overview
More informationPerformance Bounds for Detect and Avoid Signal Sensing
Perfrmance unds fr Detect and Avid Signal Sensing Sam Reisenfeld Real-ime Infrmatin etwrks, University f echnlgy, Sydney, radway, SW 007, Australia samr@uts.edu.au Abstract Detect and Avid (DAA) is a Cgnitive
More informationA Matrix Representation of Panel Data
web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins
More informationFall 2013 Physics 172 Recitation 3 Momentum and Springs
Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.
More informationTHE PARTITION OF ENERGY INTO WAVES AND CURRENTS
THE PARTITION OF ENERGY INTO WAVES AND CURRENTS W. Perrie, C. Tang, Y. Hu and B.M. DeTracy Fisheries & Oceans Canada, Bedfrd Institute f Oceangraphy, Dartmuth, Nva Sctia, Canada 1. INTRODUCTION Ocean mdels
More informationEric Klein and Ning Sa
Week 12. Statistical Appraches t Netwrks: p1 and p* Wasserman and Faust Chapter 15: Statistical Analysis f Single Relatinal Netwrks There are fur tasks in psitinal analysis: 1) Define Equivalence 2) Measure
More informationResampling Methods. Chapter 5. Chapter 5 1 / 52
Resampling Methds Chapter 5 Chapter 5 1 / 52 1 51 Validatin set apprach 2 52 Crss validatin 3 53 Btstrap Chapter 5 2 / 52 Abut Resampling An imprtant statistical tl Pretending the data as ppulatin and
More informationk-nearest Neighbor How to choose k Average of k points more reliable when: Large k: noise in attributes +o o noise in class labels
Mtivating Example Memry-Based Learning Instance-Based Learning K-earest eighbr Inductive Assumptin Similar inputs map t similar utputs If nt true => learning is impssible If true => learning reduces t
More informationDifferentiation Applications 1: Related Rates
Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm
More informationSAMPLING DYNAMICAL SYSTEMS
SAMPLING DYNAMICAL SYSTEMS Melvin J. Hinich Applied Research Labratries The University f Texas at Austin Austin, TX 78713-8029, USA (512) 835-3278 (Vice) 835-3259 (Fax) hinich@mail.la.utexas.edu ABSTRACT
More informationNGSS High School Physics Domain Model
NGSS High Schl Physics Dmain Mdel Mtin and Stability: Frces and Interactins HS-PS2-1: Students will be able t analyze data t supprt the claim that Newtn s secnd law f mtin describes the mathematical relatinship
More informationModule 4: General Formulation of Electric Circuit Theory
Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated
More informationA Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture
Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu
More informationFIELD QUALITY IN ACCELERATOR MAGNETS
FIELD QUALITY IN ACCELERATOR MAGNETS S. Russenschuck CERN, 1211 Geneva 23, Switzerland Abstract The field quality in the supercnducting magnets is expressed in terms f the cefficients f the Furier series
More informationA Scalable Recurrent Neural Network Framework for Model-free
A Scalable Recurrent Neural Netwrk Framewrk fr Mdel-free POMDPs April 3, 2007 Zhenzhen Liu, Itamar Elhanany Machine Intelligence Lab Department f Electrical and Cmputer Engineering The University f Tennessee
More informationand the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:
Algrithm fr Estimating R and R - (David Sandwell, SIO, August 4, 2006) Azimith cmpressin invlves the alignment f successive eches t be fcused n a pint target Let s be the slw time alng the satellite track
More informationChapter 11: Atmosphere
Chapter 11: Atmsphere Sectin 1: Atmspheric Basics Objectives 1. Describe the cmpsitin f the atmsphere. 2. Cmpare and cntrast the varius layers f the atmsphere. 3. Identify three methds f transferring energy
More informationLyapunov Stability Stability of Equilibrium Points
Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),
More informationWhat is Statistical Learning?
What is Statistical Learning? Sales 5 10 15 20 25 Sales 5 10 15 20 25 Sales 5 10 15 20 25 0 50 100 200 300 TV 0 10 20 30 40 50 Radi 0 20 40 60 80 100 Newspaper Shwn are Sales vs TV, Radi and Newspaper,
More informationNUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION
NUROP Chinese Pinyin T Chinese Character Cnversin NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION CHIA LI SHI 1 AND LUA KIM TENG 2 Schl f Cmputing, Natinal University f Singapre 3 Science
More informationCOMP 551 Applied Machine Learning Lecture 11: Support Vector Machines
COMP 551 Applied Machine Learning Lecture 11: Supprt Vectr Machines Instructr: (jpineau@cs.mcgill.ca) Class web page: www.cs.mcgill.ca/~jpineau/cmp551 Unless therwise nted, all material psted fr this curse
More informationPSU GISPOPSCI June 2011 Ordinary Least Squares & Spatial Linear Regression in GeoDa
There are tw parts t this lab. The first is intended t demnstrate hw t request and interpret the spatial diagnstics f a standard OLS regressin mdel using GeDa. The diagnstics prvide infrmatin abut the
More informationPerfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart
Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Sandy D. Balkin Dennis K. J. Lin y Pennsylvania State University, University Park, PA 16802 Sandy Balkin is a graduate student
More informationThermodynamics and Equilibrium
Thermdynamics and Equilibrium Thermdynamics Thermdynamics is the study f the relatinship between heat and ther frms f energy in a chemical r physical prcess. We intrduced the thermdynamic prperty f enthalpy,
More informationMaterials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals of Diffusion
Materials Engineering 272-C Fall 2001, Lecture 7 & 8 Fundamentals f Diffusin Diffusin: Transprt in a slid, liquid, r gas driven by a cncentratin gradient (r, in the case f mass transprt, a chemical ptential
More informationMath Foundations 20 Work Plan
Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant
More informationLecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff
Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeff Reading: Chapter 2 STATS 202: Data mining and analysis September 27, 2017 1 / 20 Supervised vs. unsupervised learning In unsupervised
More information(1.1) V which contains charges. If a charge density ρ, is defined as the limit of the ratio of the charge contained. 0, and if a force density f
1.0 Review f Electrmagnetic Field Thery Selected aspects f electrmagnetic thery are reviewed in this sectin, with emphasis n cncepts which are useful in understanding magnet design. Detailed, rigrus treatments
More informationReview Problems 3. Four FIR Filter Types
Review Prblems 3 Fur FIR Filter Types Fur types f FIR linear phase digital filters have cefficients h(n fr 0 n M. They are defined as fllws: Type I: h(n = h(m-n and M even. Type II: h(n = h(m-n and M dd.
More informationWRITING THE REPORT. Organizing the report. Title Page. Table of Contents
WRITING THE REPORT Organizing the reprt Mst reprts shuld be rganized in the fllwing manner. Smetime there is a valid reasn t include extra chapters in within the bdy f the reprt. 1. Title page 2. Executive
More informationEmphases in Common Core Standards for Mathematical Content Kindergarten High School
Emphases in Cmmn Cre Standards fr Mathematical Cntent Kindergarten High Schl Cntent Emphases by Cluster March 12, 2012 Describes cntent emphases in the standards at the cluster level fr each grade. These
More informationLecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff
Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeff Reading: Chapter 2 STATS 202: Data mining and analysis September 27, 2017 1 / 20 Supervised vs. unsupervised learning In unsupervised
More informationCOMP 551 Applied Machine Learning Lecture 5: Generative models for linear classification
COMP 551 Applied Machine Learning Lecture 5: Generative mdels fr linear classificatin Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Jelle Pineau Class web page: www.cs.mcgill.ca/~hvanh2/cmp551
More informationDistributions, spatial statistics and a Bayesian perspective
Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics
More informationDepartment of Electrical Engineering, University of Waterloo. Introduction
Sectin 4: Sequential Circuits Majr Tpics Types f sequential circuits Flip-flps Analysis f clcked sequential circuits Mre and Mealy machines Design f clcked sequential circuits State transitin design methd
More informationDESIGN OPTIMIZATION OF HIGH-LIFT CONFIGURATIONS USING A VISCOUS ADJOINT-BASED METHOD
DESIGN OPTIMIZATION OF HIGH-LIFT CONFIGURATIONS USING A VISCOUS ADJOINT-BASED METHOD Sangh Kim Stanfrd University Juan J. Alns Stanfrd University Antny Jamesn Stanfrd University 40th AIAA Aerspace Sciences
More information5 th grade Common Core Standards
5 th grade Cmmn Cre Standards In Grade 5, instructinal time shuld fcus n three critical areas: (1) develping fluency with additin and subtractin f fractins, and develping understanding f the multiplicatin
More informationLHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers
LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the
More informationBootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >
Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);
More information7.0 Heat Transfer in an External Laminar Boundary Layer
7.0 Heat ransfer in an Eternal Laminar Bundary Layer 7. Intrductin In this chapter, we will assume: ) hat the fluid prperties are cnstant and unaffected by temperature variatins. ) he thermal & mmentum
More informationHypothesis Tests for One Population Mean
Hypthesis Tests fr One Ppulatin Mean Chapter 9 Ala Abdelbaki Objective Objective: T estimate the value f ne ppulatin mean Inferential statistics using statistics in rder t estimate parameters We will be
More informationCAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank
CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal
More informationBuilding to Transformations on Coordinate Axis Grade 5: Geometry Graph points on the coordinate plane to solve real-world and mathematical problems.
Building t Transfrmatins n Crdinate Axis Grade 5: Gemetry Graph pints n the crdinate plane t slve real-wrld and mathematical prblems. 5.G.1. Use a pair f perpendicular number lines, called axes, t define
More informationarxiv: v1 [physics.comp-ph] 21 Feb 2018
arxiv:182.7486v1 [physics.cmp-ph] 21 Feb 218 rspa.ryalscietypublishing.rg Research Article submitted t jurnal Subject Areas: Mechanical Engineering Keywrds: Data-driven frecasting, Lng-Shrt Term Memry,
More informationChecking the resolved resonance region in EXFOR database
Checking the reslved resnance regin in EXFOR database Gttfried Bertn Sciété de Calcul Mathématique (SCM) Oscar Cabells OECD/NEA Data Bank JEFF Meetings - Sessin JEFF Experiments Nvember 0-4, 017 Bulgne-Billancurt,
More informationObservation quality control: methodology and applications
Observatin qualit cntrl: methdlg and applicatins Pierre Gauthier Department f Earth and Atmspheric Sciences Université du Québec à Mntréal CANADA) Presentatin at the 010 ESA Earth Observatin Summer Schl
More informationPhysics 2010 Motion with Constant Acceleration Experiment 1
. Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin
More informationInternal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.
Sectin 7 Mdel Assessment This sectin is based n Stck and Watsn s Chapter 9. Internal vs. external validity Internal validity refers t whether the analysis is valid fr the ppulatin and sample being studied.
More informationTechnical Bulletin. Generation Interconnection Procedures. Revisions to Cluster 4, Phase 1 Study Methodology
Technical Bulletin Generatin Intercnnectin Prcedures Revisins t Cluster 4, Phase 1 Study Methdlgy Release Date: Octber 20, 2011 (Finalizatin f the Draft Technical Bulletin released n September 19, 2011)
More informationRevision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax
.7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical
More informationPart 3 Introduction to statistical classification techniques
Part 3 Intrductin t statistical classificatin techniques Machine Learning, Part 3, March 07 Fabi Rli Preamble ØIn Part we have seen that if we knw: Psterir prbabilities P(ω i / ) Or the equivalent terms
More informationMath 302 Learning Objectives
Multivariable Calculus (Part I) 13.1 Vectrs in Three-Dimensinal Space Math 302 Learning Objectives Plt pints in three-dimensinal space. Find the distance between tw pints in three-dimensinal space. Write
More informationThe general linear model and Statistical Parametric Mapping I: Introduction to the GLM
The general linear mdel and Statistical Parametric Mapping I: Intrductin t the GLM Alexa Mrcm and Stefan Kiebel, Rik Hensn, Andrew Hlmes & J-B J Pline Overview Intrductin Essential cncepts Mdelling Design
More informationGeneral Chemistry II, Unit I: Study Guide (part I)
1 General Chemistry II, Unit I: Study Guide (part I) CDS Chapter 14: Physical Prperties f Gases Observatin 1: Pressure- Vlume Measurements n Gases The spring f air is measured as pressure, defined as the
More informationNWC SAF ENTERING A NEW PHASE
NWC SAF ENTERING A NEW PHASE Pilar Fernández Agencia Estatal de Meterlgía (AEMET) Lenard Priet Castr, 8; 28040 Madrid, Spain Phne: +34 915 819 654, Fax: +34 915 819 767 E-mail: mafernandeza@aemet.es Abstract
More information11. DUAL NATURE OF RADIATION AND MATTER
11. DUAL NATURE OF RADIATION AND MATTER Very shrt answer and shrt answer questins 1. Define wrk functin f a metal? The minimum energy required fr an electrn t escape frm the metal surface is called the
More informationEnhancing Performance of MLP/RBF Neural Classifiers via an Multivariate Data Distribution Scheme
Enhancing Perfrmance f / Neural Classifiers via an Multivariate Data Distributin Scheme Halis Altun, Gökhan Gelen Nigde University, Electrical and Electrnics Engineering Department Nigde, Turkey haltun@nigde.edu.tr
More informationEDA Engineering Design & Analysis Ltd
EDA Engineering Design & Analysis Ltd THE FINITE ELEMENT METHOD A shrt tutrial giving an verview f the histry, thery and applicatin f the finite element methd. Intrductin Value f FEM Applicatins Elements
More informationNumerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment
Presented at the COMSOL Cnference 2008 Hannver University f Parma Department f Industrial Engineering Numerical Simulatin f the Thermal Respsne Test Within the Cmsl Multiphysics Envirnment Authr : C. Crradi,
More informationx x
Mdeling the Dynamics f Life: Calculus and Prbability fr Life Scientists Frederick R. Adler cfrederick R. Adler, Department f Mathematics and Department f Bilgy, University f Utah, Salt Lake City, Utah
More informationComparing Several Means: ANOVA. Group Means and Grand Mean
STAT 511 ANOVA and Regressin 1 Cmparing Several Means: ANOVA Slide 1 Blue Lake snap beans were grwn in 12 pen-tp chambers which are subject t 4 treatments 3 each with O 3 and SO 2 present/absent. The ttal
More informationRelationship Between Amplifier Settling Time and Pole-Zero Placements for Second-Order Systems *
Relatinship Between Amplifier Settling Time and Ple-Zer Placements fr Secnd-Order Systems * Mark E. Schlarmann and Randall L. Geiger Iwa State University Electrical and Cmputer Engineering Department Ames,
More informationABSORPTION OF GAMMA RAYS
6 Sep 11 Gamma.1 ABSORPTIO OF GAMMA RAYS Gamma rays is the name given t high energy electrmagnetic radiatin riginating frm nuclear energy level transitins. (Typical wavelength, frequency, and energy ranges
More informationTree Structured Classifier
Tree Structured Classifier Reference: Classificatin and Regressin Trees by L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stne, Chapman & Hall, 98. A Medical Eample (CART): Predict high risk patients
More informationthe results to larger systems due to prop'erties of the projection algorithm. First, the number of hidden nodes must
M.E. Aggune, M.J. Dambrg, M.A. El-Sharkawi, R.J. Marks II and L.E. Atlas, "Dynamic and static security assessment f pwer systems using artificial neural netwrks", Prceedings f the NSF Wrkshp n Applicatins
More informationModelling of NOLM Demultiplexers Employing Optical Soliton Control Pulse
Micwave and Optical Technlgy Letters, Vl. 1, N. 3, 1999. pp. 05-08 Mdelling f NOLM Demultiplexers Emplying Optical Slitn Cntrl Pulse Z. Ghassemly, C. Y. Cheung & A. K. Ray Electrnics Research Grup, Schl
More informationLecture 02 CSE 40547/60547 Computing at the Nanoscale
PN Junctin Ntes: Lecture 02 CSE 40547/60547 Cmputing at the Nanscale Letʼs start with a (very) shrt review f semi-cnducting materials: - N-type material: Obtained by adding impurity with 5 valence elements
More informationStudy Group Report: Plate-fin Heat Exchangers: AEA Technology
Study Grup Reprt: Plate-fin Heat Exchangers: AEA Technlgy The prblem under study cncerned the apparent discrepancy between a series f experiments using a plate fin heat exchanger and the classical thery
More informationCHAPTER 5 THE RADIATIVE TRANSFER EQUATION (RTE)
CHAPTER 5 THE RADIATIVE TRANSFER EQUATION (RTE) 5.1 Derivatin f RTE Radiative transfer serves as a mechanism fr exchanging energy between the atmsphere and the underlying surface and between different
More information3.4 Shrinkage Methods Prostate Cancer Data Example (Continued) Ridge Regression
3.3.4 Prstate Cancer Data Example (Cntinued) 3.4 Shrinkage Methds 61 Table 3.3 shws the cefficients frm a number f different selectin and shrinkage methds. They are best-subset selectin using an all-subsets
More informationCHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India
CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion
.54 Neutrn Interactins and Applicatins (Spring 004) Chapter (3//04) Neutrn Diffusin References -- J. R. Lamarsh, Intrductin t Nuclear Reactr Thery (Addisn-Wesley, Reading, 966) T study neutrn diffusin
More informationSNOW AND AVALANCHES FORECASTING OVER THE ANDES MOUNTAlNS. Jose A. Vergara* Departamento de Geoflsica, Universidad de Chile
SNOW AND AVALANCHES FORECASTNG OVER THE ANDES MOUNTAlNS Jse A. Vergara* Departament de Geflsica, Universidad de Chile Rene len, Camins y Nieve, ANDNA-CODElCO, CHLE ABSTRACT: During the last three years
More informationWriting Guidelines. (Updated: November 25, 2009) Forwards
Writing Guidelines (Updated: Nvember 25, 2009) Frwards I have fund in my review f the manuscripts frm ur students and research assciates, as well as thse submitted t varius jurnals by thers that the majr
More informationLarge scale facilities
The Labratry f Atmspheric Pllutin and f Pllutin Cntrl Engineering f Atmspheric Pllutants, is part f the Department f Envirnmental Engineering, f Demcritus University f Thrace in Greece. The labratry clsely
More informationAerodynamic Separability in Tip Speed Ratio and Separability in Wind Speed- a Comparison
Jurnal f Physics: Cnference Series OPEN ACCESS Aerdynamic Separability in Tip Speed Rati and Separability in Wind Speed- a Cmparisn T cite this article: M L Gala Sants et al 14 J. Phys.: Cnf. Ser. 555
More informationSPH3U1 Lesson 06 Kinematics
PROJECTILE MOTION LEARNING GOALS Students will: Describe the mtin f an bject thrwn at arbitrary angles thrugh the air. Describe the hrizntal and vertical mtins f a prjectile. Slve prjectile mtin prblems.
More informationCollective effects in Beam Dynamics
Cllective effects in Beam Dynamics Givanni Ruml, Hannes Bartsik and Kevin Li USPAS, ne-week curse, 19-24 January, 2015 http://uspas.fnal.gv/index.shtml January 2015 USPAS lectures 2 Particle beam Nminal
More information