DeepDyve requires Javascript to function. Please enable Javascript on your browser to continue.

Explicit two-dimensional solutions for the ocean flow in arctic gyres
Explicit two-dimensional solutions for the ocean flow in arctic gyres
Haziot, Susanna
2018-06-06 00:00:00
Monatsh Math https://doi.org/10.1007/s00605-018-1198-3 Explicit two-dimensional solutions for the ocean ﬂow in arctic gyres Susanna V. Haziot Received: 7 May 2018 / Accepted: 29 May 2018 © The Author(s) 2018 Abstract We investigate the ocean ﬂow in an arctic gyre. For given linear or constant oceanic vorticites, we prove explicit solutions which depend on both the polar and the azimuthal angles in the spherical coordinates associated with a rotating Earth. This enables us to identify resonant modes. Keywords Geophysical ﬂow · Spherical coordinates · Mercator projection Mathematics Subject Classiﬁcation 35Q86 · 34B05 · 86A05 1 Introduction In an overly-simpliﬁed model of our planet, the major ocean currents would be solely driven by the winds. In reality however, the Coriolis force, due to the rotation of the Earth, deﬂects these practically steady currents so that they produce in the Northern Hemisphere a clockwise circular motion and in the Southern Hemisphere counter- clockwise circular paths. Furthermore, continental boundaries interrupt and divert their ﬂow as well. These large-scale linked-up ocean currents rimmed by land masses are known as “gyres” [2]. There are three types of gyres: tropical, which form near the Equator (conﬁned to either the Northern or the Southern hemisphere since gyres do not cross the Equator, which acts like a wave-guide due to the vanishing of the meridional component of the Coriolis force, thereby facilitating east-west ﬂow propagation— Communicated by A. Constantin. B Susanna V. Haziot susanna.haziot@univie.ac.at Fakultät für Mathematik, Universität Wien, Oskar-Morgenstern-Platz 1, 1090 Wien, Austria 123 S. V. Haziot see the discussion in [7,8]), subtropical (between polar and equatorial regions) and subpolar (in the polar regions). In this paper we will focus on the Arctic basin. It is worth mentioning that the Arctic is very different from the Antarctic: whereas the Arctic is a semi-enclosed ocean almost completely surrounded by landmasses, the Antarctic is almost its complete geographical opposite, a single landmass encircled by a very powerful current known as the Antarctic Circumpolar Current (ACC) [20]; we refer to [9,18,19,22] for studies of the gyre-like ﬂow of the ACC. To study gyres, one can reduce their complexity by using some of their properties: as their vertical velocity is negligeable with respect to their horizontal one, they can be considered as being two-dimensional ﬂows. This feature was exploited recently in [10] to derive a model for ocean gyres in spherical coordinates. In [3], the author uses the stereographical projection to transform this model into a planar elliptic boundary-value problem which he then reduces to a one-dimensional ODE by neglecting azimuthal variations. This permits him to provide explicit solutions for given oceanic vorticities. In this paper we replace the stereographic projection by the Mercator projection, which reduces the model in [10] to a semilinear elliptic PDE that is simpler than the equation obtained recently in [3]. The advantage of this new transformation is that for linear oceanic vorticity functions, thanks to the periodicity of the azimuthal velocity, we can use Fourier series to separate our x and y variables (in a planar geometry context, as we can take advantage of the Mercator projection). This in turn provides us with explicit solutions which take into account possible azimuthal variations. The availability of explicit linear solutions opens up possibilities for the investigation of weakly nonlinear effects by means of perturbation theory, along the lines of the developments in [4,5] of the considerations made in [3]. 2 The Mercator projection We begin by recalling the recently derived model for gyres in spherical coordinates. [10]Let θ ∈[0,π) be the polar angle, where θ = 0 corresponds to the South Pole. Our latitude angle is therefore θ − .Let ϕ ∈[0, 2π) be the azimuthal (or longitude) angle. The polar and azimuthal velocity components of the horizontal ﬂow on the Earth are given by ψ , −ψ , ϕ θ sin(θ ) respectively, where ψ(θ, ϕ) represents the stream function in spherical coordinates. The governing equation for gyres is given by + cot(θ ) + = F ( − ω cos(θ )) ϕϕ θ θθ sin (θ ) where (θ, ϕ) = ψ(θ, ϕ) + ω cos(θ ) is associated with the vorticity of motion of the ocean relative to the Earth’s surface, 2ω cos(θ ) is the spin vorticity due to the rotation of the Earth and F ( − ω cos(θ )) is the oceanic vorticity, due to the motion of the ocean and speciﬁc to a particular gyre. Given ω and F,wehavetosolve thegoverning 123 Explicit two-dimensional solutions for the ocean ﬂow… equation in a given spherical region whose boundary is ψ (ϕ). For this, we turn to the Mercator projection. The Mercator projection maps the sphere onto the plane such that the north-south direction is the horizontal direction, the east-west direction is the vertical direction with the length of the Equator preserved, and all paths of equal compass bearing on the sphere are straight lines. It is conformal but distorts areas [12]. The change of variables x =− ln tan , y = ϕ satisﬁes the above mentioned properties. With this change of variables, the North Pole (θ = π) corresponds to x =−∞ and the Equator (θ = ) to x = 0. Therefore, x < 0 in the Northern Hemisphere with −1 cos(θ ) = tanh(x ), sin(θ ) = cosh (x ) Let us set u(x , y) = ψ(θ, ϕ). We can then rewrite the governing equation as the following semilinear elliptic PDE F (u(x , y)) sinh(x ) u(x , y) = + 2ω (1) 2 3 cosh (x ) cosh (x ) with boundary condition u(x , y) = u (y), x = x < 0 0 0 0 −x if the boundary of our spherical region corresponds to the parallel θ = 2arctan(e ) ∈ ,π , situated in the Northern Hemisphere (Fig. 1). Since u(x , y) is periodic in the y-variable, we can reformulate u as the Fourier series iky u(x , y) = u (x )e , x ≤ x < 0. k 0 k∈Z and therefore, if F is linear, the linear elliptic PDE reduces to an inﬁnite number of boundary-value problems for second-order ODEs: F (u ) u (x ) − k u = cosh (x ) u (x ) = α (2) k 0 k when k ∈ Z and k = 0, and F (u ) sinh(x ) u (x ) = + 2ω 2 3 cosh (x ) cosh (x ) u (x ) = α (3) 0 0 0 123 S. V. Haziot Fig. 1 The Mercator Projection: We project the polar angle onto the y-axis and the azimuthal angle onto the x-axis. The North Pole is situated at x =−∞ where iky u(x , y) = α e . (4) 0 k k∈Z When u(x , y) is a constant the boundary is deﬁned as a streamline (i.e. a level set of ψ). For our general solution to be physically admissible, we require that our stream function be equal to some constant at the North Pole, and that the North Pole be a stagnation point. Therefore, u(x , y) needs to satisfy the following asymptotic condi- tions: lim {u(x , y)}= 0, lim { u (x , y), u (x , y) cosh(x )}= (0, 0) (5) x y x →−∞ x →−∞ uniformly in y ∈ R since the horizontal velocity is ψ , −ψ = (u , u ) ϕ θ y x sin(θ ) cosh(x ). 123 Explicit two-dimensional solutions for the ocean ﬂow… 3 Constant oceanic vorticity Let us consider the case when the oceanic vorticity is constant, F ≡ γ with γ ∈ R. Therefore, (2)–(3) takes the form: u (x ) = k u , k =0(6) γ sinh(x ) u (x ) = + 2ω (7) 2 3 cosh (x ) cosh (x ) with boundary conditions u (x ) = α , k = 0 k 0 k u (x ) = α , k =0(8) 0 0 0 The general solution to (6)is kx −kx u (x ) = c (k)e + c (k)e , k = 0. (9) k 1 2 Since (5) yields 2π −iky lim u (x ) = lim u(x , y)e dy = 0, (10) x →−∞ x →−∞ 2π we must have kx u (x ) = c (k)e , k > 0, k 1 kx u (x ) = c (k)e , k < 0. (11) k 2 Theorem 1 If u(x , y) = ψ for all y ∈[0, 2π), then u ≡ 0 for k = 0. 0 0 k iky Proof We note that u(x , y) = u (x )e for k = 0 is only equal to some 0 k 0 k∈Z constant for all y ∈[0, 2π) when u (x ) = 0 for all k = 0. Because of (11), this k 0 means u ≡ 0 for all k = 0. We are now ready to state and prove our main result. Theorem 2 Assume the boundary is not a streamline, i.e. u (y) is not a constant. If α = α = 0, then for a given γ ∈ R, 1 −1 iky k(x −x ) iky k(x −x ) 0 0 u(x , y) = α e e + α e e k k k∈Z,k≤−2 k∈Z,k≥2 + γ [x + ln(2 cosh(x ))]− ω[1 + tanh(x )] is the general solution of (2)–(3) with F = γ , satisfying the asymptotic conditions (5) and the boundary condition (4), provided α = γ [x + ln(2 cosh(x ))]− ω[1 + 0 0 0 123 S. V. Haziot tanh(x )]. There are no solutions for α = 0 or for α = 0, the modes k =±1 being 0 1 −1 resonant. Proof We begin by looking at the case k = 0. We have: γ sinh(x ) u (x ) = + 2ω , x ≤ x (12) 2 3 cosh (x ) cosh (x ) It can be easily veriﬁed that u (x ) = γ ln(cosh(x )) − ω tanh(x ) x ≤ x P 0 is a particular solution of (12) and therefore, the general solution can be written in the form: u(x ) = δx + β + γ ln(cosh(x )) − ω tanh(x ), x ≤ x ,δ,β ∈ R (13) We claim that for δ = γ and β =−ω+γ ln(2),(13) satisﬁes the asymptotic conditions (5). Indeed, setting δ = γ , we get: lim {u (x ) cosh(x )}= lim γ e − = 0 x →−∞ x →−∞ cosh(x ) For the ﬁrst condition in (5), we have lim {u(x )}= lim {γ [ln(1 + e ) − ln(2)]+ β − ω tanh(x )} x →−∞ x →−∞ =−γ ln(2) + β + ω which vanishes when β = γ ln(2) − ω. Note that these values of δ and β are the only ones that come in question. From the boundary condition (4), we get that α = γ [x + ln(2 cosh(x ))]− ω[1 + tanh(x )] 0 0 0 0 x −x We now look at the case k = 1 and k =−1. In order for u (x ) = c e + c e to 1 1 2 satisfy (5), we need to set c = c = 0. Therefore, for k =±1, u ≡ 0. Consequently, 1 2 k from the boundary condition (4) we can deduce that for α = 0or α = 0, we have 1 −1 no physically relevant solution. From (11), we can see that (5) is satisﬁed for all |k|≥ 2. In addition, using the boundary condition (4), we have iky u(x , y) = α e 0 k k∈Z 123 Explicit two-dimensional solutions for the ocean ﬂow… from which we get −kx c (k) = α e , k ≥ 2 1 k kx c (k) = α e , k ≤−2. 2 k Remark 1 The problem with no y-dependence was considered in [3]. The result in [3] corresponds to the setting u ≡ 0for k = 0. Concerning the physical relevance of the above considerations, the basic sources of oceanic vorticity are wind force and the gravitational forces due to the relative motions of the Moon, the Sun and the Earth in the form of the ﬂood and ebb tidal currents (see the discussions in [6,16]). Both oceanic vorticities can be realisticallly regarded as non-zero constants (see [11,14]), with the sign depending on the prevalent wind direction, and, respectively, on whether the ﬂow is of ebb or ﬂood type. Let us note that wave-current interactions in ﬂows with vorticity are a topic of great current interest (see the discussions in [6,11,13,15]), at the large scales that are relevant for gyre ﬂows the presence of surface waves is not relevant. 4 Some special linear oceanic vorticities Let us now consider the case of linear vorticity, setting F (u ) = au .(2) then takes k k the form: u (x ) = + k u (x ) (14) cosh (x ) for k = 0 and a sinh(x ) u (x ) = u (x ) + 2ω (15) 2 3 cosh (x ) cosh (x ) for k = 0. The case for k = 0 has been shown in [3]. To begin with, we make the change of variable z = tanh(x ) 2 −2 and thus, using the fact that 1 − z = cosh (x ),(14) becomes (1 − z )u (z) − 2zu (z) + −a − u (z) = 0. (16) k k 1 − z We rewrite √ √ 1 − 4a 1 −1 + 1 − 4a 1 + 1 − 4a − a = − = = l(l + 1). (17) 4 4 2 2 123 S. V. Haziot Setting (17)into(16) yields (1 − z )u (z) − 2zu (z) + l(l + 1) − u (z) = 0 (18) k k 1 − z which, if l is a positive integer, is known as the associated Legendre Differential equation, see [1,17]. Therefore, the general solution to (18)is x 2 cosh (t ) k k u (x ) = c P (tanh(x )) + c P (tanh(x )) dt (19) k 1 2 l l [P (tanh(t ))] −∞ where P (z) denotes the associated Legendre function. 1−(2l+1) Theorem 3 For a = , where l is any positive integer, and for any b ∈ R, k −1 k iky u(x , y) = α [P (tanh(x ))] P (tanh(x ))e k 0 l l k∈Z , |k|≤l + b[x + ln(2 cosh(x ))]− ω[1 + tanh(x )] with −1 + 1 − 4a l = is the general solution to (2) with F (u) = au + b satisfying the asymptotic conditions (5) and the boundary condition (4), provided α = b[x + ln(2 cosh(x ))]− ω[1 + 0 0 0 tanh(x )]. Proof To begin with, the associated Legendre polynomial for 0 ≤ k ≤ l is deﬁned as: k k+l (−1) 1 d 1 P (tanh(x )) = − (20) l 2 k+l 2 2 l! dx cosh (x ) cosh (x ) and we can therefore write lim P (tanh(x )) = lim (21) x →−∞ x →−∞ cosh (x ) where n ≥ 2, n ∈ N. Therefore, for (20) to satisfy the asymptotic conditions (5), we need to set c = 0. From the boundary condition (4), we have c = P (tanh(x )) thus concluding the proof for k > 0. 123 Explicit two-dimensional solutions for the ocean ﬂow… For negative k, the associated Legendre polynomial P (z) is deﬁned as (l − k)! −k k k P (z) = (−1) P (z) (l + k)! and therefore the proof follows directly from the one for k > 0. Finally, from the previous section, we know that the asymptotic conditions are satisﬁed for b sinh(x ) u (x ) = + 2ω 2 3 cosh (x ) cosh (x ) if and only if u (x ) = b[x + ln(2 cosh(x ))]− ω[1 + tanh(x )]. From the boundary condition (4), we have α = b[x + ln(2 cosh(x ))]− ω[1 + tanh(x )]. 0 0 0 0 5 Some considerations about general linear oceanic vorticity It is interesting to observe that (14) can be reduced to a Riccati’s equation, which in turn is equivalent to solving an Ermakov equation. f (x ) Indeed, assuming that the solution to (14)isofthe form u (x ) = e , where f is an arbitrary function, solving (14) is equivalent to solving the following Riccati equation (see [17]) 2 2 g (x ) =−g (x ) + + k (22) cosh (x ) with g(x ) = f (x ). Using the following substitution g(x ) = s(x )v(x ) (23) and choosing s(x ) such that s (x ) 3(s (x )) a 2 2 − + s (x ) = + k . (24) 2 2 2s(x ) 4s (x ) cosh (x ) we can write (22)as 2 2 v (x ) =−s(x )v (x ) − 2s(x )h(x )v(x ) − s(x )h (x ) + s(x ) − h (x ) (25) where s (x ) h(x ) = (26) 2s (x ) We assume for now that such an s exists. A condition for its existence will be given later on. 123 S. V. Haziot We can now rewrite (25)as v (x ) =−s(x )[(v(x ) + h(x )) − 1]− h (x ) (27) and using the substitution w(x ) = v(x ) + h(x ) we get w (x ) =−s(x )[w (x ) − 1]. (28) Thesolutionto(28)is w(x ) = tanh s(t ) dt −∞ and therefore a particular solution to (14)is cosh s(t ) dt −∞ u = √ . (29) s(x ) It remains to ﬁnd s(x ).From (24), we get the following second-order ODE 3 a 2 4 2 2 s (x )s(x ) − [s (x )] + 2s (x ) = 2s (x ) + k (30) cosh (x ) and using the substitution −1/2 γ(x ) = s (x ) we can reduce (30) to the following Ermakov equation ( see [21]) −3 2 γ (x ) = γ (x ) + γ(x ) + k (31) cosh (x ) The general solution to an Ermakov equation (see [21]) is of the form dt 2 2 2 k γ (x ) = p (x ) + p (x ) k + k (32) 1 2 1 0 0 p (t ) dt −∞ where p (x ) is a particular solution to γ˜ (x ) =˜ γ(x ) + k . (33) cosh (x ) Remark 2 If a =−l(l + 1), for some integer l ≥ 1, then from (19) we already know that a particular solution to (33) is the associated Legendre polynomial P (z) with 123 Explicit two-dimensional solutions for the ocean ﬂow… −1+ 1−4a z = tanh(x ) and l = . We can therefore conclude that s(x ) = . k 1 [P (z)] 1 + k + k 2 1 l −∞ 2 [P (z)] Acknowledgements Open access funding provided by University of Vienna. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 Interna- tional License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. References 1. Abramowitz, M., Stegun, I.: Handbook of Mathematical Functions with Formulas, Graphs and Math- ematical Tables. Dover Publications, Inc., New York (1965) 2. Apel, J.: Principles of Ocean Physics. Academic Press, London (1987) 3. Chu, J.: On a differential equation arising in geophysics. Monatsh. Math. (2017). https://doi.org/10. 1007/s00605-017-1087-1 4. Chu, J.: On a nonlinear integral equation for the ocean ﬂow in arctic gyres. Q. Appl. Math. (2017). https://doi.org/10.1090/qam/1486 5. Chu, J.: Monotone solutions of a nonlinear differential equation for geophysical ﬂuid ﬂows. Nonlinear Anal. 166, 144–153 (2018) 6. Constantin, A.: Nonlinear water waves with applications to wave-current interactions and tsunamis. In: CBMS-NSF Regional Conference Series in Applied Mathematics, vol. 81. SIAM, Philadelphia, PA (2011) 7. Constantin, A., Johnson, R.S.: The dynamics of waves interacting with the Equatorial Undercurrent. Geophys. Astrophys. Fluid Dyn. 109, 311–358 (2015) 8. Constantin, A., Johnson, R.S.: An exact, steady, purely azimuthal equatorial ﬂow with a free surface. J. Phys. Oceanogr. 46, 1935–1945 (2016) 9. Constantin, A., Johnson, R.S.: An exact, steady, purely azimuthal ﬂow as a model for the Antarctic Circumpolar Current. J. Phys. Oceanogr. 46, 3585–3594 (2016) 10. Constantin, A., Johnson, R.S.: Large gyres as a shallow-water asymptotic solution of Euler’s equation in spherical coordinates. Proc. R. Soc. Lond. A 473, 20170063 (2017) 11. Constantin, A., Strauss, W., Varvaruca, E.: Global bifurcation of steady gravity water waves with critical layers. Acta Math. 217, 195–262 (2016) 12. Daners, D.: The Mercator and stereographic projections, and many in between. Am. Math. Mon. 119, 199–210 (2012) 13. da Silva, A.F.T., Peregrine, D.H.: Steep, steady surface waves on water of ﬁnite depth with constant vorticity. J. Fluid Mech. 195, 281–302 (1988) 14. Ewing, J.A.: Wind, wave and current data for the design of ships and offshore structures. Mar. Struct. 3, 421–459 (1990) 15. Henry, D.: Large amplitude steady periodic waves for ﬁxed-depth rotational ﬂows. Commun. Partial Differ. Equ. 38, 1015–1037 (2013) 16. Jonsson, I.G.: Wave-current interactions. In: Le Méhauté, B., Hanes, D.M. (eds.) The Sea: Ocean Engineering Science, vol. 9(A), Wiley, pp. 65–120 (1990) 17. Kamke, E.: Differentialgleichungen: Lösungsmethoden und Lösungen. Akademische Verlagsge- sellschaft, Leipzig (1967) 18. Marynets, K.: A nonlinear two-point boundary-value problem in geophysics. Monatsh. Math. (2017). https://doi.org/10.1007/s00605-017-1127-x 19. Marynets, K.: A weighted Sturm-Liouville problem related to ocean ﬂows. J. Math. Fluid Mech. (2017). https://doi.org/10.1007/s00021-017-0347-0 123 S. V. Haziot 20. National Snow and Ice Data Center. All About Sea Ice. https://nsidc.org/cryosphere/seaice/index.html. Accessed 24 Nov 2017 21. Polyanin, A.D., Zaitsev, V.F.: Handbook of Exact Solutions for Ordinary Differential Equations, 2nd edn. Chapman and Hall/CRC, Boca Raton (2003) 22. Quirchmayr, R.: A steady, purely azimuthal ﬂow model for the Antarctic Circumpolar Current. Monatsh. Math. (2017). https://doi.org/10.1007/s00605-017-1097-z
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png
Monatshefte f�r Mathematik
Springer Journals
http://www.deepdyve.com/lp/springer-journals/explicit-two-dimensional-solutions-for-the-ocean-flow-in-arctic-gyres-TFjWOMP7Dd