A fourth order accurate finite difference scheme for the computation of elastic waves

A finite difference for elastic waves is introduced. The model is based on the first order system of equations for the velocities and stresses. The differencing is fourth order accurate on the spatial derivatives and second order accurate in time. The model is tested on a series of examples including the Lamb problem, scattering from plane interf aces and scattering from a fluid-elastic interface. The scheme is shown to be effective for these problems. The accuracy and stability is insensitive to the Poisson ratio. For the class of problems considered here it is found that the fourth order scheme requires for two-thirds to one-half the resolution of a typical second order scheme to give comparable accuracy

Statistics and Characteristics of Spatio-Temporally Rare Intense Events in Complex Ginzburg-Landau Models

We study the statistics and characteristics of rare intense events in two types of two dimensional Complex Ginzburg-Landau (CGL) equation based models. Our numerical simulations show finite amplitude collapse-like solutions which approach the infinite amplitude solutions of the nonlinear Schr\"{o}dinger (NLS) equation in an appropriate parameter regime. We also determine the probability distribution function (PDF) of the amplitude of the CGL solutions, which is found to be approximately described by a stretched exponential distribution, $P(|A|) \approx e^{-|A|^\eta}$, where $\eta < 1$. This non-Gaussian PDF is explained by the nonlinear characteristics of individual bursts combined with the statistics of bursts. Our results suggest a general picture in which an incoherent background of weakly interacting waves, occasionally, `by chance', initiates intense, coherent, self-reinforcing, highly nonlinear events.Comment: 7 pages, 9 figure

A system of ODEs for a Perturbation of a Minimal Mass Soliton

We study soliton solutions to a nonlinear Schrodinger equation with a saturated nonlinearity. Such nonlinearities are known to possess minimal mass soliton solutions. We consider a small perturbation of a minimal mass soliton, and identify a system of ODEs similar to those from Comech and Pelinovsky (2003), which model the behavior of the perturbation for short times. We then provide numerical evidence that under this system of ODEs there are two possible dynamical outcomes, which is in accord with the conclusions of Pelinovsky, Afanasjev, and Kivshar (1996). For initial data which supports a soliton structure, a generic initial perturbation oscillates around the stable family of solitons. For initial data which is expected to disperse, the finite dimensional dynamics follow the unstable portion of the soliton curve.Comment: Minor edit

Renormalizing Partial Differential Equations

In this review paper, we explain how to apply Renormalization Group ideas to the analysis of the long-time asymptotics of solutions of partial differential equations. We illustrate the method on several examples of nonlinear parabolic equations. We discuss many applications, including the stability of profiles and fronts in the Ginzburg-Landau equation, anomalous scaling laws in reaction-diffusion equations, and the shape of a solution near a blow-up point.Comment: 34 pages, Latex; [email protected]; [email protected]