Due to a trend towards Ultra High Bypass Ratio engines the corresponding engine/airframe interference is becoming a key aspect in aircraft design. The present economic situation increases the pressure on commercial aviation companies to reduce the Direct Operating Cost, and the environmental situation requires a new generation of aircraft with a lower environmental impact. Therefore detailed aerodynamic investigations are required to evaluate the real benefits of new technologies. The presented research activity is part of a long-term project with the main objective of generating a reliable and accurate tool to predict the performance of an aircraft over the whole flight domain. In particular the aim of this research was to perform advanced CFD in order to establish a tool able to evaluate engine installation effects for different configurations and attitudes. The developed tool can be provided with correlations of the Net Propulsive Force (NPF), the force exerted by the power-plant to the aircraft, as a function of position. This can be done in principle at cruise, hold, climb, descent, take-off and landing, to model the different integration effects at different phases. Due to the complexity of the problem it was only possible at an initial stage to determine these correlations at cruise condition. Two parametric test cases were evaluated, showing that the engine horizontal positioning can influence the mission fuel burn by up to 6.4%. According to the extensive literature review that has been done, this study can be regarded as the first open literature engine position-NPF parametric study using CFD. Even though no correlations were extracted for other conditions; a deployed high-lift wing configuration was also studied in detail, defining the main aerodynamics effects of the engine integration at high angle of attack. A topological study of the high-lift installation vortices is presented in this work and it can be considered the first in the open literature. It should be pointed out that extensive research is currently underway to correctly evaluate the high-lift aerodynamic using CFD. The Propulsive System Integration (PSI) in high-lift conditions is adding flow features to an already demanding problem, making it a real challenge for the numerical methods. Nevertheless the additional effects of a nacelle chine on the maximum lift were also evaluated. The main outcomes of this PhD research were: a coupled performance modelling tool able to handle the effects of engine-airframe integration as a function of geometry and attitude, and a topological study of the high-lift installation vortices. During the course of the work, this research was successfully suggested as an extra activity for the European NEWAC project (New Aero Engine Core Concepts), and resulted in a new deliverable for that project.

Birds and bats employ a variety of advanced wing motions in the efficient production of thrust. The purpose of this thesis is to quantify the benefit of these advanced wing motions, determine the optimal theoretical wing kinematics for a given flight condition, and to develop a methodology for applying the results in the optimal design of flapping-wing aircraft (ornithopters). To this end, a medium-fidelity, combined aero-structural model has been developed that is capable of simulating the advanced kinematics seen in bird flight, as well as the highly non-linear structural deformations typical of high-aspect ratio wings. Five unique methods of thrust production observed in natural species have been isolated, quantified and thoroughly investigated for their dependence on Reynolds number, airfoil selection, frequency, amplitude and relative phasing. A gradient-based optimization algorithm has been employed to determined the wing kinematics that result in the minimum required power for a generalized aircraft or species in any given flight condition. In addition to the theoretical work, with the help of an extended team, the methodology was applied to the design and construction of the world's first successful human-powered ornithopter. The Snowbird Human-Powered Ornithopter, is used as an example aircraft to show how additional design constraints can pose limits on the optimal kinematics. The results show significant trends that give insight into the kinematic operation of natural species. The general result is that additional complexity, whether it be larger twisting deformations or advanced wing-folding mechanisms, allows for the possibility of more efficient flight. At its theoretical optimum, the efficiency of flapping-wings exceeds that of current rotors and propellers, although these efficiencies are quite difficult to achieve in practice.

Birds and bats employ a variety of advanced wing motions in the efficient production of thrust. The purpose of this thesis is to quantify the benefit of these advanced wing motions, determine the optimal theoretical wing kinematics for a given flight condition, and to develop a methodology for applying the results in the optimal design of flapping-wing aircraft (ornithopters). To this end, a medium-fidelity, combined aero-structural model has been developed that is capable of simulating the advanced kinematics seen in bird flight, as well as the highly non-linear structural deformations typical of high-aspect ratio wings. Five unique methods of thrust production observed in natural species have been isolated, quantified and thoroughly investigated for their dependence on Reynolds number, airfoil selection, frequency, amplitude and relative phasing. A gradient-based optimization algorithm has been employed to determined the wing kinematics that result in the minimum required power for a generalized aircraft or species in any given flight condition. In addition to the theoretical work, with the help of an extended team, the methodology was applied to the design and construction of the world's first successful human-powered ornithopter. The Snowbird Human-Powered Ornithopter, is used as an example aircraft to show how additional design constraints can pose limits on the optimal kinematics. The results show significant trends that give insight into the kinematic operation of natural species. The general result is that additional complexity, whether it be larger twisting deformations or advanced wing-folding mechanisms, allows for the possibility of more efficient flight. At its theoretical optimum, the efficiency of flapping-wings exceeds that of current rotors and propellers, although these efficiencies are quite difficult to achieve in practice.

Este trabalho apresenta um procedimento simples e robusto para avaliação e análise do desempenho de instalações propulsoras de embarcações de pesca operando no Estado do Espírito Santo. O objetivo é verificar se as embarcações estão operando com baixa eficiência propulsiva, portanto, com consumo de combustível maior do que o desejável, e identificar possíveis ações de melhoria desta eficiência que poderiam ser implementadas com baixo custo. Tendo em vista as condições artesanais da construção e operação das embarcações de pesca no Espírito Santo, o procedimento foi desenvolvido de forma a poder ser aplicado com relativa facilidade pelas pessoas envolvidas na atividade pesqueira artesanal neste Estado. O procedimento envolve o levantamento de informações sobre a geometria do casco, as características da instalação propulsora existente e sobre alguns aspectos do desempenho das embarcações. Para este levantamento são propostos métodos e técnicas suficientemente simples para que possam ser facilmente aplicadas nos estaleiros. Com base nestes dados levantados implementou-se no procedimento proposto diversos algoritmos que permitem determinar o coeficiente propulsivo das embarcações de pesca e a perda de energia útil que ocorre tanto no motor, por conta de condições inadequadas de operação, quanto na transmissão desta energia para o hélice. Para a validação dos métodos e técnicas envolvidos neste procedimento e dos resultados fornecidos por ele, o mesmo foi aplicado em alguns cascos de embarcações de pesca disponíveis em estaleiros visitados no Estado. Para uma destas embarcações o procedimento foi aplicado em sua totalidade, permitindo identificar a baixa eficiência da instalação propulsora e algumas das possíveis causas que provocam a perda de energia útil no sistema. Finalmente, com base nesse exemplo de aplicação, um estudo de possíveis melhorias na instalação propulsora é apresentado ao final do trabalho. / This study presents a simple and robust procedure to the performance evaluation and analysis of propulsion systems of fishing vessels operating in the Espírito Santo State. The goals are to verify if the vessels are operating at low propulsive efficiency, therefore, with a fuel consumption greater than desired, and identify possible actions to improve this efficiency that could be implemented at a low cost. The procedure was developed in view of the conditions of construction and operation of fishing vessels in Espírito Santo, so that it can be applied with relative ease by the people involved in artisanal fishing activity in this State. The procedure involves the collection of information about the geometry of hull, the characteristics of existing propulsion plant and some aspects of the performance of the vessels. To do this data collection, this paper proposes methods and techniques simple enough to be easily applied in the shipyards. Based on these data collected, the proposed procedure implemented several algorithms to ascertain the propulsive coefficient of fishing vessels and the energy loss which may occur due to the inadequate operating conditions and in energy transmission to propeller. It has been applied in some fishing vessel hulls available on sites visited to validate the methods and techniques involved in this procedure and its results. The procedure was applied, in its entirety, in one of this vessels, that allowed to identify the low efficiency of the propulsion plant and some of the possible causes that justify this energy loss in the system. Finally, based on this application example, a study of possible improvements in propulsion plant is presented at the end of work.

Wireless applications have become a common part of daily life. Whether it is mobile phones, the Wi-Fi router at home, the keycard which has replaced the car key, a radio frequency identification access system to a building or a Bluetooth headset for your computer or phone, the means of modern wireless data exchange is an omnipresent technology. In sports, the market is in its infancy for wireless, technical applications or gadgets. Only heart rate monitors and GPS watches are currently used by recreational athletes. Even though most of the larger sports equipment companies regularly launch new products related to sports performance monitoring and mobile phone technology, product innovation leaps are rare.In this work the design of a wireless sports performance measurement platform is presented. Using the example of kayaking, this platform is configured as a paddle performance measuring system, the Kayak XL System, which can monitor propulsive paddle force, paddle kinematics and boat velocity, interalia. A common mobile phone platform has been chosen as the user interface for this system. The design approach focussing on user requests, demands and expectations in combination with the process of iterative technical development are unveiled in this thesis. An evaluation of the system is presented and the work is finalised with an overview of further systems which have been designed based on the developed measurement platform. The Kayak XL System is a flexible system designed to be mounted onto any standard kayak paddle and installed in any competition kayak. Versatility, unobtrusiveness and usability were major design concerns. The developed system consists of four modules plus a software which has been designed for Android mobile phones. The phone communicates with each of the four modules trough Bluetooth radio. These four modules are also referred to as nodes and have specific measurement purposes. Two nodes have been designed to measure paddle force and kinematics, one node has the purpose to measure foot stretcher force and boat motion data, and the fourth node enables a more convenient method of calibrating paddle force measurement. The fourth node is therefore only needed prior to performance data acquisition. Results show that paddle and foot stretcher force can be measured with a resolution below 1N after calibration. Installing the paddle nodes on a previously configured paddle without repeated calibration is facilitated with the compromise of a doubled error margin. The default sampling frequency is set to 100 Hz and can, like all system parameters, be configured on the mobile phone. Real-time computation of complex performance parameters is only limited by the phone CPU. The system adds twice 109 g to the paddle and approximately 850 g to the kayak, excluding the mass of the mobile phone / <p>QC 20120827</p>

La vectorisation de la poussée d'une tuyère propulsive supersonique axisymétrique est étudiée par le biais d'une injection fluidique transversale dans sa partie divergente. Cette étude menée dans le cadre du programme PERSEUS du CNES a été motivée par la recherche d'une solution alternative au pilotage conventionnel de la poussée par actionneurs mécaniques. Le travail de la thèse, tout en s'appuyant sur des approches expérimentale et numérique, comprend essentiellement une large étude paramétrique concernant principalement la position de l'injection, la forme de la tuyère, la nature et le débit du fluide injecté. L'analyse des résultats montre que pour certaines configurations optimales, des angles de déviation pertinents peuvent être obtenus pour des taux d'injections modérés. L'analyse numérique étendue aux écoulements chauds multi-espèces, plus proches des applications réelles, a montré que la vectorisation fluidique reste très performante lors de l'injection de produits de combustion dans le divergent.

La vectorisation de la poussée d'une tuyère propulsive supersonique axisymétrique est étudiée par le biais d'une injection fluidique transversale dans sa partie divergente. Cette étude menée dans le cadre du programme PERSEUS du CNES a été motivée par la recherche d'une solution alternative au pilotage conventionnel de la poussée par actionneurs mécaniques. Le travail de la thèse, tout en s'appuyant sur des approches expérimentale et numérique, comprend essentiellement une large étude paramétrique concernant principalement la position de l'injection, la forme de la tuyère, la nature et le débit du fluide injecté. L'analyse des résultats montre que pour certaines configurations optimales, des angles de déviation pertinents peuvent être obtenus pour des taux d'injections modérés. L'analyse numérique étendue aux écoulements chauds multi-espèces, plus proches des applications réelles, a montré que la vectorisation fluidique reste très performante lors de l'injection de produits de combustion dans le divergent. / Secondary injection into the divergent section of a supersonic rocket nozzle is investigated for the fluidic thrust vectoring effects. The study was conducted in the framework of CNES PERSEUS program and was motivated by the need for an alternative vectoring solution aimed for a small space launcher. The thesis work, based on the combined experimental and numerical approaches, essentially comprises of a wide parametric study mainly concerning the position of the injection, the shape of the primary and injection nozzles, flow regime, gas thermophysical properties and injected fluid mass-flow-rate. The analysis shows that for some optimal configurations, pertinent deflection angles can be obtained using the moderate injection rates. Furthermore, the study was extended to the hot flow multi-species investigation, simulating a case closer to the real applications. This numerical analysis indicated that the fluidic vectoring method remains effective with injection of combustion products into the divergent section of a propulsive rocket nozzle.

There is currently a great deal of interest within the aviation community in the design of small, slow-flying but manoeuvrable uninhabited vehicles for reconnaissance, surveillance, and search and rescue operations in urban environments. Inspired by observation of birds, insects, fish and cetaceans, flapping wings are being actively studied in the hope that they may provide greater propulsive efficiencies than propellers and rotors at low Reynolds numbers for such Micro-Air Vehicles (MAVs). Researchers have posited the Strouhal number (combining flapping frequency, amplitude and forward speed) as the parameter controlling flapping wing aerodynamics in cruising flight, although there is conflicting evidence. This thesis explores the effect of flapping frequency and amplitude on forces and wake structures, as well as physical mechanisms leading to optimum propulsive efficiency. Two-dimensional rigid airfoils are considered at Reynolds number 2,000 ??? 40,000. A compressible Navier-Stokes simulation is combined with numerical and analytical potential flow techniques to isolate and evaluate the effect of viscosity, leading and trailing edge vortex separation, and wake vortex dynamics. The wake structures of a plunging airfoil are shown to be sensitive to the flapping frequency independent of the Strouhal number. For a given frequency, the wake of the airfoil exhibits ???vortex lock-in??? as the amplitude of motion is increased, in a manner analogous to an oscillating circular cylinder. This is caused by interaction between the flapping frequency and the ???bluff-body??? vortex shedding frequency apparent even for streamlined airfoils at low Reynolds number. The thrust and propulsive efficiency of a plunging airfoil are also shown to be sensitive to the flapping frequency independent of Strouhal number. This dependence is the result of vortex shedding from the leading edge, and an interaction between the flapping frequency and the time for vortex formation, separation and convection over the airfoil surface. The observed propulsive efficiency peak for a pitching and plunging airfoil is shown to be the result of leading edge vortex shedding at low flapping frequencies (low Strouhal numbers), and high power requirements at large flapping amplitudes (high Strouhal numbers). The efficiency peak is governed by flapping frequency and amplitude separately, rather than the Strouhal number directly.

There is currently a great deal of interest within the aviation community in the design of small, slow-flying but manoeuvrable uninhabited vehicles for reconnaissance, surveillance, and search and rescue operations in urban environments. Inspired by observation of birds, insects, fish and cetaceans, flapping wings are being actively studied in the hope that they may provide greater propulsive efficiencies than propellers and rotors at low Reynolds numbers for such Micro-Air Vehicles (MAVs). Researchers have posited the Strouhal number (combining flapping frequency, amplitude and forward speed) as the parameter controlling flapping wing aerodynamics in cruising flight, although there is conflicting evidence. This thesis explores the effect of flapping frequency and amplitude on forces and wake structures, as well as physical mechanisms leading to optimum propulsive efficiency. Two-dimensional rigid airfoils are considered at Reynolds number 2,000 ??? 40,000. A compressible Navier-Stokes simulation is combined with numerical and analytical potential flow techniques to isolate and evaluate the effect of viscosity, leading and trailing edge vortex separation, and wake vortex dynamics. The wake structures of a plunging airfoil are shown to be sensitive to the flapping frequency independent of the Strouhal number. For a given frequency, the wake of the airfoil exhibits ???vortex lock-in??? as the amplitude of motion is increased, in a manner analogous to an oscillating circular cylinder. This is caused by interaction between the flapping frequency and the ???bluff-body??? vortex shedding frequency apparent even for streamlined airfoils at low Reynolds number. The thrust and propulsive efficiency of a plunging airfoil are also shown to be sensitive to the flapping frequency independent of Strouhal number. This dependence is the result of vortex shedding from the leading edge, and an interaction between the flapping frequency and the time for vortex formation, separation and convection over the airfoil surface. The observed propulsive efficiency peak for a pitching and plunging airfoil is shown to be the result of leading edge vortex shedding at low flapping frequencies (low Strouhal numbers), and high power requirements at large flapping amplitudes (high Strouhal numbers). The efficiency peak is governed by flapping frequency and amplitude separately, rather than the Strouhal number directly.

Stability is an important criterion in the design and performance of launch vehicles. Present day launch vehicles have become more and more flexible due to the constraints of weight reduction, necessarily imposed for enhanced performance of the vehicle. Due to higher flexibility, the launch vehicle stability becomes a concern. Instability in the launch vehicles has been noticed due to three major sources: thrust, aerodynamic forces and combustion induced instabilities. Instability in the launch vehicles may pose problem to the structural integrity leading to structural failure or it may lead to the deviation in the trajectory of the vehicle. Several structural failures of launch vehicles due to instabilities have been reported in the literature. The prediction of the structural response due to various excitations such as thrust and aerodynamic loading is essential to identify any failure scenarios and to limit the vibrations transmitted to the payload. Therefore, determination of dynamic and stability characteristics of a launch vehicle under the influence of different parameters, is of vital importance. Disciplines such as, flight mechanics (dynamics), structural dynamics, aerodynamics, propulsion, guidance and control are closely related in the design and analysis of launch vehicles. Typically, flight mechanics, guidance and control problems consider a rigid vehicle for modeling and simulation purposes. The disciplines of structural dynamics and aeroelasticity consider a flexible vehicle. In order to bring in the effect of flexibility on the flight dynamics of the launch vehicle, structural dynamics and aeroelasticity aspects need to be effected. The preliminary design of a new launch vehicle requires inputs from different disciplines and parametric studies are required to finalise the vehicle configuration. The study of the effect of different parameters on the dynamics and stability of launch vehicles is required. In this context, there is a need to develop an integrated approach that provides tools for the design and analysis of a launch vehicle. The availability of integrated modeling and simulation tools will reduce the requirement of costly prototype development and testing. In the present thesis, an attempt has been made to develop a numerical tool to conduct parametric studies for launch vehicle dynamics and stability. The developed tool is suitable for prediction of onset of instabilities under the influence of different parameters. The approach developed in this thesis is also well suited for specialized analysis of problems involving vertical launch, stage separation, engine shutdown and internal stress wave propagation related to structural integrity. Stability problems due to thrust and the aerodynamic forces (aeroelastic stability) in the launch vehicles/ missiles have been reported in the literature. Most of these works have modeled the vehicle as a beam or by using discrete degrees of freedom. In these works, the effect of thrust or aerodynamic forces on the flexible body modes is investigated and it is shown that the instability may occur in one of the bending modes due to change in the parameters such as thrust or aerodynamic forces. Traditionally, the dynamic characteristics are obtained in a body-fixed coordinate system, whereas the prediction of trajectory (rigid body dynamics) is carried out in an inertial frame of reference. Only few works have addressed the coupling of the rigid body motion and the flexible body dynamics of a vehicle. But these works also, do not consider the total derivative of displacements with respect to an inertial frame of reference. When the integrated equations of motion are derived in an inertial frame of reference, the rigid body motion and the elastic displacements are highly coupled. In this thesis, the rigid body motion and the flexible body dynamics is studied in an inertial frame of reference. The flexible body dynamics of the moving vehicle is studied in an inertial frame of reference, including velocity induced curvature effects, which have not been considered so far in the published literature. A detailed mechanics based model is developed to analyze the problem of structural instabilities in launch vehicles. Coupling among the rigid-body modes, the longitudinal vibrational modes and the transverse vibrational modes due to asymmetric lifting-body cross-section are considered. The model also incorporates the effects of aerodynamic forces and the propulsive thrust of the vehicle. The propulsive thrust is considered as a follower force. The model is one-dimensional, and it can be employed to idealized slender vehicles with complex shapes. The governing differential equations along with the boundary conditions are derived using Extended Hamilton’s principle. Subsequently, the modeling of the propulsive thrust and the aerodynamic forces are included in the formulation. In the literature, the propulsive thrust has generally been modeled as a follower force applied at the nozzle end. Few of the works in the literature have modeled the combustion process in the solid rocket motor and the liquid propellant engine in detail. This is required to understand the combustion induced instabilities. In the present thesis, the propulsive thrust is considered as a follower force and few of the combustion parameters affecting the thrust are considered. In the literature, the modeling of the aerodynamic forces acting on a launch vehicle has been carried out using general purpose computational fluid dynamics (CFD) codes or by using empirical methods. CFD codes are used to obtain the pressure and the shear stress distribution on the vehicle surface by the solution of Navier Stokes/ Euler equations. The empirical methods have been used to obtain the distributed aerodynamic forces acting on the vehicle. The aerodynamic forces are expressed in terms of distributed aerodynamic coefficients. In the present work, the modeling of the aerodynamic forces has been carried out in two different ways: using a CFD package and by using empirical methods. The stability of a system can be studied by determining the system response with time. Eigenvalue analysis is another tool to investigate the stability of a linear system. To study the stability characteristics of the system using eigenvalue analysis, a computational framework has been developed. For this purpose, the finite element discretization of the system is carried out. Further to that, two different methods are utilized for finite element discretization of the vehicle structure: Fourier Transform based Spectral Finite Element method (SFEM) and an hp Finite Element method (FEM). The conventional FEM is a versatile tool for modeling complicated structures and to obtain the solution of the system of equations for a variety of forcing functions. The SFEM is more suitable for obtaining the solution for simple 1D and 2D structures subjected to shock and transient loads, having high frequency content. In this thesis, the spectral finite element model is developed for a vehicle subjected to the propulsive thrust and the aerodynamic forces. Prediction of instability using SFEM, means solving a nonlinear eigenvalue problem. Standard computer codes or routines are not available for solving a nonlinear eigenvalue problem. A computer code has been written to solve the nonlinear eigenvalue problem using one of the algorithms available in the literature. An hp finite element model is also developed for launch vehicle. The finite element stiffness and damping matrices due to the thrust, the aerodynamic forces and the rigid body velocity and acceleration are derived using Lagrange’s equations of motion. A standard linear eigenvalue problem and a polynomial eigenvalue problem is formulated for determination of instability regimes of the vehicle. It is important to understand the influence of different parameters such as thrust, velocity, angle of attack etc. on the stability of a launch vehicle. Parametric studies are important during the preliminary design phase of a vehicle to identify the instability regimes. The design parameters can be changed to reduce the possibility of instabilities. Numerical simulations are carried out to determine the unstable regimes of a slender launch vehicle for propulsive thrust and velocity as the parameters, neglecting the aerodynamic forces. Comparison between the results based on a Fourier spectral finite element model and a hp finite element model are carried out. Phenomenon of static instability (divergence) and dynamic instability (flutter) are observed. Determination of mode shapes of the vehicle is important for deciding the placement of sensors and actuators on the vehicle. In this context, eigenvectors (mode shapes) for different end thrust and speed are analyzed. Further, numerical simulations are also carried out to determine the instabilities in a slender launch vehicle considering the combined effects of propulsive thrust, aerodynamic forces and mass variation. The finite element model simulation results for aeroelastic effects are compared with the published literature. Stability of a vehicle is analysed for velocity (free stream Mach number) as a parameter, at maximum propulsive thrust, including the effect of aerodynamic forces and mass variation. Phenomenon of static instability (divergence) and dynamic instability (flutter) are observed. With the increase in the Mach number, branching (splitting) and merging of the modes is observed. At higher Mach numbers, divergence and flutter are observed in different modes simultaneously. Numerical simulations are carried out for a typical nosecone launch vehicle configuration to analyse the aeroelastic stability at two different Mach numbers using empirical aerodynamic data. The phenomenon of flow separation and reattachment is observed at the cone-cylinder junction. The stability of a typical vehicle under propulsive thrust and aerodynamic forces is investigated using CFD derived aerodynamic data. The aerodynamic pressure and shear stress distribution for a launch vehicle are obtained from the CFD analysis. The effect of different parameters such as combustion chamber pressure, tip mass and slenderness ratio on the stability of a vehicle is studied. In the later part of the thesis, solution methodology for the time domain response for a coupled axial and transverse motion of a vehicle is developed. The axial responses (displacements and velocities) of a typical vehicle subjected to axial thrust are determined using direct integration of the equations of motion. The axial displacements due to two different thrust histories are compared. The axial velocities with time at different locations are determined. The time domain and the frequency domain responses for a representative vehicle subjected to a transverse shock force are determined using Spectral Finite Element method (SFEM). The system of equations for a coupled axial and transverse motion of a vehicle is developed. Numerical simulations are carried out to determine the coupled axial and transverse response of a vehicle subjected to axial and transverse forces. The coupling of rigid body motion with the elastic displacements is illustrated. The thesis is comprised of seven chapters. The first chapter gives a detailed introduction to launch vehicles and covers literature survey of launch vehicle dynamics and stability. The dynamics and stability related aspects of flexible structures are also discussed. In chapter 2, a detailed mathematical model of a slender launch vehicle is developed to analyze the problem of structural instabilities. Chapter 3 deals with the finite element discretization of the vehicle structure using two different methods: Fourier spectral finite element method and an hp finite element method. In chapters 4 and 5, numerical simulations are carried out to determine the instabilities in a slender launch vehicle considering the effects of propulsive thrust, aerodynamic forces and mass variation. In chapter 6, solution methodology for the time domain response for a coupled axial and transverse motion of a vehicle is developed. The last chapter gives the conclusions and the future scope of work. To summarize, this thesis is a comprehensive document, that not only describes some detailed mathematical models for launch vehicle stability studies, but also presents the effect of aerodynamic, propulsion and structural loads on the launch vehicle stability. Linear stability analysis of a representative vehicle is carried out for prediction of onset of the instabilities under the influence of different parameters such as velocity, thrust, combustion factors etc. The correlation between the stability analysis and the time domain response is established. In short, the matter presented in this thesis can serve as a useful design aide for those working in the launch vehicle design.