Asymptotic solution for the two-body problem with constant tangencial acceleration

An analytical solution of the two body problem perturbed by a constant tangential acceleration is derived with the aid of perturbation theory. The solution, which is valid for circular and elliptic orbits with generic eccentricity, describes the instantaneous time variation of all orbital elements. A comparison with high-accuracy numerical results shows that the analytical method can be effectively applied to multiple-revolution low-thrust orbit transfer around planets and in interplanetary space with negligible error.

CCSD(T) Study of CD3-O-CD3 and CH3-O-CD3 Far-Infrared Spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, 12CH3–16O–12CD3 (DME-d3) and 12CD3–16O–12CD3 (DME-d6), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d6 belongs to the G36 group, as the most abundant isotopologue 12CH3–16O–12CH3 (DME-h6), while DME-d3 is a G18 species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d6. Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N – 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d3 intensities and the band positions of the torsional overtones.

Linear Approach to the Orbiting Spacecraft Thermal Problem

A linear method is developed for solving the nonlinear differential equations of a lumped-parameter thermal model of a spacecraft moving in a closed orbit. This method, based on perturbation theory, is compared with heuristic linearizations of the same equations. The essential feature of the linear approach is that it provides a decomposition in thermal modes, like the decomposition of mechanical vibrations in normal modes. The stationary periodic solution of the linear equations can be alternately expressed as an explicit integral or as a Fourier series. This method is applied to a minimal thermal model of a satellite with ten isothermal parts (nodes), and the method is compared with direct numerical integration of the nonlinear equations. The computational complexity of this method is briefly studied for general thermal models of orbiting spacecraft, and it is concluded that it is certainly useful for reduced models and conceptual design but it can also be more efficient than the direct integration of the equations for large models. The results of the Fourier series computations for the ten-node satellite model show that the periodic solution at the second perturbative order is sufficiently accurate.

Optimization of multidimensional cross-section tables for few-group core calculations

Multigroup diffusion codes for three dimensional LWR core analysis use as input data pre-generated homogenized few group cross sections and discontinuity factors for certain combinations of state variables, such as temperatures or densities. The simplest way of compiling those data are tabulated libraries, where a grid covering the domain of state variables is defined and the homogenized cross sections are computed at the grid points. Then, during the core calculation, an interpolation algorithm is used to compute the cross sections from the table values. Since interpolation errors depend on the distance between the grid points, a determined refinement of the mesh is required to reach a target accuracy, which could lead to large data storage volume and a large number of lattice transport calculations. In this paper, a simple and effective procedure to optimize the distribution of grid points for tabulated libraries is presented. Optimality is considered in the sense of building a non-uniform point distribution with the minimum number of grid points for each state variable satisfying a given target accuracy in k-effective. The procedure consists of determining the sensitivity coefficients of k-effective to cross sections using perturbation theory; and estimating the interpolation errors committed with different mesh steps for each state variable. These results allow evaluating the influence of interpolation errors of each cross section on k-effective for any combination of state variables, and estimating the optimal distance between grid points.

An asymptotic solution for the main problem

In this paper, an analytical solution of the main problem, a satellite only perturbed by the J2 harmonic, is derived with the aid of perturbation theory and by using DROMO variables. The solution, which is valid for circular and elliptic orbits with generic eccentricity and inclination, describes the instantaneous time variation of all orbital elements, that is, the actual values of the osculating elements

Magnetic Attitude Control System for a Small Satellite. Impact on the Thermal Performance

El principal objetivo de la tesis es estudiar el acoplamiento entre los subsistemas de control de actitud y de control térmico de un pequeño satélite, con el fin de buscar la solución a los problemas relacionados con la determinación de los parámetros de diseño. Se considera la evolución de la actitud y de las temperaturas del satélite bajo la influencia de dos estrategias de orientación diferentes: 1) estabilización magnética pasiva de la orientación (PMAS, passive magnetic attitude stabilization), y 2) control de actitud magnético activo (AMAC, active magnetic attitude control). En primer lugar se presenta el modelo matemático del problema, que incluye la dinámica rotacional y el modelo térmico. En el problema térmico se considera un satélite cúbico modelizado por medio de siete nodos (seis externos y uno interno) aplicando la ecuación del balance térmico. Una vez establecido el modelo matemático del problema, se estudia la evolución que corresponde a las dos estrategias mencionadas. La estrategia PMAS se ha seleccionado por su simplicidad, fiabilidad, bajo coste, ahorrando consumo de potencia, masa coste y complejidad, comparado con otras estrategias. Se ha considerado otra estrategia de control que consigue que el satélite gire a una velocidad requerida alrededor de un eje deseado de giro, pudiendo controlar su dirección en un sistema inercial de referencia, ya que frecuentemente el subsistema térmico establece requisitos de giro alrededor de un eje del satélite orientado en una dirección perpendicular a la radiación solar incidente. En relación con el problema térmico, para estudiar la influencia de la velocidad de giro en la evolución de las temperaturas en diversos puntos del satélite, se ha empleado un modelo térmico linealizado, obtenido a partir de la formulación no lineal aplicando un método de perturbaciones. El resultado del estudio muestra que el tiempo de estabilización de la temperatura y la influencia de las cargas periódicas externas disminuye cuando aumenta la velocidad de giro. Los cambios de temperatura se reducen hasta ser muy pequeños para velocidades de rotación altas. En relación con la estrategia PMAC se ha observado que a pesar de su uso extendido entre los micro y nano satélites todavía presenta problemas que resolver. Estos problemas están relacionados con el dimensionamiento de los parámetros del sistema y la predicción del funcionamiento en órbita. Los problemas aparecen debido a la dificultad en la determinación de las características magnéticas de los cuerpos ferromagnéticos (varillas de histéresis) que se utilizan como amortiguadores de oscilaciones en los satélites. Para estudiar este problema se presenta un modelo analítico que permite estimar la eficiencia del amortiguamiento, y que se ha aplicado al estudio del comportamiento en vuelo de varios satélites, y que se ha empleado para comparar los resultados del modelo con los obtenidos en vuelo, observándose que el modelo permite explicar satisfactoriamente el comportamiento registrado. ABSTRACT The main objective of this thesis is to study the coupling between the attitude control and thermal control subsystems of a small satellite, and address the solution to some existing issues concerning the determination of their parameters. Through the thesis the attitude and temperature evolution of the satellite is studied under the influence of two independent attitude stabilization and control strategies: (1) passive magnetic attitude stabilization (PMAS), and (2) active magnetic attitude control (AMAC). In this regard the mathematical model of the problem is explained and presented. The mathematical model includes both the rotational dynamics and the thermal model. The thermal model is derived for a cubic satellite by solving the heat balance equation for 6 external and 1 internal nodes. Once established the mathematical model of the problem, the above mentioned attitude strategies were applied to the system and the temperature evolution of the 7 nodes of the satellite was studied. The PMAS technique has been selected to be studied due to its prevalent use, simplicity, reliability, and cost, as this strategy significantly saves the overall power, weight, cost, and reduces the complexity of the system compared to other attitude control strategies. In addition to that, another control law that provides the satellite with a desired spin rate along a desired axis of the satellite, whose direction can be controlled with respect to the inertial reference frame is considered, as the thermal subsystem of a satellite usually demands a spin requirement around an axis of the satellite which is positioned perpendicular to the direction of the coming solar radiation. Concerning the thermal problem, to study the influence of spin rate on temperature evolution of the satellite a linear approach of the thermal model is used, which is based on perturbation theory applied to the nonlinear differential equations of the thermal model of a spacecraft moving in a closed orbit. The results of this study showed that the temperature stabilization time and the periodic influence of the external thermal loads decreases by increasing the spin rate. However, the changes become insignificant for higher values of spin rate. Concerning the PMAS strategy, it was observed that in spite of its extended application to micro and nano satellites, still there are some issues to be solved regarding this strategy. These issues are related to the sizing of its system parameters and predicting the in-orbit performance. The problems were found to be rooted in the difficulties that exist in determining the magnetic characteristics of the ferromagnetic bodies (hysteresis rods) that are applied as damping devices on-board satellites. To address these issues an analytic model for estimating their damping efficiency is proposed and applied to several existing satellites in order to compare the results with their respective in-flight data. This model can explain the behavior showed by these satellites.

V-substituted In2S3: an intermediate band material with photocatalytic activity in the whole visible light range

We proposed in our previous work V-substituted In2S3 as an intermediate band (IB) material able to enhance the efficiency of photovoltaic cells by combining two photons to achieve a higher energy electron excitation, much like natural photosynthesis. Here this hyper-doped material is tested in a photocatalytic reaction using wavelength-controlled light. The results evidence its ability to use photons with wavelengths of up to 750 nm, i.e. with energy significantly lower than the bandgap (=2.0 eV) of non-substituted In2S3, driving with them the photocatalytic reaction at rates comparable to those of non-substituted In2S3 in its photoactivity range (λ ≤ 650 nm). Photoluminescence spectra evidence that the same bandgap excitation as in V-free In2S3 occurs in V-substituted In2S3 upon illumination with photons in the same sub-bandgap energy range which is effective in photocatalysis, and its linear dependence on light intensity proves that this is not due to a nonlinear optical property. This evidences for the first time that a two-photon process can be active in photocatalysis in a single-phase material. Quantum calculations using GW-type many-body perturbation theory suggest that the new band introduced in the In2S3 gap by V insertion is located closer to the conduction band than to the valence band, so that hot carriers produced by the two-photon process would be of electron type; they also show that the absorption coefficients of both transitions involving the IB are of significant and similar magnitude. The results imply that V-substituted In2S3, besides being photocatalytically active in the whole visible light range (a property which could be used for the production of solar fuels), could make possible photovoltaic cells of improved efficiency.

Hydrodynamic stability theory of double ablation front structures in inertial confinement fusion.

The aim of inertial confinement fusion is the production of energy by the fusion of thermonuclear fuel (deuterium-tritium) enclosed in a spherical target due to its implosion. In the direct-drive approach, the energy needed to spark fusion reactions is delivered by the irradiation of laser beams that leads to the ablation of the outer shell of the target (the so-called ablator). As a reaction to this ablation process, the target is accelerated inwards, and, provided that this implosion is sufficiently strong a symmetric, the requirements of temperature and pressure in the center of the target are achieved leading to the ignition of the target (fusion). One of the obstacles capable to prevent appropriate target implosions takes place in the ablation region where any perturbation can grow even causing the ablator shell break, due to the ablative Rayleigh-Taylor instability. The ablative Rayleigh-Taylor instability has been extensively studied throughout the last 40 years in the case where the density/temperature profiles in the ablation region present a single front (the ablation front). Single ablation fronts appear when the ablator material has a low atomic number (deuterium/tritium ice, plastic). In this case, the main mechanism of energy transport from the laser energy absorption region (low density plasma) to the ablation region is the electron thermal conduction. However, recently, the use of materials with a moderate atomic number (silica, doped plastic) as ablators, with the aim of reducing the target pre-heating caused by suprathermal electrons generated by the laser-plasma interaction, has demonstrated an ablation region composed of two ablation fronts. This fact appears due to increasing importance of radiative effects in the energy transport. The linear theory describing the Rayleigh-Taylor instability for single ablation fronts cannot be applied for the stability analysis of double ablation front structures. Therefore, the aim of this thesis is to develop, for the first time, a linear stability theory for this type of hydrodynamic structures.

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.

Cancellation zone in supersonic lifting wing theory

BASING their work on a linear theory, Evvard1 and Krasilshchikova2'3 independently developed an expression that yields the perturbation generated by a thiri lifting wing of arbitrary planform flying at supersonic speed on a point placed on the wing plane inside its planform,1 or both on and above the wing plane.2 This point must be influenced by two leading edges, one supersonic and the other partially subsonic. Although these authors followed different approaches, their methods concur in showing the existence of a perfectly defined cancellation zone. In this Note, the Evvard approach is generalized to the case solved by Krasilshchikova. Circumventing the latter's lengthy and somewhat complex approach, Evvard's simple method seems to be useful at least for educational purposes.