Tropical Cyclones: Steady State Physics

Similar documents
Adiabatic expansion Isothermal compression Adiabatic compression

The Tropical Atmosphere: Hurricane Incubator

Annual Number of Peer Reviewed Articles with Hurricane or Tropical Cyclone in their Titles, according to Meteorological and Geoastrophysical

Hurricanes are intense vortical (rotational) storms that develop over the tropical oceans in regions of very warm surface water.

Hurricanes: Their physics and relationship to climate. Kerry Emanuel Massachusetts Institute of Technology

Quasi-equilibrium Theory of Small Perturbations to Radiative- Convective Equilibrium States

Impact of Turbulence on the Intensity of Hurricanes in Numerical Models* Richard Rotunno NCAR

Self-Stratification of Tropical Cyclone Outflow. Part II: Implications for Storm Intensification

4.1 LAWS OF MECHANICS - Review

2.1 Effects of a cumulus ensemble upon the large scale temperature and moisture fields by induced subsidence and detrainment

Reconstruction of Thermodynamic Cycles in a High-Resolution Simulation of a Hurricane

Dynamics of the Zonal-Mean, Time-Mean Tropical Circulation

A Reformulation of the Logistic Growth Equation Model (LGEM) for Ensemble and Extended Range Intensity Prediction

Hurricane Science Tutorial. Kerry Emanuel Lorenz Center, MIT

Tropical Cyclone Intensification

AMERICAN METEOROLOGICAL SOCIETY

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

Self-Stratification of Tropical Cyclone Outflow. Part I: Implications for Storm Structure

PV Thinking. What is PV thinking?

PAPER 68 ACCRETION DISCS

Physical Fundamentals of Global Change Processes

1. Second Law of Thermodynamics

Section 2: Lecture 1 Integral Form of the Conservation Equations for Compressible Flow

Summarizing, Key Point: An irreversible process is either spontaneous (ΔS universe > 0) or does not occur (ΔS universe < 0)

ATMOS 5130 Lecture 9. Enthalpy Conservation Property The Second Law and Its Consequences Entropy

Thermodynamic and Flux Observations of the Tropical Cyclone Surface Layer

USING PROGRAM HURRICANE. Download all the files, including the (empty) output subdirectory into a new folder on your machine.

Lectures on Tropical Cyclones

1. Second Law of Thermodynamics

Turbomachinery & Turbulence. Lecture 2: One dimensional thermodynamics.

2 Equations of Motion

Dust devils, water spouts, tornados

Exam 1 (Chaps. 1-6 of the notes)

Exercise Set 4. D s n ds + + V. s dv = V. After using Stokes theorem, the surface integral becomes

Introduction to Chemical Engineering Thermodynamics. Chapter 7. KFUPM Housam Binous CHE 303

Lecture 35: Vapor power systems, Rankine cycle

where p oo is a reference level constant pressure (often 10 5 Pa). Since θ is conserved for adiabatic motions, a prognostic temperature equation is:

Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics

The long-lived concentric eyewall tropical cyclones with large moat and outer eyewall

X Y. Equator. Question 1

Brief Overview of the Global Atmosphere

8.21 The Physics of Energy Fall 2009

Energy Production, Frictional Dissipation, and Maximum Intensity. of a Numerically Simulated Tropical Cyclone

ρ 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

3 Hydrostatic Equilibrium

Applied Thermodynamics for Marine Systems Prof. P. K. Das Department of Mechanical Engineering Indian Institute of Technology, Kharagpur

df dz = dp dt Essentially, this is just a statement of the first law in one of the forms we derived earlier (expressed here in W m 3 ) dq p dt dp

On the validity of representing hurricanes as Carnot heat engine

Review of First and Second Law of Thermodynamics

( ) = 1005 J kg 1 K 1 ;

To study the motion of a perfect gas, the conservation equations of continuity

Classical thermodynamics

Forecasting Hurricane Intensity: Lessons from Application of the Coupled Hurricane Intensity Prediction System (CHIPS)

Gasdynamics 1-D compressible, inviscid, stationary, adiabatic flows

n v molecules will pass per unit time through the area from left to

Chapter 2 Thermodynamics

Model description of AGCM5 of GFD-Dennou-Club edition. SWAMP project, GFD-Dennou-Club

UNIVERSITY OF SOUTHAMPTON

Radiative Transfer Chapter 3, Hartmann

Practice Examinations Chem 393 Fall 2005 Time 1 hr 15 min for each set.

2σ e s (r,t) = e s (T)exp( rr v ρ l T ) = exp( ) 2σ R v ρ l Tln(e/e s (T)) e s (f H2 O,r,T) = f H2 O

Hurricane potential intensity from an energetics point of view

Isentropic Efficiency in Engineering Thermodynamics

4. Atmospheric transport. Daniel J. Jacob, Atmospheric Chemistry, Harvard University, Spring 2017

Consider a volume Ω enclosing a mass M and bounded by a surface δω. d dt. q n ds. The Work done by the body on the surroundings is

Review of fluid dynamics

Thermodynamics! for Environmentology!

CHAPTER 3 LECTURE NOTES 3.1. The Carnot Cycle Consider the following reversible cyclic process involving one mole of an ideal gas:

Outline. Property diagrams involving entropy. Heat transfer for internally reversible process

ME Thermodynamics I

Robert Rogers, Sylvie Lorsolo, Paul Reasor, John Gamache, and Frank Marks Monthly Weather Review January 2012

ATMO 551a Moist Adiabat Fall Change in internal energy: ΔU

Quantifying the Role of Tropospheric Relative Humidity on the Development of Tropical Cyclones

Chapter 7. Entropy. by Asst.Prof. Dr.Woranee Paengjuntuek and Asst. Prof. Dr.Worarattana Pattaraprakorn

5/6/ :41 PM. Chapter 6. Using Entropy. Dr. Mohammad Abuhaiba, PE

Fundamentals of compressible and viscous flow analysis - Part II

Irreversible Processes

Energy Production, Frictional Dissipation, and Maximum Intensity of a Numerically Simulated Tropical Cyclone*

MATTER TRANSPORT (CONTINUED)

Chapter 2. Quasi-Geostrophic Theory: Formulation (review) ε =U f o L <<1, β = 2Ω cosθ o R. 2.1 Introduction

First Law CML 100, IIT Delhi SS. The total energy of the system. Contribution from translation + rotation + vibrations.

Chapter 1. Governing Equations of GFD. 1.1 Mass continuity

2013/5/22. ( + ) ( ) = = momentum outflow rate. ( x) FPressure. 9.3 Nozzles. δ q= heat added into the fluid per unit mass

The Bernoulli Equation

OCN 623: Thermodynamic Laws & Gibbs Free Energy. or how to predict chemical reactions without doing experiments

Where does Bernoulli's Equation come from?

Outline. Aim. Gas law. Pressure. Scale height Mixing Column density. Temperature Lapse rate Stability. Condensation Humidity.

Chapter Two. Basic Thermodynamics, Fluid Mechanics: Definitions of Efficiency. Laith Batarseh

In this section, mathematical description of the motion of fluid elements moving in a flow field is

Using the Entropy Rate Balance to Determine the Heat Transfer and Work in an Internally Reversible, Polytropic, Steady State Flow Process

1. Thermodynamics 1.1. A macroscopic view of matter

Entropy and the Second Law of Thermodynamics

Thermodynamics of phase transitions

( ) is zero, and thus (remembering that Q H < 0),

Where is the PdV term in the first law of black-hole thermodynamics?

ME 2322 Thermodynamics I PRE-LECTURE Lesson 23 Complete the items below Name:

Inner core dynamics: Eyewall Replacement and hot towers

Irreversible Processes

From the last time, we ended with an expression for the energy equation. u = ρg u + (τ u) q (9.1)

CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION

Transcription:

Tropical Cyclones: Steady State Physics

Energy Production

Carnot Theorem: Maximum efficiency results from a particular energy cycle: Isothermal expansion Adiabatic expansion Isothermal compression Adiabatic compression Note: Last leg is not adiabatic in hurricane: Air cools radiatively. But since environmental temperature profile is moist adiabatic, the amount of radiative cooling is the same as if air were saturated and descending moist adiabatically. Maximum rate of working: W T T s o = Q T s

Total rate of heat input to hurricane: r0 ( * ) 3 0 k 0 D Q = π ρ C V k k + C V rdr Surface enthalpy flux Dissipative heating In steady state, Work is used to balance frictional dissipation: r 0 0 3 W = π ρ C V rdr D

Plug into Carnot equation: T T C rdr C ( k ) 0 k = rdr r0 r 3 0 s o * ρ 0 D V ρ 0 k V To If integrals dominated by values of integrands near radius of maximum winds, C k T s T V o * max k 0 k C D T o ( )

Problems with Energy Bound: Implicit assumption that all irreversible entropy production is by dissipation of kinetic energy. But outside of eyewall, cumuli moisten environment...accounting for almost all entropy production there Approximation of integrals dominated by high wind region is crude

Local energy balance in eyewall region:

Definition of streamfunction, ψ: ρ rw ψ =, ρru = r ψ z Flow parallel to surfaces of constant ψ, satisfies mass continuity: 1 ( ) 1 ρru + ( ρrw) = r r r z 0

Variables conserved (or else constant along streamlines) above PBL, where flow is considered reversible, adiabatic and axisymmetric: 1 E = cpt + Lvq+ gz+ V Lvq * s* = cplnt Rd ln p+ T Energy: Entropy: Angular Momentum: 1 M = rv + fr

First definition of s*: Tds* = c dt + L dq* α dp p v (1) Steady flow: p p αdp α = dr + α dz r z Substitute from momentum equations: V αdp = dz[ g + Vi w] + dr + fv Vi u r ()

Identity: 1 1 wdz+ udr= du + w + d, ρr ς ψ u w ς azimuthal vorticity z r ( ) ( ) ( Vi Vi ) where (3) Substituting (3) into () and the result into (1) gives: V 1 Tds* = de VdV + fv dr + d r ρr ς ψ (4)

One more identity: V M 1 VdV + + fv dr = f dm r r Substitute into (4): M 1 1 Tds* dm d d E fm r ρr ξ ψ + = +, (5) Note that third term on left is very small: Ignore

Integrate (5) around closed circuit: Right side vanishes; contribution to left only from end points

1 1 T T ds* + MdM = 0 r r b o ( ) b o 1 1 ds * = ( T ) b T (6) o r r MdM b o Mature storm: r >> r : o b b M ds * ( T ) b To r dm (7)

In inner core, V >> fr M rv ds* V r ( T T ) b b b o (8) dm Convective criticality: * b s = ds (9) ( ) b Vb rb Tb To dm s

ds b /dm determined by boundary layer processes:

Put (7) in differential form: ds M dm ( ) + = 0, Tb To dt r dt (10) Integrate entropy equation through depth of boundary layer: ds 1 h = C ( k k) C F dt T V + V + s * 3 k s D b, (11) where F b is the enthalpy flux through PBL top. Integrate angular momentum equation through depth of boundary layer: dm h = CDr V V, (1) dt

Substitute (11) and (1) into (10) and set F b to 0: C ( * ) k T s T V = o k 0 k (13) C D T o Same answer as from Carnot cycle. This is still not a closed expression, since we have not determined the boundary layer enthalpy, k. We can do this using boundary layer quasiequilibrium as follows. First, use moist static energy, h instead of k: C ( * ) k T s T V = o h 0 h (14) C D T o h k+ gz = c T + L q+ gz p v

Convective neutrality: h = h* b First law of thermo: dh = Tds + RTdln p, * * b b d b (15) Go back to equation (10): ds M dm Tb To dt r dt ( ) + = 0, Use definition of M and gradient wind balance: R T d ln r p V = + r fv

( ) 1 1 1 V = + + + + r 1 1 1 = d V + frv + f r + RdTln p. 4 * Tb To ds d V frv f rdr fv dr Substitute into (15): T 1 1 1 dh = d V + frv + f r + R T ln p * b b d o Tb T o 4 Define an outer radius, r a, where V=0, p=p o. Take difference between this and radius of maximum winds: T 1 p 1 h = h V + fr V + R T ln f r ( ) * * b m b a max m max d o a p Tb T o a 4

Relate p m to V max using gradient wind equation. But simpler to use an empirical relation: RT ln p bv p m d s max a Also neglect fr m in comparison to V max. and neglect difference between T s and T b : Substitute into (14): V max T T ( ) 1 T h h f r C T T C D C k 1 T s 1 b CD To s o * * s s a a k o 4 o.

Absolute upper bound on storm size: r omax = 4 T s T T s o ( * h ) s ha f romax 1000 km For r o << r omax neglect last term in numerator. C T T V h h (16) ( ) k s o * * max s a CD To

h h = h h = c T T + L q q ( H ) ( ) ( ) * * * * s a s ab p s b v s ab Lq v * sa 1 e*( Ts ) = L ε v 1 p o ( H ) Since p depends on V, this makes (16) an implicit equation for V. But write expression for wind force instead: p C T T L ρv = V = e*( T ) 1 H o k s o v max max s RT d s CD To RT v s o ( )

Otherwise, use empirical relationship between V and p: RT d p= s p ln p e a b V p a bv RT d s C T T C T T e*( T ) V h h = L 1 ( ) ε ( H ) o k s o * * k s o s max s a v CD To CD To p C T T pe V L e*( T) 1 H V b max RT d s k s o a max vε s CD To ( ) o

Taking natural log of this, the result can be written in the form: ( ) max max AV ln V B = 0

Numerical simulations

Relationship between potential intensity (PI) and intensity of real tropical cyclones

Wind speed (m/s) Potential wind speed (m/s) 80 70 60 50 V (m/s) 40 30 0 10-80 -60-40 -0 0 0 40 60 80 100 Time after maximum intensity (hours)

Evolution with respect to time of maximum intensity

Evolution with respect to time of maximum intensity, normalized by peak wind

Evolution curve of Atlantic storms whose lifetime maximum intensity is limited by declining potential intensity, but not by landfall

Evolution curve of WPAC storms whose lifetime maximum intensity is limited by declining potential intensity, but not by landfall

CDF of normalized lifetime maximum wind speeds of North Atlantic tropical cyclones of tropical storm strength (18 m s 1) or greater, for those storms whose lifetime maximum intensity was limited by landfall.

CDF of normalized lifetime maximum wind speeds of Northwest Pacific tropical cyclones of tropical storm strength (18 m s 1) or greater, for those storms whose lifetime maximum intensity was limited by landfall.

Evolution of Atlantic storms whose lifetime maximum intensity was limited by landfall

Evolution of Pacific storms whose lifetime maximum intensity was limited by landfall

Composite evolution of landfalling storms