Source author record

Kevin Zumbrun

Kevin Zumbrun 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

54works
8topics
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

54 published item(s)

preprint2022arXiv

Forward-modulated damping estimates and nonlocalized stability of periodic Lugiato-Lefever wave

In an interesting recent analysis, Haragus-Johnson-Perkins-de Rijk have shown modulational stability under localized perturbations of steady periodic solutions of the Lugiato-Lefever equation (LLE), in the process pointing out a difficulty in obtaining standard "nonlinear damping estimates" on modulated perturbation variables to control regularity of solutions. Here, we point out that in place of standard "inverse-modulated" damping estimates, one can alternatively carry out a damping estimate on the "forward-modulated" perturbation, noting that norms of forward- and inverse-modulated variables are equivalent modulo absorbable errors, thus recovering the classical argument structure of Johnson-Noble-Rodrigues-Zumbrun for parabolic systems. This observation seems of general use in situations of delicate regularity. Applied in the context of (LLE) it gives the stronger result of stability and asymptotic behavior with respect to nonlocalized perturbations.

preprint2022arXiv

Instantaneous smoothing and exponential decay of solutions for a degenerate evolution equation with application to Boltzmann's equation

We establish an instantaneous smoothing property for decaying solutions on the half-line $(0,+\infty)$ of certain degenerate Hilbert space-valued evolution equations arising in kinetic theory, including in particular the steady Boltzmann equation. Our results answer the two main open problems posed by Pogan and Zumbrun in their treatment of $H^1$ stable manifolds of such equations, showing that $L^2_{loc}$ solutions that remain sufficiently small in $L^\infty$ (i) decay exponentially, and (ii) are $C^\infty$ for $t>0$, hence lie eventually in the $H^1$ stable manifold constructed by Pogan and Zumbrun

preprint2022arXiv

Large-amplitude modulation of periodic traveling waves

We introduce a new approach to the study of modulation of high-frequency periodic wave patterns, based on pseudodifferential analysis, multi-scale expansion, and Kreiss symmetrizer estimates like those in hyperbolic and hyperbolic-parabolic boundary-value theory. Key ingredients are local Floquet transformation as a preconditioner removing large derivatives in the normal direction of background rapidly oscillating fronts and the use of the periodic Evans function of Gardner to connect spectral information on component periodic waves to block structure of the resulting approximately constant-coefficient resolvent ODEs. Our main result is bounded-time existence and valitidy to all orders of large-amplitude smooth modulations of planar periodic solutions of multi-D reaction diffusion systems in the high-frequency/small wavelength limit.

preprint2020arXiv

Stability of strong detonation waves for Majda's model with general ignition functions

For strong detonation waves of the inviscid Majda model, spectral stability was established by Jung and Yao for waves with step-type ignition functions, by a proof based largely on explicit knowledge of wave profiles. In the present work, we extend their stability results to strong detonation waves with more general ignition functions where explicit profiles are unknown. Our proof is based on reduction to a generalized Sturm-Liouville problem, similar to that used by Sukhtayev, Yang, and Zumbrun to study spectral stability of hydraulic shock profiles of the Saint-Venant equations.

preprint2020arXiv

Transverse bifurcation of viscous slow MHD shocks

We study by a combination of analytical and numerical Evans function techniques multi-D viscous and inviscid stability and associated transverse bifurcation of planar slow Lax MHD shocks in a channel with periodic boundary conditions. Notably, this includes the first multi-D numerical Evans function study for viscous MHD. Our results suggest that, rather than a planar shock, a nonplanar traveling wave with the same normal velocity is the typical mode of propagation in the slow Lax mode. Moreover, viscous and inviscid stability transitions appear to agree, answering (for this particular model and setting) an open question of Zumbrun and Serre.

preprint2016arXiv

Center manifolds for a class of degenerate evolution equations and existence of small amplitude kinetic shocks

We construct center manifolds for a class of degenerate evolution equations including the steady Boltzmann equation and related kinetic models, establishing in the process existence and behavior of small-amplitude kinetic shock and boundary layers. Notably, for Boltzmann's equation, we show that elements of the center manifold decay in velocity at near-Maxwellian rate, in accord with the formal Chapman-Enskog picture of near-equilibrium ow as evolution along the manifold of Maxwellian states, or Grad moment approximation via Hermite polynomials in velocity. Our analysis is from a classical dynamical systems point of view, with a number of interesting modifications to accommodate ill-posedness of the underlying evolution equation.

preprint2016arXiv

Derivatives of (Modified) Fredholm Determinants and Stability of Standing and Traveling Waves

Continuing a line of investigation initiated in [11] exploring the connections between Jost and Evans functions and (modified) Fredholm determinants of Birman-Schwinger type integral operators, we here examine the stability index, or sign of the first nonvanishing derivative at frequency zero of the characteristic determinant, an object that has found considerable use in the study by Evans function techniques of stability of standing and traveling wave solutions of partial differential equations (PDE) in one dimension. This leads us to the derivation of general perturbation expansions for analytically-varying modified Fredholm determinants of abstract operators. Our main conclusion, similarly in the analysis of the determinant itself, is that the derivative of the characteristic Fredholm determinant may be efficiently computed from first principles for integral operators with semi-separable integral kernels, which include in particular the general one-dimensional case, and for sums thereof, which latter possibility appears to offer applications in the multi-dimensional case. A second main result is to show that the multi-dimensional characteristic Fredholm determinant is the renormalized limit of a sequence of Evans functions defined in [23] on successive Galerkin subspaces, giving a natural extension of the one-dimensional results of [11] and answering a question of [27] whether this sequence might possibly converge (in general, no, but with renormalization, yes). Convergence is useful in practice for numerical error control and acceleration.

preprint2016arXiv

Invariant manifolds for a class of degenerate evolution equations and structure of kinetic shock layers

We describe recent results with A. Pogan developing dynamical systems tools for a class of degenerate evolution equations arising in kinetic theory, including the steady Boltzmann and BGK equations. These yield information on structure of large- and small-amplitude kinetic shocks, the first steps in a larger program toward time-evolutionary stability and asymptotic behavior

preprint2016arXiv

L infinity resolvent bounds for steady Boltzmann's equation

We derive lower bounds on the resolvent operator for the linearized steady Boltzmann equation over weighted L1 Banach spaces in velocity, comparable to those derived by Pogan & Zumbrun in an analogous weighted L2 Hilbert space setting. These show in particular that the operator norm of the resolvent kernel is unbounded in Lp(R) for all $1<p \leq \infty$, resolving an apparent discrepancy in behavior between the two settings suggested by previous work.

preprint2016arXiv

Note on the stability of viscous roll-waves

In this note, we announce a complete classification of stability of periodic roll-wave solutions of the viscous shallow-water equations, from their onset at Froude number $F\approx 2$ up to the infinite-Froude limit. For intermediate Froude numbers, we obtain numerically a particularly simple power-law relation between $F$ and the boundaries of the region of stable periods, that appears potentially useful in hydraulic engineering applications. In the asymptotic regime $F\to 2$ (onset), we provide an analytic expression of the stability boundaries whereas in the limit $F\to\infty$, we show that roll-waves are always unstable.

preprint2016arXiv

Numerical proof of stability of viscous shock profiles

We carry out the first rigorous numerical proof based on Evans function computations of stability of viscous shock profiles, for the system of isentropic gas dynamics with monatomic equation of state. We treat a selection of shock strengths ranging from the lower stability boundary of Mach number $\approx 1.86 $, below which profiles are known by energy estimates to be stable, to the upper stability boundary of $\approx 1669$, above which profiles are expected to be provable by rigorous asymptotic analysis to be stable. These results open the possibilities of: (i) automatic rigorous verification of stability or instability of individual shocks of general systems, and (ii) rigorous proof of stability of all shocks of particular systems.

preprint2016arXiv

Pointwise nonlinear stability of nonlocalized modulated periodic reaction-diffusion waves

In this paper, extending previous results of \cite{J1}, we obtain pointwise nonlinear stability of periodic traveling reaction-diffusion waves, assuming spectral linearized stability, under nonlocalized perturbations. More precisely, we establish pointwise estimate of nonlocalized modulational perturbation under a small initial perturbation consisting of a nonlocalized modulation plus a localized perturbation decaying algebraically.

preprint2016arXiv

Recent results on stability of planar detonations

We describe recent analytical and numerical results on stability and behavior of viscous and inviscid detonation waves obtained by dynamical systems/Evans function techniques like those used to study shock and reaction diffusion waves. In the first part, we give a broad description of viscous and inviscid results for 1D perturbations; in the second, we focus on inviscid high-frequency stability in multi-D and associated questions in turning point theory/WKB expansion.

preprint2016arXiv

Stable manifolds for a class of degenerate evolution equations and exponential decay of kinetic shocks

We construct stable manifolds for a class of degenerate evolution equations including the steady Boltzmann equation, establishing in the process exponential decay of associated kinetic shock and boundary layers to their limiting equilibrium states. Our analysis is from a classical dynamical systems point of view, but with a number of interesting modifications to accomodate ill-posedness with respect to the Cauchy problem of the underlying evolution equation.

preprint2015arXiv

Block-diagonalization of ODEs in the semiclassical limit and $C^ω$ vs. $C^\infty$ stationary phase

Motivated by issues in detonation stability, we study existence of block-diagonalizing transformations for ordinary differential semiclassical limit problems arising in the study of high-frequency eigenvalue problems. Our main results are to (i) establish existence of block-diagonalizing transformations in a neighborhood of infinity for analytic-coefficient ODE, and (ii) establish by a series of counterexample sharpness of hypotheses and conclusions on existence of block-diagonalizing transformations near a finite point. In particular, we show that, in general, bounded transformations exist only locally, answering a question posed by Wasow in the 1980's, and, under the minimal condition of spectral separation, for ODE with analytic rather than $C^\infty$ coefficients. The latter issue is connected with quantitative comparisons of $C^ω$ vs. $C^\infty$ stationary phase estimates

preprint2015arXiv

Periodic-coefficient damping estimates, and stability of large-amplitude roll waves in inclined thin film flow

A technical obstruction preventing the conclusion of nonlinear stability of large-Froude number roll waves of the St. Venant equations for inclined thin film flow is the "slope condition" of Johnson-Noble-Zumbrun, used to obtain pointwise symmetrizability of the linearized equations and thereby high-frequency resolvent bounds and a crucial H s nonlinear damping estimate. Numerically, this condition is seen to hold for Froude numbers 2 \textless{} F 3.5, but to fail for 3.5 F. As hydraulic engineering applications typically involve Froude number 3 F 5, this issue is indeed relevant to practical considerations. Here, we show that the pointwise slope condition can be replaced by an averaged version which holds always, thereby completing the nonlinear theory in the large-F case. The analysis has potentially larger interest as an extension to the periodic case of a type of weighted "Kawashima-type" damping estimate introduced in the asymptotically-constant coefficient case for the study of stability of large-amplitude viscous shock waves.

preprint2015arXiv

Stability of viscous St. Venant roll-waves: from onset to infinite-Froude number limit

We study the spectral stability of roll-wave solutions of the viscous St. Venant equations modeling inclined shallow-water flow, both at onset in the small-Froude number or "weakly unstable" limit $F\to 2^+$ and for general values of the Froude number $F$, including the limit $F\to +\infty$. In the former, $F\to 2^+$, limit, the shallow water equations are formally approximated by a Korteweg de Vries/Kuramoto-Sivashinsky (KdV-KS) equation that is a singular perturbation of the standard Korteweg de Vries (KdV) equation modeling horizontal shallow water flow. Our main analytical result is to rigorously validate this formal limit, showing that stability as $F\to 2^+$ is equivalent to stability of the corresponding KdV-KS waves in the KdV limit. Together with recent results obtained for KdV-KS by Johnson--Noble--Rodrigues--Zumbrun and Barker, this gives not only the first rigorous verification of stability for any single viscous St. Venant roll wave, but a complete classification of stability in the weakly unstable limit. In the remainder of the paper, we investigate numerically and analytically the evolution of the stability diagram as Froude number increases to infinity. Notably, we find transition at around $F=2.3$ from weakly unstable to different, large-$F$ behavior, with stability determined by simple power law relations. The latter stability criteria are potentially useful in hydraulic engineering applications, for which typically $2.5\leq F\leq 6.0$.

preprint2014arXiv

Entropy criteria and stability of extreme shocks: a remark on a paper of Leger and Vasseur

We show that a relative entropy condition recently shown by Leger and Vasseur to imply uniqueness and stable $L^2$ dependence on initial data of Lax 1- or $n$-shock solutions of an $n\times n$ system of hyperbolic conservation laws with convex entropy implies Lopatinski stability in the sense of Majda. This means in particular that Leger and Vasseur's relative entropy condition represents a considerable improvement over the standard entropy condition of decreasing shock strength and increasing entropy along forward Hugoniot curves, which, in a recent example exhibited by Barker, Freistühler and Zumbrun, was shown to fail to imply Lopatinski stability, even for systems with convex entropy. This observation bears also on the parallel question of existence, at least for small $BV$ or $H^s$ perturbations

preprint2014arXiv

O(2) Hopf bifurcation of viscous shock waves in a channel

Extending work of Texier and Zumbrun in the semilinear non-re ection symmetric case, we study O(2) transverse Hopf bifurcation, or \cellular instability," of viscous shock waves in a channel, for a class of quasilinear hyperbolic{parabolic systems including the equations of thermoviscoelasticity. The main difficulties are to (i) obtain Fr'echet differentiability of the time-T solution operator by appropriate hyperbolic{parabolic energy estimates, and (ii) handle O(2) symmetry in the absence of either center manifold reduction (due to lack of spectral gap) or (due to nonstandard quasilinear hyperbolic-parabolic form) the requisite framework for treatment by spatial dynamics on the space of time-periodic functions, the two standard treatments for this problem. The latter issue is resolved by Lyapunov{Schmidt reduction of the time-T map, yielding a four-dimensional problem with O(2) plus approximate S1 symmetry, which we treat \by hand" using direct Implicit Function Theorem arguments. The former is treated by balancing information obtained in Lagrangian coordinates with that from an augmented system. Interestingly, this argument does not apply to gas dynamics or magnetohydrodynamics (MHD), due to the infinite-dimensional family of Lagrangian symmetries corresponding to invariance under arbitrary volume-preserving diffeomorphisms.

preprint2013arXiv

High-frequency stability of multidimensional ZND detonations

The rigorous study of spectral stability of strong detonations was begun by Erpenbeck in the 1960s. Working with the Zeldovitch-von Neumann-Döring (ZND) model, he identified two fundamental classes of detonation profiles, referred to as those of decreasing (D) and increasing (I) type, which appeared to exhibit very different behavior with respect to high-frequency perturbations. Using a combination of rigorous and non-rigorous arguments, Erpenbeck concluded that type I detonations were unstable to some oscillatory perturbations for which the (vector) frequency was of arbitrarily large magnitude, while type D detonations were stable provided the frequency magnitude was sufficiently high. For type D detonations Erpenbeck's methods did not allow him to obtain a cutoff magnitude for stability that was \emph{uniform} with respect to frequency direction. Thus, he left open the question whether the cutoff magnitude for stability might approach $+\infty$ as certain frequency directions were approached. In this paper we show by quite different methods that for type D detonations there exists a uniform cutoff magnitude for stability independent of frequency direction. By reducing the search for unstable frequencies to a bounded frequency set, the uniform cutoff obtained here is a key step toward the rigorous validation of a number of results in the computational detonation literature. The main difficulty in the analysis is to treat "turning points at infinity".

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.

preprint2013arXiv

Nonlocalized modulation of periodic reaction diffusion waves: The Whitham equation

In a companion paper, we established nonlinear stability with detailed diffusive rates of decay of spectrally stable periodic traveling-wave solutions of reaction diffusion systems under small perturbations consisting of a nonlocalized modulation plus a localized perturbation. Here, we determine time-asymptotic behavior under such perturbations, showing that solutions consist to leading order of a modulation whose parameter evolution is governed by an associated Whitham averaged equation.

preprint2012arXiv

Behavior of periodic solutions of viscous conservation laws under localized and nonlocalized perturbations

We establish nonlinear stability and asymptotic behavior of traveling periodic waves of viscous conservation laws under localized perturbations or nonlocalized perturbations asymptotic to constant shifts in phase, showing that long-time behavior is governed by an associated second-order formal Whitham modulation system. A key point is to identify the way in which initial perturbations translate to initial data for this formal system, a task accomplished by detailed estimates on the linearized solution operator about the background wave. Notably, our approach gives both a common theoretical treatment and a complete classification in terms of "phase-coupling" or "-decoupling" of general systems of conservation or balance laws, encompassing cases that had previously been studied separately or not at all. At the same time, our refined description of solutions gives the new result of nonlinear asymptotic stability with respect to localized perturbations in the phase-decoupled case, further distinguishing behavior in the different cases. An interesting technical aspect of our analysis is that for systems of conservation laws the Whitham modulation description is of system rather than scalar form, as a consequence of which renormalization methods such as have been used to treat the reaction-diffusion case in general do not seem to apply.

preprint2012arXiv

Convex entropy, Hopf bifurcation, and viscous and inviscid shock stability

We consider by a combination of analytical and numerical techniques some basic questions regarding the relations between inviscid and viscous stability and existence of a convex entropy. Specifically, for a system possessing a convex entropy, in particular for the equations of gas dynamics with a convex equation of state, we ask: (i) can inviscid instability occur? (ii) can there occur viscous instability not detected by inviscid theory? (iii) can there occur the ---necessarily viscous--- effect of Hopf bifurcation, or "galloping instability"? and, perhaps most important from a practical point of view, (iv) as shock amplitude is increased from the (stable) weak-amplitude limit, can there occur a first transition from viscous stability to instability that is not detected by inviscid theory? We show that (i) does occur for strictly hyperbolic, genuinely nonlinear gas dynamics with certain convex equations of state, while (ii) and (iii) do occur for an artifically constructed system with convex viscosity-compatible entropy. We do not know of an example for which (iv) occurs, leaving this as a key open question in viscous shock theory, related to the principal eigenvalue property of Sturm Liouville and related operators. In analogy with, and partly proceeding close to, the analysis of Smith on (non-)uniqueness of the Riemann problem, we obtain convenient criteria for shock (in)stability in the form of necessary and sufficient conditions on the equation of state.

preprint2012arXiv

Nonlinear modulational stability of periodic traveling-wave solutions of the generalized Kuramoto-Sivashinsky equation

In this paper we consider the spectral and nonlinear stability of periodic traveling wave solutions of a generalized Kuramoto-Sivashinsky equation. In particular, we resolve the long-standing question of nonlinear modulational stability by demonstrating that spectrally stable waves are nonlinearly stable when subject to small localized (integrable) perturbations. Our analysis is based upon detailed estimates of the linearized solution operator, which are complicated by the fact that the (necessarily essential) spectrum of the associated linearization intersects the imaginary axis at the origin. We carry out a numerical Evans function study of the spectral problem and find bands of spectrally stable periodic traveling waves, in close agreement with previous numerical studies of Frisch-She-Thual, Bar-Nepomnyashchy, Chang-Demekhin-Kopelevich, and others carried out by other techniques. We also compare predictions of the associated Whitham modulation equations, which formally describe the dynamics of weak large scale perturbations of a periodic wave train, with numerical time evolution studies, demonstrating their effectiveness at a practical level. For the reader's convenience, we include in an appendix the corresponding treatment of the Swift-Hohenberg equation, a nonconservative counterpart of the generalized Kuramoto-Sivashinsky equation for which the nonlinear stability analysis is considerably simpler, together with numerical Evans function analyses extending spectral stability analyses of Mielke and Schneider.

preprint2012arXiv

Quasi-gradient systems, modulational dichotomies, and stability of spatially periodic patterns

Extending the approach of Grillakis-Shatah-Strauss, Bronski-Johnson-Kapitula, and others for Hamiltonian systems, we explore relations between the constrained variational problem $\min_{X:C(X)=c_0} \mathcal{E}(X)$, $c_0\in \RM^r$, and stability of solutions of a class of degenerate "quasi-gradient" systems $dX/dt=-M(X)\nabla \mathcal{E}(X)$ admitting constraints, including Cahn-Hilliard equations, one- and multi-dimensional viscoelasticity, and coupled conservation law-reaction diffusion systems arising in chemotaxis and related settings. Using the relation between variational stability and the signature of $\partial c/\partial ω\in \R^{r\times r}$, where $c(ω)=C(X^*_ω)\in \RM^r$ denote the values of the imposed constraints and $ω\in \RM^r$ the associated Lagrange multipliers at a critical point $X^*_ω$, we obtain as in the Hamiltonian case a general criterion for co-periodic stability of periodic waves, illuminating and extending a number of previous results obtained by direct Evans function techniques. More interestingly, comparing the form of the Jacobian arising in the co-periodic theory to Jacobians arising in the formal Whitham equations associated with modulation, we recover and substantially generalize a previously mysterious "modulational dichotomy" observed in special cases by Oh-Zumbrun and Howard, showing that co-periodic and sideband stability are incompatible. In particular, we both illuminate and extend to general viscosity/strain-gradient effects and multidimensional deformations the result of Oh-Zumbrun of universal modulational instability of periodic solutions of the equations of viscoelasticity with strain-gradient effects, considered as functions on the whole line. Likewise, we generalize to multi-dimensions corresponding results of Howard on periodic solutions of Cahn-Hilliard equations.

preprint2012arXiv

Spectral stability of periodic wave trains of the Korteweg-de Vries/Kuramoto-Sivashinsky equation in the Korteweg-de Vries limit

We study the spectral stability of a family of periodic wave trains of the Korteweg-de Vries/Kuramoto-Sivashinsky equation $ \partial_t v+v\partial_x v+\partial_x^3 v+δ(\partial_x^2 v +\partial_x^4 v)=0$, $δ>0$, in the Korteweg-de Vries limit $δ\to 0$, a canonical limit describing small-amplitude weakly unstable thin film flow. More precisely, we carry out a rigorous singular perturbation analysis reducing the problem to the evaluation for each Bloch parameter $ξ\in [0,2π]$ of certain elliptic integrals derived formally (on an incomplete set of frequencies/Bloch parameters, hence as necessary conditions for stability) and numerically evaluated by Bar and Nepomnyashchy \cite{BN}, thus obtaining, up to machine error, complete conclusions about stability. The main technical difficulty is in treating the large-frequency and small Bloch-parameter regimes not studied by Bar and Nepomnyashchy \cite{BN}, which requires techniques rather different from classical Fenichel-type analysis. The passage from small-$δ$ to small-$ξ$ behavior is particularly interesting, using in an essential way an analogy with hyperbolic relaxation at the level of the Whitham modulation equations.

preprint2012arXiv

Viscous boundary layers in hyperbolic-parabolic systems with Neumann boundary conditions

We initiate the study of noncharacteristic boundary layers in hyperbolic-parabolic problems with Neumann boundary conditions. More generally, we study boundary layers with mixed Dirichlet--Neumann boundary conditions where the number of Dirichlet conditions is fewer than the number of hyperbolic characteristic modes entering the domain, that is, the number of boundary conditions needed to specify an outer hyperbolic solution. We have shown previously that this situation prevents the usual WKB approximation involving an outer solution with pure Dirichlet conditions. It also rules out the usual maximal estimates for the linearization of the hyperbolic-parabolic problem about the boundary layer. Here we show that for linear, constant-coefficient, hyperbolic-parabolic problems one obtains a reduced hyperbolic problem satisfying Neumann or mixed Dirichlet--Neumann rather than Dirichlet boundary conditions. When this hyperbolic problem can be solved, a unique formal boundary-layer expansion can be constructed. In the extreme case of pure Neumann conditions and totally incoming characteristics, we carry out a full analysis of the quasilinear case, obtaining a boundary-layer approximation to all orders with a rigorous error analysis. As a corollary we characterize the small viscosity limit for this problem. The analysis shows that although the associated linearized hyperbolic and hyperbolic--parabolic problems do not satisfy the usual maximal estimates for Dirichlet conditions, they do satisfy analogous versions with losses.

preprint2011arXiv

Nash--Moser iteration and singular perturbations

We present a simple and easy-to-use Nash--Moser iteration theorem tailored for singular perturbation problems admitting a formal asymptotic expansion or other family of approximate solutions depending on a parameter $\eps\to 0.$ The novel feature is to allow loss of powers of $\eps$ as well as the usual loss of derivatives in the solution operator for the associated linearized problem. We indicate the utility of this theorem by describing sample applications to (i) systems of quasilinear Schrödinger equations, and (ii) existence of small-amplitude profiles of quasilinear relaxation systems.

preprint2011arXiv

Nonlocalized modulation of periodic reaction diffusion waves: Nonlinear stability

By a refinement of the technique used by Johnson and Zumbrun to show stability under localized perturbations, we show that spectral stability implies nonlinear modulational stability of periodic traveling-wave solutions of reaction diffusion systems under small perturbations consisting of a nonlocalized modulation plus a localized perturbation. The main new ingredient is a detailed analysis of linear behavior under modulational data $\bar u'(x)h_0(x)$, where $\bar u$ is the background profile and $h_0$ is the initial modulation

preprint2011arXiv

The Erpenbeck high frequency instability theorem for ZND detonations

The rigorous study of spectral stability for strong detonations was begun by J.J. Erpenbeck in [Er1]. Working with the Zeldovitch-von Neumann-Döring (ZND) model, which assumes a finite reaction rate but ignores effects like viscosity corresponding to second order derivatives, he used a normal mode analysis to define a stability function $V(τ,\eps)$ whose zeros in $\Re τ>0$ correspond to multidimensional perturbations of a steady detonation profile that grow exponentially in time. Later in a remarkable paper [Er3] he provided strong evidence, by a combination of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of $V$ exist for perturbations of sufficiently large transverse wavenumber $\eps$, even when the von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly twenty years later) by Majda. In spite of a great deal of later numerical work devoted to computing the zeros of $V(τ,\eps)$, the paper \cite{Er3} remains the only work we know of that presents a detailed and convincing theoretical argument for detecting them. The analysis in [Er3] points the way toward, but does not constitute, a mathematical proof that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved in [Er3] and provide proofs, together with certain simplifications and extensions, of the main conclusions about stability and instability of detonations contained in that paper. The main mathematical problem, and our principal focus here, is to determine the precise asymptotic behavior as $\eps\to \infty$ of solutions to a linear system of ODEs in $x$, depending on $\eps$ and a complex frequency $τ$ as parameters, with turning points $x_*$ on the half-line $[0,\infty)$.

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.

preprint2010arXiv

$2$-modified characteristic Fredholm determinants, Hill's method, and the periodic Evans function of Gardne

Using the relation established by Johnson--Zumbrun between Hill's method of aproximating spectra of periodic-coefficient ordinary differential operators and a generalized periodic Evans function given by the $2$-modified characteristic Fredholm determinant of an associated Birman--Schwinger system, together with a Volterra integral computation introduced by Gesztesy--Makarov, we give an explicit connection between the generalized Birman--Schwinger-type periodic Evans function and the standard Jost function-type periodic Evans function defined by Gardner in terms of the fundamental solution of the eigenvalue equation written as a first-order system. This extends to a large family of operators the results of Gesztesy--Makarov for scalar Schrödinger operators and of Gardner for vector-valued second-order elliptic operators, in particular recovering by independent argument the fundamental result of Gardner that the zeros of the Evans function agree in location and (algebraic) multiplicity with the periodic eigenvalues of the associated operator

preprint2010arXiv

A numerical stability investigation of strong ZND detonations for Majda's model

We carry out a systematic numerical stability analysis of ZND detonations of Majda's model with Arrhenius-type ignition function, a simplified model for reacting flow, as heat release and activation energy are varied. Our purpose is, first, to answer a question of Majda whether oscillatory instabilities can occur for high activation energies as in the full reacting Euler equations, and, second, to test the efficiency of various versions of a numerical eigenvalue-finding scheme suggested by Humpherys and Zumbrun against the standard method of Lee and Stewart. Our results suggest that instabilities do not occur for Majda's model with Arrhenius-type ignition function, nor with a modified Arrhenius-type ignition function suggested by Lyng--Zumbrun, even in the high-activation energy limit. We find that the algorithm of Humpherys--Zumbrun is in the context of Majda's model $100$-$1,000$ times faster than the one described in the classical work of Lee and Stewart and $1$-$10$ times faster than an optimized version of the Lee--Stewart algorithm using an adaptive-mesh ODE solver

preprint2010arXiv

Convergence of Hill's method for nonselfadjoint operators

By the introduction of a generalized Evans function defined by an appropriate 2-modified Fredholm determinant, we give a simple proof of convergence in location and multiplicity of Hill's method for numerical approximation of spectra of periodic-coefficient ordinary differential operators. Our results apply to operators of nondegenerate type, under the condition that the principal coefficient matrix be symmetric positive definite (automatically satisfied in the scalar case). Notably, this includes a large class of nonselfadjoint operators, which were previously not treated. The case of general coefficients depends on an interesting operator-theoretic question regarding properties of Toeplitz matrices.

preprint2010arXiv

Existence and stability of steady states of a reaction convection diffusion equation modeling microtubule formation

We generalize the Dogterom-Leibler model for microtubule dynamics [DL] to the case where the rates of elongation as well as the lifetimes of the elongating and shortening phases are a function of GTP-tubulin concentration. We study also the effect of nucleation rate in the form of a damping term which leads to new steady-states. For this model, we study existence and stability of steady states satisfying the boundary conditions at x = 0. Our stability analysis introduces numerical and analytical Evans function computations as a new mathematical tool in the study of microtubule dynamics.

preprint2010arXiv

Existence and stability of viscoelastic shock profiles

We investigate existence and stability of viscoelastic shock profiles for a class of planar models including the incompressible shear case studied by Antman and Malek-Madani. We establish that the resulting equations fall into the class of symmetrizable hyperbolic--parabolic systems, hence spectral stability implies linearized and nonlinear stability with sharp rates of decay. The new contributions are treatment of the compressible case, formulation of a rigorous nonlinear stability theory, including verification of stability of small-amplitude Lax shocks, and the systematic incorporation in our investigations of numerical Evans function computations determining stability of large-amplitude and or nonclassical type shock profiles.

preprint2010arXiv

High-frequency asymptotics and 1-D stability of ZND detonations in the small-heat release and high-overdrive limits

We establish one-dimensional spectral, or "normal modes", stability of ZND detonations in the small heat release limit and the related high overdrive limit with heat release and activation energy held fixed, verifying numerical observations of Erpenbeck in the 1960s. The key technical points are a strategic rescaling of parameters converting the infinite overdrive limit to a finite, regular perturbation problem, and a careful high-frequency analysis depending uniformly on model parameters. The latter recovers the important result of high-frequency stability established by Erpenbeck by somewhat different techniques. Notably, the techniques used here yield quantitative estimates well suited for numerical stability investigation.

preprint2010arXiv

Metastability of solitary roll wave solutions of the St. Venant equations with viscosity

We study by a combination of numerical and analytical Evans function techniques the stability of solitary wave solutions of the St. Venant equations for viscous shallow-water flow down an incline, and related models. Our main result is to exhibit examples of metastable solitary waves for the St. Venant equations, with stable point spectrum indicating coherence of the wave profile but unstable essential spectrum indicating oscillatory convective instabilities shed in its wake. We propose a mechanism based on ``dynamic spectrum'' of the wave profile, by which a wave train of solitary pulses can stabilize each other by de-amplification of convective instabilities as they pass through successive waves. We present numerical time evolution studies supporting these conclusions, which bear also on the possibility of stable periodic solutions close to the homoclinic. For the closely related viscous Jin-Xin model, by contrast, for which the essential spectrum is stable, we show using the stability index of Gardner--Zumbrun that solitary wave pulses are always exponentially unstable, possessing point spectra with positive real part.

preprint2010arXiv

Nonlinear stability of periodic traveling wave solutions of systems of viscous conservation laws in the generic case

Extending previous results of Oh--Zumbrun and Johnson--Zumbrun, we show that spectral stability implies linearized and nonlinear stability of spatially periodic traveling-wave solutions of viscous systems of conservation laws for systems of generic type, removing a restrictive assumption that wave speed be constant to first order along the manifold of nearby periodic solutions.

preprint2010arXiv

Nonlinear stability of periodic traveling wave solutions of viscous conservation laws in dimensions one and two

Extending results of Oh and Zumbrun in dimensions $d\ge 3$, we establish nonlinear stability and asymptotic behavior of spatially-periodic traveling-wave solutions of viscous systems of conservation laws in critical dimensions $d=1,2$, under a natural set of spectral stability assumptions introduced by Schneider in the setting of reaction diffusion equations. The key new steps in the analysis beyond that in dimensions $d\ge 3$ are a refined Green function estimate separating off translation as the slowest decaying linear mode and a novel scheme for detecting cancellation at the level of the nonlinear iteration in the Duhamel representation of a modulated periodic wave.

preprint2010arXiv

Nonlinear stability of spatially-periodic traveling-wave solutions of systems of reaction diffusion equations

Using spatial domain techniques developed by the authors and Myunghyun Oh in the context of parabolic conservation laws, we establish under a natural set of spectral stability conditions nonlinear asymptotic stability with decay at Gaussian rate of spatially periodic traveling-waves of systems of reaction diffusion equations. In the case that wave-speed is identically zero for all periodic solutions, we recover and slightly sharpen a well-known result of Schneider obtained by renormalization/Bloch transform techniques; by the same arguments, we are able to treat the open case of nonzero wave-speeds to which Schneider's renormalization techniques do not appear to apply

preprint2010arXiv

Nonlinear stability of viscous roll waves

Extending results of Oh--Zumbrun and Johnson--Zumbrun for parabolic conservation laws, we show that spectral stability implies nonlinear stability for spatially periodic viscous roll wave solutions of the one-dimensional St. Venant equations for shallow water flow down an inclined ramp. The main new issues to be overcome are incomplete parabolicity and the nonconservative form of the equations, which leads to undifferentiated quadratic source terms that cannot be handled using the estimates of the conservative case. The first is resolved by treating the equations in the more favorable Lagrangian coordinates, for which one can obtain large-amplitude nonlinear damping estimates similar to those carried out by Mascia--Zumbrun in the related shock wave case, assuming only symmetrizability of the hyperbolic part. The second is resolved by the observation that, similarly as in the relaxation and detonation cases, sources occurring in nonconservative components experience greater than expected decay, comparable to that experienced by a differentiated source.

preprint2010arXiv

On the Modulation Equations and Stability of Periodic GKdV Waves via Bloch Decompositions

In this paper, we complement recent results of Bronski and Johnson and of Johnson and Zumbrun concerning the modulational stability of spatially periodic traveling wave solutions of the generalized Korteweg-de Vries equation. In this previous work it was shown by rigorous Evans function calculations that the formal slow modulation approximation resulting in the Whitham system accurately describes the spectral stability to long wavelength perturbations. Here, we reproduce this result without reference to the Evans function by using direct Bloch-expansion methods and spectral perturbation analysis. This approach has the advantage of applying also in the more general multi-periodic setting where no conveniently computable Evans function is yet devised. In particular, we complement the picture of modulational stability described by Bronski and Johnson by analyzing the projectors onto the total eigenspace bifurcating from the origin in a neighborhood of the origin and zero Floquet parameter. We show the resulting linear system is equivalent, to leading order and up to conjugation, to the Whitham system and that, consequently, the characteristic polynomial of this system agrees (to leading order) with the linearized dispersion relation derived through Evans function calculation.

preprint2010arXiv

Stability of Periodic Kuramoto-Sivashinsky Waves

In this note, we announce a general result resolving the long-standing question of nonlinear modulational stability, or stability with respect to localized perturbations, of periodic traveling-wave solutions of the generalized Kuramoto-Sivashinsky equation, establishing that spectral modulational stability, defined in the standard way, implies nonlinear modulational stability with sharp rates of decay. The approach extends readily to other second- and higher-order parabolic equations, for example, the Cahn-Hilliard equation or more general thin film models.

preprint2010arXiv

Transition to longitudinal instability of detonation waves is generically associated with Hopf bifurcation to time-periodic galloping solutions

We show that transition to longitudinal instability of strong detonation solutions of reactive compressible Navier--Stokes equations is generically associated with Hopf bifurcation to nearby time-periodic "galloping", or "pulsating", solutions, in agreement with physical and numerical observation. In the process, we determine readily numerically verifiable stability and bifurcation conditions in terms of an associated Evans function, and obtain the first complete nonlinear stability result for strong detonations of the reacting Navier--Stokes equations, in the limit as amplitude (hence also heat release) goes to zero. The analysis is by pointwise semigroup techniques introduced by the authors and collaborators in previous works.

preprint2010arXiv

Transverse Instability of Periodic Traveling Waves in the Generalized Kadomtsev-Petviashvili Equation

In this paper, we investigate the spectral instability of periodic traveling wave solutions of the generalized Korteweg-de Vries equation to long wavelength transverse perturbations in the generalized Kadomtsev-Petviashvili equation. By analyzing high and low frequency limits of the appropriate periodic Evans function, we derive an orientation index which yields sufficient conditions for such an instability to occur. This index is geometric in nature and applies to arbitrary periodic traveling waves with minor smoothness and convexity assumptions on the nonlinearity. Using the integrable structure of the ordinary differential equation governing the traveling wave profiles, we are then able to calculate the resulting orientation index for the elliptic function solutions of the Korteweg-de Vries and modified Korteweg-de Vries equations.

preprint2010arXiv

Whitham Averaged Equations and Modulational Stability of Periodic Traveling Waves of a Hyperbolic-Parabolic Balance Law

In this note, we report on recent findings concerning the spectral and nonlinear stability of periodic traveling wave solutions of hyperbolic-parabolic systems of balance laws, as applied to the St. Venant equations of shallow water flow down an incline. We begin by introducing a natural set of spectral stability assumptions, motivated by considerations from the Whitham averaged equations, and outline the recent proof yielding nonlinear stability under these conditions. We then turn to an analytical and numerical investigation of the verification of these spectral stability assumptions. While spectral instability is shown analytically to hold in both the Hopf and homoclinic limits, our numerical studies indicates spectrally stable periodic solutions of intermediate period. A mechanism for this moderate-amplitude stabilization is proposed in terms of numerically observed "metastability" of the the limiting homoclinic orbits.

preprint2009arXiv

Conditional stability of unstable viscous shock waves in compressible gas dynamics and MHD

Extending our previous work in the strictly parabolic case, we show that a linearly unstable Lax-type viscous shock solution of a general quasilinear hyperbolic--parabolic system of conservation laws possesses a translation-invariant center stable manifold within which it is nonlinearly orbitally stable with respect to small $L^1\cap H^3$ perturbations, converging time-asymptotically to a translate of the unperturbed wave. That is, for a shock with $p$ unstable eigenvalues, we establish conditional stability on a codimension-$p$ manifold of initial data, with sharp rates of decay in all $L^p$. For $p=0$, we recover the result of unconditional stability obtained by Mascia and Zumbrun. The main new difficulty in the hyperbolic--parabolic case is to construct an invariant manifold in the absence of parabolic smoothing.

preprint2009arXiv

Stability of radiative shock profiles for hyperbolic-elliptic coupled systems

Extending previous work with Lattanzio and Mascia on the scalar (in fluid-dynamical variables) Hamer model for a radiative gas, we show nonlinear orbital asymptotic stability of small-amplitude shock profiles of general systems of coupled hyperbolic--eliptic equations of the type modeling a radiative gas, that is, systems of conservation laws coupled with an elliptic equation for the radiation flux, including in particular the standard Euler--Poisson model for a radiating gas. The method is based on the derivation of pointwise Green function bounds and description of the linearized solution operator, with the main difficulty being the construction of the resolvent kernel in the case of an eigenvalue system of equations of degenerate type. Nonlinear stability then follows in standard fashion by linear estimates derived from these pointwise bounds, combined with nonlinear-damping type energy estimates

preprint2009arXiv

The refined inviscid stability condition and cellular instability of viscous shock waves

Combining work of Serre and Zumbrun, Benzoni-Gavage, Serre, and Zumbrun, and Texier and Zumbrun, we propose as a mechanism for the onset of cellular instability of viscous shock and detonation waves in a finite-cross-section duct the violation of the refined planar stability condition of Zumbrun--Serre, a viscous correction of the inviscid planar stability condition of Majda. More precisely, we show for a model problem involving flow in a rectangular duct with artificial periodic boundary conditions that transition to multidimensional instability through violation of the refined stability condition of planar viscous shock waves on the whole space generically implies for a duct of sufficiently large cross-section a cascade of Hopf bifurcations involving more and more complicated cellular instabilities. The refined condition is numerically calculable as described in Benzoni-Gavage--Serre-Zumbrun.

preprint2008arXiv

Existence and stability of noncharacteristic boundary-layers for the compressible Navier-Stokes and viscous MHD equations

For a general class of hyperbolic-parabolic systems including the compressible Navier-Stokes and compressible MHD equations, we prove existence and stability of noncharacteristic viscous boundary layers for a variety of boundary conditions including classical Navier-Stokes boundary conditions. Our first main result, using the abstract framework established by the authors in the companion work \cite{GMWZ6}, is to show that existence and stability of arbitrary amplitude exact boundary-layer solutions follow from a uniform spectral stability condition on layer profiles that is expressible in terms of an Evans function (uniform Evans stability). Whenever this condition holds we give a rigorous description of the small viscosity limit as the solution of a hyperbolic problem with "residual" boundary conditions. Our second is to show that uniform Evans stability for small-amplitude layers is equivalent to Evans stability of the limiting constant layer, which in turn can be checked by a linear-algebraic computation. Finally, for a class of symmetric-dissipative systems including the physical examples mentioned above, we carry out energy estimates showing that constant (and thus small-amplitude) layers always satisfy uniform Evans stability. This yields existence of small-amplitude multi-dimensional boundary layers for the compressible Navier-Stokes and MHD equations. For both equations these appear to be the first such results in the compressible case.

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.

preprint2003arXiv

A stability index for detonation waves in Majda's model for reacting flow

Using Evans function techniques, we develop a stability index for weak and strong detonation waves analogous to that developed for shock waves in [GZ,BSZ], yielding useful necessary conditions for stability. Here, we carry out the analysis in the context of the Majda model, a simplified model for reacting flow; the method is extended to the full Navier-Stokes equations of reacting flow in [Ly,LyZ]. The resulting stability condition is satisfied for all nondegenerate, i.e., spatially exponentially decaying, weak and strong detonations of the Majda model in agreement with numerical experiments of [CMR] and analytical results of [Sz,LY] for a related model of Majda and Rosales. We discuss also the role in the ZND limit of degenerate, subalgebraically decaying weak detonation and (for a modified, ``bump-type'' ignition function) deflagration profiles, as discussed in [GS.1-2] for the full equations.