3,635 research outputs found
Solitary waves in the Nonlinear Dirac Equation
In the present work, we consider the existence, stability, and dynamics of
solitary waves in the nonlinear Dirac equation. We start by introducing the
Soler model of self-interacting spinors, and discuss its localized waveforms in
one, two, and three spatial dimensions and the equations they satisfy. We
present the associated explicit solutions in one dimension and numerically
obtain their analogues in higher dimensions. The stability is subsequently
discussed from a theoretical perspective and then complemented with numerical
computations. Finally, the dynamics of the solutions is explored and compared
to its non-relativistic analogue, which is the nonlinear Schr{\"o}dinger
equation. A few special topics are also explored, including the discrete
variant of the nonlinear Dirac equation and its solitary wave properties, as
well as the PT-symmetric variant of the model
Uncertainty quantification for kinetic models in socio-economic and life sciences
Kinetic equations play a major rule in modeling large systems of interacting
particles. Recently the legacy of classical kinetic theory found novel
applications in socio-economic and life sciences, where processes characterized
by large groups of agents exhibit spontaneous emergence of social structures.
Well-known examples are the formation of clusters in opinion dynamics, the
appearance of inequalities in wealth distributions, flocking and milling
behaviors in swarming models, synchronization phenomena in biological systems
and lane formation in pedestrian traffic. The construction of kinetic models
describing the above processes, however, has to face the difficulty of the lack
of fundamental principles since physical forces are replaced by empirical
social forces. These empirical forces are typically constructed with the aim to
reproduce qualitatively the observed system behaviors, like the emergence of
social structures, and are at best known in terms of statistical information of
the modeling parameters. For this reason the presence of random inputs
characterizing the parameters uncertainty should be considered as an essential
feature in the modeling process. In this survey we introduce several examples
of such kinetic models, that are mathematically described by nonlinear Vlasov
and Fokker--Planck equations, and present different numerical approaches for
uncertainty quantification which preserve the main features of the kinetic
solution.Comment: To appear in "Uncertainty Quantification for Hyperbolic and Kinetic
Equations
The low dimensional dynamical system approach in General Relativity: an example
In this paper we explore one of the most important features of the Galerkin
method, which is to achieve high accuracy with a relatively modest
computational effort, in the dynamics of Robinson-Trautman spacetimes.Comment: 7 pages, 5 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
A time-splitting pseudospectral method for the solution of the Gross-Pitaevskii equations using spherical harmonics with generalised-Laguerre basis functions
We present a method for numerically solving a Gross-Pitaevskii system of equations with a harmonic and a toroidal external potential that governs the dynamics of one- and two-component Bose-Einstein condensates. The method we develop maintains spectral accuracy by employing Fourier or spherical harmonics in the angular coordinates combined with generalised-Laguerre basis functions in the radial direction. Using an error analysis, we show that the method presented leads to more accurate results than one based on a sine transform in the radial direction when combined with a time-splitting method for integrating the equations forward in time. In contrast to a number of previous studies, no assumptions of radial or cylindrical symmetry is assumed allowing the method to be applied to 2D and 3D time-dependent simulations. This is accomplished by developing an efficient algorithm that accurately performs the generalised-Laguerre transforms of rotating Bose-Einstein condensates for different orders of the Laguerre polynomials. Using this spatial discretisation together with a second order Strang time-splitting method, we illustrate the scheme on a number of 2D and 3D computations of the ground state of a non-rotating and rotating condensate. Comparisons between previously derived theoretical results for these ground state solutions and our numerical computations show excellent agreement for these benchmark problems. The method is further applied to simulate a number of time-dependent problems including the Kelvin-Helmholtz instability in a two-component rotating condensate and the motion of quantised vortices in a 3D condensate
Structure preserving schemes for the continuum Kuramoto model: phase transitions
The construction of numerical schemes for the Kuramoto model is challenging
due to the structural properties of the system which are essential in order to
capture the correct physical behavior, like the description of stationary
states and phase transitions. Additional difficulties are represented by the
high dimensionality of the problem in presence of multiple frequencies. In this
paper, we develop numerical methods which are capable to preserve these
structural properties of the Kuramoto equation in the presence of diffusion and
to solve efficiently the multiple frequencies case. The novel schemes are then
used to numerically investigate the phase transitions in the case of identical
and non identical oscillators
- …