Source author record

Margaret Beck

Margaret Beck appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Localized radial roll patterns in higher space dimensions

Localized roll patterns are structures that exhibit a spatially periodic profile in their center. When following such patterns in a system parameter in one space dimension, the length of the spatial interval over which these patterns resemble a periodic profile stays either bounded, in which case branches form closed bounded curves ("isolas"), or the length increases to infinity so that branches are unbounded in function space ("snaking"). In two space dimensions, numerical computations show that branches of localized rolls exhibit a more complicated structure in which both isolas and snaking occur. In this paper, we analyse the structure of branches of localized radial roll solutions in dimension 1+$\varepsilon$, with $0<\varepsilon\ll1$, through a perturbation analysis. Our analysis sheds light on some of the features visible in the planar case.

preprint2020arXiv

A dynamical approach to semilinear elliptic equations

A characterization of a semilinear elliptic partial differential equation (PDE) on a bounded domain in $\mathbb{R}^n$ is given in terms of an infinite-dimensional dynamical system. The dynamical system is on the space of boundary data for the PDE. This is a novel approach to elliptic problems that enables the use of dynamical systems tools in studying the corresponding PDE. The dynamical system is ill-posed, meaning solutions do not exist forwards or backwards in time for generic initial data. We offer a framework in which this ill-posed system can be analyzed. This can be viewed as generalizing the theory of spatial dynamics, which applies to the case of an infinite cylindrical domain.

preprint2020arXiv

Exponential dichotomies for elliptic PDE on radial domains

It was recently shown by the authors that a semilinear elliptic equation can be represented as an infinite-dimensional dynamical system in terms of boundary data on a shrinking one-parameter family of domains. The resulting system is ill-posed, in the sense that solutions do not typically exist forward or backward in time. In this paper we consider a radial family of domains and prove that the linearized system admits an exponential dichotomy, with the unstable subspace corresponding to the boundary data of weak solutions to the linear PDE. This generalizes the spatial dynamics approach, which applies to infinite cylindrical (channel) domains, and also generalizes previous work on radial domains as we impose no symmetry assumptions on the equation or its solutions.

preprint2020arXiv

Selection of quasi-stationary states in the stochastically forced Navier-Stokes equation on the torus

The stochastically forced vorticity equation associated with the two dimensional incompressible Navier-Stokes equation on $D_δ:=[0,2πδ]\times [0,2π]$ is considered for $δ\approx 1$, periodic boundary conditions, and viscocity $0<ν\ll 1$. An explicit family of quasi-stationary states of the deterministic vorticity equation is known to play an important role in the long-time evolution of solutions both in the presence of and without noise. Recent results show the parameter $δ$ plays a central role in selecting which of the quasi-stationary states is most important. In this paper, we aim to develop a finite dimensional model that captures this selection mechanism for the stochastic vorticity equation. This is done by projecting the vorticity equation in Fourier space onto a center manifold corresponding to the lowest eight Fourier modes. Through Monte Carlo simulation, the vorticity equation and the model are shown to be in agreement regarding key aspects of the long-time dynamics. Following this comparison, perturbation analysis is performed on the model via averaging and homogenization techniques to determine the leading order dynamics for statistics of interest for $δ\approx1$.

preprint2015arXiv

Analysis of enhanced diffusion in Taylor dispersion via a model problem

We consider a simple model of the evolution of the concentration of a tracer, subject to a background shear flow by a fluid with viscosity $ν\ll 1$ in an infinite channel. Taylor observed in the 1950's that, in such a setting, the tracer diffuses at a rate proportional to $1/ν$, rather than the expected rate proportional to $ν$. We provide a mathematical explanation for this enhanced diffusion using a combination of Fourier analysis and center manifold theory. More precisely, we show that, while the high modes of the concentration decay exponentially, the low modes decay algebraically, but at an enhanced rate. Moreover, the behavior of the low modes is governed by finite-dimensional dynamics on an appropriate center manifold, which corresponds exactly to diffusion by a fluid with viscosity proportional to $1/ν$.

preprint2013arXiv

Computing the Maslov index for large systems

We address the problem of computing the Maslov index for large linear symplectic systems on the real line. The Maslov index measures the signed intersections (with a given reference plane) of a path of Lagrangian planes. The natural chart parameterization for the Grassmannian of Lagrangian planes is the space of real symmetric matrices. Linear system evolution induces a Riccati evolution in the chart. For large order systems this is a practical approach as the computational complexity is quadratic in the order. The Riccati solutions, however, also exhibit singularites (which are traversed by changing charts). Our new results involve characterizing these Riccati singularities and two trace formulae for the Maslov index as follows. First, we show that the number of singular eigenvalues of the symmetric chart representation equals the dimension of intersection with the reference plane. Second, the Cayley map is a diffeomorphism from the space of real symmetric matrices to the manifold of unitary symmetric matrices. We show the logarithm of the Cayley map equals the arctan map (modulo 2i) and its trace measures the angle of the Langrangian plane to the reference plane. Third, the Riccati flow under the Cayley map induces a flow in the manifold of unitary symmetric matrices. Using the natural unitary action on this manifold, we pullback the flow to the unitary Lie algebra and monitor its trace. This avoids singularities, and is a natural robust procedure. We demonstrate the effectiveness of these approaches by applying them to a large eigenvalue problem. We also discuss the extension of the Maslov index to the infinite dimensional case.

preprint2013arXiv

Nonlinear stability of source defects in the complex Ginzburg-Landau equation

In an appropriate moving coordinate frame, source defects are time-periodic solutions to reaction-diffusion equations that are spatially asymptotic to spatially periodic wave trains whose group velocities point away from the core of the defect. In this paper, we rigorously establish nonlinear stability of spectrally stable source defects in the complex Ginzburg-Landau equation. Due to the outward transport at the far field, localized perturbations may lead to a highly non-localized response even on the linear level. To overcome this, we first investigate in detail the dynamics of the solution to the linearized equation. This allows us to determine an approximate solution that satisfies the full equation up to and including quadratic terms in the nonlinearity. This approximation utilizes the fact that the non-localized phase response, resulting from the embedded zero eigenvalues, can be captured, to leading order, by the nonlinear Burgers equation. The analysis is completed by obtaining detailed estimates for the resolvent kernel and pointwise estimates for the Green's function, which allow one to close a nonlinear iteration scheme.

preprint2012arXiv

Metastability and rapid convergence to quasi-stationary bar states for the 2D Navier-Stokes Equations

Quasi-stationary, or metastable, states play an important role in two-dimensional turbulent fluid flows where they often emerge on time-scales much shorter than the viscous time scale, and then dominate the dynamics for very long time intervals. In this paper we propose a dynamical systems explanation of the metastability of an explicit family of solutions, referred to as bar states, of the two-dimensional incompressible Navier-Stokes equation on the torus. These states are physically relevant because they are associated with certain maximum entropy solutions of the Euler equations, and they have been observed as one type of metastable state in numerical studies of two-dimensional turbulence. For small viscosity (high Reynolds number), these states are quasi-stationary in the sense that they decay on the slow, viscous timescale. Linearization about these states leads to a time-dependent operator. We show that if we approximate this operator by dropping a higher-order, non-local term, it produces a decay rate much faster than the viscous decay rate. We also provide numerical evidence that the same result holds for the full linear operator, and that our theoretical results give the optimal decay rate in this setting.

preprint2011arXiv

Toward nonlinear stability of sources via a modified Burgers equation

Coherent structures are solutions to reaction-diffusion systems that are time-periodic in an appropriate moving frame and spatially asymptotic at $x=\pm\infty$ to spatially periodic travelling waves. This paper is concerned with sources which are coherent structures for which the group velocities in the far field point away from the core. Sources actively select wave numbers and therefore often organize the overall dynamics in a spatially extended system. Determining their nonlinear stability properties is challenging as localized perturbations may lead to a non-localized response even on the linear level due to the outward transport. Using a modified Burgers equation as a model problem that captures some of the essential features of coherent structures, we show how this phenomenon can be analysed and nonlinear stability be established in this simpler context.

preprint2008arXiv

Nonlinear stability of time-periodic viscous shocks

In order to understand the nonlinear stability of many types of time-periodic travelling waves on unbounded domains, one must overcome two main difficulties: the presence of embedded neutral eigenvalues and the time-dependence of the associated linear operator. This problem is studied in the context of time-periodic Lax shocks in systems of viscous conservation laws. Using spatial dynamics and a decomposition into separate Floquet eigenmodes, it is shown that the linear evolution for the time-dependent operator can be represented using a contour integral similar to that of the standard time-independent case. By decomposing the resulting Green's distribution, the leading order behavior associated with the embedded eigenvalues is extracted. Sharp pointwise bounds are then obtained, which are used to prove that time-periodic Lax shocks are linearly and nonlinearly stable under the necessary conditions of spectral stability and minimal multiplicity of the translational eigenvalues. The latter conditions hold, for example, for small-oscillation time-periodic waves that emerge through a supercritical Hopf bifurcation from a family of time-independent Lax shocks of possibly large amplitude.

preprint2008arXiv

Using global invariant manifolds to understand metastability in Burgers equation with small viscosity

The large-time behavior of solutions to Burgers equation with small viscosity is described using invariant manifolds. In particular, a geometric explanation is provided for a phenomenon known as metastability, which in the present context means that solutions spend a very long time near the family of solutions known as diffusive N-waves before finally converging to a stable self-similar diffusion wave. More precisely, it is shown that in terms of similarity, or scaling, variables in an algebraically weighted $L^2$ space, the self-similar diffusion waves correspond to a one-dimensional global center manifold of stationary solutions. Through each of these fixed points there exists a one-dimensional, global, attractive, invariant manifold corresponding to the diffusive N-waves. Thus, metastability corresponds to a fast transient in which solutions approach this "metastable" manifold of diffusive N-waves, followed by a slow decay along this manifold, and, finally, convergence to the self-similar diffusion wave.