Numerical assessment of mass conservation on a MAC-type method for viscoelastic free surface flows

A numerical study of mass conservation of MAC-type methods is presented, for viscoelastic free-surface flows. We use an implicit formulation which allows for greater time steps, and therefore time marching schemes for advecting the free surface marker particles have to be accurate in order to preserve the good mass conservation properties of this methodology. We then present an improvement by using a Runge-Kutta scheme coupled with a local linear extrapolation on the free surface. A thorough study of the viscoelastic impacting drop problem, for both Oldroyd-B and XPP fluid models, is presented, investigating the influence of timestep, grid spacing and other model parameters to the overall mass conservation of the method. Furthermore, an unsteady fountain flow is also simulated to illustrate the low mass conservation error obtained.

DYNAMICS AT INFINITY of A CUBIC CHUA'S SYSTEM

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

Attitude stability of artificial satellites subject to gravity gradient torque

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

Third-body perturbation in the case of elliptic orbits for the disturbing body

This work presents a semi-analytical and numerical study of the perturbation caused in a spacecraft by a third-body using a double averaged analytical model with the disturbing function expanded in Legendre polynomials up to the second order. The important reason for this procedure is to eliminate terms due to the short periodic motion of the spacecraft and to show smooth curves for the evolution of the mean orbital elements for a long-time period. The aim of this study is to calculate the effect of lunar perturbations on the orbits of spacecrafts that are traveling around the Earth. An analysis of the stability of near-circular orbits is made, and a study to know under which conditions this orbit remains near circular completes this analysis. A study of the equatorial orbits is also performed. Copyright (C) 2008 R. C. Domingos et al.

Dynamic critical behavior of percolation observables in the 2d Ising model

We present preliminary results of our numerical study of the critical dynamics of percolation observables for the two-dimensional Ising model. We consider the (Monte-Carlo) short-time evolution of the system obtained with a local heat-bath method and with the global Swendsen-Wang algorithm. In both cases, we find qualitatively different dynamic behaviors for the magnetization and Omega, the order parameter of the percolation transition. This may have implications for the recent attempts to describe the dynamics of the QCD phase transition using cluster observables.

Localization of the optimal solution and a posteriori bounds for aggregation

After an aggregated problem has been solved, it is often desirable to estimate the accuracy loss due to the fact that a simpler problem than the original one has been solved. One way of measuring this loss in accuracy is the difference in objective function values. To get the bounds for this difference, Zipkin (Operations Research 1980;28:406) has assumed, that a simple (knapsack-type) localization of an original optimal solution is known. Since then various extensions of Zipkin's bound have been proposed, but under the same assumption. A method to compute the bounds for variable aggregation for convex problems, based on general localization of the original solution is proposed. For some classes of the original problem it is shown how to construct the localization. Examples are given to illustrate the main constructions and a small numerical study is presented.

Resonance and chaos .1. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted three-body problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincare surfaces of section with a mass ratio of 10(-3) (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. on theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

The defocusing mechanism of isochronous resonances

We present a numerical study concerning the defocusing mechanism of isochronous resonance island chains in the presence of two permanent robust tori. The process is initialized and concluded through bifurcations of fixed points located on the robust tori. Our approach is based on a Hamiltonian system derived from the resonant normal form. Choosing a convenient parameter in this system, we are able to depict a comprehensive analysis of the dynamics of the problem. (c) 2004 Elsevier B.V. All rights reserved.

Lagrangian heuristic for a class of the generalized assignment problems

A Lagrangian based heuristic is proposed for many-to-many assignment problems taking into account capacity limits for task and agents. A modified Lagrangian bound studied earlier by the authors is presented and a greedy heuristic is then applied to get a feasible Lagrangian-based solution. The latter is also used to speed up the subgradient scheme to solve the modified Lagrangian dual problem. A numerical study is presented to demonstrate the efficiency of the proposed approach. (C) 2010 Elsevier Ltd. All rights reserved.

Sphere of influence and gravitational capture radius: a dynamical approach

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

Resonance and chaos: I. First-order interior resonances

Analytical models for studying the dynamical behaviour of objects near interior, mean motion resonances are reviewed in the context of the planar, circular, restricted threebody problem. The predicted widths of the resonances are compared with the results of numerical integrations using Poincaré surfaces of section with a mass ratio of 10-3 (similar to the Jupiter-Sun case). It is shown that for very low eccentricities the phase space between the 2:1 and 3:2 resonances is predominantly regular, contrary to simple theoretical predictions based on overlapping resonance. A numerical study of the 'evolution' of the stable equilibrium point of the 3:2 resonance as a function of the Jacobi constant shows how apocentric libration at the 2:1 resonance arises; there is evidence of a similar mechanism being responsible for the centre of the 4:3 resonance evolving towards 3:2 apocentric libration. This effect is due to perturbations from other resonances and demonstrates that resonances cannot be considered in isolation. On theoretical grounds the maximum libration width of first-order resonances should increase as the orbit of the perturbing secondary is approached. However, in reality the width decreases due to the chaotic effect of nearby resonances.

Free expansion of attractive and repulsive Bose-Einstein condensed vortex states

A numerical study of the time-dependent Gross-Pitaevskii equation for an axially symmetric trap to obtain insight into the free expansion of vortex states of BEC is presented. As such, the ratio of vortex-core radius to radia rms radius xc/xrms(<1) is found to play an interesting role in the free expansion of condensed vortex states. the larger this ratio, the more prominent is the vortex core and the easier is the possibility of experimental detection of vortex states.

Preheating in quintessential inflation

We perform a numerical study of the preheating mechanism of particle production in models of quintessential inflation and compare it with the usual gravitational production mechanism. We find that even for a very small coupling between the inflaton field and a massless scalar field, g ≳ 10 -6, preheating dominates over gravitational particle production. Reheating temperatures in the range 10 4 ≲ T rh ≲ 10 15 GeV can be easily obtained. © 2003 Published by Elsevier B.V.

Quintessential inflation and non-oscillatory reheating model

A numerical study of the non-oscillatory reheating mechanism in a quintessential inflation context shows that high reheating temperature can be achieved compared with the usual reheating mechanism in which particles are produced gravitationally. We find that even for a very small coupling between the inflaton field and a massless scalar field, the non-oscillatory reheating production of particles dominates over the gravitational production mechanism. © 2004 Published by Elsevier B.V.

Mean-field model of jet formation in a collapsing Bose-Einstein condensate

We perform a systematic numerical study, based on the time-dependent Gross-Pitaevskii equation, of jet formation in collapsing and exploding Bose-Einstein condensates as in the experiment by Donley et al (2001 Nature 412 295). In the actual experiment, via a Feshbach resonance, the scattering length of atomic interaction was suddenly changed from positive to negative on a pre-formed condensate. Consequently, the condensate collapsed and ejected atoms via explosion. On disruption of the collapse by suddenly changing the scattering length to zero, a radial jet of atoms was formed in the experiment. We present a satisfactory account of jet formation under the experimental conditions and also make predictions beyond experimental conditions which can be verified in future experiments.