AdR11 lezione 02.ppt [modalità compatibilità] -...

Post on 22-Oct-2018

215 views 0 download

Transcript of AdR11 lezione 02.ppt [modalità compatibilità] -...

UNIVERSITA’ DEGLI STUDI DI NAPOLI “FEDERICO II”FACOLTA’ DI INGEGNERIA

DIPARTIMENTO DI INGEGNERIA AEROSPAZIALE

AERODINAMICA del ROTORELezione 2

A.A. 2011-2012

Paolo CACCAVALE

9 novembre 2011

2

Soluzioni 3D a potenziale ‘semplici’U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

+

doppietta puntiforme

corrente uniforme

3

Pannelli 3D: campi di moto (versori)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

FILE TECPLOT

4

Confronto vortone-doppiettaU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

Vortex ring Vorton),()( trdst v

5

Pannello 3D: sorgente (rif. Katz-Plotkin ‘Low Speed Aerodynamics’)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

formule di Hess e Smith

14

,43

,32

,21

11

14

,43

,32

,21

14

,43

,32

,21

tantan4

),,(

lnˆ

4),,(

lnˆ

4),,(

ji j

jjij

i

iiij

ijji

ijji

ji ij

ij

ijji

ijji

ji ij

ij

zrhem

zrhem

zyxw

drrdrr

dx

zyxv

drrdrr

dy

zyxu

14

,43

,32

,21

ˆˆ

4,3,2,1ˆˆ

4,3,2,1ˆ

ˆˆ

22

22

222

ji

xxyy

m

yyxxd

yyyxxx

jyyyxxx

i

yyxxhzxxe

zyyxxryyyxxx

ij

ijij

ijijij

ijijijij

jjjj

iii

ii

iii

iiii

0,0

2,0

panelPzw

panelPzw

J. L. HESS and A. M. O. SMITH (Douglas Aircraft Company)

CALCULATION OF POTENTIAL FLOW ABOUT ARBITRARY BODIES

6

Pannello 3D: doppietta (rif. Katz-Plotkin ‘Low Speed Aerodynamics’)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

14

,43

,32

,21

ˆˆ

4,3,2,1ˆˆ

4,3,2,1ˆ

ˆˆ

22

22

222

ji

xxyy

m

yyxxd

yyyxxx

jyyyxxx

i

yyxxhzxxe

zyyxxryyyxxx

ij

ijij

ijijij

ijijijij

jjjj

iii

ii

iii

iiii

14

,43

,32

,21

2

14

,43

,32

,21

2

14

,43

,32

,21

2

ˆˆˆˆ4),,(

ˆˆˆˆˆ

4),,(

ˆˆˆˆˆ

4),,(

ji jijijiji

jijiij

ji jijijiji

jiij

ji jijijiji

jiij

zyyxxrrrrrryyxxyyxx

zyxw

zyyxxrrrrrrxz

zyxv

zyyxxrrrrrryz

zyxu

formule di Hess e Smith

2

,0 panelPz

7

Pannello 3D: doppietta (rif. ‘Program VSAERO Theory Doc.’ NASA CR 4023)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

Biot-Savart equation:

8

SingolaritàU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

9

CertificazioneU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

10

Induzione di un segmentoU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore c vettori A e B (P-A, P-B)

AX = PX-XAAY = PY-YAAZ = PZ-ZAA = SQRT(AX*AX+AY*AY+AZ*AZ)BX = PX-XBBY = PY-YBBZ = PZ-ZBB = SQRT(BX*BX+BY*BY+BZ*BZ)

c vettore S (edge B-A)SX=XB-XASY=YB-YASZ=ZB-ZAS=SQRT(SX*SX+SY*SY+SZ*SZ)

c prodotto vettoriale AxBAVBX = BZ*AY-BY*AZAVBY = AZ*BX-BZ*AXAVBZ = AX*BY-BX*AYAVBS = SQRT(AVBX*AVBX+AVBY*AVBY+AVBZ*AVBZ)

c prodotto scalare A.BADB = AX*BX+AY*BY+AZ*BZ

c vettori S e A lungo gli assi del pannello L e MSM = SX*CMX+SY*CMY+SZ*CMZSL = SX*CLX+SY*CLY+SZ*CLZAM = AX*CMX+AY*CMY+AZ*CMZAL = AX*CLX+AY*CLY+AZ*CLZ

A

B

P

y

x

PA

PB

11

Induzione di un segmentoU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore c doppietta unitaria

SCALE=0.

IF (ABS(A*B*(A*B+ADB)).GT.EPS) THENSCALE = (A+B)/(A*B*(A*B+ADB))

END IF

VXD = SCALE*AVBX/(4.*PI)VYD = SCALE*AVBY/(4.*PI)VZD = SCALE*AVBZ/(4.*PI)

c potenziale

IF (ABS(PNN).LT.EPS) CJK1 = 0.0IF ( PNS .LT.EPS) CJK1 = -TWOPIIF ( PNS .LT.EPS) UNS = -TWOPI

DIC = (CJK1 )*/(4.*PI)SIC = (RJ31-PNN*CJK1)*/(4.*PI)

c sorgente unitaria

RJ3=0.IF (ABS(A+B-S).GT.EPS.AND.S.GT.EPS) THEN

RJ3=LOG((A+B+S)/(A+B-S))/SEND IF

ALL = AM*SL-AL*SM

PA = PNN*PNN*SL+ALL*AMPB = PA-ALL*SMRNUM = SM*PNN*(B*PA-A*PB)DNOM = PA*PB+PNN*PNN*A*B*SM*SM

CJK=0.IF (RNUM.NE.0.) CJK =ATAN2(RNUM,DNOM)

CJK1=CJK1+CJK

ULS = RJ3*SM/(4.*PI)UMS = RJ3*SL/(4.*PI)UNS = CJK /(4.*PI)

VXS = ULS*CLX+UMS*CMX+UNS*CNXVYS = ULS*CLY+UMS*CMY+UNS*CNYVZS = ULS*CLZ+UMS*CMZ+UNS*CNZ

12

Uni

vers

ità d

egli

Stud

i di N

apol

i “Fe

deric

o II”

Dip

artim

ento

di I

ngeg

neria

Aer

ospa

zial

eA

erod

inam

ica

del R

otor

e

FILE TECPLOT

13

Core RadiusU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

FILE EXCEL

14

Core RadiusU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

FILE TECPLOT

15

Influence coeff. of a rectilinear panel (Katz-Plotkin)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

16

Influence coeff. of a rectilinear panel (Katz-Plotkin)U

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

17

Dal filamento vorticoso al vortoneU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

18

Induzione di un vortoneU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore c vettore S (edge B-A)

SX=XB-XASY=YB-YASZ=ZB-ZAS=SQRT(SX*SX+SY*SY+SZ*SZ)

c vorticità edge/vortoneGAMMAVOR=GAMMAEDGE*S

c versori edge/vortoneCAX=SX/SCAY=SY/SCAZ=SZ/S

c induzione vortonePX=PXI-PXVPY=PYI-PYVPZ=PZI-PZV

PNS=SQRT(PX*PX+PY*PY+PZ*PZ)

VXV=-(+PY*CAZ-PZ*CAY)/(PNS**3)/4./PIVYV=-(-PX*CAZ+PZ*CAX)/(PNS**3)/4./PIVZV=-(+PX*CAY-PY*CAX)/(PNS**3)/4./PI

A

B

P

y

x

PNS

19

Conversione di un pannello in vortoniU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

20

Uni

vers

ità d

egli

Stud

i di N

apol

i “Fe

deric

o II”

Dip

artim

ento

di I

ngeg

neria

Aer

ospa

zial

eA

erod

inam

ica

del R

otor

e

FILE TECPLOT

21

Sistema lineareU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

FILE EXCEL

22

Core RadiusU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore

FILE TECPLOT

23

Riferimenti bibliograficiU

nive

rsità

deg

li St

udi d

i Nap

oli “

Fede

rico

II”D

ipar

timen

to d

i Ing

egne

ria A

eros

pazi

ale

Aer

odin

amic

a de

l Rot

ore J. Katz, A. Plotkin

‘Low Speed Aerodynamics’

McGraw-Hill

B. Maskew

‘Program VSAERO Theory Document’

NASA Contractor Report 4023

L.E. Browne – D.L. Ashby

‘Study of the Integration of Wind Tunnel and Computational Methods for Aerodynamic Configurations’ (PMARC)

NASA Technical Memorandum 102196

D.J. Willis, J. Peraire, J.K. White (MIT)

A Combined pFFT-Multipole Tree Code, Unsteady Panel Method with Vortex Particle Wakes

43rd AIAA Aerospace Sciences Meeting and Exhibit

FINE LEZIONE 2