Pedagogical Mathematica Platform Visualizing the Coriolis Effects in 3Cell Atmospheric Circulation Model
 Author: Kim Bogyeong, Yun HeeJoong
 Organization: Kim Bogyeong; Yun HeeJoong
 Publish: Journal of Astronomy and Space Sciences Volume 31, Issue1, p91~99, 15 March 2014

ABSTRACT
The atmospheric flow in the 3Cell model of global atmosphere circulation is described by the Lagrange's equation of the noninertial 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 noninertial 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.

KEYWORD
Mathematica , fictitious forces , Coriolis force , noninertial frame

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 noninertial frame of reference according to the 3Cell general circulation of the atmosphere (GCA) model rather than Hadley Cell model (Ahrens 2001). Longterm recording data from the satellites approve of the 3Cell GCA model (“AMNHWeather 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 largescale 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 GoRound table (“MerryGoRound” 2014) or BathTub 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,
Mathematica is helpful for the convenient function of symbolic calculations and graphic manipulations. TheMathematica simulation of the atmospheric circulation matching the 3Cell GCA model has been presented in our previous work (Yun 2006), however, the 3Cell GCA model was shown with the illustration rather than with active platform. In this paper, we present anewMathematica platform presenting the GCA simulation according to the 3Cell 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.2. LAGRANGE’S EQUATION IN A NONINERTIAL 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 noninertial 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_{o} for the Lagrange's equation in the noninertial 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 noninertial frame is well described in other books (Symon 1971. Landau & Lifshitz 1976).
3.
MATHEMATICA PROGRAMING FOR GERNERAL CIRCULATION OF THE ATMOSPHERE3.1 Vector calculation of the effective force
To analyze the deflection effect of the atmospheric flow based on the 3Cell GCA model at any point on the Earth's surface, we create a position function
point[t, θ] in solving the Eq. (4) inMathematica (Zimmerman & Olness 1995, Tam 1997). If a latitude (λ = π /2 θ) is given, the wind vector function [θ] is determined by the position functionpoint[t, θ] . We assume that the initial wind blows along the meridian and its velocity is determined according to the 3Cell GCA model.Mathematica coding for solving Eq. (4) is shown in In[1]  In[10] and itsMathematica solution is Out[11]:point[t, θ] ={t v0x, t^2ω0x Cos[θ], 1/2 gt^2} while n Order = 0 for simplicity. Oncepoint[t, θ] is given, theParametric Plot draw the path of the wind with a time domain vector array of thepoint[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 3Cell 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 noninertial 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.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 thevector product ortriple vector product is coding as cross[a, b] or cross[a, cross[b, c]] inMathematica 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[3133]), 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, AMNHWeather and Climate Events 2014).3.2
Mathematica platform for the 3cell atmospheric general circulation modelUnderstanding 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 3Cell GCA model. The fundamental concept of the inertial force effects in the noninertial 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 3Cell 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
Parametric Plot of Mathematica. On starting the program, platform will show Fig. 3a. Simulation performs the builtin program when you click the ► appearing while you spread the ⊕ oft panel of the platform. While the program is executing, you can change the parameters ofWind Speed, Earth's Rotation , andTheta , 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 3Cell 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 2dimensional deflections but also 3dimensional 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 functionpoint[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). InMathematica , we can save the snapshots of the simulations and print it. This program also executes well on later version of 8.0Mathematica .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): 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 noninertial 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 longrange 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 3Cell GCA model in
Mathematica. The platform draws the realistic path of the wind, throughpoint[t, θ] with the time domain of the vector array of the solution of Eq. (4) using theParametric Plot inMathematica. We expect this platform will be a helpful tool for the physicist and the scientist to analyze the atmospheric general circulation dynamics.

[]

[]

[]

[]

[]

[]

[]

[FIG. 1.] 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.

[FIG. 2.] The table of the deflections at eight points on the globe, The points are selected In the different atmospheric zone of the 3Cell GCA model.

[]

[]

[]

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

[FIG. 3.] Snapshots of the Mathematica platform of the general atmospheric circulation to the 3Cell 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.

[FIG. 4.] 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 3Cell model from the North pole: θ = 1/8 π: NorthPolar Easterlies, 2/8 π: SouthWest Westerlies, 3/8 π: NorthEast Trade Winds, 5/8 π: SouthEast Trade Winds, 6/8 π: NorthWest Westerlies, 7/8 π: SouthPolar Easterlies in both columns.

[FIG. 5.] Wind deflections at the equator in 2D (left column) graphics and 3D graphics (right column) for the different wind velocities.