Supersonic regime of the Hall-magnetohydrodynamics resistive tearing instability

An earlier analysis of the Hall-magnetohydrodynamics (MHD) tearing instability [E. Ahedo and J. J. Ramos, Plasma Phys. Controlled Fusion 51, 055018 (2009)] is extended to cover the regime where the growth rate becomes comparable or exceeds the sound frequency. Like in the previous subsonic work, a resistive, two-fluid Hall-MHD model with massless electrons and zero-Larmor-radius ions is adopted and a linear stability analysis about a force-free equilibrium in slab geometry is carried out. A salient feature of this supersonic regime is that the mode eigenfunctions become intrinsically complex, but the growth rate remains purely real. Even more interestingly, the dispersion relation remains of the same form as in the subsonic regime for any value of the instability Mach number, provided only that the ion skin depth is sufficiently small for the mode ion inertial layer width to be smaller than the macroscopic lengths, a generous bound that scales like a positive power of the Lundquist number

Similarities between ion waves in plasmas and gravity waves in incompressible fluids

Some similarities between ion waves in plasmas and gravity waves in incompressible fluids are investigated. It is shown that for zero ion temperature the ion-wave dispersion relation is similar to that of gravity waves in a stratified liquid between rigid, horizontal walls; for large wavelength the ion waves behave as the surface gravity waves of shallow-water theory. The general character of the pattern of ion waves arising in steady plasma flows is analyzed for arbitrary ion temperature, wavelength, and acoustic mach number (which is based on the ion-acoustic speed), and is compared to the pattern of surface gravity waves in steady water flows when surface tension is taken into account.

End effect in Langimur probe reponse under ionospheric satellite conditions

A theory is presented for an end effect in the current response of a highly negative, cylindrical Langmuir probe in a collisionless plasma flow. Under conditions where the ratio of probe radius to debye length is small and the ion-acoustic Mach number is large, the current exhibits a strong peak when the probe axis is brought into alignment with the flow direction. Closed formulas are given for the maximum and angular half-width of the peak, and universal graphical results are presented for the entire peak structure. The theory shows very good agreement with experimental data. The use of the end effect for diagnostic purposes, in particular, for the determination of the ion temperature, is discussed.

A new basic effect in retarding potential analyzers

The Retarding Potential Analyzer (RPA) is the standard instrument for in situ measurement of ion temperature and other ionospheric parameters. The fraction of incoming ions rejected by a RPA produces perturbations that reach well ahead of a thin Debye sheath, a feature common to all collisionless, hypersonic flows past ion-rejecting bodies. This phenomenon is here found to result in a correction to Whipple’s classical law for the current characteristic of an ideal RPA sheath thin; inverse ram ion Mach number M-1, and ram angle of RPA aperture u, small or moderately small.

Advances in global instability computations: from incompressible to hypersonic flow

Esta tesis constituye un gran avance en el conocimiento del estudio y análisis de inestabilidades hidrodinámicas desde un punto de vista físico y teórico, como consecuencia de haber desarrollado innovadoras técnicas para la resolución computacional eficiente y precisa de la parte principal del espectro correspondiente a los problemas de autovalores (EVP) multidimensionales que gobiernan la inestabilidad de flujos con dos o tres direcciones espaciales inhomogéneas, denominados problemas de estabilidad global lineal. En el contexto del trabajo de desarrollo de herramientas computacionales presentado en la tesis, la discretización mediante métodos de diferencias finitas estables de alto orden de los EVP bidimensionales y tridimensionales que se derivan de las ecuaciones de Navier-Stokes linealizadas sobre flujos con dos o tres direcciones espaciales inhomogéneas, ha permitido una aceleración de cuatro órdenes de magnitud en su resolución. Esta mejora de eficiencia numérica se ha conseguido gracias al hecho de que usando estos esquemas de diferencias finitas, técnicas eficientes de resolución de problemas lineales son utilizables, explotando el alto nivel de dispersión o alto número de elementos nulos en las matrices involucradas en los problemas tratados. Como más notable consecuencia cabe destacar que la resolución de EVPs multidimensionales de inestabilidad global, que hasta la fecha necesitaban de superordenadores, se ha podido realizar en ordenadores de sobremesa. Además de la solución de problemas de estabilidad global lineal, el mencionado desarrollo numérico facilitó la extensión de las ecuaciones de estabilidad parabolizadas (PSE) lineales y no lineales para analizar la inestabilidad de flujos que dependen fuertemente en dos direcciones espaciales y suavemente en la tercera con las ecuaciones de estabilidad parabolizadas tridimensionales (PSE-3D). Precisamente la capacidad de extensión del novedoso algoritmo PSE-3D para el estudio de interacciones no lineales de los modos de estabilidad, desarrollado íntegramente en esta tesis, permite la predicción de transición en flujos complejos de gran interés industrial y por lo tanto extiende el concepto clásico de PSE, el cuál ha sido empleado exitosamente durante las pasadas tres décadas en el mismo contexto para problemas de capa límite bidimensional. Típicos ejemplos de flujos incompresibles se han analizado en este trabajo sin la necesidad de recurrir a restrictivas presuposiciones usadas en el pasado. Se han estudiado problemas vorticales como es el caso de un vórtice aislado o sistemas de vórtices simulando la estela de alas, en los que la homogeneidad axial no se impone y así se puede considerar la difusión viscosa del flujo. Además, se ha estudiado el chorro giratorio turbulento, cuya inestabilidad se utiliza para mejorar las características de funcionamiento de combustores. En la tesis se abarcan adicionalmente problemas de flujos compresibles. Se presenta el estudio de inestabilidad de flujos de borde de ataque a diferentes velocidades de vuelo. También se analiza la estela formada por un elemento rugoso aislado en capa límite supersónica e hipersónica, mostrando excelentes comparaciones con resultados obtenidos mediante simulación numérica directa. Finalmente, nuevas inestabilidades se han identificado en el flujo hipersónico a Mach 7 alrededor de un cono elíptico que modela el vehículo de pruebas en vuelo HIFiRE-5. Los resultados comparan favorablemente con experimentos en vuelo, lo que subraya aún más el potencial de las metodologías de análisis de estabilidad desarrolladas en esta tesis. ABSTRACT The present thesis constitutes a step forward in advancing the frontiers of knowledge of fluid flow instability from a physical point of view, as a consequence of having been successful in developing groundbreaking methodologies for the efficient and accurate computation of the leading part of the spectrum pertinent to multi-dimensional eigenvalue problems (EVP) governing instability of flows with two or three inhomogeneous spatial directions. In the context of the numerical work presented in this thesis, the discretization of the spatial operator resulting from linearization of the Navier-Stokes equations around flows with two or three inhomogeneous spatial directions by variable-high-order stable finite-difference methods has permitted a speedup of four orders of magnitude in the solution of the corresponding two- and three-dimensional EVPs. This improvement of numerical performance has been achieved thanks to the high-sparsity level offered by the high-order finite-difference schemes employed for the discretization of the operators. This permitted use of efficient sparse linear algebra techniques without sacrificing accuracy and, consequently, solutions being obtained on typical workstations, as opposed to the previously employed supercomputers. Besides solution of the two- and three-dimensional EVPs of global linear instability, this development paved the way for the extension of the (linear and nonlinear) Parabolized Stability Equations (PSE) to analyze instability of flows which depend in a strongly-coupled inhomogeneous manner on two spatial directions and weakly on the third. Precisely the extensibility of the novel PSE-3D algorithm developed in the framework of the present thesis to study nonlinear flow instability permits transition prediction in flows of industrial interest, thus extending the classic PSE concept which has been successfully employed in the same context to boundary-layer type of flows over the last three decades. Typical examples of incompressible flows, the instability of which was analyzed in the present thesis without the need to resort to the restrictive assumptions used in the past, range from isolated vortices, and systems thereof, in which axial homogeneity is relaxed to consider viscous diffusion, as well as turbulent swirling jets, the instability of which is exploited in order to improve flame-holding properties of combustors. The instability of compressible subsonic and supersonic leading edge flows has been solved, and the wake of an isolated roughness element in a supersonic and hypersonic boundary-layer has also been analyzed with respect to its instability: excellent agreement with direct numerical simulation results has been obtained in all cases. Finally, instability analysis of Mach number 7 ow around an elliptic cone modeling the HIFiRE-5 flight test vehicle has unraveled flow instabilities near the minor-axis centerline, results comparing favorably with flight test predictions.

Laminar-turbulent transition induced by a discrete roughness element in a supersonic boundary layer

The linear instability and breakdown to turbulence induced by an isolated roughness element in a boundary layer at Mach 2:5, over an isothermal flat plate with laminar adiabatic wall temperature, have been analysed by means of direct numerical simulations, aided by spatial BiGlobal and three-dimensional parabolized (PSE-3D) stability analyses. It is important to understand transition in this flow regime since the process can be slower than in incompressible flow and is crucial to prediction of local heat loads on next-generation flight vehicles. The results show that the roughness element, with a height of the order of the boundary layer displacement thickness, generates a highly unstable wake, which is composed of a low-velocity streak surrounded by a three-dimensional high-shear layer and is able to sustain the rapid growth of a number of instability modes. The most unstable of these modes are associated with varicose or sinuous deformations of the low-velocity streak; they are a consequence of the instability developing in the three-dimensional shear layer as a whole (the varicose mode) or in the lateral shear layers (the sinuous mode). The most unstable wake mode is of the varicose type and grows on average 17% faster tan the most unstable sinuous mode and 30 times faster than the most unstable boundary layer mode occurring in the absence of a roughness element. Due to the high growthrates registered in the presence of the roughness element, an amplification factor of N D 9 is reached within 50 roughness heights from the roughness trailing edge. The independently performed Navier–Stokes, spatial BiGlobal and PSE-3D stability results are in excellent agreement with each other, validating the use of simplified theories for roughness-induced transition involving wake instabilities. Following the linear stages of the laminar–turbulent transition process, the roll-up of the three-dimensional shear layer leads to the formation of a wedge of turbulence, which spreads laterally at a rate similar to that observed in the case of compressible turbulent spots for the same Mach number.

Physics of Vibrating Airfoils at Low Reduced Frequency

The unsteady aerodynamics of low pressure turbine vibrating airfoils in flap mode is studied in detail using a frequency domain linearized Navier-Stokes solver. Both the travelling-wave and influence coefficient formulations of the problem are used to highlight key aspects of the physics and understand different trends such as the effect of reduced frequency and Mach number. The study is focused in the low-reduced frequency regime which is of paramount relevance for the design of aeronautical low-pressure turbines and compressors. It is concluded that the effect of the Mach number on the unsteady pressure phase can be neglected in first approximation and that the unsteadiness of the vibrating and adjacent airfoils is driven by vortex shedding mechanisms. Finally a simple model to estimate the work-per-cycle as a function of the reduced frequency and Mach Number is provided. The edge-wise and torsion modes are presented in less detail but it is shown that acoustic waves are essential to explain its behaviour. The non-dimensional work-per-cycle of the edge-wise mode shows a large dependence with the Mach number while in the torsion mode a large number of airfoils is needed to reconstruct the work-per-cycle departing from the influence coefficients.

Effect of the pulse trajectory on ultrasonic fluid velocity measurement

In a general situation a non-uniform velocity field gives rise to a shift of the otherwise straight acoustic pulse trajectory between the transmitter and receiver transducers of a sonic anemometer. The aim of this paper is to determine the effects of trajectory shifts on the velocity as measured by the sonic anemometer. This determination has been accomplished by developing a mathematical model of the measuring process carried out by sonic anemometers; a model which includes the non-straight trajectory effect. The problem is solved by small perturbation techniques, based on the relevant small parameter of the problem, the Mach number of the reference flow, M. As part of the solution, a general analytical expression for the deviations of the computed measured speed from the nominal speed has been obtained. The correction terms of both the transit time and of the measured speed are of M 2 order in rotational velocity field. The method has been applied to three simple, paradigmatic flows: one-directional horizontal and vertical shear flows, and mixed with a uniform horizontal flow.

Sonic anemometry of planetary atmospheres

The application of a recently developed model of sonic anemometers measuring process has revealed that these sensors cannot be considered as absolute ones when measuring spectral characteristics of turbulent wind speed since it is demonstrated that the ratios of measured to real spectral density functions depend on the composition and temperature of the considered planetary atmosphere. The new model of the measuring process of sonic anemometers is applied to describe the measuring characteristics of these sensors as fluid/flow dependent (against the traditional hypothesis of fluid/flow independence) and hence dependent on the considered planetary atmosphere. The influence of fluid and flow characteristics (quantified via the Mach number of the flow) and the influence of the design parameters of sonic anemometers (mainly represented by time delay between pulses shots and geometry) on turbulence measurement are quantified for the atmospheres of Mars, Jupiter, and Earth. Important differences between the behavior of these sensors for the same averaged wind speed in the three considered atmospheres are detected in terms of characteristics of turbulence measurement as well as in terms of optimum values of anemometer design parameters for application on the different considered planetary atmospheres. These differences cannot be detected by traditional models of sonic anemometer measuring process based on line averaging along the sonic acoustic paths.

On multiple-path sonic anemometer measurement theory

In this paper, a model of the measuring process of sonic anemometers with more than one measuring path is presented. The main hypothesis of the work is that the time variation of the turbulent speed field during the sequence of pulses that produces a measure of the wind speed vector affects the measurement. Therefore, the previously considered frozen flow, or instantaneous averaging, condition is relaxed. This time variation, quantified by the mean Mach number of the flow and the time delay between consecutive pulses firings, in combination with both the full geometry of sensors (acoustic path location and orientation) and the incidence angles of the mean with speed vector, give rise to significant errors in the measurement of turbulence which are not considered by models based on the hypothesis of instantaneous line averaging. The additional corrections (relative to the ones proposed by instantaneous line-averaging models) are strongly dependent on the wave number component parallel to the mean wind speed, the time delay between consecutive pulses, the Mach number of the flow, the geometry of the sensor and the incidence angles of mean wind speed vector. Kaimal´s limit k W1=1/l (where k W1 is the wave number component parallel to mean wind speed and l is the path length) for the maximum wave numbers from which the sonic process affects the measurement of turbulence is here generalized as k W1=C l /l, where C l is usually lesser than unity and depends on all the new parameters taken into account by the present model.

On sonic anemometer measurement theory

In this paper a model for the measuring process of sonic anemometers (ultrasound pulse based) is presented. The differential equations that describe the travel of ultrasound pulses are solved in the general case of non-steady, non-uniform atmospheric flow field. The concepts of instantaneous line-average and travelling pulse-referenced average are established and employed to explain and calculate the differences between the measured turbulent speed (travelling pulse-referenced average) and the line-averaged one. The limit k1l=1 established by Kaimal in 1968, as the maximum value which permits the neglect of the influence of the sonic measuring process on the measurement of turbulent components is reviewed here. Three particular measurement cases are analysed: A non-steady, uniform flow speed field, a steady, non-uniform flow speed field and finally an atmospheric flow speed field. In the first case, for a harmonic time-dependent flow field, Mach number, M (flow speed to sound speed ratio) and time delay between pulses have revealed themselves to be important parameters in the behaviour of sonic anemometers, within the range of operation. The second case demonstrates how the spatial non-uniformity of the flow speed field leads to an influence of the finite transit time of the pulses (M≠0) even in the absence of non-steady behaviour of the wind speed. In the last case, a model of the influence of the sonic anemometer processes on the measurement of wind speed spectral characteristics is presented. The new solution is compared to the line-averaging models existing in the literature. Mach number and time delay significantly distort the measurement in the normal operational range. Classical line averaging solutions are recovered when Mach number and time delay between pulses go to zero in the new proposed model. The results obtained from the mathematical model have been applied to the calculation of errors in different configurations of practical interest, such as an anemometer located on a meteorological mast and the transfer function of a sensor in an atmospheric wind. The expressions obtained can be also applied to determine the quality requirements of the flow in a wind tunnel used for ultrasonic anemometer calibrations.

Projectile Combustion Effects on Ram Accelerator Performance

Thesis (Master's)--University of Washington, 2016-06

Effects of hydrogen-air non-equilibrium chemistry on the performance of a model scramjet thrust nozzle

Two aspects of hydrogen-air non-equilibrium chemistry related to scramjets are nozzle freezing and a process called 'kinetic afterburning' which involves continuation of combustion after expansion in the nozzle. These effects were investigated numerically and experimentally with a model scramjet combustion chamber and thrust nozzle combination. The overall model length was 0.5m, while precombustion Mach numbers of 3.1 +/- 0.3 and precombustion temperatures ranging from 740K to 1,400K were involved. Nozzle freezing was investigated at precombustion pressures of 190kPa and higher, and it was found that the nozzle thrusts were within 6% of values obtained from finite rate numerical calculations, which were within 7% of equilibrium calculations. When precombustion pressures of 70kPa or less were used, kinetic afterburning was found to be partly responsible for thrust production, in both the numerical calculations and the experiments. Kinetic afterburning offers a means of extending the operating Mach number range of a fixed geometry scramjet.

Dynamic force balances for short-duration hypersonic testing facilities

Two force balance techniques for use in hypersonic impulse facilities are compared by measuring the drag force on a 30° semi-apex-angle blunt cone model in a hypersonic shock tunnel at a free stream Mach number of 5.75. An accelerometer-based balance and a stress-wave force balance were tested simultaneously on the same model to measure the drag force. It was found that drag force measurements could be made using both techniques in a flow with a 450-μ s test period. The measured drag forces compared well with the theoretical values estimated using Newtonian theory.

Performance of a scramjet combustor with combined normal and tangential fuel injection

The performance of a scramjet combustor with combined normal and tangential injection was experimentally investigated. Experiments were performed on a 500-mm cylindrical scramjet combustor at a freestream Mach number of 4.5, a nozzle supply pressure of 35.8 MPa, and a nozzle supply enthalpy of 5.8 MJ/kg. Hydrogen fuel was injected normally through portholes to promote combustion and tangentially through a slot to reduce viscous drag. A series of fuel injectors were used to vary the proportion of tangential to normal fuel between 45 and 100%. Reductions in the viscous drag of up to 25% were observed with the greatest reductions occurring at the lowest total equivalence ratio tested for each injector. However, the average pressure produced by combustion with combined normal and tangential injection was approximately 50% less than that produced by normal injection alone. An analysis of the change in specific impulse of the scramjet combustor indicated that the best overall performance was produced by 100% normal injection.