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.

Probing the properties of Be star discs with spectroastrometry and NLTE radiative transfer modelling: beta CMi

While the presence of discs around classical Be stars is well established, their origin is still uncertain. To understand what processes result in the creation of these discs and how angular momentum is transported within them, their physical properties must be constrained. This requires comparing high spatial and spectral resolution data with detailed radiative transfer modelling. We present a high spectral resolution, R similar to 80 000, sub-milliarcsecond precision, spectroastrometric study of the circumstellar disc around the Be star beta CMi. The data are confronted with 3D, non-local thermodynamic equilibrium radiative transfer calculations to directly constrain the properties of the disc. Furthermore, we compare the data to disc models featuring two velocity laws: Keplerian, the prediction of the viscous disc model, and angular momentum conserving rotation. It is shown that the observations of beta CMi can only be reproduced using Keplerian rotation. The agreement between the model and the observed spectral energy distribution, polarization and spectroastrometric signature of beta CMi confirms that the discs around Be stars are well modelled as viscous decretion discs.

Beyond the diffraction limit of optical/IR interferometers I. Angular diameter and rotation parameters of Achernar from differential phases

Context. Spectrally resolved long-baseline optical/IR interferometry of rotating stars opens perspectives to investigate their fundamental parameters and the physical mechanisms that govern their interior, photosphere, and circumstellar envelope structures. Aims. Based on the signatures of stellar rotation on observed interferometric wavelength-differential phases, we aim to measure angular diameters, rotation velocities, and orientation of stellar rotation axes. Methods. We used the AMBER focal instrument at ESO-VLTI in its high-spectral resolution mode to record interferometric data on the fast rotator Achernar. Differential phases centered on the hydrogen Br gamma line (K band) were obtained during four almost consecutive nights with a continuous Earth-rotation synthesis during similar to 5h/night, corresponding to similar to 60 degrees position angle coverage per baseline. These observations were interpreted with our numerical code dedicated to long-baseline interferometry of rotating stars. Results. By fitting our model to Achernar's differential phases from AMBER, we could measure its equatorial radius R-eq = 11.6 +/- 0.3 R-circle dot, equatorial rotation velocity V-eq = 298 +/- 9 km s(-1), rotation axis inclination angle i = 101.5 +/- 5.2 degrees, and rotation axis position angle (from North to East) PA(rot) = 34.9 +/- 1.6 degrees. From these parameters and the stellar distance, the equatorial angular diameter circle divide(eq) of Achernar is found to be 2.45 +/- 0.09 mas, which is compatible with previous values derived from the commonly used visibility amplitude. In particular, circle divide(eq) and PA(rot) measured in this work with VLTI/AMBER are compatible with the values previously obtained with VLTI/VINCI. Conclusions. The present paper, based on real data, demonstrates the super-resolution potential of differential interferometry for measuring sizes, rotation velocities, and orientation of rotating stars in cases where visibility amplitudes are unavailable and/or when the star is partially or poorly resolved. In particular, we showed that differential phases allow the measurement of sizes up to similar to 4 times smaller than the diffraction-limited angular resolution of the interferometer.

A regional multirnodulus algorithm for blind equalization of QAM signals: Introduction and steady-state analysis

It is well known that constant-modulus-based algorithms present a large mean-square error for high-order quadrature amplitude modulation (QAM) signals, which may damage the switching to decision-directed-based algorithms. In this paper, we introduce a regional multimodulus algorithm for blind equalization of QAM signals that performs similar to the supervised normalized least-mean-squares (NLMS) algorithm, independently of the QAM order. We find a theoretical relation between the coefficient vector of the proposed algorithm and the Wiener solution and also provide theoretical models for the steady-state excess mean-square error in a nonstationary environment. The proposed algorithm in conjunction with strategies to speed up its convergence and to avoid divergence can bypass the switching mechanism between the blind mode and the decision-directed mode. (c) 2012 Elsevier B.V. All rights reserved.

Planetary nebulae and determination of the bulge-disk boundary

In this paper, a sample of planetary nebulae in the Galaxy's inner-disk and bulge is used to find the galactocentric distance that optimally separates these two populations in terms of their abundances. Statistical distance scales were used to investigate the distribution of abundances across the disk–bulge interface, while a Kolmogorov–Smirnov test was used to find the distance at which the chemical properties of these regions separate optimally. The statistical analysis indicates that, on average, the inner population is characterized by lower abundances than the outer component. Additionally, for the α-element abundances, the inner population does not follow the disk's radial gradient toward the Galactic Center. Based on our results, we suggest a bulge–disk interface at 1.5 kpc, marking the transition between the bulge and the inner disk of the Galaxy as defined by the intermediate-mass population.

Ventilation-perfusion ratio in perflubron during partial liquid ventilation

BACKGROUND: Functional magnetic resonance imaging (fMRI) of fluorine-19 allows for the mapping of oxygen partial pressure within perfluorocarbons in the alveolar space (Pao(2)). Theoretically, fMRI-detected Pao(2) can be combined with the Fick principle approach, i.e., a mass balance of oxygen uptake by ventilation and delivery by perfusion, to quantify the ventilation-perfusion ratio (Va/Q) of a lung region: The mixed venous blood and the inspiratory oxygen fraction, which are equal for all lung regions, are measured. In addition, the local expiratory oxygen fraction and the end capillary oxygen content, both of which may differ between the lung regions, are calculated using the fMRI-detected Pao(2). We investigated this approach by numerical simulations and applied it to quantify local Va/Q in the perfluorocarbons during partial liquid ventilation. METHODS: Numerical simulations were performed to analyze the sensitivity of the Va/Q calculation and to compare this approach with another one proposed by Rizi et al. in 2004 (Magn Reson Med 2004;52:65-72). Experimentally, the method was used during partial liquid ventilation in 7 anesthetized pigs. The Pao(2) distribution in intraalveolar perflubron was measured by fluorine-19 MRI. Respiratory gas fractions together with arterial and mixed venous blood samples were taken to quantify oxygen partial pressure and content. Using the Fick principle, the local Va/Q was estimated. The impact of gravity (nondependent versus dependent) of perflubron dose (10 vs 20 mL/kg body weight) and of inspired oxygen fraction (Fio(2)) (0.4-1.0) on Va/Q was examined. RESULTS: In numerical simulations, the Fick principle proved to be appropriate over the Va/Q range from 0.02 to 2.5. Va/Q values were in acceptable agreement with the method published by Rizi et al. In the experimental setting, low mean Va/Q values were found in perflubron (confidence interval [CI] 0.08-0.29 with 20 mL/kg perflubron). At this dose, Va/Q in the nondependent lung was higher (CI 0.18-0.39) than in the dependent lung regions (CI 0.06-0.16; P = 0.006; Student t test). Differences depending on Fio(2) or perflubron dose were, however, small. CONCLUSION: The results show that derivation of Va/Q from local Po(2) measurements using fMRI in perflubron is feasible. The low detected Va/Q suggests that oxygen transport into the perflubron-filled alveolar space is significantly restrained.

On the analytical approach for modeling photovoltaic systems behavior

The 1-diode/2-resistors electric circuit equivalent to a photovoltaic system is analyzed. The equations at particular points of the I–V curve are studied considering the maximum number of terms. The maximum power point as a boundary condition is given special attention. A new analytical method is developed based on a reduced amount of information, consisting in the normal manufacturer data. Results indicate that this new method is faster than numerical methods and has similar (or better) accuracy than other existing methods, numerical or analytical.

Comportamento tribo-mecânico e desgaste adesivo de materiais em nanoescala: análises por dinâmica molecular e mecânica do contínuo.

As formulações baseadas na mecânica do contínuo, embora precisas até certo ponto, por vezes não podem ser utilizadas, ou não são conceitualmente corretas para o entendimento de fenômenos em escalas reduzidas. Estas limitações podem aparecer no estudo dos fenômenos tribológicos em escala nanométrica, que passam a necessitar de novos métodos experimentais, teóricos e computacionais que permitam explorar estes fenômenos com a resolução necessária. Simulações atomísticas são capazes de descrever fenômenos em pequena escala, porém, o número necessário de átomos modelados e, portanto, o custo computacional - geralmente torna-se bastante elevado. Por outro lado, os métodos de simulação associados à mecânica do contínuo são mais interessantes em relação ao custo computacional, mas não são precisos na escala atômica. A combinação entre essas duas abordagens pode, então, permitir uma compreensão mais realista dos fenômenos da tribologia. Neste trabalho, discutem-se os conceitos básicos e modelos de atrito em escala atômica e apresentam-se estudos, por meio de simulação numérica, para a análise e compreensão dos mecanismos de atrito e desgaste no contato entre materiais. O problema é abordado em diferentes escalas, e propõe-se uma abordagem conjunta entre a Mecânica do Contínuo e a Dinâmica Molecular. Para tanto, foram executadas simulações numéricas, com complexidade crescente, do contato entre superfícies, partindo-se de um primeiro modelo que simula o efeito de defeitos cristalinos no fenômeno de escorregamento puro, considerando a Dinâmica Molecular. Posteriormente, inseriu-se, nos modelos da mecânica do contínuo, considerações sobre o fenômeno de adesão. A validação dos resultados é feita pela comparação entre as duas abordagens e com a literatura.

Simulated magnetic field expulsion in neutron star cores

The study of long-term evolution of neutron star (NS) magnetic fields is key to understanding the rich diversity of NS observations, and to unifying their nature despite the different emission mechanisms and observed properties. Such studies in principle permit a deeper understanding of the most important parameters driving their apparent variety, e.g. radio pulsars, magnetars, X-ray dim isolated NSs, gamma-ray pulsars. We describe, for the first time, the results from self-consistent magnetothermal simulations considering not only the effects of the Hall-driven field dissipation in the crust, but also adding a complete set of proposed driving forces in a superconducting core. We emphasize how each of these core-field processes drive magnetic evolution and affect observables, and show that when all forces are considered together in vectorial form, the net expulsion of core magnetic flux is negligible, and will have no observable effect in the crust (consequently in the observed surface emission) on megayear time-scales. Our new simulations suggest that strong magnetic fields in NS cores (and the signatures on the NS surface) will persist long after the crustal magnetic field has evolved and decayed, due to the weak combined effects of dissipation and expulsion in the stellar core.

Adaptive stepsize based on control theory for stochastic differential equations

The numerical solution of stochastic differential equations (SDEs) has been focussed recently on the development of numerical methods with good stability and order properties. These numerical implementations have been made with fixed stepsize, but there are many situations when a fixed stepsize is not appropriate. In the numerical solution of ordinary differential equations, much work has been carried out on developing robust implementation techniques using variable stepsize. It has been necessary, in the deterministic case, to consider the best choice for an initial stepsize, as well as developing effective strategies for stepsize control-the same, of course, must be carried out in the stochastic case. In this paper, proportional integral (PI) control is applied to a variable stepsize implementation of an embedded pair of stochastic Runge-Kutta methods used to obtain numerical solutions of nonstiff SDEs. For stiff SDEs, the embedded pair of the balanced Milstein and balanced implicit method is implemented in variable stepsize mode using a predictive controller for the stepsize change. The extension of these stepsize controllers from a digital filter theory point of view via PI with derivative (PID) control will also be implemented. The implementations show the improvement in efficiency that can be attained when using these control theory approaches compared with the regular stepsize change strategy. (C) 2004 Elsevier B.V. All rights reserved.