60 research outputs found
Discrete embedded solitons
We address the existence and properties of discrete embedded solitons (ESs),
i.e., localized waves existing inside the phonon band in a nonlinear
dynamical-lattice model. The model describes a one-dimensional array of optical
waveguides with both the quadratic (second-harmonic generation) and cubic
nonlinearities. A rich family of ESs was previously known in the continuum
limit of the model. First, a simple motivating problem is considered, in which
the cubic nonlinearity acts in a single waveguide. An explicit solution is
constructed asymptotically in the large-wavenumber limit. The general problem
is then shown to be equivalent to the existence of a homoclinic orbit in a
four-dimensional reversible map. From properties of such maps, it is shown that
(unlike ordinary gap solitons), discrete ESs have the same codimension as their
continuum counterparts. A specific numerical method is developed to compute
homoclinic solutions of the map, that are symmetric under a specific reversing
transformation. Existence is then studied in the full parameter space of the
problem. Numerical results agree with the asymptotic results in the appropriate
limit and suggest that the discrete ESs may be semi-stable as in the continuous
case.Comment: A revtex4 text file and 51 eps figure files. To appear in
Nonlinearit
Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow
Steady states and traveling waves play a fundamental role in understanding
hydrodynamic problems. Even when unstable, these states provide the
bifurcation-theoretic explanation for the origin of the observed states. In
turbulent wall-bounded shear flows, these states have been hypothesized to be
saddle points organizing the trajectories within a chaotic attractor. These
states must be computed with Newton's method or one of its generalizations,
since time-integration cannot converge to unstable equilibria. The bottleneck
is the solution of linear systems involving the Jacobian of the Navier-Stokes
or Boussinesq equations. Originally such computations were carried out by
constructing and directly inverting the Jacobian, but this is unfeasible for
the matrices arising from three-dimensional hydrodynamic configurations in
large domains. A popular method is to seek states that are invariant under
numerical time integration. Surprisingly, equilibria may also be found by
seeking flows that are invariant under a single very large Backwards-Euler
Forwards-Euler timestep. We show that this method, called Stokes
preconditioning, is 10 to 50 times faster at computing steady states in plane
Couette flow and traveling waves in pipe flow. Moreover, it can be carried out
using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes
to these popular spectral codes. We explain the convergence rate as a function
of the integration period and Reynolds number by computing the full spectra of
the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid
Dynamics, ed. Alexander Gelfgat (Springer, 2018
Algebraic Methods for Determining Hamiltonian Hopf Bifurcations in Three-Degree-of-Freedom Systems
Nonlinear tapping dynamics of multi-walled carbon nanotube tipped atomic force microcantilevers
- …