Three decades of Space Activities at the Instituto Universitario de Microgravedad "Ignacio Da Riva" / Universidad Politécnica de Madrid (IDR/UPM) en History of Rocketry and Astronautics

Strictly speaking, space related activities at the Escuela Técnica Superior de Ingenieros Aeronáuticos (ETSIA) begun in 1973, when Prof. Ignacio Da Riva got a contract from the European Space Agency (ESA) to compile a handbook on spacecraft thermal control. By the same time, ESA issued an announcement of opportunities offering to the European scientific community the possibility of perform microgravity relevant experiments on board space platform like the European orbital laboratory Spacelab. Prof. Da Riva proposed one of the few selected experiments dealing with fluid physics under microgravity conditions, later flown on Spacelab-1 mission in 1983. These two events were the starting point where Prof. Da Riva, full professor of Aerodynamics at ETSIA, nucleated a small group of young professors and students located at the Laboratorio de Aerodinámica y Mecánica de Fluidos (LAMF) of ETSIA. Such group was leaded by Prof. Da Riva since its creation till 1991, when Prof. Da Riva died, and it was the seed of the more recently created research institute for aerospace science and technology named "Ignacio Da Riva" (IDR) in his honour. In this communication space related activities performed either at LAMF or IDR during the last three decades are briefly described.

Asymptotic Description of Flutter Amplitude Saturation by Nonlinear Friction Forces

The computation of the non-linear vibration dynamics of an aerodynamically unstable bladed-disk is a formidable numerical task, even for the simplified case of aerodynamic forces assumed to be linear. The nonlinear friction forces effectively couple dif- ferent travelling waves modes and, in order to properly elucidate the dynamics of the system, large time simulations are typically required to reach a final, saturated state. Despite of all the above complications, the output of the system (in the friction microslip regime) is basically a superposition of the linear aeroelastic un- stable travelling waves, which exhibit a slow time modulation that is much longer than the elastic oscillation period. This slow time modulation is due to both, the small aerodynamic effects and the small nonlinear friction forces, and it is crucial to deter- mine the final amplitude of the flutter vibration. In this presenta- tion we apply asymptotic techniques to obtain a new simplified model that captures the slow time dynamics of the amplitudes of the travelling waves. The resulting asymptotic model is very re- duced and extremely cheap to simulate, and it has the advantage that it gives precise information about the characteristics of the nonlinear friction models that actually play a role in the satura- tion of the vibration amplitude.

O Botafumeiro: parametric pumping in the Middle Age

The pendular motion of a giant censer (O Botafumeiro) that hangs in the transept of the cathedral of Santiago de Compostela, and is cyclically pumped by men who pull at the supporting rope, is analyzed. Maximum angular amplitude attainable, and number of cycles and time needed to attain it, are calculated; the results agree with observed values (~ 82°, ~ 17 cycles, ~ 80 seconds) to the few percent accuracy of both the analysis and the observations and parameter measurements. The energy gain in a pumping cycle is obtained for an arbitrary pumping procedure to two orders in the small fractional change of pendular length; the relevance of the ratio (characteristic radial acceleration during pumping)/g to the gain is discussed- Effects due to rope mass, air drag on both Censer and rope, and the fact that the Censer is not a point mass, are considered. If the pumping cycle is inverted once the maximum amplitude has been attained, the Censer could be swiftly brought to rest, avoiding the usual violent stop. Historically recorded accidents, rope shape, and the influence of relevant parameters on the motion are discussed.

Bare-tether cathodic contact through thermoionic emission by low work-function materials

A new material, C12A7 : electride, which might present a work function as low as 0.6 eV and moderately high temperature stability, was recently proposed as coating for floating bare tethers. Arising from heating under space operation, current is emitted by thermionic emission along a thus coated cathodic segment. A preliminary study on the space-charge-limited (SCL) double layer in front of the cathodic segment is presented using Langmuir’s SCL electron current between cylindrical electrodes and orbital-motion-limited ion-collection sheath. A detailed calculation of current and bias profiles along the entire tether length is carried out with ohmic effects and the transition from SCL to full Richardson-Dushman emission included. Analysis shows that in the simplest drag mode, under typical orbital and tether conditions, thermionic emission leads to a short cathodic section and may eliminate the need for an active cathodic device and its corresponding gas feed requirements and power subsystem, which results in a truly “propellant-less” tether system for such basic applications as de-orbiting low earth orbit satellites.

Electrodynamic tethers for exploration of Jupiter and its icy moons

Use of electrodynamic bare tethers in exploring the Jovian system by tapping its rotational energy for power and propulsion is studied. The position of perijove and apojove in elliptical orbits, relative to the synchronous orbit at 2.24 times Jupiter’s radius, is exploited to conveniently make the induced Lorentz force to be drag or thrust, while generating power, and navigating the system. Capture and evolution to a low elliptical orbit near Jupiter, and capture into low circular orbits at moons Io and Europa are discussed.

Electrodynamic tether at Jupiter II:Fast moon tour after capture

An electrodynamic bare-tether mission to Jupiter,following the capture of a spacecraft (SC) into an equatorial highly elliptical orbit with perijove at about 1.3 times the Jovian radius, is discussed. Repeated applications of the propellantless Lorentz drag on a spinning tether, at the perijove vicinity, can progressively lower the apojove at constant perijove, for a tour of Galilean moons. Electrical energy is generated and stored as the SC moves from an orbit at 1 : 1 resonance with a moon, down to resonance with the next moon; switching tether current off, stored power is then used as the SC makes a number of flybys of each moon. Radiation dose is calculated throughout the mission,during capture, flybys and moves between moons. The tour mission is limited by both power needs and accumulated dose. The three-stage apojove lowering down to Ganymede, Io, and Europa resonances would total less than 14 weeks, while 4 Ganymede, 20 Europa, and 16 Io flybys would add up to 18 weeks, with the entire mission taking just over seven months and the accumulated radiation dose keeping under 3 Mrad (Si) at 10-mm Al shield thickness.

A proposed two-stage two-tether scientific mission at Jupiter

A two-stage mission to place a spacecraft (SC) below the Jovian radiation belts, using a spinning bare tether with plasma contactors at both ends to provide propulsion and power,is proposed. Capture by Lorentz drag on the tether, at the periapsis of a barely hyperbolic equatorial orbit, is followed by a sequence of orbits at near-constant periapsis, drag finally bringing the SC down to a circular orbit below the halo ring. Although increasing both tether heating and bowing, retrograde motion can substantially reduce accumulated dose as compared with prograde motion, at equal tether-to-SC mass ratio. In the second stage,the tether is cut to a segment one order of magnitude smaller, with a single plasma contactor, making the SC to slowly spiral inward over severalmonths while generating large onboard power, which would allow multiple scientific applications, including in situ study of Jovian grains, auroral sounding of upper atmosphere, and space- and time-resolved observations of surface and subsurface.

Efficient adaptive methods based on adjoint equations and truncation error estimation for unstructured grids

El propósito de esta tesis es la implementación de métodos eficientes de adaptación de mallas basados en ecuaciones adjuntas en el marco de discretizaciones de volúmenes finitos para mallas no estructuradas. La metodología basada en ecuaciones adjuntas optimiza la malla refinándola adecuadamente con el objetivo de mejorar la precisión de cálculo de un funcional de salida dado. El funcional suele ser una magnitud escalar de interés ingenieril obtenida por post-proceso de la solución, como por ejemplo, la resistencia o la sustentación aerodinámica. Usualmente, el método de adaptación adjunta está basado en una estimación a posteriori del error del funcional de salida mediante un promediado del residuo numérico con las variables adjuntas, “Dual Weighted Residual method” (DWR). Estas variables se obtienen de la solución del problema adjunto para el funcional seleccionado. El procedimiento habitual para introducir este método en códigos basados en discretizaciones de volúmenes finitos involucra la utilización de una malla auxiliar embebida obtenida por refinamiento uniforme de la malla inicial. El uso de esta malla implica un aumento significativo de los recursos computacionales (por ejemplo, en casos 3D el aumento de memoria requerida respecto a la que necesita el problema fluido inicial puede llegar a ser de un orden de magnitud). En esta tesis se propone un método alternativo basado en reformular la estimación del error del funcional en una malla auxiliar más basta y utilizar una técnica de estimación del error de truncación, denominada _ -estimation, para estimar los residuos que intervienen en el método DWR. Utilizando esta estimación del error se diseña un algoritmo de adaptación de mallas que conserva los ingredientes básicos de la adaptación adjunta estándar pero con un coste computacional asociado sensiblemente menor. La metodología de adaptación adjunta estándar y la propuesta en la tesis han sido introducidas en un código de volúmenes finitos utilizado habitualmente en la industria aeronáutica Europea. Se ha investigado la influencia de distintos parámetros numéricos que intervienen en el algoritmo. Finalmente, el método propuesto se compara con otras metodologías de adaptación de mallas y su eficiencia computacional se demuestra en una serie de casos representativos de interés aeronáutico. ABSTRACT The purpose of this thesis is the implementation of efficient grid adaptation methods based on the adjoint equations within the framework of finite volume methods (FVM) for unstructured grid solvers. The adjoint-based methodology aims at adapting grids to improve the accuracy of a functional output of interest, as for example, the aerodynamic drag or lift. The adjoint methodology is based on the a posteriori functional error estimation using the adjoint/dual-weighted residual method (DWR). In this method the error in a functional output can be directly related to local residual errors of the primal solution through the adjoint variables. These variables are obtained by solving the corresponding adjoint problem for the chosen functional. The common approach to introduce the DWR method within the FVM framework involves the use of an auxiliary embedded grid. The storage of this mesh demands high computational resources, i.e. over one order of magnitude increase in memory relative to the initial problem for 3D cases. In this thesis, an alternative methodology for adapting the grid is proposed. Specifically, the DWR approach for error estimation is re-formulated on a coarser mesh level using the _ -estimation method to approximate the truncation error. Then, an output-based adaptive algorithm is designed in such way that the basic ingredients of the standard adjoint method are retained but the computational cost is significantly reduced. The standard and the new proposed adjoint-based adaptive methodologies have been incorporated into a flow solver commonly used in the EU aeronautical industry. The influence of different numerical settings has been investigated. The proposed method has been compared against different grid adaptation approaches and the computational efficiency of the new method has been demonstrated on some representative aeronautical test cases.

Magnetic self-field effects on bare-tether current collection

It has been recently suggested that the magnetic field created by the current in a bare tether could sensibly reduce its electron collection capability in the magnetised ionosphere, a region of closed magnetic surfaces disconnecting the cylinder from infinity. In this paper, the ohmic voltage drop along the tether is taken into account in considering self-field effects. Separate analyses are carried out for the thrust and power generation and drag modes of operation, which are affected in different ways. In the power generation and drag modes, bias decreases as current increases along the tether, starting at the anodic, positively-biased end (upper end in the usual, eastward-flying spacecraft); in the thrust mode of operation, bias increases as current increases along the tether, starting at the lower end. When the ohmic voltage drop is considered, self-field effects are shown to be weak, in all cases, for tape tethers, and for circular cross-section tethers just conductive in a thin outer layer. Self-field effects might become important, in the drag case only, for tethers with fully conductive cross sections that are unrealistically heavy.

Short electrodynamic tethers

ED bare theters are best systems to deorbit S/C at end of service. For near polar orbits, usual tethers kept vertical by the gravity gradient, yield too weak magnetic drag. Here we propose keeping tethers perpendicular to the orbital plane. they mus be rigid and short for structural reasons, requiring power supply like Ion thrusters. terher tube-booms that can be rolled up on a drum would lie on each side of the S/C. One boom, carying in idle Hollow Cathode, collects electrons; the opposite boom's HC ejects electrons.

A "Free-Lunch" tour of the Jovian System

An ED-tether mission to Jupiter is presented. A bare tether carrying cathodic devices at both ends but no power supply, and using no propellant, could move 'freely' among Jupiter's 4 great moons. The tour scheme would have current naturally driven throughout by the motional electric field, the Lorentz force switching direction with current around a 'drag' radius of 160,00 kms, where the speed of the jovian ionosphere equals the speed of a spacecraft in circular orbit. With plasma density and magnetic field decreasing rapidly with distance from Jupiter, drag/thrust would only be operated in the inner plasmasphere, current being near shut off conveniently in orbit by disconnecting cathodes or plugging in a very large resistance; the tether could serve as its own power supply by plugging in an electric load where convenient, with just some reduction in thrust or drag. The periapsis of the spacecraft in a heliocentric transfer orbit from Earth would lie inside the drag sphere; with tether deployed and current on around periapsis, magnetic drag allows Jupiter to capture the spacecraft into an elliptic orbit of high eccentricity. Current would be on at succesive perijove passes and off elsewhere, reducing the eccentricity by lowering the apoapsis progressively to allow visits of the giant moons. In a second phase, current is on around apoapsis outside the drag sphere, rising the periapsis until the full orbit lies outside that sphere. In a third phase, current is on at periapsis, increasing the eccentricity until a last push makes the orbit hyperbolic to escape Jupiter. Dynamical issues such as low gravity-gradient at Jupiter and tether orientation in elliptic orbits of high eccentricity are discussed.

Electrodynamic tether at Jupiter 2. Tour missions after capture

Three separate scenarios of an electrodynamic tether mission at Jupiter following capture of a spacecraft (SC) into an equatorial, highly elliptical orbit around the planet, with perijove at about 1.5 times the Jovian radius, are discussed. Repeated application of Lorentz drag on the spinning tether, at the perijove vicinity, can progressively lower the apojove. One mission involves the tethered-SC rapidly and frequently visiting Galilean moons; elliptical orbits with apojove down at the Ganymede, Europa, and Io orbits are in 2:5, 4:9, and 1:2 resonances with the respective moons. About 20 slow flybys of Io would take place before the accumulated radiation dose exceeds 3 Mrad (Si) at 10 mm Al shield thickness, with a total duration of 5 months after capture (4 months for lowering the apojove to Io and one month for the flybys). The respective number of flybys for Ganymede would be 10 with a total duration of about 9 months. An alternative mission would have the SC acquire a low circular orbit around Jupiter, below the radiation belts, and manoeuvre to get an optimal altitude, with no major radiation effects, in less than 5 months after capture. In a third mission, repeated thrusting at the apojove vicinity, once down at the Io torus, would raise the perijove itself to the torus to acquire a low circular orbit around Io in about 4 months, for a total of 8 months after capture; this corresponds, however, to over 100 apojove passes with an accumulated dose, of about 8.5 Mrad (Si), that poses a critical issue.

Propulsive small expendable deployer system experiment

Relatively short electrodynamic tethers can extract orbital energy to "push" against a planetary magnetic field to achieve propulsion without the expenditure of propellant. The Propulsive Small Expendable Deployer System experiment will use the flight-proven Small Expendable Deployer System to deploy a 5-km bare aluminum tether from a Delta II upper stage to achieve ~0.4-N drag thrust, thus lowering the altitude of the stage. The experiment will use a predominantly bare tether for current collection in lieu of the endmass collector and insulated tether used on previous missions. The flight experiment is a precursor to a more ambitious electrodynamic tether upper-stage demonstration mission that will be capable of orbit-raising,lowering, and inclination changes, all using electrodynamic thrust. The expected performance of the tether propulsion system during the experiment is described.

Bare Tethers for Electrodynamic Spacecraft Propulsion

Electrodynamic tether thrusters can use the power provided by solar panels to drive a current in the tether and then the Lorentz force to push against the Earth's magnetic field, thereby achieving propulsion without the expenditure of onboard energy sources or propellant. Practical tether propulsion depends critically on being able to extract multiamp electron currents from the ionosphere with relatively short tethers (10 km or less) and reasonably low power. We describe a new anodic design that uses an uninsulated portion of the metallic tether itself to collect electrons. Because of the efficient collection of this type of anode, electrodynamic thrusters for reboost of the International Space Station and for an upper stage capable of orbit raising, lowering, and inclination changes appear to be feasible. Specifically, a 10-km-long bare tether, utilizing 10 kW of the space station power could save most of the propellant required for the station reboost over its 10-year lifetime. The propulsive small expendable deployer system experiment is planned to test the bare-tether design in space in the year 2000 by deploying a 5-km bare aluminum tether from a Delta II upper stage to achieve up to 0.5-N drag thrust, thus deorbiting the stage.

Multi-fluid Eulerian model of an electrospray in a host gas.

An Eulerian multifluid model is used to describe the evolution of an electrospray plume and the flow induced in the surrounding gas by the drag of the electrically charged spray droplets in the space between an injection electrode containing the electrospray source and a collector electrode. The spray is driven by the voltage applied between the two electrodes. Numerical computations and order-of-magnitude estimates for a quiescent gas show that the droplets begin to fly back toward the injection electrode at a certain critical value of the flux of droplets in the spray, which depends very much on the electrical conditions at the injection electrode. As the flux is increased toward its critical value, the electric field induced by the charge of the droplets partially balances the field due to the applied voltage in the vicinity of the injection electrode, leading to a spray that rapidly broadens at a distance from its origin of the order of the stopping distance at which the droplets lose their initial momentum and the effect of their inertia becomes negligible. The axial component of the electric field first changes sign in this region, causing the fly back. The flow induced in the gas significantly changes this picture in the conditions of typical experiments. A gas plume is induced by the drag of the droplets whose entrainment makes the radius of the spray away from the injection electrode smaller than in a quiescent gas, and convects the droplets across the region of negative axial electric field that appears around the origin of the spray when the flux of droplets is increased. This suppresses fly back and allows much higher fluxes to be reached than are possible in a quiescent gas. The limit of large droplet-to-gas mass ratio is discussed. Migration of satellite droplets to the shroud of the spray is reproduced by the Eulerian model, but this process is also affected by the motion of the gas. The gas flow preferentially pushes satellite droplets from the shroud to the core of the spray when the effect of the inertia of the droplets becomes negligible, and thus opposes the well-established electrostatic/inertial mechanism of segregation and may end up concentrating satellite droplets in an intermediate radial region of the spray.