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

Similar documents
Steady waves in compressible flow

DRAG PREDICTION AND VALIDATION OF STANDARD M549, 155mm PROJECTILE

AERODYNAMIC CHARACTERIZATION OF A CANARD GUIDED ARTILLERY PROJECTILE

LEE-SIDE FLOW SIMULATIONS OF CRUCIFORM WING- BODY CONFIGURATIONS AT INCOMPRESSIBLE MACH NUMBERS

Lecture-4. Flow Past Immersed Bodies

Brenda M. Kulfan, John E. Bussoletti, and Craig L. Hilmes Boeing Commercial Airplane Group, Seattle, Washington, 98124

Department of Mechanical Engineering

Module3: Waves in Supersonic Flow Lecture14: Waves in Supersonic Flow (Contd.)

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

Wall treatments and wall functions

COMPUTATIONAL SIMULATION OF THE FLOW PAST AN AIRFOIL FOR AN UNMANNED AERIAL VEHICLE

Given a stream function for a cylinder in a uniform flow with circulation: a) Sketch the flow pattern in terms of streamlines.

WALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS

Numerical study of battle damaged two-dimensional wings

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

RECENT near-sonic and low-sonic boom transport aircraft

MDTS 5705 : Aerodynamics & Propulsion Lecture 2 : Missile lift and drag. G. Leng, MDTS, NUS

A numeric far field model for support interference studies in a slotted wall wind tunnel (European Transonic Windtunnel)

SPC Aerodynamics Course Assignment Due Date Monday 28 May 2018 at 11:30

Given the water behaves as shown above, which direction will the cylinder rotate?

Aerodynamic Investigation of a 2D Wing and Flows in Ground Effect

NUMERICAL OPTIMIZATION OF THE SHAPE OF A HOLLOW PROJECTILE

Analysis of Missile Bodies with Various Cross sections and Enhancement of Aerodynamic Performance

Studies on the Transition of the Flow Oscillations over an Axisymmetric Open Cavity Model

In which of the following scenarios is applying the following form of Bernoulli s equation: steady, inviscid, uniform stream of water. Ma = 0.

Validation 3. Laminar Flow Around a Circular Cylinder

Numerical Investigation of the Transonic Base Flow of A Generic Rocket Configuration

CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH

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

High Speed Aerodynamics. Copyright 2009 Narayanan Komerath

CFD COMPUTATION OF THE GROUND EFFECT ON AIRPLANE WITH HIGH ASPECT RATIO WING

Simulation of Aeroelastic System with Aerodynamic Nonlinearity

CALIFORNIA STATE UNIVERSITY, NORTHRIDGE TRANSONIC FLOW PROBLEMS AS CFD INSTRUCTIONAL TOOLS. Daniel John Reyna

Department of Energy Sciences, LTH

TRANSONIC AERODYNAMIC AND SCALING ISSUES FOR LATTICE FIN PROJECTILES TESTED IN A BALLISTICS RANGE

Flight Vehicle Terminology

Investigation potential flow about swept back wing using panel method

Coupled Fluid and Heat Flow Analysis Around NACA Aerofoil Profiles at Various Mach Numbers

Correction of the wall interference effects in wind tunnel experiments

Surrogate Representation for the Aerodynamic Coefficients of Supersonic Airfoils

AEROSPACE ENGINEERING DEPARTMENT. Second Year - Second Term ( ) Fluid Mechanics & Gas Dynamics

Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza

Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2

ACTIVE CONTROL OF BASE PRESSURE IN SUDDENLY EXPANDED FLOW FOR AREA RATIO 4.84

Numerical Validation of Flow Through an S-shaped Diffuser

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

TURBULENT FLOW ACROSS A ROTATING CYLINDER WITH SURFACE ROUGHNESS

Wings and Bodies in Compressible Flows

A Numerical Study of Circulation Control on a Flapless UAV

Active Control of Separated Cascade Flow

Experimental investigation of flow control devices for the reduction of transonic buffeting on rocket afterbodies

Theoretical Gas Flow through Gaps in Screw-type Machines

Introduction to Aerospace Engineering

DRAG REDUCTION FOR HYPERSONIC RE- ENTRY VEHICLES

ADVANCES in NATURAL and APPLIED SCIENCES

Fluctuating Pressure Inside/Outside the Flow Separation Region in High Speed Flowfield

CFD Study of Flow Over Parallel Ridges with Varying Height and Spacing

International Journal of Modern Trends in Engineering and Research e-issn No.: , Date: 2-4 July, 2015

Experimental Evaluation of Aerodynamics Characteristics of a Baseline Airfoil

The Simulation of Wraparound Fins Aerodynamic Characteristics

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

Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration

CFD DESIGN OF A GENERIC CONTROLLER FOR VORTEX-INDUCED RESONANCE

Fluid Mechanics Prof. T. I. Eldho Department of Civil Engineering Indian Institute of Technology, Bombay

FUNDAMENTALS OF AERODYNAMICS

Applied Fluid Mechanics

A Multi-Dimensional Limiter for Hybrid Grid

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

Paper: ASAT PP

Numerical Computations of Supersonic Flow Over a Square Cross-Section Missile

6.1 According to Handbook of Chemistry and Physics the composition of air is

An Experimental Validation of Numerical Post-Stall Aerodynamic Characteristics of a Wing

Design of a Droopnose Configuration for a Coanda Active Flap Application. Marco Burnazzi and Rolf Radespiel

EFFECT OF FLARE ANGLE AND NATURAL VENTILATION ON THE AERODYNAMIC CHARACTERISTICS OF A TYPICAL RE-ENTRY BODY AT SUBSONIC AND TRANSONIC MACH NUMBERS

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

AN ENGINEERING LEVEL PREDICTION METHOD FOR NORMAL-FORCE INCREASE DUE TO WEDGE SECTIONS

Separable warhead mathematical model of Supersonic & Hypersonic Re-entry Vehicles

Transonic Aerodynamics Wind Tunnel Testing Considerations. W.H. Mason Configuration Aerodynamics Class

CFD Time Evolution of Heat Transfer Around A Bundle of Tubes In Staggered Configuration. G.S.T.A. Bangga 1*, W.A. Widodo 2

LONGITUDINAL STABILITY AND TRIM OF AN ARIANE 5 FLY-BACK BOOSTER

Evaluating the Side Force in Oseen Flow Slender Body Theory and Comparison with CFD Simulation Results

MODELLING AND CFD ANALYSIS OVER A MISSILE BLUNT BODY WITH SPHERICAL AEROSPIKE

Fundamentals of Aerodynamits

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

OpenFOAM Simulations for MAV Applications

Wind Tunnel Measurements of Transonic Aerodynamic Loads on Mine Clearing Darts

Numerical Simulations of the Mars Science! Laboratory Supersonic Parachute!

SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE CAVITY FLOWS

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

Experimental and Numerical Investigation of Flow over a Cylinder at Reynolds Number 10 5

Development of Multi-Disciplinary Simulation Codes and their Application for the Study of Future Space Transport Systems

Aerodynamics of Wedge-Shaped Deflectors for Jet Noise Reduction

58:160 Intermediate Fluid Mechanics Bluff Body Professor Fred Stern Fall 2014

Wind Tunnel Study of a Large Aerostat, CFD Validation

ABSTRACT I. INTRODUCTION

THE EFFECT OF LEADING EDGE ROUNDING ON A CARET WAVERIDER CONFIGURATION AT MACH 8.2

Aerodynamics. Lecture 1: Introduction - Equations of Motion G. Dimitriadis

A comparison of turbulence models for an impinging jet in a crossflow

Elliptic Trailing Edge for a High Subsonic Turbine Cascade

Separation in three-dimensional steady flow. Part 3: TOPOLOGY OF SOME REMARKABLE THREE-DIMENSIONAL FLOWS

Transcription:

16 th International Conference on AEROSPACE SCIENCES & AVIATION TECHNOLOGY, ASAT - 16 May 26-28, 2015, E-Mail: asat@mtc.edu.eg Military Technical College, Kobry Elkobbah, Cairo, Egypt Tel : +(202) 24025292 24036138, Fax: +(202) 22621908 Paper: ASAT-16-172-AE Numerical Investigation of Wind Tunnel Wall Effects on a Supersonic Finned Missile Loai A. Elmahdy *, M.Y.M Ahmed *, O.K. Mahmoud * and O.E. Abdel Hamid * Abstract: One major aspect of wind tunnel measurements is the influence of test section walls on the measured data. In this paper, a numerical investigation of the supersonic flow around a conventional missile configuration is conducted. Wall effects are addressed by comparing two computational domains, unbounded and bounded. The calculated aerodynamic coefficients are compared with own corresponding measured values. The results show a good agreement between the measured parameters and those derived from the bounded computational domain. The impact of the walls on the flow field structure is investigated by comparing the bounded and unbounded numerical simulations. Keywords: missile aerodynamics, CFD, wall effect C A C df C m C N X cp axial force coefficient. fore body drag coefficient. pitching moment coefficient. normal force coefficient. normalized center of pressure position. Nomenclature: 1. Introduction Generally, there are three possible methods to estimate aerodynamic coefficients for a flying body namely, experimental, computational, and analytical. Experimental methods in which the aerodynamic coefficients are measured on an aerodynamically similar model in a wind tunnel, should be the closest to real flights. However, they are expensive, time consuming, and need measurement facilities that may not be available. Moreover, some real flight conditions may not be replicable in wind tunnels. In the Computational methods, based on computational fluid dynamics, CFD approach, the flow around the body is simulated by solving the governing equations of the flow numerically in the computational domain around the body. Once the simulation is properly set, the computation cost and time can be more economic compared with experimental methods. The processes for the preparation of CFD simulations can be complicated however, they are indeed more flexible compared with experimental measurements. The advantage of this method is the ability for visualizing flow features. Finally in the analytical method, simplified, closed form, and empirical equations are used. These equations are the simplest and least expensive and time consuming. However, they are basically suitable for simple flight problems in terms of body geometry and flight conditions. They yield poor results for high flight speeds or incidences. * Egyptian armed forces

Wind tunnel measurements should be as close as possible to real flight. However, it has been shown that these measurements are subject to various types of errors. They are listed briefly as[1]: (1) Solid Blockage: the free stream velocity around the model increases due to the confinement of the flow. (2) Streamline curvature: The normal curvature of the flow around a model is straightened by the tunnel walls. For subsonic airfoils in the tunnel, more lift and moment will be measured at a given angle of attack. (3) Wake blockage: The blockage due to the wake of the model. This may influence the upstream flow especially in subsonic conditions (4) Scaling error: the scaled model gives values that are different from those of the full scale. (5) Shock reflections: the source of error is witnessed in transonic and supersonic free stream conditions. Shock waves that are created around the model are reflected on the tunnel walls. The reflected waves may impinge the model causing changes in the measured parameters. (6) Horizontal buoyancy: the boundary layer on the test section walls has the effect of reducing the flow area of the free stream. All the above mentioned errors are commonly referred as wall interference errors. Estimating and correcting the wind tunnel wall effects has drawn the attention of researchers for decades. For example, In 1944 Allen and Vincbnti[2] derived theoretical tunnel-wall corrections for an airfoil of finite thickness and camber in a two dimensional flow wind tunnel. In their study, they considered the effect of the wake of the airfoil and of the compressibility of the fluid with the restriction that the chord of the wing is small compared with the tunnel height. They presented correction equations for the different parameters measured in the wind tunnel at low Mach numbers, however it has only limited applications at higher Mach numbers. Thom and Klanfer[3] studied the effect of the presence of walls on an airfoil in a subsonic flow. They found that: the walls of the wind tunnel increase the surface velocity on the airfoil. To get an accurate results and a closer simulation for the real flow they suggested using a thinner airfoil, which leads to another wind tunnel error namely, scaling error. Maskell[4] studied the effect of blockage on a wing and presented empirical correction equations depending on the dynamic pressure and the value of base pressure coefficient which showed good results with the experimental work. Awbi[5] experimentally investigated the effect of wall constraint on rectangular cross-section prism in a subsonic wind tunnel for different blockage and depth-to-width ratios of the test section. He showed that the effect of confinement of the flow on the mean surface pressure on the model was not uniform in the stream wise direction. West and Apelt[6] studied the blockage effects on the base pressure and on the drag coefficients of circular cylinders using different blockage and aspect ratios. They found that for blockage ratios above 6% the pressure distribution differed from that with smaller ratios while the change in drag coefficient was small compared with that with no blockage. They also concluded that reduction in aspect ratio affects drag coefficient and base pressure coefficient which were similar to those associated with the increase in blockage ratio. Duraisamy[7] studied the effect of wall interference on steady and oscillating airfoils in a subsonic wind tunnel using various theoretical methods, compressible, inviscid and viscous computations, and compared his work with experimental data. He presented a methodology

that could be used to predict wall interference effects and to extrapolate linear and nonlinear wind tunnel data to free-air conditions. The above-mentioned studies focused on bodies exposed to subsonic free stream conditions. In the open literature, the wind tunnel wall effects in the supersonic free-stream conditions were treated only recently. Ivanov and Streltsov[8] in 2014 introduced a new method to the boundary interference problem for subsonic and transonic wind tunnel which is based on the idea of using controlled boundary layer on the test section walls. Doig[9] studied the aerodynamics of a projectile in proximity to a solid wall flying at Mach 2.4 at zero incidence using the computational fluid dynamics approach. Considering the spin of the projectile, he studied the effect of changing the distance from the wall. It was concluded that the proximity increased the drag and generated a normal force on the projectile for wall clearance less than one caliber. For values of wall clearance between 0.4 and 1 caliber, the pitching moment produced was nose-down, for lower clearances, a more distinct nose-up trend was produced. Finally, Roychowdhury and Unnikrishnan[10] investigated the wall interference and solid blockage on a supersonic blunt model representing the payload of a launch vehicle. The simulation was axisymmetric (for zero angle of attack), at Mach numbers 1.5 and 3 for various blockage ratios. This approach implies cylindrical test section. It was found that when the blockage ratio is small, the detached shock extends to the tunnel wall reflecting back to the body various times which makes oscillations in the pressure distribution on the surface of the bod. The study presented recommendations for the blockage ratio and the model length to get smooth pressure distribution un-affected by the wall presence. However the study lacks for validation with experimental measurements. The present paper aims to shed more light on the wall effects in supersonic wind tunnel measurements. The supersonic flow around a model of a finned missile at incidence is simulated numerically. The free stream Mach number varies from 1.5 to 4 at 4⁰ incidence angle. To address the impact of test section walls, two 3-dimensional computational domains are investigated namely, open (unbounded) and closed (bounded with test section walls). The flow field features in both domains are explored and compared, and the overall aerodynamic properties are compared with own wind tunnel measurements. The remainder of the paper is organized as follows, details of the case study are presented next followed by the simulation setup details. The main results are then discussed and paper finalizes with the key findings. 2. Case study 2.1 Model The model investigated is a 1:16 scaled missile body configuration composed of conecylinder body, two outer conduits, and four fins. The model has a total length and caliber of 551.25 mm and 34 mm, respectively. The sketch of the model configuration is shown in Fig. 1

Fig. 1 Full model configuration 2.2 Wind tunnel measurements The model is tested in a tri-sonic wind tunnel which test section dimensions are 0.6 0.6m 2 and the length is 1.575m. The test conditions are variable from Mach 0.4 to 4.45, the corresponding Reynolds number of the flow at the test section inlet varies from 8.7e06 to 26.5e6. A set of conventional attack angle mechanism is installed in the test section of the wind tunnel, which can change the attack angle in the range of 15 ~38. The aerodynamic loads on the model are directly measured using a sting balance that is internally fit in the model base. Fig. 2 shows the model fitted inside the test section. Fig. 2 model installation in the wind tunnel test section

The accuracy of recorded data in measurements of the aerodynamic loads are listed in Table 1. Table 1 Recorded measurement errors for different parameters CN CAF CAB 0.0016 0.0042 0.0007 Where CN is the error in measured normal force coefficient. CAF is the error in measured axial force coefficient on the forebody. CAB is the error in measured axial force coefficient on the base. 2.3 Computational domain As stated earlier, CFD simulation is used to illustrate the effect of the presence of test section walls on the aerodynamic parameters of the model and the flow features around the model. To do so, bounded and unbounded flow simulations are conducted in two different domains each of which should have the proper definitions of boundaries. Both domains are threedimensional. Since only the incidence angle is considered, the flow around the model is pitchplane symmetric. Thus, only half 3D domains are constructed. The bounded domain is simply identical to the test section of the wind tunnel. The domain is defined by the dimensions of the test section, whereas the model is placed inclined in the test section orientation in the middle of the test section yielding the angle of attack of the model with respect to the flow through the test section. In Fig. 3 the bounded computational domain is illustrated. To replicate the wind tunnel experiments as closely as possible, the sting balance is added to the model. It is represented as cylindrical wall boundary that extends from the model base to the end of the domain. Outlet (pressure outlet) Body (solid wall) Inlet (pressure inlet) Sting (wall) Pitch plane (symmetry) Exit (pressure outlet) Test section walls (slip wall) Fig. 3 Bounded computational domain Test section walls (slip wall) As shown in Fig. 3, the pitch plane is defined as symmetric boundary whereas the model surface is defined as solid wall with no-slip condition specified. All side walls of the test sections are represented as slip wall boundaries on which no boundary layer is created numerically. The inlet and exit of the test section are represented as pressure inlet and pressure outlet boundaries, respectively.

In contrast, the unbounded domain is designed so that it captures all possible flow features that would happen around the model in the unbounded flight. The unbounded domain has the form of a spherical domain in the middle of which the model is positioned. Domain extents are 4 calibers, 6 calibers, 20 calibers in the upstream, downstream, and the lateral directions, respectively. Such dimensions are found to adequately accommodate all flow features of the supersonic flow conditions in concern, in the unbounded domain, the model is drawn aligned with the coordinate system while the angle of attack is set by defining an inclined incoming free stream flow. Definitions of model surface, pitch plane, sting, and exit are similar to those in the bounded domain. The outer boundary of the domain is defined as pressure far field. Fig. 4 shows the unbounded computational domain. Inlet (pressure far-field) Outlet (pressure outlet) Model (wall) Sting (wall) Pitch plane (symmetry) Fig. 4 Unbounded computational domain 2.4 Domain discretization The two domains are then discretized. It s intended to generate grids that yield equal solution in both domains despite the difference in their configuration. the scoping methods available in the CFD solver[11] are thus adopted. Initially, a coarse unstructured mesh is generated in the domain. Then, the grid is refined at user-specified zones. These zones are chosen so as to enable resolving the key flow features. Two scoping methods are successively implemented; the element sizing then the sphere of influence sizing. Body element yields a clustered fine grid around the body whereas the sphere of influence is drawn in the areas of the domain where the shock and expansion waves take place. Five spheres of influence are drawn (marked by the thick circles in Fig. 5). Eventually, a nearly feature aligned grid is achieved. Fig. 5 shows the grid on the symmetry plane for both bounded and unbounded domains.

(a) Unbounded domain (b) Bounded domain Fig. 5 Domain discretization using sphere of influence scoping method A grid sensitivity check is then conducted to address how refining the grid improves the solution. This is achieved by continuously refining the grid and comparing the successive simulation results until no significant improvement is achieved. The solution is then said to be grid-independent solution. Three flow simulations are run based on these grids. The features of these grids and the resulting fore drag coefficients are compared with the experimental data Table 2. The improvement between the second and the third grid is minor, so the coarser grid (the second grid) is chosen. Table 2 Grid sensitivity check results Grid No. of cells C df % improvement 1 602061 0.376 ----------- 2 1825562 0.404 7.4 3 2216425 0.4189 3.6 2.5 Turbulence model The criterion for choosing the appropriate turbulence model for the present case is the nature of the flow. Spalart Almaras is designed specifically for applications where wall bounded flows are involved. [11]. 3. Results and discussion 3.1 features of the unbounded and bounded flow fields Mach contours and the stream lines along the symmetry plane for both bounded and unbounded domains at Mach 1.5, 2.5, and 4 are shown in Fig. 6 Contours and stream lines in both domains at different Mach numbers a, b, and c, respectively. Bounded domain is shown on the lift whereas unbounded domain is shown on the right.

(a)m= 1.5 (b) M=2.5 (c) M=4 Fig. 6 Contours and stream lines in both domains at different Mach numbers As shown in Fig. 6a, at Mach=1.5 a shock wave is created ahead of the model nose. The windward side of this conical shock is stronger that it s leeward side due to the incidence of the model. At the nose shoulder, expansion fan is created. At the beginning of the upper and lower ducts, an oblique shock is created. Another shock wave is created at the wing leading edge. At the base another expansion occurs ahead of the wake behind the model. In addition, the boundary layer around the model is thinner in the windward side than in the leeward side of the model. The wake is not fully developed due to the existence of the sting at the model base.

As the Mach number increases to 2.5 in Fig. 6b, the shock wave ahead of the model nose begins to be more oblique and expansion at the nose shoulder is weaker than that of Mach 1.5. Finally at Mach 4 Fig. 6c, the shock ahead of the model nose becomes closer to the body wall and the oblique shocks at the ducts diminish especially at the upper duct. Moreover, expansion at the lower nose shoulder is clearer than that of the upper. Also, expansion at the lower base is not as significant as that at lower Mach. In all Mach values, it seems that the curvature of streamlines around the model is confined in a small area away from the model. This can be addressed in both bounded and bounded domain. Away from the model and in the vicinity of the bounded domain boundaries, the streamlines are straight and are, hence, not affected by the presence of the tunnel walls. This may be owed to the nature of supersonic flow as well as the slenderness of the model. At Mach 1.5 in the bounded domain, the shock generated ahead of the model nose impinges the upper and lower surfaces in the bounded domain. The reflected shock (which is stronger in the windward side) interferes with the model closer to the model fins and base area. As Mach increases to 2.5, reflection migrates downstream such that the reflected wave interferes with the model wave far downstream of the model base. Pressure contours on the model surface Fig. 7 pressure contours on the model surfaceon the windward and leeward sides of the model in Mach 2. Windward Leeward Fig. 7 pressure contours on the model surface As shown in the above figure, the pressure at the nose tip is the highest value as it is nearly the stagnation pressure. The pressure on the lower side of the nose is much higher than that on the upper side that is the result for the incidence, the same feature for that at the lower an upper sides of the wing is shown in the figure. At the base, the pressure is the lowest value where the wake zone is generated. 3.2 Variation of aerodynamic coefficients with Mach number a. Normal force coefficient The change in the normal force coefficient with the Mach number is shown in Fig. 8 Change in normal force coefficient with the Mach number (α=4⁰) as predicted by both bounded and unbounded solutions along with the experimental measurements.

CA CN Paper: ASAT-16-172-AE 0.65 0.6 0.55 0.5 0.45 0.4 1 1.5 2 2.5 3 3.5 4 4.5 MACH NUMBER CFD BOUNDED CFD OPEN EXPERIMENTAL Fig. 8 Change in normal force coefficient with the Mach number (α=4⁰) The measured normal force coefficient decreases as the Mach number increases. The simulation results based on both of the computational domains follows the same trend as the experimental ones, but with some difference of about 10% at Mach number 1.5, which is being currently investigated. At higher Mach numbers, the values of the normal force coefficients based on both domains are almost equal to each other and to the measured values. This may indicate that the wind tunnel walls have insignificant effect on the measured model aerodynamic parameters at high Mach numbers. b. Axial force coefficient The change in axial force coefficient with the Mach number is shown in Fig. 9 0.65 0.6 0.55 0.5 0.45 0.4 0.35 0.3 0.25 0.2 1 1.5 2 2.5 3 3.5 4 4.5 MACH NUMBER CFD BOUNDED CFD OPEN EXPERIMENTAL Fig. 9 Change in axial force coefficient with the Mach number

Xcp Paper: ASAT-16-172-AE The measured axial force coefficient decreases as the Mach number increases. Both of the two computational domains are following the same trend, which means that the axial force coefficient is not affected by the presence of the test section walls. c. Normalized center of pressure position The variation of the normalized center of pressure position measured from the tip of body nose with the Mach number is shown in Fig. 10. 0.7 0.65 0.6 0.55 0.5 0.45 0.4 1 1.5 2 2.5 3 3.5 4 4.5 MACH NUMBER CFD BOUNDED CFD OPEN EXPERIMENTAL Fig. 10 change in normalized center of pressure with Mach number (α=4⁰) The experimental data shows that the center of pressure shifts towards the body nose as the Mach number increases. The bounded computational domain shows the same behavior with an exception at M=4. Similarly the unbounded computational domain results behavior agrees with that of the experimental ones. The definition of the simulation results at Mach 4 requires a deeper investigation. d. Drag polar To construct a drag polar based on simulation results, three additional simulations namely, at 10⁰, 14⁰, and 18⁰ incidences, are conducted at every Mach number for the bounded domain. Fig. 11 shows the drag polar for Mach numbers 2 and 4. The drag polar points conducted from the computational domains, bounded or unbounded generally obey the trend of those conducted from experimental measurements

CL Paper: ASAT-16-172-AE 4.5 4 3.5 α=18 3 α=14 2.5 EXPERIMENTAL M=4 2 1.5 1 0.5 α=4 α=10 CFD BOUNDED M=4 CFD unbounded M=4 EXPERIMENTAL M=2 CFD BOUNDED M=2 CFD UNBOUNDED M=2 0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2-0.5 CD Fig. 11 Drag polar for Mach numbers 2 and 4 4. Conclusion Numerical simulations are conducted to investigate the impact of the wind tunnel walls on the measured data in supersonic conditions. A model of a supersonic finned missile is taken as the case study and the free-stream conditions are taken corresponding to Mach 1.5 to 4 at 4 o incidence. The simulations are carried out using two three-dimensional computational domains, bounded and unbounded ones. Results from the simulations are mutually compared versus experimental measurements. Moreover, the flow field associated with both domains are analyzed to show the considered effects. For the case considered, numerical simulations have shown a fairly accurate and reliable tool in predicting the aerodynamic coefficients, as compared with the experimental measurements. It was found that, the presence of wind tunnel walls has insignificant impacts on the measured data in the tunnel. However, at low supersonic Mach numbers, slight impact can be addressed possibly due to the shock reflections on the walls of the test section. This should be taken into consideration when correcting the tunnel data to resemble the real flight conditions. The present study considered a fixed incidence angle and specific test section blockage. The study is currently extended to investigate the impact of test section blockage and incidence angle on the quality of the tunnel measurements and the associated flow field features. 5. References [1] D. Althaus, "Tunnel-Wall Corrections at the Laminar Wind Tunnel," Institute report, IAG, http://www. iag. uni-stuttgart. de/laminarwindkanal, 2003. [2] H. J. Allen and W. G. Vincenti, "Wall interference in a two-dimensional-flow wind tunnel, with consideration of the effect of compressibility," 1944. [3] A. Thom and L. Klanfer, Tunnel-Wall Effect on an Aerofoil at Subsonic Speeds: HM Stationery Office, 1957. [4] E. Maskell, "A theory of the blockage effects on bluff bodies and stalled wings in a closed wind tunnel," DTIC Document1963.

[5] H. Awbi, "Wind-tunnel-wall constraint on two-dimensional rectangular-section prisms," Journal of Wind Engineering and Industrial Aerodynamics, vol. 3, pp. 285-306, 1978. [6] G. S. West and C. J. Apelt, "The effects of tunnel blockage and aspect ratio on the mean flow past a circular cylinder with Reynolds numbers between 104 and 105," Journal of Fluid Mechanics, vol. 114, pp. 361-377, 1982. [7] K. Duraisamy, W. J. McCroskey, and J. D. Baeder, "Analysis of wind tunnel wall interference effects on subsonic unsteady airfoil flows," Journal of aircraft, vol. 44, pp. 1683-1690, 2007. [8] A. Ivanov and E. Streltsov, "Controlled Boundary Layer on the Solid Walls of Wind Tunnels: New Approach to The Boundary Interference Problem." [9] G. Doig, T. J. Barber, E. Leonardi, A. J. Neely, and H. Kleine, "Aerodynamics of a supersonic projectile in proximity to a solid surface," AIAA journal, vol. 48, pp. 2916-2930, 2010. [10] A. Roychowdhury, C. Unnikrishnan, and R. Aerodynamics, "Numerical Simulation of Wind tunnel Blockage and Wall Interference Effects at Supersonic Mach Number Flows," International Journal of Innovative Research in Technology & Science(IJIRTS) vol. 2, pp. 63-72, 2013. [11] I. ANSYS/FLUENT, "ANSYS/FLUENT, Inc. "Fluent 14.5 user guide,"", ed, 2014. [12] DSSC Technicl Report TR70120, 2012.