Determination of the critical coupling for oscillators in a ring

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

Chaotic thermalization in Yang-Mills-Higgs theory on a spacial lattice

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

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.

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.

Effect of Jupiter's mass growth on satellite capture

We study the effects of Jupiter mass growth in order to permanently capture prograde satellites. Adopting the restricted three-body problem, Sun-Jupiter-Particle, we performed numerical simulations backward in time while considering the decrease in Jupiter's mass. We considered the particle's initial conditions to be prograde, at pericenter, in the region 100R(4) <= a <= 400R(4) and 0 <= e <= 0.5. The results give Jupiter's mass at the moment when the particle escapes from the planet. Such values give an indication of the conditions that are necessary for capture. An analysis of these results shows that prograde satellite capture is more complex than a retrograde one. It occurs in a two-step process. First, when the particles get inside about 0.85R(Hill) (Hills' radius), they become weakly bound to Jupiter. Then, they keep migrating toward the planet with a strong decrease in eccentricity, while the planet is growing. The radial oscillation of the particles reduces significantly when they reach a radial distance that is less than about 0.45R(Hill) from the planet. Three-dimensional simulations for the known prograde satellites of Jupiter were performed. The results indicate that Leda, Himalia, Lysithea, and Elara could have been permanently captured when Jupiter had between 50% and 60% of its present mass.

The effect of Jupiter's mass growth on satellite capture - Retrograde case

Gravitational capture can be used to explain the existence of the irregular satellites of giants planets. However, it is only the first step since the gravitational capture is temporary. Therefore, some kind of non-conservative effect is necessary to to turn the temporary capture into a permanent one. In the present work we study the effects of Jupiter mass growth for the permanent capture of retrograde satellites. An analysis of the zero velocity curves at the Lagrangian point L-1 indicates that mass accretion provides an increase of the confinement region ( delimited by the zero velocity curve, where particles cannot escape from the planet) favoring permanent captures. Adopting the restricted three-body problem, Sun-Jupiter-Particle, we performed numerical simulations backward in time considering the decrease of M-4. We considered initial conditions of the particles to be retrograde, at pericenter, in the region 100 R-4 less than or equal to a less than or equal to 400 R-4 and 0 less than or equal to e < 0.5. The results give Jupiter's mass at the moment when the particle escapes from the planet. Such values are an indication of the necessary conditions that could provide capture. An analysis of these results shows that retrograde satellites would be captured as soon as they get inside the Hills' radius and after that they keep migrating toward the planet while it is growing. For the region where the orbits of the four old retrograde satellites of Jupiter ( Ananke, Carme, Pasiphae and Sinope) are located we found that such satellites could have been permanently captured when Jupiter had between 62% and 93% of its present mass.

On energy transfer between linear and nonlinear oscillators

In this paper energy transfer in a dissipative mechanical system is analysed. Such system is composed of a linear and a nonlinear oscillator with a nonlinearizable cubic stiffness. Depending on initial conditions, we find energy transfer either from linear to nonlinear oscillator (energy pumping) or from nonlinear to linear. Such results are valid for two different potentials. However, under resonance and absence of external excitation, if the mass of the nonlinear oscillator is adequately small then the linear oscillator always loses energy. Our approach uses rigorous Regular Perturbation Theory. Besides, we have included the case of two linear oscillators under linear or cubic interactions. Comparisons with the earlier case are made. (c) 2008 Elsevier Ltd. All rights reserved.

In-flight and collisional dissipation as a mechanism to suppress Fermi acceleration in a breathing Lorentz gas

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

Avaliação da silagem de diferentes híbridos de milho (Zea mays, L.) por meio do desempenho de bezerros confinados em fase de crescimento

Objetivou-se avaliar, por meio do desempenho animal, silagens de três híbridos de milho (Zea mays, L.): Agroceres AG-5011, Braskalb XL-344 e Cargill C-806. Foram utilizados 27 machos, não-castrados, mestiços Nelore x Charolês, com idade média de 10 meses e peso médio inicial de 252 kg, confinados durante 84 dias, subdivididos em quatro períodos de 21 dias. As silagens dos híbridos testados corresponderam a 70% da matéria seca ofertada. Não houve interação do híbrido de milho testado com o período de avaliação em confinamento, para os parâmetros estudados. A dieta que incluiu a silagem do híbrido AG-5011 resultou em maior consumo voluntário médio diário de matéria seca, expresso por 100 kg de peso vivo e tamanho metabólico, e maior consumo médio diário de fibra em detergente ácido, expresso nas diferentes formas, em relação aos híbridos XL-344 e C-806, que, por sua vez, foram semelhantes entre si. As dietas que incluíram silagens dos híbridos AG-5011, XL-344 e C-806 apresentaram semelhanças quanto ao consumo voluntário de matéria seca, expresso em kg/animal/dia (6,73; 6,05; e 6,19), consumo médio diário de FDN, em kg/dia (3,054; 2,837; e 2,836), e energia digestível, expressa nas diferentes formas, ganho de peso médio diário, em kg/dia (1,283; 1,219; e 1,249), ganho em estado corporal, conversão alimentar (5,26; 4,98; e 4,97) e eficiência energética (14,02; 13,68; e 13,21). Os três híbridos são indicados como materiais genéticos de muito boa qualidade para produção de silagem, promovendo ganhos de peso diário superior a 1,2 kg/animal/dia, quando utilizados em dietas com relação volumoso:concentrado de 70:30.

A SIMPLE NUMERICAL-METHOD FOR HYDROSTATIC INCOMPRESSIBLE MODELS WITH RIGID LIDS

A simple and easily implemented method is developed to keep the vertical velocity equal to zero at the bottom and top of hydrostatic incompressible numerical models. The pressure is computed at the top by correcting its value given in the previous time step so that the vertical integral of the horizontal divergence is zero at each column. Numerical experiments that exhibit small time variations of pressure at the top are able to simplify the algorithm and save computer time. Numerical simulations illustrate the method effectiveness for a horizontal deformation-induced frontogenesis.

Solitons in Bose-Einstein condensates trapped in a double-well potential

We investigate, analytically and numerically, families of bright solitons in a system of two linearly coupled nonlinear Schrodinger/Gross-Pitaevskii equations, describing two Bose-Einstein condensates trapped in an asymmetric double-well potential, in particular, when the scattering lengths in the condensates have arbitrary magnitudes and opposite signs. The solitons are found to exist everywhere where they are permitted by the dispersion law. Using the Vakhitov-Kolokolov criterion and numerical methods, we show that, except for small regions in the parameter space, the solitons are stable to small perturbations. Some of them feature self-trapping of almost all the atoms in the condensate with no atomic interaction or weak repulsion is coupled to the self-attractive condensate. An unusual bifurcation is found, when the soliton bifurcates from the zero solution with vanishing amplitude and width simultaneously diverging but at a finite number of atoms in the soliton. By means of numerical simulations, it is found that, depending on values of the parameters and the initial perturbation, unstable solitons either give rise to breathers or completely break down into incoherent waves (radiation). A version of the model with the self-attraction in both components, which applies to the description of dual-core fibers in nonlinear optics, is considered too, and new results are obtained for this much studied system. (C) 2003 Elsevier B.V. All rights reserved.

A critical investigation of the Tanford-Kirkwood scheme by means of Monte Carlo simulations

Monte Carlo simulations are used to assess the adequacy of the Tanford-Kirkwood prescription for electrostatic interactions in macromolecules. Within a continuum dielectric framework, the approach accurately describes salt screening of electrostatic interactions for moderately charged systems consistent with common proteins at physiological conditions. The limitations of the Debye-Huckel theory, which forms the statistical mechanical basis for the Tanford-Kirkwood result, become apparent for highly charged systems. It is shown, both by an analysis of the Debye-Huckel theory and by numerical simulations, that the difference in dielectric permittivity between macromolecule and surrounding solvent does not play a significant role for salt effects if the macromolecule is highly charged. By comparison to experimental data, the continuum dielectric model (combined with either an approximate effective Hamiltonian as in the Tanford-Kirkwood treatment or with exact Monte Carlo simulations) satisfactorily predicts the effects of charge mutation on metal ion binding constants, but only if the macromolecule and solvent are assigned the same or similar permittivities.

The strands of the F ring disturbed by its closest satellites

Some Voyager images showed that the F ring of Saturn is composed of at least four separate, non-intersecting, strands covering about 45 degrees in longitude. According to Murray et al. [Murray, C.D., Gordon, M., Giuliatti Winter, S.M. Unraveling the strands of Saturn's F ring. Icarus 129, 304, 1997.] this structure may be caused by undetected satellites embedded in the gaps.Due to precession, the satellites Prometheus and Pandora and the ring particles can experience periodic close encounters. Giuliatti Winter et al. [Giuliatti Winter, S.M, Murray, C.D., Gordon, M. Perturbations to Saturn's F-ring strands at their closest approach to Prometheus. Plan. Space Sciences, 48, 817, 2000.] analysed the behaviour of these four strands at closest approach with the satellite Prometheus. Their work suggests that Prometheus can induce the ring particles to scatter in the direction of the planet, thus increasing the population of small bodies in this region.In this work we analysed the effects of Prometheus on the radial structure of Saturn's F ring during the Voyager and early Cassini epochs. Our results show that at Voyager epoch Prometheus, and also Pandora, had a negligible influence in the strands. However, during the Cassini encounter Prometheus could affect the strands significantly, scattering particles of the inner strand in the direction of the planet. This process can contribute to the replenishment of material in the region between the F ring and the A ring, where two rings have recently been discovered [Porco, C. et al. Cassini imaging science. Initial results on Saturn's rings and small Satellites. Science, 307, 1226, 2005].We also analyse the behaviour of undetected satellites under the effects of these two satellites by computing the Lyapunov Characteristic Exponent. Our results show that these satellites have a chaotic behaviour which leads to a much more complex scenario. The new satellite S/2004 S6 also presents a chaotic behaviour with can alter the dynamic of the system, since this satellite crosses the orbit of the strands. (C) 2006 COSPAR. Published by Elsevier Ltd. All rights reserved.

Third body perturbation using a single averaged model considering elliptic orbits for the disturbing body

In the present work, we expanded the study done by Solorzanol(1) including the eccentricity of the perturbing body. The assumptions used to develop the single-averaged analytical model are the same ones of the restricted elliptic three-body problem. The disturbing function was expanded in Legendre polynomials up to fourth-order. After that, the equations of motion are obtained from the planetary equations and we performed a set of numerical simulations. Different initial eccentricities for the perturbing and perturbed body are considered. The results obtained perform an analysis of the stability of a near-circular orbits and investigate under which conditions this orbit remain near-circular.

Sphere of influence and gravitational capture radius: a dynamical approach

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