Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

In this paper we present a finite difference method for solving two-dimensional viscoelastic unsteady free surface flows governed by the single equation version of the eXtended Pom-Pom (XPP) model. The momentum equations are solved by a projection method which uncouples the velocity and pressure fields. We are interested in low Reynolds number flows and, to enhance the stability of the numerical method, an implicit technique for computing the pressure condition on the free surface is employed. This strategy is invoked to solve the governing equations within a Marker-and-Cell type approach while simultaneously calculating the correct normal stress condition on the free surface. The numerical code is validated by performing mesh refinement on a two-dimensional channel flow. Numerical results include an investigation of the influence of the parameters of the XPP equation on the extrudate swelling ratio and the simulation of the Barus effect for XPP fluids. (C) 2010 Elsevier B.V. All rights reserved.

A Numerical method for the semilinear stochastic transport equation

Trabalho apresentado no XXXV CNMAC, Natal-RN, 2014.

Stability of multistep second derivative methods for integro-differential equations

The stability of multistep second derivative methods for integro-differential equations is examined through a test equation which allows for the construction of the associated characteristic polynomial and its region of stability (roots in the unit circle) at a proper parameter space. (c) 2004 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Numerical solution of the eXtended Pom-Pom model for viscoelastic free surface flows

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

Theoretical study of the structure and stability of NbxOy, and NbxOy+ (x = 1-3; y = 2-5, 7, 8) clusters

Geometric, thermodynamic and electronic properties of cluster neutrals NbxOy and cations NbxOy+ (x = 1-3; y = 2-5, 7, 8) have been characterized theoretically. A DFT calculation using a hybrid combination of B3LYP with contracted Huzinaga basis sets. Numerical results of the relative stabilities, ionization potentials and band gaps of different clusters are in agreement with experiment. Analysis of dissociation channels supports the more stable building blocks as formed by NbO2, NbO2+ NbO3 and NbO3+ stoichiometries. The net atomic charges suggest that oxygen donor molecules can interact more favorably on central niobium atoms of cluster cations, while the interaction with oxygen acceptor molecules is more favorable on the terminal oxygen atoms of neutral clusters. A topological analysis of the electron localization function gradient field indicates that the clusters may be described as having a strong ionic interaction between Nb and O atoms. Published by Elsevier B.V. B.V.

Attitude stability of artificial satellites subject to gravity gradient torque

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

On the stability of the satellites of asteroid 87 Sylvia

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

Study of Stability of Rotational Motion of Spacecraft with Canonical Variables

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

Superfluid Fermi-Fermi mixture: Phase diagram, stability, and soliton formation

We study the phase diagram for a dilute Bardeen-Cooper-Schrieffer superfluid Fermi-Fermi mixture (of distinct mass) at zero temperature using energy densities for the superfluid fermions in one (1D), two (2D), and three (3D) dimensions. We also derive the dynamical time-dependent nonlinear Euler-Lagrange equation satisfied by the mixture in one dimension using this energy density. We obtain the linear stability conditions for the mixture in terms of fermion densities of the components and the interspecies Fermi-Fermi interaction. In equilibrium there are two possibilities. The first is that of a uniform mixture of the two components, the second is that of two pure phases of two components without any overlap between them. In addition, a mixed and a pure phase, impossible in 1D and 2D, can be created in 3D. We also obtain the conditions under which the uniform mixture is stable from an energetic consideration. The same conditions are obtained from a modulational instability analysis of the dynamical equations in 1D. Finally, the 1D dynamical equations for the system are solved numerically and by variational approximation (VA) to study the bright solitons of the system for attractive interspecies Fermi-Fermi interaction in 1D. The VA is found to yield good agreement to the numerical result for the density profile and chemical potential of the bright solitons. The bright solitons are demonstrated to be dynamically stable. The experimental realization of these Fermi-Fermi bright solitons seems possible with present setups.

Synchronization: Stability and duration time

We consider the problem of stability and duration of the synchronization process between self-excited oscillators, both in their regular and chaotic states. Making use of the properties of Hill equation describing the deviation between the slave and the master, we derive the stability conditions and expressions of the synchronization time. A fairly good agreement is obtained between the analytical and numerical results.

Stability analysis of the D-dimensional nonlinear Schrodinger equation with trap and two- and three-body interactions

Considering the static solutions of the D-dimensional nonlinear Schrodinger equation with trap and attractive two-body interactions, the existence of stable solutions is limited to a maximum critical number of particles, when D greater than or equal to 2. In case D = 2, we compare the variational approach with the exact numerical calculations. We show that, the addition of a positive three-body interaction allows stable solutions beyond the critical number. In this case, we also introduce a dynamical analysis of the conditions for the collapse. (C) 2000 Published by Elsevier B.V. B.V. All rights reserved.

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)

On the stability of hypothetical satellites coorbital to Mimas or Enceladus

In the present work, we study the stability of hypothetical satellites that are coorbital with Enceladus and Mimas. We performed numerical simulations of 50 particles around the triangular Lagrangian equilibrium points of Enceladus and Mimas taking into account the perturbation of Mimas, Enceladus, Tethys, Dione, Titan and the oblateness of Saturn. All particles remain on tadpole orbits after 10 000 yr of integration. Since in the past the orbit of Enceladus and Mimas expanded due to the tidal perturbation, we also simulated the system with Enceladus and Mimas at several different values of semimajor axes. The results show that in general the particles remain on tadpole orbits. The exceptions occur when Enceladus is at semimajor axes that correspond to 6:7, 5:6 and 4:5 resonances with Mimas. Therefore, if Enceladus and Mimas had satellites librating around their Lagrangian triangular points in the past, they would have been removed if Enceladus crossed one of these first-order resonances with Mimas.

ITERATIVE NUMERICAL-SOLUTION OF SCATTERING PROBLEMS

An iterative Neumann series method, employing a real auxiliary scattering integral equation, is used to calculate scattering lengths and phase shifts for the atomic Yukawa and exponential potentials. For these potentials the original Neumann series diverges. The present iterative method yields results that are far better, in convergence, stability and precision, than other momentum space methods. Accurate result is obtained in both cases with an estimated error of about 1 in 10(10) after some 8-10 iterations.

Numerical study of low-cost alternative orbits around the Moon

In this paper, we have investigated a region of direct stable orbits around the Moon, whose stability is related to the H2 Family of periodic orbits and to the quasi-periodic orbits that oscillate around them. The stability criteria adopted was that the path did not escape from the Moon during an integration period of 1000 days (remaining with negative two-body Moon-probe orbital energy during this period). Considering the three-dimensional four-body Sun-Earth-Moon-probe problem, we investigated the evolution of the size of the stability region, taking into account the eccentricity of the Earth's orbit, the eccentricity and inclination of the Moon's orbit, and the solar radiation pressure on the probe. We also investigated the evolution of the region's size and its location by varying the inclination of the probe's initial osculating orbit relative to the Moon's orbital plane between 0 degrees and 180 degrees. The size of the stability region diminishes; nevertheless, it remains significant for 0 <= i <= 25 degrees and 35 degrees <= i <= 45 degrees. The orbits of this region could be useful for missions by space vehicles that must remain in orbit around the Moon for periods of up to 1000 days, requiring low maintenance costs. (c) 2005 COSPAR. Published by Elsevier Ltd. All rights reserved.