Elliptic solutions and solitary waves of a higher order KdV--BBM long wave equation

We provide conditions for existence of hyperbolic, unbounded periodic and elliptic solutions in terms of Weierstrass $\wp$ functions of both third and fifth-order KdV--BBM (Korteweg-de Vries--Benjamin, Bona \& Mahony) regularized long wave equation. An analysis for the initial value problem is developed together with a local and global well-posedness theory for the third-order KdV--BBM equation. Traveling wave reduction is used together with zero boundary conditions to yield solitons and periodic unbounded solutions, while for nonzero boundary conditions we find solutions in terms of Weierstrass elliptic $\wp$ functions. For the fifth-order KdV--BBM equation we show that a parameter $\gamma=\frac {1}{12}$, for which the equation has a Hamiltonian, represents a restriction for which there are constraint curves that never intersect a region of unbounded solitary waves, which in turn shows that only dark or bright solitons and no unbounded solutions exist. Motivated by the lack of a Hamiltonian structure for $\gamma\neq\frac{1}{12}$ we develop $H^k$ bounds, and we show for the non Hamiltonian system that dark and bright solitons coexist together with unbounded periodic solutions. For nonzero boundary conditions, due to the complexity of the nonlinear algebraic system of coefficients of the elliptic equation we construct Weierstrass solutions for a particular set of parameters only.Comment: 13 pages, 6 figure

A new dissipation term for finite-difference simulations in Relativity

We present a new numerical dissipation algorithm, which can be efficiently used in combination with centered finite-difference methods. We start from a formulation of centered finite-volume methods for Numerical Relativity, in which third-order space accuracy can be obtained by employing just piecewise-linear reconstruction. We obtain a simplified version of the algorithm, which can be viewed as a centered finite-difference method plus some 'adaptive dissipation'. The performance of this algorithm is confirmed by numerical results obtained from 3D black hole simulations.Comment: Talk presented at the Spanish Relativity Meeting (Tenerife 2007

Constraint-preserving boundary conditions in the 3+1 first-order approach

A set of energy-momentum constraint-preserving boundary conditions is proposed for the first-order Z4 case. The stability of a simple numerical implementation is tested in the linear regime (robust stability test), both with the standard corner and vertex treatment and with a modified finite-differences stencil for boundary points which avoids corners and vertices even in cartesian-like grids. Moreover, the proposed boundary conditions are tested in a strong field scenario, the Gowdy waves metric, showing the expected rate of convergence. The accumulated amount of energy-momentum constraint violations is similar or even smaller than the one generated by either periodic or reflection conditions, which are exact in the Gowdy waves case. As a side theoretical result, a new symmetrizer is explicitly given, which extends the parametric domain of symmetric hyperbolicity for the Z4 formalism. The application of these results to first-order BSSN-like formalisms is also considered.Comment: Revised version, with conclusive numerical evidence. 23 pages, 12 figure

The limit of a Stanley-Wilf sequence is not always rational, and layered patterns beat monotone patterns

We show the first known example for a pattern $q$ for which $\lim_{n\to \infty} \sqrt[n]{S_n(q)}$ is not an integer. We find the exact value of the limit and show that it is irrational. Then we generalize our results to an infinite sequence of patterns. Finally, we provide further generalizations that start explaining why certain patterns are easier to avoid than others. Finally, we show that if $q$ is a layered pattern of length $k$, then $L(q)\geq (k-1)^2$ holds.Comment: 10 pages, 1 figur