Showing 1–19 of 19 results for author: Stechmann, S N

Beyond Linear Decomposition: a Nonlinear Eigenspace Decomposition for a Moist Atmosphere with Clouds

Authors: Antoine Remond-Tiedrez, Leslie M. Smith, Samuel N. Stechmann

Abstract: A linear decomposition of states underpins many classical systems. This is the case of the Helmholtz decomposition, used to split vector fields into divergence-free and potential components, and of the dry Boussinesq system in atmospheric dynamics, where identifying the slow and fast components of the flow can be viewed as a decomposition. The dry Boussinesq system incorporates two leading ingredi… ▽ More A linear decomposition of states underpins many classical systems. This is the case of the Helmholtz decomposition, used to split vector fields into divergence-free and potential components, and of the dry Boussinesq system in atmospheric dynamics, where identifying the slow and fast components of the flow can be viewed as a decomposition. The dry Boussinesq system incorporates two leading ingredients of mid-latitude atmospheric motion: rotation and stratification. In both cases the leading order dynamics are linear so we can rely on an eigendecomposition to decompose states. Here we study the extension of dry Boussinesq to incorporate another important ingredient in the atmosphere: moisture and clouds. The key challenge with this system is that nonlinearities are present at leading order due to phase boundaries at cloud edge. Therefore standard tools of linear algebra, relying on eigenvalues and eigenvectors, are not applicable. The question we address in this paper is this: in spite of the nonlinearities, can we find a decomposition for this moist Boussinesq system? We identify such a decomposition adapted to the nonlinear balances arising from water phase boundaries. This decomposition combines perspectives from partial differential equations (PDEs), the geometry, and the conserved energy. Moreover it sheds light on two aspects of previous work. First, this decomposition shows that the nonlinear elliptic PDE used for potential vorticity and moisture inversion can be used outside the limiting system where it was first derived. Second, we are able to rigorously justify, and interpret geometrically, an existing numerical method for this elliptic PDE. This decomposition may be important in applications because, like its linear counterparts, it may be used to analyze observational data. Moreover, by contrast with previous decompositions, it may be used even in the presence of clouds. △ Less

Submitted 17 May, 2024; originally announced May 2024.

MSC Class: Primary: 86A10; 76U60; Secondary: 37N10; 53Z05; 65N99

LEMDA: A Lagrangian-Eulerian Multiscale Data Assimilation Framework

Authors: Quanling Deng, Nan Chen, Samuel N. Stechmann, Jiuhua Hu

Abstract: Lagrangian trajectories are widely used as observations for recovering the underlying flow field via Lagrangian data assimilation (DA). However, the strong nonlinearity in the observational process and the high dimensionality of the problems often cause challenges in applying standard Lagrangian DA. In this paper, a Lagrangian-Eulerian multiscale DA (LEMDA) framework is developed. It starts with e… ▽ More Lagrangian trajectories are widely used as observations for recovering the underlying flow field via Lagrangian data assimilation (DA). However, the strong nonlinearity in the observational process and the high dimensionality of the problems often cause challenges in applying standard Lagrangian DA. In this paper, a Lagrangian-Eulerian multiscale DA (LEMDA) framework is developed. It starts with exploiting the Boltzmann kinetic description of the particle dynamics to derive a set of continuum equations, which characterize the statistical quantities of particle motions at fixed grids and serve as Eulerian observations. Despite the nonlinearity in the continuum equations and the processes of Lagrangian observations, the time evolutions of the posterior distribution from LEMDA can be written down using closed analytic formulae. This offers an exact and efficient way of carrying out DA, which avoids using ensemble approximations and the associated tunings. The analytically solvable properties also facilitate the derivation of an effective reduced-order Lagrangian DA scheme that further enhances computational efficiency. The Lagrangian DA within the framework has advantages when a moderate number of particles is used, while the Eulerian DA can effectively save computational costs when the number of particle observations becomes large. The Eulerian DA is also valuable when particles collide, such as using sea ice floe trajectories as observations. LEMDA naturally applies to multiscale turbulent flow fields, where the Eulerian DA recovers the large-scale structures, and the Lagrangian DA efficiently resolves the small-scale features in each grid cell via parallel computing. Numerical experiments demonstrate the skilful results of LEMDA and its two components. △ Less

Submitted 4 February, 2024; v1 submitted 31 January, 2024; originally announced January 2024.

Estimating the time-evolving refractivity of a turbulent medium using optical beam measurements: a data assimilation approach

Authors: Anjali Nair, Qin Li, Samuel N. Stechmann

Abstract: In applications such as free-space optical communication, a signal is often recovered after propagation through a turbulent medium. In this setting, it is common to assume that limited information is known about the turbulent medium, such as a space- and time-averaged statistic (e.g., root-mean-square), but without information about the state of the spatial variations. It could be helpful to gain… ▽ More In applications such as free-space optical communication, a signal is often recovered after propagation through a turbulent medium. In this setting, it is common to assume that limited information is known about the turbulent medium, such as a space- and time-averaged statistic (e.g., root-mean-square), but without information about the state of the spatial variations. It could be helpful to gain more information if the state of the turbulent medium can be characterized with the spatial variations and evolution in time described. Here, we propose to investigate the use of data assimilation techniques for this purpose. A computational setting is used with the paraxial wave equation, and the extended Kalman filter is used to conduct data assimilation using intensity measurements. To reduce computational cost, the evolution of the turbulent medium is modeled as a stochastic process. Following some past studies, the process has only a small number of Fourier wavelengths for spatial variations. The results show that the spatial and temporal variations of the medium are recovered accurately in many cases. In some time windows in some cases, the error is larger for the recovery. Finally we discuss the potential use of the spatial variation information for aiding the recovery of the transmitted signal or beam source. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Non-conservation and conservation for different formulations of moist potential vorticity

Authors: Parvathi Kooloth, Leslie M. Smith, Samuel N. Stechmann

Abstract: Potential vorticity (PV) is one of the most important quantities in atmospheric science. The PV of each fluid parcel is known to be conserved in the case of a dry atmosphere. However, a parcel's PV is not conserved if clouds or phase changes of water occur. Recently, PV conservation laws were derived for a cloudy atmosphere, where each parcel's PV is not conserved but parcel-integrated PV is conse… ▽ More Potential vorticity (PV) is one of the most important quantities in atmospheric science. The PV of each fluid parcel is known to be conserved in the case of a dry atmosphere. However, a parcel's PV is not conserved if clouds or phase changes of water occur. Recently, PV conservation laws were derived for a cloudy atmosphere, where each parcel's PV is not conserved but parcel-integrated PV is conserved, for integrals over certain volumes that move with the flow. Hence a variety of different statements are now possible for moist PV conservation and non-conservation, and in comparison to the case of a dry atmosphere, the situation for moist PV is more complex. Here, in light of this complexity, several different definitions of moist PV are compared for a cloudy atmosphere. Numerical simulations are shown for a rising thermal, both before and after the formation of a cloud. These simulations include the first computational illustration of the parcel-integrated, moist PV conservation laws. The comparisons, both theoretical and numerical, serve to clarify and highlight the different statements of conservation and non-conservation that arise for different definitions of moist PV. △ Less

Submitted 15 November, 2023; originally announced November 2023.

Element learning: a systematic approach of accelerating finite element-type methods via machine learning, with applications to radiative transfer

Authors: Shukai Du, Samuel N. Stechmann

Abstract: In this paper, we propose a systematic approach for accelerating finite element-type methods by machine learning for the numerical solution of partial differential equations (PDEs). The main idea is to use a neural network to learn the solution map of the PDEs and to do so in an element-wise fashion. This map takes input of the element geometry and the PDEs' parameters on that element, and gives o… ▽ More In this paper, we propose a systematic approach for accelerating finite element-type methods by machine learning for the numerical solution of partial differential equations (PDEs). The main idea is to use a neural network to learn the solution map of the PDEs and to do so in an element-wise fashion. This map takes input of the element geometry and the PDEs' parameters on that element, and gives output of two operators -- (1) the in2out operator for inter-element communication, and (2) the in2sol operator (Green's function) for element-wise solution recovery. A significant advantage of this approach is that, once trained, this network can be used for the numerical solution of the PDE for any domain geometry and any parameter distribution without retraining. Also, the training is significantly simpler since it is done on the element level instead on the entire domain. We call this approach element learning. This method is closely related to hybridizbale discontinuous Galerkin (HDG) methods in the sense that the local solvers of HDG are replaced by machine learning approaches. Numerical tests are presented for an example PDE, the radiative transfer equation, in a variety of scenarios with idealized or realistic cloud fields, with smooth or sharp gradient in the cloud boundary transition. Under a fixed accuracy level of $10^{-3}$ in the relative $L^2$ error, and polynomial degree $p=6$ in each element, we observe an approximately 5 to 10 times speed-up by element learning compared to a classical finite element-type method. △ Less

Submitted 4 August, 2023; originally announced August 2023.

Particle-Continuum Multiscale Modeling of Sea Ice Floes

Authors: Quanling Deng, Samuel N. Stechmann, Nan Chen

Abstract: Sea ice profoundly influences the polar environment and the global climate. Traditionally, Sea ice has been modeled as a continuum under Eulerian coordinates to describe its large-scale features, using, for instance, viscous-plastic rheology. Recently, Lagrangian particle models, also known as the discrete element method (DEM) models, have been utilized for characterizing the motion of individual… ▽ More Sea ice profoundly influences the polar environment and the global climate. Traditionally, Sea ice has been modeled as a continuum under Eulerian coordinates to describe its large-scale features, using, for instance, viscous-plastic rheology. Recently, Lagrangian particle models, also known as the discrete element method (DEM) models, have been utilized for characterizing the motion of individual sea ice fragments (called floes) at scales of 10 km and smaller, especially in marginal ice zones. This paper develops a multiscale model that couples the particle and the continuum systems to facilitate an effective representation of the dynamical and statistical features of sea ice across different scales. The multiscale model exploits a Boltzmann-type system that links the particle movement with the continuum equations. For the small-scale dynamics, it describes the motion of each sea ice floe. Then, as the large-scale continuum component, it treats the statistical moments of mass density and linear and angular velocities. The evolution of these statistics affects the motion of individual floes, which in turn provides bulk feedback that adjusts the large-scale dynamics. Notably, the particle model characterizing the sea ice floes is localized and fully parallelized, in a framework that is sometimes called superparameterization, which significantly improves computation efficiency. Numerical examples demonstrate the effective performance of the multiscale model. Additionally, the study demonstrates that the multiscale model has a linear-order approximation to the truth model. △ Less

Submitted 15 August, 2023; v1 submitted 14 March, 2023; originally announced March 2023.

A nonlinear elliptic PDE from atmospheric science: well-posedness and regularity at cloud edge

Authors: Antoine Remond-Tiedrez, Leslie M. Smith, Samuel N. Stechmann

Abstract: The precipitating quasi-geostrophic equations go beyond the (dry) quasi-geostrophic equations by incorporating the effects of moisture. This means that both precipitation and phase changes between a water-vapour phase (outside a cloud) and a water-vapour-plus-liquid phase (inside a cloud) are taken into account. In the dry case, provided that a Laplace equation is inverted, the quasi-geostrophic e… ▽ More The precipitating quasi-geostrophic equations go beyond the (dry) quasi-geostrophic equations by incorporating the effects of moisture. This means that both precipitation and phase changes between a water-vapour phase (outside a cloud) and a water-vapour-plus-liquid phase (inside a cloud) are taken into account. In the dry case, provided that a Laplace equation is inverted, the quasi-geostrophic equations may be formulated as a nonlocal transport equation for a single scalar variable (the potential vorticity). In the case of the precipitating quasi-geostrophic equations, inverting the Laplacian is replaced by a more challenging adversary known as potential-vorticity-and-moisture inversion. The PDE to invert is nonlinear and piecewise elliptic with jumps in its coefficients across the cloud edge. However, its global ellipticity is a priori unclear due to the dependence of the phase boundary on the unknown itself. This is a free boundary problem where the location of the cloud edge is one of the unknowns. Potential vorticity-and-moisture inversion differs from oft-studied free-boundary problems in the literature, and here we present the first rigorous analysis of this PDE. Introducing a variational formulation of the inversion, we use tools from the calculus of variations to deduce existence and uniqueness of solutions and de Giorgi techniques to obtain sharp regularity results. In particular, we show that the gradient of the unknown pressure or streamfunction is Hölder continuous at cloud edge. △ Less

Submitted 9 April, 2024; v1 submitted 18 January, 2023; originally announced January 2023.

MSC Class: 49J10; 35B65; 86A10 (Primary) 35J20; 35Q86; 76U60 (Secondary)

Scintillation Minimization versus Intensity Maximization in Optimal Beams

Authors: Qin Li, Anjali Nair, Samuel N Stechmann

Abstract: In free-space optical communications and other applications, it is desirable to design optical beams that have reduced or even minimal scintillation. However, the optimization problem for minimizing scintillation is challenging, and few optimal solutions have been found. Here we investigate the general optimization problem of minimizing scintillation and formulate it as a convex optimization probl… ▽ More In free-space optical communications and other applications, it is desirable to design optical beams that have reduced or even minimal scintillation. However, the optimization problem for minimizing scintillation is challenging, and few optimal solutions have been found. Here we investigate the general optimization problem of minimizing scintillation and formulate it as a convex optimization problem. An analytical solution is found and demonstrates that a beam that minimizes scintillation is incoherent light (i.e., spatially uncorrelated). Furthermore, numerical solutions show that beams minimizing scintillation give very low intensity at the receiver. To counteract this effect, we study a new convex cost function that balances both scintillation and intensity. We show through numerical experiments that the minimizers of this cost function reduce scintillation while preserving a significantly higher level of intensity at the receiver. △ Less

Submitted 28 April, 2023; v1 submitted 16 November, 2022; originally announced November 2022.

Systematic assessment of the effects of space averaging and time averaging on weather forecast skill

Authors: Ying Li, Samuel N. Stechmann

Abstract: Intuitively, one would expect a more skillful forecast if predicting weather averaged over one week instead of the weather averaged over one day, and similarly for different spatial averaging areas. However, there are few systematic studies of averaging and forecast skill with modern forecasts, and it is therefore not clear how much improvement in forecast performance is produced via averaging. He… ▽ More Intuitively, one would expect a more skillful forecast if predicting weather averaged over one week instead of the weather averaged over one day, and similarly for different spatial averaging areas. However, there are few systematic studies of averaging and forecast skill with modern forecasts, and it is therefore not clear how much improvement in forecast performance is produced via averaging. Here we present a direct investigation of averaging effects, based on data from operational numerical weather forecasts. Data is analyzed for precipitation and surface temperature, for lead times of roughly 1 to 7 days, and for time- and space-averaging diameters of 1 to 7 days and 100 to 4500 km, respectively. For different geographic locations, the effects of time- or space-averaging can be different, and while no clear geographical pattern is seen for precipitation, a clear spatial pattern is seen for temperature. For temperature, in general, time averaging is most effective near coastlines, also effective over land, and least effective over oceans. Based on all locations globally, time averaging was less effective than one might expect. To help understand why time averaging may sometimes be minimally effective, a stochastic model is analyzed as a synthetic weather time series, and analytical formulas are presented for the decorrelation time. In effect, while time averaging creates a time series that is visually smoother, it does not necessarily cause a substantial increase in the predictability of the time series. △ Less

Submitted 16 July, 2022; originally announced July 2022.

Conservation laws for potential vorticity in a salty ocean or cloudy atmosphere

Authors: Parvathi Kooloth, Leslie M. Smith, Samuel N. Stechmann

Abstract: One of the most important conservation laws in atmospheric and oceanic science is conservation of potential vorticity. The original derivation is approximately a century old, in the work of Rossby and Ertel, and it is related to the celebrated circulation theorems of Kelvin and Bjerknes. However, the laws apply to idealized fluids, and extensions to more realistic scenarios have been problematic.… ▽ More One of the most important conservation laws in atmospheric and oceanic science is conservation of potential vorticity. The original derivation is approximately a century old, in the work of Rossby and Ertel, and it is related to the celebrated circulation theorems of Kelvin and Bjerknes. However, the laws apply to idealized fluids, and extensions to more realistic scenarios have been problematic. Here, these laws are extended to hold with additional fundamental complexities, including salinity in the ocean, or moisture and clouds in the atmosphere. In the absence of these additional complexities, it is known that potential vorticity is conserved following each fluid parcel; here, for a salty ocean or cloudy atmosphere, the general conserved quantity is potential vorticity integrated over certain pancake-shaped volumes. Furthermore, the conservation laws are also related to a symmetry in the Lagrangian, which brings a connection to the symmetry-conservation relationships seen in other areas of physics. △ Less

Submitted 16 June, 2022; originally announced June 2022.

Comments: 11 pages, 3 figures

Hamilton's Principle with Phase Changes and Conservation Principles for Moist Potential Vorticity

Authors: Parvathi Kooloth, Leslie M. Smith, Samuel N. Stechmann

Abstract: Many definitions of moist potential vorticity (PV) have been proposed to extend the dry theory of Ertel PV. None of the moist PV definitions seem to have all of the desirable properties of the dry Ertel PV. For instance, dry PV is not only a globally conserved quantity, but also a material invariant that is conserved along fluid parcel trajectories. Therefore, an open question remains: is there a… ▽ More Many definitions of moist potential vorticity (PV) have been proposed to extend the dry theory of Ertel PV. None of the moist PV definitions seem to have all of the desirable properties of the dry Ertel PV. For instance, dry PV is not only a globally conserved quantity, but also a material invariant that is conserved along fluid parcel trajectories. Therefore, an open question remains: is there a moist PV that is a material invariant, if clouds and phase changes of water are present? In prior studies, definitions of moist PV have been proposed based on physical and mathematical intuition. Here, a systematic approach is used. In particular, a particle relabeling symmetry is devised for a moist atmosphere and then Noether's theorem is employed to arrive at the associated conservation laws for moist PV. A priori, it is not clear whether this systematic approach will be viable, since it relies on variational derivatives in Hamilton's principle, and phase changes introduce singularities that could potentially prevent derivatives at cloud edge. However, it is shown that the energy and the Lagrangian density are sufficiently smooth to allow variational derivatives, in a moist Boussinesq system with reversible phase transitions between water vapor and liquid cloud water.. From the particle relabeling symmetry, a moist Kelvin circulation theorem is found, along with a moist PV conservation law that applies not for each individual parcel but for parcel-integrated PV, integrated over certain local volumes. △ Less

Submitted 15 June, 2022; originally announced June 2022.

Comments: 24 pages, 3 figures

Computation of optimal beams in weak turbulence

Authors: Qin Li, Anjali Nair, Samuel N Stechmann

Abstract: When an optical beam propagates through a turbulent medium such as the atmosphere or ocean, the beam will become distorted. It is then natural to seek the best or optimal beam that is distorted least, under some metric such as intensity or scintillation. We seek to maximize the light intensity at the receiver using the paraxial wave equation with weak-fluctuation as the model. In contrast to class… ▽ More When an optical beam propagates through a turbulent medium such as the atmosphere or ocean, the beam will become distorted. It is then natural to seek the best or optimal beam that is distorted least, under some metric such as intensity or scintillation. We seek to maximize the light intensity at the receiver using the paraxial wave equation with weak-fluctuation as the model. In contrast to classical results that typically confine original laser beams to be from a special class, we allow the beam to be general, which leads to an eigenvalue problem of a large-sized matrix with each entry being a multi-dimensional integral. This is an expensive and sometimes infeasible computational task in many practically reasonable settings. To overcome this, we utilize an asymptotic expansion and transform the derivation to Fourier space, which allows us to incorporate some optional turbulence assumptions, such as homogeneous-statistics assumption, small-length-scale cutoff assumption, and Markov assumption, to reduce the dimension of the numerical integral. The proposed methods provide a computational strategy that is numerically feasible, and results are demonstrated in several numerical examples. △ Less

Submitted 17 August, 2022; v1 submitted 1 April, 2022; originally announced April 2022.

Superfloe Parameterization with Physics Constraints for Uncertainty Quantification of Sea Ice Floes

Authors: Nan Chen, Quanling Deng, Samuel N. Stechmann

Abstract: The discrete element method (DEM) is providing a new modeling approach for describing sea ice dynamics. It exploits particle-based methods to characterize the physical quantities of each sea ice floe along its trajectory under Lagrangian coordinates. One major challenge in applying the DEM models is the heavy computational cost when the number of floes becomes large. In this paper, an efficient La… ▽ More The discrete element method (DEM) is providing a new modeling approach for describing sea ice dynamics. It exploits particle-based methods to characterize the physical quantities of each sea ice floe along its trajectory under Lagrangian coordinates. One major challenge in applying the DEM models is the heavy computational cost when the number of floes becomes large. In this paper, an efficient Lagrangian parameterization algorithm is developed, which aims at reducing the computational cost of simulating the DEM models while preserving the key features of the sea ice. The new parameterization takes advantage of a small number of artificial ice floes, named the superfloes, to effectively approximate a considerable number of the floes, where the parameterization scheme satisfies several important physics constraints. The physics constraints guarantee the superfloe parameterized system will have similar short-term dynamical behavior as the full system. These constraints also allow the superfloe parameterized system to accurately quantify the long-range uncertainty, especially the non-Gaussian statistical features, of the full system. In addition, the superfloe parameterization facilitates a systematic noise inflation strategy that significantly advances an ensemble-based data assimilation algorithm for recovering the unobserved ocean field underneath the sea ice. Such a new noise inflation method avoids ad hoc tunings as in many traditional algorithms and is computationally extremely efficient. Numerical experiments based on an idealized DEM model with multiscale features illustrate the success of the superfloe parameterization in quantifying the uncertainty and assimilating both the sea ice and the associated ocean field. △ Less

Submitted 7 March, 2022; v1 submitted 27 May, 2021; originally announced May 2021.

Convergence of Rain Process Models to Point Processes

Authors: Scott Hottovy, Samuel N. Stechmann

Abstract: A moisture process with dynamics that switch after hitting a threshold gives rise to a rainfall process. This rainfall process is characterized by its random holding times for dry and wet periods. On average, the holding times for the wet periods are much shorter than the dry. Here convergence is shown for the rain fall process to a point process that is a spike train. The underlying moisture proc… ▽ More A moisture process with dynamics that switch after hitting a threshold gives rise to a rainfall process. This rainfall process is characterized by its random holding times for dry and wet periods. On average, the holding times for the wet periods are much shorter than the dry. Here convergence is shown for the rain fall process to a point process that is a spike train. The underlying moisture process for the point process is a threshold model with a teleporting boundary condition. This approximation allows simplification of the model with many exact formulas for statistics. The convergence is shown by a Fokker-Planck derivation, convergence in mean-square with respect to continuous functions, of the moisture process, and convergence in mean-square with respect to generalized functions, of the rain process. △ Less

Submitted 28 November, 2020; originally announced November 2020.

Comments: 11 pages, 2 figures

Multiscale Asymptotics for the Skeleton of the Madden-Julian Oscillation and Tropical-Extratropical Interactions

Authors: Shengqian Chen, Andrew J. Majda, Samuel N. Stechmann

Abstract: A new model is derived and analyzed for tropical-extratropical interactions involving the Madden-Julian oscillation (MJO). The model combines (i) the tropical dynamics of the MJO and equatorial baroclinic waves and (ii) the dynamics of barotropic Rossby waves with significant extratropical structure, and the combined system has a conserved energy. The method of multiscale asymptotics is applied to… ▽ More A new model is derived and analyzed for tropical-extratropical interactions involving the Madden-Julian oscillation (MJO). The model combines (i) the tropical dynamics of the MJO and equatorial baroclinic waves and (ii) the dynamics of barotropic Rossby waves with significant extratropical structure, and the combined system has a conserved energy. The method of multiscale asymptotics is applied to systematically derive a system of ordinary differential equations (ODEs) for three-wave resonant interactions. Two novel features are (i) a degenerate auxiliary problem with overdetermined equations due to a compatibility condition (meridional geostrophic balance) and (ii) cubic self-interaction terms that are not typically found in three-wave resonance ODEs. Several examples illustrate applications to MJO initiation and termination, including cases of (i) the MJO, equatorial baroclinic Rossby waves, and barotropic Rossby waves interacting, and (ii) the MJO, baroclinic Kelvin waves, and barotropic Rossby waves interacting. Resonance with the Kelvin wave is not possible here if only dry variables are considered, but it occurs in the moist model here through interactions with water vapor and convective activity. △ Less

Submitted 28 October, 2015; v1 submitted 3 October, 2015; originally announced October 2015.

Nonlinear traveling waves for the skeleton of the Madden-Julian oscillation

Authors: Shengqian Chen, Samuel N. Stechmann

Abstract: The Madden-Julian Oscillation (MJO) is the dominant component of intraseasonal (30-90 days) variability in the tropical atmosphere. Here, traveling wave solutions are presented for the MJO skeleton model of Majda and Stechmann. The model is a system of nonlinear partial differential equations that describe the evolution of the tropical atmosphere on planetary (10,000-40,000 km) spatial scales. The… ▽ More The Madden-Julian Oscillation (MJO) is the dominant component of intraseasonal (30-90 days) variability in the tropical atmosphere. Here, traveling wave solutions are presented for the MJO skeleton model of Majda and Stechmann. The model is a system of nonlinear partial differential equations that describe the evolution of the tropical atmosphere on planetary (10,000-40,000 km) spatial scales. The nonlinear traveling waves come in four types, corresponding to the four types of linear wave solutions, one of which has the properties of the MJO. In the MJO traveling wave, the convective activity has a pulse-like shape, with a narrow region of enhanced convection and a wide region of suppressed convection. Furthermore, an amplitude-dependent dispersion relation is derived, and it shows that the nonlinear MJO has a lower frequency and slower propagation speed than the linear MJO. By taking the small-amplitude limit, an analytic formula is also derived for the dispersion relation of linear waves. To derive all of these results, a key aspect is the model's conservation of energy, which holds even in the presence of forcing. In the limit of weak forcing, it is shown that the nonlinear traveling waves have a simple sech-squared waveform. △ Less

Submitted 28 October, 2015; v1 submitted 3 October, 2015; originally announced October 2015.

Hasimoto Transformation and Vortex Soliton Motion Driven by Fluid Helicity

Authors: Darryl D. Holm, Samuel N. Stechmann

Abstract: Vorticity filament motions with respect to the Dirac bracket of Rasetti and Regge [1975] are known to be related to the nonlinear Schrödinger equation by the Hasimoto transformation (HT), when the Hamiltonian is the Local Induction Approximation (LIA) of the kinetic energy. We show that when the Hamiltonian is the LIA of Euler-fluid helicity $\int\mathbf{u}\cdot {\rm curl} \mathbf{u}$, the vorti… ▽ More Vorticity filament motions with respect to the Dirac bracket of Rasetti and Regge [1975] are known to be related to the nonlinear Schrödinger equation by the Hasimoto transformation (HT), when the Hamiltonian is the Local Induction Approximation (LIA) of the kinetic energy. We show that when the Hamiltonian is the LIA of Euler-fluid helicity $\int\mathbf{u}\cdot {\rm curl} \mathbf{u}$, the vorticity filament equation of motion under the Rasetti-Regge Dirac bracket is mapped by HT to the integrable complex modified Korteweg-de Vries (cmKdV) equation, the second equation in the nonlinear Schrödinger hierarchy. △ Less

Submitted 20 September, 2004; originally announced September 2004.

Comments: 15 pages

Reduced Singular Solutions of EPDiff Equations on Manifolds with Symmetry

Authors: Darryl D. Holm, Jonathan Munn, Samuel N. Stechmann

Abstract: The EPDiff equation governs geodesic flow on the diffeomorphisms with respect to a chosen metric, which is typically a Sobolev norm on the tangent space of vector fields. EPDiff admits a remarkable ansatz for its singular solutions, called ``diffeons,'' whose momenta are supported on embedded subspaces of the ambient space. Diffeons are true solitons for some choices of the norm. The diffeon sol… ▽ More The EPDiff equation governs geodesic flow on the diffeomorphisms with respect to a chosen metric, which is typically a Sobolev norm on the tangent space of vector fields. EPDiff admits a remarkable ansatz for its singular solutions, called ``diffeons,'' whose momenta are supported on embedded subspaces of the ambient space. Diffeons are true solitons for some choices of the norm. The diffeon solution ansatz is a momentum map. Consequently. the diffeons evolve according to canonical Hamiltonian equations. We examine diffeon solutions on Einstein spaces that are "mostly" symmetric, i.e., whose quotient by a subgroup of the isometry group is 1-dimensional. An example is the two-sphere, whose isometry group $\SO{3}$ contains $S^1$. In this situation, the singular diffeons (called ``Puckons'') are supported on latitudes (``girdles'') of the sphere. For this $S^1$ symmetry of the two-sphere, the canonical Hamiltonian dynamics for Puckons reduces from integral partial differential equations to a dynamical system of ordinary differential equations for their colatitudes. Explicit examples are computed numerically for the motion and interaction of the Puckons on the sphere with respect to the $H^1$ norm. We analyse this case and several other 2-dimensional examples. From consideration of these 2-dimensional spaces, we outline the theory for reduction of diffeons on a general manifold possessing a metric equivalent to the warped product of the line with the bi-invariant metric of a Lie group. △ Less

Submitted 27 September, 2004; v1 submitted 24 February, 2004; originally announced February 2004.

Comments: 46 pages, 5 figures, numerical results added

Rotating Concentric Circular Peakons

Authors: Darryl D. Holm, Vakhtang Putkaradze, Samuel N. Stechmann

Abstract: We study invariant manifolds of measure-valued solutions of the partial differential equation for geodesic flow of a pressureless fluid. These solutions describe interaction dynamics on lower-dimensional support sets; for example, curves, or filaments, of momentum in the plane. The 2+1 solutions we study are planar generalizations of the 1+1 peakon solutions of Camassa & Holm [1993] for shallow… ▽ More We study invariant manifolds of measure-valued solutions of the partial differential equation for geodesic flow of a pressureless fluid. These solutions describe interaction dynamics on lower-dimensional support sets; for example, curves, or filaments, of momentum in the plane. The 2+1 solutions we study are planar generalizations of the 1+1 peakon solutions of Camassa & Holm [1993] for shallow water solitons. As an example, we study the canonical Hamiltonian interaction dynamics of $N$ rotating concentric circles of peakons, whose solution manifold is $2N-$dimensional. Thus, the problem is reduced from infinite dimensions to a finite-dimensional, canonical, invariant manifold. The existence of this reduced solution manifold and many of its properties may be understood, by noticing that it is also the momentum map for the action of diffeomorphisms on the space of curves in the plane. We show both analytical and numerical results. △ Less

Submitted 5 December, 2003; originally announced December 2003.

Comments: 29 pages, 12 figures