The Destabilization of Rossby Normal Modes by Meridional Baroclinic Shear
|
|
- Milton Hubbard
- 5 years ago
- Views:
Transcription
1 The Destabilizatin f Rssby Nrmal Mdes by Meridinal Barclinic Shear by Jseph Pedlsky Wds Hle Oceangraphic Institutin Wds Hle, MA 0543 Abstract The Rssby nrmal mdes f a tw-layer fluid in a meridinal channel f width L * are altered by the presence f a meridinal flw with a small vertical shear. The stability f the mdes in the presence f the weak shear is cnsidered. It is fund that the jint presence f the Rssby mdes and the vertical shear leads t barclinic instability even fr arbitrarily small values f the shear. The results are used t explain previus numerical calculatins f the persistent instability f meridinal flws when the rati βl D / V >1, where V is the magnitude f the shear, β is the planetary vrticity gradient and L D is the defrmatin radius. If the flw were znal it wuld be stable fr such weak shears. The grwth rates are weak when βl D / V >>1 and each unstable mde exists in a narrw range f meridinal wavenumber. The asympttic results qualitatively agree with the earlier numerical results at mderate values f the same parameter. 1. Intrductin In a recent paper Walker and Pedlsky (00), hereafter WP, examined the instability, within the tw-layer mdel, f a meridinal flw n the beta plane in a meridinal channel f width L *. If the channel were unbunded and the meridinal flw were f infinite lateral extent earlier results (Pedlsky, 1987) shw that the flw wuld always be unstable t a wave-like perturbatin independent f the znal directin (x). Such a perturbatin wuld thus be insensitive t the stabilizing presence f β. The finite width f the channel frces an x variatin in the perturbatin stream functin field and thus a meridinal velcity which senses the planetary vrticity. Nevertheless, WP fund that the meridinal shear flw was unstable fr all values f β examined, even fr thse values f βl D / V >1 (symbls have cnventinal meanings and are defined belw) fr which the flw wuld be stable if it were znal instead f meridinal. In additin, the instability extended t wave lengths shrter than the classical shrt- jpedlsky@whi.edu
2 wave cut-ff f the tw-layer mdel. WP speculated that the extended range f instability in shear and wave number was due t the destabilizatin by the presence f weak shear f the Rssby nrmal mdes present in the channel. (Nte that fr a plane wave in an unbunded regin a stability threshld des exist if the wave vectr is nt purely meridinal). Hwever, the difficulty f the numerical analysis when the shear is weak and the grwth rates are small precluded a satisfying verificatin f this hypthesis. In this paper I present an asympttic perturbatin analysis valid fr very weak shear t demnstrate the persistence f at least weak instabilities fr small values f the shear. In additin t its explanatry quality with regard t the earlier results in WP it is suggestive that the ptential energy present in the weak shear can be tapped by the Rssby mde. Althugh the basic flw cnsidered here is cnsiderably simpler than the circulatin in a cmplete subtrpical gyre it is pssible that the instability utlined here can prvide a mechanism t maintain Rssby mdes in such gyres against the inevitable effects f dissipatin whse presence has ften been cited as a reasn fr the unlikelihd fr the existence f Rssby nrmal mdes. In sectin I frmulate the basic prblem and exhibit the perturbatin analysis. Sectin 3 is a presentatin f results and sectin 4 is a brief discussin f the results and their significance.. Frmulatin Cnsider a tw layer, quasi-gestrphic mdel n the beta plane (Pedlsky, 1987). Fr simplicity we will take each f the layers t have the same basic thickness, H, in the absence f mtin. Imagine a channel f width L *, riented nrth suth, as shwn in Figure 1. There is a unifrm nrthward flw, V, in nly the upper layer f the tw layers. In rder fr such a flw t be a cnsistent slutin f the ptential vrticity equatin there must be a vrticity surce n the beta plane t maintain the flw. One can imagine a unifrm wind stress curl being respnsible fr maintaining V. Small perturbatins that are wave-like in y disturb the basic flw and the quasi-gestrphic ptential vrticity equatin is linearized in rder t describe the initial evlutin f the perturbatin field. In dimensinless units the perturbatin equatins are: (v s) [ φ 1xx l φ 1 + (φ φ 1 )]+ 1 il φ 1x + vφ 1 = 0 ( s) [ φ xx l φ + (φ 1 φ )]+ 1 il φ x vφ = 0 (.1 a,b)
3 Figure 1. The meridinal channel cntaining the flw. The upper layer basic state velcity is V and the channel width is L. In (.1) the crss channel crdinate x has been scaled with the defrmatin radius, (g ' H) 1/ / f where the reduced gravity and Crilis parameter are in standard ntatin. Subscripts x dentes derivatives with respect t x. The alng channel wave number l is similarly scaled with the inverse defrmatin radius. Bth the meridinal basic state velcity and the cmplex phase speed f the wave perturbatin, s, are scaled with the characteristic Rssby lng wave speed βl D. Fr details f the derivatin the reader is referred t WP. If the basic state shear is weak s that v = V / βl D << 1 an expansin in an asympttic series in v is suggested. First, hwever, there is a suggestin frm the numerical results f WP that the unstable mdes, if they exist, will have scales that can be shrt cmpared t a defrmatin radius and wavelengths in the y-directin that are als shrt. This suggests that the fllwing transfrmatins fr x,l and c are useful: 3
4 ξ = x / v 1/, a = lv 1/ (. a,b,c) s = v(s + vs ) where we have expanded the phase speed in a series in the small parameter v, the perturbatin equatins are nw: [ ] i a φ 1ξ + vφ 1 = 0 (1 s vs ) φ 1ξξ a φ 1 + v(φ φ 1 ) (.3a,b) ( s vs ) [ φ ξξ a φ + v(φ 1 φ )] i a φ ξ vφ = 0 Nte that in these stretched crdinates the bundaries f the channel are at ξ = 0 and ξ = L * (L D v 1/ ) L. We als nte that t lwest rder in the shear, v, the tw layers are decupled. This is cnsistent with the results f WP which shwed that unstable mdes exist fr βl D / V >1, with strikingly different crss channel scales in the x-directin. The perturbatin streamfunctin is als expanded in a series in v, as is a, thus: φ n = φ n () + vφ n (1) +..., n =1, a = a (1 + vα +...) (.4) At lwest rder this yields the prblem d (1 s ) dξ φ 1 () a () φ 1 i φ a 1ξ = 0, d ( s ) dξ φ () a φ () i φ a ξ = 0 (.5a,b) 4
5 whse slutins, subject t the bundary cnditins f vanishing stream functin at the channel walls are: φ 1 () = A 1 e iξ /[a (1 s )] sinmξ, (.6a,b) φ () = A e iξ /[a s ] sin nξ where m Mπ / L, M =1,,3... n = Nπ / L, N = 1,,3... (.7 a,b) Since L is the channel width scaled inversely with v 1/ it fllws that fr cnsistency M 1/ and N, each an integer, must be large, O ( v ) s that m and n are O(1) in accrdance with the assumptins f the asympttics. We nte that the slutins (.6) are each the Rssby nrmal mdes f the upper and lwer layers respectively in which a carrier wave is mdulated by a slutin f the Helmhltz equatin (Pedlsky, 1987). The nrmal mde structure in the upper layer is altered nly by the Dppler shift f the meridinal flw in the upper layer. The slutins (.7) reflect the fact that the hrizntal structure f the mde may be quite different frm ne layer t the ther if M and N are different. The cnditin that bth (.6 a,b) are slutins crrespnding t the same frequency r phase speed is simply that the tw slutins fr s, frm (.6b) and s = s =1 1 a (a + n ) 1/ (.8) 1 a (a + m ) 1/ (.9) frm (.6a), must yield identical results. The cnditin that s be the same in (.8) and (.9), i.e. that we are dealing with a single cmpsite mde, yields the cnditin, 5
6 a = 1 ( a + m ) 1/ + 1 ( a + n 1/. (.10) ) Fr any pair (m,n) there is a single slutin fr the (scaled) y- wave number a. Figure shws a map in the M, N, plane f a fr the case v =.0, L * =10. Althugh M and N are integers the figure treats the variables as cntinuus fr graphical clarity. It is imprtant t nte that t this rder the phase speed is strictly real s that instability will be apparent nly at the next rder in the expansin in v. At next rder the perturbatin equatins are: (1 s ) (1) φ 1ξξ a (1) [ φ 1 ] i (1) φ a 1ξ = (1 s ) () φ () φ1 [ ] () +s 1 φ 1ξξ a () [ φ 1 ] φ 1 () (.11a) +(1 s ) α a () [ φ 1 ] i αφ a 1 ξ () (1) s φ ξξ a (1) φ i (1) φ a ξ = s () (φ () 1 φ ) () +φ +s 1 φ () a () ξξ φ s a () αφ (.11b) i a αφ ξ (). It is nly at this rder that the cupling between the tw layers enters the perturbatin equatins and the ptential vrticity gradient assciated with the vertical shear is explicitly included in the analysis. If barclinic instability is t ccur these must be essential ingredients. The instability prperties will be cntained in behavir f s 1 whse imaginary part will yield the grwth rate fr the perturbatin (after multiplicatin by the y wave-number). 6
7 Figure. Cnturs f the critical wave number a in the M, N plane fr v =.0, L =10. T find s 1 it is nly necessary t remve resnant terms frm the right-hand sides f (.11a,b). This is easily dne by multiplying (.11a) by the functin f 1 = e iξ /(a [1 s]) sin mξ which has the frm f the cmplex cnjugate f the nrmal mde f the upper layer, and multiplying (.11b) by f = e i ξ /(a s) sinnξ and integrating ver the width f the channel. After sme algebra this leads t tw equatins relating the amplitudes f the nrmal mdes f the tw layers, viz.: 7
8 and A = A 1 ( ) s s 1 4a (1 s ) + 1 (1 s )α 4a + m γ 1 (1 s ) L (.1a) A 1 = A s 1 4a s (1 s ) + s α (4a + n ) L γ 1 * s (.1b) where the cupling cnstant γ 1 is given by: γ 1 = i mn(k m + k n ) e i(k m+ kn)l ( 1) M +N [ 1 ] [(m n) (k m + k n ) ][ (m + n) (k m + k n ) ] (.13) where: k m = m + a ( ) 1/ k n = n + a ( ) 1/. Nte that the asterisk in (.1b) dentes the cmplex cnjugate f the cupling cefficient γ 1. Eliminating A 1 and A between (.1 a,b) leads t a quadratic equatin fr s s + s B+ C = 0, (1 ) s (1 s ) B = a s s s (1 s) α( a ) ( ) + n α a + m (1 s ) s (1 ) / C = γ s s L 4 4 a s (1 s ) (1 s) sα( a + n ) s (1 s) α( a + m ) (.14 a,b,c) whse slutin determines the stability f the mde whse structure t first-rder is given by (.6 a,b). 8
9 3. Results In terms f ur riginal nn-dimensinal variables the grwth rate f an unstable mde is given by lc i where c i is the imaginary part f the wave s phase speed. We have already nted that t lwest rder in ur expansin in pwers f v the phase speed is real and an imaginary part, if it exists, will depend n the imaginary part f s 1. It then fllws that the grwth rate σ will be σ = lc i = av 3/ imag (s 1 ) (3.1) where the scaled wave number a is itself given by the expansin (.4) s that the grwth rate is a functin f α. In the fllwing figures the grwth rate is given in terms f its riginal variables, i.e. l and in terms f v r its inverse β =1/v. Figure 3. Grwth rate versus wave number fr v =.0, M = 50, N = 0, L = 5. Nte the narrw windw in l fr which σ I > 0. 9
10 Figure 3 shws the grwth rate as a functin f l fr v=.0( β = 50 ) fr the mde crrespnding t M=50 and N =0 fr a channel f width 5 L D. There is a narrw range f wavenumbers in which the grwth rate is psitive and the magnitude f the grwth rate is quite small, fr these values f the rder f 10-4 (in units βl D ). There is a symmetry in M,N s that the mde with M and N reversed has the same grwth rate. Any cmbinatin (M,N) will have sme range ver which its grwth rate will be nn zer s that a plt f grwth rate versus l fr small v wuld shw a bewildering superpsitin f intervals f unstable wave numbers. This is precisely what WP fund numerically but the smallness f the grwth rates precluded a clear picture fr v less than 0.5. Fr example, the grwth rate curve fr M = 5, N = 5 is shwn in Figure 4. The maximum grwth rate is slightly larger than in the previus example. Indeed, as shwn belw, the maximum crrespnds t the diagnal slutin M=N but the differences are nt great enugh t truly privilege that mde. Nte t that the interval f wave number is different than the previus example. As WP fund, these mdes are unstable fr values f l which exceed the classical shrt wave cut-ff f 1/4. Figure 4. As in Figure 3 except fr M =5, N = 5. Nte the change in the l interval f instability. 10
11 An alternative representatin represents the maximum grwth rate fr each M,N in the M,N plane. One can identify the crrespnding y-wavenumber l with the aid f Figure. Figure 5 shws just such a cntur plt. The maximum ccurs very clse t the values chsen fr Figure 4 and the maximum value is abut and ccurs n the M,N diagnal. Nte that n the diagnal s = 0.5, a result that culd be anticipated frm symmetry, s that fr such diagnal mdes the structure f the eigenfunctins is, t lwest rder, the same in each layer. There will be a departure at higher rder and the ff-diagnal unstable mdes, whse grwth rates are cmmensurate with the diagnal mdes will have very different structures in tw layers. Figure 5. Cnturs f grwth rate in the M N plane. Nte that the maximum ccurs alng the line M = N. 4. Discussin The examinatin f the instability f meridinal barclinic shear flws prduces sme nvel instability characteristics when cmpared with the classical prblem f the instability f znal flws. Perhaps nne is mre surprising than the absence f a critical threshld fr instabil- 11
12 ity fr the shear even in thse cases, as studied in this paper, where the gemetry f the flw frces the disturbance t sense the effect f β. In spite f cnsiderable effrt a necessary cnditin fr instability has nt been prven fr the system (.1) and althugh a negative is nt a prf it was thught t suggest that all meridinal flws wuld be unstable. That interesting hypthesis is put n firmer grund by the asympttic result f this paper which shws that very weak meridinal shears can destabilize therwise neutral Rssby nrmal mdes fr the channel. The preexisting mde has its structure slightly altered by the shear allwing the release f the available ptential energy in the shear flw. Since the shears we have cnsidered are weak the crrespnding energy surce is rather feeble and the resulting grwth rates are small. One culd stretch the asympttics t values f v which are nt very small t btain larger grwth rates but that range is already cvered by the detailed numerical analysis f WP which qualitatively agrees with the present results. Frm the pint f view f ceanic applicatins that may be unnecessary. Althugh the gemetry f the channel is simple cmpared t the gemetry f the subtrpical gyre it is still true that instabilities f the flw in the eastern regins f the gyre can be interpreted in terms f the instability f meridinal flw (see fr example Spall, 000) and emphasizes the imprtant rle f znal bundaries (here taken as simple meridians) in affecting the variability f the mid-basin flw. This leads t the interesting pssibility that instabilities f the type described in this paper can extract sufficient energy frm the gyre circulatin t maintain Rssby nrmal mdes against dissipatin as lng as the weak grwth rates pertinent t such mdes exceed the dissipatin rates fr such mdes. Cessi and Primeau (001) and LaCasce (000) have already pinted ut the existence f such special Rssby nrmal mdes with exceedingly lng dissipative times under the influence f scale selective dissipatin. The cupled instability f the Rssby nrmal mdes and the barclinic shear culd prvide a mechanism fr the generatin and maintenance f the nrmal mdes rather that relying n the persistent ringing f the system by repeated external frcing. Acknwledgments. This research was supprted in part by a grant frm the Natinal Science Fundatin, OCE References Cessi, P. and F. Primeau, 001. Dissipative selectin f lw-frequency mdes in a reduced gravity basin. Jurnal f Physical Oceangraphy, 31, LaCasce, J. H., 000. Barclinic Rssby waves in a square basin. Jurnal f Physical Oceangraphy, 30, Pedlsky, J., Gephysical Fluid Dynamics. Springer-Verlag. New Yrk, pp Spall, M., 000. Generatin f strng messcale eddies by weak cean gyres. Jurnal f Marine Research, 58, Walker, A. and J. Pedlsky, 00. On the instability f meridinal barclinic currents. Jurnal f Physical Oceangraphy (t appear). 1
Module 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 informationKinematic transformation of mechanical behavior Neville Hogan
inematic transfrmatin f mechanical behavir Neville Hgan Generalized crdinates are fundamental If we assume that a linkage may accurately be described as a cllectin f linked rigid bdies, their generalized
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 information, which yields. where z1. and z2
The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin
More informationHomology groups of disks with holes
Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.
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 informationChapter 3 Kinematics in Two Dimensions; Vectors
Chapter 3 Kinematics in Tw Dimensins; Vectrs Vectrs and Scalars Additin f Vectrs Graphical Methds (One and Tw- Dimensin) Multiplicatin f a Vectr b a Scalar Subtractin f Vectrs Graphical Methds Adding Vectrs
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 informationPressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects
Pressure And Entrpy Variatins Acrss The Weak Shck Wave Due T Viscsity Effects OSTAFA A. A. AHOUD Department f athematics Faculty f Science Benha University 13518 Benha EGYPT Abstract:-The nnlinear differential
More informationOF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION
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 informationElectric Current and Resistance
Electric Current and Resistance Electric Current Electric current is the rate f flw f charge thrugh sme regin f space The SI unit f current is the ampere (A) 1 A = 1 C / s The symbl fr electric current
More informationAdmissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs
Admissibility Cnditins and Asympttic Behavir f Strngly Regular Graphs VASCO MOÇO MANO Department f Mathematics University f Prt Oprt PORTUGAL vascmcman@gmailcm LUÍS ANTÓNIO DE ALMEIDA VIEIRA Department
More informationLecture 17: Free Energy of Multi-phase Solutions at Equilibrium
Lecture 17: 11.07.05 Free Energy f Multi-phase Slutins at Equilibrium Tday: LAST TIME...2 FREE ENERGY DIAGRAMS OF MULTI-PHASE SOLUTIONS 1...3 The cmmn tangent cnstructin and the lever rule...3 Practical
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 information1 The limitations of Hartree Fock approximation
Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants
More information[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )
(Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well
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 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 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 informationIntroductory Thoughts
Flw Similarity By using the Buckingham pi therem, we have reduced the number f independent variables frm five t tw If we wish t run a series f wind-tunnel tests fr a given bdy at a given angle f attack,
More informationCompressibility and collisional effects on thermal instability of a partially ionized medium
Pram~na, Vl. I0, N. 3, March 978, pp. 267-272, printed in India. Cmpressibility and cllisinal effects n thermal instability f a partially inized medium R C SHARMA and K C SHARMA Department f Mathematics,
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More information13. PO TREATMENT OF DEPT (DISTORTIONLESS ENHANCEMENT POLARIZATION TRANSFER)
94 Prduct Operatr Treatment 3. PO TREATMENT OF DEPT (DISTORTIONLESS ENHANCEMENT POLARIZATION TRANSFER) DEPT is a ne-dimensinal sequence used as a tl fr unambiguus identificatin f the CH, CH, and CH 3 peaks
More information^YawataR&D Laboratory, Nippon Steel Corporation, Tobata, Kitakyushu, Japan
Detectin f fatigue crack initiatin frm a ntch under a randm lad C. Makabe," S. Nishida^C. Urashima,' H. Kaneshir* "Department f Mechanical Systems Engineering, University f the Ryukyus, Nishihara, kinawa,
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 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 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 informationNUMBERS, MATHEMATICS AND EQUATIONS
AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t
More informationCourse Stabilty of Structures
Curse Stabilty f Structures Lecture ntes 2015.03.06 abut 3D beams, sme preliminaries (1:st rder thery) Trsin, 1:st rder thery 3D beams 2:nd rder thery Trsinal buckling Cupled buckling mdes, eamples Numerical
More informationSections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.
Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage
More informationEXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE
EXPERIMENTAL STUD ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE Tmnbu Gt, Masaaki Ohba, Takashi Kurabuchi 2, Tmyuki End 3, shihik Akamine 4, and Tshihir Nnaka 2
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 informationWe can see from the graph above that the intersection is, i.e., [ ).
MTH 111 Cllege Algebra Lecture Ntes July 2, 2014 Functin Arithmetic: With nt t much difficulty, we ntice that inputs f functins are numbers, and utputs f functins are numbers. S whatever we can d with
More informationGAUSS' LAW E. A. surface
Prf. Dr. I. M. A. Nasser GAUSS' LAW 08.11.017 GAUSS' LAW Intrductin: The electric field f a given charge distributin can in principle be calculated using Culmb's law. The examples discussed in electric
More informationLecture 5: Equilibrium and Oscillations
Lecture 5: Equilibrium and Oscillatins Energy and Mtin Last time, we fund that fr a system with energy cnserved, v = ± E U m ( ) ( ) One result we see immediately is that there is n slutin fr velcity if
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 informationIntroduction to Smith Charts
Intrductin t Smith Charts Dr. Russell P. Jedlicka Klipsch Schl f Electrical and Cmputer Engineering New Mexic State University as Cruces, NM 88003 September 2002 EE521 ecture 3 08/22/02 Smith Chart Summary
More informationSupporting information
Electrnic Supplementary Material (ESI) fr Physical Chemistry Chemical Physics This jurnal is The wner Scieties 01 ydrgen perxide electrchemistry n platinum: twards understanding the xygen reductin reactin
More informationFlipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System
Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed
More informationChapter 2 GAUSS LAW Recommended Problems:
Chapter GAUSS LAW Recmmended Prblems: 1,4,5,6,7,9,11,13,15,18,19,1,7,9,31,35,37,39,41,43,45,47,49,51,55,57,61,6,69. LCTRIC FLUX lectric flux is a measure f the number f electric filed lines penetrating
More informationKinetics of Particles. Chapter 3
Kinetics f Particles Chapter 3 1 Kinetics f Particles It is the study f the relatins existing between the frces acting n bdy, the mass f the bdy, and the mtin f the bdy. It is the study f the relatin between
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 informationEquilibrium of Stress
Equilibrium f Stress Cnsider tw perpendicular planes passing thrugh a pint p. The stress cmpnents acting n these planes are as shwn in ig. 3.4.1a. These stresses are usuall shwn tgether acting n a small
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 informationarxiv:hep-ph/ v1 2 Jun 1995
WIS-95//May-PH The rati F n /F p frm the analysis f data using a new scaling variable S. A. Gurvitz arxiv:hep-ph/95063v1 Jun 1995 Department f Particle Physics, Weizmann Institute f Science, Rehvt 76100,
More informationSurface and Contact Stress
Surface and Cntact Stress The cncept f the frce is fundamental t mechanics and many imprtant prblems can be cast in terms f frces nly, fr example the prblems cnsidered in Chapter. Hwever, mre sphisticated
More informationAn Introduction to Complex Numbers - A Complex Solution to a Simple Problem ( If i didn t exist, it would be necessary invent me.
An Intrductin t Cmple Numbers - A Cmple Slutin t a Simple Prblem ( If i didn t eist, it wuld be necessary invent me. ) Our Prblem. The rules fr multiplying real numbers tell us that the prduct f tw negative
More informationMODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b
. REVIEW OF SOME BASIC ALGEBRA MODULE () Slving Equatins Yu shuld be able t slve fr x: a + b = c a d + e x + c and get x = e(ba +) b(c a) d(ba +) c Cmmn mistakes and strategies:. a b + c a b + a c, but
More informationCHAPTER 8b Static Equilibrium Units
CHAPTER 8b Static Equilibrium Units The Cnditins fr Equilibrium Slving Statics Prblems Stability and Balance Elasticity; Stress and Strain The Cnditins fr Equilibrium An bject with frces acting n it, but
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 informationAP Statistics Notes Unit Two: The Normal Distributions
AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).
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 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 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 informationOn Boussinesq's problem
Internatinal Jurnal f Engineering Science 39 (2001) 317±322 www.elsevier.cm/lcate/ijengsci On Bussinesq's prblem A.P.S. Selvadurai * Department f Civil Engineering and Applied Mechanics, McGill University,
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 information( ) (1) ρ c crustal density 2600 kg m -3 ρ w water density 1000 kg m -3. HEAT FLOW PARADOX (Copyright 2001, David T. Sandwell)
1 HEAT FLOW PARADOX (Cpyright 2001, David T. Sandwell) (See Special Issue f J. Gephys. Res., v.85, 1980: A) Turctte, Tag, and Cper, A Steady-State mdel fr the distributin f stress and temperature n the
More informationLecture 6: Phase Space and Damped Oscillations
Lecture 6: Phase Space and Damped Oscillatins Oscillatins in Multiple Dimensins The preius discussin was fine fr scillatin in a single dimensin In general, thugh, we want t deal with the situatin where:
More informationDepartment of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets
Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0
More informationInterference is when two (or more) sets of waves meet and combine to produce a new pattern.
Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme
More informationPhysics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018
Michael Faraday lived in the Lndn area frm 1791 t 1867. He was 29 years ld when Hand Oersted, in 1820, accidentally discvered that electric current creates magnetic field. Thrugh empirical bservatin and
More informationEngineering Approach to Modelling Metal THz Structures
Terahertz Science and Technlgy, ISSN 1941-7411 Vl.4, N.1, March 11 Invited Paper ngineering Apprach t Mdelling Metal THz Structures Stepan Lucyszyn * and Yun Zhu Department f, Imperial Cllege Lndn, xhibitin
More informationPreparation work for A2 Mathematics [2018]
Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time
More informationTHE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS
THE TOPOLOGY OF SURFACE SKIN FRICTION AND VORTICITY FIELDS IN WALL-BOUNDED FLOWS M.S. Chng Department f Mechanical Engineering The University f Melburne Victria 3010 AUSTRALIA min@unimelb.edu.au J.P. Mnty
More informationFundamental Concepts in Structural Plasticity
Lecture Fundamental Cncepts in Structural Plasticit Prblem -: Stress ield cnditin Cnsider the plane stress ield cnditin in the principal crdinate sstem, a) Calculate the maximum difference between the
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 informationAircraft Performance - Drag
Aircraft Perfrmance - Drag Classificatin f Drag Ntes: Drag Frce and Drag Cefficient Drag is the enemy f flight and its cst. One f the primary functins f aerdynamicists and aircraft designers is t reduce
More informationENGI 4430 Parametric Vector Functions Page 2-01
ENGI 4430 Parametric Vectr Functins Page -01. Parametric Vectr Functins (cntinued) Any nn-zer vectr r can be decmpsed int its magnitude r and its directin: r rrˆ, where r r 0 Tangent Vectr: dx dy dz dr
More informationDispersion Ref Feynman Vol-I, Ch-31
Dispersin Ref Feynman Vl-I, Ch-31 n () = 1 + q N q /m 2 2 2 0 i ( b/m) We have learned that the index f refractin is nt just a simple number, but a quantity that varies with the frequency f the light.
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 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 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 informationMedium Scale Integrated (MSI) devices [Sections 2.9 and 2.10]
EECS 270, Winter 2017, Lecture 3 Page 1 f 6 Medium Scale Integrated (MSI) devices [Sectins 2.9 and 2.10] As we ve seen, it s smetimes nt reasnable t d all the design wrk at the gate-level smetimes we just
More informationSuggested reading: Lackmann (2011), Sections
QG Thery and Applicatins: Apprximatins and Equatins Atms 5110 Synptic Dynamic Meterlgy I Instructr: Jim Steenburgh jim.steenburgh@utah.edu 801-581-8727 Suite 480/Office 488 INSCC Suggested reading: Lackmann
More informationB. Definition of an exponential
Expnents and Lgarithms Chapter IV - Expnents and Lgarithms A. Intrductin Starting with additin and defining the ntatins fr subtractin, multiplicatin and divisin, we discvered negative numbers and fractins.
More informationPreparation work for A2 Mathematics [2017]
Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will
More informationGENERAL FORMULAS FOR FLAT-TOPPED WAVEFORMS. J.e. Sprott. Plasma Studies. University of Wisconsin
GENERAL FORMULAS FOR FLAT-TOPPED WAVEFORMS J.e. Sprtt PLP 924 September 1984 Plasma Studies University f Wiscnsin These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated.
More informationQ1. A string of length L is fixed at both ends. Which one of the following is NOT a possible wavelength for standing waves on this string?
Term: 111 Thursday, January 05, 2012 Page: 1 Q1. A string f length L is fixed at bth ends. Which ne f the fllwing is NOT a pssible wavelength fr standing waves n this string? Q2. λ n = 2L n = A) 4L B)
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 informationMATHEMATICS SYLLABUS SECONDARY 5th YEAR
Eurpean Schls Office f the Secretary-General Pedaggical Develpment Unit Ref. : 011-01-D-8-en- Orig. : EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 6 perid/week curse APPROVED BY THE JOINT TEACHING COMMITTEE
More informationI. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is
Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,
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 informationFree Vibrations of Catenary Risers with Internal Fluid
Prceeding Series f the Brazilian Sciety f Applied and Cmputatinal Mathematics, Vl. 4, N. 1, 216. Trabalh apresentad n DINCON, Natal - RN, 215. Prceeding Series f the Brazilian Sciety f Cmputatinal and
More informationBASD HIGH SCHOOL FORMAL LAB REPORT
BASD HIGH SCHOOL FORMAL LAB REPORT *WARNING: After an explanatin f what t include in each sectin, there is an example f hw the sectin might lk using a sample experiment Keep in mind, the sample lab used
More informationChE 471: LECTURE 4 Fall 2003
ChE 47: LECTURE 4 Fall 003 IDEL RECTORS One f the key gals f chemical reactin engineering is t quantify the relatinship between prductin rate, reactr size, reactin kinetics and selected perating cnditins.
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 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 informationChapter 9 Vector Differential Calculus, Grad, Div, Curl
Chapter 9 Vectr Differential Calculus, Grad, Div, Curl 9.1 Vectrs in 2-Space and 3-Space 9.2 Inner Prduct (Dt Prduct) 9.3 Vectr Prduct (Crss Prduct, Outer Prduct) 9.4 Vectr and Scalar Functins and Fields
More informationIntroduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem
A Generalized apprach fr cmputing the trajectries assciated with the Newtnian N Bdy Prblem AbuBar Mehmd, Syed Umer Abbas Shah and Ghulam Shabbir Faculty f Engineering Sciences, GIK Institute f Engineering
More informationAdvanced Heat and Mass Transfer by Amir Faghri, Yuwen Zhang, and John R. Howell
6.5 Natural Cnvectin in Enclsures Enclsures are finite spaces bunded by walls and filled with fluid. Natural cnvectin in enclsures, als knwn as internal cnvectin, takes place in rms and buildings, furnaces,
More informationDead-beat controller design
J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable
More informationCHAPTER Read Chapter 17, sections 1,2,3. End of Chapter problems: 25
CHAPTER 17 1. Read Chapter 17, sectins 1,2,3. End f Chapter prblems: 25 2. Suppse yu are playing a game that uses tw dice. If yu cunt the dts n the dice, yu culd have anywhere frm 2 t 12. The ways f prducing
More informationSimple Coupled Midlatitude Climate Models
016 JOURNAL OF PHYSICAL OCEANOGRAPHY VOLUME 9 Simple Cupled Midlatitude Climate Mdels LYNNE D. TALLEY Scripps Institutin f Oceangraphy, University f Califrnia, San Dieg, La Jlla, Califrnia (Manuscript
More informationDr M. BROUARD. 5. Thermodynamic formulation of Transition State Theory Entropy of activation. Thermochemical kinetics. CHEMICAL REACTION RATES
CHEMICAL REACTION RATES Dr M. BROUARD Trinity Term 2003 A. Bimlecular Reactins 5 Lectures 1. Intrductin Simple cllisin thery. Ptential energy curves and surfaces. The reactin crdinate and barriers t reactin.
More informationPhys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1
Crdinatr: Nasser Wednesday, January 4, 007 Page: 1 Q1. Tw transmitters, S 1 and S shwn in the figure, emit identical sund waves f wavelength λ. The transmitters are separated by a distance λ /. Cnsider
More informationL a) Calculate the maximum allowable midspan deflection (w o ) critical under which the beam will slide off its support.
ecture 6 Mderately arge Deflectin Thery f Beams Prblem 6-1: Part A: The department f Highways and Public Wrks f the state f Califrnia is in the prcess f imprving the design f bridge verpasses t meet earthquake
More informationA Correlation of. to the. South Carolina Academic Standards for Mathematics Precalculus
A Crrelatin f Suth Carlina Academic Standards fr Mathematics Precalculus INTRODUCTION This dcument demnstrates hw Precalculus (Blitzer), 4 th Editin 010, meets the indicatrs f the. Crrelatin page references
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 informationHow do scientists measure trees? What is DBH?
Hw d scientists measure trees? What is DBH? Purpse Students develp an understanding f tree size and hw scientists measure trees. Students bserve and measure tree ckies and explre the relatinship between
More information