Non-Equilibrium Reaction Rates in the Macroscopic Chemistry Method for DSMC Calculations

The Direct Simulation Monte Carlo (DSMC) method is used to simulate the flow of rarefied gases. In the Macroscopic Chemistry Method (MCM) for DSMC, chemical reaction rates calculated from local macroscopic flow properties are enforced in each cell. Unlike the standard total collision energy (TCE) chemistry model for DSMC, the new method is not restricted to an Arrhenius form of the reaction rate coefficient, nor is it restricted to a collision cross-section which yields a simple power-law viscosity. For reaction rates of interest in aerospace applications, chemically reacting collisions are generally infrequent events and, as such, local equilibrium conditions are established before a significant number of chemical reactions occur. Hence, the reaction rates which have been used in MCM have been calculated from the reaction rate data which are expected to be correct only for conditions of thermal equilibrium. Here we consider artificially high reaction rates so that the fraction of reacting collisions is not small and propose a simple method of estimating the rates of chemical reactions which can be used in the Macroscopic Chemistry Method in both equilibrium and non-equilibrium conditions. Two tests are presented: (1) The dissociation rates under conditions of thermal non-equilibrium are determined from a zero-dimensional Monte-Carlo sampling procedure which simulates ‘intra-modal’ non-equilibrium; that is, equilibrium distributions in each of the translational, rotational and vibrational modes but with different temperatures for each mode; (2) The 2-D hypersonic flow of molecular oxygen over a vertical plate at Mach 30 is calculated. In both cases the new method produces results in close agreement with those given by the standard TCE model in the same highly nonequilibrium conditions. We conclude that the general method of estimating the non-equilibrium reaction rate is a simple means by which information contained within non-equilibrium distribution functions predicted by the DSMC method can be included in the Macroscopic Chemistry Method.

Flush air data system calibration using numerical simulation

The use of computational fluid dynamics simulations for calibrating a flush air data system is described, In particular, the flush air data system of the HYFLEX hypersonic vehicle is used as a case study. The HYFLEX air data system consists of nine pressure ports located flush with the vehicle nose surface, connected to onboard pressure transducers, After appropriate processing, surface pressure measurements can he converted into useful air data parameters. The processing algorithm requires an accurate pressure model, which relates air data parameters to the measured pressures. In the past, such pressure models have been calibrated using combinations of flight data, ground-based experimental results, and numerical simulation. We perform a calibration of the HYFLEX flush air data system using computational fluid dynamics simulations exclusively, The simulations are used to build an empirical pressure model that accurately describes the HYFLEX nose pressure distribution ol cr a range of flight conditions. We believe that computational fluid dynamics provides a quick and inexpensive way to calibrate the air data system and is applicable to a broad range of flight conditions, When tested with HYFLEX flight data, the calibrated system is found to work well. It predicts vehicle angle of attack and angle of sideslip to accuracy levels that generally satisfy flight control requirements. Dynamic pressure is predicted to within the resolution of the onboard inertial measurement unit. We find that wind-tunnel experiments and flight data are not necessary to accurately calibrate the HYFLEX flush air data system for hypersonic flight.

Design, modelling and analysis of a six component force balance for hypervelocity wind tunnel testing

A combination of modelling and analysis techniques was used to design a six component force balance. The balance was designed specifically for the measurement of impulsive aerodynamic forces and moments characteristic of hypervelocity shock tunnel testing using the stress wave force measurement technique. Aerodynamic modelling was used to estimate the magnitude and distribution of forces and finite element modelling to determine the mechanical response of proposed balance designs. Simulation of balance performance was based on aerodynamic loads and mechanical responses using convolution techniques. Deconvolution was then used to assess balance performance and to guide further design modifications leading to the final balance design. (C) 2001 Elsevier Science Ltd. All rights reserved.

Heat transfer and flow behind a step in high enthalpy superorbital flow

Heat transfer levels have been investigated behind a rearward-facing step in a superorbital expansion tube. The heat transfer was measured along a flat plate and behind 2 and 3mm steps with the same length to step height ratio. Results were obtained with air as the test gas at speeds of 6.76kms(-1) and 9-60kms(-1) corresponding to stagnation enthalpies of 26MJ/kg and 48MJ/kg respectively. A laminar boundary layer was established on the flat plate and measured heat transfer levels were consistent with classical empirical correlations. In the case of flow behind a step, the measurements showed a gradual rise in heat transfer from the rear of the step to a plateau several step heights downstream for both flow conditions. Reattachment distance was estimated to be approximately 1.6 step heights downstream of the 2mm step at the low enthalpy condition through the use of flow visualisation.

Molecular Physics of Elementary Processes relevant to Hypersonics: atom-molecule, molecule-molecule and atom-surface processes.

In the present chapter some prototype gas and gas-surface processes occurring within the hypersonic flow layer surrounding spacecrafts at planetary entry are discussed. The discussion is based on microscopic dynamical calculations of the detailed cross sections and rate coefficients performed using classical mechanics treatments for atoms, molecules and surfaces. Such treatment allows the evaluation of the efficiency of thermal processes (both at equilibrium and nonequilibrium distributions) based on state-to-state and state specific calculations properly averaged over the population of the initial states. The dependence of the efficiency of the considered processes on the initial partitioning of energy among the various degrees of freedom is discussed.

On the development of an unstructured grid solver for inert and reactive high speed flow simulations

An unstructured grid Euler solver for reactive compressible flow applications is presented. The method is implemented in a cell centered, finite volume context for unstructured triangular grids. Three different schemes for spatial discretization are implemented and analyzed. Time march is implemented in a time-split fashion with independent integrators for the flow and chemistry equations. The capability implemented is tested for inert flows in a hypersonic inlet and for inert and reactive supersonic flows over a 2-D wedge. The results of the different schemes are compared with each other and with independent calculations using a structured grid code. The strengths and the possible weaknesses of the proposed methods are discussed.

Solution of aerospace problems using structured and unstructured strategies

Products developed at industries, institutes and research centers are expected to have high level of quality and performance, having a minimum waste, which require efficient and robust tools to numerically simulate stringent project conditions with great reliability. In this context, Computational Fluid Dynamics (CFD) plays an important role and the present work shows two numerical algorithms that are used in the CFD community to solve the Euler and Navier-Stokes equations applied to typical aerospace and aeronautical problems. Particularly, unstructured discretization of the spatial domain has gained special attention by the international community due to its ease in discretizing complex spatial domains. This work has the main objective of illustrating some advantages and disadvantages of numerical algorithms using structured and unstructured spatial discretization of the flow governing equations. Numerical methods include a finite volume formulation and the Euler and Navier-Stokes equations are applied to solve a transonic nozzle problem, a low supersonic airfoil problem and a hypersonic inlet problem. In a structured context, these problems are solved using MacCormack’s implicit algorithm with Steger and Warming’s flux vector splitting technique, while, in an unstructured context, Jameson and Mavriplis’ explicit algorithm is used. Convergence acceleration is obtained using a spatially variable time stepping procedure.

Brillouin spectroscopy : measurement of elastic and photoelastic constants of some alkali holide crystals

The assembly and testing of apparatus for the measurement of elastic and photoelastic constants by Brillouin scattering, using a Fabry-Perot interferometer and with argon ion laser excitation is described. Such measurements are performed on NaCI, KBr and LiF using the A = 488.0 nm laser line. The elastic constants obtained here are in very good agreement with the ultrasonic data for all three materials. The discrepancy between ultrasonic and hypersonic sound velocities which was reported by some authors for KBr and LiF is not confirmed, and the elastic constants obtained for LiF are the most accurate to date. Also, the present photoelastic constants are in good agreement with the data obtained by ultrasonic techniques for all three crystals. The results for the KBr and LiF crystals constitute the first set of photoelastic constants obtained for these materials by Brillouin spectroscopy. Our results for LiF are the best available to date.

Simulação numérica de escoamentos hipersônicos sobre corpos rombudos pelo método de elementos finitos

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Determinação de temperatura de chamas por espectroscopia de espécies quimioluminescentes

In a combustion process involving fossil fuels, there is the formation of species Chemiluminescent, especially CH*, C2* and OH*, whose spontaneous emission can be used as a diagnostic tool. In the present work, mapping and determination of the rotational temperature of the species CH* produced in flames on a burner fueled by Liquefied Petroleum Gas (LPG) was carried out. This study is part of a project involving the characterization of supersonic combustion in scramjets engines, whose study has been conducted in the hypersonic shock tunnel IEAv laboratories. The technique used was the natural emission spectroscopy, which has as main advantage of being non-intrusive. The rotational temperature determination was made using the Boltzmann method, whose principle is to relate the emission intensity of the species to the temperature by means of spectroscopic constants established.The temperature values were determined from the analysis of electronic bands AX and BX of the radical CH*. In order to confirm the results of flame temperatures obtained by the natural emission technique, was also used the technique of line reversal sodium. The results of both techniques showed that the temperature of the flames investigated is about 2500K a 2700K

Hirota method for oblique solitons in two-dimensional supersonic nonlinear Schrodinger flow

In a previous work El et al. (2006) [1] exact stable oblique soliton solutions were revealed in two-dimensional nonlinear Schrodinger flow. In this work we show that single soliton solution can be expressed within the Hirota bilinear formalism. An attempt to build two-soliton solutions shows that the system is "close" to integrability provided that the angle between the solitons is small and/or we are in the hypersonic limit. (C) 2012 Elsevier B.V. All rights reserved.

Diagnostic techniques for EHD and MHD interaction

The impact of plasma technologies is growing both in the academic and in the industrial fields. Nowadays, a great interest is focused in plasma applications in aeronautics and astronautics domains. Plasma actuators based on the Magneto-Hydro-Dynamic (MHD) and Electro- Hydro-Dynamic (EHD) interactions are potentially able to suitably modify the fluid-dynamics characteristics around a flying body without utilizing moving parts. This could lead to the control of an aircraft with negligible response time, more reliability and improvements of the performance. In order to study the aforementioned interactions, a series of experiments and a wide number of diagnostic techniques have been utilized. The EHD interaction, realized by means of a Dielectric Barrier Discharge (DBD) actuator, and its impact on the boundary layer have been evaluated by means of two different experiments. In the first one a three phase multi-electrode flat panel actuator is used. Different external flow velocities (from 1 to 20m/s) and different values of the supplied voltage and frequency have been considered. Moreover a change of the phase sequence has been done to verify the influence of the electric field existing between successive phases. Measurements of the induced speed had shown the effect of the supply voltage and the frequency, and the phase order in the momentum transfer phenomenon. Gains in velocity, inside the boundary layer, of about 5m/s have been obtained. Spectroscopic measurements allowed to determine the rotational and the vibrational temperature of the plasma which lie in the range of 320 ÷ 440°K and of 3000 ÷ 3900°K respectively. A deviation from thermodynamic equilibrium had been found. The second EHD experiment is realized on a single electrode pair DBD actuator driven by nano-pulses superimposed to a DC or an AC bias. This new supply system separates the plasma formation mechanism from the acceleration action on the fluid, leading to an higher degree of the control of the process. Both the voltage and the frequency of the nano-pulses and the amplitude and the waveform of the bias have been varied during the experiment. Plasma jets and vortex behavior had been observed by means of fast Schlieren imaging. This allowed a deeper understanding of the EHD interaction process. A velocity increase in the boundary layer of about 2m/s had been measured. Thrust measurements have been performed by means of a scales and compared with experimental data reported in the literature. For similar voltage amplitudes thrust larger than those of the literature, had been observed. Surface charge measurements led to realize a modified DBD actuator able to obtain similar performances when compared with that of other experiments. However in this case a DC bias replacing the AC bias had been used. MHD interaction experiments had been carried out in a hypersonic wind tunnel in argon with a flow of Mach 6. Before the MHD experiments a thermal, fluid-dynamic and plasma characterization of the hypersonic argon plasma flow have been done. The electron temperature and the electron number density had been determined by means of emission spectroscopy and microwave absorption measurements. A deviation from thermodynamic equilibrium had been observed. The electron number density showed to be frozen at the stagnation region condition in the expansion through the nozzle. MHD experiments have been performed using two axial symmetric test bodies. Similar magnetic configurations were used. Permanent magnets inserted into the test body allowed to generate inside the plasma azimuthal currents around the conical shape of the body. These Faraday currents are responsible of the MHD body force which acts against the flow. The MHD interaction process has been observed by means of fast imaging, pressure and electrical measurements. Images showed bright rings due to the Faraday currents heating and exciting the plasma particles. Pressure measurements showed increases of the pressure in the regions where the MHD interaction is large. The pressure is 10 to 15% larger than when the MHD interaction process is silent. Finally by means of electrostatic probes mounted flush on the test body lateral surface Hall fields of about 500V/m had been measured. These results have been used for the validation of a numerical MHD code.

Development and characterization of non-oxide ceramic composites for mechanical and tribological applications

The main reasons for the attention focused on ceramics as possible structural materials are their wear resistance and the ability to operate with limited oxidation and ablation at temperatures above 2000°C. Hence, this work is devoted to the study of two classes of materials which can satisfy these requirements: silicon carbide -based ceramics (SiC) for wear applications and borides and carbides of transition metals for ultra-high temperatures applications (UHTCs). SiC-based materials: Silicon carbide is a hard ceramic, which finds applications in many industrial sectors, from heat production, to automotive engineering and metals processing. In view of new fields of uses, SiC-based ceramics were produced with addition of 10-30 vol% of MoSi2, in order to obtain electro conductive ceramics. MoSi2, indeed, is an intermetallic compound which possesses high temperature oxidation resistance, high electrical conductivity (21·10-6 Ω·cm), relatively low density (6.31 g/cm3), high melting point (2030°C) and high stiffness (440 GPa). The SiC-based ceramics were hot pressed at 1900°C with addition of Al2O3-Y2O3 or Y2O3-AlN as sintering additives. The microstructure of the composites and of the reference materials, SiC and MoSi2, were studied by means of conventional analytical techniques, such as X-ray diffraction (XRD), scanning electron microscopy (SEM) and energy dispersive spectroscopy (SEM-EDS). The composites showed a homogeneous microstructure, with good dispersion of the secondary phases and low residual porosity. The following thermo-mechanical properties of the SiC-based materials were measured: Vickers hardness (HV), Young’s modulus (E), fracture toughness (KIc) and room to high temperature flexural strength (σ). The mechanical properties of the composites were compared to those of two monolithic SiC and MoSi2 materials and resulted in a higher stiffness, fracture toughness and slightly higher flexural resistance. Tribological tests were also performed in two configurations disco-on-pin and slideron cylinder, aiming at studying the wear behaviour of SiC-MoSi2 composites with Al2O3 as counterfacing materials. The tests pointed out that the addition of MoSi2 was detrimental owing to a lower hardness in comparison with the pure SiC matrix. On the contrary, electrical measurements revealed that the addition of 30 vol% of MoSi2, rendered the composite electroconductive, lowering the electrical resistance of three orders of magnitude. Ultra High Temperature Ceramics: Carbides, borides and nitrides of transition metals (Ti, Zr, Hf, Ta, Nb, Mo) possess very high melting points and interesting engineering properties, such as high hardness (20-25 GPa), high stiffness (400-500 GPa), flexural strengths which remain unaltered from room temperature to 1500°C and excellent corrosion resistance in aggressive environment. All these properties place the UHTCs as potential candidates for the development of manoeuvrable hypersonic flight vehicles with sharp leading edges. To this scope Zr- and Hf- carbide and boride materials were produced with addition of 5-20 vol% of MoSi2. This secondary phase enabled the achievement of full dense composites at temperature lower than 2000°C and without the application of pressure. Besides the conventional microstructure analyses XRD and SEM-EDS, transmission electron microscopy (TEM) was employed to explore the microstructure on a small length scale to disclose the effective densification mechanisms. A thorough literature analysis revealed that neither detailed TEM work nor reports on densification mechanisms are available for this class of materials, which however are essential to optimize the sintering aids utilized and the processing parameters applied. Microstructural analyses, along with thermodynamics and crystallographic considerations, led to disclose of the effective role of MoSi2 during sintering of Zrand Hf- carbides and borides. Among the investigated mechanical properties (HV, E, KIc, σ from room temperature to 1500°C), the high temperature flexural strength was improved due to the protective and sealing effect of a silica-based glassy phase, especially for the borides. Nanoindentation tests were also performed on HfC-MoSi2 composites in order to extract hardness and elastic modulus of the single phases. Finally, arc jet tests on HfC- and HfB2-based composites confirmed the excellent oxidation behaviour of these materials under temperature exceeding 2000°C; no cracking or spallation occurred and the modified layer was only 80-90 μm thick.

High frequency acoustics in colloid-based meso- and nanostructures by spontaneous Brillouin light scattering

Materials that can mold the flow of elastic waves of certain energy in certain directions are called phononic materials. The present thesis deals essentially with such phononic systems, which are structured in the mesoscale (<1 µm), and with their individual components. Such systems show interesting phononic properties in the hypersonic region, i.e., at frequencies in the GHz range. It is shown that colloidal systems are excellent model systems for the realization of such phononic materials. Therefore, different structures and particle architectures are investigated by Brillouin light scattering, the inelastic scattering of light by phonons.rnThe experimental part of this work is divided into three chapters: Chapter 4 is concerned with the localized mechanical waves in the individual spherical colloidal particles, i.e., with their resonance- or eigenvibrations. The investigation of these vibrations with regard to the environment of the particles, their chemical composition, and the influence of temperature on nanoscopically structured colloids allows novel insights into the physical properties of colloids at small length scales. Furthermore, some general questions concerning light scattering on such systems, in dispute so far, are convincingly addressed.rnChapter 5 is a study of the traveling of mechanical waves in colloidal systems, consisting of ordered and disordered colloids in liquid or elastic matrix. Such systems show acoustic band gaps, which can be explained geometrically (Bragg gap) or by the interaction of the acoustic band with the eigenvibrations of the individual spheres (hybridization gap).rnWhile the latter has no analogue in photonics, the presence of strong phonon scatterers, when a large elastic mismatch between the composite components exists, can largely impact phonon propagation in analogy to strong multiple light scattering systems. The former is exemplified in silica based phononic structures that opens the door to new ways of sound propagation manipulation.rnChapter 6 describes the first measurement of the elastic moduli in newly fabricated by physical vapor deposition so-called ‘stable organic glasses’. rnIn brief, this thesis explores novel phenomena in colloid-based hypersonic phononic structures, utilizing a versatile microfabrication technique along with different colloid architectures provided by material science, and applying a non-destructive optical experimental tool to record dispersion diagrams.rn

[The role of thoracoscopy in paediatric spontaneous pneumothorax]

Spontaneous pneumothorax in children is an emergency. A rapid and structured treatment is necessary. Indicating clinical signs are dyspnoea with unilateral reduced breath sounds and hypersonic percussion note. Chest x-ray confirms diagnosis. CT scan shows detailed information's of pathological changes in lung tissue and offers important considerations for the therapy. Today video assisted thoracoscopic surgery (VATS) takes a central role in the treatment of spontaneous pneumothorax. Thoracoscopy completes diagnosis detecting the pathologies, which were missed in imaging procedures. Definitive surgical treatment is performed simultaneously in the same anaesthesia.