In this thesis the characteristics of pitch and roll damping coefficients of missile models are investigated by using Computational Fluid Dynamics (CFD) techniques. Experimental data of NACA0012 airfoil, Basic Finner (BF) and Modified Basic Finner (MBF) models are used for validation and verification studies. Numerical computations are performed from subsonic to supersonic flow regimes. Grid refinement and turbulence model selection studies are conducted before starting the dynamic motion simulations. Numerical method of dynamic motion simulation is validated with a 2D NACA0012 airfoil. After the validation of numerical method, forced-oscillation motion is given to the BF and MBF models. In order to get deeper understandings about the characteristics of dynamic pitching and rolling motions, parametric studies are performed. The amplitude and frequency of forced-oscillation motions are investigated one by one. The effects of angle of attacks are also investigated for both pitching and rolling motions. The results of CFD simulations are compared with experimental data obtained from different wind tunnel and free flight tests. It is seen from these comparisons that experimental and numerical results are in good agreement throughout the whole flow regime. In conclusion, the numerical method presented in this study is validated and can be used for the prediction of pitch and roll damping coefficient of any missile configurations.

Reliability is one of the most critical performance measures of guided missile systems. It is directly related to missile mission success. In order to have a high reliability value, reliability analysis should be carried out at all phases of the system design. Carrying out reliability analysis at all the phases of system design helps the designer to make reliability related design decisions in time and update the system design. In this study, reliability analysis process performed during the conceptual design phase of a Medium Range Anti-Tank Missile System Inertial Measurement Unit (IMU) was introduced. From the reliability requirement desired for the system, an expected IMU reliability value was derived by using reliability allocation methods. Then, reliability prediction for the IMU was calculated by using Relex Software. After that, allocated and predicted reliability values of the IMU were compared. It was seen that the predicted reliability value of the IMU did not meet the required reliability value. Therefore, reliability improvement analysis was carried out.

In this thesis, a simple moving boundary technique has been suggested, implemented and verified. The technique may be considered as a generalization of the well-known "ghost" cell approach for boundary condition implementation. According to the proposed idea, the moving body does not appear on the computational grid and is allowed to move over the grid. The impermeable wall boundary condition is enforced by assigning proper gasdynamic values at the grid nodes located inside the moving body close to its boundaries (ghost nodes). The reflection principle taking into account the velocity of the boundaries assigns values at the ghost nodes. The new method does not impose any particular restrictions on the geometry, deformation and law of motion of the moving body. / The developed technique is rather general and can be used with virtually any finite-volume or finite-difference scheme, since the modifications of the schemes themselves are not required. In the present study the proposed technique has been incorporated into a one-dimensional non-adaptive Euler code and a two-dimensional locally adaptive unstructured Euler code. / It is shown that the new approach is conservative with the order of approximation near the moving boundaries. To reduce the conservation error, it is beneficial to use the method in conjunction with local grid adaptation. / The technique is verified for a number of one and two dimensional test cases with analytical solutions. It is applied to the problem of supersonic inlet starting via variable geometry approach. At first, a classical starting technique of changing exit area by a moving wedge is numerically simulated. Then, the feasibility of some novel ideas such as a collapsing frontal body and "tractor-rocket" are explored.

A three-dimensional, parallel, finite volume solver which uses Roe&#039 / s upwind flux differencing scheme for spatial and Runge-Kutta explicit multistage time stepping scheme for temporal discretization on unstructured meshes is developed for the unsteady solution of external viscous flow around rotating bodies. The main aim of this study is to evaluate the aerodynamic dynamic stability derivative coefficients for rotating missile configurations. Arbitrary Lagrangian Eulerian (ALE) formulation is adapted to the solver for the simulation of the rotation of the body. Eigenvalues of the Euler equations in ALE form has been derived. Body rotation is simply performed by rotating the entire computational domain including the body of the projectile by means of rotation matrices. Spalart-Allmaras one-euqation turbulence model is implemented to the solver. The solver developed is first verified in 3-D for inviscid flow over two missile configurations. Then inviscid flow over a rotating missile is tested. Viscous flux computation algorithms and Spalarat-Allmaras turbulence model implementation are validated in 2-D by performing calculations for viscous flow over flat plate, NACA0012 airfoil and NLR 7301 airfoil with trailing edge flap. The ALE formulation is validated in 2-D on a rapidly pitching NACA0012 airfoil. Afterwards three-dimensional validation studies for viscous, laminar and turbulent flow calculations are performed on 3-D flat plate problem. At last, as a validation test case, unsteady laminar and turbulent viscous flow calculations over a spinning M910 projectile configuration are performed. Results are qualitatively in agreement with the analytical solutions, experimental measurements and previous studies for steady and unsteady flow calculations.

No description available.

This study has conducted an OLS-regression to examine the relationship between gasoline and crude oil prices in the Swedish petroleum market, focusing on potential asymmetry, where gasoline prices respond more quickly to increases in crude oil prices compared to decreases. Additionally, we examine the impact of individual petroleum stations' competitiveness on this asymmetry, applying the distance to the nearest station as a measure of competitiveness. To explore this relationship, we utilise two datasets: one comprises unique user-generated data for individual gas stations spanning the period from 2019 to 2022, while the other includes recommended prices covering the period from 2001 to 2020. The findings provide some evidence supporting the existence of asymmetry, indicating the presence of inefficiencies within the market. However, no evidence suggesting larger asymmetry concerning individual competitiveness was discovered. Overall, this research offers novel insights into the dynamics of the Swedish fuel market in recent years.

The initial stage of liquid phase sintering, involving liquid penetration, agglomerate separation, particle spreading and growth has been investigated in experiments using tungsten heavy alloys. The particle composites used were produced by hot isostatic pressing (HIP) of pure powder mixtures of W-Ni-Fe-(Co). By using different HIP temperatures, volume fractions of tungsten, alloying elements like Cobalt and Sulphur or excluding Iron from the matrix, liquid penetration, agglomerate separation and particle growth conditions were affected. The investigations were performed mainly under microgravity (sounding rockets or parabolic trajectories by airplanes) but at high tungsten particle fractions, short sintering times or at infiltration of solid pure tungsten, they were performed at normal gravity. The liquid penetration of the tungsten agglomerates is explained by initial wetting under non-equilibrium conditions, due to the reaction between the liquid matrix and the particles, and a decrease of interfacial energy. The dissolving of tungsten gives a pressure drop in the penetrating liquid and a driving force for the liquid movement by a suggested parabolic penetration model. For cold worked tungsten, a penetration theory was proposed, where an internal stress release in the penetrated tungsten grains creates space for the advancing liquid. The spreading of the tungsten agglomerates is explained by an interagglomerate melt swelling due to a Kirkendall effect. The liquid matrix undergoes a volume increase since the diffusion rates of Ni-Fe are higher than for W and initial concentration gradients of W and Ni, Fe exists. The suggested model by Kirkendall are also used for an analysis of the interaction behaviour between solid particles and a solidification front and inclusion behaviour in iron base alloys during teeming and deoxidation. The average tungsten particles size decrease initially since part of the tungsten particles is dissolved when the non-equilibrium matrix phase is melting. When equilibrium is reached, the tungsten particles grow in accordance with the Ostwald ripening process by an approximately 1/3 power law. Larger particle fraction of particles showed a higher growth rate, due to shorter diffusion distances between the particles. Cobalt, Sulphur and absence of iron in the matrix were found to increase the growth rate of the tungsten particles due to a higher surface tension between the solid tungsten particles and the matrix melt. / QC 20100528

Controlling unsteady combustion induced gas flow fluctuations and the resultant motor vibrations is a very significant step in rocket motor design. It occurs when the unsteady heat release due to combustion happens to feed the acoustic oscillations of the closed duct forming a feed-back system. The resultant vibrations concerned may even lead to total failure of the rocket system unless analysed and tested thoroughly. This thesis aims developing a linear numerical analysis method for the growth rate of instabilities and possible mode shape of a liquid-propelled chamber geometry. In particular, A 3-D Helmholtz code, utilizing Culicks spatial averaging linear iterative method, is developed to find the form of deformed mode shapes iteratively to obtain possible effects of heat source and impedance boundary conditions. The natural mode shape phase is solved through finite volume discretization and the open-source eigenvalue extractor, ARPACK, and its parallel implementation PARPACK. The iterative method is particularly used for analyzing the geometries with complex shapes and essentially for disturbances of small magnitudes to natural mode shapes. The developed tools are tested via two simple cases, a duct with inactive flame and a Rijke tube, used as validation cases for the code particularly with only boundary contribution and heat contribution respectively. A sample 2-D and 3-D liquid-propelled combustion chamber is also analysed with heat sources. After comparing with the expected values, it is eventually proved that the method should be only used for determining the modes instability analysis, as to whether it keeps vibrating or decays. The methodology described can be used as a preliminary design tool for the design of liquid-propellant rocket engine combustors, rapidly revealing only the onset of instabilities.

Characteristics of impinging-jet atomizers in the context of application in liquid propulsion systems are studied in this thesis. A review of past studies on impinging jets revealed the necessity of a correlation in terms of injector parameters for predicting Sauter Mean Diameter (SMD) of a spray. So, an experimental study of atomization in doublet and triplet impinging jet injectors is conducted using water as the stimulant? The major injector parameters considered are orifice diameter, impingement angle and jet velocity. Relative influences of these parameters are explained in terms of a single parameter, specific normal momentum. SMD of the spray reduces as specific normal momentum is increased. A universal expression between non-dimensional SMD and specific normal momentum is obtained, which satisfactorily predicts SMD in doublets as well as triplets. Noting that practical impinging injectors are likely to have skewness (partial impingement), the study is extended to understand the behavior of such jets. In perfectly impinging doublet, a high aspect ratio ellipse-like mass distribution pattern is obtained with major axis normal to the plane of two jets whereas in skewed jets the major axis turns from its normal position. A simple correlation is obtained, which shows that this angle of turn is a function of skewness fraction and impingement angle only and is independent of injection velocity. Experimental data from both mass distribution and photographic technique validate this prediction. SMD is found to decrease as skewness is increased. This may be the combined effect of shearing of liquid sheet at the point of impingement and more sheet elongation. Hence, skewness turns out to be an important parameter in controlling drop size.

Made available in DSpace on 2014-12-17T14:55:41Z (GMT). No. of bitstreams: 1 GlaubertoLAA.pdf: 3487440 bytes, checksum: bcd5cef1c854f4d01f9c73419e1d7d42 (MD5) Previous issue date: 2009-11-20 / Global Positioning System, or simply GPS, it is a radionavigation system developed by United States for military applications, but it becames very useful for civilian using. In the last decades Brazil has developed sounding rockets and today many projects to build micro and nanosatellites has appeared. This kind of vehicles named spacecrafts or high dynamic vehicles, can use GPS for its autonome location and trajectories controls. Despite of a huge number of GPS receivers available for civilian applications, they cannot used in high dynamic vehicles due environmental issues (vibrations, temperatures, etc.) or imposed dynamic working limits. Only a few nations have the technology to build GPS receivers for spacecrafts or high dynamic vehicles is available and they imposes rules who difficult the access to this receivers. This project intends to build a GPS receiver, to install them in a payload of a sounding rocket and data collecting to verify its correct operation when at the flight conditions. The inner software to this receiver was available in source code and it was tested in a software development platform named GPS Architect. Many organizations cooperated to support this project: AEB, UFRN, IAE, INPE e CLBI. After many phases: defining working conditions, choice and searching electronic, the making of the printed boards, assembling and assembling tests; the receiver was installed in a VS30 sounding rocket launched at Centro de Lan?amento da Barreira do Inferno in Natal/RN. Despite of the fact the locations data from the receiver were collected only the first 70 seconds of flight, this data confirms the correct operation of the receiver by the comparison between its positioning data and the the trajectory data from CLBI s tracking radar named ADOUR / O Sistema de Posicionamento Global, conhecido mundialmente pala sigla GPS, ? um sistema de radionavega??o constru?do pelos norte-americanos com inten??es militares, mas que encontraram, com o passar do tempo, muitas aplica??es de uso civil. No Brasil, al?m do desenvolvimento de foguetes de sondagem, come?am a aparecer projetos de constru??o de micro e nanosat?lites. Estes ve?culos denominados espaciais ou de alta din?mica podem, quando em voo, usufruir do sistema GPS para localiza??o aut?noma e verifica??o/controle das suas trajet?rias. Apesar da enorme disponibilidade de receptores GPS no mercado civil, estes n?o podem ser utilizados em ve?culos de alta din?mica, seja por quest?es ambientais (vibra??es, temperaturas elevadas, etc.) ou por prote??o l?gica (via software). Os receptores para uso em ve?culos de alta din?mica, ou ve?culos espaciais, fazem parte de uma tecnologia restrita a poucos pa?ses, que estabelecem regras muito r?gidas para suas aquisi??es. O presente projeto objetiva construir e validar funcionamento b?sico deste receptor ao instal?-lo num foguete de sondagem e coleta de dados em voo. O software a ser utilizado no receptor j? estava dispon?vel em c?digo fonte e testado em uma plataforma de desenvolvimento denominada GPS Architect. V?rios organismos cooperaram para realiza??o projeto: AEB, UFRN, IAE, INPE e CLBI. Ap?s v?rios passos para realiza??o do projeto: defini??o das condi??es de funcionamento, escolha e aquisi??o dos componentes eletr?nicos, fabrica??o das placas de circuito impresso, montagem e testes de integra??o; o mesmo foi instalado num foguete de sondagem VS30 lan?ado a partir do Centro de Lan?amento da Barreira do Inferno em Natal/RN. Apesar da coleta parcial dos dados do receptor, por falha t?cnica do sistema de telemetria do foguete, os resultados obtidos foram suficientes para validar o funcionamento do receptor a partir da compara??o entre os dados de trajetografia fornecidos pelo receptor GPS e o radar de trajetografia do CLBI conhecido como Radar ADOUR