6 research outputs found
Gradient Flow Finite Element Discretizations with Energy-Based Adaptivity for the Gross-Pitaevskii Equation
We present an effective adaptive procedure for the numerical approximation of
the steady-state Gross-Pitaevskii equation. Our approach is solely based on
energy minimization, and consists of a combination of gradient flow iterations
and adaptive finite element mesh refinements. Numerical tests show that this
strategy is able to provide highly accurate results, with optimal convergence
rates with respect to the number of freedom
Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency
We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii
eigenvalue problem based on an energy inner product that depends on time
through the density of the flow itself. The gradient flow is well-defined and
converges to an eigenfunction. For ground states we can quantify the
convergence speed as exponentially fast where the rate depends on spectral gaps
of a linearized operator. The forward Euler time discretization of the flow
yields a numerical method which generalizes the inverse iteration for the
nonlinear eigenvalue problem. For sufficiently small time steps, the method
reduces the energy in every step and converges globally in to an
eigenfunction. In particular, for any nonnegative starting value, the ground
state is obtained. A series of numerical experiments demonstrates the
computational efficiency of the method and its competitiveness with established
discretizations arising from other gradient flows for this problem
An inverse iteration method for eigenvalue problems with eigenvector nonlinearities
© 2014 Society for Industrial and Applied Mathematics. Consider a symmetric matrix A(v) â ânxn depending on a vector v â ân and satisfying the property A(αv) = A(v) for any α â â\{0}. We will here study the problem of finding (λ, v) â â x ân\{0} such that (â, v) is an eigenpair of the matrix A(v) and we propose a generalization of inverse iteration for eigenvalue problems with this type of eigenvector nonlinearity. The convergence of the proposed method is studied and several convergence properties are shown to be analogous to inverse iteration for standard eigenvalue problems, including local convergence properties. The algorithm is also shown to be equivalent to a particular discretization of an associated ordinary differential equation, if the shift is chosen in a particular way. The algorithm is adapted to a variant of the Schrödinger equation known as the Gross-Pitaevskii equation. We use numerical simulations to illustrate the convergence properties, as well as the efficiency of the algorithm and the adaption.status: publishe