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.

Numerical correction for the diffuse solar irradiance by the melo-escobedo shadowring measuring method

The present paper deals with numerical corrections factors proposed as a function of the clearness index in order to correct the diffuse solar irradiance measured with the Melo-Escobedo Shadowring Measuring Method (ME shadowring). The global irradiance was measured by an Eppley - PSP pyranometer ; direct normal irradiance by an Eppley-NIP pyrheliometer fitted to a ST-3 sun tracking device and the diffuse irradiance by an Eppley-PSP pyranometer fitted to a ME shadowring. The validations were performed by the MBE and RMSE statistical indicators. The results showed that the numerical correction factors were appropriate to correct the shadowring diffuse irradiance.

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.

A Numerical Study of the FENE-CR Model Applied to a Jet Flow Problem

The FENE-CR model is investigated through a numerical algorithm to simulate the time-dependent moving free surface flow produced by a jet impinging on a flat surface. The objective is to demonstrate that by increasing the extensibility parameter L, the numerical solutions converge to the solutions obtained with the Oldroyd-B model. The governing equations are solved by an established free surface flow solver based on the finite difference and marker-and-cell methods. Numerical predictions of the extensional viscosity obtained with several values of the parameter L are presented. The results show that if the extensibility parameter L is sufficiently large then the extensional viscosities obtained with the FENE-CR model approximate the corresponding Oldroyd-B viscosity. Moreover, the flow from a jet impinging on a flat surface is simulated with various values of the extensibility parameter L and the fluid flow visualizations display convergence to the Oldroyd-B jet flow results.

Short Lyapunov time: a method for identifying confined chaos

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

SPH simulations of clumps formation by dissipative collisions of molecular clouds - II. Magnetic case

We performed computer simulations of interstellar cloud-cloud collisions using the three-dimensional smoothed particle magnetohydrodynamics method. In order to study the role of the magnetic field on the process of collision-triggered fragmentation, we focused our attention on head-on supersonic collisions between two identical spherical molecular-clouds. Two extreme configurations of the magnetic field were adopted: parallel and perpendicular to the initial clouds motion. The initial magnetic field strength was approximately 12.0 muG. In the parallel case, much more of the collision debris were retained in the shocking region than in the non-magnetic case where gas escaped freely throughout the symmetry plane. Differently from the non-magnetic case, eddy-like vortices were formed. The regions of highest vorticity and the the regions of highest density are offset. We found clumps formation only in the parallel case, however, they were larger, hotter and less dense than in the analogous non-magnetic case. In the perpendicular case, the compressed field works as a magnetic wall, preventing a stronger compression of the colliding clouds. This last effect inhibits direct contact of the two clouds. In both cases, we found that the field lines show a chaotic aspect in large scales. Also, the field magnitude is considerably amplified in the shock layer. However, the field distribution is almost coherent in the higher density regions.

SPH simulations of clumps formation by dissipative collision of molecular clouds - I. Non magnetic case

Computer experiments of interstellar cloud collisions were performed with a new smoothed-particle-hydrodynamics (SPH) code. The SPH quantities were calculated by using spatially adaptive smoothing lengths and the SPH fluid equations of motion were solved by means of a hierarchical multiple time-scale leapfrog. Such a combination of methods allows the code to deal with a large range of hydrodynamic quantities. A careful treatment of gas cooling by H, H(2), CO and H II, as well as a heating mechanism by cosmic rays and by H(2) production on grains surface, were also included in the code. The gas model reproduces approximately the typical environment of dark molecular clouds. The experiments were performed by impinging two dynamically identical spherical clouds onto each other with a relative velocity of 10 km s(-1) but with a different impact parameter for each case. Each object has an initial density profile obeying an r(-1)-law with a cutoff radius of 10 pc and with an initial temperature of 20 K. As a main result, cloud-cloud collision triggers fragmentation but in expense of a large amount of energy dissipated, which occurred in the head-on case only. Off-center collision did not allow remnants to fragment along the considered time (similar to 6 Myr). However, it dissipated a considerable amount of orbital energy. Structures as small as 0.1 pc, with densities of similar to 10(4) cm(-3), were observed in the more energetic collision.

Moonlets wandering on a leash-ring

Since the Voyager flybys, embedded moonlets have been proposed to explain some of the surprising structures observed in Saturn's narrow F ring. Experiments conducted with the Cassini spacecraft support this suggestion. Images of the F ring show bright compact spots, and seven occultations of stars by the F ring, monitored by ultraviolet and infrared experiments, revealed nine events of high optical depth. These results point to a large number of such objects, but it is not clear whether they are solid moonlets or rather loose particle aggregates. Subsequent images suggested an irregular motion of these objects so that a determination of their orbits consistent with the F ring failed. Some of these features seem to cross the whole ring. Here we show that these observations are explained by chaos in the F ring driven mainly by the 'shepherd' moons Prometheus and Pandora. It is characterized by a rather short Lyapunov time of about a few hundred orbital periods. Despite this chaotic diffusion, more than 93 per cent of the F-ring bodies remain confined within the F ring because of the shepherding, but also because of a weak radial mobility contrasted by an effective longitudinal diffusion. This chaotic stirring of all bodies involved prevents the formation of 'propellers' typical of moonlets, but their frequent ring crossings explain the multiple radial 'streaks' seen in the F ring. The related 'thermal' motion causes more frequent collisions between all bodies which steadily replenish F-ring dust and allow for ongoing fragmentation and re-accretion processes (ring recycling).

Sphere of influence and gravitational capture radius: a dynamical approach

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

The role of Mab as a source for the mu ring of Uranus (Research Note)

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

A new solution to the optimal power flow problem

A new approach to solving the Optimal Power Flow problem is described, making use of some recent findings, especially in the area of primal-dual methods for complex programming. In this approach, equality constraints are handled by Newton's method inequality constraints for voltage and transformer taps by the logarithmic barrier method and the other inequality constraints by the augmented Lagrangian method. Numerical test results are presented, showing the effective performance of this algorithm. © 2001 IEEE.

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)

TERRESTRIAL PLANET FORMATION IN A PROTOPLANETARY DISK WITH A LOCAL MASS DEPLETION: A SUCCESSFUL SCENARIO FOR THE FORMATION OF MARS

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

Long-term dynamics of Methone, Anthe and Pallene

We numerically investigate the long-term dynamics of the Saturn's small satellites Methone (S/2004 S1), Anthe (S/2007 S4) and Pallene (S/2004 S2). In our numerical integrations, these satellites are disturbed by non-spherical shape of Saturn and the six nearest regular satellites. The stability of the small bodies is studied here by analyzing long-term evolution of their orbital elements.We show that long-term evolution of Pallene is dictated by a quasi secular resonance involving the ascending nodes (12) and longitudes of pericentric distances (pi) of Mimas (subscript 1) and Pallene (subscript 2), which critical argument is pi(2) - pi(1) - Omega(1) + Omega(2) Long-term orbital evolution of Methone and Anthe are probably chaotic since: i) their orbits randomly cross the orbit of Mimas in time scales of thousands years); ii) long-term numerical simulations involving both small satellites are strongly affected by small changes in the initial conditions.

Comparison of methodologies for degree-day estimation using numerical methods

O desenvolvimento de projetos relacionados ao desempenho de diversas culturas tem recebido aperfeiçoamento cada vez maior, incorporado a modelos matemáticos sendo indispensável à utilização de equações cada vez mais consistentes que possibilitem previsão e maior aproximação do comportamento real, diminuindo o erro na obtenção das estimativas. Entre as operações unitárias que demandam maior estudo estão aquelas relacionadas com o crescimento da cultura, caracterizadas pela temperatura ideal para o acréscimo de matéria seca. Pelo amplo uso dos métodos matemáticos na representação, análise e obtenção de estimativas de graus-dia, juntamente com a grande importância que a cultura da cana-de-açúcar tem para a economia brasileira, foi realizada uma avaliação dos modelos matemáticos comumente usados e dos métodos numéricos de integração na estimativa da disponibilidade de graus-dia para essa cultura, na região de Botucatu, Estado de São Paulo. Os modelos de integração, com discretização de 6 em 6 h, apresentaram resultados satisfatórios na estimativa de graus-dia. As metodologias tradicionais apresentaram desempenhos satisfatórios quanto à estimativa de grausdia com base na curva de temperatura horária para cada dia e para os agrupamentos de três, sete, 15 e 30 dias. Pelo método numérico de integração, a região de Botucatu, Estado de São Paulo, apresentou disponibilidade térmica anual média de 1.070,6 GD para a cultura da cana-de-açúcar.