Research Plan: Transport and Mixing in Aperiodic Systems
Background and Significance
Over the past 10 years, much work has been done in applying the approach and methods of dynamical theory to the study of transport and mixing in fluids. For immersed point mass particles, the phase space is the physical space in which the fluid flow takes place. Lagrangian Coherent Structures (LCS's), such as KAM tori and the invariant manifolds of hyperbolic fixed points and periodic orbits divide the phase space of these systems into regions of qualitatively different dynamics [Yuan et al., 2004; Babiano et al., 1994].
The LCS's and their interactions govern the entanglement of the trajectories and generate chaotic mixing and transport phenomena [Haller and Poje, 1998; Wiggins, 1992].
In recent years, Lagrangian transport in the context of two-dimensional, time-periodic flows has been extensively developed. For example, the persistence of invariant curves in the Poincar ´e map (KAM curves) gives rise to barriers in the flow; Chaos and Smale horseshoes are mechanisms for the "randomization" of fluid particle trajectories [Babiano et al., 1994]; Melnikov's method provides flux estimates as well as a description of the parameter regimes where chaotic particle motion occurs [Melnikov, 1963; Guckenheimer and Holmes, 1983]. Lobe dynamics enables efficient computations of transport between qualitatively different flow regimes [Rom-Kedar, 1990; Rom-Kedar and Wiggins, 1991]. The central theme of my research program is the development and the application of this framework for time-chaotic systems [Haller, 2000; Haller, 2001; Ide et al., 2002]. This includes a major effort to study the properties of "moving" LCS's and investigate how deterministic chaos is generated in aperiodic flows.
There are a variety of important applied problems that can be addressed with this methodology. These include space missions, molecular fragmentation, and turbulence in fluids. In recent years, numerous experiments have revealed the presence of LCS's in geophysical fluid dynamics settings [Pierrehumbert, 1991; Samelson, 1992; Voth et al., 2002]. In coastal areas and for global ocean circulation, the identification of dynamical barriers and alleyways provides a better understanding of the mixing processes [Coulliette and Wiggins, 2001; Yuan et al., 2002]. Such studies also permits various optimization problems such as reducing the impact of coastal pollutants, finding energy-efficient routes, or devising adaptive sampling strategies in the ocean.
The following sections detail the challenges and approaches for each aspect of my research program: theoretical developments, numerical implementations, and applications.
A Observation and Control of Exponential Dichotomies

Figure 1. An analytical system corresonding to two
rotating eigenvectors (eigenvalues: +1 and 1). For
slow rotation rate (left), a parcel of particles (green
dots) stretches along a moving unstable manifold
(red line). For fast rotations (right), there is not any
exponential dichotomy or unstable manifold.
Figure 2. LCS's (red and blue curves) during the
break-up of the Antarctic polar vortex. Particles
(green parcels) reveal the hyperbolic structure at the
pole and stretch along the attracting LCS.
Autonomous and periodic systems have hyperbolic fixed points (or periodic orbits) and associated stable and unstable manifolds [ Hirsch and Smale, 1974]. An equivalent notion for aperiodic systems are trajectories and invariant manifolds with exponential dichotomies [ Coppel, 1978]. In the vicinity of these structures, particles experience exponential stretching and exponential shrinking [ Ide et al., 2002]. The stretching causes sensitivity to initial conditions, i.e., the exponential separation of two particles initially close to each other. The shrinking is responsible for the folding of the invariant manifolds in a compact region and the existence of a compact chaotic invariant set [ Smale, 1980; Wiggins, 1992]. This description of chaotic transport is based on the observation of particle trajectories (i.e., the Lagrangian dynamics) near the exponential dichotomies [ Ide et al., 2002]. My research plan aims at understanding, detecting and controlling these structures.
As shown in Fig. 1, the relationship between Lagrangian dynamics and observed Eulerian features, i.e., the "frozen"
velocity field, is not straightforward for time-varying systems. The objective is to develop Eulerian criteria for Lagrangian hyperbolicity, i.e., the existence of stable and unstable manifolds. Such criteria are based only on the velocity field and do not require the integration of particles. A necessary criterion [ Haller and Poje, 1998; Poje and Haller, 1999] indicates that fast rotation of Eulerian hyperbolic structures can induce non-hyperbolic flows. This phenomenon is the major source of difficulties in establishing Eulerian rules for Lagrangian hyperbolicity.
Simplified analytical models such as the one depicted in
Fig. 1 reveals typical behaviors and provides approximate
criteria for hyperbolicity that are applicable in more general
settings [ Lekien, 2003].
One such example is the break-up of the Antarctic Polar
vortex. During September 2002, the polar vortex became
unstable and split into two vortices, each holding a smaller
ozone hole and stretching outside the polar area. Figure 2
shows the hyperbolic trajectory at the South pole and its
associated rotating invariant manifolds [ Lekien et al., 2005d].
The 4 point vortex model of the splitting event [ Newton, 2001] is quasi-analytical and has a fixed point at the South pole,
where the dynamics is equivalent to the rotating saddle of
Fig. 1. Similar developments lead to the identification of
vortices [ Haller, 2005].
A.1 Past Accomplishments
Figure 3. Left panel. DLE field and LCS
in high-frequency radar data collected on
the coast of Florida [Lekien et al., 2005b].
Right panel. The estimated flux along the
LCS is negligible with respect to typical currents in the area.
Figure 4. Lagrangian separation at the boundary
of Rayleigh-Bénard convection cells (upside down).
Red triangle: Separation point, Blue triangle: Re-attachment point, Green dots : Particles released in
the vicinity of the boundary.
Lagrangian Coherent Structures. The Lagrangian behavior of the flow during a finite interval of time can be extracted from Finite-Time Lyapunov Exponents (FTLE), a measure of the separation rate. The Direct Lyapunov Exponent (DLE) is a robust method to compute these exponents [ Haller, 2000; Haller, 2001]. Ridges in the DLE field approximate Lagrangian Coherent Structures (LCS's) [ Shadden et al., 2005]. We have shown that these ridges render a geometric description of the flow [ Lekien et al., 2005b].
Figure 3 shows an LCS based on coastal high-frequency
radar data. The ridge is the boundary between two regions
of qualitatively different dynamics: a recirculation zone and
fast Northward current.
Reduced analysis on boundaries: Separation. Hyperbolic stretching near a boundary is best studied in terms of the existence of a separation profile emanating from the boundary [ Sears and Telionis, 1975; Van Dommelen and Cowley, 1990]. A necessary criterion for the occurrence of separation at the boundary of time-chaotic two-dimensional flows is given in [ Haller, 2004]. In recent
work [ Lekien and Haller, 2005], we extended these results to boundaries with a free-slip boundary condition (e.g., geophysical flows, nano-fluids). The separation criterion developed is necessary (based on [ [Mané, 1978]) and sufficient (based on [ Fenichel, 1971]).
A.2 Present and Future Work
Figure 5. LCS's for a cylindrical Rayleigh-Bénard convection cell. The two LCS's intersect and form lobes. These lobes carry fluid to and from the core cylindrical cell.
Figure 6 describes transport in a three-layer quasigeostrophic model of the North Atlantic [ Rowley, 1996]. Intergyre transport is governed by the motion of two sequences of lobes moving along the Gulf stream. In particular, the pinching of the invariant manifolds is responsible for the formation of large rings in the ocean [ Haller and Poje, 1998]. I seek global dynamical constraints that provide a geometrical understanding of transport. For example, the uniqueness of trajectory and the continuity of the flow can be used to prove that only anti-cyclonic rings are involved in transport across a meandering jet [ Lekien, 2003]. I am also investigating how lobe dynamics in higher dimensional spaces (based on [ Beigie, 1995]) can be adapted to time-chaotic systems such as the three-dimensional Rayleigh-Bénard cells in Fig. 5.
Figure 6. Lobe dynamics and transport in a layered quasigeostrophic (QG) model of the North Atlantic.
Finite-Size Lyapunov Exponent. The coherent structures are defined using Finite-Time Lyapunov Exponents (FTLE) in [ Shadden et al., 2005]. The works of [ Joseph and Legras, 2002] and [ Koh and Legras, 2002] demonstrate the use of Finite-Size Lyapunov Exponent (FSLE) for locating LCS's in atmospheric data. G. Haller and I are currently investigating the relationship between FSLE and FTLE. In particular we analyze the curvature of these fields, a critical quantity determining the properties of the LCS's [ Shadden et al., 2005; Lekien, 2005].
Physical constraints. Early works in dynamical systems theory were mainly applied to kinematically defined velocity fields. Using dynamically evolving fields brings more physical relevance to these transport studies [ Del-Castillo-Negrete and Morrison, 1993; Ngan and Shepherd, 1997]. A next step is to investigate the effect of these dynamical constraints on the existence and properties of LCS's [ Rogerson et al., 1999]. For example, the separation criterion of [ Haller, 2004] is derived using the assumption that the density of the fluid remains bounded. Physical constraints facilitate the development of theorems and methods applicable to physical systems. One of my objectives is to show that Navier-Stokes equations cannot sustain a sufficiently high rotation rate, hence non-hyperbolic rotating saddles (see Fig. 1) do not exist in real fluids.
I am also interested in investigating LCS's for particles with a non-zero volume [ Cartwright et al., 2002]. The motion of neutrally buoyant particles in a two-dimensional aperiodic flow can be analyzed in a 4 dimensional phase space [ Babiano et al., 2000]. Due to the fast action of on-board controllers, the dynamics of complex (non-spherical) underwater vehicles can also be constrained to low dimensional spaces. The differences between the LCS's for point mass particles and immersed solids is an important theoretical question. Furthermore, it is a critical to investigate the effect on the accuracy of the barriers to contaminants and the alleyways for autonomous vehicles.
B Numerical Methods
B.1 Past Accomplishments
Interpolation. Dynamical systems theory typically applies to velocity fields that are sufficiently smooth. This requires a careful interpolation of the data. In [ Lekien and Marsden, 2005], we provide C 1 interpolation for two-dimensional time-dependent flows. Modal analysis [ Lipphardt et al., 2000; Chu et al., 2003] is used to interpolate, extrapolate and filter incomplete datasets. In comparison to orthogonal modes [ Rowley et al., 2004], normal modes do not depend on time or on the measured data and can be used during long real-time experiments. In [ Lekien et al., 2004], we develop Open-Boundary Modes (OMA) that are suitable for assimilating radar data in coastal regions such as Monterey Bay [ Paduan and Rosenfeld, 1996; Paduan and Cook, 1997].
Figure 7. MANGEN: Graphical User Interface.
Coherent Packages. MANGEN (MANifold GENerator) is a collection of tools that I wrote to study mixing and transport in experimental fluid flows. The interface is flexible and permits the computation of Lagrangian coherent structures (LCS's) and separation profiles for a wide range of problems. The abstraction layer between the dynamical systems methods and the physical systems (e.g., high-frequency radar data, model data, Hamiltonian) permits the development of new tools and their direct application to all the systems configured.
MANGEN has been successfully exported to other institutions, including the Monterey Bay Aquarium Research Institute (for Lagrangian studies during the Adaptive Ocean Sampling Network (AOSN-II) project), the Jet Propulsion Laboratory (as a means to offer end-user products based on the ROMS ocean model), the Rosentiel School of Marine and Atmospheric Science (to investigate the properties of the solutions of the parabolic wave equation) and the French Research Institute for Exploitation of the Sea (to study Lagrangian residual circulation in the Iroise Sea).
Concurrent Programming and Real-Time Applications. Real-time applications and high resolution simulations require the use of concurrent algorithms for the computation of LCS's and separation profiles. Parallel code is currently available in MANGEN and is used on clusters at Caltech, JPL, and Princeton University. The real-time capabilities of MANGEN were during the Adaptive Ocean Sampling Network (AOSN-II) experiment in August 2003. High-frequency radar data and ocean model output were assimilated in real-time to display LCS's in the Monterey Bay Aquarium Research Institute's control room.
B.2 Present and Future Work
Lagrangian Coherent Structures. For the purpose of extracting LCS's (as opposed to observing them on a FTLE plot), a precise definition is developed in [ Shadden et al., 2005; Lekien et al., 2005e]. Ongoing work [ Lekien, 2005] reduces the definition and shows that LCS's can be extracted as the level sets of a computable function. The extraction of the LCS's is based on a finite-element method. Consequently, I am also developing more efficient and accurate FTLE computations on finite-element meshes. In particular, the algorithms of MANGEN have been adapted to work on arbitrary manifolds, instead of only Cartesian subspaces. For example, LCS's are computed on a spherical Earth in [ Lekien et al., 2005d].
I am also investigating the relationship between LCS's and almost invariant sets [ Dellnitz and Preis, 2003]. In [ Dellnitz et al., 2005], we combine and compare these two techniques for the restricted three-body problem. The formulation of lobe dynamics and transport problems in terms of almost invariant sets permits the use of highly efficient algorithms based on graph theory [ Dellnitz et al., 2005]. Conversely, I plan to study how the extraction of LCS's can be cast as a set-oriented problem and solved as such. The objective is to replace the stochastic Frobenius-Perron operator in [ Dellnitz and Preis, 2003] by a deterministic transfer function satisfying the objective of [ Shadden et al., 2005].
Optimal Filtering. The recent use of OMA modes [ Lekien et al., 2004] for biological studies in Bodega Bay [ Kaplan et al., 2005] suggested that OMA modes can be generalized and optimized. In particular, a new class of modes are currently developed [ Kaplan and Lekien, 2005] to avoid the artificial vanishing vorticity near the coastline.
Integration with Ocean and Biological Models. Interfaces for MANGEN are developed to interact directly with experimental data and advanced ocean models. In particular, the parallel code of MANGEN has been installed at the Jet Propulsion Laboratory to function directly with the Regional Ocean Modeling System (ROMS). In addition, interfaces are developed for the Harvard Ocean Model and Prediction System (HOPS) allowing the computation of LCS for real-time forecasts of the ocean as well as the direct integration of modeled uncertainties [ Lermusiaux and Lekien, 2005].
C Applications
C.1 Past Accomplishments
Figure 8. Recirculating zone near Fort-Lauderdale, FL..
Pollution Control. Coherent structures in the ocean mark the boundary between regions of qualitatively different dynamics. The impact of contaminants in a coastal environment can be reduced by the use of a release schedule based on the motion of nearby LCS's.
In [ Lekien et al., 2005b], we show that the average concentration of pollutants and the peaks of maximum concentration can be significantly reduced using the LCS's observed in high-frequency radar data collected on the coast of Florida (Fig. 8).
In [ Lekien et al., 2005a], we show the existence of highly sensitive intervals of time in the coastal dynamics of Monterey Bay. Avoiding dangerous operations such as the dismantlement of obsolete oil tanks during these periods lessens the disastrous outcome in the event of an accidental oil spill.
Fuel Efficient Routes. During the Autonomous Ocean Sampling Network ( AOSN-II) experiment in August 2003, autonomous underwater vehicles (SIO and Slocum gliders) were used to sample geophysical processes in Monterey Bay, CA [ Fiorelli et al., 2003]. Coordinated control maintained either specific formations (small-scale sampling) or coordinated arrays (large-scale coverage). Real-time LCS computations permit the optimization of the paths followed by the weakly actuated vehicles ( 25 cm/s) and take advantage of the currents (as high as 100 cm/s). LCS's indicate regions of high stretching and are avoided by vehicles formations [ Bhatta et al., 2005]; Isolated vehicles follow the structures during transient periods to reach their assigned position quickly [ Inanc et al., 2004].
Similar methodologies can be applied to multi-body problems in dynamical astronomy, molecular dynamics, and chemical reaction kinetics. In [ Dellnitz et al., 2005], we focus on the transport rates between two resonance regions for the three body system consisting of the Sun, Jupiter and a third body (such as an asteroid or an autonomous vehicle). The dynamics of the lobes formed by the intersection of the LCS's provide fuel-efficient routes in the solar system.
Control of fluid flows. The recent development of exact criteria for separation in fluid flows [ Haller, 2004; Lekien and Haller, 2005] opens the door to a variety of approaches to actuate the flow and induce separation at specified positions (and with a specified strength or separation angle). Such a method is described in [ Insperger et al., 2004; Lekien et al., 2005c] and permits a rapid control of the lift on a wing or the transport of fuel within the fluid.
|
|
|
Figure 9. Two wall jets used as actuators and inducing fixed separation in the presence of an unknown background flow.
|
Figure 10. A delayed dead-beat controller [Lekien et al., 2005c] produces Lagrangian patterns (green lines) revealing separation at the desired position (blue circle).
|
C.2 Present and Future Work
Fragmentation Dynamics. As an example of a time-dependent 4-dimensional system, I am investigating the Van Der Waals complex He-I 2 using the potential of [ Gaspard and Rice, 1989]. An alternative to statistical studies [ MacKay et al., 1984a; MacKay et al., 1984b] of dissociation of the compound is the geometrical approach based on the motion of the lobes between the LCS's. The first 6 steps of the dynamics of the lobes involved in the decay of a restricted two-dimensional (T-shaped) model of the molecule can be found on Fig. 11. Figure 12 shows a section of the LCS's of the full problem where the lobes are bounded by pieces of three-dimensional LCS's [ Gillilan and Ezra, 1991; Beigie, 1995].
|
|
|
Figure 11. Geometrical description of the fragmentation of He-I2 for the reduced two-dimensional model. An initial ensemble of molecules (shaded area) dissociates progressively as a result of lobes invading the interior of the the quasi-invariant region.
|
Figure 12. LCS and lobes in the full 4-dimensional phase space of HeI2
|
Figure 13. Atmospheric LCS's based on NCEP/NCAR reanalysis.
Mixing and Transport in the Atmosphere. In collaboration with G. Haller and R.T. Pierrehumbert, I am computing LCS's for the atmosphere using the NCEP/NCAR Reanalysis Project [ Kalnay et al., 1996], an analysis (and forecast) system performing data assimilation using past data from 1948 to the present. The objective is to identify Lagrangian processes and quasi-invariant regions in the reanalysis data (Fig. 13).
Ocean Sampling. The Adaptive Sampling And Prediction ( ASAP) project aim to develop a sustainable, portable, adaptive ocean observing and prediction system for use in coastal environments [ Leonard et al., 2005]. This project employs, among other observation platforms, autonomous underwater vehicles that carry sensors to measure physical and biological signals in the ocean. The measurements from all sensing platforms are assimilated in real-time into advanced ocean models. The objective is to coordinate the mobile assets in order to collect data of highest possible value. To this purpose, optimal sampling arrays are designed based on minimizing the residual error of an objective analysis scheme or advanced ocean assimilation models [ Leonard et al., 2005; Bhatta et al., 2005]. Symmetries in the optimal solutions are broken by properly aligning the arrays with the LCS's.
Figure 14. Left panel: FTLE and LCS based on HOPS ensemble average.
Right panel: Corresponding error based on ensemble statistics.
Uncertainties. The impact of uncertainties on the LCS's is an important question for practical applications. Realistic ocean modeling with dense data assimilation provides accurate nowcasts and prediction of ocean flow fields as well as associated uncertainties [ [Jazwinski, 1970; Lermusiaux, 2005]. P.F.J. Lermusiaux and I are investigating the transfer of these uncertainties from the ocean state to the LCS's [ Lermusiaux and Lekien, 2005]. Figure 14 shows that the "upwelling barrier" in Monterey Bay is a robust LCS with minimum uncertainty [ Lermusiaux et al., 2005].
Long Term Objectives
A central objective of my research program is to enhance the understanding and use of Lagrangian Coherent Structures (LCS's) in time-chaotic dynamical systems. Classical results are concerned with autonomous and time-periodic velocity fields. For turbulent flows (or chaotic systems), the velocity field (or state of the system) does not reach an equilibrium and does not follow a periodic pattern in phase space. Following D. Ruelle and F. Takens, I assume that the state of the turbulent flow is wrapping around a strange attractor [ Ruelle and Takens, 1971]. Experimental measurements and analytical examples have shown that dissipative systems such as viscous fluids will not have quasi-periodic motions. The understanding and control of turbulence with dynamical systems methods requires the generalization of these methods to arbitrary time-dependence [ Ruelle and Takens, 1971].
The study of the relationships between the dynamics in phase space and the LCS's in the physical space will provide a deeper understanding of the mechanisms and implications of turbulence, as well as new methods for the control of transport and mixing in turbulent flows. Lagrangian patterns in the physical space are the meaningful observables. The Eulerian velocity field (i.e., the state in phase space) is not measured directly but derived from particle trajectories (e.g., derivative of drifter paths, Digital Particle Image Velocimetry (DPIV) or Doppler effect). In Engineering and Science, the knowledge or the prediction of the velocity is only used as a means to forecast or control the path of actual particles. For example, in a search and rescue mission, the path followed by the drifting body must be forecast, not the currents. My research program approaches turbulence by considering instabilities in the Navier-Stokes equations and observing them from the point of view of a drifting body, not the velocity field. The motivation for this approach is twofold: first, turbulence is studied from a physically relevant coordinate system; chaotic disturbances that do not influence the motion of passive drifters are ignored. Second, it opens the door to a possibly simpler description of turbulence. The apparition of a strange attractor in phase space [ Ruelle and Takens, 1971] is interpreted by using (smooth) moving and bifurcating LCS's in physical space.
My approach is based exclusively on time-dependent dynamical models and is well suited for the direct application to geophysical flows and to experimentally measured velocity fields. In the future, the theory and numerical methods will adapt to particles with a finite volume (pollutants or vehicles) and incorporate statistical models. This will increase the utility of these methods in Science and Engineering and provide a clearer mathematical and physical understanding of hyperbolicty, mixing and separation.
Literature Cited
-
Babiano, A., Cartwright, J. H. E., Piro, O., and Provenzale, A. (2000). Dynamics of a small neutrally buoyant sphere in a fluid and targeting in Hamiltonian systems. Phys. Rev. Let., 84(25):5764-5767.
-
Babiano, A., Provenzale, A., and Vulpiani, A. (1994). Chaotic advection, tracer dynamics, and turbulent dispersion Proceedings of the NATO Advanced Research Workshop and EGS Topical Workshop on Chaotic Advection, Conference Center Sereno di Gavo, Italy, 24-28 May 1993. Physica D, 76(1-3):1-329.
-
Beigie, D. (1995). Multiple separatrix crossing in multi-degree-of-freedom Hamiltonian flows. Journal of Nonlinear Science,5(1):57-103.
-
Bhatta, P., Fiorelli, E., Lekien, F., Leonard, N. E., Paley, D. A., Zhang, F., Bachmayer, R., Davis, R. E., Fratantoni, D., and Sepulchre, R. (2005). Coordination of an underwater fleet for adaptive sampling. In International Workshop on Underwater Robotics, pages 61-69.
-
Cartwright, J. H. E., Magnasco, M. O., Piro, O., and Tuval, I. (2002). Bailout embeddings and neutrally buoyant particles in three--dimensional flows. Phys. Rev. Let., 89(26):264501.1-4.
-
Chu, P. C., Ivanov, L. M., Korzhova, T. P., Margolina, T. M., and Melnichenko, V. M. (2003). Analysis of sparse and noisy ocean current data using flow decomposition. J. of Atmospheric and Oceanic Tech., 20:478-491.
-
Coppel, W. A. (1978). Dichotomies in Stability Theory. Springer-Verlag, Berlin, Heidelberg, New-York.
-
Coulliette, C. and Wiggins, S. (2001). Intergyre transport in a wind-driven, quasigeostrophic double gyre: An application of lobe dynamics. Nonlinear Proc. Geophys, 8:69-94.
-
Del-Castillo-Negrete, D. and Morrison, P. J. (1993). Chaotic transport of Rossby waves in shear flow. Phys. Fluids A, 5(4):948-965.
-
Dellnitz, M., Junge, O., Koon, W.-S., Lekien, F., Lo, M., Marsden, J., Padberg, K., Preisand, R., Ross, S., and Thiere, B. (2005). Transport in dynamical astronomy and multibody problems. Int. J. of Bifurcation and Chaos, 15(3):699-727.
-
Dellnitz, M. and Preis, R. (2003). Congestion and almost invariant sets in dynamical systems. Symbolic and Numerical Scientific Computation, 2630:183-209.
-
Fenichel, N. (1971). Persistence and smoothness of invariant manifolds for flows. Indiana University Mathematics Journal, 21(3):193-226.
-
Fiorelli, E., Bhatta, P., Leonard, N. E., and Schulman, I. (2003). Adaptive sampling using feedback control of an autonomous underwater glider fleet. In Proceedings 13th. Int. Symp. on Unmanned Untethered Submersible Technology.
-
Gaspard, P. and Rice, S. A. (1989). Hamiltonian mapping models of molecular fragmentation. Journal of Physical Chemistry, 93(19):6947-6957.
-
Gillilan, R. E. and Ezra, G. S. (1991). Transport and turnstiles in multidimensional Hamiltonian mappings for unimolecular fragmentation - application to Vanderwaals predissociation. Journal of Chemical Physics, 94(4):2648-2668.
-
Guckenheimer, J. and Holmes, P. J. (1983). Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields. Springer-Verlag, New York, Heidelberg, Berlin.
-
Haller, G. (2000). Finding finite-time invariant manifolds in two-dimensional velocity fields. Chaos, 10:99-108.
-
Haller, G. (2001). Lagrangian structures and the rate of strain in a partition of two-dimensional turbulence. Phys. Fluids A, 13:3368-3385.
-
Haller, G. (2004). Exact theory of unsteady separation for two-dimensional flows. Journal of Fluid Mechanics, 512:257-311.
-
Haller, G. (2005). An objective definition of a vortex. J. of Fluid Mechanics, 525:1-26.
-
Haller, G. and Poje, A. (1998). Finite-time transport in aperiodic flows. Physica D, 119(3-4):352-380.
-
Hirsch, M. W. and Smale, S. (1974). Differential Equations, Dynamical Systems and Linear Algebra. Academic Press, New York.
-
Ide, K., Small, D., and Wiggins, S. (2002). Distinguished hyperbolic trajectories in time-dependent fluid flows: analytical and computational approach for velocity fields defined as data sets. Nonlinear Processes in Geophysics, 9(3-4):237-263.
-
Inanc, T., Shadden, S., and Marsden, J. (2004). Optimal trajectory generation in ocean flows. In Proceedings 43rd IEEE Conf. Decision and Control.
-
Insperger, T., Lekien, F., Salman, H., Haller, G., and Stepan, G. (2004). Control of separation in periodic flows including delay effects. In Proc. of IFAC Workshop on Adaptation and Learning in Control and Signa Processing, and IFAC Workshop on Periodic Control Systems, Yokomaha, Japan, volume ALCOSP/PSYCO 04, pages 451-455.
-
Jazwinski, A. H. (1970). Stochastic Processes and Filtering Theory. Academic Press.
-
Joseph, B. and Legras, B. (2002). Relation between kinematic boundaries, stirring, and barriers for the Antarctic polar vortex. J. Atmos. Sci., 59:1198-1212.
-
Kalnay, E., Kanamitsu, M., Kistler, R., Collins, W., Deaven, D., Gandin, L., Iredell, M., Saha, S., White, G., Woollen, J., Zhu, Y., Chelliah, M., Ebisuzaki, W., Higgins, W., Janowiak, J., Mo, K. C., Ropelewski, C., Wang, J., Leetmaa, A., Reynolds, R., Jenne, R., and Joseph, D. (1996). The NCEP/NCAR 40-year reanalysis project. Bulletin of the American Meteorological Society, 77(3):437-471.
-
Kaplan, D. M., Largier, J. L., and Botsford, L. W. (in press, 2005). HF radar observations of surface circulation off Bodega Bay (Northern California, USA). J. Geophys. Res. - Oceans.
-
Kaplan, D. M. and Lekien, F. (in preparation, 2005). Optimal boundary modal decomposition. J. of Geophys. Res. - Oceans.
-
Koh, T.-Y. and Legras, B. (2002). Hyperbolic lines and the stratospheric polar vortex. Chaos, 12:382-394.
-
Lekien, F. (2003). Time-Dependent Dynamical Systems and Geophysical Flows. PhD thesis, California Institute of Technology.
-
Lekien, F. (in preparation, 2005). Computation and extraction of Lagrangian coherent structures in time-chaotic flows. Physica D.
-
Lekien, F., Coulliette, C., Bank, R., and Marsden, J. E. (2004). Open-boundary modal analysis: Interpolation, extrapolation, and filtering. J. Geophys. Res. Oceans, 109(Art. No. C12004):1-13.
-
Lekien, F., Coulliette, C., Haller, G., Paduan, J., and Marsden, J. E. (submitted, 2005a). Optimal polution release in Monterey Bay based on nonlinear analysis of coastal radar data. Environmental Science and Technology.
-
Lekien, F., Coulliette, C., Mariano, A. J., Ryan, E. H., Shay, L. K., Haller, G., and Marsden, J. E. (2005b). Pollution release tied to invariant manifolds: A case study for the coast of Florida. Physica D, 210(12):1-20.
-
Lekien, F. and Haller, G. (submitted, 2005). Unsteady separation criterion for free-slip walls. Phys. of Fluids.
-
Lekien, F., Insperger, T., Stepan, G., Salman, H., and Haller, G. (submitted, 2005c). Control of separation in two-dimensional flows: Stability and robustness for delayed oscillatory flows. Automatica.
-
Lekien, F. and Marsden, J. (2005). Tricubic interpolation in three dimensions. Int. J. for Numerical Methods in Engineering, 63(3):455-471.
-
Lekien, F., Ross, S. D., and Newton, P. K. (in preparation, 2005d).
-
Lekien, F., Shadden, S. C., and Marsden, J. E. (submitted, 2005e). Lagrangian coherent structures in n-dimensional systems. Physica D.
-
Leonard, N. E., Paley, D. A., Lekien, F., Sepulchre, R., Fratantoni, M., and Davis, R. E. (submitted, 2005). Collective motion, sensor networks and ocean sampling. In Proceedings of the IEEE, Special Issue on Networked Control Systems.
-
Lermusiaux, P. F. J. (submitted, 2005). Uncertainty estimation and prediction for the interdisciplinary ocean. J. Comp. Phys.
-
Lermusiaux, P. F. J., Chiu, C.-S., Gawarkiewicz, G. G., Abbot, P., Robinson, A. R., Miller, R. N., Haley, P. J., Leslie, W. G., Majumdar, S. J., Pang, A., and Lekien, F. (submitted, 2005). Quantifying uncertainties in ocean predictions. Oceanography Magazine.
-
Lermusiaux, P. F. J. and Lekien, F. (2005). Dynamics and Lagrangian coherent structures in the ocean and their uncertainty. In Marsden, J. and Scheurle, J., editors, Dynamical Systems Methods in Fluid Dynamics, volume 531, pages 19-20. Mathematisches Forschungsinstitut Oberwolfach.
-
Lipphardt, Jr., B. L., Kirwan, Jr., A. D., Grosch, C. E., Lewis, J. K., and Paduan, J. D. (2000). Blending HF radar and model velocities in Monterey Bay through normal mode analysis. J. of Geophys. Res. Oceans, 105(C2), 3:425-450.
-
Mané, R. (1978). Persistent manifolds are normally hyperbolic. Transactions of the American Mathematical Society, 246(DEC):261-283.
-
MacKay, R. S., Meiss, J. D., and Percival, I. C. (1984a). Stochasticity and transport in Hamiltonian-systems. Phys. Rev. Letters, 52(9):697-700.
-
MacKay, R. S., Meiss, J. D., and Percival, I. C. (1984b). Transport in Hamiltonian-systems. Physica D, 13(1-2):55-81.
-
Melnikov, V. K. (1963). On the stability of center for time-periodic perturbations. Trans. Moscow Math. Soc., 12:1-57.
-
Newton, P. (2001). The N-Vortex Problem: Analytical Techniques, volume 145 of Applied Math. Sci. Series. Springer-Verlag.
-
Ngan, K. and Shepherd, T. G. (1997). Chaotic mixing and transport in Rossby-wave critical layers. Physica D, 110:105-122.
-
Paduan, J. D. and Cook, M. S. (1997). Mapping surface currents in Monterey Bay with CODAR-type HR radar. Oceanography, 10:49-52.
-
Paduan, J. D. and Rosenfeld, L. K. (1996). Remotely sensed surface currents in Monterey Bay from shore-based HF radar (coastal ocean dynamics application radar). J. of Geophys. Res. Oceans, 101(C9):20669-20686.
-
Pierrehumbert, R. T. (1991). Chaotic mixing of tracer and vorticity by modulated traveling rossby waves. Geophys. Astrophys. Fluid Dyn., 58:285-319.
-
Poje, A. C. and Haller, G. (1999). Geometry of cross-stream mixing in a double-gyre ocean model. J. Phys. Oceanogr, 29:1649-1665.
-
Rogerson, A., Miller, P. D., Pratt, L. J., and Jones, C. K. R. T. (1999). Lagrangian motion and fluid exchange in a barotropic meandering jet. J. Phys. Oceanogr., 29(10):2635-2655.
-
Rom-Kedar, V. (1990). Transport rates of a class of two-dimensional maps and flows. Physica D, 43(2-3):229-268.
-
Rom-Kedar, V. and Wiggins, S. (1991). Transport in two-dimensional maps: Concepts, examples, and a comparison of the theory of Rom-Kedar and Wiggins with the Markov Model of MacKay, Meiss, Ott, and Percival. Physica D, 51(13):248-266.
-
Rowley, C. A. (1996). A modelling Study of the North Atlantic Current. PhD thesis, University of Rhode Island, Kingston, RI.
-
Rowley, C. W., Colonius, T., and Murray, R. M. (2004). Model reduction for compressible flows using POD and Galerkin projection. Physica D, 189(12):115-129.
-
Ruelle, D. and Takens, F. (1971). On the nature of turbulence. Commun. Math. Phys., 20:167-192.
-
Samelson, R. M. (1992). Fluid exchange across a meandering jet. J. Phys. Oceanogr., 22:431-440.
-
Sandstede, B., Balasuriya, S., Jones, C. K. R. T., and Miller, P. (2000). Melnikov theory for finite-time vector fields. Nonlinearity, 13(4):1357-1377.
-
Sears, W. R. and Telionis, D. P. (1975). Boundary-layer separation in unsteady flows. SIAM Journal of Applied Mathematics, 28:215.
-
Shadden, S. C., Lekien, F., and Marsden, J. E. (2005). Definition and properties of Lagrangian coherent structures from finite-time Lyapunov exponents in two-dimensional aperiodic flows. Physica D, 212(3-4):271-304.
-
Smale, S. (1980). The Mathematics of Time: Essays on Dynamical Systems, Economic Processes, and Related Topics. Springer-Verlag, New-York, Heidelberg, Berlin.
-
Solomon, T. H. and Gollub, J. P. (1988a). Chaotic particle-transport in time-dependent rayleigh-b´enard convection. Physical Review A, 38(12):6280-6286.
-
Solomon, T. H. and Gollub, J. P. (1988b). Passive transport in steady Rayleigh-B´enard convection. Phys. of Fluids, 31(6):1372-1379.
-
Van Dommelen, L. L. and Cowley, S. J. (1990). On the Lagrangian description of unsteady boundary layer separation. Part 1. General theory. J. of Fluid Mech., 210:593-626.
-
Voth, G. A., Haller, G., and Gollub, J. P. (2002). Experimental measurements of stretching fields in fluid mixing. Phys. Rev. Letters, 88(25):254501.1 - 254501.4.
-
Wiggins, S. (1992). Chaotic Transport in Dynamical Systems. Springer-Verlag, New York.
-
Yuan, G. C., Pratt, L. J., and Jones, C. K. R. T. (2002). Barrier destruction and Lagrangian predictability at depth in a meandering jet. Dynamics of Atmospheres and Oceans, 35(1):41-61.
-
Yuan, G. C., Pratt, L. J., and Jones, C. K. R. T. (2004). Cross-jet Lagrangian transport and mixing in a 2 1/2-layer model. Journal of Physical Oceanography, 34(9):1991-2005.
|