Dinâmica do sistema binário Plutão-Caronte

Pós-graduação em Física - FEG

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.

On the effects of each term of the geopotential perturbation along the time I: Quasi-circular orbits

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

Satélites estabilizados por rotação e torque de radiação solar direta

The aims of this work are to analyze the direct solar radiation pressure torque (TPRS) in the rotational motion of spin-stabilized artificial satellites, to numerically implement these solutions and to compare the results with real data of the Brazilian Satellite Data Collection – SCD1 and SCD2, supplied by INPE. The mathematical model for this torque is determined for a cylindrical satellite, and the components of this torque are determined in a fixed system in the satellite. An analytical solution for the spin motion equations is proposed, in which TPRSD does not affect the spin velocity of the satellite. Two approaches are adopted in the numerical implementation of the developed theory: the first one considers the proposed theory and the second introduces a variation in the spin velocity based on its real variation. The results obtained indicate that the solar radiation pressure torque has little influence in the right ascension and declination axis of rotation due to the small dimension of the satellite and altitude in which it is found. To better validate the application of the presented theory, the angular deviation of the spin axis and solar aspect angle were also analyzed. The comparison of the results of the approaches conducted with real data show good precision in the theory, which can be applied in the prediction of the rotational motion of the spin-stabilized artificial satellites, when others external torques are considered besides the direct solar radiation pressure torque

Study of some strategies for disposal of the gnss satellites

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

Regiões de estabilidade no ambiente de um sistema de asteroide triplo

This work was developed from the study by Araujo, R.A.N. et al. Stability regions around the components of the triple system 2001 SN263. (Monthly Notices Of The Royal Astronomical Society, 2012, v. 423(4), 3058-3073 p.) where it was studied the stable and unstable regions system (2001 SN263), which is a triple asteroid system, and these are celestial orbiting our sun. Being close to the Earth is characterized as NEA (Near-Earth Asteroids), asteroids and which periodically approach the Earth's orbit, given that there is great interest in the study and exploitation of these objects, it is the key can carry features that contribute to better understand the process of formation of our solar system. Study the dynamics of bodies that govern those systems proves to be greatly attractive because of the mutual gravitational perturbation of bodies and also by external disturbances. Recently, NEA 2001 SN263 was chosen as a target of Aster mission where a probe is sent for this triple system, appearing therefore the need for obtaining information for characterizing stable regions internal and external to the system, with respect to the effects of radiation pressure. First, this study demonstrated that the integrator used showed satisfactory results of the orbital evolution of bodies in accordance with previous studies and also the characterization of stable and unstable regions brought similar results to the study by Araujo et al. (2012). From these results it was possible to carry out the implementation of the radiation pressure in the system in 2001 SN263, in a region close to the central body, where the simulations were carried out, which brought as a result that the regions before being characterized as stable in unstable true for small particles size from 1 to 5 micrometers. So the next orbital region to the central body and the ... ( Complete abstract click electronic access below)

Regiões de estabilidade no ambiente de um sistema de asteroide triplo

This work was developed from the study by Araujo, R.A.N. et al. Stability regions around the components of the triple system 2001 SN263. (Monthly Notices Of The Royal Astronomical Society, 2012, v. 423(4), 3058-3073 p.) where it was studied the stable and unstable regions system (2001 SN263), which is a triple asteroid system, and these are celestial orbiting our sun. Being close to the Earth is characterized as NEA (Near-Earth Asteroids), asteroids and which periodically approach the Earth's orbit, given that there is great interest in the study and exploitation of these objects, it is the key can carry features that contribute to better understand the process of formation of our solar system. Study the dynamics of bodies that govern those systems proves to be greatly attractive because of the mutual gravitational perturbation of bodies and also by external disturbances. Recently, NEA 2001 SN263 was chosen as a target of Aster mission where a probe is sent for this triple system, appearing therefore the need for obtaining information for characterizing stable regions internal and external to the system, with respect to the effects of radiation pressure. First, this study demonstrated that the integrator used showed satisfactory results of the orbital evolution of bodies in accordance with previous studies and also the characterization of stable and unstable regions brought similar results to the study by Araujo et al. (2012). From these results it was possible to carry out the implementation of the radiation pressure in the system in 2001 SN263, in a region close to the central body, where the simulations were carried out, which brought as a result that the regions before being characterized as stable in unstable true for small particles size from 1 to 5 micrometers. So the next orbital region to the central body and the ... ( Complete abstract click electronic access below)

Resonances in Mie scattering by an inhomogeneous atomic cloud

Despite the quantum nature of the process, collective scattering by dense cold samples of two-level atoms can be interpreted classically describing the sample as a macroscopic object with a complex refractive index. We demonstrate that resonances in Mie theory can be easily observable in the cooperative scattering by tuning the frequency of the incident laser field or the atomic number. The solution of the scattering problem is obtained for spherical atomic clouds who have the parabolic density characteristic of BECs, and the cooperative radiation pressure force calculated exhibits resonances in the cloud displacement for dense clouds. At odds with uniform clouds which show a complex structure including narrow peaks, these densities show resonances, yet only under the form of quite regular and contrasted oscillations. Copyright (C) EPLA, 2012

Untersuchungen zu Mechanismen der Laser-Teilchenbeschleunigung

Auf dem Gebiet der Teilchenbeschleunigung mittels Hochintensitäts-Lasern wurden in der letzten Dekade viele erfolgreiche Entwicklungen hin zu immer höheren Energien und größeren Teilchenzahlen veröffentlicht. In den meisten Fällen wurde der sogenannte TNSA-Prozess (engl. Target-Normal-Sheath-Acceleration (TNSA)) untersucht. Bei diesem Prozess erfolgt die Beschleunigung in dem an der Oberfläche durch Ladungstrennung erzeugten Potential. Ein kaum vermeidbares Problem ist hierbei das resultierende breite Energie-Spektrum der beschleunigten Teilchen. Diese Situation konnte in den letzten Jahren zwar verbessert, aber nicht vollständig gelöst werden. Für Intensitäten größer 10^(20..21) W/cm^2 sagen theoretische Modellrechnungen eine auf dem Lichtdruck basierende Beschleunigung (engl. Radiation-Pressure-Acceleration (RPA)) mit deutlich eingegrenztem, fast monoenergetischem Spektrum voraus. Im Rahmen dieser Arbeit wurde ein Experiment zur Untersuchung dieses Prozesses bei Intensitäten von einigen 10^19 W/cm^2 durchgeführt. Dazu wurden zunächst spezielle Targets entwickelt und als Patent angemeldet, welche den Experimentbedingungen angepasst sind. Durch die Auslegung des experimentellen Aufbaus und der Diagnostiken auf hohe Repetitionsraten, in Verbindung mit einem geeigneten Lasersystem, konnte auf Basis einer Statistik von mehreren Tausend Schüssen ein großer Parameterraum untersucht werden. Untersucht wurden unter anderem die Abhängigkeit von Targetmaterial und Dicke, Intensität, Laserpolarisation und Vorplasmabedingungen. Aus den gewonnenen Daten und Vergleichen mit 2-dimensionalen numerischen Simulationen konnte ein Modell des Beschleunigungsprozesses aufgestellt und durch Vergleich mit den experimentellen Ergebnissen geprüft werden. Dabei wurden klare Indizien für die Existenz eines neuen, nicht feldinduzierten, Beschleunigungsprozesses gefunden. Darüber hinaus wurde zur Polarisationsbeeinflussung ein optisches System entwickelt, das ausschließlich mit reflexiven Elementen arbeitet. Damit konnten viele Nachteile bestehender, auf Verzögerungsplatten beruhender Elemente vermieden, und die Anwendbarkeit bei hohen Laserenergien erreicht werden.

Modeling the heterogeneous ice and gas coma of Comet 103P/Hartley 2

The spectacular images of Comet 103P/Hartley 2 recorded by the Medium Resolution Instrument (MRI) and High Resolution Instrument (HRI) on board of the Extrasolar Planet Observation and Deep Impact Extended Investigation (EPOXI) spacecraft, as the Deep Impact extended mission, revealed that its bi-lobed very active nucleus outgasses volatiles heterogeneously. Indeed, CO2 is the primary driver of activity by dragging out chunks of pure ice out of the nucleus from the sub-solar lobe that appear to be the main source of water in Hartley 2's coma by sublimating slowly as they go away from the nucleus. However, water vapor is released by direct sublimation of the nucleus at the waist without any significant amount of either CO2 or icy grains. The coma structure for a comet with such areas of diverse chemistry differs from the usual models where gases are produced in a homogeneous way from the surface. We use the fully kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J. 685, 659-677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J. 732, 104-120) applied to Comet 103P/Hartley 2 including sublimating icy grains to reproduce the observations made by EPOXI and ground-based measurements. A realistic bi-lobed nucleus with a succession of active areas with different chemistry was included in the model enabling us to study in details the coma of Hartley 2. The different gas production rates from each area were found by fitting the spectra computed using a line-by-line non-LTE radiative transfer model to the HRI observations. The presence of icy grains with long lifetimes, which are pushed anti-sunward by radiation pressure, explains the observed OH asymmetry with enhancement on the night side of the coma.

Assessment of detectability of neutral interstellar deuterium by IBEX observations

Context. The abundance of deuterium in the interstellar gas in front of the Sun gives insight into the processes of filtration of neutral interstellar species through the heliospheric interface and potentially into the chemical evolution of the Galactic gas. Aims: We investigate the possibility of detection of neutral interstellar deuterium at 1 AU from the Sun by direct sampling by the Interstellar Boundary Explorer (IBEX). Methods: Using both previous and the most recent determinations of the flow parameters of neutral gas in the local interstellar cloud (LIC) and an observation-based model of solar radiation pressure and ionization in the heliosphere, we simulated the flux of neutral interstellar D at IBEX for the actual measurement conditions. We assessed the number of interstellar D atom counts expected during the first three years of IBEX operation. We also simulated the observations expected during an epoch of high solar activity. In addition, we calculated the expected counts of D atoms from the thin terrestrial water layer covering the IBEX-Lo conversion surface, sputtered by neutral interstellar He atoms. Results: Most D counts registered by IBEX-Lo are expected to come from the water layer, exceeding the interstellar signal by 2 orders of magnitude. However, the sputtering should stop once the Earth leaves the portion of orbit traversed by interstellar He atoms. We identify seasons during the year when mostly the genuine interstellar D atoms are expected in the signal. During the first 3 years of IBEX operations about 2 detectable interstellar D atoms are expected. This number is comparable to the expected number of sputtered D atoms registered during the same time intervals. Conclusions: The most favorable conditions for the detection occur during low solar activity, in an interval including March and April each year. The detection chances could be improved by extending the instrument duty cycle, say, by making observations in the special deuterium mode of IBEX-Lo.

Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Contribution of Starlette, Stella, and AJISAI to the SLR-derived global reference frame

The contribution of Starlette, Stella, and AJI-SAI is currently neglected when defining the International Terrestrial Reference Frame, despite a long time series of precise SLR observations and a huge amount of available data. The inferior accuracy of the orbits of low orbiting geodetic satellites is the main reason for this neglect. The Analysis Centers of the International Laser Ranging Service (ILRS ACs) do, however, consider including low orbiting geodetic satellites for deriving the standard ILRS products based on LAGEOS and Etalon satellites, instead of the sparsely observed, and thus, virtually negligible Etalons. We process ten years of SLR observations to Starlette, Stella, AJISAI, and LAGEOS and we assess the impact of these Low Earth Orbiting (LEO) SLR satellites on the SLR-derived parameters. We study different orbit parameterizations, in particular different arc lengths and the impact of pseudo-stochastic pulses and dynamical orbit parameters on the quality of the solutions. We found that the repeatability of the East and North components of station coordinates, the quality of polar coordinates, and the scale estimates of the reference are improved when combining LAGEOS with low orbiting SLR satellites. In the multi-SLR solutions, the scale and the Z component of geocenter coordinates are less affected by deficiencies in solar radiation pressure modeling than in the LAGEOS-1/2 solutions, due to substantially reduced correlations between the Z geocenter coordinate and empirical orbit parameters. Eventually, we found that the standard values of Center-of-mass corrections (CoM) for geodetic LEO satellites are not valid for the currently operating SLR systems. The variations of station-dependent differential range biases reach 52 and 25 mm for AJISAI and Starlette/Stella, respectively, which is why estimating station dependent range biases or using station-dependent CoM, instead of one value for all SLR stations, is strongly recommended.This clearly indicates that the ILRS effort to produce CoM corrections for each satellite, which are site-specific and depend on the system characteristics at the time of tracking,is very important and needs to be implemented in the SLR data analysis.

A comment on the article “A collinearity diagnosis of the GNSS geocenter determination” by P. Rebischung, Z. Altamimi, and T. Springer

Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) showed that the geocenter z -component estimated from observations of global navigation satellite systems (GNSS) is strongly correlated to a particular parameter of the solar radiation pressure (SRP) model developed by Beutler et al. (Manuscr Geod 19:367–386, 1994). They analyzed the forces caused by SRP and the impact on the satellites’ orbits. The authors achieved their results using perturbation theory and celestial mechanics. Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026, 2013) also deal with the geocenter determination with GNSS. The authors carried out a collinearity diagnosis of the associated parameter estimation problem. They conclude “without much exaggerating that current GNSS are insensitive to any component of geocenter motion”. They explain this inability by the high degree of collinearity of the geocenter coordinates mainly with satellite clock corrections. Based on these results and additional experiments, they state that the conclusions drawn by Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) are questionable. We do not agree with these conclusions and present our arguments in this article. In the first part, we review and highlight the main characteristics of the studies performed by Meindl et al. (Adv Space Res 51(7):1047–1064, 2013) to show that the experiments are quite different from those performed by Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026,2013) . In the second part, we show that normal equation (NEQ) systems are regular when estimating geocenter coordinates, implying that the covariance matrices associated with the NEQ systems may be used to assess the sensitivity to geocenter coordinates in a standard way. The sensitivity of GNSS to the components of the geocenter is discussed. Finally, we comment on the arguments raised by Rebischung et al. (J Geod doi:10.1016/j.asr.2012.10.026, 2013) against the results of Meindl et al. (Adv Space Res 51(7):1047–1064, 2013).

Galileo Orbit and Clock Quality of the IGS Multi-GNSS Experiment

The Multi-GNSS Experiment (MGEX) of the International GNSS Service (IGS) aims at the data collection and analysis of all available satellite navigation systems. In particular the new global and regional satellite navigation systems are of interest, i.e., the European Galileo, the Chinese BeiDou, the Japanese QZSS as well as satellite based augmentation systems. This article analyzes the orbit and clock quality of the Galileo products of four MGEX analysis centers for a common time period of 20 weeks. Orbit comparisons of the individual analysis centers have a consistency at the 5–30 cm level. Day boundary discontinuities range from 4 to 28 cm whereas 2-day orbit fit RMS values vary between 1 and 7 cm. The accuracy evaluated by satellite laser ranging residuals is on the one decimeter level with a systematic bias of about −5 cm for all analysis centers. In addition, systematic errors on the decimeter level related to solar radiation pressure mismodeling are present in all orbit products. Due to the correlation of radial orbit errors with the clock parameters, these errors are also visible as a bump in the Allan deviation of the Galileo satellite clocks at the orbital frequency.