Researcher profile

Mark Williams

Mark Williams contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2022arXiv

Reflection of conormal pulse solutions to large variable-coefficient semilinear hyperbolic systems

We provide a rigorous justication of nonlinear geometric optics expansions for reflecting \emph{pulses} in space dimensions $n>1$. The pulses arise as solutions to variable coefficient semilinear first-order hyperbolic systems. The justification applies to $N\times N$ systems with $N$ interacting pulses which depend on phases that may be nonlinear. The \emph{coherence} assumption made in a number of earlier works is dropped. We consider problems in which incoming pulses are generated from pulse boundary data as well as problems in which a single outgoing pulse reflects off a possibly curved boundary to produce a number of incoming pulses. Although we focus here on boundary problems, it is clear that similar results hold by similar methods for the Cauchy problem for $N\times N$ systems in free space.

preprint2016arXiv

Geometric optics for surface waves in nonlinear elasticity

This work is devoted to the analysis of high frequency solutions to the equations of nonlinear elasticity in a half-space. We consider surface waves (or more precisely, Rayleigh waves) arising in the general class of isotropic hyperelastic models, which includes in particular the Saint Venant-Kirchhoff system. Work has been done by a number of authors since the 1980s on the formulation and well-posedness of a nonlinear evolution equation whose (exact) solution gives the leading term of an \emph{approximate} Rayleigh wave solution to the underlying elasticity equations. This evolution equation, which we refer to as "the amplitude equation", is an integrodifferential equation of nonlocal Burgers type. We begin by reviewing and providing some extensions of the theory of the amplitude equation. The remainder of the paper is devoted to a rigorous proof in 2D that exact, highly oscillatory, Rayleigh wave solutions $u^\eps$ to the nonlinear elasticity equations exist on a fixed time interval independent of the wavelength $\eps$, and that the approximate Rayleigh wave solution provided by the analysis of the amplitude equation is indeed close in a precise sense to $u^\eps$ on a time interval independent of $\eps$. The paper focuses mainly on the case of Rayleigh waves that are \emph{pulses}, which have profiles with continuous Fourier spectrum, but our method applies equally well to the case of wavetrains, whose Fourier spectrum is discrete.

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

The Mach stem equation and amplification in strongly nonlinear geometric optics

We study highly oscillating solutions to a class of weakly well-posed hyperbolic initial boundary value problems. Weak well-posedness is associated with an amplification phenomenon of oscillating waves on the boundary. In the previous works [CGW14, CW14], we have rigorously justified a weakly nonlinear regime for semilinear problems. In that case, the forcing term on the boundary has amplitude O($ε$^2) and oscillates at a frequency O(1/$ε$). The corresponding exact solution, which has been shown to exist on a time interval that is independent of $ε$ $\in$ (0,1], has amplitude O($ε$). In this paper, we deal with the exact same scaling, namely O($ε$^2) forcing term on the boundary and O($ε$) solution, for quasilinear problems. In analogy with [CGM03], this corresponds to a strongly nonlinear regime, and our main result proves solvability for the corresponding WKB cascade of equations, which yields existence of approximate solutions on a time interval that is independent of $ε$ $\in$ (0,1]. Existence of exact solutions close to approximate ones is a stability issue which, as shown in [CGM03], highly depends on the hyperbolic system and on the boundary conditions; we do not address that question here. This work encompasses previous formal expansions in the case of weakly stable shock waves [MR83] and two-dimensional compressible vortex sheets [AM87]. In particular, we prove well-posedness for the leading amplitude equation (the "Mach stem equation") of [MR83] and generalize its derivation to a large class of hyperbolic boundary value problems and to periodic forcing terms. The latter case is solved under a crucial nonresonant assumption and a small divisor condition.

preprint2013arXiv

Amplification of pulses in nonlinear geometric optics

In this companion paper to our study of amplification of wavetrains, we study weakly stable semilinear hyperbolic boundary value problems with pulse data. Here weak stability means that exponentially growing modes are absent, but the so-called uniform Lopatinskii condition fails at some boundary frequency in the hyperbolic region. As a consequence of this degeneracy there is again an amplification phenomenon: outgoing pulses of amplitude $O(\eps^2)$ and wavelength $\eps$ give rise to reflected pulses of amplitude $O(\eps)$, so the overall solution has amplitude $O(\eps)$. Moreover, the reflecting pulses emanate from a radiating pulse that propagates in the boundary along a characteristic of the Lopatinskii determinant. In the case of N*N systems considered here, a single outgoing pulse produces on reflection a family of incoming pulses traveling at different group velocities. Unlike wavetrains, pulses do not interact to produce resonances that affect the leading order profiles. However, pulse interactions do affect lower order profiles and so these interactions have to be estimated carefully in the error analysis. Whereas the error analysis in the wavetrain case dealt with small divisor problems by approximating periodic profiles by trigonometric polynomials (which amounts to using a high frequency cutoff), in the pulse case we approximate decaying profiles with nonzero moments by profiles with zero moments (a low frequency cutoff). Unlike the wavetrain case, we are now able to obtain a rate of convergence in the limit describing convergence of approximate to exact solutions.

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 geometric optics for reflecting uniformly stable pulses

We provide a justification with rigorous error estimates showing that the leading term in weakly nonlinear geometric optics expansions of highly oscillatory reflecting pulses is close to the uniquely determined exact solution for small wavelengths. Pulses reflecting off fixed noncharacteristic boundaries are considered under the assumption that the underlying boundary problem is uniformly spectrally stable in the sense of Kreiss. There are two respects in which these results make rigorous earlier formal treatments of pulses. First, we give a rigorous construction of leading pulse profiles in problems where pulses traveling with many distinct group velocities are, unavoidably, present; and second, we provide a rigorous error analysis which yields a rate of convergence of approximate to exact solutions as the wavelength approaches zero. Unlike wavetrains, interacting pulses do not produce resonances that affect leading order profiles. However, our error analysis shows the importance of estimating pulse interactions in the construction and estimation of correctors. Our results apply to a general class of systems that includes quasilinear problems like the compressible Euler equations; moreover, the same methods yield a stability result for uniformly stable Euler shocks perturbed by highly oscillatory pulses.

preprint2012arXiv

Semilinear geometric optics with boundary amplification

We study weakly stable semilinear hyperbolic boundary value problems with highly oscillatory data. Here weak stability means that exponentially growing modes are absent, but the so-called uniform Lopatinskii condition fails at some boundary frequency $β$ in the hyperbolic region. As a consequence of this degeneracy there is an amplification phenomenon: outgoing waves of amplitude $O(\eps^2)$ and wavelength $\eps$ give rise to reflected waves of amplitude $O(\eps)$, so the overall solution has amplitude $O(\eps)$. Moreover, the reflecting waves emanate from a radiating wave that propagates in the boundary along a characteristic of the Lopatinskii determinant. An approximate solution that displays the qualitative behavior just described is constructed by solving suitable profile equations that exhibit a loss of derivatives, so we solve the profile equations by a Nash-Moser iteration. The exact solution is constructed by solving an associated singular problem involving singular derivatives of the form $\partial_{x'}+β\frac{\partial_{θ_0}}{\eps}$, $x'$ being the tangential variables with respect to the boundary. Tame estimates for the linearization of that problem are proved using a first-order calculus of singular pseudodifferential operators constructed in the companion article \cite{CGW2}. These estimates exhibit a loss of one singular derivative and force us to construct the exact solution by a separate Nash-Moser iteration. The same estimates are used in the error analysis, which shows that the exact and approximate solutions are close in $L^\infty$ on a fixed time interval independent of the (small) wavelength $\eps$. The approach using singular systems allows us to avoid constructing high order expansions and making small divisor assumptions.

preprint2012arXiv

Singular pseudodifferential calculus for wavetrains and pulses

We develop a singular pseudodifferential calculus. The symbols that we consider do not satisfy the standard decay with respect to the frequency variables. We thus adopt a strategy based on the Calderon-Vaillancourt Theorem. The remainders in the symbolic calculus are bounded operators on $L^2$, whose norm is measured with respect to some small parameter. Our main improvement with respect to an earlier work by Williams consists in showing a regularization effect for the remainders. Due to a nonstandard decay in the frequency variables, the regularization takes place in a scale of anisotropic, and singular, Sobolev spaces. Our analysis allows to extend previous results on the existence of highly oscillatory solutions to nonlinear hyperbolic problems. The results are also used in a companion work to justify nonlinear geometric optics with boundary amplification, which corresponds to a more singular regime than any other one considered before. The analysis is carried out with either an additional real or periodic variable in order to cover problems for pulses or wavetrains in geometric optics.

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

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)$.

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.