Pedagogical Mathematica Platform Visualizing the Coriolis Effects in 3-Cell Atmospheric Circulation Model

Journal of Astronomy and Space Sciences.
2014.
Mar,
31(1):
91-99

This is an open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/) which permits unrestrictednon-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

- Received : January 27, 2013
- Accepted : February 20, 2014
- Published : March 15, 2014

Download

PDF

e-PUB

PubReader

PPT

Export by style

Share

Article

Metrics

Cited by

TagCloud

The atmospheric flow in the 3-Cell model of global atmosphere circulation is described by the Lagrange's equation of the non-inertial frame where pressure force, frictional force and fictitious force are mixed in complex form. The Coriolis force is an important factor which requires calculation of fictitious force effects on atmospheric flow viewed from the rotating Earth. We make new Mathematica platform to solve Lagrange's equation by numerical analysis in order to analyze dynamics of atmospheric general circulation in the non-inertial frame. It can simulate atmospheric circulation process anywhere on the earth. It is expected that this pedagogical platform can be utilized to help students studying the atmospheric flow understand the mechanisms of atmospheric global circulation.
Mathematica
is helpful for the convenient function of symbolic calculations and graphic manipulations. The
Mathematica
simulation of the atmospheric circulation matching the 3-Cell GCA model has been presented in our previous work (
Yun 2006
), however, the 3-Cell GCA model was shown with the illustration rather than with active platform. In this paper, we present anew
Mathematica
platform presenting the GCA simulation according to the 3-Cell GCA model in 2D or 3D graphics. In the platform, we can simulate the Coriolis effects at any point of the globe automatically and confirm the atmospheric circulation dynamics interactively.
_{o}
for the Lagrange's equation in the non-inertial frame of reference. This transformation is done in two steps. Firstly, we consider a frame of reference K' which moves with a translational velocity
relative to the ineltlal frame K
_{o}
. Next, we bring a new frame of K which rotates relative to K' with angular velocity,
. As a result, K executes both translational and rotational transformation to the inertial frame K
_{o}
(fix star):
. The Lagrangian in K
_{o}
frame is (
Landau & Lifshitz 1976
)
where, V
_{o}
is the velocity of a particle in K
_{o}
frame and U is a potential. The velocities
come from a transformanon
and the
come from a transfonnation
Finally, me Lagrangian in K frame is
where,
Then Lagrange's equation, which satisfy the principle ofleast action,
If we suppose that angular velocity is constant and neglect dle translational velocity, we omit
and
. Then Eq. (3) is
We rewrite this again as
This is just a Newton's equation, which shows the motion of an object by the force
Here
includes the friction and the pressure gradient force. and the
includes the CorioHs force 2m
wand the centrifugal force m(
) ×(
) which deflect the atmospheric flows on the rotating Earth. As shown in
Fig.1
, the CorioJis acceleration vector
direct to ±Y direction. Deriving process of the Lagrange's equation in the non-inertial frame is well described in other books (
Symon 1971
.
Landau & Lifshitz 1976
).
Actual and fictitious forces on the rotating Earth's surface at the latitude of λ = π/2-θ in a topocentric frame of XYZ axis system. The angular velocity of the Earth is ω_{0} and θ is the polar angle in the geocentric reference frame.
point[t, θ]
in solving the Eq. (4) in
Mathematica
(
Zimmerman & Olness 1995
,
Tam 1997
). If a latitude (λ = π /2 -θ) is given, the wind vector function
[θ] is determined by the position function
point[t, θ]
. We assume that the initial wind blows along the meridian and its velocity is determined according to the 3-Cell GCA model.
Mathematica
coding for solving Eq. (4) is shown in In[1] - In[10] and its
Mathematica
solution is Out[11]:
point[t, θ]
=
{t v0x, -t^2ω0x Cos[θ], -1/2 gt^2}
while n Order = 0 for simplicity. Once
point[t, θ]
is given, the
Parametric Plot
draw the path of the wind with a time domain vector array of the
point[t, θ]
in Mathematica such as that shown in
Fig. 2
. The Out[22] (
Fig. 2
) is a list of plot in deflection on the eight points on the globe selected respectively at the different atmospheric zone of the 3-Cell GCA model. As shown in
Fig. 2
, the winds deflect right in the Northern Hemisphere and deflect left in the Southern Hemisphere regardless of the wind direction. No deflection occurs at the equator (θ = π /2) because the Coriolis force m
= 0 at the equator such as that shown in
Fig. 1
. It was not until comprehension of the vectorial nature of the effective force in the non-inertial frame of reference that we could analyze the deflection effects of the atmospheric circulations effectively; we present some parts of the Mathematica coding to show the vector calculations and the results with their vector components.
The table of the deflections at eight points on the globe, The points are selected In the different atmospheric zone of the 3-Cell GCA model.
However,
Fig. 2
does not show which fictitious forces cause the deflections unless we calculate those with their vector components respectively. Because the deflections of the winds over the globe come from the resultant effective force
in Eq. (5), we must calculate accurately the ingredients of the effective force for the probable cause of deflections. For the examination of the deflection effects of the effective force, we calculate the constituent parts of the effective force accelerations with its vector components at the eight cities respectively, and summarize those in
Table 1
. The vector calculation of the fictitious forces with those components is easy in
Mathematica
as shown coding In[31]. The vector command of the
vector product
or
triple vector product
is coding as cross[a, b] or cross[a, cross[b, c]] in
Mathematica
as if we write down the formula in text. Resultant accelerations occur in X, Y, Z directions, however, only the Coriolis force acceleration occur in ±Y direction (In[31], Out[31]) which acts as a deflections force of the wind flowing ±X direction. The effective force accelerations in ±X, Z directions reduce or enhance of wind speed and gravity (In[31], Out[31-33]), those are not forces deflecting the wind direction for the wind speed in ±X, gravity in Z. These calculations confirm that the Coriolis force is the unique deflection force in the atmospheric circulation dynamics even though its relative intensity is so weak;
Coriolis force: centrifugal force: gravity
= 10
^{-4}
: 10
^{-2}
: 1 as shown in
Table 1
. Cyclone is an area of low pressure around which the winds blow parallel to the pressure gradient force balanced to the Coriolis force. A hurricane is a severe tropical cyclone which occurs over the northern Atlantic and eastern North Pacific oceans, and it is called the typhoon in the western North Pacific. This same type of cyclone has a different name in a different region of the world, however, all the climate events occur by the Coriolis force effect (
Ahrens 2001
, AMNH-Weather and Climate Events 2014).
^{a}deflection to East (+) or West (+), ^{b}South (+) or North (-), ^{c}no factor
Parametric Plot
of Mathematica. On starting the program, platform will show Fig. 3a. Simulation performs the built-in program when you click the ► appearing while you spread the ⊕ of
t
panel of the platform. While the program is executing, you can change the parameters of
Wind Speed, Earth's Rotation
, and
Theta
, then the change will be effective immediately by the platform. When you stop the simulation, then the platform present the trace so far, and you can change the parameters of the platform to see another simulation. We look at the different deflection in a simulation if we only click the different point in the status of pausing platform.
Fig. 4
shows the deflections at six points both 2D and 3D Graphics at once by clicking the points on the platform panel. It shows promptly winds of the GCA 3-Cell model on the one platform: two polar easterlies, two westerlies, and two trade winds. The performance mode of the platform will change to the 2D Graphics or 3D Graphics anytime you want. Manipulation of both modes enables you to analyze not only 2-dimensional deflections but also 3-dimensional deflections. For example, at the equator we can analyze the X, Z deflection as the wind velocity is varying along the ±X direction and confirm the deflection effects such that as shown in
Fig. 5
. Because the program simulate the position function
point[t, θ]
of solution of the vector differential equation with effective force, the simulation performs the physics behavior of three vectors – wind velocity
angular velocity of the Earth (
), and position vector
– and their associated effects accurately. Hence, the program differs from the animation that animate the assign functions with the Animate function from the convenient graphic tool (
Zeleny 2010
). In
Mathematica
, we can save the snapshots of the simulations and print it. This program also executes well on later version of 8.0
Mathematica
.
Snapshots of the Mathematica platform of the general atmospheric circulation to the 3-Cell model. (a) is the 2D Graphics presentation and (b) is that of 3D Graphics. Simulations begins on click ► in the popup when you click the ⊕ right of t panel. Anytime you may change the Graphics dimension and Graphics menu by selecting the panel menu. In addition, you may change the Graphics mode even if the simulation is in pause, then the platform shows the modified simulation.
The wind deflections at the six points of the globe in both 2D Graphics (a) and 3D Graphics (b). The six points are selected for the GCA 3-Cell model from the North pole: θ = 1/8 π: North-Polar Easterlies, 2/8 π: South-West Westerlies, 3/8 π: North-East Trade Winds, 5/8 π: South-East Trade Winds, 6/8 π: North-West Westerlies, 7/8 π: South-Polar Easterlies in both columns.
Wind deflections at the equator in 2D (left column) graphics and 3D graphics (right column) for the different wind velocities.
Here the Coriolis force 2 m
wand centrifugal force m
× (
) are the most important fictitious forces which play a significant role in a variety of natural processes, most prominently the atmospheric circulation dynamics. In particular, the Coriolis force is also responsible for the circular motion of the tropical cyclone, the tropical typhoon and the hurricane. Physical comprehension and manifesting ability for the non-inertial frame of reference on the rotating Earth becomes a merit of asset to physicist and natural scientist. The conception of the inertial forces on the Earth has been recognized recently because the Coriolis force is not only the most important effective force in the GCA dynamics but also an indispensable task for the scientist in long-range missile launching, satellite operation and GPS position sensors in modern technology. We demonstrate the Coriolis effects platform to simulate the wind progressing on the rotating Earth's surface matching the 3-Cell GCA model in
Mathematica.
The platform draws the realistic path of the wind, through
point[t, θ]
with the time domain of the vector array of the solution of Eq. (4) using the
Parametric Plot
in
Mathematica.
We expect this platform will be a helpful tool for the physicist and the scientist to analyze the atmospheric general circulation dynamics.

1. INTRODUCTION

Motion observed on the rotating Earth is generally explained by invoking inertial forces described in the noninertial frame of reference (
Symon 1971
,
Landau & Lifshitz 1976
). Because viewing the motion from accelerating or rotating frame of reference introduces fictitious forces added to actual forces. On account of the inequality in energy absorbed at a spherical surface of the rotating Earth, seven zones of atmospheric pressure are formed over the Earth surface: an intertropical convergence zone (ITCZ) near the equator, two subtropical highs in both Hemispheres at the latitude of 30°, two subpolar lows on both Hemispheres at the latitude of 60°, two polar highs on both poles. Accordingly, global scale circulation of the atmosphere is described by the Lagrange's equation in the non-inertial frame of reference according to the 3-Cell general circulation of the atmosphere (GCA) model rather than Hadley Cell model (
Ahrens 2001
). Long-term recording data from the satellites approve of the 3-Cell GCA model (“AMNH-Weather and Climate Events” 2014, “Global Climate Animation” 2014). On the rotating Earth frame, the Coriolis force acts as a most important force to change the direction of surface airflows on the Earth. The deflection is not only instrumental in large-scale atmospheric circulations, the development of tropical cyclones, hurricanes and typhoons, also it can affect missile launching, satellite operation, and GPS position sensors (
Bikonis & Demkovicz 2013
) in the modern sciences. Effects upon the weather, ocean currents, rivers and projectile motions are well documented (
Graney 2011
,
Mclntyre 2000
), but the motions over very long distance are required for discernible effects. Common pedagogical tools are helpful to explain the Coriolis effects. For example, Merry- Go-Round table (“Merry-Go-Round” 2014) or Bath-Tub Vortex (
Trefethen et al. 1965
) is helpful for explaining the Coriolis effects, but it cannot be distinguished whether its effect is from Coriolis force or centrifugal force unless we calculate the forces with their vector components. While this approach simplifies some problems, there is often little physical insight into the motion, in particular, into the fictitious force of the vectorial characteristic.
Recently, efficient visualization programs are utilized for the Coriolis force effects (
Zimmerman & Olness 1995
,
Tam 1997
,
Yun 2005
,
Zeleny 2010
). In particular,
2. LAGRANGE’S EQUATION IN A NON-INERTIAL FRAME OF REFERENCE

In the inertial frame, space should be homogeneous and lime is isotropic to assert the invariance of the mechanical sysrem. If we were to choose an arbirrary frame of reference, space would be inhomogeneous and anisotropic. Therefore, the equation of motion in the rotating Earth system should be described in the non-inertial frame ofreference (
Landau & Lifshitz 1976
). For the validity of the principle of least action in the mechanical system independent of the frame of reference chosen, we must carry out the necessary transformation of the Lagrangian L
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

3.MATHEMATICAPROGRAMING FOR GERNERAL CIRCULATION OF THE ATMOSPHERE

- 3.1 Vector calculation of the effective force

To analyze the deflection effect of the atmospheric flow based on the 3-Cell GCA model at any point on the Earth's surface, we create a position function
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

The accelerations of the effective force calculated with those vector components at the eight cities. The calculation parameters: ω0= 7.292 × 10-5sec-1, g=9.8 m/sec2, and the unit of the acceleration is m/sec2. Minus sign of thestands for the wind direction to the north unit of m/sec.

PPT Slide

Lager Image

- 3.2Mathematicaplatform for the 3-cell atmospheric general circulation model

Understanding the vectorial nature of the effective forces on the Earth's surface is a keyword of the atmosphere circulation dynamics. However, it is not easy to evaluate the wind deflection effects at any point on the Earth's surface, because the deflection is varying on the resultant wind vector and associated effective force vectors at the point belong to the atmospheric pressure zone matching to 3-Cell GCA model. The fundamental concept of the inertial force effects in the non-inertial frame of reference is essential to analyze the GCA dynamics. Visual representation of such vectorial nature of the GCA will be in valuable to the researchers working in this field and help teach physics or meteorology students.
We provide Coriolis effects platform to simulate the wind progressing on the rotating Earth's surface matching the 3-Cell GCA model using the function Manipulate in Mathematica.
Fig. 3
a is a 2D Graphics snapshot of the platform and
Fig. 3b
is a 3D Graphics snapshot of the platform. The platform draws the path of the wind through the point[t, θ] with a time domain of vector array of the solution of Eq. (4) using the
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

4.SUMMURY

The Earth is a planet of a sphere which revolves the fix star while rotating itself with a constant angular velocity. The frame of reference on the rotating Earth is a non· inertial frame of reference since the frame is in acceleration continuously. Therefore, the equation of motion should be modified through the coordinate's lransformarions by the least action principle in the mechanical system of Eq. (4):
PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

PPT Slide

Lager Image

Acknowledgements

This research was partially supported by the ReSEAT program fund by the Korean Ministry of Education, Science and Technology through the National Research Foundation of Korea and the Korea Lottery Commission grants. REFERENCES

Ahrens CD
2001
Essentials of Meteorology
3rd
BROOKS-COLE
Pacific Grove
180 -

AMNH-Weather and Climate Events [Internet]
http://www.amnh.org/sciencebulletins/climate

Global Climate Animations [Internet]
Global Climate Animations
http://geogrphy.uoregon.edu/envchange/climate

Landau LD
,
Lifshitz EM
1976
Mechanics
3rd
Pergamon Press
Oxford

Merry-go-round table [Internet]
http://www.nasa.gov/audience/forstudents/brainbites/nonflash/bb_home_corioliseffect.html

Symon KR
1971
Mechanics
3rd
Addison-Wesley
Boston
279 -

Tam PT
1997
A Physicist's Guide to Mathematica
Academic Press
San Diego
349 -

Yun HJ
2005
Interactive Vector Analysis and Coordinates Manipulation with Mathematica
New Physics
http://koix.ksci.re.kr/KISTI1.1003/JNL.JAKO200510103443436
50
(3)
134 -
141

Yun HJ
2006
Mathematica Simulations for the Analysis to the Coriolis Force in a Non-Inertial Frame of Reference
New Physics
http://koix.ksci.re.kr/KISTI1.1003/JNL.JAKO200610103418019
52
(2)
87 -
95

Motion on surface of the Earth: Coriolis Force
Zeleny
http://demonstrations.wolfram.com/MotionOnSurfaceOfTheEarthCoriolisForce

Zimmerman RL
,
Olness FI
1995
Mathematica for Physics
Addison-Wesley
New York
44 -

Citing 'Pedagogical Mathematica Platform Visualizing the Coriolis Effects in 3-Cell Atmospheric Circulation Model
'

@article{ OJOOBS_2014_v31n1_91}
,title={Pedagogical Mathematica Platform Visualizing the Coriolis Effects in 3-Cell Atmospheric Circulation Model}
,volume={1}
, url={http://dx.doi.org/10.5140/JASS.2014.31.1.91}, DOI={10.5140/JASS.2014.31.1.91}
, number= {1}
, journal={Journal of Astronomy and Space Sciences}
, publisher={The Korean Space Science Society}
, author={Kim, Bogyeong
and
Yun, Hee-Joong}
, year={2014}
, month={Mar}