Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms

We present a numerical scheme for solving the time-independent nonlinear Gross-Pitaevskii equation in two dimensions describing the Bose-Einstein condensate of trapped interacting neutral atoms at zero temperature. The trap potential is taken to be of the harmonic-oscillator type and the interaction both attractive and repulsive. The Gross-Pitaevskii equation is numerically integrated consistent with the correct boundary conditions at the origin and in the asymptotic region. Rapid convergence is obtained in all cases studied. In the attractive case there is a limit Co the maximum number of atoms in the condensate. (C) 2000 Published by Elsevier B.V. B.V. All rights reserved.

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical approximation of the Ginzburg-Landau equation with memory effects in the dynamics of phase transitions

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Physical approximations for the nonlinear evolution of perturbations in inhomogeneous dark energy scenarios

The abundance and distribution of collapsed objects such as galaxy clusters will become an important tool to investigate the nature of dark energy and dark matter. Number counts of very massive objects are sensitive not only to the equation of state of dark energy, which parametrizes the smooth component of its pressure, but also to the sound speed of dark energy, which determines the amount of pressure in inhomogeneous and collapsed structures. Since the evolution of these structures must be followed well into the nonlinear regime, and a fully relativistic framework for this regime does not exist yet, we compare two approximate schemes: the widely used spherical collapse model and the pseudo-Newtonian approach. We show that both approximation schemes convey identical equations for the density contrast, when the pressure perturbation of dark energy is parametrized in terms of an effective sound speed. We also make a comparison of these approximate approaches to general relativity in the linearized regime, which lends some support to the approximations.

Numerical simulation of Ginzburg-Landau-Langevin equations

This work is concerned with non-equilibrium phenomena, with focus on the numerical simulation of the relaxation of non-conserved order parameters described by stochastic kinetic equations known as Ginzburg-Landau-Langevin (GLL) equations. We propose methods for solving numerically these type of equations, with additive and multiplicative noises. Illustrative applications of the methods are presented for different GLL equations, with emphasis on equations incorporating memory effects.

High-order upwinding and the hydraulic jump

This work is concerned with the computation of incompressible axisymmetric and fall three-dimensional free-surface flows. In particular, the circular-hydraulic jump is simulated and compared with approximate analytic solutions. However, the principal thrust of this paper is to provide a real problem as a test bed for comparing the many existing convective approximations. Their performance is compared; SMART, HLPA and VONOS emerge as acceptable upwinding methods for this problem. Copyright (C) 2002 John Wiley Sons, Ltd.

Some comments on the numerical simulation of self-synchronization of four non-ideal exciters

In this paper, self-synchronization of four non-ideal exciters is examined via numerical simulation. The mathematical model consists of four unbalanced direct Current motors with limited power supply mounted on a flexible Structural frame support. (c) 2004 Elsevier B.V. All rights reserved.

Numerical simulation of a radial diffuser turbulent airflow

In the present work are presented results from numerical simulations performed with the ANSYS-CFX (R) code. We have studied a radial diffuser flow case, which is the main academic problem used to study the flow behavior on flat plate valves. The radial flow inside the diffuser has important behavior such as the turbulence decay downstream and recirculation regions inside the valve flow channel due to boundary layer detachment. These flow structures are present in compressor reed valve configurations, influencing to a greater extent the compressor efficiency. The main target of the present paper was finding the simulation set-up (computational domain, boundary conditions and turbulence model) that better fits with experimental data published by Tabatabai and Pollard. The local flow turbulence and velocity profiles were investigated using four different turbulence models, two different boundary conditions set-up, two different computational domains and three different flow conditions (Re-in - Reynolds number at the diffuser inlet). We used the Reynolds stress (BSL); the k-epsilon; the RNG k-epsilon; and the shear stress transport (SST) k-omega turbulence models. The performed analysis and comparison of the computational results with experimental data show that the choice of the turbulence model, as well as the choice of the other computational conditions, plays an important role in the results physical quality and accuracy. (c) 2007 Elsevier B.V. All rights reserved.

Comments on nonlinear dynamics of a non-ideal Duffing-Rayleigh oscillator: Numerical and analytical approaches

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical exploration of resonant dynamics in the system of Saturnian major satellites

We numerically investigate the long-term dynamics of the Saturnian system by analyzing the Fourier spectra of ensembles of orbits taken around the current orbits of Mimas, Enceladus, Tethys, Rhea and Hyperion. We construct dynamical maps around the current position of these satellites in their respective phase spaces. The maps are the result of a great deal of numerical simulations where we adopt dense sets of initial conditions and different satellite configurations. Several structures associated to the current two-body mean-motion resonances, unstable regions associated to close approaches between the satellites, and three-body mean-motion resonances in the system, are identified in the map. (C) 2010 Elsevier Ltd. All rights reserved.

Rational approximations of the Arrhenius integral using Jacobi fractions and gaussian quadrature

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

An introduction to numerical methods in low-dimensional quantum systems

This is an introductory course to the Lanczos Method and Density Matrix Renormalization Group Algorithms (DMRG), two among the leading numerical techniques applied in studies of low-dimensional quantum models. The idea of studying the models on clusters of a finite size in order to extract their physical properties is briefly discussed. The important role played by the model symmetries is also examined. Special emphasis is given to the DMRG.

Numerical chromosome polymorphism in Mikania cordifolia Willd. (Asteraceae)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Numerical study about natural escape and capture routes by the Moon via Lagrangian points L1 and L2

The lunar sphere of influence, whose radius is some 66,300 km, has regions of stable orbits around the Moon and also regions that contain trajectories which, after spending some time around the Moon, escape and are later recaptured by lunar gravity. Both the escape and the capture occur along the Lagrangian equilibrium points L1 and L2. In this study, we mapped out the region of lunar influence considering the restricted three-body Earth-Moon-particle problem and the four-body Sun-Earth-Moon-particle (probe) problem. We identified the stable trajectories, and the escape and capture trajectories through the L I and L2 in plots of the eccentricity versus the semi-major axis as a function of the time that the energy of the osculating lunar trajectory in the two-body Moon-particle problem remains negative. We also investigated the properties of these routes, giving special attention to the fact that they supply a natural mechanism for performing low-energy transfers between the Earth and the Moon, and can thus be useful on a great number of future missions. (C) 2007 Published by Elsevier Ltd on behalf of COSPAR.