A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows

Similar documents
Calculation of The Vehicle Drag and Heating Reduction at Hypervelocities with Laser- Induced Air Spike

NAPC Numerical investigation of axisymmetric underexpanded supersonic jets. Pratikkumar Raje. Bijaylakshmi Saikia. Krishnendu Sinha 1

Introduction to Aerodynamics. Dr. Guven Aerospace Engineer (P.hD)

V (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)

Numerical Analysis of Swirl Effects on Conical Diffuser Flows

Taylor Dispersion Created by Robert P. Hesketh, Chemical Engineering, Rowan University Fall 2005

WALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS

To study the motion of a perfect gas, the conservation equations of continuity

Simulation of unsteady muzzle flow of a small-caliber gun

Validation 3. Laminar Flow Around a Circular Cylinder

Pressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers

FUNDAMENTALS OF AERODYNAMICS

Numerical investigation on the effect of inlet conditions on the oblique shock system in a high-speed wind tunnel

Hypersonic flow and flight

BRAZILIAN 14-X HYPERSONIC WAVERIDER SCRAMJET AEROSPACE VEHICLE DIMENSIONAL DESIGN AT MACH NUMBER 10

AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS

Effect of Applied Magnetic Field on Shock Boundary Layer Interaction

MIRAMARE - TRIESTE June 2001

Experimental Investigation of Ionization Sensors under Shock Tube Stagnation Conditions

Carbon Science and Technology

Flow field in the compressor blade cascade NACA

Numerical Investigation of Geometrical Influence On Isolator Performance

Numerical investigation of swirl flow inside a supersonic nozzle

TAU Extensions for High Enthalpy Flows. Sebastian Karl AS-RF

The Simulation of Wraparound Fins Aerodynamic Characteristics

CFD Simulation of Internal Flowfield of Dual-mode Scramjet

Uncertainty in airflow field parameters in a study of shock waves on flat plate in transonic wind tunnel

TABLE OF CONTENTS CHAPTER TITLE PAGE

A Plasma Torch Model. 1. Introduction

Analysis of the Cooling Design in Electrical Transformer

Shock and Expansion Waves

NUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE

Follow this and additional works at:

MAE 598 Project #1 Jeremiah Dwight

Development of Flow over Blunt-Nosed Slender Bodies at Transonic Mach Numbers

RE-ENTRY TRAJECTORY SIMULATION OF A SMALL BALLISTIC RECOVERABLE SATELLITE

CFD Analysis of the Effect of Material Properties of Nose Cone on the Heat Flux and Thermal Field during Re-entry of Space Vehicle

AERODYNAMIC SHAPING OF PAYLOAD FAIRING FOR A LAUNCH VEHICLE Irish Angelin S* 1, Senthilkumar S 2

CALCULATION OF SHOCK STAND-OFF DISTANCE FOR A SPHERE

Transient growth of a Mach 5.92 flat-plate boundary layer

A finite-volume algorithm for all speed flows

DRAG REDUCTION FOR HYPERSONIC RE- ENTRY VEHICLES

Active Control of Separated Cascade Flow

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

ADVANCES in NATURAL and APPLIED SCIENCES

The Computations of Jet Interaction on a Generic Supersonic Missile

Analysis of a Dual-Mode Scramjet Engine Isolator Operating From Mach 3.5 to Mach 6

Numerical Investigation of Shock wave Turbulent Boundary Layer Interaction over a 2D Compression Ramp

Compressible Flow. Professor Ugur GUVEN Aerospace Engineer Spacecraft Propulsion Specialist

Steady waves in compressible flow

INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---

Introduction to Aerospace Engineering

Shape Optimisation of Axisymmetric Scramjets

Comparison of Turbulence Models in the Flow over a Backward-Facing Step Priscila Pires Araujo 1, André Luiz Tenório Rezende 2

2 Navier-Stokes Equations

The deposition efficiency and spatial thickness distribution of films created by Directed

CFD ANALYSIS OF AERODYNAMIC HEATING FOR HYFLEX HIGH ENTHALPY FLOW TESTS AND FLIGHT CONDITIONS

COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL

All-Particle Multiscale Computation of Hypersonic Rarefied Flow

Computational Investigations of High-Speed Dual-Stream Jets

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Performance. 5. More Aerodynamic Considerations

Numerical Study of Boundary-Layer Receptivity on Blunt Compression-Cones in Mach-6 Flow with Localized Freestream Hot-Spot Perturbations

Numerical Investigation of Wind Tunnel Wall Effects on a Supersonic Finned Missile

HIGH SPEED GAS DYNAMICS HINCHEY

CFD ANALYSIS OF HYPERSONIC NOZZLE THROAT ANALYSIS

A Computational Study on the Thrust Performance of a Supersonic Pintle Nozzle

Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow

Numerical Analysis of a Helical Coiled Heat Exchanger using CFD

Detached Eddy Simulation on Hypersonic Base Flow Structure of Reentry-F Vehicle

DIRECT SIMULATION CALCULATIONS OF THE RAREFIED HYPERSONIC FLOW PAST A BACKWARD-FACING STEP

2. NUMERICAL INVESTIGATION OF HYPERSONIC INVISCID AND VISCOUS FLOW DECELERATION BY MAGNETIC FIELD

Theoretical analysis on the effect of divergent section with laminar boundary layer of sonic nozzles

Application of a Modular Particle-Continuum Method to Partially Rarefied, Hypersonic Flows

REDUCTION OF AERODYNAMIC HEATING AND DRAG WITH OPPOSING JET THROUGH EXTENDED NOZZLE IN HIGH ENTHALPY FLOW

The effect of momentum flux ratio and turbulence model on the numerical prediction of atomization characteristics of air assisted liquid jets

NUMERICAL INVESTIGATION ON THE FLOW CHARACTERISTICS OF A SUPERSONIC JET IMPINGING ON AN AXI-SYMMETRIC DEFLECTOR

FEDSM COMPUTATIONAL AEROACOUSTIC ANALYSIS OF OVEREXPANDED SUPERSONIC JET IMPINGEMENT ON A FLAT PLATE WITH/WITHOUT HOLE

NUMERICAL SIMULATION OF FLUID FLOW BEHAVIOUR ON SCALE UP OF OSCILLATORY BAFFLED COLUMN

CHARACTERIZATION OF A RADIATION DETECTOR FOR AIRCRAFT MEASUREMENTS

Application of Computational Fluid Dynamics in Discontinuous Unsteady Flow with Large Amplitude Changes; the Shock Tube Problem

EXTERNAL-JET (FLUID) PROPULSION ANALOGY FOR PHOTONIC (LASER) PROPULSION By John R. Cipolla, Copyright February 21, 2017

Numerical study on performance of supersonic inlets with various three-dimensional bumps

Concept: AERODYNAMICS

Aerodynamics Simulation of Hypersonic Waverider Vehicle

High Pressure Zone Capture Wing Configuration for High Speed Air Vehicles

MULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544

HEAT TRANSFER COEFFICIENT CHARACTERIZATION AT THE SOLAR COLLECTOR WALL-FLUID INTERFACE

Lecture1: Characteristics of Hypersonic Atmosphere

EVALUATION OF EFFECT OF SHAPE AND LENGTH OF SPIKE ON AERODYNAMICS PERFORMANCE OF SUPERSONIC AXI-SYMMETRIC BODIES

IX. COMPRESSIBLE FLOW. ρ = P

Review of Fundamentals - Fluid Mechanics

A Non-Intrusive Polynomial Chaos Method For Uncertainty Propagation in CFD Simulations

CFD Analysis of Micro-Ramps for Hypersonic Flows Mogrekar Ashish 1, a, Sivakumar, R. 2, b

[Prasanna m a*et al., 5(6): July, 2016] ISSN: IC Value: 3.00 Impact Factor: 4.116

DESIGN & COMPUTATIONAL FLUID DYNAMICS ANALYSES OF AN AXISYMMETRIC NOZZLE AT TRANSONIC FREE STREAM CONDITIONS

NUMERICAL STUDY OF THE AERODYNAMIC PERFORMANCE OF NACA 0012 IN THE PRESENCE OF AN UNSTEADY HEAT SOURCE

Analysis of Shock Motion in STBLI Induced by a Compression Ramp Configuration Using DNS Data

3. FORMS OF GOVERNING EQUATIONS IN CFD

CFD ANALYSIS OF CD NOZZLE AND EFFECT OF NOZZLE PRESSURE RATIO ON PRESSURE AND VELOCITY FOR SUDDENLY EXPANDED FLOWS. Kuala Lumpur, Malaysia

Transcription:

doi: 152/jatmv5i227 A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows André Carlos Fraile Jr 1, Mauricio Antoniazzi Pinheiro Rosa 1 ABSTRACT: This work presents a numerical analysis of the gas dynamic problem of localized and steady energy addition to uniform high-speed airflows Firstly, the general effects caused by a localized energy deposition on the flow were investigated, then an extensive parametric analysis concerning the effects of energy deposition rate and dispersion for different free stream flow speeds was performed As a general result, localized and steady energy deposition generates compression waves and constant property flow stream tube downstream to the source The parametric analysis results have shown that either increasing the rate or decreasing dispersion in the orthogonal direction to the flow of the deposited energy has the effect of enhancing property flow changes, which is even more pronounced for lower flow speeds In contrast, energy dispersion in the flow direction has presented very little effect on the flow changes The results of this numerical analysis should be very helpful in studies of energy addition applications in hypersonic KEYWORDS: High speed flow, Energy addition, Flow control, Numerical analysis INTRODUCTION The need of aerospace vehicles for high speeds discloses a series of aerodynamics problems, such as large pressure drag forces and strong shock waves It is known that, for a blunt body at supersonic speeds, a physical spike placed at its nose modifies the structure of the strong shock wave, which can reduce wave drag in an effective way, although this structure also requires an undesirable additional cooling system (Riggins et al, 1999, Riggins and Nelson, 1999, Knight, 23) Another way to control high speed flows, which has also been considered for the reduction of aerodynamic drag of aerospace vehicles, corresponds to energy deposition in a small region of the airflow upstream to the vehicle This is similarly capable of modifying the properties and path of the fluid elements and consequently of reducing the strength of the vehicle shock wave by modifying its shape (Riggins et al, 1999, Riggins and Nelson, 1999, Knight, 23) Besides this application, energetic techniques for high-speed flow control have been recently widely studied with several other purposes, such as inlet mass capture increase by scramjet engines operating at off design speeds, attitude control of supersonic vehicles, scramjet isolator shortening, among others Deposition of energy to airflows has been experimentally accomplished by several different means, such as plasma arcs, laser pulse, microwave, electron beam, glow discharges, and so on (Oliveira, 2c) This work is part of some development activities that are being carried out at the Institute for Advanced Studies (IEAv) regarding high-speed airflow control The advance of 1Instituto de Estudos Avançados São José dos Campos/SP Brazil Author for correspondence: André Carlos Fraile Jr Instituto de Estudos Avançados Divisão de Aerotermodinâmica e Hipersônica Trevo Coronel Aviador José Alberto Albano do Amarante, 1 CEP 1222-1 São José dos Campos/SP Brazil E-mail: fraile@ieavctabr Received: 3/12/12 Accepted: 1//13 J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

17 Fraile Jr, AC and Rosa, MAP experimental studies previously performed at the IEAv, with respect to the addition of laser energy pulses to investigate vehicle drag and heat transfer (Minnuci et al, 25, Oliveira, 2a, 2b, 2c, Salvador et al, 25, 27, 2), has shown the need of numerical analyses to obtain a better understanding of the phenomena generally observed and also to support new experiments Therefore, in this work, which focused on understanding the properties of an energized flow without the presence of bodies, it is presented a numerical investigation of a localized steady energy addition to highspeed airflows, by varying the energy deposition rate and the spatial dispersion for different free stream airflow speeds (Fraile Jr, 211) While this subject has already been studied, this paper presents a series of results showing airflow variables such as velocity, density, pressure, and temperature when several conditions were modified, and these results can provide information to another ongoing numerical study at IEAv related to wave drag reduction of blunt bodies in supersonic flows using energy addition Herein, energy addition to high-speed flows is studied from thermodynamics and gas dynamics standpoints by considering that the energy addition effectively represents the direct energy transfer to the air translational (thermal) mode Therefore, no plasma effects are considered and localized energy addition to the flow is simply treated as a volumetric heat source into account a two-dimensional representation in space by considering axial symmetry The physical (computational) domain, also shown in Fig, 1, contains a cylindrical volumetric heat source representing a localized steady energy deposition in a supersonic uniform airflow This cylindrical heat source is aligned with the flow The Navier-Stokes equations (Eq 1), with a volumetric source term representing the energy addition, were used in this work (FLUENT INC, 26, Fraile Jr, 211): дρ + (ρv) = дt д (ρv) + (ρvv) = p + ( = τ ) дt д (ρe) + [v (ρe + p)] = [k T + ( = τ v)] + Q дt where ρ is the flow density, p is the pressure,v is the velocity, τ represents the influence of viscosity on momentum and energy equations, k is the thermal conductivity, E is the energy per mass unit, and Q is the energy source term (per time and volume unit) Also, a small region of the domain receives energy and the term Q is the mathematical way employed to deposit energy to the airflow The tensor τ can be written as in Eq 2 (FLUENT INC, 26, Fraile Jr, 211): (1) METHODOLOGY = 2 τ = µ v + v T 3 vi (2) The schematic shown in Fig 1 represents the numerical problem treated in this work The analysis performed takes y pressure far field where µ is the molecular viscosity and I is the unit tensor The gas is considered ideal, therefore it can be written as in Eq 3: p = ρ RT (3) pressure far field pressure outlet where R is the gas constant source axis x Figure 1 Physical (computational) domain employed for calculation of a steady heat source in a supersonic flow The energy per mass can be defined as Eq (FLUENT INC, 26; Fraile Jr, 211): E = h p ρ where dh = c p dt + v 2 () 2 J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows 171 The numerical solution to this problem is obtained by using the Fluent solver (FLUENT INC, 26) for the Navier-Stokes equations, with the conditions of laminar flow and ideal gas, in a square grid that is well suitable for the geometry shown in Fig 1 Since the numerical domain is axisymmetric, the Navier-Stokes equations are solved in cylindrical coordinates (FLUENT INC, 26, Tannehill et al, 1997) The equations are discretized in accordance with an implicit method and a first order upwind spatial scheme A pressurebased numerical method is used with an algorithm of coupling between the pressure and velocity named coupled, for which the equations of continuity and momentum are solved simultaneously, and then the other ones of the model are solved separately The boundary conditions employed to solve the problem are also presented in Fig 1 The addition of energy to the flow in a region of the computational domain is modeled through a user-defined function (UDF), which is an available tool of the solver The steady-state solution for the problem is sought since only steady energy depositions in the flow are considered (FLUENT INC, 26) The supersonic free stream flows from left to right with a static pressure of p = 1,x1 5 Pa and a temperature of T = 3K While it is intended to use the results of this work for better understanding the effects of energy addition on supersonic airflow, they have been also used to study the effects of energy addition on blunt body wave drag reduction This fact justifies some of the variables used herein, like the heat source dimensions, which have the same magnitude (millimeters) of those used in shock tunnels, and also of the heat power presented in the following sections, which have been used in recent works at the IEAv RESULTS AND DISCUSSION This section presents the results concerning the modification of a uniform flow due to the presence of an energy source confined to a small flow region in order to understand how the source affects the flow properties General characteristics of the flow with localized steady energy addition In order to identify the key changes of the properties of a uniform flow caused by the addition of energy to a small region in space, the physical domain shown in Fig 1 is used to perform numerical calculations It is considered that energy is steadily deposited uniformly in a small cylindrical region (length l = 1 mm and radius r = 5 mm), which is centered at the origin of the reference axes, with total power of 71 W, in a uniform Mach free stream flow The square grid used in all calculations is 2 cells/mm, which has yielded good accuracy results The data used to validate this grid will be presented in this section after some basic effects of energy addition to the airflow have been discussed The results corresponding to the scalar fields of the flow properties, such as pressure (p), temperature (T), density (ρ), as well as the flow movement quantities such as Mach number (M), axial (V x ) and radial velocities (V y ) are presented in Fig 2 Figure 3 presents the outcomes for the flow variables along the symmetry axis (lower boundary in Fig 1) and in a cross-sectional plane downstream to the heat source perpendicular to the free stream flow direction (right-hand boundary in Fig 1) As can be seen in the scalar fields of Fig 2 and in the plots of Fig 3, inside the heat source the flow is heated up with the effect of compressing the flow and diverting some part obliquely with respect to the symmetry axis The localized energy deposition in the flow generates compression waves that propagate radially with the local sound speed and axially with the flow one Inside the heat source, the flow density increases only slightly while the Mach number is reduced mainly because of the increase in temperature The results in Fig 3 show that there are two distinct flow regions around the symmetry axis downstream to the source: one just behind it, of about six source radius length, where all flow properties change considerably while the flow expands (pressure and density are reduced) and accelerates accompanied by a certain reduction in the temperature; and the other region following the first one, where flow properties are basically constants along the source axis, ie the pressure returns to the free stream value while the other flow variables reach values that are different from the free stream ones, even for distances far downstream from the source Actually, this latter region, as can be observed in Figs 2 and 3b, is characterized by a long constant property stream tube, aligned to the symmetry axis, of about the same source radius Also, the wave fronts propagate radially in both regions, but only in the second one the flow between the constant property stream tube and the wave fronts is very close to the free stream one As can be noticed, in the constant J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

Fraile Jr, AC and Rosa, MAP 172 2-2 - 12 16 x/r 2 2 2 - - - - - 2 2 2 (e) - 12 16 x/r 12 16 x/r 2 2 2 32 2 2 2 32 2 2 2 12 32 2 2 2 32 - Vy(m/s) 16 131 12 73 2 1-15 16 375 35 325 3 275 25 12 32 y/r 12 12 16 x/r M 16 2 Vx(m/s) 16 131 12 1373 133 131 125 16 y/r - 2 2 12 32 19 132 97 62 5 12 72 71 639 567 95 23 352 2 ρ (kg/m3) 16 y/r y/r 12 T (K) 16 y/r y/r 16 2 p (Pa) 25x1+5 22x1+5 19x1+5 17x1+5 1x1+5 x1+5 - - 12 16 x/r (f) 12 16 x/r Figure 2 Scalar fields for a Mach free stream and a 71 W heat source, of: pressure; temperature; density; Mach; (e) axial velocity; and (f) radial velocity property stream tube, the temperature has doubled whereas the Mach number and density have reduced about 3 and 5%, respectively, of their free stream values However, it should be noted that the laminar model used in this work does not take into account diffusion and convection effects in flows far from walls, therefore it behaves as an inviscid one that justifies the later results On the other hand, in real life flows are always at least somewhat turbulent, thus it is expected that the influence of the source on the flow in terms of the distance from the source that all flow properties return to their free stream values is smaller than for laminar ones The results shown in this work were obtained with a square grid of 2 cells/mm (each millimeter along the source axis or the direction perpendicular to it is divided into 2 parts), which can be validated through comparison with others of different resolutions Figure presents some flow variables (pressure, temperature, Mach number, and density) along the source symmetry axis when the resolution of the square grid is modified Even with the 1 cells/mm mesh, it is already possible to capture the main effects discussed herein, and it is also clear that the use of more refined meshes provides more accurate results Therefore, it was chosen the 2 cells/mm mesh, instead of the higher resolution one, because this yields enough accurate results, for the present purpose, without being so computationally costly J Aerosp Technol Manag, São José dos Campos, Vol5, No 2, pp169-1, Apr-Jun, 213

A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows 173,T/T, ()$ 3 &)* 25 T/T &)$ 2 %)* 15 %)$ 1 ρ/ρ /, M/,, V/V %)% 1 %)$ $), 9 $)* 5 $)( 3 -!" -!# $ # " 12 %& 16 %' &$ 2 &# 2 &" 2 (& 32 -!" -!# $ # " %& 12 %' 16 &$ 2 &# 2 &" 2 (& 32 $)" $)+ 7 $)' 6 $)* 5 $)# - ρ/ρ/ -M/M -V/V p/p ρ/ρ, T/T, T/T 22 "'" 2 "'! 1 &'% 16 &'$ 1 &'# ρ/ρ T/T 12 &'" 1 &'!!'%! 2" # 6$ % y/r 1 &! 12 &" 1 &# ρ/ρ /, M/,, V/V &'& 1 &'! 9!'*!'% 7!') 6!'$ 5!'( + ρ/ρ / + M/ + V/V!'#! " 2 # $ 6 % &! 1 &" 12 &# 1 y/r Figure 3 Normalized flow properties for a Mach free stream and a 71 W heat source, along the lower (symmetry axis) and right-hand boundaries of the domain in Fig 1 In order to show that the boundary condition of free stream near the energy source is not affecting the results, it is possible to analyze the airflow when the distance between the source center and the left pressure far field boundary shown in Fig 1 is increased Figure 5 presents the pressure distribution along the symmetry axis for a Mach free stream and a 71 W heat source when it is centered in the same position of all results of this paper, =, and when it is centered 12 units of r away from its original position When the pressure distribution for = 12 is graphically shifted 12 units to the left along the axis (Fig 5b), both distributions are approximately coincident, which indicates that the far field boundary condition is not affecting the solution The same analysis was performed for other flow properties (temperature, density, and Mach number) and they showed the same behavior presented in Fig 5 In the next sections, the effects of energy deposition rate and dispersion as well as the free stream velocity in supersonic flows are analyzed based on the fluid behavior nearby, inside the source and also in both regions downstream to the source mentioned before, which actually are of more interest for practical applications Effects of energy deposition rate (heat source power) When the effects of energy addition to high speed flows are analyzed, it is also of interest to study how the rate of energy deposition changes the flow properties Figure 6 presents results of the flow properties along the symmetry axis whereas Fig 7 shows results in a downstream perpendicular plane to the symmetry axis for three different values of energy deposition rate (heat source power) J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

17 Fraile Jr, AC and Rosa, MAP p/p M M 3 ()$ 25 &)* 2 &)$ 15 %)* 1 %)$ 5 ρ/ρ T/T T/T $)* $)" -!" -!# $ # " 12 %& 16 %' &$ 2 &# 2 &" (&!"!# $ # " %& %' &$ &# &" (& #)& 2 #)$ ()" 3 ()' 36 ()# 3 ()& 32 ()$ 3 &)" 2 &)' 26 &)# 2 &)& 22 2 2 &)" Mesh: 26 &)' 1 cells/mm &)# 2 cells/mm 2 cells/mm 22 &)& 2 &)$ 1 %)" 16 %)' 1 %)# Mesh: 12 %)& 1 cells/mm 2 cells/mm 1 %)$ cells/mm 2 32 - - 12 16 2 2 2 32 Mesh: 1 cells/mm 2 cells/mm cells/mm / &)$ $)( 3 -!" -!# $ # " 12 %& 16 %' &$ 2 &# 2 &" 2 (& 32 -!" -!# $ # " 12 %& 16 %' 2 &$ &# 2 &" 2 (& 32 %)% %)$ 1 $), 9 $)" $)+ 7 $)' 6 $)* 5 $)# Mesh: 1 cells/mm 2 cells/mm cells/mm Figure Changes in flow properties for different grids (Mach free stream and 71 W heat source): ; T/T ; M; ρ/ρ ()$ 3 &)* 25 &)$ 2 Source center position: = = 12 ()$ 3 &)* 25 &)$ 2 Source center position: = = 12 ( ) %)* 15 p/p ρ/ρ %)* 15 %)$ 1 %)$ 1 5 $)* $)* -!" -!# $ # " 12 %& 16 %' 2 &$ 2 &# 2 &" 32 (& -!" -!# $ # " 12 %& 16 %' 2 &$ 2 &# 2 &" 32 (& 5 Figure 5 Pressure along the symmetry axis (Mach free stream and 71 W) for different source center positions: original distributions; distribution when source center in =12 is shifted 12 units of to the left J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows 175 M/ 1 9 7 6 5 3 2 9 11 W 11 W 7 71 W 71 W 7 5 5 3 2 3 1 2 1 - # " $ 12 %& 16 %' 2 &# 2 &" 2 &$ 32 - # " $ 12 %& 16 %' 2 &# 2 &" &$ 2 32 (& 11 W 71 W ρ/ρ T/T - 12 16 2 2 2 32-12 16 2 2 2 32 1 12 1 6 2 11 W 71 W Figure 6 Changes in flow properties with the power ( =,) along the symmetry axis: ; T/T ; M/ ; ρ/ρ It is observed in Fig 6 that the higher the heat source power the more pronounced the effects on the flow properties, as previously discussed In other words, increasing the source power, the pressure, temperature, and density also increase inside the heat source, making the flow compression higher and consequently diverting more flow obliquely to the symmetry axis Furthermore, by increasing the source power, the expansion of the flow is more pronounced (higher variable changes) in the region immediately downstream to the source, although its length is not much affected by the source power for the same free stream velocity Considering the effects on the flow for distances beyond the downstream expanding region, that is, where the pressure along the symmetry axis has already basically returned to its free stream value, it can be inferred from Figs 6 and 7 that by increasing the power the axial length of the first region behind the source, where the flow properties are still changing, is not much affected (Fig 6) The changes in temperature, Mach number and density get more pronounced inside the flow stream tube (Figs 6 and 7), the radius of the stream tube gets only slightly bigger (Fig 7), and the wave fronts get stronger and slightly faster (Figs 7a and d) Effect of free stream flow speed Another aspect to be considered in this study is the effect of the free stream flow speed on the flow properties with localized steady energy addition to the flow Figures and 9 show their results along, respectively, the symmetry axis and the right-hand boundary in Fig 1, corresponding to three values for (for the same free stream sound speed) and 71 W power added to the flow J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

176 Fraile Jr, AC and Rosa, MAP M/ 13 125 7 11 W 71 W 12 5 15 1 95 9 5 3 2 1 2 6 1 12 1 y/r y/r 1 9 ρ/ρ T/T 12 1 6 11 W 71 W 2 6 1 12 1 7 6 11 W 71 W 2 11 W 71 W 5 2 6 1 12 1 2 6 1 12 1 y/r y/r Figure 7 Change of flow properties ( = ) along the cross section plane at the right-hand boundary as a function of power: ; T/T ; M/ ; ρ/ρ The results in Fig confirm the trend already observed that the most significant variations of the flow properties occur basically inside the heat source and in a small distance downstream from it Also, in Fig, it can be seen that the expansion region is elongated for higher Mach numbers, due to the increase in the flow speed, which corresponds to the transport of changes in flow property to greater distances from the source before reaching almost constant values When the free stream Mach number is higher, due to an increase in the free stream speed, the fluid residence time inside the heat source decreases and so does the rate of temperature change inside the source (lower peak) resulting in smaller alterations to the other flow properties not only inside the source but also downstream of it For instance (Fig ), the Mach number and density reductions in the flow stream tube for =7 are about 1% smaller than for = Therefore, the higher the free stream velocity the higher should be the rate of energy addition to the flow to keep its property percent changes in the stream tube Still, as observed in Fig 9, the stream tube radius is basically not affected by the free stream velocity However, as shown in Figs 9a and c, the wave fronts reach the right-hand boundary in Fig 1 at lower coordinates for higher free stream Mach numbers This happens mainly because of the higher free stream flow velocity, since the local sound speed is basically the same for all cases As a consequence, the diverted flow angle is smaller for higher free stream flow velocity It is also interesting to mention that the strength of the wave fronts is not very affected by the free stream velocity, as can be seen in Fig 9a J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows 177 M/ 3 275 = 25 25 225 = 7 2 2 175 15 15 = 125 1 1 = 7 5 75-12 16 2 2 2 32-12 16 2 2 2 32 1 9 7 6 5 ρ/ρ T/T = 1 = M 9 = 7 = 7 7 6 5 3-12 16 2 2 2 32-12 16 2 2 2 32 Figure Flow properties (71 W heat source) along the symmetry axis as a function of free stream flow velocity: ; T/T ; M/ ; ρ/ρ Effects of energy deposition dispersion (heat source size) This section presents the overall results of a study about the effects of the energy deposition dispersion on a supersonic flow This will be treated here by considering different heat source sizes, that is, by varying the radius (r) and length (l) of the cylindrical source, taking as reference the heat source used in the previous sections (l =1 mm and r =5 mm; 71 W total power) For the sake of clarity, the cylindrical source size (length and radius) will be written in the normalized form (l/l r/r ) Four different cases have been considered: the reference source (1 1); the double length (2 1); the double radius (1 2); and the double length and radius (2 2) In terms of the source volumetric heat rate, considering that the reference source is q, then the others are q /2, q / and q /, respectively Figures 1 and 11 show the flow property results along, respectively, the symmetry axis and the right-hand boundary in Fig 1, corresponding to the four cases mentioned for a Mach free stream and 71 W source total power As can be observed in these figures, the source length has very little effect on the main aspects of the flow behavior, although it seems that as the volumetric heat rate gets higher (reducing source radius), the source length begins to show some influence on the flow On the other hand, the source radius has major effects on the flow behavior, for example, decreasing the source radius: the compression inside the source increases; the expanding region behind the source gets shorter but still about six times the source radius; the differences between the constant flow properties inside the flow stream tube and their respective free stream ones get even bigger; the J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

17 Fraile Jr, AC and Rosa, MAP M/ 5 15 1 22 2 1 16 1 = 95 = 7 12 1 9 2 6 1 12 1 y/r y/r ρ/ρ T/T = = 7 2 6 1 12 1 15 12 1 95 9 1 9 5 7 = M 6 = M 75 = 7 5 = 7 7 2 6 1 12 1 2 6 y/r y/r 1 12 1 Figure 9 Flow properties (71 W heat source) along the cross-sectional plane at the right-hand boundary in Fig 1 as a function of free stream flow velocity: ; T/T ; M/ ; ρ/ρ flow stream tube radius is still about the source radius; and the wave fronts gain in strength but with unchangeable wave speed In conclusion, energy deposition dispersion in the flow direction (x-direction) has little effect on the flow, except for small radial dispersion whereas dispersion in the radial direction (y-direction) has major impact on the flow COMMENTS AND CONCLUSIONS This work has shown the effects on the flow structure of a localized and steady energy addition to uniform supersonic airflows, based only on thermodynamics and gas dynamics standpoints The knowledge obtained here about the highspeed flow changes due to localized energy addition certainly will be of great importance for the subsequent studies related to the applications of energetics in hypersonics at the IEAv The results presented here show that a localized energy source in a uniform high-speed airflow can modify its path lines at oblique angles, generating compression wave fronts whose strength increases with the rate of energy deposition and also with the reduction in radial dispersion Free stream velocity as well as axial dispersion has shown little effect on the wave front strength, although the oblique wave angle gets smaller for higher free stream speeds The results have shown that the flow downstream to the heat source may be separated into two distinct regions: one just behind it, which is characterized by considerable flow property variations; and the other, that follows the first one, where the flow properties are basically constants but different from those of the free stream one, with the exception of the pressure that returns to its original value J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

A Numerical Investigation of Localized and Steady Energy Addition to High Speed Airflows 179 M/ 3 25 2 15 1 ρ/ρ T/T 5-12 16 2 2 2 32-12 16 2 2 2 32 1 9 7 6 5 2 x 1 2 2 2 16 12 2 x 1 2 x 1 1 2 x 1 9 7 6 5 3-12 16 2 2 2 32-12 16 2 2 2 32 Figure 1 Flow properties ( =,) along the symmetry axis of the energy source as a function of source geometry: ; T/T ; M/ ; ρ/ρ Actually, this latter region is characterized by a long constant flow property stream tube aligned to the symmetry axis, of about the same source radius, with lower density and Mach number and higher temperature than their respective free stream values Between these two regions and the wave fronts, the flow remains very close to the free stream one The outcomes have also shown that either increasing energy deposition rate, or reducing free stream speed as well as decreasing radial dispersion, has the effect of increasing the changes in the flow properties inside the source and in both mentioned regions around the symmetry axis In contrast, axial dispersion has presented little effect on the flow, mainly when the radial dispersion is comparable or higher than the axial one Another important parameter is the first region length, which increases substantially by increasing the free stream speed or source radius Another interesting result is that for a fixed free stream speed, the length of this region is always given by a fixed number of source radius, for instance, for a Mach free stream, this region length is approximately six times the source radius Also, it is important to mention that the constant property stream tube radius has always approximately the source radius and is not much affected by source power, free stream speed, or axial dispersion It could be observed that the source radius (radial dispersion) has a major effect on the flow properties, whereas its length (axial dispersion) has a negligible one Therefore, highly focused energy deposition, mainly in the radial direction, will always need less energy deposition rate As the free stream speed increases, higher energy addition rate, or lower radial dispersion, is required to achieve about the same relative flow changes J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213

1 Fraile Jr, AC and Rosa, MAP M/ 5 2 22 2 x 1 2 15 1 95 9 2 6 1 12 1 1 16 1 12 1 y/r y/r ρ/ρ T/T 15 1 1 95 9 9 5 7 2 x 1 2 x 1 75 6 7 5 65 2 6 1 12 1 2 6 1 12 1 y/r y/r Figure 11 Flow properties ( =,) in a cross-sectional plane at the right-hand boundary in Fig 1 as a function of source geometry: ; T/T ; M/ ; ρ/ρ 2 x 1 2 6 1 12 1 REFERENCES FLUENT INC, 26, FLUENT 63 User s Guide, Retrieved in May 3, 213, from http://cdlab2fluidtuwienacat/lehre/turb/ FluentInc/fluent6326/help/html/ug/main_prehtm Fraile Jr, AC, 211, Um estudo numérico da redução de arrasto em corpos rombudos por adição de energia em escoamentos de altas velocidades (In Portuguese), Thesis, Instituto Tecnológico de Aeronáutica, 12 p Knight, D, 23, Survey of Aerodynamic Flow Control at High Speed by Energy Deposition, 1st AIAA Aerospace Sciences Meeting and Exhibit, 19 p Minucci, MAS et al, 25, Laser-supported directed-energy airspike in hypersonic flow, Journal of Spacecraft and Rockets, Vol 2, No 1, pp 51-57 Oliveira, AC et al, 2a, Bow shock wave mitigation by laserplasma energy addition in hypersonic flow, Journal of Spacecraft and Rockets, Vol 5, No 5, pp 921-927 Oliveira, AC et al, 2b, Drag Reduction by Laser-Plasma Energy Addition in Hypersonic Flow, Fifth International Symposium on Beamed Energy Propulsion, Vol 997, pp 379-39 Oliveira, AC, 2c, Investigação experimental da adição de energia por laser em escoamento hipersônico de baixa densidade (In Portuguese), Thesis, Instituto Nacional de Pesquisas Espaciais, 22 p Riggins, DW et al, 1999, Blunt-body wave drag reduction using focused energy deposition, AIAA Journal, Vol 37, No, pp 6-67 Riggins, DW, Nelson, HF, 1999, Hypersonic flow control using upstream focused energy deposition, Aerospace Sciences Meeting & Exhibit, 37, Reno Salvador, II et al, 25, Experimental Analysis of Heat Flux to a Blunt Body in Hypersonic Flow with Upstream Laser Energy Deposition Preliminary Results, th International Symposium on Beamed Energy Propulsion, Vol 3, pp 163-171 Salvador, II et al, 27, Surface Heat Flux and Pressure Distribution on a Hypersonic Blunt Body with DEAS, Fifth International Symposium on Beamed Energy Propulsion, Vol 997, pp 367-37 Salvador, II et al, 2, Experimental Analysis of Heat Flux to a Blunt Body in Hypersonic Flow with Upstream Laser Energy Deposition, Fifth International Symposium on Beamed Energy Propulsion, Vol 997, pp 379-39 Tannehill, JC et al, 1997, Computational fluid mechanics and heat transfer, Taylor and Francis, London, England, 792 p J Aerosp Technol Manag, São José dos Campos, Vol5, N o 2, pp169-1, Apr-Jun, 213