Source author record

Gregory Lyng

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

ResearcherUnclaimed source record

nlin.SI math.AP math.CA math.NA

Catalog footprint

What is connected

4works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2014arXiv

Analysis and development of compact finite difference schemes with optimized numerical dispersion relations

Finite difference approximation, in addition to Taylor truncation errors, introduces numerical dispersion-and-dissipation errors into numerical solutions of partial differential equations. We analyze a class of finite difference schemes which are designed to minimize these errors (at the expense of formal order of accuracy), and we analyze the interplay between the Taylor truncation errors and the dispersion-and-dissipation errors during mesh refinement. In particular, we study the numerical dispersion relation of the fully discretized non-dispersive transport equation in one and two space dimensions. We derive the numerical phase error and the $L^2$-norm error of the solution in terms of the dispersion-and-dissipation error. Based on our analysis, we investigate the error dynamics among various optimized compact schemes and the unoptimized higher-order generalized Padé compact schemes. The dynamics shed light on the principles of designing suitable optimized compact schemes for a given problem. Using these principles as guidelines, we then propose an optimized scheme that prescribes the numerical dispersion relation before finding the corresponding discretization. This approach produces smaller numerical dispersion-and-dissipation errors for linear and nonlinear problems, compared with the unoptimized higher-order compact schemes and other optimized schemes developed in the literature. Finally, we discuss the difficulty of developing an optimized composite boundary scheme for problems with non-trivial boundary conditions. We propose a composite scheme that introduces a buffer zone to connect an optimized interior scheme and an unoptimized boundary scheme. Our numerical experiments show that this strategy produces small $L^2$-norm error when a wave packet passes through the non-periodic boundary.

preprint2012arXiv

The Gaussian semiclassical soliton ensemble and numerical methods for the focusing nonlinear Schrödinger equation

We report on a number of careful numerical experiments motivated by the semiclassical (zero-dispersion, ε\downarrow 0) limit of the focusing nonlinear Schrödinger equation. Our experiments are designed to study the evolution of a particular family of perturbations of the initial data. These asymptotically small perturbations are precisely those that result from modifying the initial-data by using formal approximations to the spectrum of the associated spectral problem; such modified data has always been a standard part of the analysis of zero-dispersion limits of integrable systems. However, in the context of the focusing nonlinear Schrödinger equation, the ellipticity of the Whitham equations casts some doubt on the validity of this procedure. To carry out our experiments, we introduce an implicit finite difference scheme for the partial differential equation, and we validate both the proposed scheme and the standard split-step scheme against a numerical implementation of the inverse scattering transform for a special case in which the scattering data is known exactly. As part of this validation, we also investigate the use of the Krasny filter which is sometimes suggested as appropriate for nearly ill-posed problems such as we consider here. Our experiments show that that the O(ε) rate of convergence of the modified data to the true data is propagated to positive times including times after wave breaking.

preprint2005arXiv

The N-soliton of the focusing nonlinear Schroedinger equation for N large

We present a detailed analysis of the solution of the focusing nonlinear Schrödinger equation with initial condition $ψ(x,0)=N {\rm sech}(x)$ in the limit $N\to\infty$. We begin by presenting new and more accurate numerical reconstructions of the $N$-soliton by inverse scattering (numerical linear algebra) for N=5, 10, 20, and 40. We then recast the inverse-scattering problem as a Riemann-Hilbert problem and provide a rigorous asymptotic analysis of this problem in the large-N limit. For those $(x,t)$ where results have been obtained by other authors, we improve the error estimates from $O(N^{-1/3})$ to $O(N^{-1})$. We also analyze the Fourier power spectrum in this regime and relate the results to the optical phenomenon of supercontinuum generation. We then study the $N$-soliton for values of $(x,t)$ where analysis has not been carried out before, and we uncover new phenomena. The main discovery of this paper is the mathematical mechanism for a secondary caustic (phase transition), which turns out to differ from the mechanism that generates the primary caustic. The mechanism for the generation of the secondary caustic depends essentially on the discrete nature of the spectrum for the $N$-soliton, and more significantly, cannot be recovered from an analysis of an ostensibly similar Riemann-Hilbert problem in the conditions of which a certain formal continuum limit is taken on an {\em ad hoc} basis

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.