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.

Adaptive scheme for Accurate orbit propagation

Two extensions of the fast and accurate special perturbation method recently developed by Peláez et al. are presented for respectively elliptic and hyperbolic motion. A comparison with Peláez?s method and with the very efficient Stiefel- Scheifele?s method, for the problems of oblate Earth plus Moon and continuous radial thrust, shows that the new formulations can appreciably improve the accuracy of Peláez?s method and have a better performance of Stiefel-Scheifele?s method. Future work will be to include the two new formulations and the original one due to Peláez into an adaptive scheme for highly accurate orbit propagation

The theory of electrostatic probes in strong magnetic fields : Thesis submitted to the Faculty of The Graduate School of the University of Colorado for the Degree Doctor of Philosophy Department of Aerospace Engineering Science

A theory is developed of an electrostatic probe in a fully-ionized plasma in the presence of a strong magnetic field. The ratio of electron Larmor radius to probe transverse dimension is assumed to be small. Poisson's equation, together with kinetic equations for ions and electrons are considered. An asymptotic perturbation method of multiple scales is used by considering the characteristic lengths appearing in the problem. The leading behavior of the solution is found. The results obtained appear to apply to weaker fields also, agreeing with the solutions known in the limit of no magnetic field. The range of potentials for wich results are presented is limited. The basic effects produced by the field are a depletion of the plasma near the probe and a non-monotonic potential surrounding the probe. The ion saturation current is not changed but changes appear in both the floating potential Vf and the slope of the current-voltage diagram at Vf. The transition region extends beyond the space potential Vs,at wich point the current is largely reduced. The diagram does not have an exponential form in this region as commonly assumed. There exists saturation in electron collection. The extent to which the plasma is disturbed is determined. A cylindrical probe has no solution because of a logarithmic singularity at infinity. Extensions of the theory are considered.

Zero-G Gauging systems. CONTRACT No 17/8/4

Various systems for measuring propellant content in spacecrafts under weightlessness conditions are reviewed. The cavity resonator method is found to be the most suitable measurement; technique. This method is analyzed in detail. A determination of errors intrinsec to the method is carried out.

A new set of integrals of motion to propagate the perturbed two-body problem

A formulation of the perturbed two-body problem that relies on a new set of orbital elements is presented. The proposed method represents a generalization of the special perturbation method published by Peláez et al. (Celest Mech Dyn Astron 97(2):131?150,2007) for the case of a perturbing force that is partially or totally derivable from a potential. We accomplish this result by employing a generalized Sundman time transformation in the framework of the projective decomposition, which is a known approach for transforming the two-body problem into a set of linear and regular differential equations of motion. Numerical tests, carried out with examples extensively used in the literature, show the remarkable improvement of the performance of the new method for different kinds of perturbations and eccentricities. In particular, one notable result is that the quadratic dependence of the position error on the time-like argument exhibited by Peláez?s method for near-circular motion under the J2 perturbation is transformed into linear.Moreover, themethod reveals to be competitive with two very popular elementmethods derived from theKustaanheimo-Stiefel and Sperling-Burdet regularizations.

Influencia de un fluido en las características dinámicas de placas circulares y casicirculares

La influencia de un fluido en las características dinámicas de estructuras se ha estudiado desde hace tiempo. Sin embargo muchos estudios se refieren a aplicaciones bajo el agua, como es el caso del sonar de un submarino por lo que el fluido circundante se considera líquido (sin efectos de compresibilidad). Más recientemente en aplicaciones acústicas y espaciales tales como antenas o paneles muy ligeros, ha sido estudiada la influencia en las características dinámicas de una estructura rodeada por un fluido de baja densidad. Por ejemplo se ha mostrado que el efecto del aire en el transmisor-reflector del Intelsat VI C-B con un diámetro de 3,2 metros y con un peso de sólo 34,7 kg disminuye la primera frecuencia en torno a un 20% con respecto a su valor en vacío. Por tanto es importante en el desarrollo de estas grandes y ligeras estructuras disponer de un método con el que estimar el efecto del fluido circundante sobre las frecuencias naturales de éstas. De esta manera se puede evitar el ensayo de la estructura en una cámara de vacío que para el caso de una gran antena o panel puede ser difícil y costoso. Se ha desarrollado un método de elementos de contorno (BEM) para la determinación del efecto del fluido en las características dinámicas de una placa circular. Una vez calculados analíticamente los modos de vibración de la placa en vacío, la matriz de masa añadida debido a la carga del fluido se determina por el método de elementos de contorno. Este método utiliza anillos circulares de manera que el número de elementos para obtener unos resultados precisos es muy bajo. Se utiliza un procedimiento de iteración para el cálculo de las frecuencias naturales del acoplamiento fluido-estructura para el caso de fluido compresible. Los resultados del método se comparan con datos experimentales y otros modelos teóricos mostrando la precisión y exactitud para distintas condiciones de contorno de la placa. Por otro lado, a veces la geometría de la placa no es circular sino casi-circular y se ha desarrollado un método de perturbaciones para determinar la influencia de un fluido incompresible en las características dinámicas de placas casi-circulares. El método se aplica a placas con forma elíptica y pequeña excentricidad. Por una parte se obtienen las frecuencias naturales y los modos de deformación de la placa vibrando en vacío. A continuación, se calculan los coeficientes adimensionales de masa virtual añadida (factores NAVMI). Se presentan los resultados de estos factores y el efecto del fluido en las frecuencias naturales. ABSTRACT The influence of the surrounding fluid on the dynamic characteristics of structures has been well known for many years. However most of these works were more concerned with underwater applications, such as the sonar of a submarine and therefore the surrounding fluid was considered a liquid (negligible compressibility effects). Recently for acoustical and spatial applications such as antennas or very light panels the influence on the dynamic characteristics of a structure surrounded by a fluid of low density has been studied. Thus it has been shown that the air effect for the Intelsat VI C-B transmit reflector with a diameter of 3,2 meters and weighting only 34,7 kg decreases the first modal frequency by 20% with respect to the value in vacuum. It is important then, in the development of these light and large structures to have a method that estimates the effect that the surrounding fluid will have on the natural frequencies of the structure. In this way it can be avoided to test the structure in a vacuum chamber which for a large antenna or panel can be difficult and expensive A BEM method for the determination of the effect of the surrounding fluid on the dynamic characteristics of a circular plate has been developed. After the modes of the plate in vacuum are calculated in an analytical form, the added mass matrix due to the fluid loading is determined by a boundary element method. This method uses circular rings so the number of elements to obtain an accurate result is very low. An iteration procedure for the computation of the natural frequencies of the couple fluid-structure system is presented for the case of the compressibility effect of air. Comparisons of the present method with various experimental data and other theories show the efficiency and accuracy of the method for any support condition of the plate. On the other hand, sometimes the geometry of the plate is not circular but almost-circular, so a perturbation method is developed to determine the influence of an incompressible fluid on the dynamic characteristics of almost-circular plates. The method is applied to plates of elliptical shape with low eccentricity. First, the natural frequencies and the mode shapes of the plate vibrating in vacuum are obtained. Next, the nondimensional added virtual mass coefficients (NAVMI factors) are calculated. Results of this factors and the effect of the fluid on the natural frequencies are presented.

One-dimensional self-similar solution of the dynamics of axisymmetric slender liquid bridges

In this paper the dynamics of axisymmetric, slender, viscous liquid bridges having volume close to the cylindrical one, and subjected to a small gravitational field parallel to the axis of the liquid bridge, is considered within the context of one-dimensional theories. Although the dynamics of liquid bridges has been treated through a numerical analysis in the inviscid case, numerical methods become inappropriate to study configurations close to the static stability limit because the evolution time, and thence the computing time, increases excessively. To avoid this difficulty, the problem of the evolution of these liquid bridges has been attacked through a nonlinear analysis based on the singular perturbation method and, whenever possible, the results obtained are compared with the numerical ones.

DROMO formulation for planar motions: Solution to the Tsien Problem

The two-body problem subject to a constant radial thrust is analyzed as a planar motion. The description of the problem is performed in terms of three perturbation methods: DROMO and two others due to Deprit. All of them rely on Hansen?s ideal frame concept. An explicit, analytic, closed-form solution is obtained for this problem when the initial orbit is circular (Tsien problem), based on the DROMO special perturbation method, and expressed in terms of elliptic integral functions. The analytical solution to the Tsien problem is later used as a reference to test the numerical performance of various orbit propagation methods, including DROMO and Deprit methods, as well as Cowell and Kustaanheimo?Stiefel methods.

EDROMO: An accurate propagator for elliptical orbits in the perturbed two-body problem

EDROMO is a special perturbation method for the propagation of elliptical orbits in the perturbed two-body problem. The state vector consists of a time-element and seven spatial elements, and the independent variable is a generalized eccentric anomaly introduced through a Sundman time transformation. The key role in the derivation of the method is played by an intermediate reference frame which enjoys the property of remaining fixed in space as long as perturbations are absent. Three elements of EDROMO characterize the dynamics in the orbital frame and its orientation with respect to the intermediate frame, and the Euler parameters associated to the intermediate frame represent the other four spatial elements. The performance of EDromo has been analyzed by considering some typical problems in astrodynamics. In almost all our tests the method is the best among other popular formulations based on elements.

A De-Embedding Method to Characterize Differential

This paper analyses the noise and gain measurement of microwave differential amplifiers using two passive baluns. A model of the baluns that includes potential losses and unbalances has been considered. This analysis allows to de-embed the actual performance of the differential device from the single-ended measurements of the two-port cascaded system and the baluns. The method has been validated with measured results from a fully-differential amplifier prototype.

Effect of Aspect Ratio on the Three-Dimensional Global Instability Analysis of Incompressible Open Cavity Flows

The stability analysis of open cavity flows is a problem of great interest in the aeronautical industry. This type of flow can appear, for example, in landing gears or auxiliary power unit configurations. Open cavity flows is very sensitive to any change in the configuration, either physical (incoming boundary layer, Reynolds or Mach numbers) or geometrical (length to depth and length to width ratio). In this work, we have focused on the effect of geometry and of the Reynolds number on the stability properties of a threedimensional spanwise periodic cavity flow in the incompressible limit. To that end, BiGlobal analysis is used to investigate the instabilities in this configuration. The basic flow is obtained by the numerical integration of the Navier-Stokes equations with laminar boundary layers imposed upstream. The 3D perturbation, assumed to be periodic in the spanwise direction, is obtained as the solution of the global eigenvalue problem. A parametric study has been performed, analyzing the stability of the flow under variation of the Reynolds number, the L/D ratio of the cavity, and the spanwise wavenumber β. For consistency, multidomain high order numerical schemes have been used in all the computations, either basic flow or eigenvalue problems. The results allow to define the neutral curves in the range of L/D = 1 to L/D = 3. A scaling relating the frequency of the eigenmodes and the length to depth ratio is provided, based on the analysis results.

A three-dimensional self-adaptive hp finite element method for the characterization of waveguide discontinuities

We propose the use of a highly-accurate three-dimensional (3D) fully automatic hp-adaptive finite element method (FEM) for the characterization of rectangular waveguide discontinuities. These discontinuities are either the unavoidable result of mechanical/electrical transitions or deliberately introduced in order to perform certain electrical functions in modern communication systems. The proposed numerical method combines the geometrical flexibility of finite elements with an accuracy that is often superior to that provided by semi-analytical methods. It supports anisotropic refinements on irregular meshes with hanging nodes, and isoparametric elements. It makes use of hexahedral elements compatible with high-order H(curl)H(curl) discretizations. The 3D hp-adaptive FEM is applied for the first time to solve a wide range of 3D waveguide discontinuity problems of microwave communication systems in which exponential convergence of the error is observed.

Instability Analysis of Incompressible Open Cavity Flows

El problema del flujo sobre una cavidad abierta ha sido estudiado en profundidad en la literatura, tanto por el interés académico del problema como por sus aplicaciones prácticas en gran variedad de problemas ingenieriles, como puede ser el alojamiento del tren de aterrizaje de aeronaves, o el depósito de agua de aviones contraincendios. Desde hace muchos a˜nos se estudian los distintos tipos de inestabilidades asociadas a este problema: los modos bidimensionales en la capa de cortadura, y los modos tridimensionales en el torbellino de recirculación principal dentro de la cavidad. En esta tesis se presenta un estudio paramétrico completo del límite incompresible del problema, empleando la herramienta de estabilidad lineal conocida como BiGlobal. Esta aproximación permite contemplar la estabilidad global del flujo, y obtener tanto la forma como las características de los modos propios del problema físico, sean estables o inestables. El estudio realizado permite caracterizar con gran detalle todos los modos relevantes, así como la envolvente de estabilidad en el espacio paramétrico del problema incompresible (Mach nulo, variación de Reynolds, espesor de capa límite incidente, relación altura/profundidad de la cavidad, y longitud característica de la perturbación en la dirección transversal). A la luz de los resultados obtenidos se proponen una serie de relaciones entre los parámetros y características de los modos principales, como por ejemplo entre el Reynolds crítico de un modo, y la longitud característica del mismo. Los resultados numéricos se contrastan con una campaña experimental, siendo la principal conclusión de dicha comparación que los modos lineales están presentes en el flujo real saturado, pero que existen diferencias notables en frecuencia entre las predicciones teóricas y los experimentos. Para intentar determinar la naturaleza de dichas diferencias se realiza una simulación numérica directa tridimensional, y se utiliza un algoritmo de DMD (descomposición dinámica de modos) para describir el proceso de saturación. ABSTRACT The problem of the flow over an open cavity has been studied in depth in the literature, both for being an interesting academical problem and due to the multitude of industrial applications, like the landing gear of aircraft, or the water deposit of firefighter airplanes. The different types of instabilities appearing in this flow studied in the literature are two: the two-dimensional shear layer modes, and the three-dimensional modes that appear in the main recirculating vortex inside the cavity. In this thesis a parametric study in the incompressible limit of the problem is presented, using the linear stability analysis known as BiGlobal. This approximation allows to obtain the global stability behaviour of the flow, and to capture both the morphological features and the characteristics of the eigenmodes of the physical problem, whether they are stable or unstable. The study presented here characterizes with great detail all the relevant eigenmodes, as well as the hypersurface of instability on the parameter space of the incompressible problem (Mach equal to zero, and variation of the Reynolds number, the incoming boundary layer thickness, the length to depth aspect ratio of the cavity and the spanwise length of the perturbation). The results allow to construct parametric relations between the characteristics of the leading eigenmodes and the parameters of the problem, like for example the one existing between the critical Reynolds number and its characteristic length. The numerical results presented here are compared with those of an experimental campaign, with the main conclusion of said comparison being that the linear eigenmode are present in the real saturated flow, albeit with some significant differences in the frequencies of the experiments and those predicted by the theory. To try to determine the nature of those differences a three-dimensional direct numerical simulation, analyzed with Dynamic Mode Decomposition algorithm, was used to describe the process of saturation.

A new regularization method for fast and accurate propagation of roto-translationally coupled asteroids

There is a well-distinguished group of asteroids for which the roto-translational cou-pling is known to have a non-negligible e�ect in the long-term. The study of such asteroids suggests the use of specialized propagation techniques, where perturbation methods make their best. The techniques from which the special regularization method DROMO is derived, have now been extended to the attitude dynamics, with equally remarkable results in terms of speed and accuracy, thus making the combination of these algorithms specially. well-suited to deal with the propagation of bodies with strong attitude coupling.