Position Error Studies for Air Data Boom of Flying Wing Configuration using RANS CFD Simulation

Similar documents
Estimation of Flow Field & Drag for Aerofoil Wing

Introduction to ANSYS CFX

Application of Wray-Agarwal Turbulence Model for Accurate Numerical Simulation of Flow Past a Three-Dimensional Wing-body

CFD Simulation of flow around angle of attack and sideslip angle vanes on a BAe Jetstream Part 1

Introduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.

Aerodynamic Analysis of Forward Swept Wing Using Prandtl-D Wing Concept

Debojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering

Modeling External Compressible Flow

Missile External Aerodynamics Using Star-CCM+ Star European Conference 03/22-23/2011

Research Article A Computational Investigation of Unsteady Aerodynamics of Insect-Inspired Fixed Wing Micro Aerial Vehicle s 2D Airfoil

Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil

A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS

NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING

ANSYS FLUENT. Airfoil Analysis and Tutorial

Keywords: CFD, aerofoil, URANS modeling, flapping, reciprocating movement

NUMERICAL INVESTIGATION OF THE FLOW BEHAVIOR INTO THE INLET GUIDE VANE SYSTEM (IGV)

AIR LOAD CALCULATION FOR ISTANBUL TECHNICAL UNIVERSITY (ITU), LIGHT COMMERCIAL HELICOPTER (LCH) DESIGN ABSTRACT

Estimating Vertical Drag on Helicopter Fuselage during Hovering

Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit

Ashwin Shridhar et al. Int. Journal of Engineering Research and Applications ISSN : , Vol. 5, Issue 6, ( Part - 5) June 2015, pp.

Research and Design working characteristics of orthogonal turbine Nguyen Quoc Tuan (1), Chu Dinh Do (2), Quach Thi Son (2)

Faculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia

Usage of CFX for Aeronautical Simulations

SHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY

Non-Newtonian Transitional Flow in an Eccentric Annulus

CFD ANALYSIS OF ENGINE AIR INTAKE FOR A PUSHER TYPE LIGHT TRANSPORT AIRCRAFT

Aerodynamic Analyses of Aircraft-Blended Winglet Performance

Investigation of a General Aviation Differential Pressure Angle of Attack Probe

FAR-Wake Workshop, Marseille, May 2008

Grid Dependence Study of Transonic/Supersonic Flow Past NACA Air-foil using CFD Hemanth Kotaru, B.Tech (Civil Engineering)

THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS

Analysis of an airfoil

Detached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads

Experimental study of UTM-LST generic half model transport aircraft

Computational Flow Analysis of Para-rec Bluff Body at Various Reynold s Number

Study of Swept Angle Effects on Grid Fins Aerodynamics Performance

Available online at ScienceDirect. Procedia Engineering 99 (2015 )

Effect of Position of Wall Mounted Surface Protrusion in Drag Characteristics At Low Reynolds Number

AERODYNAMIC DESIGN AND OPTIMIZATION TOOLS ACCELERATED BY PARAMETRIC GEOMETRY PREPROCESSING

Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit

CFD ANALYSIS OF AN RC AIRCRAFT WING

The Spalart Allmaras turbulence model

Team 194: Aerodynamic Study of Airflow around an Airfoil in the EGI Cloud

UNSTEADY RANS BASED IMPULSE RESPONSE STUDIES OF AGARD WING FOR AEROELASTIC AND FLUTTER ANALYSIS

Grid. Apr 09, 1998 FLUENT 5.0 (2d, segregated, lam) Grid. Jul 31, 1998 FLUENT 5.0 (2d, segregated, lam)

COMPUTATIONAL FLUID DYNAMICS ANALYSIS OF ORIFICE PLATE METERING SITUATIONS UNDER ABNORMAL CONFIGURATIONS

AERODYNAMIC DESIGN OF THE STRAKE FOR THE ROCKET PLANE IN TAILLESS CONFIGURATION.

AERODYNAMIC DESIGN OF FLYING WING WITH EMPHASIS ON HIGH WING LOADING

AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET

MASSACHUSETTS INSTITUTE OF TECHNOLOGY. Analyzing wind flow around the square plate using ADINA Project. Ankur Bajoria

High-Lift Aerodynamics: STAR-CCM+ Applied to AIAA HiLiftWS1 D. Snyder

An efficient method for predicting zero-lift or boundary-layer drag including aeroelastic effects for the design environment

Verification and Validation of Turbulent Flow around a Clark-Y Airfoil

Accepted Manuscript. Jetstream 31 national flying laboratory: Lift and drag measurement and modelling

DYNAMICS OF A VORTEX RING AROUND A MAIN ROTOR HELICOPTER

AERODYNAMIC SHAPES DESIGN ON THE BASE OF DIRECT NEWTON TYPE OPTIMIZATION METHOD

Simulation of Flow Development in a Pipe

CFD Analysis of conceptual Aircraft body

Design and Optimization of SUAV Empennage

COMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY. Abstract. I. Introduction

SONIC BOOM MINIMIZATION OF AIRFOILS THROUGH COMPUTATIONAL FLUID DYNAMICS AND COMPUTATIONAL ACOUSTICS

Numerical and theoretical analysis of shock waves interaction and reflection

Post Stall Behavior of a Lifting Line Algorithm

The viscous forces on the cylinder are proportional to the gradient of the velocity field at the

McNair Scholars Research Journal

Comparison of Results from Turbulence Models for the Nomad Flaperon-Configured Aerofoil

NUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF TEST MODELS AERODYNAMICS

September 17, 2003 Rev 2.0 (321-06) SUBJECT: Analysis of a Cylindrical pitot-static device for use in Air Flow Measurement

Summary and additional notes for Static ti pressure measurements:

Design and Development of Unmanned Tilt T-Tri Rotor Aerial Vehicle

The Numerical Simulation of Civil Transportation High-lift Configuration

Optimization of under-relaxation factors. and Courant numbers for the simulation of. sloshing in the oil pan of an automobile

Supersonic Flow Over a Wedge

MODELLING OF DYNAMIC STABILITY DERIVATIVES USING CFD

Compressible Flow in a Nozzle

Design Optimization of a Subsonic Diffuser. for a Supersonic Aircraft

INVESTIGATION OF HYDRAULIC PERFORMANCE OF A FLAP TYPE CHECK VALVE USING CFD AND EXPERIMENTAL TECHNIQUE

Recent developments for the multigrid scheme of the DLR TAU-Code

A B C D E. Settings Choose height, H, free stream velocity, U, and fluid (dynamic viscosity and density ) so that: Reynolds number

COMPUTATIONAL MODELING OF WIND ENERGY SYSTEMS

Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow

Viscous/Potential Flow Coupling Study for Podded Propulsors

Simulation of Turbulent Flow in an Asymmetric Diffuser

Numerical Simulation of Heat Transfer by Natural Convection in Horizontal Finned Channels

MSC Software Aeroelastic Tools. Mike Coleman and Fausto Gill di Vincenzo

ISSN(PRINT): ,(ONLINE): ,VOLUME-1,ISSUE-1,

AIAA A Computational Study of Subsonic Flows over A Medium Range Cargo Aircraft

Aerodynamic Design Optimization of UAV Rotor Blades using a Genetic Algorithm

Fluid-Structure Interaction Over an Aircraft Wing

Shock Wave Reduction via Wing-Strut Geometry Design

Driven Cavity Example

CFD Simulations of Flow over Airfoils:

Shape Optimization of Supersonic Bodies to Reduce Sonic Boom Signature

Numerical Simulations of Fluid-Structure Interaction Problems using MpCCI

APPLICATION OF A NAVIER-STOKES SOLVER TO THE ANALYSIS OF MULTIELEMENT AIRFOILS AND WINGS USING MULTIZONAL GRID TECHNIQUES

LES Applications in Aerodynamics

RV. Contextual Spelling 10. Grammar 24. Punctuation 2. Sentence Structure 1. Style 37. Vocabulary enhancement 21. of 100

Computational Fluid Dynamics (CFD) Simulation in Air Duct Channels Using STAR CCM+

CFD MODELING FOR PNEUMATIC CONVEYING

Algorithmic Developments in TAU

Transcription:

20 th Annual CFD Symposium August 9-10, 2018, NAL, Bangalore, India Position Error Studies for Air Data Boom of Flying Wing Configuration using RANS CFD Simulation G Prabhakaran*, Arvind Bobby Alphonse, B Saravanan, K Selvaraj ῼ Aerodynamics Division Aeronautical Development Establishment, Bangalore ABSTRACT: Accuracy in air speed measurements is critical for safe operation and optimal performance of the aircraft. If the measured air speed is less than the actual speed, aircraft may enter the stall regime under certain conditions. On the other hand if the measured air speed is more, then performance of the aircraft will be affected. Air speeds are measured by air data boom located ahead of the aircraft. The static pressure and total pressure measured by air data boom are translated into air speed and altitude. The static pressure data measured by air data boom are subject to instrument errors and position errors. Instruments error are due to errors in the standalone probe, plumbing errors and errors associated with calibration which are typically very small. The major portion of the error is the position error which is due to the location of static pressure port. This position error depends on the configuration of aircraft and the stand-off distance from nose of the aircraft. In this paper the effect of flying wing airframe on the position error is presented. A commercially available air data boom with pitot static probe, angle of attack vane and angle of side slip vane is opted for analysis. Studies are carried out for two different boom lengths and the effect of increase in standoff distance is brought out. This work also presents the computational fluid dynamics requirements and methodologies to estimate the position error for alpha vane, beta vane and static probe, which are located on the boom. The analysis is carried out at M=0.2 for the angle of attack range from -8 to 16.The measured velocity is less than the free stream velocity as expected in the lower range of angle of attack and it decreases with increasing angle of attack. At higher angle of attack, the probe measures a higher velocity as compared to the free stream due to local flow acceleration around the boom. Based on the computational results, the error in velocity and flow angles are quantified. Keywords: Alpha and Beta Vane, Air data Boom, Position Error, Flying Wing Nomenclature IRA = Indian Reference Atmosphere. AOA = Angle of Attack in degree. AOSS = Angle of Side Slip in degree. V = Free stream velocity in m/s. P o = Free stream total pressure in Pascal. P = Pressure Measured by static probe in Pascal. a = Speed of sound in m/s. V measured = Velocity measured by Vane in m/s; V measured = a ( 2 γ 1 ) γ 1 [(P o ) γ 1] P ρ = Free stream density in Kg/m 3. V error = Velocity error in Percentage; [(V measured -V ) X100 ] / (V ) P = Free Stream Static Pressure in Pascal. V x = Velocity in X-direction in m/s. V Y = Velocity in Y-direction in m/s. *Research Fellow Scientist Scientist ῼ Scientist

V z = Velocity in Z-direction in m/s. α vane = AOA indicated by AOA vane; tan -1 (V Z/V X). α actual = Free stream AOA. α error = AOA error; α vane - α actual vane = AOSS indicated by AOSS vane; tan -1 (V Y/V X) actual = Free Stream AOSS. error = AOSS error; vane - actual. C P = Pressure Coefficient; (P-P )/(0.5*ρ * V * V ) 1. Introduction: Pitot-static tube is an important air data sensor and is designed to measure the ambient static pressure and total pressure corresponding to the aircraft flight conditions. The pressure measured by tube is translated into aircraft speed, altitude, Mach number, vertical velocity information etc. that is necessary for various aircraft subsystems. The measurement of correct ambient pressure by the Pitot tube over the whole Mach number of aircraft is critical 1. Commonly, Pitot tube is placed ahead of the aircraft where the local static pressure in the subsonic flight condition is higher than the ambient static pressure, i.e. positive Cp value 1. The presence of the aircraft in the air stream causes error to the measuring instruments. The local flow angles also differ from the free-stream flow direction. Accurate measurements of flow angles, static pressure, and total pressure are necessary for safe flight operations. Several flight subsystems, such as auto flight controls, engine controls, and cockpit and cabin environmental control depend on accurate air data. References 5 and 6 specify the accuracy level required of pitot static tube by civil and military organizations. Flow angle vanes are small, mass balanced winglets designed to align with free stream of air during flight and it is the most common method for measuring the angle of attack and side slip in civil aviation. If the vane is fitted ahead of the airframe, it will allow the vanes to read the angles near to the free stream conditions. The vane is mounted behind the airframe and is prone to more error because of the influence of wing and propeller etc. Ideally, a system measuring flow angle such as vanes and differential pressure sensors or probes, should be fitted to a boom usually extended out from the aircraft i.e well away from influence of airframe on the flow field 2. All systems for measuring angle of attack and sideslip mentioned above are susceptible to error due to a number of factors 3. In this paper, performance characteristics of flow vanes and static probe are integrated on the air data booms and it is analyzed for various flight conditions. There are two types of boom configurations used for analysis, one is called as an original boom and another one is an extended boom (265mm longer than original boom). Computational Fluid dynamics (CFD) techniques have been used in the analysis to estimate the recommended extension length of the nose boom which is to be installed in the first flight of flying wing. 2. Computational Methodology: Commercial CFD software ICEM CFD, HiFUN and Tecplot are used as Pre-processor, Solver and Post processor respectively. The details of Pre-Processor, solver and post-processor are given below 2.1 Preprocessing The geometry is cleaned using ICEM CFD software. The air data boom is modeled without vanes as they will have negligible influence over the local flow field. Spherical computational domain is used for generating volume mesh around the model. The computational domain radius is taken as 10 times of the fuselage length. The mesh size of 17 million is used to capture all the geometry features and to keep the wall Y + value below 5. The typical grid has been shown in the Fig.1.1 and 1.2.

Fig.1.1 Flying Wing with Computational Domain Fig.1.2 Nose Boom 2.2 Solution The flow solver is a Reynolds Averaged NS (RANS) equation solver. Second order discretization scheme is used to solve for all flow variables. Inviscid and viscous flux calculation is carried out using Roe and Green Gauss schemes respectively. The Spalart Allmaras turbulence model is used for turbulent viscosity calculations. IRA sea level conditions are assumed for CFD analysis. Pressure far-field conditions have been used for the computational domain boundary. No-slip adiabatic wall condition is applied to all solid surfaces. The solution was deemed converged when the aerodynamic coefficients were changing less than 0.5% with iteration. The Free stream conditions and reference values used in the solver are given in Table 1. Table 1: Free stream Conditions and references Values Parameter Mach Number 0.2 Values Density(kg/m 3 ) 1.165 Pressure(N/m 2 ) 101325 Temperature(K) 303.15 Viscosity (Kg/m-s) 1.789 X 10-5 Ratio of Specific Heats 1.4 2.3 Post-Processing Tecplot is used for post processing the results of HiFUN solver. Data have been analyzed and presented in the subsequent section. The velocity components are measured in three directions using probe option and average pressure is calculated around static pressure probe. In Fig.2, the locations of alpha vane, beta vane, and static probe are shown.

Location of Alpha Vane Location of Static Pressure Probe Location of Beta Vane Fig.2 location of Vanes and Static pressure Probe 3. Boom Geometry Description The air data boom with pitot static probe is shown in Fig-3. Fig-3: Straight Nose Air Data Boom-Heated 4. Results and Discussion Because of the flow field created by the flying wing, the flow angle at any given location in the vicinity of the flying wing will generally differ from the true angle of attack of flying wing. At subsonic speeds the effects of the flow field extend in all directions. As the flow field around the each configuration is unique, in order to analyze the CFD results, first a vane center point was defined in Tecplot, then three velocity components are calculated at this point. The Computational simulation matrix is given in Table 3. Table 3: Summary of CFD Simulation matrix Parameter Values Mach Number 0.2 Angle of attack (deg) -8 to 16 in steps of 2 Beta(deg) 2 to 10 in steps of 2

The static pressure is measured in two ways during flight, with Pitot-static tube or static probe on the fuselage. The static pressure distribution on the air data boom is influenced by various factors like aircraft angle of attack, sideslip angle, speed etc. The surface distribution of the static pressure error (ΔP/q ) along the boom length is used to reveal the local variations due to major elements of aircraft, control surfaces, and antenna (see Fig.4).There is a gradual decrease in velocity error upto 10degree. Due to favorable pressure gradients around the air data boom there is an increase in error at higher angle of attack, which is shown in Fig-9. Favorable pressure gradients are observed along the boom length, at 12degree angle of attack due to local flow acceleration as shown in Fig-5 and 6. From Fig-10, it is observed that velocity error reduces with increasing sideslip. The variation of pressure coefficient across the boom at various locations in longitudinal direction for both booms are shown in Fig-7 and 8. It is observed that favorable pressure gradient increases along the boom length from the aircraft nose. There will be difference in local AOA measured by the vane and free stream AOA due to the location of the AOA vane. The local streamline direction is estimated by estimating the local velocity components. It follows from Fig-11 and Fig-12 that angle of attack error increases with angle of attack and angle of sideslip error increases with sideslip. In all cases, it is observed that extended boom is associated with lesser position error. 1.250 1.000 0.750 Aircraft Nose 0.500 ΔP/q 0.250 0.000-0.250-0.500-1.2-1 -0.8-0.6-0.4-0.2 0 X/D Fig-4: CFD simulated static pressure-error (ΔP/q ) distribution on the original boom plotted against boom centerline axis

Fig-5: Static Pressure Variation along the original boom (Side View) Fig-6: Static Pressure Variation along the extended boom (Side View) Fig-7: C P Variation across the original boom from aircraft nose for12deg angle of attack (Front View)

V error (%) Fig-8: C P Variation across the extended boom from aircraft nose for 12deg angle of attack (Front View) 3 2 1 0-1 -2-3 Extended Boom Original Boom -4-10 -5 0 5 10 15 20 AOA (deg) Fig-9: Variation of V error with AOA

AOA error (deg) V error (%) 0-0.5-1 -1.5-2 -2.5-3 Extended Boom Original Boom -3.5-4 0 2 4 6 8 10 12 AOSS (deg) Fig-10: Variation of V error with AOSS 8 6 4 2 0 Extended Boom Original Boom -10-5 0 5 10 15 20-2 -4 AOA (deg) Fig-11: Variation of AOA error with free stream AOA

AOSS error (deg) 3 2.5 2 1.5 1 0.5 Extended Boom Original Boom 0 0 2 4 6 8 10 12 AOSS (deg) Fig-12: Variation of AOSS error with free stream AOSS 5. Conclusion: There is a gradual decrease in velocity error upto 10deg for both the booms, the measured velocity is said to lesser than free stream velocity. An increase in velocity error is observed after 10 degree due to local acceleration around the boom, the measured velocity is said to higher than free stream velocity. Different length of air data boom are taken to estimate the error values such as velocity error, angle of attack and angle of sideslip error. Marginal reducation in error can be seen when the length of the boom is increased. The effect of flying wing airframe on air data boom is less significant when compared to other configurations.

6. References 1. Gracy,W., Measurement of Aircraft Speed and Altitude, NASA Reference Publication 1046,1980. 2. G.M Sakamoto. Aerodynamic characteristics of vane flow angularity sensor system capable of measuring flight path accelerations for Mach number range from 0.40 to 2.54. NASA technical Note,(TN D_8242),1976. 3. C.J.Bennett, N.J.Lawson, J.E. Gautery, A.Cooke. CFD simulation of flow around angle of attack and sideslip angle vanes on a bae jetstream 3102-part.Submitted to aerospace Science and technology,2016 4. Gracey, W., Measurement of Static Pressure on Aircraft, NACA Rept. 5. Wusest, W., Pressure and Flow Measurement, AGARD AG-160, Vo1.11, AGARD Flight Test Instrumentation Series, 1980. 6. Federal Aviation Administration, Federal Aviation Regulations, Volume 3,Transmittal m11,part25, Airworthiness standrads: Transport Category Airplanes U.S. Government Printing office,washington,dc,june.1974 7. U.S.Nay, Bureau of Aeronautics, Instrument Systems, Pitot Tube and Flush Static Port Operated, Installation of, MIL-I-6115A,Mar.1951. 8. Haering, E.A., Airdata Measurement and calibration, NASA TM-104316,1995 9. Installation and Inspection of Pitot-Static Pressure Systems U.S. Air Force MIL-P_26292 MIL-P-26292C, Wright-Patterson AFB, 1969. 10. J.A Lawford and K.R Nippress., Calibration of air-data systems and flow direction sensors. Advisory Group for Aerospace Research and Development flight test techniques series (AG-300), 1984. 11. J.D.Anderson. Computation Fluid Dynamics-First Edition McGraw-Hill Education, New York, USA, 1995.