Analytical models for the power subsystem and the attitude control subsystem of a microsatellite

El trabajo realizado en la presente tesis doctoral se debe considerar parte del proyecto UPMSat-2, que se enmarca dentro del ámbito de la tecnología aeroespacial. El UPMSat-2 es un microsatélite (de bajo coste y pequeño tamaño) diseñado, construido, probado e integrado por la Universidad Politécnica de Madrid (España), para fines de demostración tecnológica y educación. El objetivo de la presente tesis doctoral es presentar nuevos modelos analíticos para estudiar la interdependencia energética entre los subsistemas de potencia y de control de actitud de un satélite. En primer lugar, se estudia la simulación del subsistema de potencia de un microsatélite, prestando especial atención a la simulación de la fuente de potencia, esto es, los paneles solares. En la tesis se presentan métodos sencillos pero precisos para simular la producción de energía de los paneles en condiciones ambientales variables a través de su circuito equivalente. Los métodos propuestos para el cálculo de los parámetros del circuito equivalente son explícitos (o al menos, con las variables desacopladas), no iterativos y directos; no se necesitan iteraciones o valores iniciales para calcular los parámetros. La precisión de este método se prueba y se compara con métodos similares de la literatura disponible, demostrando una precisión similar para mayor simplicidad. En segundo lugar, se presenta la simulación del subsistema de control de actitud de un microsatélite, prestando especial atención a la nueva ley de control propuesta. La tesis presenta un nuevo tipo de control magnético es aplicable a la órbita baja terrestre (LEO). La ley de control propuesta es capaz de ajustar la velocidad de rotación del satélite alrededor de su eje principal de inercia máximo o mínimo. Además, en el caso de órbitas de alta inclinación, la ley de control favorece la alineación del eje de rotación con la dirección normal al plano orbital. El algoritmo de control propuesto es simple, sólo se requieren magnetopares como actuadores; sólo se requieren magnetómetros como sensores; no hace falta estimar la velocidad angular; no incluye un modelo de campo magnético de la Tierra; no tiene por qué ser externamente activado con información sobre las características orbitales y permite el rearme automático después de un apagado total del subsistema de control de actitud. La viabilidad teórica de la citada ley de control se demuestra a través de análisis de Monte Carlo. Por último, en términos de producción de energía, se demuestra que la actitud propuesto (en eje principal perpendicular al plano de la órbita, y el satélite que gira alrededor de ella con una velocidad controlada) es muy adecuado para la misión UPMSat-2, ya que permite una área superior de los paneles apuntando hacia el sol cuando se compara con otras actitudes estudiadas. En comparación con el control de actitud anterior propuesto para el UPMSat-2 resulta en un incremento de 25% en la potencia disponible. Además, la actitud propuesto mostró mejoras significativas, en comparación con otros, en términos de control térmico, como la tasa de rotación angular por satélite puede seleccionarse para conseguir una homogeneización de la temperatura más alta que apunta satélite y la antena. ABSTRACT The work carried out in the present doctoral dissertation should be considered part of the UPMSat-2 project, falling within the scope of the aerospace technology. The UPMSat-2 is a microsatellite (low cost and small size) designed, constructed integrated and tested for educational and technology demonstration purposes at the Universidad Politécnica de Madrid (Spain). The aim of the present doctoral dissertation is to present new analytical models to study the energy interdependence between the power and the attitude control subsystems of a satellite. First, the simulation of the power subsystem of a microsatellite is studied, paying particular attention to the simulation of the power supply, i.e. the solar panels. Simple but accurate methods for simulate the power production under variable ambient conditions using its equivalent circuit are presented. The proposed methods for calculate the equivalent circuit parameters are explicit (or at least, with decoupled variables), non-iterative and straight forward; no iterations or initial values for the parameters are needed. The accuracy of this method is tested and compared with similar methods from the available literature demonstrating similar precision but higher simplicity. Second, the simulation of the control subsystem of a microsatellite is presented, paying particular attention to the new control law proposed. A new type of magnetic control applied to Low Earth Orbit (LEO) satellites has been presented. The proposed control law is able to set the satellite rotation speed around its maximum or minimum inertia principal axis. Besides, the proposed control law favors the alignment of this axis with the normal direction to the orbital plane for high inclination orbits. The proposed control algorithm is simples, only magnetorquers are required as actuators; only magnetometers are required as sensors; no estimation of the angular velocity is needed; it does not include an in-orbit Earth magnetic field model; it does not need to be externally activated with information about the orbital characteristics and it allows automatic reset after a total shutdown of attitude control subsystem. The theoretical viability of the control law is demonstrated through Monte Carlo analysis. Finally, in terms of power production, it is demonstrated that the proposed attitude (on principal axis perpendicular to the orbit plane, and the satellite rotating around it with a controlled rate) is quite suitable for the UPMSat-2 mission, as it allows a higher area of the panels pointing towards the sun when compared to other studied attitudes. Compared with the previous attitude control proposed for the UPMSat-2 it results in a 25% increment in available power. Besides, the proposed attitude showed significant improvements, when compared to others, in terms of thermal control, as the satellite angular rotation rate can be selected to achieve a higher temperature homogenization of the satellite and antenna pointing.

Outreach Testing of Ancient Astronomy

This work is an outreach approach to an ubiquitous recent problem in secondary-school education: how to face back the decreasing interest in natural sciences shown by students under ‘pressure’ of convenient resources in digital devices/applications. The approach rests on two features. First, empowering of teen-age students to understand regular natural events around, as very few educated people they meet could do. Secondly, an understanding that rests on personal capability to test and verify experimental results from the oldest science, astronomy, with simple instruments as used from antiquity down to the Renaissance (a capability restricted to just solar and lunar motions). Because lengths in astronomy and daily life are so disparate, astronomy basically involved observing and registering values of angles (along with times), measurements being of two types, of angles on the ground and of angles in space, from the ground. First, the gnomon, a simple vertical stick introduced in Babylonia and Egypt, and then in Greece, is used to understand solar motion. The gnomon shadow turns around during any given day, varying in length and thus angle between solar ray and vertical as it turns, going through a minimum (noon time, at a meridian direction) while sweeping some angular range from sunrise to sunset. Further, the shadow minimum length varies through the year, with times when shortest and sun closest to vertical, at summer solstice, and times when longest, at winter solstice six months later. The extreme directions at sunset and sunrise correspond to the solstices, swept angular range greatest at summer, over 180 degrees, and the opposite at winter, with less daytime hours; in between, spring and fall equinoxes occur, marked by collinear shadow directions at sunrise and sunset. The gnomon allows students to determine, in addition to latitude (about 40.4° North at Madrid, say), the inclination of earth equator to plane of its orbit around the sun (ecliptic), this fundamental quantity being given by half the difference between solar distances to vertical at winter and summer solstices, with value about 23.5°. Day and year periods greatly differing by about 2 ½ orders of magnitude, 1 day against 365 days, helps students to correctly visualize and interpret the experimental measurements. Since the gnomon serves to observe at night the moon shadow too, students can also determine the inclination of the lunar orbital plane, as about 5 degrees away from the ecliptic, thus explaining why eclipses are infrequent. Independently, earth taking longer between spring and fall equinoxes than from fall to spring (the solar anomaly), as again verified by the students, was explained in ancient Greek science, which posited orbits universally as circles or their combination, by introducing the eccentric circle, with earth placed some distance away from the orbital centre when considering the relative motion of the sun, which would be closer to the earth in winter. In a sense, this can be seen as hint and approximation of the elliptic orbit proposed by Kepler many centuries later.

Estimation of satellite antenna phase center offsets for Galileo

Satellite antenna phase center offsets for the GalileoInOrbitValidation(IOV) and FullOperationalCapability (FOC) satellites are estimated by two different analysiscenters based on tracking data of a global GNSS network. The mean x- and y-offsets could be determined with a precision of a few centimeters. However, daily estimates of thex-offsets of the IOV satellites show pronounced systematic effects with a peak-to-peak amplitude of up to 70 cm that depend on the orbit model and the elevation of the Sun above the orbital plane. For the IOV y-offsets, no dependence on the orbit model exists but the scatter strongly depends on the elevation of the Sun above the orbital plane. In general, these systematic effects are significantly smaller for the FOC satellites. The z-offsets of the two analysis centers agree within the 10–15 cm level, and the time series do not show systematic effects. The application of an averaged Galileo satellite antenna model obtained from the two solutions results in a reduction of orbit day boundary discontinuities by up to one third—even if an independent software package is used.

Type Ia supernovae from violent mergers of carbon-oxygen white dwarfs: polarization signatures

The violent merger of two carbon-oxygen white dwarfs has been proposed as a viable progenitor for some Type Ia supernovae. However, it has been argued that the strong ejecta asymmetries produced by this model might be inconsistent with the low degree of polarization typically observed in Type Ia supernova explosions. Here, we test this claim by carrying out a spectropolarimetric analysis for the model proposed by Pakmor et al. for an explosion triggered during the merger of a 1.1 and 0.9 M⊙ carbon-oxygen white dwarf binary system. Owing to the asymmetries of the ejecta, the polarization signal varies significantly with viewing angle. We find that polarization levels for observers in the equatorial plane are modest (≲1 per cent) and show clear evidence for a dominant axis, as a consequence of the ejecta symmetry about the orbital plane. In contrast, orientations out of the plane are associated with higher degrees of polarization and departures from a dominant axis. While the particular model studied here gives a good match to highly polarized events such as SN 2004dt, it has difficulties in reproducing the low polarization levels commonly observed in normal Type Ia supernovae. Specifically, we find that significant asymmetries in the element distribution result in a wealth of strong polarization features that are not observed in the majority of currently available spectropolarimetric data of Type Ia supernovae. Future studies will map out the parameter space of the merger scenario to investigate if alternative models can provide better agreement with observations.

Magnetic structure and orbital ordering in tetragonal and monoclinic KCrF3 from first-principles calculations

KCrF3 has been systematically investigated by using the full-potential linearized augmented plane wave plus local orbital method within the generalized gradient approximation and the local spin density approximation plus the on-site Coulomb repulsion approach. The total energies for ferromagnetic and three different antiferromagnetic configurations are calculated in the high-temperature tetragonal and low-temperature monoclinic phases, respectively.

Charge, orbital, and magnetic ordering in YBaFe2O5 from first-principles calculations

First principles calculations using the augmented plane wave plus local orbitals method, as implemented in the WIEN2k code, have been used to investigate the electronic and magnetic properties of YBaFe2O5, especially as regards the charge-orbital ordering. Although the total 3d charge disproportion is rather small, an orbital order parameter defined as the difference between t(2g) orbital occupations of Fe2+ and Fe3+ cations is large (0,73) and gives unambiguous evidence for charge and orbital ordering: Strong hybridization between O 2p and Fe e(g) states results in the nearly complete loss of the separation between the total charges at the Fe2+ and Fe3+ atoms.

Orbital ordering in Cs2AgF4 from first principles

First-principles calculations using the augmented plane wave plus local orbital method, as implemented in the WIEN2K code, have been used to investigate the structural, electronic, and magnetic properties of the layered perovskite Cs2AgF4. Our calculations indicate that an orthorhombic ground state for Cs2AgF4 is energetically favored over tetragonal. We also find that Cs2AgF4 should be a strong two-dimensional ferromagnet, with very weak antiferromagnetic coupling between the layers, in agreement with the experiment. More importantly, an antiferrodistortive ordering of z(2)-x(2) and z(2)-y(2) orbitals is inferred from the density of states and from a spin density isosurface analysis.

Two-dimensional electron states bound to an off-plane donor in a magnetic field

The states of an electron confined in a two-dimensional (2D) plane and bound to an off-plane donor impurity center, in the presence of a magnetic field, are investigated. The energy levels of the ground state and the first three excited states are calculated variationally. The binding energy and the mean orbital radius of these states are obtained as a function of the donor center position and the magnetic field strength. The limiting cases are discussed for an in-plane donor impurity (i.e. a 2D hydrogen atom) as well as for the donor center far away from the 2D plane in strong magnetic fields, which corresponds to a 2D harmonic oscillator.

DETERMINATION of THE INDIRECT ORIENTATION of ORBITAL PUSHBROOM IMAGES USING CONTROL STRAIGHT LINES

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

Two-dimensional electron states bound to an off-plane donor in a magnetic field

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

An analytical and numerical study of plane change maneuvers using aerodynamic force

The goal of the present work is to analyze space missions that use the terrestrial atmosphere to accomplish orbital maneuvers that involve a plane change. A set of analytical solutions is presented for the variation of the orbital elements due to a single passage through the atmosphere, assuming that the interval the spacecraft travels through the atmosphere is not too large. The study considers both the lift influence on the spacecraft orbit as well as drag. The final equations are tested with numerical integration and can be considered in accordance with the numerical results whenever the perigee height is larger than a critical value. Next, a numerical study of the ratio between the velocity increment required to correct the semimajor axis decay due to the atmospheric passage and the velocity variation required to obtain the change in the inclination is also presented. This analysis can be used to decide if a maneuver passing through the atmosphere can decrease the fuel consumption of the mission and, in the cases where this technique can be used, if a multiple passage is more efficient than a single passage.

An analysis of the Hygiea asteroid family orbital region

(10) Hygiea is the fourth largest asteroid of the main belt, by volume and mass, and it is the largest member of its family, that is made mostly by low-albedo, C-type asteroids, typical of the outer main belt. Like many other large families, it is associated with a 'halo' of objects, that extends far beyond the boundary of the core family, as detected by traditional hierarchical clustering methods (HCM) in proper element domains. Numerical simulations of the orbital evolution of family members may help in estimating the family and halo family age, and the original ejection velocity field. But, in order to minimize the errors associated with including too many interlopers, it is important to have good estimates of family membership that include available data on local asteroid taxonomy, geometrical albedo and local dynamics. For this purpose, we obtained synthetic proper elements and frequencies of asteroids in the Hygiea orbital region, with their errors. We revised the current knowledge on asteroid taxonomy, including Sloan Digital Sky Survey-Moving Object Catalog 4th release (SDSS-MOC 4) data, and geometric albedo data from Wide-field Infrared Survey Explorer (WISE) and Near-Earth Object WISE (NEOWISE). We identified asteroid family members using HCM in the domain of proper elements (a, e, sin (i)) and in the domains of proper frequencies most appropriate to study diffusion in the local web of secular resonances, and eliminated possible interlopers based on taxonomic and geometrical albedo considerations. To identify the family halo, we devised a new hierarchical clustering method in an extended domain that includes proper elements, principal components PC1, PC2 obtained based on SDSS photometric data and, for the first time, WISE and NEOWISE geometric albedo. Data on asteroid size distribution, light curves and rotations were also revised for the Hygiea family. The Hygiea family is the largest group in its region, with two smaller families in proper element domain and 18 families in various frequencies domains identified in this work for the first time. Frequency groups tend to extend vertically in the (a, sin (i)) plane and cross not only the Hygiea family but also the near C-type families of Themis and Veritas, causing a mixture of objects all of relatively low albedo in the Hygiea family area. A few high-albedo asteroids, most likely associated with the Eos family, are also present in the region. Finally, the new multidomains hierarchical clustering method allowed us to obtain a good and robust estimate of the membership of the Hygiea family halo, quite separated from other asteroids families halo in the region, and with a very limited (about 3 per cent) presence of likely interlopers. © 2013 The Author Published by Oxford University Press on behalf of the Royal Astronomical Society.

Controle orbital de satélites artificias com propulsão e uso de gravidade lunar

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

Merging of low-mass systems and the origin of the Fundamental Plane

We present a new set of dissipationless N-body simulations to examine the feasibility of creating bright ellipticals (following the Kormendy relation, hereafter KR) by hierarchically merging present-day early-type dwarf galaxies, and to study how the encounter parameters affect the location of the end product in the (mu(e))-R-e plane. We investigate the merging of one-component galaxies of both equal and different masses, the merging of two-component galaxy models to explore the effect of dark haloes on the final galaxy characteristics, and the merging of ultracompact dwarf galaxies. We find that the increase of (mu(e)) with R-e is attributable to an increase in the initial orbital energy. The merger remnants shift down in the (mu(e))-R-e plane and fail to reach the KR. Thus, the KR is not reproducible by mergers of dwarf early-type systems, rendering untenable the theory that present-day dwarfs are responsible for even a small fraction of the present-day ellipticals, unless a considerable amount of dissipation is invoked. However, we do find that present-day dwarfs can be formed by the merger of ultracompact dwarfs.

Why H Atom Prefers the On-Top Site and Alkali Metals Favor the Middle Hollow Site on the Basal Plane of Graphite

In this work, the different adsorption properties of H and alkali metal atoms on the basal plane of graphite are studied and compared using a density functional method on the same model chemistry level. The results show that H prefers the on-top site while alkali metals favor the middle hollow site of graphite basal plane due to the unique electronic structures of H, alkali metals, and graphite. H has a higher electronegativity than carbon, preferring to form a covalent bond with C atoms, whereas alkaline metals have lower electronegativity, tending to adsorb on the highest electrostatic potential sites. During adsorption, there are more charges transferred from alkali metal to graphite than from H to graphite.