AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET

Similar documents
AERODYNAMIC DESIGN OF FLYING WING WITH EMPHASIS ON HIGH WING LOADING

(c)2002 American Institute of Aeronautics & Astronautics or Published with Permission of Author(s) and/or Author(s)' Sponsoring Organization.

NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING

Introduction to ANSYS CFX

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

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

GRID PATTERN EFFECTS ON AERODYNAMIC CHARACTERISTICS OF GRID FINS

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

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

Estimation of Flow Field & Drag for Aerofoil Wing

Shock Wave Reduction via Wing-Strut Geometry Design

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

AN INVERSE DESIGN METHOD FOR ENGINE NACELLES AND WINGS

AERODYNAMIC DESIGN AND OPTIMIZATION TOOLS ACCELERATED BY PARAMETRIC GEOMETRY PREPROCESSING

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

INVERSE METHODS FOR AERODYNAMIC DESIGN USING THE NAVIER-STOKES EQUATIONS

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

Usage of CFX for Aeronautical Simulations

CFD Analysis of conceptual Aircraft body

State of the art at DLR in solving aerodynamic shape optimization problems using the discrete viscous adjoint method

Shape optimisation using breakthrough technologies

The Spalart Allmaras turbulence model

AIRFOIL SHAPE OPTIMIZATION USING EVOLUTIONARY ALGORITHMS

How to Enter and Analyze a Wing

Subsonic Airfoils. W.H. Mason Configuration Aerodynamics Class

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

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

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

Design Optimization of a Subsonic Diffuser. for a Supersonic Aircraft

Impact of Computational Aerodynamics on Aircraft Design

Second Symposium on Hybrid RANS-LES Methods, 17/18 June 2007

The Numerical Simulation of Civil Transportation High-lift Configuration

CFD ANALYSIS OF AN RC AIRCRAFT WING

LESSONS FROM WIND TUNNEL MODELS MADE BY RAPID PROTOTYPING

Performance improvement of a wind turbine blade using a developed inverse design method

Subsonic Airfoils. W.H. Mason Configuration Aerodynamics Class

A DRAG PREDICTION VALIDATION STUDY FOR AIRCRAFT AERODYNAMIC ANALYSIS

Study of Swept Angle Effects on Grid Fins Aerodynamics Performance

WAVERIDER DESIGH WITH PARAMETRIC FLOW QUALITY CONTROL BY INVERSE METHOD OF CHARACTERISTICS

4. RHEOELECTRIC ANALOGY

39th AIAA Aerospace Sciences Meeting and Exhibit January 8 11, 2001/Reno, NV

ACTIVE SEPARATION CONTROL WITH LONGITUDINAL VORTICES GENERATED BY THREE TYPES OF JET ORIFICE SHAPE

German Aerospace Center, Institute of Aerodynamics and Flow Technology, Numerical Methods

Progress and Future Prospect of CFD in Aerospace

OPTIMIZATIONS OF AIRFOIL AND WING USING GENETIC ALGORITHM

CONFIGURATION TEST CASES FOR AIRCRAFT WING ROOT DESIGN AND OPTIMIZATION

Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization

NUMERICAL AND EXPERIMENTAL INVESTIGATIONS OF TEST MODELS AERODYNAMICS

Abstract. Introduction

Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit

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

Estimating Vertical Drag on Helicopter Fuselage during Hovering

LOW-SPEED HIGH-LIFT PERFORMANCE IMPROVEMENTS OBTAINED AND VALIDATED BY THE EC-PROJECT EPISTLE

AERODYNAMIC OPTIMIZATION OF NEAR-SONIC PLANE BASED ON NEXST-1 SST MODEL

CAD-BASED WORKFLOWS. VSP Workshop 2017

Modeling External Compressible Flow

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

LAMDES User s Manual VLMpc

Supersonic Wing Design Method Using an Inverse Problem for Practical Application

Experimental study of UTM-LST generic half model transport aircraft

NUMERICAL SIMULATION OF 3D FLAPPING WING BASED ON CHIMERA METHOD

An advanced RBF Morph application: coupled CFD-CSM Aeroelastic Analysis of a Full Aircraft Model and Comparison to Experimental Data

Aerodynamic Design of a Tailless Aeroplan J. Friedl

Computational shock and Mach waves visualization aiding the development of aerodynamic design techniques

Introduction to Aerodynamic Shape Optimization

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

Optimization of Laminar Wings for Pro-Green Aircrafts

INTERACTIVE AERODYNAMICS ANALYSIS AND DESIGN PROGRAMS FOR USE IN THE UNDERGRADUATE ENGINEERING CURRICULUM

SONIC-BOOM PREDICTION USING EULER CFD CODES WITH STRUCTURED/UNSTRUCTURED OVERSET METHOD

SMALL HEIGHT DUCT DESIGN FOR MULTICOPTER FAN

Small Height Duct Design for 17 Multicopter Fan Considering Its Interference on Quad-copter

OPTIMISATION OF THE HELICOPTER FUSELAGE WITH SIMULATION OF MAIN AND TAIL ROTOR INFLUENCE

Computational Study of Unsteady Flows around Dragonfly and Smooth Airfoils at Low Reynolds Numbers

Available online at ScienceDirect. Procedia Engineering 99 (2015 )

Numerical Investigation of Transonic Shock Oscillations on Stationary Aerofoils

Simulation of Turbulent Flow around an Airfoil

Analysis of an airfoil

LES Analysis on Shock-Vortex Ring Interaction

What s New in AAA? Design Analysis Research. Version 3.3. February 2011

GAs for aerodynamic shape design II: multiobjective optimization and multi-criteria design

CFD++ APPLICATION ON WIND TUNNEL DATA ANALYSIS

ANSYS FLUENT. Airfoil Analysis and Tutorial

AIAA Using Flow Control Devices in Small Wind Tunnels

OVER-THE-WING-NACELLE-MOUNT CONFIGURATION FOR NOISE REDUCTION

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

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

FAR-Wake Workshop, Marseille, May 2008

Aerodynamic Inverse Design Framework using Discrete Adjoint Method

CFD Analysis of 2-D Unsteady Flow Past a Square Cylinder at an Angle of Incidence

Single and multi-point aerodynamic optimizations of a supersonic transport aircraft using strategies involving adjoint equations and genetic algorithm

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

NUMERICAL RESEARCH OF THE NACELLE STRAKE ON A CIVIL JET

Keisuke Sawada. Department of Aerospace Engineering Tohoku University

Design and Computational Fluid Dynamics Analysis of an Idealized Modern Wingsuit

DYNAMICS OF A VORTEX RING AROUND A MAIN ROTOR HELICOPTER

Comparison of CFD Simulation of a Hyundai I20 Model with Four Different Turbulence Models

Transition Flow and Aeroacoustic Analysis of NACA0018 Satish Kumar B, Fred Mendonç a, Ghuiyeon Kim, Hogeon Kim

A Surface Parameterization Method for Airfoil Optimization and High Lift 2D Geometries Utilizing the CST Methodology

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

COMPARISON OF SHOCK WAVE INTERACTION FOR THE THREE-DIMENSIONAL SUPERSONIC BIPLANE WITH DIFFERENT PLANAR SHAPES

Transcription:

25 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET Qingzhen YANG*,Yong ZHENG* & Thomas Streit** *Northwestern Polytechincal University, 772,Xi an, China, ** Institute of Aerodynamic and Flow Technology, DLR, 388 Braunschweig,Germany Keywords: Aerodynamic, Design, Blended Wing-Body, Inlet Abstract A wing-body blended configuration for transport and three wing-body blended flying wings for unmanned air vehicle are designed using box complex search approach optimization method and inverse design method. The twists of wing sections are optimized using composite method, and the airfoils on wing sections are designed using iterative residual correction method. The designed blended wing-body is analyzed using CFD method and its aerodynamic characteristic is obtained in detail. The embedding inlet of engine is designed and integrated into the wing-body blended configuration to investigate the influence of inlet on the aerodynamic performance. Numerical simulation shows that the inlet results in a favorable effect on the lift and pitching moment but unfavorable effect on the total drag. The stall characteristic becomes better due to the embedding inlet. Introduction The wing-body blended configuration is widely investigated in recent years due to its advantage of higher cruise efficiency [][2][3], larger loading space for civil aircraft, and also better electromagnet-ism stealth performance for military aircraft. The European transport solution VELA is a good example for civil usage. B2 is a well-known success in military application. Nevertheless, the stability and maneuverability difficulty is the disadvantage of thus configuration. In this paper twists of wing sections are optimized using box complex search approach to get maximum ratio of lift to drag and desired pitching moment characteristic. The iterative residual correction concept is used as inverse design method that modifies the profile of wing section according to the target pressure desired for special purpose. The natural laminar flow on the wing surface is desired. The embedding inlet and power system of blended wing-body () is designed to enhance drag performance and stealth characteristics. 2 Methods 2. Optimization Method Box complex search (M.J. Box, 965) approach is a direct search procedure to solve constrained nonlinear problems. The procedure is initiated by randomly placing a set of N+2 <=K <= 2N search points in the feasible region, where N is the number of search variables. A worst point is identified and replaced by a reflection point according to the given relation. This procedure is iterated, in each step discarding the least desirable point and replacing it with a new, and hopefully superior, search point, within the constrained bounds. The search is terminated and the best solution is taken as the optimal solution. The target function and constraints are aerodynamic characteristics computed using a CFD solver. Due to the large number of search variables, a fast working CFD analysis code is required, thus the full potential solver with three-dimensional boundary layer correction is chosen as CFD tool. The maximum ratio of lift to drag and minimum pitching moment at cruise situation is desired target. 2.2 Inverse Design Method The Takanashi [4] iterative residual correction concept is used for geometry modification according to the difference between initial

QINGZHEN YANG,YONG ZHENG & THOMAS STREIT pressure and target pressure on the wing surface. This method is based on the assumption that the correction is small, and the residuals of pressure between initial and target surfaces are also small, thus the correction flow satisfies the small perturbation equation. This assumption is reasonable for instance of subsonic and transonic flow with weak shocks. If the residuals are not small enough, they can be relaxed. In order to reach the target pressure the design procedure is iterated several times. As improvement, the geometry smooth procedure is applied and upper wind equation is introduced for instance of supersonic flow occurs on wing surface [5][6]. For () configuration, original airfoils for all section are laminar airfoil, but the analysis shows that the transition takes place much earlier since the three-dimensional flow effect on wing surface. So the target pressure is specified to keep a constant pressure along span-wise in middle part of the wings so that span-wise flow on middle sections is limited and the laminar boundary layer last longer in chord-wise. 2.3 Inlet Design Normally there are more than one engines distributed aside the body of configuration which looks like a bomber. For smaller unmanned aircrafts, the inlet of single engine can be installed inside the body like a fighter, or back on the body at symmetry plane. An embedding inlet of single engine flying wing is designed as a S-shaped channel for air breath jet engine. The geometry of the inlet is also optimized for the best pressure recovery on geometrically confined condition. 3 Results and Analysis 3. Case The VELA is designed by DLR(see Fig. (a)). The detailed investigation shows that the pressure distribution on inner sections is not so good that lift coefficient of inner is quite small, and pressure distribution show that complex intersected shocks appear on up surface. One design purpose is to change the target pressure distribution on inner sections, and let design code change the airfoils. The design results are shown in Fig.3, where the geometry and pressure distribution on inner sections is modified slightly. 3.2 Case 2 A flying wing configuration of is analysed using CFD and the section geometry is redesigned to improve its pitching moment performance, the initial flying wing has good lift-drag performance, but separation flow occurs early on up surface. Thus the pitching moment performance becomes worse at attack angle 8 degree and Mach number.2, which is exactly in taking off/landing case. The detailed analysis shows that the initial airfoil is a typical higher subsonic laminar supercritical airfoil. Although it has good performance at Mach number.6, at the higher angle of attack and/or the lower speed case the cross flow on middle sections is much strong, the boundary layer transition and separation occur very early. In order to solve this problem, the pressure distribution on upper surface is modified to decrease the pressure peak, and the design code corrects the airfoils of the wing sections. The pressure distribution on upper surface becomes more uniform along span-wise (see Fig.3). Finally, the new and initial one are analysed in detail with RANS code. Fig.4 gives comparison of aerodynamic characteristics between initial, designed and the configuration of new with the embedding inlet. It shows that the new configuration has higher maximum lift coefficient but poorer linearity. Fig.5 shows that the new has smaller drag coefficient at lower lift coefficient. The drag coefficient of the configuration with embedding inlet is larger at lower lift coefficient but becomes smaller at higher lift coefficient. It means that the flow separation situation at higher lift coefficient is improved. 2

AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET Fig.6 shows the new configuration of with embedding inlet has the better linearity of pitching moment coefficient, it validates further that flow separation is improved. Comparison of Fig.(a) and Fig.(b) shows that the embedding inlet changes Mach number contour distribution on up surface of, it decreases the Mach number of up surface, specially on inner wing sections. The pressure distribution on up surface is also changed, lower pressure region become smaller near the body (see Fig.3 (a) and (b)). Moreover the lift coefficient of the with embedding inlet is increased compared with alone (see Fig.4), and the drag coefficient is increased at lower lift coefficient and decreased near maximum lift coefficient. To analyse the influence of embedding inlet on aerodynamic characteristics of, the pressure contour on symmetry plane is shown if Fig.7 (a) and (b). It is shown that the embedding inlet increases the pressure of up and low side on symmetry plane section, i.e. decreases the flow velocity on this section. The same results can be observed from streamline distributions (see Fig.8 (a) and (b)). Fig.9 and Fig. show that the streamlines are displaced less due to the inlet. 3.3 Case 3 Although the stall angle becomes greater and linearity looks better than the old due to the embedding inlet, the absolute value of pitching moment coefficient is larger and the static stability is relative too larger for flying wing configuration. Generally pitching moment value is a key design parameter for the flying wing. It must be small enough because thus type aircraft can provide only small control efficiency of pitching moment. In order to decrease the value of pitching moment coefficient, a newer configuration is designed. Its swept angle is decreased so that the pitching moment can be smaller. Meantime the aerodynamic focus may change more slowly with the angle of attack. The geometry and surface pressure distribution of the newer are show in Fig.. The airfoil of every section is the same as case 3, but the platform is different. The design purpose of the newer is to get better performance of pitching moment. The detailed RANS analysis results are given out. Fig. shows the pressure contour and oil stream picture on upper surface in this case. The oil stream is strongly influenced by the embedding inlet. It is evident that the flow separation is improved due to the embedding inlet influence. The aerodynamic characteristic of the newer is given out in Fig.2-fig.4. They show that the embedding inlet increases the maximum lift coefficient and stall angle of the. The drag characteristic become much better at higher lift coefficient, i.e. higher angle of attack, but the ratio of lift to drag is decreased due the inlet influence. 4 Conclusion The numerical analysis is down for 3 configurations in this paper, it shows some interesting results. Body embedding inlet can improve lift characteristic of configurations, it improves separation flow at higher lift coefficient and increases the maximum lift coefficient. Body embedding inlet can increase the stall angle of about 2 degree, thus the pitching moment coefficient curve can last long linearity region. The friction drag of body embedding increases the total drag of +inlet configuration, thus the ratio of lift to drag is decreased. The influence of planform shape on pitching moment performance is quite complex. The stall angle can be enlarged and linearity become better when swept angle become smaller. References [] B. Mialon, ONERA Chatillon, France; M. Hepperle. Flying Wing Aerodynamics Studies at ONERA and DLR Proc of CEAS/KATnet Conference on Key 3

QINGZHEN YANG,YONG ZHENG & THOMAS STREIT Aerodynamic Technologies. Bremen Germany, 2-22 June 25. [2] Richard L. Campbell, Melissa B. Carter and Odis C. Pendergraft, Jr. Design and Testing of a Blended Wing Body with Boundary Layer Ingestion Nacelles at High Reynolds Numbers. Proc Conference of the 43rd AIAA Aerospace Sciences Meeting, AIAA- 25-459. [3] JOSÉ JIMÉNEZ-VARONA.Design of Blendedwing-body Configurations using a Constrained Numerical Optimisation Method. Proc of CEAS/ KATnet Conference on Key Aerodynamic Technologies. Bremen Germany, 2-22 June 25. [4] S. Takanashi: An Iterative Procedure for Three- Dimensional Transonic Wing Design by the Integral Equation Method. AIAA-84 255,984. [5] J. Hua, F. M. Kong, J. J. Wang, H. Zhan, R. Z. Sun, D. Qiu: Recent Examples on Design Aerodynamics for Transport Aircraft. ICAS 22 Congress. [6] Q.Z. Yang, Th. Streit and G. Wichmann. Aerodynamic Design Study for Multi-lifting Surfaces. Proc of CEAS/KATnet Conference on Key Aerodynamic Technologies. Bremen Germany, 2-22 June 25. - -.5 Cp.5 Comparison of Cp () Z/b=. Designed Initial Target.2.4.6.8 x/c.5 Comparison of Cp () Z/b=.263 Designed Initial Target.2.4.6.8 -.2 Fig.2 Pressure distribution on inner sections of VELA (a) Initial fly wing x/c.8.6.4 y/c.2 Figures (b) Designed flying wing Fig.3 Comparison of pressure distribution(ma=.6) (a) (VELA, α=4,ma=.85) (b) (Blended wing body α=4,ma=.6).5 (Ma=.6) old new New + inlet 5 AOA Fig.4 Comparison of lift coefficient (c) (Blended wing body + embedding inlet α=4,ma=.6) Fig. Mach number contour of surface 4

AERODYNAMIC DESIGN FOR WING-BODY BLENDED AND INLET.5 (Ma=.6) old new New + inlet..5.5..5 Cd Fig.5 Comparison of polar curves old new new + inlet (Ma=.6) (b) +Inlet Fig.7 Pressure contour on symmetry plane Cm -.5 -. -.5 -.2 5 AOA Fig.6 Comparison of pitching moment coefficient. Fig.8 Comparison of streamline near the symmetry ( M =.6, α = ) (b) +Inlet 5

QINGZHEN YANG,YONG ZHENG & THOMAS STREIT 2.5.5 -.5 Mach=.2 +Inlet (b) +Inlet - -5 - -5 5 5 2 AOA(Deg) Fig.2 Comparision of lift coefficient 2 Fig.9 The streamline near inner sections.5.5 -.5 Mach=.2 +Inlet, M.2, α = = - -...2.3.4 CD Fig.3 Comparision of drag coefficient.4.2 Mach=.2 +Inlet (b) +Inlet, M.2, α = = Fig. The geometry and pressure contour of case 3 CMZ -.2 -.4 -.6-5 - -5 5 5 2 AOA(Deg) Fig.4 Comparison of pitching moment coefficient 25 2 5 5 (b)+inlet Fig Comparison of pressure contour and oil stream picture on up surface( M =.2, α = 6 ) -5 - Mach=.2 +Inlet -5-5 - -5 5 5 2 AOA(Deg) Fig.5 Comparison of ratio of lift to drag 6