arxiv:math/ v1 [math.ds] 22 Jan 1999

Transcription

1 Escape Probabilit, Mean Residence Time and Geophsical Fluid Particle Dnamics arxiv:math/9999v [math.ds] 22 Jan 999 James R. Brannan, Jinqiao Duan and Vincent J. Ervin Clemson Universit Department of Mathematical Sciences Clemson, South Carolina 29634, USA. December 25, 998 (Final Version) Abstract Stochastic dnamical sstems arise as models for fluid particle motion in geophsical flows with random velocit fields. Escape probabilit (from a fluid domain) and mean residence time (in a fluid domain) quantif fluid transport between flow regimes of different characteristic motion. We consider a quasigeostrophic meandering jet model with random perturbations. This jet is parameterized b the parameter β = 2Ω cos(θ), where Ω is the rotation rate of the earth, r the earth s r radius and θ the latitude. Note that Ω and r are fied, so β is a monotonic decreasing function of the latitude. The unperturbed jet (for < β < 2 ) consists of a basic flow with attached eddies. With 3 random perturbations, there is fluid echange between regimes of different characteristic motion. We quantif the echange b escape probabilit and mean residence time. For an edd, the average escape probabilit for fluid particles (initiall inside the edd) escape into the eterior retrograde region is smaller than escape into the jet core for < β <.3333, while for.3333 < β < 2, the opposite holds. 3 For a unit jet core near the jet troughs, the average escape probabilit for fluid particles (initiall inside the jet core) escape into the northern recirculating region is greater than escape into the southern recirculating region for < β <.5, while for.385 < β < 2, the opposite holds. Moreover, for 3.5 < β <.385, fluid particles are about equall likel to escape into either recirculating regions. Furthermore, for a unit jet core near the jet crests, the situation is the opposite as for near the jet troughs. The maimal mean residence time of fluid particles initiall in an edd increases as β increases from to.432 (or as latitude decreases accordingl), then decreases as β increases from.432 to 2 (or as 3 latitude decreases accordingl). However, the maimal mean residence time of fluid particles initiall in a unit jet core alwas increases as β increases (or as latitude decreases). Ke Words: stochastic dnamics, transport and miing, geophsical flows, escape probabilit, mean residence time 273 Author for correspondence Dr. Jinqiao Duan, duan@math.clemson.edu Fa: (864) Tel: (864)656-

2 Stochastic dnamics: Escape probabilit and mean residence time Stochastic dnamical sstems are used as models for various scientific and engineering problems. We consider the following class of stochastic dnamical sstems ẋ = a (, ) + b (, )ẇ, () ẏ = a 2 (, ) + b 2 (, )ẇ 2, (2) where w (t), w 2 (t) are two real independent Wiener processes (white noises) and a, a 2, b, b 2 are given deterministic functions. More complicated stochastic sstems also occur in applications ([], [25], []), [], [34], [7]). For a planar bounded domain D, we can consider the eit problem of random solution trajectories of ()-(2) from D. To this end, let D denote the boundar of D and let Γ be a part of the boundar D. The escape probabilit p(, ) is the probabilit that the trajector of a particle starting at (, ) in D first hits D (or escapes from D) at some point in Γ, and p(, ) is known to satisf ([9], [32] and references therein) 2 b2 (, )p + 2 b2 2(, )p + a (, )p + a 2 (, )p =, (3) p Γ =, (4) p D Γ =. (5) Suppose that initial conditions (or initial particles) are uniforml distributed over D. The average escape probabilit P that a trajector will leave D along the subboundar Γ, before leaving the rest of the boundar, is given b (e.g., [9], [32]) P = p(, )dd, (6) D D where D is the area of domain D. The residence time of a particle initiall at (, ) inside D is the time until the particle first hits D (or escapes from D). The mean residence time u(, ) is given b (e.g., [32], [24], [28] and references therein) 2 b2 (, )u + 2 b2 2(, )u + a (, )u + a 2 (, )u =, (7) 2 A quasigeostrophic jet model u D =. (8) The Lagrangian view of fluid motion is particularl important in geophsical flows since onl Lagrangian data can be obtained in man situations. It is essential to understand fluid particle trajectories in man fluid problems. Escape probabilit (from a fluid domain) and mean residence time (in a fluid domain) quantif fluid transport between flow regimes of different characteristic motion. Deterministic quantities 2

3 like escape probabilit and mean residence time can be computed b solving Fokker-Planck tpe partial differential equations. We now use the these ideas in the investigation of meandering oceanic jets. Meandering oceanic jets such as the Gulf Stream are strong currents dividing different bodies of water. Recentl, del-castillo-negrete and Morrison ([6]), and Pratt et al. ([27], [22]) have studied models for oceanic jets. These models are dnamicall consistent to within a linear approimation, i.e., the potential vorticit is approimatel conserved. Del-castillo-Negrete and Morrison s model consists of the basic flow plus time-periodic linear neutral modes. In this paper, we consider an oceanic jet consisting of the basic flow as in del-castillo-negrete and Morrison ([6]), plus random-in-time noise. This model incorporates small-scale oceanic motions such as the molecular diffusion ([6]), which is an important factor in the Gulf Stream ([9], [2]). The irregularit of RAFOS floats ([3], [33], [2]) also suggests the inclusion of random effects in Gulf Stream modeling. This random jet ma also be viewed as satisfing, approimatel in the spirit of del-castillo-negrete and Morrison ([6]), the randoml wind forced quasigeostrophic model. Several authors have considered the randoml wind forced quasigeostrophic model in order to incorporate the impact of uncertain geophsical forces ([3], [3], [5], [23], [7], [4]). The studied statistical issues such as estimating correlation coefficients for the linearized quasigeostrophic equation with random forcing. There is also recent work which investigates the impact of the uncertaint of the ocean bottom topograph on quasigeostrophic dnamics ([8]). The randoml forced quasigeostrophic equation takes the form ([23]) ψ t + J(ψ, ψ) + βψ = dw dt, (9) where W(,, t) is a space-time Wiener process (white noise). The stream function would have a random or noise component ([2], [9]). Note that β is the meridional derivative of the Coriolis parameter [26] [5], i.e., β = 2Ω r cos(θ), where Ω is the rotation rate of the earth, r the earth s radius and θ the latitude. Since Ω and r are fied, β is a monotonic decreasing function of the latitude. The deterministic meandering jet derived in del-castillo-negrete and Morrison ([6]) is Ψ(, ) = tanh() + a sech 2 ()cos(k) + c, where a =., c = 3 ( + 2( 3 2 β), k = β), β 2 3. This Ψ(, ) is an approimate solution of the usual quasigeostrophic model ψ t + J(ψ, ψ) + βψ =. () With random wind forcing or molecular diffusive forcing in the stochastic quasigeostrophic model (9), the stream function would have a random or noise component. We approimate this noise component b adding a noise term to the above deterministic stream function Ψ(, ), that is, in the rest of this paper, we consider the following random stream function as a model for a quasigeostrophic meandering jet, Ψ(, ) = tanh() + a sech 2 ()cos(k) + c + noise. 3

4 Figure : Unperturbed jet: ǫ = and β = /3. The equations of motion for fluid particles in this jet then have noise terms. We further approimate them as white noises (or Wiener processes) or more specificall, d = Ψ dt + ǫdw, () d = Ψ dt + ǫdw 2, (2) d = [sech 2 () + 2a sech 2 ()tanh()cos(k) c]dt + ǫdw, (3) d = ak sech 2 ()sin(k)dt + ǫdw 2, (4) where < ǫ <, and w (t), w 2 (t) are two real independent Wiener processes (in time onl). The calculations below are for ǫ =.. Note that β is now the onl parameter 3, 4, as a is given and c, and k depend onl on β [, 2 3 ]. When ǫ =, the deterministic jet consists of the jet core and two rows of recirculating eddies, which are called the northern and southern recirculating regions. Outside the recirculating regions are the eterior retrograde regions; see Figure. 3 Escape probabilit We take D to be either an edd or a piece of jet core (see Figures 2, 3). This piece of jet core has the same horizontal length scale as an edd, and it is one period of the deterministic jet core (note that the deterministic velocit field is periodic in ). Thus we call this piece of jet core a unit jet core. 4

5 Figure 2: An edd: β = / Figure 3: A unit jet core near a trough: β = /3 5

6 .8.6 p Figure 4: Escape probabilit of fluid particles (initiall in an edd in Figure 2) eiting into the eterior retrograde region: β = /3 From (3), (4), (5), the escape probabilit p(, ) that a fluid particle, initiall at (, ), crosses the subboundar Γ of the domain D satisfies ǫ p + [sech 2 () + 2a sech 2 ()tanh()cos(k) c]p ak sech 2 ()sin(k)p =, (5) p Γ =, (6) p D Γ =. (7) We take Γ to be either top or bottom boundar of an edd or a unit jet core (see Figures 2, 3). We numericall solve this elliptic sstem for various values of β between and 2 3. In the unit jet core case, we take periodic boundar condition in horizontal (meridional) direction, with period 2π k. A piecewise linear, finite element approimation scheme was used for the numerical solutions of the escape probabilit p(, ), and the mean residence time u(, ), described b the elliptic equations (3), and (7), respectivel. Using a collection of points ling on the boundar, piecewise cubic splines were constructed to define the boundar of the edd, and the top and bottom boundaries of the jet core. Computational (triangular) grids for the edd and the jet core were then obtained b deforming regular grids constructed for an ellipse and a rectangular region, respectivel. The computed escape probabilit crossing the upper or lower boundar of an edd or a unit jet core, for the case of β = 3, are shown in Figures 4, 5, 6 and 7. Suppose that the fluid particles are initiall uniforml distributed in D (an edd or a unit jet core). We can also compute the average escape probabilit P that a particle will leave D along the upper or lower subboundar Γ, using the formula (6); see Figures 8 and 9. For an edd, the average escape probabilit for fluid particles (initiall inside the edd) escape into the eterior retrograde region is smaller than escape into the jet core for < β <.3333, while for 6

7 .8.6 p Figure 5: Escape probabilit of fluid particles (initiall in an edd in Figure 2) eiting into the jet core: β = /3.8.6 p Figure 6: Escape probabilit of fluid particles (initiall in the unit jet core in Figure 3) eiting into the northern recirculating region: β = /3 7

8 p Figure 7: Escape probabilit of fluid particles (initiall in the unit jet core in Figure 3) eiting into the southern recirculating region: β = / < β < 2 3, the opposite holds (Figure 8). Thus β =.3333 is a bifurcation point. Also, the average escape probabilit for fluid particles escape into the eterior retrograde region increases as β increases from to.54 (or as latitude decreases accordingl), and then decreases as β increases from.54 to 2 3 (or as latitude decreases accordingl). Thus β =.54 is another bifurcation point. The opposite holds for the average escape probabilit for fluid particles escape into the jet core. For a unit jet core near the jet troughs, the average escape probabilit for fluid particles (initiall inside the jet core) escape into the northern recirculating region is greater than escape into the southern recirculating region for < β <.5, while for.385 < β < 2 3, the opposite holds. Moreover, for.5 < β <.385, fluid particles are about equall likel to escape into either recirculating regions (Figure 9). Furthermore, for a unit jet core near the jet crests, the situation is the opposite as for near the jet troughs. 4 Mean residence time The mean residence time u(, ) of a fluid particle, initiall at (, ) in either an edd or a piece of jet core (see Figures 2, 3), satisfies ǫ u + [sech 2 () + 2a sech 2 ()tanh()cos(k) c]u ak sech 2 ()sin(k)u =, (8) u D = (9) The mean residence times of fluid particles in an edd or a unit jet core are shown, for the case of β = 3, in Figures,. 8

9 P β Figure 8: Average escape probabilit for a fluid particle in an edd: for eiting into the eterior retrograde region; for eiting into the jet core region P β Figure 9: Average escape probabilit for a fluid particle in the unit jet core near a trough: for eiting into the northern recirculating region; for eiting into the southern recirculating region. 9

10 4 2 8 u Figure : Mean residence time in an edd: β = u Figure : Mean residence time in a unit jet core: β = 3

11 u ma β Figure 2: Maimal value of mean residence time, as a function of β, in an edd. The maimal mean residence time of fluid particles initiall in an edd increases as β increases from to.432 (or as latitude decreases accordingl), then decreases as β increases from.432 to 2 3 (or as latitude decreases accordingl); see Figure 2. However, the maimal mean residence time of fluid particles initiall in a unit jet core alwas increases as β increases (or as latitude decreases accordingl); see Figure 3. 5 Discussions The present work on fluid particle motion in random flows takes into account of fluid particle diffusive as well as advective motion. There has been recent work on fluid particle advective motion (molecular diffusion ignored) in time-periodic, quasi-periodic and aperiodic flows (periodic quasi-periodic aperiodic random); see, for eample, [29], [35], [3], [2], [8], [2] and [22]. Our work on random particle motion does not require that the random part is small. However it does require that the flow can be decomposed into stead or unstead deterministic (drift) and random (diffusion) parts; Otherwise the Fokker-Planck formalism does not hold. Acknowledgement. J. Duan would like to acknowledge the hospitalit of the Center for Nonlinear Studies, Los Alamos National Laborator, the Institute for Mathematics and its Applications (IMA), Minnesota, and the Woods Hole Oceanographic Institution, Massachusetts. He is also grateful for the discussions on this work with Diego del-castillo-negrete (Scripps Institution of Oceanograph), Greg Hollowa (Institute of Ocean Sciences, Canada), Julian Hunt (Cambridge Universit), Peter Kiessler (Clemson Universit), Pat Miller (Stevens Institute of Technolog), Larr Pratt (Woods Hole Oceanographic Institution), and Roger Samelson (Oregan State Universit). This work was supported b the National Science Foundation Grant DMS

12 u ma β References Figure 3: Maimal value of mean residence time, as a function of β, in a unit jet core. [] L. Arnold, Stochastic differential equations, John Wile & Sons, New York, 974. [2] D. Beigie, A. Leonard and S. Wiggins, Chaotic Transport in the Homoclini c and Heteroclinic Tangle Regions of Quasiperiodicall Forced Two Dimensional Dnamical Sstems, Non linearit, 4 (99), [3] A. S. Bower and T. Rossb, Evidence of cross-frontal echange processes in the Gulf Stream based on isopcnal RAFOS float data, J. Phs. Oceanograph 9 (989), [4] J. Brannan, J. Duan and T. Wanner, Dissipative quasigeostrophic dnamics under random forcing, J. Math. Anal. Appl. 228 (998), [5] B. Cushman-Roisin, Introduction to Geophsical Fluid Dnamics, Prentice Hall, 994. [6] D. del-castillo-negrete and P. J. Morrison, Chaotic transport b Rossb waves in shear flow, Phs. Fluids A 5 (993), [7] T. DelSole and B. F. Farrell, A stochasticall ecited linear sstem as a model for quasigeostrophic turbulence: Analtic results for one- and two-laer fluids, J. Atmos. Sci. 52 (995) [8] J. Duan and S. Wiggins, Fluid echange across a meandering jet with quasiperiodic variabilit, J. Phs. Oceanograph 26 (996), [9] S. Dutkiewicz, A. Griffa and D. B. Olson, Particle diffusion in a meandering jet, J. Geophs. Res. 98 (993), C9,

13 [] M. I. Freidlin and A. D. Wentzell, Random Perturbations of Dnamical Sstems, Springer-Verlag, New York, 984. [] M. I. Freidlin and A. D. Wentzell, Random Perturbations of Hamiltonian Sstems, American Mathematical Societ, Providence, 994. [2] A. Griffa, Applications of stochastic particle models to oceanographic problems, In Stochastic Modelling and Phsical Oceanograph, R. Adler, P. Müller and B. Rozovskii eds., Birkhauser, 996. [3] A. Griffa and S. Castellari, Nonlinear general circulation of an ocean model driven b wind with a stochastic component, J. Marine Research 49 (99), [4] D. Gurarie and V. Klatskin, Turbulent transport of passive tracers and the onset of diffusivit, In Two-dimensional turbulence in plasmas and fluids, AIP Conference Proceedings, , 997. [5] G. Hollowa, Ocean circulation: Flow in probabilit under statistical dnamical forcing, In Stochastic Models in Geosstems, S. Molchanov and W. Wocznski (eds.), Springer-Verlag, New York, 996. [6] M. B. Isichenko, Percolation, statistical topograph and transport in random medis, Rev. Modern Phs. 64 (992), [7] P. E. Kloeden and E. Platen, Numerical solution of stochastic differential equations, Springer-Verlag, 992; second corrected printing 995. [8] V. Klatskin and D. Gurarie, Random topograph in geophsical models, In Stochastic Models in Geosstems, S. Molchanov and W. Wocznski (eds.), Springer-Verlag, New York, 996. [9] C. C. Lin and L. A. Segel, Mathematics Applied to Deterministic Problems in the Natural Sciences, SIAM, Philadelphia, 988. [2] M. S. Lozier, L. J. Pratt, A. M. Rogerson and P. Miller, Echange geometr revealed b float trajectories in the Gulf Stream, J. Phs. Oceanograph 27 (997), [2] N. Malhotra and S. Wiggins, Geometric Structures, Lobe Dnamics, and Lagrangian Transport in Flows with Aperiodic Time-Dependence, with Applications to Rossb Wave Flow, J. Nonlinear Sci. 8 (998), [22] P. Miller, C. K. R. T. Jones, G. Haller and L. Pratt, Chaotic miing across oceanic jets, In Chaotic, Fractal and Nonlinear Signal Processing, R. Katz (eds.), AIP Conference Proceedings 375, 59-64, 996. [23] P. Müller, Stochastic forcing of quasi-geostrophic eddies, in Stochastic Modelling in Phsical Oceanograph, R. J. Adler, P. Müller and B. Rozovskii (eds.), Birkhäuser, 996. [24] T. Naeh, M. M. Klosek, B. J. Matkowsk and Z. Schuss, A direct approach to the eit problem, SIAM J. Appl. Math. 5 (99), [25] B. Oksendal, Stochastic Differential Equations: An Introduction with Applications, Springer-Verlag, New York, 4th ed., 995. [26] J. Pedlosk, Geophsical Fluid Dnamics, Springer-Verlag, 2nd edition,

14 [27] L. J. Pratt, M. S. Lozier and N. Beliakova, Parcel trajectories in quasigeostrophic jets: Neutral Modes, J. Phs. Oceanograph 25 (995), [28] H. Risken, The Fokker-Planck Equation, Springer-Verlag, New York, 984. [29] V. Rom-Kedar, A. Leonard and S. Wiggins, An analtical stud of transport, miing, and chaos in an unstead vortical flow, F. Fluid Mech. 24 (99), [3] R. M. Samelson, Stochasticall forced current fluctuations in vertical shear and over topograph, J. Geophs. Res. 94 (989), [3] R. M. Samelson, Fluid echange across a meandering jet, J. Phs. Oceanogr., 22(4) (992), [32] Z. Schuss, Theor and Applications of Stochastic Differential Equations, Wile & Sons, New York, 98. [33] T. Song, T. Rossb and E. Carter, Lagrangian studies of fluid echange between the Gulf Stream and surrounding waters, J. Phs. Oceanograph 25 (995), [34] E. Stone and P. Holmes, Random perturbations of heteroclinic attractors, SIAM J. Appl. Math. 5 (99), [35] S. Wiggins, Chaotic Transport in Dnamical Sstems. Springer-Verlag: New York, Heidelberg, Berlin,