Chapter 4: Fundamental Forces Newton s Second Law: F=ma In atmospheric science it is typical to consider the force per unit mass acting on the atmosphere: Force mass = a In order to understand atmospheric motion accelerations) we need to know what forces act on the atmosphere. Body and Surface Forces Body or volume) force a force that acts on a volume of the atmosphere Surface force a force that acts locally upon a part of a fluid Forces in an Inertial Reference Frame Inertial reference frame a non-accelerating frame of reference Gravity a force that arises from the mutual attraction between two objects g # = GM &# r % % & $ $ r G gravitational constant = 6.673x10-11 N m 2 kg -2 ) M mass of the Earth kg) r distance between objects For a unit mass of air at the surface of the Earth: M = 5.9742x10 24 kg mass of Earth) r = 6.37x10 6 m radius of Earth) g 0 = 9.83 m s -2 r 2
Geopotential Φ): work required to raise a unit mass to height z z Φz) = g z)dz 0 Note that g varies with z since it depends on the distance from the center of the Earth. Geopotential height Z): Z = Φ g 0 = 1 g 0 z 0 g dz Pressure Gradient Force Example: Real-time weather map The force exerted on the left face of this air parcel due to pressure is: pa = pδyδz The force exerted on the right face of this air parcel due to pressure is: % p + δp)δyδz = p + p & x δx δyδz )
The net force exerted by pressure on this air parcel is: p x δxδyδz By dividing by the mass of the air parcel ρδxδyδz) we get the force per unit mass due to changes in pressure i.e. the pressure gradient force): P g i = 1 p ρ x We can write all three components of the pressure gradient force as: P g = 1 % p i + p j + p k ρ & x y z ) In what direction does this force act relative to locations with high and low pressure? Example: Calculation of the pressure gradient force from a weather map Viscous Force
Each layer in this fluid is experiencing a shear stress due to differences in the speed of motion of adjacent layers of fluid. τ i = µ u z µ - coefficient of shear viscosity depends on fluid being considered) How does the shear stress change between adjacent layers of the fluid? What does this imply about the net force acting on a layer of this fluid? For a net force to arise there must be a gradient in the shear stress. The force that arises due to a gradient in the shear stress is called a viscous force and can be represented as: F r i = µ $ 2 u ρ x + 2 u 2 y + 2 u & % 2 z 2 ) Kinematic viscosity coefficient: ν = µ ρ ~ 1.5x10-5 m 2 s -1 in the atmosphere) In three dimensions this force is written as: $ 2 u x + 2 u 2 y + 2 u & 2 z 2 ) -, i + /,% /, $ 2 v F r = ν x + 2 v 2 y + 2 v & % 2 z 2 ) /, j + /, /, $ 2 w x + 2 w 2 y + 2 w, & % 2 z 2 ) k / / +,./
Hydrostatic Balance For an atmosphere that is at rest, what forces act in the vertical direction? 1 dp ρ dz = g dp dz = ρg This equation is known as the hydrostatic equation. Hydrostatic balance a balance between the vertical pressure gradient force and the gravitational force Combining the hydrostatic equation with the ideal gas law gives: " ln p % $ = # p 0 & z 0 g R d T dz Using a layer average temperature, T, we can write: where H = R d T g p p 0 e z H This equation is known as the hypsometric equation and relates pressure and height.
Using the definition of geopotential we can also write: dp p = g R d T dz dφ = R d T dp p ΔΦ = R d p Td ln p p 0 ΔZ = R d g 0 p Td ln p p 0 Thickness ΔZ) the difference in geopotential height Z) between two pressure levels What causes the thickness to increase or decrease)? Example: Thickness and atmospheric soundings Forces in a Rotating Reference Frame Non-inertial frame of reference a frame of reference which is undergoing an acceleration Why is a frame of reference fixed to the Earth a non-inertial frame of reference? Angular velocity of rotation of the Earth: Ω = 2π 1 sidereal day = 2π 86164 sec = 7.292 10 5 s 1 Sidereal day the time required for the earth to make one complete revolution in an absolute fixed) coordinate system, that is, with respect to the stars.
Centrifugal Force Rotating table experiment: What happens to the object when it is placed on the rotating table? What force is required to keep the object stationary on the rotating table? mrω 2 Centrifugal force outward radial force Centripetal force inward radial force Is this object experiencing an acceleration when viewed in a frame of reference rotating with the table? Is this object experiencing an acceleration when viewed from a fixed frame of reference? Which forces are acting in each frame of reference? The centrifugal force arises only in observations taken in a rotating frame of reference and is due to the acceleration of the frame of reference.
The centrifugal force and the Earth: What is the direction of the centrifugal force for an object on the Earth? The centrifugal force is a body force and can be combined with the gravity force to give an effective gravity. Effective gravity: g = 9.81 m s -2 at sea level For atmospheric science applications we use effective gravity g) rather than gravity g ). Therefore effective gravity should be used in place of gravity in all of the previous equations in this chapter. The addition of the centrifugal force to the gravity force results in the effective gravity force not being directed towards the center of the Earth. In reality the Earth is not a perfect sphere and the effective gravity is exactly normal to the surface of the Earth at all locations neglecting topography). Therefore, the effective gravity force acts only in the z- direction. This discussion has only considered an object at rest in a rotating frame of reference. What happens when an object is in motion relative to a rotating frame of reference?
Coriolis Force For an object in motion relative to a rotating frame of reference we need to consider: - an additional centrifugal force - changes in the relative angular momentum of the object What is the centrifugal force acting on an object moving towards the east? CF total = m Ω + u 2 $ & ) R φ % R φ = mω 2 R φ + 2mΩu R φ R φ + m u2 R φ 2 The first term on the RHS of this equation is the same as the centrifugal force for an object at rest, and is incorporated into the effective gravity. The last term on the RHS of this equation is assumed to be small since R φ >> u. The remaining term on the RHS of this equation can be written as: R φ CF motion = 2mΩusin φ j + 2mΩucosφk and is the centrifugal force due to the motion of the object that we need to consider in our study of the atmosphere. This component of the centrifugal force arises only for zonal east/west) motion.
Conservation of angular momentum For a rotating object the angular momentum of the object is given by: Iω = 1 mr 2 2 φ Ω The angular momentum Iω) of this object will remain constant. What happens to the terms in this equation if the object is moving towards the North Pole v>0)? ) 2 Ω + 1 mr 2 2 φ Ω = 1 m R 2 φ + δr % & δu R φ + δr) R φ 2 Ω = R φ 2 Ω + 2R φ δrω + R 2 φδu R φ + δr As the object moves towards the pole its zonal velocity u) increases. Noting that δr = vsin φδt we get: As δt 0 we then have: δu = 2Ωvsin φδt du dt = 2Ωvsin φ Conservation of angular momentum results in a zonal acceleration for objects moving in a meridional north/south) direction. Similarly, conservation of angular momentum also results in objects moving in the vertical direction experiencing a zonal acceleration equal to: du dt = 2Ωwcosφ
Conservation of angular momentum does not lead to an acceleration for objects moving in a purely zonal direction. The accelerations due to the additional component of the centrifugal force and due to conservation of angular momentum for an object in motion in a rotating frame of reference can be combined to form a single force known as the Coriolis force. F Coriolis m = 2Ωvsin φ 2Ωwcosφ ) i 2Ωusin φ j + 2Ωucosφk Using vector notation, and noting that Ω = Ωcosφ j + Ωsin φk the Coriolis force can also be expressed as: F Coriolis m = 2 Ω u In what direction does the Coriolis force act relative to the motion of the air? Does this differ in the Northern and Southern hemispheres? Example: The Coriolis force on the roof of Duane Physics The Navier-Stokes Equation The sum of the forces acting in each direction gives the acceleration experienced by an air parcel. Du Dt = Dv Dt = 1 p ρ x + ν & 2 u x + 2 u 2 y + 2 u) 2 z 2 + + 2Ωvsin φ 2Ωwcosφ) 1 p ρ y + ν & 2 v x + 2 v 2 y + 2 v) 2 z 2 + 2Ωusin φ) Dw Dt = g 1 p ρ z + ν & 2 w x + 2 w 2 y + 2 w) 2 z 2 + + 2Ωucosφ) What does each term in these equations represent?
These equations are known as the Navier-Stokes equations. We can rewrite the Navier-Stokes equations in an Eulerian framework as: u t + u u x + v u y + w u z = v t + u v x + v v y + w v z = 1 p ρ x + ν & 2 u x + 2 u 2 y + 2 u) 2 z 2 + + 2Ωvsin φ 2Ωwcosφ) 1 p ρ y + ν & 2 v x + 2 v 2 y + 2 v) 2 z 2 + 2Ωusin φ) w t + u w x + v w y + w w z = g 1 p ρ z + ν & 2 w x + 2 w 2 y + 2 w) 2 z 2 + + 2Ωucosφ) Perturbation Pressure We can rewrite the vertical momentum equation by defining a reference pressure, p 0, that is in hydrostatic balance with a reference density, ρ 0, as: The total pressure is then given by: dp 0 dz = ρ 0g p = p o z) + p d Substituting this definition into the vertical component of the Navier-Stokes equation gives: ρ Dw Dt = ρg z p 0 z ) + p d ) + ρν 2 w & x + 2 w 2 y 2 + 2 w z 2 ) + + ρ 2Ωucosφ) = ρg + ρ 0 g p d z + ρν & 2 w x + 2 w 2 y + 2 w) 2 z 2 + + ρ 2Ωucosφ) ) Dw Dt = ρ ρ 0 ρ g 1 p d ρ z + ν & 2 w x + 2 w 2 y + 2 w) 2 z 2 + + 2Ωucosφ )