Literatura científica selecionada sobre o tema "Reusable Launchers"

The increasing market of small and affordable space systems requires fast and reliablelaunch capabilities to cover the current and future demand. This project aims to studyand implement guidance and control schemes for vertical ascent and descent phases ofa reusable launcher. Specifically, the thesis focuses on developing and applying ModelPredictive Control (MPC) and optimisation techniques to several kino-dynamic modelsof rockets. Moreover, the classical MPC method has been modified to include a decreasingfactor for the horizon used, enhancing the performance of the guidance and control.Multiple scenarios of vertical launch, landing and full fligth guidance on Earth and Mars,along with Monte Carlo analysis, were carried out to demonstrate the robustness of thealgorithm against different initial conditions. The results were promising and invite tofurther research.

Dans le cadre du développement de moteurs de fusée réutilisables, les exigences de fonctionnement des différents éléments composant un moteur ont connu de grandes évolutions. Alors qu'un moteur classique était conçu pour un nombre restreint de points de fonctionnement, un moteur réutilisable doit répondre à des exigences sur une large plage de points, afin d'effectuer des manoeuvres plus complexes. En conséquence, les lois de commande des moteurs fusées ont subi une évolution similaire, rendant nécessaire la loi de commande en boucle fermée. Bien que de nombreuses études aient été réalisées sur des lois de commande, peu de travaux portent sur la stabilité du moteur en boucle fermée. Dans cette optique, l'objectif de ces travaux est de proposer une démonstration de stabilité d'un modèle de moteur fusée, ainsi qu'un contrôleur permettant d'obtenir des garanties de stabilité du modèle. En premier lieu, un modèle typique de moteur de fusée à ergols liquide est développé, sous forme d'espace d'états. Ce type de modèle, bien que plus courant, se révèle peu adapté à l'étude de la stabilité, de par sa formulation hautement non-linéaire. Dans ce cadre, l'utilisation d'une fonction de Lyapunov se révèle complexe, et une reformulation du modèle est envisagée, sous forme d'un modèle Hamiltonien à ports. Un second chapitre permet d'introduire la notion de modèle Hamiltonien à ports. Ce type de modèle met en valeur les transferts énergétiques qui ont lieu entre les différents éléments d'un système, et sont construits avec une structure géométrique fixe. Ces différentes caractéristiques permettent une étude directe de la passivité d'un système, un outil d'analyse de la stabilité d'un système. La reformulation permet de trouver une fonction caractéristique d'un système Hamiltonien à ports, l'Hamiltonien, qui prouve la passivité d'un système et peut être formulé comme une fonction de Lyapunov. Cette démonstration donne des conditions de stabilité sur la modélisation du système, ainsi que sur le contrôleur appliqué en boucle fermée. Dans le cas où la démonstration directe de passivité n'est pas réalisable, un contrôleur peut être construit pour assurer la passivité de la boucle fermée. Pour conférer les propriétés de la passivité au modèle de moteur utilisé, la théorie du contrôle par passivité est présentée. Le principe d'un tel contrôleur est d'assurer la stabilité d'un système en rendant la boucle fermée passive. Avec la théorie des systèmes Hamiltonien à ports cependant, ce contrôleur permet aussi de modifier la structure géométrique hamiltonienne, afin de reformuler un système sous forme Hamiltonienne à ports. Ce contrôleur permet de rendre le système passif autour d'un point de fonctionnement désiré par l'utilisateur, qui peut être changé au cours du temps. Ainsi, ce contrôleur permet un suivi de trajectoire avec des garanties de passivité du système au cours du temps. Le quatrième chapitre propose une approche différente pour établir un contrôleur stabilisant, à l'aide de la théorie de la contraction. La propriété de contraction d'un système dénote sa capacité à converger rapidement vers une trajectoire de référence. Cette propriété constitue une forme de stabilité exponentielle, plus puissante que la stabilité par passivation. Le contrôleur peut de plus être réalisé aisément, en résolvant des inégalités linéaires matricielles. Enfin, les résultats de ces travaux sont présentés à l'aide de simulations sur MATLAB Simulink, et permettent de conclure sur les différents contrôleurs présentés. Un contrôleur simple proportionnel intégral dérivé (PID) est construit pour permettre une comparaison. Les résultats montrent que les contrôleurs réalisés proposent des propriétés stabilisantes, alors que le contrôleur PID est instable dans certaines zones de fonctionnement. Le contrôleur par passivité étend le domaine de stabilité du système, et le contrôleur par contraction empêche le système de quitter le domaine de stabilité du système original<br>With the development of reusable rocket engines, the operating requirements of the various components in an engine have significantly increased. While a non-reusable engine was designed for a limited number of operating points, a reusable engine must meet requirements over a wide range of points to perform complex maneuvers. Consequently, rocket engine control laws have evolved similarly, with the introduction of closed-loop control laws. Although many studies have been conducted on control laws, few works focus on the stability of the engine in closed-loop control. In this context, the objective of this work is to propose a demonstration of the stability of a rocket engine model, as well as a controller that guarantees the stability of the model. First, a model of a liquid propelled rocket engine is proposed under a state-space form. Although more common, this type of modeling does not allow for an easy stability analysis due to its highly nonlinear terms. In this context, the use of a Lyapunov function proves to be cumbersome, and a reformulation of the model is considered, in the form of a Port-Hamiltonian model, more suited for stability analysis of the system. A second chapter introduces the concept of the Port-Hamiltonian model. This type of model highlights the energy transfers that occur between the various components of a system and is built with a fixed geometric structure. These characteristics allow for a direct study of the passivity of a system, a tool for stability analysis the stability. The reformulation allows for the identification of a characteristic function of a Port-Hamiltonian system, the Hamiltonian function, which can be used to prove the passivity of a system and can be formulated as a Lyapunov function. This demonstration provides stability conditions for the system as well as the controller applied in the closed-loop system. In cases where a direct demonstration of passivity is not possible, a controller can be constructed to ensure the passivity of the closed-loop system. To endow the rocket engine model with passivity properties, the third chapter presents passivity-based control (PBC) theory. The principle of such a controller is to ensure the stability of a system by making the closed-loop system passive. Coupled with Port-Hamiltonian systems theory, however, this controller also allows for modification of the Hamiltonian geometric structure to reformulate a system into Port-Hamiltonian form. This controller makes the system passive around a desired operating point, which can be changed over time. Thus, this controller enables trajectory tracking with passivity guarantees over time. The fourth chapter proposes a different approach to establish a stabilizing controller using contraction theory. The contraction property of a system indicates its ability to rapidly converge towards a reference trajectory. This property represents a form of exponential stability, which is more robust than stability through passivation. Moreover, the controller can be easily implemented by solving linear matrix inequalities. Finally, the results of this work are presented through simulations on MATLAB Simulink, allowing for conclusions on the various controllers presented. A simple proportional-integralderivative (PID) controller is constructed for comparison. The results show that the designed controllers offer stabilizing properties, while the PID controller is unstable in certain operating regions. The passivity-based controller extends the stability domain of the system, and the contraction-based controller prevents the system from leaving the stability domain of the original system

The design of complex systems such as a launcher, or subsystems like its pressurization system, is known to be fastidious and expensive especially in the space domain. However, the recent emergence of new actors in this domain has been a game changer, binding these systems to fulfil more and more requirements (cheap, efficient, rapidly produced, aesthetic, environmentally friendly…) in order to compete with the market. For these reasons, engineers now need more than ever to consider the full picture of a system or subsystem in order to optimize it. In this context, this document aims to present the method and the results obtained in the optimization of the pressurization system of a reusable rocket first stage named Themis 3. The modelling of the pressurization system has been realized through a software called Geeglee. In order to cover as many impacts of this system as possible, an important part of the rocket stage has been considered and modelled (from the pressurization gases to the propellant tanks, passing by the pressurant and propellant feeding lines). Three High Level Requirements (HLR) have been identified as of major importance for the trade-off in the design of the pressurization system: the total mass impact, the total Recursive Cost (RC) impact and the total Non-Recursive Cost (NRC) impact. This optimization has in particular permitted to confirm some well-known results, so to say that exogenous pressurization systems represent a smaller mass impact on the vehicle, at the expense of a higher RC compared with autogenous systems.<br>Utformningen av komplexa system, som till exempel en bärraket, eller delsystem som dess trycksättningssystem, är känd för att vara krävande och kostsam, särskilt inom rymdteknikområdet. Den senaste tidens uppkomst av nya aktörer på detta område har dock förändrat spelregler då dessa system tvingas uppfylla allt fler krav (billiga, effektiva, snabbt producerade, estetiska, miljövänliga etc.) kunna vara konkurrenskraftiga på marknaden. Av dessa skäl måste ingenjörer nu mer än någonsin beakta hela bilden av ett system eller delsystem för att kunna optimera det. I detta sammanhang syftar detta dokument till att presentera metoden och resultaten från optimeringen av trycksättningssystemet för en återanvändbar rakets första steg som heter Themis 3. Modelleringen av trycksättningssystemet har genomförts med hjälp av ett systemmodelleringsverktyg som kallas Geeglee. För att täcka in så många effekter av detta system som möjligt har en viktig del av raketsteget beaktats och modellerats (från trycksättningsgaserna till drivmedelstankarna, via tryck- och drivmedelsmatningsledningarna). Tre krav på hög nivå har identifierats som mycket viktiga för avvägningen vid utformningen av trycksättningssystem konstruktion: den totala masspåverkan, den totala icke-rekursiva kostnaden och den totala rekursiva kostnaden. Denna optimering har framför allt gjort det möjligt att bekräfta vissa välkända resultat, dvs. att exogena trycksättningssystem ger en mindre masspåverkan på fordonet, på bekostnad av en högre RC jämfört med autogena system.

Starting with a brief overview of thrust generation for launchers, this study focuses on the design process of the demonstrator aerospike engine, DEMOP-1, of the Pangea Aerospace's commercial grade engine and its flow field analysis. The primary goal of the study is to obtain the plug nozzle design delivers 30 kN thrust using cryogenic liquid oxygen (LOX) as the oxidizer and cryogenic liquid methane (LCH4) as the fuel, with the mixture ratio of 3.4. Design parameters considered as 30 bar of combustion chamber pressure (Po) and expansion ratio as 15 for an optimum expanded nozzle. On the basis of decided design characteristics, Angelino's method is used to design the nozzle contour through MATLAB. The flow field over the aerospike analyzed using commercial CFD program FLUENT for sea level, optimum expansion and vacuum conditions. Flow simulations are carried out for air (specific heat ratio, gamma= 1.4), and afterwards based on the obtained thrust values at each altitude for air, expected thrust values for the real propellant, LOX/LCH4 (specific heat ratio, gamma = 1.1664), are calculated. Finally, the study is concluded with the comparison of trend in thrust and specific impulse for conventional bell nozzle and aerospike. For the conventional bell engine the values obtained in commercial computational simulation of chemical rocket propulsion and combustion software RPA for bell nozzle with same characteristics with aerospike, Po = 30 bar and expansion ratio = 15, are taken as reference for sea level, optimum expansion level and vacuum condition performance. Due to its ability to adopt the altitude, aerospike delivers higher performance at the low altitudes with respect to the conventional bell nozzle which has the same expansion ratio and combustion chamber pressure. Last in order but not in importance, after obtaining the flow field on plug of the aerospike, the shock wave impingement on the nozzle surface at sea level has been investigated.

Cette étude consiste à développer une méthode de prévision de la distribution du flux de chaleur surfacique pour des lanceurs réutilisables qui soit adaptée à une phase de conception préliminaire. Les phases de vol hypersonique les plus critiques lors de la réentrée atmosphérique du lanceur sont celles pour lesquelles les flux de chaleurs les plus importants apparaissent à la surface. Ainsi, pour les études de conception et de préparation de mission, il est indispensable de posséder des outils élaborés et rapides qui permettent d’obtenir des estimations de ces flux de chaleur surfaciques. Cela permet ainsi un dimensionnement rapide du design étudié car le flux de chaleur est une données d’entrée pour une multitude de spécialités: l’aérothermodynamique, l’optimisation de la trajectoire, la conception structurelle, etc. L’objectif principal de cette thèse est d’évaluer si l’utilisation des méthodes de « Machine Learning » apporte une plus-value pour le développent d’un modèle de prévision de la distribution du flux de chaleur convectif-diffusif surfacique. Pour cela, nous proposons une méthode qui consiste à développer deux modèles de substitution consécutifs. Le premier modèle de substitution est établi pour effectuer un adimensionnement des variables. Cela nous permet par la suite d’utiliser le modèle indépendamment des conditions de vol et d’altitude. Le deuxième modèle de substitution basé sur les réseaux de neurones établi la corrélation entre les variables d’entrée – topologie de pression surfacique et variables géométriques – et le flux de chaleur en sortie. Dans le cadre de cette thèse la méthode a été testée et développée pour un écoulement hypersonique laminaire en régime continue. Pour analyser ses capacités et limites, nous avons appliqué cette méthode sur l’étage orbital duconcept de lanceur « SpaceLiner ». Une base de données est établie avec les calculs de Navier-Stokes sur un ensemble de points de vols et d’orientations choisi. Les variables d’entrées sont adimensionnées à partir d’équations existantes ou de transformations mathématiques simples. Pour l’adimensionnement de flux de chaleur un modèle de substitution pour la prévision du flux de chaleur au point d’arrêt pour notre domaine d’intérêt est développé. Ce modèle est basé sur le modèle de Lepage-Vérant et du Krigeage. Cela nous permet d’obtenir une prévision avec une erreur moyenne relative de 1.9% par rapport aux résultats de CFD. Le développement puis l’évaluation du modèle de substitution pour la prévision du flux de chaleur en tout point de la paroi sont effectués. Nous avons comparé nos résultats à ceux issus de calculs CFD à la fois dans des zones convexes et –partiellement– dans des zones planes de façon satisfaisante avec une précision supérieure ou égale à celle des méthodes existantes. Des réseaux de neurones entrainés sur un petit nombre de points de vol sont capables d’interpoler et d’extrapoler les résultats à d’autres points de vol à proximité. Cependant, dans des zones à flux de chaleur adimensionné très faible (inférieur à 0.05) le modèle reste peu fiable mais cet aspect n’est pas essentiel dans une étude de conception préliminaire. De plus, la base de données n’est pas adaptée au développement d’un modèle de substitution pour la zone concave. Comme les modèles existants, la méthode proposée ne réussit pas à prédire les phénomènes physiques complexes comme des interactions choc-choc. En conclusion la méthode proposée a démontré son potentiel d’intégration dans des phases de conception préliminaire<br>For reusable launch vehicles critical aerothermal loads occur on the descent trajectory in the hypersonic continuum regime. Thus, for the pre-design phase of the mission preparation, it is essential to have advanced models with low response time (CPU time) at one’s disposal to predict the heat fluxes on the vehicle’s surface. This is necessary to make informed decisions on trade-offs required between system, aerothermal and trajectory optimisation aspects. The objective of this thesis is to investigate whether machine learning techniques offer additional value in the development of such a model predicting the wall heat flux distribution. For this, we propose a method which consists of two consecutive surrogate models. The first one provides a method of nondimensionalisation, to enable the prediction of heat flux independent of the flight points at different altitudes and freestream conditions. The second surrogate model maps four geometric and pressure based dimensionless input variables to the dimensionless heat flux output variable. This correlation is modelled with neural networks. In order to analyse the prediction accuracy, capabilities and limitations, we apply the proposed method to the SpaceLiner concept. For this, we build a database with Computational fluid dynamics (CFD) simulations which contains data from different flight points and orientations of the SpaceLiner orbiter stage. The nondimensionalisation model of the stagnation point heat flux is build based on the Lepage-Vérant model and Kriging and permits its prediction with an accuracy of a mean relative error of 1.9% compared to the results of CFD equations. The surrogate model based on neural networks for the prediction of the heat flux distribution at the current development state complies to our objective for the convex zones and partially for the flat zones if we disregard the very low heat flux zones and complex physical phenomena such as shock-shock interactions. These results ensure better or equivalent prediction accuracy as existing pre-design tools. The database was unsuited to developed a reliable model for the concave zones. Overall interpolation and close proximity extrapolation to different flight points for the given vehicle is possible. However further generalisation towards different vehicle shapes remain unreachable with the current setup. The proposed method has demonstrated potential to be integrated in the pre-design process