Orbital debris mitigation through deorbiting with passive electrodynamic drag

The increase of orbital debris and the consequent proliferation of smaller objects through fragmentation are driving the need for mitigation strategies. The issue is how to deorbit the satellite with an efficient system that does not impair drastically the propellant budget of the satellite and, consequently, reduces its operating life. We have been investigating, in the framework of a European-Community-funded project, a passive system that makes use of an electrodynamics tether to deorbit a satellite through Lorentz forces. The deorbiting system will be carried by the satellite itself at launch and deployed from the satellite at the end of its life. From that moment onward the system operates passively without requiring any intervention from the satellite itself. The paper summarizes the results of the analysis carried out to show the deorbiting performance of the system starting from different orbital altitudes and inclinations for a reference satellite mass. Results can be easily scaled to other satellite masses. The results have been obtained by using a high-fidelity computer model that uses the latest environmental routines for magnetic field, ionospheric density, atmospheric density and a gravity field model. The tether dynamics is modelled by considering all the main aspects of a real system as the tether flexibility and its temperature-dependent electrical conductivity. Temperature variations are computed by including all the major external and internal input fluxes and the thermal flux emitted from the tether. The results shows that a relatively compact and light system can carry out the complete deorbit of a relatively large satellite in a time ranging from a month to less than a year starting from high LEO with the best performance occurring at low orbital inclinations.

The orbital-motion-limited regime of cylindrical Langmuir probes

An asymptotic analysis of electron collection at high bias Fp serves to determine the domain of validity of the orbital-motion-limited regime of cylindrical Langmuir probes, which is basic for the workings of conductive bare tethers. The radius of a wire collecting OML current in an unmagnetized plasma at rest cannot exceed a value, Rmax , which is found to exhibit a minimum as a function of Fp ; atFp values of interest, Rmax is already increasing and is larger than the electron Debye length lDe . The breakdown of the regime relates to conditions far fromthe probe, at electron energies comparable to the ion thermal energy, kTi ; Rmax is found to increase with Ti . It is also found that ~1! the maximumwidth of a thin tape, if used instead of a wire, is 4Rmax ; ~2! the electron thermal gyroradius must be larger than both R and lDe for magnetic effects to be negligible; and ~3! conditions applying to the tether case are such that trapped-orbit effects are negligible.

Cylindrical Langmuir probes beyond the orbital-motion-limited regime

The current I to a cylindrical probe at rest in an unmagnetized plasma, with probe bias highly positive, is determined. The way I lags behind the orbital-motion-limited OMLcurrent, 1 OML R, as the radius R exceeds the maximum radius for the OML regime to hold, is of interest for space-tether applications. The ratio I/I OML is roughly a decreasing function of R/lD R max /lDe , which is independent of bias, with lDe the electron Debye length and Rmax /l De roughly an increasing function of the temperature ratio, Ti /Te. The dependence of current on ion energy is used to discuss the effect of probe motion through the plasma, a case applying to tethers in low orbit.

Jupiter power generation with electrodynamic tethers at constant orbital energy

An electrodynamic tether system for power generation at Jupiter is presented that allows extracting energy from Jupiter's corotating plasmasphere while leaving the system orbital energy unaltered to first order. The spacecraft is placed in a polar orbit with the tether spinning in the orbital plane so that the resulting Lorentz force, neglecting Jupiter's magnetic dipole tilt, is orthogonal to the instantaneous velocity vector and orbital radius, hence affecting orbital inclination rather than orbital energy. In addition, the electrodynamic tether subsystem, which consists of two radial tether arms deployed from the main central spacecraft, is designed in such a way as to extract maximum power while keeping the resulting Lorentz torque constantly null. The power-generation performance of the system and the effect on the orbit inclination is evaluated analytically for different orbital conditions and verified numerically. Finally, a thruster-based inclination-compensation maneuver at apoapsis is added, resulting in an efficient scheme to extract energy from the plasmasphere of the planet with minimum propellant consumption and no inclination change. A tradeoff analysis is conducted showing that, depending on tether size and orbit characteristics, the system performance can be considerably higher than conventional power-generation methods.

Analysis of tape tether survival in LEO against orbital debris

The low earth orbit (LEO) environment contains a large number of artificial debris, of which a significant portion is due to dead satellites and fragments of satellites resulted from explosions and in-orbit collisions. Deorbiting defunct satellites at the end of their life can be achieved by a successful operation of an Electrodynamic Tether (EDT) system. The effectiveness of an EDT greatly depends on the survivability of the tether, which can become debris itself if cut by debris particles; a tether can be completely cut by debris having some minimal diameter. The objective of this paper is to develop an accurate model using power laws for debris-size ranges, in both ORDEM2000 and MASTER2009 debris flux models, to calculate tape tether survivability. The analytical model, which depends on tape dimensions (width, thickness) and orbital parameters (inclinations, altitudes) is then verified with fully numerical results to compare for different orbit inclinations, altitudes and tape width for both ORDEM2000 and MASTER2009 flux data.

A Special Perturbation Method in Orbital Dynamics

Bead models are used in dynamical simulation of tethers. These models discretize a cable using beads distributed along its length. The time evolution is obtained nu- merically. Typically the number of particles ranges between 5 and 50, depending on the required accuracy. Sometimes the simulation is extended over long periods (several years). The complex interactions between the cable and its spatial environment require to optimize the propagators —both in runtime and precisión that constitute the central core of the process. The special perturbation method treated on this article conjugates simpleness of computer implementation, speediness and precision, and is capable to propagate the orbit of whichever material particle. The paper describes the evolution of some orbital elements, which are constants in a non-perturbed problem, but which evolve in the time scale imposed by the perturbation. It can be used with any kind of orbit and it is free of sin- gularities related to small inclination and/or small eccentricity. The use of Euler parameters makes it robust.

Low Work-Function Thermionic Emission and Orbital-Motion-Limited Ion Collection at Bare-Tether Cathodic Contact

Una amarra electrodinámica (electrodynamic tether) opera sobre principios electromagnéticos intercambiando momento con la magnetosfera planetaria e interactuando con su ionosfera. Es un subsistema pasivo fiable para desorbitar etapas de cohetes agotadas y satélites al final de su misión, mitigando el crecimiento de la basura espacial. Una amarra sin aislamiento captura electrones del plasma ambiente a lo largo de su segmento polarizado positivamente, el cual puede alcanzar varios kilómetros de longitud, mientras que emite electrones de vuelta al plasma mediante un contactor de plasma activo de baja impedancia en su extremo catódico, tal como un cátodo hueco (hollow cathode). En ausencia de un contactor catódico activo, la corriente que circula por una amarra desnuda en órbita es nula en ambos extremos de la amarra y se dice que ésta está flotando eléctricamente. Para emisión termoiónica despreciable y captura de corriente en condiciones limitadas por movimiento orbital (orbital-motion-limited, OML), el cociente entre las longitudes de los segmentos anódico y catódico es muy pequeño debido a la disparidad de masas entre iones y electrones. Tal modo de operación resulta en una corriente media y fuerza de Lorentz bajas en la amarra, la cual es poco eficiente como dispositivo para desorbitar. El electride C12A7 : e−, que podría presentar una función de trabajo (work function) tan baja como W = 0.6 eV y un comportamiento estable a temperaturas relativamente altas, ha sido propuesto como recubrimiento para amarras desnudas. La emisión termoiónica a lo largo de un segmento así recubierto y bajo el calentamiento de la operación espacial, puede ser más eficiente que la captura iónica. En el modo más simple de fuerza de frenado, podría eliminar la necesidad de un contactor catódico activo y su correspondientes requisitos de alimentación de gas y subsistema de potencia, lo que resultaría en un sistema real de amarra “sin combustible”. Con este recubrimiento de bajo W, cada segmento elemental del segmento catódico de una amarra desnuda de kilómetros de longitud emitiría corriente como si fuese parte de una sonda cilíndrica, caliente y uniformemente polarizada al potencial local de la amarra. La operación es similar a la de una sonda de Langmuir 2D tanto en los segmentos catódico como anódico. Sin embargo, en presencia de emisión, los electrones emitidos resultan en carga espacial (space charge) negativa, la cual reduce el campo eléctrico que los acelera hacia fuera, o incluso puede desacelerarlos y hacerlos volver a la sonda. Se forma una doble vainas (double sheath) estable con electrones emitidos desde la sonda e iones provenientes del plasma ambiente. La densidad de corriente termoiónica, variando a lo largo del segmento catódico, podría seguir dos leyes distintas bajo diferentes condiciones: (i) la ley de corriente limitada por la carga espacial (space-charge-limited, SCL) o (ii) la ley de Richardson-Dushman (RDS). Se presenta un estudio preliminar sobre la corriente SCL frente a una sonda emisora usando la teoría de vainas (sheath) formada por la captura iónica en condiciones OML, y la corriente electrónica SCL entre los electrodos cilíndricos según Langmuir. El modelo, que incluye efectos óhmicos y el efecto de transición de emisión SCL a emisión RDS, proporciona los perfiles de corriente y potencial a lo largo de la longitud completa de la amarra. El análisis muestra que en el modo más simple de fuerza de frenado, bajo condiciones orbitales y de amarras típicas, la emisión termoiónica proporciona un contacto catódico eficiente y resulta en una sección catódica pequeña. En el análisis anterior, tanto la transición de emisión SCL a RD como la propia ley de emisión SCL consiste en un modelo muy simplificado. Por ello, a continuación se ha estudiado con detalle la solución de vaina estacionaria de una sonda con emisión termoiónica polarizada negativamente respecto a un plasma isotrópico, no colisional y sin campo magnético. La existencia de posibles partículas atrapadas ha sido ignorada y el estudio incluye tanto un estudio semi-analítico mediante técnica asintóticas como soluciones numéricas completas del problema. Bajo las tres condiciones (i) alto potencial, (ii) R = Rmax para la validez de la captura iónica OML, y (iii) potencial monotónico, se desarrolla un análisis asintótico auto-consistente para la estructura de plasma compleja que contiene las tres especies de cargas (electrones e iones del plasma, electrones emitidos), y cuatro regiones espaciales distintas, utilizando teorías de movimiento orbital y modelos cinéticos de las especies. Aunque los electrones emitidos presentan carga espacial despreciable muy lejos de la sonda, su efecto no se puede despreciar en el análisis global de la estructura de la vaina y de dos capas finas entre la vaina y la región cuasi-neutra. El análisis proporciona las condiciones paramétricas para que la corriente sea SCL. También muestra que la emisión termoiónica aumenta el radio máximo de la sonda para operar dentro del régimen OML y que la emisión de electrones es mucho más eficiente que la captura iónica para el segmento catódico de la amarra. En el código numérico, los movimientos orbitales de las tres especies son modelados para potenciales tanto monotónico como no-monotónico, y sonda de radio R arbitrario (dentro o más allá del régimen de OML para la captura iónica). Aprovechando la existencia de dos invariante, el sistema de ecuaciones Poisson-Vlasov se escribe como una ecuación integro-diferencial, la cual se discretiza mediante un método de diferencias finitas. El sistema de ecuaciones algebraicas no lineal resultante se ha resuelto de con un método Newton-Raphson paralelizado. Los resultados, comparados satisfactoriamente con el análisis analítico, proporcionan la emisión de corriente y la estructura del plasma y del potencial electrostático. ABSTRACT An electrodynamic tether operates on electromagnetic principles and exchanges momentum through the planetary magnetosphere, by continuously interacting with the ionosphere. It is a reliable passive subsystem to deorbit spent rocket stages and satellites at its end of mission, mitigating the growth of orbital debris. A tether left bare of insulation collects electrons by its own uninsulated and positively biased segment with kilometer range, while electrons are emitted by a low-impedance active device at the cathodic end, such as a hollow cathode, to emit the full electron current. In the absence of an active cathodic device, the current flowing along an orbiting bare tether vanishes at both ends and the tether is said to be electrically floating. For negligible thermionic emission and orbital-motion-limited (OML) collection throughout the entire tether (electron/ion collection at anodic/cathodic segment, respectively), the anodic-to-cathodic length ratio is very small due to ions being much heavier, which results in low average current and Lorentz drag. The electride C12A7 : e−, which might present a possible work function as low as W = 0.6 eV and moderately high temperature stability, has been proposed as coating for floating bare tethers. Thermionic emission along a thus coated cathodic segment, under heating in space operation, can be more efficient than ion collection and, in the simplest drag mode, may eliminate the need for an active cathodic device and its corresponding gas-feed requirements and power subsystem, which would result in a truly “propellant-less” tether system. With this low-W coating, each elemental segment on the cathodic segment of a kilometers-long floating bare-tether would emit current as if it were part of a hot cylindrical probe uniformly polarized at the local tether bias, under 2D probe conditions that are also applied to the anodic-segment analysis. In the presence of emission, emitted electrons result in negative space charge, which decreases the electric field that accelerates them outwards, or even reverses it, decelerating electrons near the emitting probe. A double sheath would be established with electrons being emitted from the probe and ions coming from the ambient plasma. The thermionic current density, varying along the cathodic segment, might follow two distinct laws under different con ditions: i) space-charge-limited (SCL) emission or ii) full Richardson-Dushman (RDS) emission. A preliminary study on the SCL current in front of an emissive probe is presented using the orbital-motion-limited (OML) ion-collection sheath and Langmuir’s SCL electron current between cylindrical electrodes. A detailed calculation of current and bias profiles along the entire tether length is carried out with ohmic effects considered and the transition from SCL to full RDS emission is included. Analysis shows that in the simplest drag mode, under typical orbital and tether conditions, thermionic emission provides efficient cathodic contact and leads to a short cathodic section. In the previous analysis, both the transition between SCL and RDS emission and the current law for SCL condition have used a very simple model. To continue, considering an isotropic, unmagnetized, colissionless plasma and a stationary sheath, the probe-plasma contact is studied in detail for a negatively biased probe with thermionic emission. The possible trapped particles are ignored and this study includes both semianalytical solutions using asymptotic analysis and complete numerical solutions. Under conditions of i) high bias, ii) R = Rmax for ion OML collection validity, and iii) monotonic potential, a self-consistent asymptotic analysis is carried out for the complex plasma structure involving all three charge species (plasma electrons and ions, and emitted electrons) and four distinct spatial regions using orbital motion theories and kinetic modeling of the species. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between the sheath and the quasineutral region. The parametric conditions for the current to be space-chargelimited are obtained. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers. In the numerical code, the orbital motions of all three species are modeled for both monotonic and non-monotonic potential, and for any probe radius R (within or beyond OML regime for ion collection). Taking advantage of two constants of motion (energy and angular momentum), the Poisson-Vlasov equation is described by an integro differential equation, which is discretized using finite difference method. The non-linear algebraic equations are solved using a parallel implementation of the Newton-Raphson method. The results, which show good agreement with the analytical results, provide the results for thermionic current, the sheath structure, and the electrostatic potential.

Low work-function thermionic emission and orbital-motion-limited ion collection at bare-tether cathodic contact

With a thin coating of low-work-function material, thermionic emission in the cathodic segment of bare tethers might be much greater than orbital-motion-limited (OML) ion collection current. The space charge of the emitted electrons decreases the electric field that accelerates them outwards, and could even reverse it for high enough emission, producing a potential hollow. In this work, at the conditions of high bias and relatively low emission that make the potential monotonic, an asymptotic analysis is carried out, extending the OML ion-collection analysis to investigate the probe response due to electrons emitted by the negatively biased cylindrical probe. At given emission, the space charge effect from emitted electrons increases with decreasing magnitude of negative probe bias. Although emitted electrons present negligible space charge far away from the probe, their effect cannot be neglected in the global analysis for the sheath structure and two thin layers in between sheath and the quasineutral region. The space-charge-limited condition is located. It is found that thermionic emission increases the range of probe radius for OML validity and is greatly more effective than ion collection for cathodic contact of tethers.

Spectrum of 100-kyr glacial cycle: Orbital inclination, not eccentricity

Spectral analysis of climate data shows a strong narrow peak with period ≈100 kyr, attributed by the Milankovitch theory to changes in the eccentricity of the earth’s orbit. The narrowness of the peak does suggest an astronomical origin; however the shape of the peak is incompatible with both linear and nonlinear models that attribute the cycle to eccentricity or (equivalently) to the envelope of the precession. In contrast, the orbital inclination parameter gives a good match to both the spectrum and bispectrum of the climate data. Extraterrestrial accretion from meteoroids or interplanetary dust is proposed as a mechanism that could link inclination to climate, and experimental tests are described that could prove or disprove this hypothesis.

A long marine history of carbon cycle modulation by orbital-climatic changes

Pacing of the marine carbon cycle by orbital forcing during the Pliocene and Pleistocene Ice Ages [past 2.5 million years (Myr)] is well known. As older deep-sea sediment records are being studied at greater temporal resolution, it is becoming clear that similar fluctuations in the marine carbon system have occurred throughout the late Mesozoic and Tertiary, despite the absence of large continental ice sheets over much of this time. Variations in both the organic and the calcium carbonate components of the marine carbon system seem to have varied cyclically in response to climate forcing, and carbon and carbonate time series appear to accurately characterize the frequency spectrum of ancient climatic change. For the past 35 Myr, much of the variance in carbonate content carries the “polar” signal of obliquity [41,000 years (41 kyr)] forcing. Over the past 125 Myr, there is evidence from marine sediments of the continued role of precessional (≈21 kyr) climatic cycles. Repeat patterns of sedimentation at about 100, 400, and 2,400 kyr, the modulation periods of precession, persistently enter into marine carbon cycle records as well. These patterns suggest a nonlinear response of climate and/or the sedimentation of organic carbon and carbonates to precessional orbital perturbations. Nonlinear responses of the carbon system may help to amplify relatively weak orbital insolation anomalies into more significant climatic perturbations through positive feedback effects. Nonlinearities in the carbon cycle may have transformed orbital-climatic cycles into long-wavelength features on time scales comparable to the residence times of carbon and nutrient elements in the ocean.

Propriedades magnéticas de trímeros de FexCo1-x depositados em Pt(111)

Quais propriedades magnéticas são modificadas quando se agrupam átomos de Fe/Co para formar estruturas quasi-2D, se comparadas aos nanofios (quasi-1D) de FexCo1-x? E como estas propriedades reagem com a variação da proporção de Fe/Co nos aglomerados? A fim de responder a estas questões, trímeros de FexCo1-x depositados em Pt(111) são investigados utilizando o método de primeiros princípios Real Space-Linear Muffin-Tin Orbital-Atomic Sphere Approximation (RS-LMTO-ASA) no âmbito da Teoria do Funcional da Densidade (DFT). Diferentes configurações de trímeros triangulares são consideradas, variando-se as posições e a concentração dos átomos de Fe/Co. Neste trabalho, demonstra-se a ocorrência de uma tendência não-linear estritamente decrescente dos momentos orbitais médios como função da concentração de Fe, distinta do encontrado tanto para os nanofios de FexCo1-x (dependência linear) quanto para a monocamada correspondente (dependência não-linear). Os resultados obtidos mostram ainda que os momentos orbitais variam com o ambiente local e com a direção de magnetização, especialmente quando associados aos átomos de Co, em concordância com publicações anteriores. A mudança de dimensionalidade quasi-1D (nanofios) para quasi-2D (trímeros compactos) não afeta o comportamento dos momentos de spin, que permanecem descritos por uma função linear com respeito à proporção de Fe/Co. Ambos o formato e a concentração de Fe nos sistemas apresentam um papel importante nos valores de energia de anisotropia magnética. Em adição, observou-se que o subtrato de Pt opera ativamente na definição das propriedades magnéticas dos aglomerados. Embora todas as configurações lineares e compactas dos aglomerados de FexCo1-x sejam estáveis e exibam interações fortemente ferromagnéticas entre os primeiros vizinhos, nem todas revelaram o ordenamento colinear como estado fundamental, apresentando uma interação de Dzyaloshinskii-Moriya não-desprezível induzida pelo acoplamento spin-órbita. Estes casos específicos são: o trímero triangular de Co puro e o trímero linear (nanofio) de Fe puro, para o qual foi verificado o acoplamento do tipo Ruderman-Kittel-Kasuya-Yosida entre os átomos de Fe constituintes. Os resultados obtidos contribuem para o entendimento de quais mecanismos definem o magnetismo nos trímeros de FexCo1-x/Pt(111), e discutem as questões presentes atualmente na literatura no contexto destes sistemas.

Magnetic and orbital blocking in Ni nanocontacts

We address the fundamental question of whether magnetoresistance (MR) of atomic-sized contacts of nickel is very large because of the formation of a domain wall (DW) at the neck. Using ab initio transport calculations we find that, as in the case of nonmagnetic electrodes, transport in Ni nanocontacts depends very much on the orbital nature of the electrons. Our results are in agreement with several experiments in the average value of the conductance. On the other hand, contrary to existing claims, DW scattering does not account for large MR in Ni nanocontacts.

Vacancy-induced magnetism in graphene and graphene ribbons

We address the electronic structure and magnetic properties of vacancies and voids both in graphene and graphene ribbons. By using a mean-field Hubbard model, we study the appearance of magnetic textures associated with removing a single atom (vacancy) and multiple adjacent atoms (voids) as well as the magnetic interactions between them. A simple set of rules, based on the Lieb theorem, link the atomic structure and the spatial arrangement of the defects to the emerging magnetic order. The total spin S of a given defect depends on its sublattice imbalance, but some defects with S=0 can still have local magnetic moments. The sublattice imbalance also determines whether the defects interact ferromagnetically or antiferromagnetically with one another and the range of these magnetic interactions is studied in some simple cases. We find that in semiconducting armchair ribbons and two-dimensional graphene without global sublattice imbalance, there is a maximum defect density above which local magnetization disappears. Interestingly, the electronic properties of semiconducting graphene ribbons with uncoupled local moments are very similar to those of diluted magnetic semiconductors, presenting giant Zeeman splitting.

Astrophysical parameters and orbital solution of the peculiar X-ray transient IGR J00370+6122

Context. BD + 60° 73 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. Aims. We aim to characterise the binary system IGR J00370+6122 by deriving its orbital and physical parameters. Methods. We obtained high-resolution spectra of BD + 60° 73 at different epochs. We used the fastwind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. Results. BD + 60° 73 is a BN0.7 Ib low-luminosity supergiant located at a distance ~3.1 kpc, in the Cas OB4 association. We derive Teff = 24 000 K and log gc = 3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1-σ level with a mass M∗ ≈ 15 M⊙. The recurrence time of the X-ray flares is the orbital period of the system. The neutron star is in a high-eccentricity (e = 0.56 ± 0.07) orbit, and the X-ray emission is strongly peaked around orbital phase φ = 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. Conclusions. The X-ray behaviour of IGR J00370+6122 is reminiscent of “intermediate” supergiant X-ray transients, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which when combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will very likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.