2.3. PBL Equations for Mean Flow and Their Applications
|
|
- Loreen Flynn
- 5 years ago
- Views:
Transcription
1 .3. PBL Equations for Mean Flow and Their Applications Read Holton Section 5.3!.3.1. The PBL Momentum Equations We have derived the Reynolds averaed equations in the previous section, and they describe the mean flow, takin into account the effect of turbulence motion. Generally, inside the boundary layer, vertical turbulent flux diverence is much bier than horiontal one, because vertical shear/radient is much bier. We often consider the so-called horiontally homoeneous turbulence for which the horiontal turbulent flux diverence is nelected. We therefore have, from (.5), du 1 p wu ' ' = + fv, (.7a) dt ρ x 0 dv 1 p wv ' ' = fu. (.7b) dt ρ y 0 Consider further that the mean variables are for synoptic-scale horiontal flows, the acceleration terms are much smaller than the PGF and Coriolis terms, they can be nelected. 1 p 1 p Further, because fv =, and fu =, we replace the PGF ρ0 x ρ0 y terms usin u and v : wu ' ' f( v v ) = 0, (.8a) wv ' ' f( u u ) = 0. (.8a) 1
2 They represent the balance amon Coriolis force, pressure radient force and the vertical momentum flux diverence (a net stress applied to the air parcel). We will use the above two equations to solve for boundary-layer wind in two different models..3.. The Mixed Layer Model (Section of Holton) As we saw earlier, in the well-mixed boundary layer, the mean wind and mean potential temperature are nearly constant with heiht. If we look at the flux profiles, they decrease linearly with heiht, and to ero at the top of boundary layer (they often o to neative values before becomin ero) where turbulence becomes weak. Accordin to observations, the momentum fluxes at the surface can be accurately represented by the bulk aerodynamic formula (or dra law): ( wu ' ') sfc = Cd V u, (.9) ( wv ' ') = C V v, sfc d where C d = non-dimensional dra coefficient and V = u + v is the wind speed at the surface, or more precisely, at the anemometer heiht (10 m). The above formulae are very important and are used commonly in numerical weather prediction models to model the surface fluxes. Similar formulae are available for heat and moisture. The dra coefficient C d depends on the rouhness of the surface and the stability of the surface layer air. Over ocean, it is on the order of 10-3 and over land on the order of Now that we know the fluxes at the surface (~=0) and at the top of the boundary layer (=h) and know that the fluxes are linear function of heiht. We also assumin that the mean wind is constant with heiht, we can easily interate Equations (.8a,b) vertically, i.e.,
3 h h wu ' ' f ( v v ) d = d 0 0 h h wu ' ' f ( v v ) d = d = w' u ' w' u ' = Cd V u 0 0 h 0 Similarly, u v f ( v v ) h =. (.30a) C V d f ( u u ) h =. (.30b) C V d Without loss of enerality, we can choose our coordinate system so that the x-axis is parallel to u so that v = 0. Therefore, u = u K V v v = K V s where K C /( fh) and has units of s/m. s d s (.31a) (.31b) Since u is the wind speed (positive), we find that inside the mixed layer, the wind speed is subeostrophic (we now see why), i.e., less than the eostrophic wind u. Also we have a component of motion directin towards the low-pressure side (to the left of eostrophic wind in northern hemisphere and the riht in the S.H.) that is proportional to K s which is ero without surface dra. The v component is completely aeostrophic (perpendicular to the isobars), and arises due to the dra/fiction/momentum flux. The stroner is the surface dra, the larer is the anle of the wind to the eostrophic wind / pressure contours. 3
4 v u u Example, u = 10 m/s K s = / ( ) = 0.05 s m -1 u =8.8 m/s, v =3.77 m/s, V = 9.10m/s at all heiht within this well mixed layer. All of this is consistent with what you have learned before about the three-force balance. In a eostrophic flow, the PGF and the Coriolis force are in balance. When friction is present, friction acts in the opposite direction as the velocity vector and is therefore also perpendicular to the Coriolis force. The balance between the three forces are illustrated as follows: P- PGF P 0 φ V F r Cor P + As can be seen above, the wind vector is no loner parallel to the pressure contours, the flow is crossin the isobars to the lower pressure side! In vector form, the balance is expressed as 4
5 ˆ 1 C d fk V = p V sfc V. (.3) ρ h 0 The last term is sometimes called friction but keep in mind that this is due to turbulent mixin, not molecular viscosity, the latter is much smaller K-Theory: First-order Turbulence Closure (Section 5.3. of Holton) When the boundary layer is neutral or stable, the constant wind speed assumption used in the previous slab model is no loner ood (because of the absence of stron vertical mixin). To solve the equations for the mean variables, we have to express the turbulent fluxes in terms of the mean variables, to close the equation set (so that we have the same number of unknowns as the number of equations). This is called the closure problem. One of the commonly used methods is to express the fluxes in terms of the radient of the mean variables, in a way analoous to the molecular diffusion. The flux is assumed to be proportional to the local radient of the mean: u uw ' ' = K m, v vw ' ' = K m, (.33a) (.33b) and similarly to the potential temperature: θ θ ' w' = Kh. (.33c) 5
6 We call K m the eddy viscosity (has units of m /s) and K h the eddy diffusivity of heat. Both are positive. This closure scheme is often called the K-theory. This closure says that flux is always down-radient. With the closure iven in (.33), we can write equations in (.7) as du 1 p fv K u = + + m dt ρ0 x, (.34a) dv 1 p fu K v = + m dt ρ y. (.34b) 0 We see that the perturbation variables no loner appear in the equations. Given p we can solve for uandv. For the potential temperature, the equation is d θ θ = Kh + Qθ dt. (.34c) We will use this set of equations to obtain the Ekman layer solution in the next section. Note also that the K-theory actually has a number of limitations. Firstly K should be flow-dependent and it is often not easy to determine its value. Secondly, often the eddies that cause most of the mixin have sies that are comparable to the PBL depth and are therefore not so much related to local radient therefore the closure assumption does not work well. For this reason, certain non-local closure schemes have been developed. 6
7 .3.4. The Ekman Layer (Section of Holton) The Ekman layer is the part of the atmosphere that extends from the anemometer level up to the heiht where the radient wind balance holds (PGF + COR + CENTRIF). In this layer, we make the followin assumptions within the context of K-theory: a) The force balance is amon PGF, Coriolis and turbulence. Centrifual force is nelected. b) The flow is steady-state, i.e., / t = 0. c) The atmosphere is horiontally homoeneous, i.e., ()/ x+ ()/ y = 0. this implies that V is independent of heiht, via the thermal wind law. d) Lapse rate is neutral, i.e., θ / = 0. So, the PBL is well mixed. e) K m is constant f) Earth is flat with w =0 at round. ) Flow is incompressible. With these assumptions, our mean flow momentum equations in (.34) become u + f( v v ) 0 = v f( u u ) 0 = (.35a) (.35b) The above two equations are coupled. One way to simplify the solution procedure is to combine the two equations into one, by usin complex variables. Define complex velocities w= u + iv and w = u + iv. (.35a) + i (.35b) ( u + iv) m K + f( v ui v + u i) = 0 7
8 ( u + iv) m K + fi( iv u + iv + u ) = 0 w fi( w w ) 0 = (.36) Since we assumed that w is constant with heiht, (.36) can be rewritten as ( w w ) if ( w w ) 0 = K m (.37) which is a linear, nd-order ODE for ( w w ) with constant coefficients a basic ODE problem. The characteristic equation for (.37) is if λ = 0 if therefore λ =± where λ's are the eienvalues of the PDE. K m The solutions to (.37) are therefore (verify for yourself): if if w w = Ae + Be. (.38) Here A and B are arbitrary constants to be determined by the boundary conditions. The boundary conditions are w =0 at = 0 0 w = A+ B w= w as 0 = Ae A = 0. if Therefore B = w. Solution (.38) becomes 8
9 w= w w e if K m Without loss of enerality, we choose our x coordinate axis to be parallel to the eostrophic wind, so that w = u, v =0. w u u e if K m = (.39) Let's see what i is equal to. i= 1, 1+ i= 1+ 1, (1 + i) = (1 + 1) 1+ i (1 + i) = = i = i Our solution is therefore f K (1 i) m w= u (1 e + ) γ γ u + iv = u [1 e cos( γ ) + ie sin( γ )] γ u = u [1 e cos( γ )] (.40a) γ v = u e sin( γ ) (.40b) When plotted the solution on a hodoraph, which is a plot of the wind vector on a polar coordinate diaram, we obtain what is called the Ekman Spiral solution as shown below: 9
10 We see that with this solution, there is an aeostrophic wind component pointin to the left of eostrophic wind, i.e., towards the low-pressure side, this arees with the results of the earlier mixed-layer model, and the solution represents a three-way balance, i.e., a balance between the PGF, Coriolis and frictional forces. Compared to the mixed layer model, we don t assume that the boundary wind is constant with heiht, and we actually obtained the solution for the vertical wind profile. The idealied Ekman spiral wind profile is rarely observed in the atmosphere, however, because of other effects that are not considered in the idealied model. In the surface layer, the equations used to obtain this solution are not valid. One way to take it into account is to match the Ekman layer solution with that in the surface layer, doin so ives us the so called modified Ekman layer solution is that often a better match of the observations. The followin fiure ives such an example. 10
11 .3.5. Secondary Circulation and Spin-Down (Section 5.4 of Holton) From both the Ekman layer solution and the mix-layer solution obtained in section.3. (Equation.31), we know that the friction reduces the boundary layer wind speed below eostrophic wind, and causes it to cross the isobars from hih towards low pressure. In a synoptic situation where the isobars are curved, such as in low and hih pressure systems, the cross-isobaric component of flow near the surface causes risin air in low-pressure systems and descendin air in hihs. The processin of inducin vertical motions by boundary layer friction is called Ekman pumpin. 11
12 Usin Ekman spiral solution, we can quantitatively estimate the vertical motion at the top of boundary layer. Aain assume v = 0. From mass continuity equation and the Ekman spiral solution, the vertical velocity at the top of the boundary layer is D w d = 0 0 D v d y D v D u γ γ 0 (.41) 0 w( D) = d = e sin( ) d y y where D is the boundary layer depth. The mass flux across the isobars in a column of unit width spannin the depth of the boundary layer is D D (.4) γ M = ρ vd = ρ u e sin( γ) d where the depth D can be defined as D = π / γ, the level at which the wind (as iven by the spiral solution) becomes parallel to the eostrophic wind. Comparin (.4) to (.41), we see that M ρ0 wd ( ) = (.43) y therefore the vertical mass flux out of the PBL is equal to the converence of the cross-isobaric mass transport in the layer. u Further, = ξ, the eostrophic vorticity, and substitutin the Ekman y spiral solution into (.41) (notin the eostrophic wind does not chane with heiht), we obtain 1
13 1/ f wd ( ) = ξ f f (.44) hence the vertical mass flux is proportional to the eostrophic vorticity and the square root of the boundary layer eddy viscocity. In the northern hemisphere (f>0), w at Z=D is positive inside a cyclone ( ξ > 0) and neative inside an anticyclone. This w induces a secondary (vertical) circulation as shown below. Above the PBL inside a cyclone, there exists diverence, accordin to anular momentum conservation, the horiontal rotation rate (vorticity) will decrease with time, hence the spindown of the cyclone. Therefore, the spin down is not by simply surface friction and vertical momentum exchane, but throuh a most effective way that involves the vertical (secondary) circulation! A quantitative discussion of this usin a vorticity equation is iven in Holton section 5.4. Read it yourself! 13
14 Tea Cup Example. After the tea in a cup is stirred, the bottom friction slows down the spinnin tea in the cup near the bottom surface, breakin the balance between the PGF and centrifual force. The PGF cause flow towards the center of the cup near the bottom, creatin vertical motion alon the central axis and induces compensatin flow above this layer away from the central axis. To conserve anular momentum, the spin has to slow down. This is analoous to the spin-down of a cyclone. 14
Chapter 2. Turbulence and the Planetary Boundary Layer
Chapter 2. Turbulence and the Planetary Boundary Layer In the chapter we will first have a qualitative overview of the PBL then learn the concept of Reynolds averain and derive the Reynolds averaed equations.
More informationESCI 343 Atmospheric Dynamics II Lesson 1 Ageostrophic Wind
ESCI 343 Atmospheric Dynamics II Lesson 1 Aeostrophic Wind References: An Introduction to Dynamic Meteoroloy (3 rd edition), J.R. Holton THE QG MOMENTUM EQUATIONS The QG momentum equations are derived
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Energy Cycle
Course.8, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Enery Cycle Enery Forms: As we saw in our discussion of the heat budet, the enery content of the atmosphere per
More informationDevelopment Theory. Chapter 10
Chapter 1 Development Theory Development Theory In this section I will discuss: - further aspects of the structure and dynamics of synoptic-scale disturbances, and - derive Sutcliffe's development theory,
More informationEffective Depth of Ekman Layer.
5.5: Ekman Pumping Effective Depth of Ekman Layer. 2 Effective Depth of Ekman Layer. Defining γ = f/2k, we derived the solution u = u g (1 e γz cos γz) v = u g e γz sin γz corresponding to the Ekman spiral.
More information5 Shallow water Q-G theory.
5 Shallow water Q-G theory. So far we have discussed the fact that lare scale motions in the extra-tropical atmosphere are close to eostrophic balance i.e. the Rossby number is small. We have examined
More informationReynolds Averaging. We separate the dynamical fields into slowly varying mean fields and rapidly varying turbulent components.
Reynolds Averaging Reynolds Averaging We separate the dynamical fields into sloly varying mean fields and rapidly varying turbulent components. Reynolds Averaging We separate the dynamical fields into
More informationGeopotential tendency and vertical motion
Geopotential tendency and vertical motion Recall PV inversion Knowin the PV, we can estimate everythin else! (Temperature, wind, eopotential ) In QG, since the flow is eostrophic, we can obtain the wind
More informationBalanced Flow Geostrophic, Inertial, Gradient, and Cyclostrophic Flow
Balanced Flow Geostrophic, Inertial, Gradient, and Cyclostrophic Flow The types of atmospheric flows describe here have the following characteristics: 1) Steady state (meaning that the flows do not change
More informationMotion in Two Dimensions Sections Covered in the Text: Chapters 6 & 7, except 7.5 & 7.6
Motion in Two Dimensions Sections Covered in the Tet: Chapters 6 & 7, ecept 7.5 & 7.6 It is time to etend the definitions we developed in Note 03 to describe motion in 2D space. In doin so we shall find
More informationThe broadest definition of a low-level jet (LLJ) is simply any lower-tropospheric maximum in the vertical profile of the horizontal winds.
2.4.3. Low-level (especially nocturnal) Jet The broadest definition of a low-level jet (LLJ) is simply any lower-tropospheric maximum in the vertical profile of the horizontal winds. A LLJ can occur under
More informationESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr. DeCaria
ESCI 485 Air/Sea Interaction Lesson 1 Stresses and Fluxes Dr DeCaria References: An Introduction to Dynamic Meteorology, Holton MOMENTUM EQUATIONS The momentum equations governing the ocean or atmosphere
More informationThe dynamics of high and low pressure systems
The dynamics of high and low pressure systems Newton s second law for a parcel of air in an inertial coordinate system (a coordinate system in which the coordinate axes do not change direction and are
More informationANALYZE In all three cases (a) (c), the reading on the scale is. w = mg = (11.0 kg) (9.8 m/s 2 ) = 108 N.
Chapter 5 1. We are only concerned with horizontal forces in this problem (ravity plays no direct role). We take East as the +x direction and North as +y. This calculation is efficiently implemented on
More informationESS Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1
ESS5203.03 - Turbulence and Diffusion in the Atmospheric Boundary-Layer : Winter 2017: Notes 1 Text: J.R.Garratt, The Atmospheric Boundary Layer, 1994. Cambridge Also some material from J.C. Kaimal and
More informationA Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere
A Mathematical Model for the Fire-extinuishin Rocket Fliht in a Turbulent Atmosphere CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tiroviste, Km 8 ROMANIA crismihailescu@yahoo.com http://www.elmec.ro
More informationThe homework problem (pucks_on_ice) illustrates several important points:
Ekman laers, friction & eostrohic flow The homework roblem (ucks_on_ice) illustrates several imortant oints: 1. Particles move erendicular to the alied force, to the riht in the northern hemishere, to
More informationAorangi (Mt. Cook) we typically view the state of the troposphere on "isobaric surfaces" (surfaces of constant pressure) rather
Ch3 (ctd.) EAS270_Ch3_BehaviourAtmos_B.odp JDW, EAS UAlberta, last mod. 19 Sept. 2016 Atmospheric "behaviour" Aorani (Mt. Cook) "ride" we typically view the state of the troposphere on "isobaric surfaces"
More informationThe atmospheric boundary layer: Where the atmosphere meets the surface. The atmospheric boundary layer:
The atmospheric boundary layer: Utrecht Summer School on Physics of the Climate System Carleen Tijm-Reijmer IMAU The atmospheric boundary layer: Where the atmosphere meets the surface Photo: Mark Wolvenne:
More informationOCN/ATM/ESS 587. The wind-driven ocean circulation. Friction and stress. The Ekman layer, top and bottom. Ekman pumping, Ekman suction
OCN/ATM/ESS 587 The wind-driven ocean circulation. Friction and stress The Ekman layer, top and bottom Ekman pumping, Ekman suction Westward intensification The wind-driven ocean. The major ocean gyres
More information(! g. + f ) Dt. $% (" g. + f ) + f
The QG Vorticity equation The complete derivation of the QG vorticity equation can be found in Chapter 6.2.. Please read this section and keep in mind that extra approximations have been made besides those
More informationDust devils, water spouts, tornados
Balanced flow Things we know Primitive equations are very comprehensive, but there may be a number of vast simplifications that may be relevant (e.g., geostrophic balance). Seems that there are things
More informationMidterm Feb. 17, 2009 Physics 110B Secret No.=
Midterm Feb. 17, 29 Physics 11B Secret No.= PROBLEM (1) (4 points) The radient operator = x i ê i transforms like a vector. Use ɛ ijk to prove that if B( r) = A( r), then B( r) =. B i = x i x i = x j =
More information(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2
11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity
More informationWire antenna model of the vertical grounding electrode
Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna
More informationESCI 344 Tropical Meteorology Lesson 12 Tropical Cyclones: Motion and Analysis
ESCI 344 Tropical Meteoroloy Lesson 12 Tropical Cyclones: Motion and Analysis References: A Global View of Tropical Cyclones, Elsberry (ed.) The Hurricane, Pielke Forecasters Guide to Tropical Meteoroloy,
More informationConservation of absolute vorticity. MET 171A: Barotropic Midlatitude Waves. Draw a picture of planetary vorticity
Conservation of absolute vorticity : Barotropic Midlatitude Waves Recall the important terms in the vorticity equation in the middle troposphere, near the level of non-diverence Lecture Outline 1. Conservation
More informationExam 2A Solution. 1. A baseball is thrown vertically upward and feels no air resistance. As it is rising
Exam 2A Solution 1. A baseball is thrown vertically upward and feels no air resistance. As it is risin Solution: Possible answers: A) both its momentum and its mechanical enery are conserved - incorrect.
More informationFundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationEART164: PLANETARY ATMOSPHERES
EART164: PLANETARY ATMOSPHERES Francis Nimmo Last Week Radiative Transfer Black body radiation, Planck function, Wien s law Absorption, emission, opacity, optical depth Intensity, flux Radiative diffusion,
More informationChapter K. Oscillatory Motion. Blinn College - Physics Terry Honan. Interactive Figure
K. - Simple Harmonic Motion Chapter K Oscillatory Motion Blinn Collee - Physics 2425 - Terry Honan The Mass-Sprin System Interactive Fiure Consider a mass slidin without friction on a horizontal surface.
More informationProblem #1: The Gradient Wind in Natural Coordinates (Due Friday, Feb. 28; 20 pts total)
METR 50: Atmospheric Dynamics II Dr. Dave Dempsey Spring 014 Problem #1: The Gradient Wind in Natural Coordinates (Due Friday, Feb. 8; 0 pts total) In natural (s,n,p) coordinates, the synoptic-scaled,
More information11 Free vibrations: one degree of freedom
11 Free vibrations: one deree of freedom 11.1 A uniform riid disk of radius r and mass m rolls without slippin inside a circular track of radius R, as shown in the fiure. The centroidal moment of inertia
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 informationProblem Set 5 Solutions
Problem Set 5 Solutions Dorian Abbot APM47 0/8/04. (a) We consider π θ π. The pendulum pointin down corresponds to θ=0 and the pendulum pointin up corresponds to θ π. Define ν θ. The system can be rewritten:
More information1 Introduction to Governing Equations 2 1a Methodology... 2
Contents 1 Introduction to Governing Equations 2 1a Methodology............................ 2 2 Equation of State 2 2a Mean and Turbulent Parts...................... 3 2b Reynolds Averaging.........................
More informationChapter 1. Introduction
Chapter 1. Introduction In this class, we will examine atmospheric phenomena that occurs at the mesoscale, including some boundary layer processes, convective storms, and hurricanes. We will emphasize
More informationPhysics 111 P 2 A = P 1. A + mg = P 1. A + ρ( AΔh)g. Wednesday, 8-9 pm in NSC 118/119 Sunday, 6:30-8 pm in CCLIR 468.
ics Announcements day, ember 11, 011 C5: Fluids Pascal s Principle Archimede s Principle Fluid Flows Continuity Equation Bernoulli s Equation Toricelli s Theorem Wednesday, 8-9 pm in NSC 118/119 Sunday,
More informationREVIEW: Going from ONE to TWO Dimensions with Kinematics. Review of one dimension, constant acceleration kinematics. v x (t) = v x0 + a x t
Lecture 5: Projectile motion, uniform circular motion 1 REVIEW: Goin from ONE to TWO Dimensions with Kinematics In Lecture 2, we studied the motion of a particle in just one dimension. The concepts of
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 informationAltitude measurement for model rocketry
Altitude measurement for model rocketry David A. Cauhey Sibley School of Mechanical Aerospace Enineerin, Cornell University, Ithaca, New York 14853 I. INTRODUCTION In his book, Rocket Boys, 1 Homer Hickam
More informationQuasi-geostrophic motion
Capter 8 Quasi-eostropic motion Scale analysis for synoptic-scale motions Simplification of te basic equations can be obtained for synoptic scale motions. Consider te Boussinesq system ρ is assumed to
More informationg (z) = 1 (1 + z/a) = 1 1 ( km/10 4 km) 2
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More informationOcean dynamics: the wind-driven circulation
Ocean dynamics: the wind-driven circulation Weston Anderson March 13, 2017 Contents 1 Introduction 1 2 The wind driven circulation (Ekman Transport) 3 3 Sverdrup flow 5 4 Western boundary currents (western
More informationEnergizing Math with Engineering Applications
Enerizin Math with Enineerin Applications Understandin the Math behind Launchin a Straw-Rocket throuh the use of Simulations. Activity created by Ira Rosenthal (rosenthi@palmbeachstate.edu) as part of
More informationONLINE: MATHEMATICS EXTENSION 2 Topic 6 MECHANICS 6.3 HARMONIC MOTION
ONINE: MATHEMATICS EXTENSION Topic 6 MECHANICS 6.3 HARMONIC MOTION Vibrations or oscillations are motions that repeated more or less reularly in time. The topic is very broad and diverse and covers phenomena
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science
Massachusetts nstitute of Technoloy Department of Electrical Enineerin and Computer Science 6.685 Electric Machines Class Notes 2 Manetic Circuit Basics September 5, 2005 c 2003 James L. Kirtley Jr. 1
More information1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance
CHAPTER 7 PROJECTILES 7 No Air Resistance We suppose that a particle is projected from a point O at the oriin of a coordinate system, the y-axis bein vertical and the x-axis directed alon the round The
More information2.2 Differentiation and Integration of Vector-Valued Functions
.. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin
More informationg (z) = 1 (1 + z/a) = 1
1.4.2 Gravitational Force g is the gravitational force. It always points towards the center of mass, and it is proportional to the inverse square of the distance above the center of mass: g (z) = GM (a
More informationConvective Fluxes: Sensible and Latent Heat Convective Fluxes Convective fluxes require Vertical gradient of temperature / water AND Turbulence ( mixing ) Vertical gradient, but no turbulence: only very
More informationLecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2)
Lecture 3. Turbulent fluxes and TKE budgets (Garratt, Ch 2) The ABL, though turbulent, is not homogeneous, and a critical role of turbulence is transport and mixing of air properties, especially in the
More informationPneumatic Conveying in Horizontal Pipes: Eulerian Modeling and Pressure Drop Characteristics
American Journal of Mathematical and Computational Sciences 08; (): 0-6 http://www.aascit.or/journal/ajmcs Pneumatic Conveyin in Horizontal Pipes: Eulerian Modelin and Pressure Drop Characteristics Pandaba
More informationGeneral Comment on Lab Reports: v. good + corresponds to a lab report that: has structure (Intro., Method, Results, Discussion, an Abstract would be
General Comment on Lab Reports: v. good + corresponds to a lab report that: has structure (Intro., Method, Results, Discussion, an Abstract would be a bonus) is well written (take your time to edit) shows
More informationOn side wall labeled A: we can express the pressure in a Taylor s series expansion: x 2. + higher order terms,
Chapter 1 Notes A Note About Coordinates We nearly always use a coordinate system in this class where the vertical, ˆk, is normal to the Earth s surface and the x-direction, î, points to the east and the
More information( ) = 1005 J kg 1 K 1 ;
Problem Set 3 1. A parcel of water is added to the ocean surface that is denser (heavier) than any of the waters in the ocean. Suppose the parcel sinks to the ocean bottom; estimate the change in temperature
More informationPhysics 18 Spring 2011 Homework 2 - Solutions Wednesday January 26, 2011
Physics 18 Sprin 011 Homework - s Wednesday January 6, 011 Make sure your name is on your homework, and please box your final answer. Because we will be ivin partial credit, be sure to attempt all the
More informationExperiment 1: Simple Pendulum
COMSATS Institute of Information Technoloy, Islamabad Campus PHY-108 : Physics Lab 1 (Mechanics of Particles) Experiment 1: Simple Pendulum A simple pendulum consists of a small object (known as the bob)
More information[ ][ ] mg = R GM g = R GM = g R 2. Definition of G: From (1) we have, 2 NEWTON S LAW OF GRAVITATION:
NEWTON S LAW OF GAVITATION: Statement: Every particle in the universe attracts every other particle with a force which is directly proportional to the product of the masses and inversely proportional to
More informationPhysics 20 Homework 1 SIMS 2016
Physics 20 Homework 1 SIMS 2016 Due: Wednesday, Auust 17 th Problem 1 The idea of this problem is to et some practice in approachin a situation where you miht not initially know how to proceed, and need
More informationρ x + fv f 'w + F x ρ y fu + F y Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ Dv Dt + u2 tanφ + vw a a = 1 p Dw Dt u2 + v 2
Fundamental Equation in z coordinate p = ρrt or pα = RT Du uv tanφ + uw Dt a a = 1 p ρ x + fv f 'w + F x Dv Dt + u2 tanφ + vw a a = 1 p ρ y fu + F y Dw Dt u2 + v 2 = 1 p a ρ z g + f 'u + F z Dρ Dt + ρ
More informationFrictional Damping of Baroclinic Waves
Frictional Damping of Baroclinic Waves HHH, 26th May 2006 Bob Plant With thanks to Stephen Belcher, Brian Hoskins, Dan Adamson Motivation Control simulation, T+60 Simulation with no boundary layer turbulence,
More information(a) Find the function that describes the fraction of light bulbs failing by time t e (0.1)x dx = [ e (0.1)x ] t 0 = 1 e (0.1)t.
1 M 13-Lecture March 8, 216 Contents: 1) Differential Equations 2) Unlimited Population Growth 3) Terminal velocity and stea states Voluntary Quiz: The probability density function of a liht bulb failin
More informationHouze sections 7.4, 8.3, 8.5, Refer back to equations in Section 2.3 when necessary.
Thunderstorm Dynamics Houze sections 7.4, 8.3, 8.5, Refer back to equations in Section.3 when necessary. Bluestein Vol. II section 3.4.6. Review article "Dynamics of Tornadic Thunderstorms" by Klemp handout.
More informationSynoptic Meteorology I: Other Force Balances
Synoptic Meteorology I: Other Force Balances For Further Reading Section.1.3 of Mid-Latitude Atmospheric Dynamics by J. Martin provides a discussion of the frictional force and considerations related to
More informationHurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.
Hurricanes: Observations and Dynamics Houze Section 10.1. Holton Section 9.7. Emanuel, K. A., 1988: Toward a general theory of hurricanes. American Scientist, 76, 371-379 (web link). http://ww2010.atmos.uiuc.edu/(gh)/guides/mtr/hurr/home.rxml
More informationCirculation and Vorticity
Circulation and Vorticity Example: Rotation in the atmosphere water vapor satellite animation Circulation a macroscopic measure of rotation for a finite area of a fluid Vorticity a microscopic measure
More information7 Balanced Motion. 7.1 Return of the...scale analysis for hydrostatic balance! CSU ATS601 Fall 2015
7 Balanced Motion We previously discussed the concept of balance earlier, in the context of hydrostatic balance. Recall that the balanced condition means no accelerations (balance of forces). That is,
More information4. Atmospheric transport. Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017
4. Atmospheric transport Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017 Forces in the atmosphere: Gravity g Pressure-gradient ap = ( 1/ ρ ) dp / dx for x-direction (also y, z directions)
More information1/25/2010. Circulation and vorticity are the two primary
Lecture 4: Circulation and Vorticity Measurement of Rotation Circulation Bjerknes Circulation Theorem Vorticity Potential Vorticity Conservation of Potential Vorticity Circulation and vorticity are the
More informationActual bathymetry (with vertical exaggeration) Geometry of the ocean 1/17/2018. Patterns and observations? Patterns and observations?
Patterns and observations? Patterns and observations? Observations? Patterns? Observations? Patterns? Geometry of the ocean Actual bathymetry (with vertical exaggeration) Continental Continental Basin
More informationGround Rules. PC1221 Fundamentals of Physics I. Position and Displacement. Average Velocity. Lectures 7 and 8 Motion in Two Dimensions
PC11 Fundamentals of Physics I Lectures 7 and 8 Motion in Two Dimensions Dr Tay Sen Chuan 1 Ground Rules Switch off your handphone and paer Switch off your laptop computer and keep it No talkin while lecture
More informationLecture 2 Atmospheric Boundary Layer
Lecture Atmospheric Boundary Layer H. J. Fernando Ariona State niversity Region of the lower atmosphere where effects of the Earth surface are felt Surface fluxes of momentum, buoyancy.. Neutral, Convective,
More informationXI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com
XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K affan_414@live.com https://promotephysics.wordpress.com [MOTION IN TWO DIMENSIONS] CHAPTER NO. 4 In this chapter we are oin to discuss motion in projectile
More informationThe General Circulation of the Atmosphere: A Numerical Experiment
The General Circulation of the Atmosphere: A Numerical Experiment Norman A. Phillips (1956) Presentation by Lukas Strebel and Fabian Thüring Goal of the Model Numerically predict the mean state of the
More informationEddy viscosity. AdOc 4060/5060 Spring 2013 Chris Jenkins. Turbulence (video 1hr):
AdOc 4060/5060 Spring 2013 Chris Jenkins Eddy viscosity Turbulence (video 1hr): http://cosee.umaine.edu/programs/webinars/turbulence/?cfid=8452711&cftoken=36780601 Part B Surface wind stress Wind stress
More informationVorticity in natural coordinates
Vorticity in natural coordinates (see Holton pg 95, section 4.2.) Let s consider the vertical vorticity component only, i.e. ζ kˆ ω, we have ω u dl kˆ ω lim --- lim ----------------- curve is in xy plane
More informationLecture 1. Equations of motion - Newton s second law in three dimensions. Pressure gradient + force force
Lecture 3 Lecture 1 Basic dynamics Equations of motion - Newton s second law in three dimensions Acceleration = Pressure Coriolis + gravity + friction gradient + force force This set of equations is the
More informationd v 2 v = d v d t i n where "in" and "rot" denote the inertial (absolute) and rotating frames. Equation of motion F =
Governing equations of fluid dynamics under the influence of Earth rotation (Navier-Stokes Equations in rotating frame) Recap: From kinematic consideration, d v i n d t i n = d v rot d t r o t 2 v rot
More informationLecture 2. Lecture 1. Forces on a rotating planet. We will describe the atmosphere and ocean in terms of their:
Lecture 2 Lecture 1 Forces on a rotating planet We will describe the atmosphere and ocean in terms of their: velocity u = (u,v,w) pressure P density ρ temperature T salinity S up For convenience, we will
More informationExperiment 3 The Simple Pendulum
PHY191 Fall003 Experiment 3: The Simple Pendulum 10/7/004 Pae 1 Suested Readin for this lab Experiment 3 The Simple Pendulum Read Taylor chapter 5. (You can skip section 5.6.IV if you aren't comfortable
More informationThe Equations of Motion in a Rotating Coordinate System. Chapter 3
The Equations of Motion in a Rotating Coordinate System Chapter 3 Since the earth is rotating about its axis and since it is convenient to adopt a frame of reference fixed in the earth, we need to study
More informationChapter (3) TURBULENCE KINETIC ENERGY
Chapter (3) TURBULENCE KINETIC ENERGY 3.1 The TKE budget Derivation : The definition of TKE presented is TKE/m= e = 0.5 ( u 2 + v 2 + w 2 ). we recognize immediately that TKE/m is nothing more than the
More informationHomework 2: Solutions GFD I Winter 2007
Homework : Solutions GFD I Winter 007 1.a. Part One The goal is to find the height that the free surface at the edge of a spinning beaker rises from its resting position. The first step of this process
More informationConvergence of DFT eigenvalues with cell volume and vacuum level
Converence of DFT eienvalues with cell volume and vacuum level Sohrab Ismail-Beii October 4, 2013 Computin work functions or absolute DFT eienvalues (e.. ionization potentials) requires some care. Obviously,
More informationULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor
ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor Ade afar Abdullah Electrical Enineerin Department, Faculty of Technoloy and Vocational Education Indonesia University of Education
More information(A) (B) (C) (D) None of these
Exercise OBJECTIVE PROBLEMS. Action and reaction (A) act on two different objects (C) have opposite directions. Which fiure represents the correct F.B.D. of rod of mass m as shown in fiure : (B) have equal
More informationDirection of Hurricane Beta Drift in Horizontally Sheared Flows
146 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 54 Direction of Hurricane Beta Drift in Horizontally Sheared Flows BIN WANG, XIAOFAN LI,* AND LIGUANG WU Department of Meteoroloy, School of Ocean and Earth
More informationHomework # 2. SOLUTION - We start writing Newton s second law for x and y components: F x = 0, (1) F y = mg (2) x (t) = 0 v x (t) = v 0x (3)
Physics 411 Homework # Due:..18 Mechanics I 1. A projectile is fired from the oriin of a coordinate system, in the x-y plane (x is the horizontal displacement; y, the vertical with initial velocity v =
More informationGoal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations
Scale analysis relevant to the tropics [large-scale synoptic systems]* Goal: Use understanding of physically-relevant scales to reduce the complexity of the governing equations *Reminder: Midlatitude scale
More informationPhysics 11 Fall 2012 Practice Problems 2 - Solutions
Physics 11 Fall 01 Practice Problems - s 1. True or false (inore any effects due to air resistance): (a) When a projectile is fired horizontally, it takes the same amount of time to reach the round as
More informationOn the effect of forward shear and reversed shear baroclinic flows for polar low developments. Thor Erik Nordeng Norwegian Meteorological Institute
On the effect of forward shear and reversed shear baroclinic flows for polar low developments Thor Erik Nordeng Norwegian Meteorological Institute Outline Baroclinic growth a) Normal mode solution b) Initial
More informationPAPER 333 FLUID DYNAMICS OF CLIMATE
MATHEMATICAL TRIPOS Part III Wednesday, 1 June, 2016 1:30 pm to 4:30 pm Draft 21 June, 2016 PAPER 333 FLUID DYNAMICS OF CLIMATE Attempt no more than THREE questions. There are FOUR questions in total.
More informationControl Volume. Dynamics and Kinematics. Basic Conservation Laws. Lecture 1: Introduction and Review 1/24/2017
Lecture 1: Introduction and Review Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study of motion without regard for the cause. Dynamics: On the other hand, dynamics
More informationLecture 1: Introduction and Review
Lecture 1: Introduction and Review Review of fundamental mathematical tools Fundamental and apparent forces Dynamics and Kinematics Kinematics: The term kinematics means motion. Kinematics is the study
More informationPHYSICS 218 SOLUTION TO HW 8. Created: November 20, :15 pm Last updated: November 21, 2004
Created: November 20, 2004 7:5 pm Last updated: November 2, 2004. Schroeder.6 (a) The three forces actin on the slab of thickness dz are ravity and the pressure from above and below. To achieve equilibrium
More informationSMAST Technical Report The Performance of a Coupled 1-D Circulation and Bottom Boundary Layer Model with Surface Wave Forcing
1 SMAST Technical Report 01-03-20 The Performance of a Coupled 1-D Circulation and Bottom Boundary Layer Model with Surface Wave Forcing Y. Fan and W. S. Brown Ocean Process Analysis Laboratory Institute
More informationInternal boundary layers in the ocean circulation
Internal boundary layers in the ocean circulation Lecture 9 by Andrew Wells We have so far considered boundary layers adjacent to physical boundaries. However, it is also possible to find boundary layers
More informationNote that Rossby waves are tranverse waves, that is the particles move perpendicular to the direction of propagation. f up, down (clockwise)
Ocean 423 Rossby waves 1 Rossby waves: Restoring force is the north-south gradient of background potential vorticity (f/h). That gradient can be due to either the variation in f with latitude, or to a
More informationJ.L. Kirtley Jr. September 4, 2010
Massachusetts Institute of Technoloy Department of Electrical Enineerin and Computer Science 6.007 Electromanetic Enery: From Motors to Lasers Supplemental Class Notes Manetic Circuit Analo to Electric
More information