2,202 research outputs found
Compact high order finite difference schemes for linear Schrödinger problems on non-uniform meshes
In the present paper a general technique is developed for construction
of compact high-order finite difference schemes to approximate Schrödinger
problems on nonuniform meshes. Conservation of the finite difference schemes
is investigated. Discrete transparent boundary conditions are constructed for
the given high-order finite difference scheme. The same technique is applied
to construct compact high-order approximations of the Robin and Szeftel type
boundary conditions. Results of computational experiments are presente
Reduced-Order Modeling based on Approximated Lax Pairs
A reduced-order model algorithm, based on approximations of Lax pairs, is
proposed to solve nonlinear evolution partial differential equations. Contrary
to other reduced-order methods, like Proper Orthogonal Decomposition, the space
where the solution is searched for evolves according to a dynamics specific to
the problem. It is therefore well-suited to solving problems with progressive
waves or front propagation. Numerical examples are shown for the KdV and FKPP
(nonlinear reaction diffusion) equations, in one and two dimensions
Approximated Lax Pairs for the Reduced Order Integration of Nonlinear Evolution Equations
A reduced-order model algorithm, called ALP, is proposed to solve nonlinear
evolution partial differential equations. It is based on approximations of
generalized Lax pairs. Contrary to other reduced-order methods, like Proper
Orthogonal Decomposition, the basis on which the solution is searched for
evolves in time according to a dynamics specific to the problem. It is
therefore well-suited to solving problems with progressive front or wave
propagation. Another difference with other reduced-order methods is that it is
not based on an off-line / on-line strategy. Numerical examples are shown for
the linear advection, KdV and FKPP equations, in one and two dimensions
Recommended from our members
Compact high order finite difference schemes for linear Schrödinger problems on non-uniform meshes
In the present paper a general technique is developed for construction
of compact high-order finite difference schemes to approximate Schrödinger
problems on nonuniform meshes. Conservation of the finite difference schemes
is investigated. Discrete transparent boundary conditions are constructed for
the given high-order finite difference scheme. The same technique is applied
to construct compact high-order approximations of the Robin and Szeftel type
boundary conditions. Results of computational experiments are presente
High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension
The nonlinear Helmholtz equation (NLH) models the propagation of
electromagnetic waves in Kerr media, and describes a range of important
phenomena in nonlinear optics and in other areas. In our previous work, we
developed a fourth order method for its numerical solution that involved an
iterative solver based on freezing the nonlinearity. The method enabled a
direct simulation of nonlinear self-focusing in the nonparaxial regime, and a
quantitative prediction of backscattering. However, our simulations showed that
there is a threshold value for the magnitude of the nonlinearity, above which
the iterations diverge. In this study, we numerically solve the one-dimensional
NLH using a Newton-type nonlinear solver. Because the Kerr nonlinearity
contains absolute values of the field, the NLH has to be recast as a system of
two real equations in order to apply Newton's method. Our numerical simulations
show that Newton's method converges rapidly and, in contradistinction with the
iterations based on freezing the nonlinearity, enables computations for very
high levels of nonlinearity. In addition, we introduce a novel compact
finite-volume fourth order discretization for the NLH with material
discontinuities.The one-dimensional results of the current paper create a
foundation for the analysis of multi-dimensional problems in the future.Comment: 47 pages, 8 figure
On Form Factors in nested Bethe Ansatz systems
We investigate form factors of local operators in the multi-component Quantum
Non-linear Schr\"odinger model, a prototype theory solvable by the so-called
nested Bethe Ansatz. We determine the analytic properties of the infinite
volume form factors using the coordinate Bethe Ansatz solution and we establish
a connection with the finite volume matrix elements. In the two-component
models we derive a set of recursion relations for the "magnonic form factors",
which are the matrix elements on the nested Bethe Ansatz states. In certain
simple cases (involving states with only one spin-impurity) we obtain explicit
solutions for the recursion relations.Comment: 34 pages, v2 (minor modifications
- …