Küchemann, D.

Weber, J.

German

Institute of Transport Research:Publications

~~~----~---------------(' • l ( i Aerodynamische ~ Versud\san-stol~ Göttingen e. V. · A~tell11n~: Vfi S S ,1 B•richt: :s 4~/ A/ 04 . A'uftt99JLA 113 Detli!ft: 30. '· 42 •' . . _.,. Üb;e,:r die · St:römu.l,).s~ g .:rt~PSfpt:m:i&en Ve:r~le~p~aen Vl:.M:ltteilmtg~ .Wet t.~r~ Mese,un~ep.~ a». ,!l,itlJ;@iutae:rä'lrea · Über$ioht:Der vorliegende :Seri:oht'bringt~ als Ergänzung zur !V .Mitteilung. (FB 1236/4) wei te1·e · Windkanalmessungen .• an Einlaufg~räten, die vor allem <den Einfluß de.s N'asen:rad:i,ua der:Ve:rkleidungauf d:Le St:römungaver-häl tnisae zeigen eoll.en. Aus~e;rdam wil:d~ eil'le ein-. faohe Regel zum Entwurf soloher An:ordrtungen, b.ei denen sich e.ine Nabe itn . Innenraum befinde·t,vorge-schlagen. Gl;iederu,ng: !.Vorbemerkungen. II .• Ergebnisse .. l.Eine Verkl.eidu.ng mit großem Nasen.radius,. 2.Verklei.dungen mitve+acbiedener Naeenfortn. II I. ~use.mmenfasaung. . .Der Bericht mnfasst: l5 Seiten , Der Leiter: Die )3earb~i ter: -AVA Nr. 101. 5. 39. 5000 ,, - 2-I.Vorqemerkungen. Der vorliegende Bericht· bringt wei·tere Windkanalmessungen an Einlaufgeräten als Nachtrag zur IV .. Mitteilu.ng (FB 1236/4). Bei den früheren Messungen waren noch einige Fragen offen ge-blieben .. So war der Einfluß des Nasenradius und der Nasenform der Verkleidung nicht eingehender behandeltrdieser sollte in einer neuen MefJreihe näher untersucht werden. Wei terhdin soll ein ... , .. einfaches Entwnrfsverfahren beschrieben und geprüft werden. Die friiher untersuchten Verkleidungen hatten alle den gleichen Krümmungsradius 8N = 0,067 Ra l).Auf Grund dieses kleinen Nasenradius hatten sich bei einigen Verkleidungen vor. allem am Stand ·innen am Einlauf lcleine lokale Ablöseerscheinungen gezeigt"'~d es war vermutet,daß durch größeren Nasenradius die Strömung im Innenrau.mnoo:Q,;ge:ringere Verluste zeigen Würde.Da-her wu..rde bei der :neuen Verkleidung 5 ~ N = O, 095 Ra gewä.hl t. Um einen besseren ttberblio:tr über den Einflu.13 von Nasenradius und·N'asenf'orm zu bekommen,v.rurde der Nasenteil von Verkleidung 5 noch abgeändert (Verkleidung 6 und 7 ;vgl.Abb.lO) ;.sodaß sioh kleinere Nasenradien 9N = 0,040 Ra 'bzw. 9N ::: Ot014 Ra er-gaben. Bei. den früheren Messungen hat sich wmiterhin gezeigt,daß es zur Vermeidung vonStrömungsverlusten sehr wichtig ist,daß sich die Querschnitte des Innenraumes zumindest nicht erweitern. l)Die Bezeichnungen sind die gleichen wie in der IV .. Mitteilung. '•.·. - 3-Die Aufgabe besteht also da.rin;eine Anordnung zu finden,bei der die Außense.i te im Schnellflug die gewünschten niedrigen Übergeschwindigkeiten besitzt und der Innenraum sich nicht erweitert.Bei den Anwendungen ist der Innenraum stets tei;t.-weise oder ganz ausgefüllt z.B. durch eine Nabe oder einen Kühler (Sternmotor).zurL(;)sung der strömungsphysikalischen Auf'g~be kann man die Verdrängungsktlrper.duroh Singu.laritä-ten ~-rsetzen.Ma.n kann sich jedoch die Aufgabe wesentlich er-leiohtern,wenn man die aus den Messungen hervorgehende Tat-sa.che·lieranzieht,daß die Strömung auf der Außenseite in der Hauptsache nur von der Durohflu.ßmenge (Staupunktslage) ab-hängt und nicht von den Verhältnissen im llUlenraum,sodaß man die 'Gestaltung des Innenraumes unabhängig von der Außenseite behandeln kann.Der.Einbau von Kühlern ist~iesem Sinne bereits untersucht 2 ) .Der Fall einer im Inneren eingebauten IIJ'abe soll im folgenden besproohen.werden. Bei der Umst:römung einer Nabe liegt die Schwieri_gkei t bei der Konstruktion eines z"B.konstanten Que:rsohriittsver-laufes vor allem darin,die nstrtlm-imgsquerschnitten selbst zu bestinunen,'besonders in dem Ge'b.iet,wo die Nabe beginnt. Wir betrachten dazu die freie Umstrtlmung der Nabe,welohe für den einfachsten Fall des räumlichen Halbkö.rpers ( Einzelquelle) .in A'bb.l dargestellt ist. (Anders geformte Naben kann man durch geeignete Quellsenkenanordnungen ereetzen,wa.s jedoch 2lvgl.D.Küohemann und J.Weber:Der Einbau von Kühlern.Erscheint demnächst als FB .. 1 -·'-Nr. ~; :'~ / 13 / te ,/ it / 10 -:: g / -8 / 7 -:"' 6 ---5 _,"' ---/ --1 - -'\. .,, '\. .. 1\ \. \. ~ ·\ - f.-----\ ---I--- l---\ :----~ -_.-:. ~ f..--- ~ ~ -- V - ~ ~ [-\ / 'V I J!1 1,8 L: /. ·L ' / / / --/ 1--~ / - ~ ---1---- / -+-----/ -~ :---7' -J-----' ~ ---d---:2 r:=- ' ~ -,_ -------'-, 1 v---iL_ q6 -r RN t qE 1;,. ~.12 ~10 -08 -06 ·-O'r. 0 --oo- I . I. I. -cz I. 0 . 02 04. Ofi 08- .1{_1 c 00 I . Nr.O 1 ·,: 3 4-g f'r ~-J 1/ ~k ~ _.-:' ~ !"-' -I I / ~ r I I :' .. . -:. -V ~: \ "\ '\. ""\ .\"':'.·· ·;, :; . , .... \ ~ ~ -- ~\ ·--::::--~ ~· ,_" 1.0 ---~-- ·- ··- ---'--'-- -'---p-pc . T . f -----1--~-- --~----qs of -· 0.6 0,5---~-- -----0,4-0,3 0,2 . 0,1 X 0,2 0/t 0/) 08 :- RN1, 1,'r -~2 -~0 -0,8 '-0,6 -0,4- 1\~ ~ t.::::: -_\ ~ ~ - :::::-:::. [}: I'--~ ~ -\ I'----~ ~ ~ -q ........ I---·Oir Abb. 1. Stromlinienbild der räumlichen Halbkörperströmung und DruckverteiLungen aüf · . den Stromlinien. Linie I: Beginn des Druckabfalls auf den Stromlinien; I-inie j[: p•p0 ; Linie fil: Ende des Druckabfallgebietes. - 5 -für den freifahrenden Körper keine große Reohena:rbei t bedeu-tet. )Die Druckverteilungen auf den Stromlinien· (Abb.l) zei-gen,daß in·dem.Gebiet zwischen den Linien I und I!I Druckab-fall hel"rsoht ,d .h. do.ß sich die Ströraungsquerschni tte in die-sem Gebiet verengen"S.troma.ufwärts vo:il Linie I und stromab-wärts von Linie III ergibt si~h eine Erweiteru..ng.M.acht.ma.n eine dieser Stromlinien zur Innenwand der Verltleidung,wie es ~ in Abb,2 an eip.em Bel-spiel gezeigt isttso kann man vermuten, daß $1oh die Strömung an dieser Wand nicht ablöstoWill man .je~ooh ~.:a. einen ko~stanten oder sichstetig verengenden ' Querschnitt er.balten;so sind n~rnoch geringe Abänderungen . . . .. voti der vo~gesohlägenen Innenseite nötig, die leicht· durohzu-fühxeil s.ind. · II.E~rsebn.isse. 4tl4!ai Yerl$le;t§.ung m;L t g:rgßem !IQ.Senradi)l.s •• Die V.erkleidung 5 (AblJ~3J .. wurd~ ·als une:n.dlich langer Ringkörper ohne Nabe, waloher.}1in Geschwindiglteitsverhä.ltnia u1/U0 ~ o,; herstellt (Q/'ii'Ra2l.l0 = Otl6),bereehnet .• Die zugeh5rige Skelettlinie er-gab sich aus einer längs der Tiefe konstanten Wirbelbelegung, . wobe:i mit Hilfe eines Iterationsverfahrens erreicht wurde.daß Belegungafläche und Stromfläche der Gesamtströmung zusa.mmen-fallen.Duroh Überlagerung einer Quellringströmung ergab sich die in .A.'bb.:; eingezeichnete Kontur.(Das Hechenverfahren soll ausführlich 1n einer demnächst erscheinenden. Mitteilung be-schrieben werden.)In den Innenra.um wurde die vorhandene Nabe 1 1,0 Stromlinie der Halbkörperströmung 0,+ q2 0 0,1 02 : I . 6-0,5 0,7 0,8 A.bb.a. ·ZurFormgebu~g des ln~enra.umes . . 0 O,f O,t;. 0,5 q7 0,8 0,9 verkteia1,mg 3 Halbkörper 1,1 12 •I (2 Abb.3. Vergleich der vermessenen und gerechneten Anordnungen. 1,3 1,5 ~..L ·. Ra 1,3 --..5.. f5 Ra - 7-gebracht und die Innenseite der Verkleidung in der oben an-gegebenen Weise abgeänd~rt.Dabei war zu berüoksiohttgen,daß die I~abe l nicht genau die Form eines Halbkörpers hat.Es_ wu:r-de erstrebt •. daß die Quersohni tte gleichmäßig abnehmen, um Strömungsverluste im Innenraum zu ve:r;meiden.Bei Anwesenheit der Nabe wird das Geschwindigkeitsverhältnis tr1/U0 _.bei wel.:.. chem die Verkleidung stoßfrei angestr8nrt wird,von o, 3 ver-schieden sein. Wir erhalten dießelbe D.urob.flußmenge ~/'Jf Ra 2u0 = 0,16 wiebei der leeren Verkleidung"wenn wir bei der neu-en Anordnung 51 bei x = 2,5 Ra (vgl,d:ie Zal:ilentafel) das Ge-. schwindigkeitsverhältnis u1/U0 = 0,_85 haben .• Zahlent~f'el x;a·.; .a 0 0,05 O;l· 0,2 o,; 0,4 0,5 r/R aus~en 0,600 0,695 o,12e 0,772 ö 802 t .· .. 0 826' ' . 0,845 r/Rän o,6oO 0,516 O;c503. 0,501 0,515 ._0,539 0,565 inn , /R x a 0 6 , . 0,7 0 8 ' 0,9 liO· l;l 1,2 r/R aust!en 0,863 .0,879 0,892 o,go:; 0·913 ·' ~,922 0,931 r/R~ -o, 59"3. 0,621 0;645 0,66? 0,682 0,696 Ot~708 inn x/Ra 1,3 1,4 1,5 2-o 2.-.? ;.,o r/H. . aus~en 0,939 o;946 0,953 0,978 0,992 1;,000 r/R · inn~n Ot718 0,726 0,733 0,758, 0 776 , . 0,821 -8-. ';.' ;- -_;' ' . . - . - '< .. - :· ..... ·.. -.' .. J. .-·: ;-·. tn Abl.h4 'sind die g$messenen ·. :;) und gerechneten Wandga ....... . ·• eohwi:ndigkeltsverteilungen ~Ur den Entwurfszustand miteinander. '• .< I .·• ' • • • • • sei,~e ##d ~m beaonderen;daß die Theorie die ma:1rimale tlberge-· .. schwin~igltei t hinrei()h~nd genau wiedergibt.Dieae. beträgt im ·: ' -, . . . .'- ~- -- . : . . . _·. . . - . . . . vor1J..eg$nden FE:tlle 22% • Geringe Abweichungen zwischen Theorie 0 ' .' > • • I •, :' • •o • 0 ·~ :• ', ' ' . l.m~ 1\'leesung in dem Debiet 0•1 ' x/Ra ·~ 0,3 sind a.~f ei11e He;r;... . stell~saU:ngenauigkeit beim Modell zurüokzu:fü.hren.Da der Innen-·. ratinl durch. d.as Hinzukommen der Na.be beträchtlich abgeändert ist,lln.terscheidet sieh nun auch O-ie gemessene Wa.ridgesohwindig-'5) Die l\IIessu.ng~n wurdenwie in der IV.Mitteilung durchgeführt. Aussenseite 'C'a:O [j- Gw flw Drosselstellung J 7I; a;; 1R:U0 ® 1 0 - iJ 0 3 0.35 0.95 ao65. • 4 0.60 095 0,t2 . 0 s {?84- . ä95 ots ·. G) 6 .1.08 0.95 0.21 e 7 !27 0,9+ OZit-. () 8 1.78 093 Q54 () g 3.12 ({92 0.59 Anordnung 51 - tO ~--'-------'-------'----'-------'--~-'------' Abb.S. Druckverteilungen auf der Aussenseite bei verschiedenen Betriebszuständen [ui gemessen bei x·Zß Ra]. Abb.fi. Druckverteilungen auf lnne'!seite und Nabe bei verschiedenen Betriebszuständen. ~0 P-Po T ~ I 'f q. -·~···· ' .. (). ' ' q. ? · .. _I E'. I·. l ~-. r _,_,_-r~-l\ -1. ~ \ .. .· _"i r -=-1 ~ ~ li ~\ .. ·. _.",_ ' . Ir\ l\ ' .. t ~ I I~ • 1. I I . I I ' -·--1 ---= =-- J --... qc ..... t\6 qa ~0 12 f/r tG 1,8 2fJ _....,...- ....K. 2,~ ~- 'I '! 0 I I , R. I l . . i I i a I ,. i . ! J I i -~ ~ I~ I i I i J l ~~ ~o--e---e. ~ j i l ' i .,~ ~ I I ~ I i I ' ! . ~ I ; \ I j I T-~ I I \ I I r ' I I I i I ~ Ci I Anord~ung 51 I ! .~ I l .. ,, ~ I •. I \ .-..:. "· ........... ' l ·-,q. -o. -.o.. -o. ~tl -1, -f. 1\ -~ ,, ~ k:: 1--"' \ --'- ., r ~ ---~, ~~ I -· Innenseite ........ r ..... t...... ~'-----'---Nabe . , ..... ........ ~ r. ..... -c•o ....... I'>· ~- lii flw Gw ~ Drosselstellung lio ()tb rR:uo .:. @ 1 0 - 0 •[-'---0 5 0.8+ 0.95 0,16 e 7 t.Z7 o,g9- O.Z'fr. j () 8 1.78 0,93 C?31r -2,6 1 ((:) I ; ,; - 10-keitsverteilung auf der Innenseite entsprechend von der theo-retisch (ohne Nabe) berechneten.Auf der Innenseite ergibt. sich sehr bald die gewüneohte,nur mehr schwaoh·zunehmende Geschwin-digkeitsverteilung.Auch die Geschwindigkeit an de:r Nabenwand / fügt sich diesem Gang ein,d.h.die Geschwindigkeitsverteilung über den Strömungsquerschnitt ist für x ~ 01 5 Ra praktisch ·konstant .Damit hat sich die vorgesehlag.ene Entwul:'fsregel als brauchbar erwiesen,au.oh in diesem Falle·,wo die Abänderung des Innenraumes groß ist. Die VJan.ddruekverteilungen bei nioh'tstoßfi'eiem Eintritt sind in den Abb~ 5 und 6 aufgetragen.Abb.7 zeigt die Druokve:r-; ,·. J-: .... <·-:c.:~, Abb.7, . Wanat;Jtuc~verteilung~haf(! st~~~~·.vir9Je/cn·"c ····••··. zweier. Vetirteicito1gen mit. v.~rsifhi~denen' ·. < .-•.. "' · · · .. f'lasenra.dlen /vgt.Abb.-3. · ' ' · · , ,. . · ·. ·". -:··-',! - ll-teilungen am Stand. Von be$onderem Interesse ist nun der Ver-gleich mit Verkleidung l ,v.,~lohe sioh von der Verkle;tdung5'in deJ· Haup:haaohe nur durch den l'fasenradiua u.nterach$idet (vgl. Abbil :3) .Der Einfluß au;f die U'bergeschwindigkei·hen auf d~r · · ·. Außenseite ist nur· gering (vgl • .A.bb.s),Einen we.sentliohen Un-, teraohi.$d bemerkt man jedooh bei den Druckverteilungen .· a:m. St.and {Ah.b. 7) .. Trotz des etwas kleineren Eintri ttsquersehnittes .· . . -. '. . b~:t. Ver~leid1lng 5 ist die .Unte~druckspi·tze ·auf der Innenseit·e . . . . bei. weitem niob.t niehr so groß.AUoh d~e lokale Stramungah'bltia'l!ng . im .. vord.er~li Teil d.er Innenseite bei Ve·rkle.idung l ist bei Ver-. klei(}Ul').g· 5 kaUm no.ch feetßt.ellbitr. (vgl.auoh Abb.l$·). • Um· die Güte. d.er Umström®g 'zu kentw@iqhnen •. wurde in der IV,l\!litteilung .~die Qrtiße Q./Qth · .ge:brauoht.piese stellt das v~rhä.ltnis zwisoh~n· . der .Wirld.:iohe~· und d~r :potentialtheoretischen Durchflußmenge .· · . da.r.pa.bet ist ·jEHiocb · ~u bedenken,da.ß die unvermeidliche Wand.;. . reibungm?.c.hieht bei kleinen st:römungaquerschnitten. stärker ins . · · · Gewia:nt· fäll-t· a.lp bei gr~ß~:ren. so ·nimmt die verdrängungsdicke · det {}:ren.zeohiQht bei .der Alto:rd~1mg 51 einen größeren !Ce:Ll des · · .. Queraqhn.ittes ein als bei der Ano:rdn'Ullg 11 .• :Bei unseren Mod.e11 .... . : . · abinessun$en ist· die _Grösse . Qv/Qth, a.~so kein. ~indell,tig~s .· . . . . ~tlternJ~;"kma.l,wa.s ~.B. beim Vergleich der in Abb,7 (und i;) .an- · gegeb.enen Werte zu beachten ist. In Ab.l5~9 sind die durch Integration aus den Wa.nddruokvet• · teilu.ngentgewonnen~n Axialkräfte mit d.en theoretischen. We:rten ;, - 12-(vgl. die IV.Mitteilung) vergliohan.Die Anzahl der Druokm~ß ... stellen reichte ·j edooh im rückwärtigen Teil des Innenraumes · nicht aus »um in allen Fällen einwa.ndfr.eie E:rge'brlase zu erhal-ten. Dieser Mangel macht .eich aber erst bei gr<ißeren Durch~ flüssen (u1 > U0 )in nennenswertem Mli.ße bemerkbar.Bei gerin-gen Durchflüsse~1 ist die Übereinstimmung· zwischen Messung und Theo;r-ie recht befriedigend.Es zeigt sioh,daß die bei der theo,... retischen Aufteilung der Axialkräfte -in Nabenschub und Verklei-dungsschub gemachte Annahm.e,die Nabe befinde sich sehr weit im Inneren der Verkleidung selbst in unserem Falle brauchbar ist • . t:-· . 0,8 ; -·. ··,, . ::-.-·· .. · ' .. -, Abb.!i · iJem;$.r;eh{iJnri gere~h#'et~·~;rldtlftlift~.· ·· j - 13-Das Verhalten bei Schräganströmung i$t gegenüber An-ordnung 11 praktisch nicht geändert,weshalb nicht näher darauf eingegangen zu werden braucht. ·· a._,yerk;t.eid:tmßen pti t ver.s_ghiedener· Nasenform.Die Ver-. kleid·ttngen 6 und 7 ,.deren Form man aus Abb.lO entnehmen ) kann,wurden durch Aufsetzen eines Holzringes a'Uf Verklei-dun;g 5 hergestellt ,Da diese ~~essungen nur zur Obersicht dienen sollteh;wurde auf die vorderen Druokmeßatelle.n ver-zichtet.Zunächst zeigt aioh,daß die·ttbergeschwindigkeiten auf der Au~enseite wahrscheinliöh verrin~ert werden (Abb. · ll).Andererseits ergibt sich offenbar bei Anordnung 71 für Durchfluß Null ein lokales Abreißgebiet an der :Nase,wie aus der Grenzschichtmessung an der Außenseite (Abb.l2) hervor• . . . . •. geht.Wesentlioh verschlechtert werden jedoch die Strömungs-verhältnisse im Innenraum bei größeren Durchflüssen (Stau.-. pu.n.kt auf der Aui~eneei te), besonders am Stand. Strömungsbe-obachtungen mit Hilf'e.elner Fadensonde zeigten bei Ariord..;. nung 61 und 71 ,daß die Strömung au.f der mnnenseite tinmit-telbar hinter der Nase abriß.Dieses Verhalten findet man auch in den Geschwindigkeitsverteilungen über einen Innen-qu.erechnitt .(A'bb .. l3) wieder. Beim Vergleich der Anordnungen 51,61 und 71 geben die Werte von QwfQth einen Anb.alt über die Strömungsverluste.Aus Abb.l3 geht weiter du.rch Vergleich der Anordnungen 51 und 4 leer hervor,daß sich die v_erluste aucl1. bei 51 nur auf die Wandreibung beschränlten.Daraus folgt, daß bei Anordnungen,bei denen sich unmittelbar hinter dem o. 'Ahb.tO. OleVer/aeidungen 6und7ats Abänd,erungen des Nasenteiles der Verkleidi.mg s ... 1,02 ·. Ano.ranung · • ' 11··. . 0. '51 0 61 Ä . 71 ... -• Abb. 12. Geschw;nqigkeitspro(ile im Aussenraum für verschiedene Anordnungen bei x -.c;s Ra für · Durchfluss Null. \ ~Or---,---,---~~-----,r---.---. P-P, Aussenseite ~ . ~-o 'lo --------- Anorcinung ' 0 51 Tq6 o · et r[mm] .:1 71 -'\111!0- J!. . tO Ra. flr Abb. · ff. Einfluss der Nasenform auf die Druck-verteilungen auf der Aussenseite. ~f~~~-~~--+---r•~~ j ~~~~-~~~-+~ I am.Sti:md · aw . a,h 6 51 q92 0 61 v.89 A 71. Q84-• 4-teer 0,98 q6 0.8 tO lJw -.;...-uth Abb.1.3. Geschwindigkeitsprofile im Innen-raum für verschiedene Anordnungen bei x • 2,5 Ra am Stand. - 15..:. Einlauf. eine Nabe befindet,der gewählte Krümmungsradiu.s aus-reieht,um auch am Stand eine verlustlose Strömung zu gewähr-leisten. II I ... Zusa.mmenf'?-ssu.ns. Während die in der IV.Mitteilung untersuchten Ringkörper den gleichen Nasen11 adius hatten, wird in diesem Bericht eine Reihe von Messungen an Einlaufgeräten mit verschiedener Na.-senform vorgelegt .Die Messungen zeigen,, daß es möglich ist, durch geeignete Wahl des Nasenradius und der nasenform bei Anordnungen mit einer Nabe im Innenraum auch am Stand eine vollkomll\.en verlustlose Strömung ,(abgesehen von der Wandrei...: bung)zu erreichen. Eine spitze N·ase verschlechtert die Strö-mungsverhältnisse im Innenraum bei großen Durchflüssen erheb-lich. Es wird ·weiter vo:r•geschlagen, den Entwurf' von Anordnungen mit Nabe dahingehend zu vereinfaöhen.daß die Gestaltv.ng des Innenraumes unabhängig VOll der Außenseite beha.ndel t wird.Di.e Wände des Innenraumes si,n.d so zu formen,daß die Strömungs .... querschnitte konstant bleiben,bzw.leioht; abnehmen.Zur Lösung dieser Aufgäbe wird das Strömungsbild der frei unströmten Nabe herangezogen .. Die Brauchbarkeit der vorgeschlagene11. Ent-vvurfsregel wird durch die Maßergebnisse bestätigt. Die tibergeschwi.tidigkeiten auf der Außenseite im. Schnell-flug sind hinreichend klein und 'Werden ftir sto!Jfreien Ein-tritt durch die Theorie gut wiedergegeben. I 

Über die Strömung an ringförmigen Verkleidungen. Teil 6: Mitteilung: Weitere Messungen an Einlaufgeräten

Abstract

Similar works

Full text

Available Versions

Institute of Transport Research:Publications