Stability regions around the components of the triple system 2001 SN263

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The stability evolution of a family of simply periodic lunar orbits

The present work deals with a family of simply periodic orbits around the Moon in the rotating Earth Moon-particle system. Taking the framework of the planar, circular, restricted three-body problem, we follow the evolution of this family of periodic orbits using the numerical technique of Poincaré surface of section. The maximum amplitude of oscillation about the periodic orbits are determined and can be used as a parameter to measure the degree of stability in the phase space for such orbits. Despite the fact that the whole family of periodic orbits remain stable, there is a dichotomy in the quasi-periodic ones at the Jacobi constant Cj = 2.85. The quasi-periodic orbits with Cj < 2.85 oscillate around the periodic orbits in a different way from those with Cj > 2.85. © 1999 Elsevier Science Ltd. All rights reserved.

Stability of trapped Bose-Einstein condensates

In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a Gaussian variational approach and exact numerical simulations. We also discuss the validity of the criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions for the existence of two phases in the condensate. In this case, the magnitude of the oscillations between the two phases are studied considering sufficient large initial chirps. The occurrence of collapse in a BEC with repulsive two-body interaction is also shown to be possible.

Improving the stability conditions of TS fuzzy systems with fuzzy Lyapunov functions

In this paper, the fuzzy Lyapunov function approach is considered for stabilizing continuous-time Takagi-Sugeno fuzzy systems. Previous linear matrix inequality (LMI) stability conditions are relaxed by exploring further the properties of the time derivatives of premise membership functions and by introducing a slack LMI variable into the problem formulation. The stability results are thus used in the state feedback design which is also solved in terms of LMIs. Numerical examples illustrate the efficiency of the new stabilizing conditions presented. © 2011 IFAC.

Stochastic exponential stability of singular linear systems with Markov jump parameters

This paper deals with exponential stability of discrete-time singular systems with Markov jump parameters. We propose a set of coupled generalized Lyapunov equations (CGLE) that provides sufficient conditions to check this property for this class of systems. A method for solving the obtained CGLE is also presented, based on iterations of standard singular Lyapunov equations. We present also a numerical example to illustrate the effectiveness of the approach we are proposing.

Three-dimensional transient complex free surface flows: Numerical simulation of XPP fluid

In this paper we present a finite difference MAC-type approach for solving three-dimensional viscoelastic incompressible free surface flows governed by the eXtended Pom-Pom (XPP) model, considering a wide range of parameters. The numerical formulation presented in this work is an extension to three-dimensions of our implicit technique [Journal of Non-Newtonian Fluid Mechanics 166 (2011) 165-179] for solving two-dimensional viscoelastic free surface flows. To enhance the stability of the numerical method, we employ a combination of the projection method with an implicit technique for treating the pressure on the free surfaces. The differential constitutive equation of the fluid is solved using a second-order Runge-Kutta scheme. The numerical technique is validated by performing a mesh refinement study on a pipe flow, and the numerical results presented include the simulation of two complex viscoelastic free surface flows: extrudate-swell problem and jet buckling phenomenon. © 2013 Elsevier B.V.

Stability of trapped degenerate dipolar Bose and Fermi gases

Trapped degenerate dipolar Bose and Fermi gases of the cylindrical symmetry with the polarization vector along the symmetry axis are only stable for the strength of dipolar interaction below a critical value. In the case of bosons, the stability of such a dipolar Bose-Einstein condensate (BEC) is investigated for different strengths of contact and dipolar interactions using a variational approximation and a numerical solution of a mean-field model. In the disc shape, with the polarization vector perpendicular to the plane of the disc, the atoms experience an overall dipolar repulsion and this fact should contribute to the stability. However, a complete numerical solution of the dynamics leads to the collapse of a strongly disc-shaped dipolar BEC due to the long-range anisotropic dipolar interaction. In the case of fermions, the stability of a trapped single-component degenerate dipolar Fermi gas is studied including the Hartree-Fock exchange and Brueckner-Goldstone correlation energies in the local-density approximation valid for a large number of atoms. Estimates for the maximum allowed number of polar Bose and Fermi molecules in the BEC and degenerate Fermi gas are given. © 2013 IOP Publishing Ltd.

Semi-analytical study of the rotational motion stability of artificial satellites using quaternions

This study at aims performing the stability analysis of the rotational motion to artificial satellites using quaternions to describe the satellite attitude (orientation on the space). In the system of rotational motion equations, which is composed by four kinematic equations of the quaternions and by the three Euler equations in terms of the rotational spin components. The influence of the gravity gradient and the direct solar radiation pressure torques have been considered. Equilibrium points were obtained through numerical simulations using the softwares Matlab and Octave, which are then analyzed by the Routh-Hurwitz Stability Criterion.

3D stability orbits close to 433 Eros using an effective polyhedral model method

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

Stability analysis of Crank-Nicolson and Euler schemes for time-dependent diffusion equations

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

MHD numerical simulations of colliding winds in massive binary systems - I. Thermal versus non-thermal radio emission

In the past few decades detailed observations of radio and X-ray emission from massive binary systems revealed a whole new physics present in such systems. Both thermal and non-thermal components of this emission indicate that most of the radiation at these bands originates in shocks. O and B-type stars and WolfRayet (WR) stars present supersonic and massive winds that, when colliding, emit largely due to the freefree radiation. The non-thermal radio and X-ray emissions are due to synchrotron and inverse Compton processes, respectively. In this case, magnetic fields are expected to play an important role in the emission distribution. In the past few years the modelling of the freefree and synchrotron emissions from massive binary systems have been based on purely hydrodynamical simulations, and ad hoc assumptions regarding the distribution of magnetic energy and the field geometry. In this work we provide the first full magnetohydrodynamic numerical simulations of windwind collision in massive binary systems. We study the freefree emission characterizing its dependence on the stellar and orbital parameters. We also study self-consistently the evolution of the magnetic field at the shock region, obtaining also the synchrotron energy distribution integrated along different lines of sight. We show that the magnetic field in the shocks is larger than that obtained when the proportionality between B and the plasma density is assumed. Also, we show that the role of the synchrotron emission relative to the total radio emission has been underestimated.

Architectures, stability and optimization for clock distribution networks

Synchronous telecommunication networks, distributed control systems and integrated circuits have its accuracy of operation dependent on the existence of a reliable time basis signal extracted from the line data stream and acquirable to each node. In this sense, the existence of a sub-network (inside the main network) dedicated to the distribution of the clock signals is crucially important. There are different solutions for the architecture of the time distribution sub-network and choosing one of them depends on cost, precision, reliability and operational security. In this work we expose: (i) the possible time distribution networks and their usual topologies and arrangements. (ii) How parameters of the network nodes can affect the reachability and stability of the synchronous state of a network. (iii) Optimizations methods for synchronous networks which can provide low cost architectures with operational precision, reliability and security. (C) 2011 Elsevier B. V. All rights reserved.

Modeling the emergence of HIV-1 drug resistance resulting from antiretroviral therapy: Insights from theoretical and numerical studies

The use of antiretroviral therapy has proven to be remarkably effective in controlling the progression of human immunodeficiency virus (HIV) infection and prolonging patient's survival. Therapy however may fail and therefore these benefits can be compromised by the emergence of HIV strains that are resistant to the therapy. In view of these facts, the question of finding the reason for which drug-resistant strains emerge during therapy has become a worldwide problem of great interest. This paper presents a deterministic HIV-1 model to examine the mechanisms underlying the emergence of drug-resistance during therapy. The aim of this study is to determine whether, and how fast, antiretroviral therapy may determine the emergence of drug resistance by calculating the basic reproductive numbers. The existence, feasibility and local stability of the equilibriums are also analyzed. By performing numerical simulations we show that Hopf bifurcation may occur. The model suggests that the individuals with drug-resistant infection may play an important role in the epidemic of HIV. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

Stability analysis and area spectrum of three-dimensional Lifshitz black holes

In this work, we probe the stability of a z = 3 three-dimensional Lifshitz black hole by using scalar and spinorial perturbations. We found an analytical expression for the quasinormal frequencies of the scalar probe field, which perfectly agree with the behavior of the quasinormal modes obtained numerically. The results for the numerical analysis of the spinorial perturbations reinforce the conclusion of the scalar analysis, i.e., the model is stable under scalar and spinor perturbations. As an application we found the area spectrum of the Lifshitz black hole, which turns out to be equally spaced.

Structural stability of flexible lines in catenary configuration under torsion

Catenary risers can present during installation a very low tension close to seabed, which combined with torsion moment can lead to a structural instability, resulting in a loop. This is undesirable once it is possible that the loop turns into a kink, creating damage. This work presents a numerical methodology to analyze the conditions of loop formation in catenary risers. Stability criteria were applied to finite element models, including geometric nonlinearities and contact constraint due to riser-seabed interaction. The classical Greenhill's formula was used to predict the phenomenon and parametric analysis shows a “universal plot” able to predict instability in catenaries using a simple equation that can be applied for typical risers installation conditions and, generically, for catenary lines under torsion.