Pervanenin girdap kafes yöntemiyle aerodinamik analizi
Aerodynamic analysis of a propeller by vortex lattice method
- Tez No: 389183
- Danışmanlar: PROF. DR. MAHMUT ADİL YÜKSELEN
- Tez Türü: Yüksek Lisans
- Konular: Uçak Mühendisliği, Aircraft Engineering
- Anahtar Kelimeler: Belirtilmemiş.
- Yıl: 2015
- Dil: Türkçe
- Üniversite: İstanbul Teknik Üniversitesi
- Enstitü: Fen Bilimleri Enstitüsü
- Ana Bilim Dalı: Uçak ve Uzay Mühendisliği Ana Bilim Dalı
- Bilim Dalı: Uçak ve Uzay Mühendisliği Bilim Dalı
- Sayfa Sayısı: 91
Özet
Bu çalışmada pervanelerin aerodinamik performanslarının nümerik yöntem ile belirlenmesi amacıyla girdap kafes yöntemini temel alan ve FORTRAN dilinde yazılmış bir bilgisayar programı geliştirilmiştir. Program girdi olarak verilen pervane geometri bilgisini, pervane ilerleme oranını ve pervane palalarının kesit profil bilgilerini kullanarak hesap yapmaktadır. Literatürde girdap kafes yöntemi ile yapılmış çok sayıda çalışma mevcuttur. Bu yöntem özellikle tasarım aşamasında oldukça hızlı ve doğru sonuç verdiği için tercih edilmektedir. Ancak girdap kafes yöntemleri taşıyıcı elemanı kamburluk eğrisinden ibaretmiş gibi varsaydığından viskoz etkileri hesaba katamamaktadır. Çeşitli kanat profilleri ve test pervanesi için programda kullanılan yöntemin güvenilirliğini ve doğruluğunu göstermek için birçok uygulama yapılmıştır. Çalışma sonuçlarını değerlendirmek amacıyla daha önce geliştirilmiş olan taşıyıcı çizgi teorisi ve pala elemanı-momentum teorisi programları kullanılmıştır. Kanat üzerindeki doğrulama çalışmaları taşıyıcı çizgi yöntemi programı (MAYwingPLL) ile gerçekleştirilmiştir. İlk uygulamada düz levha şeklindeki kanat ve taşıma eğrisi eğimi düz levhadan daha küçük olan kamburluksuz hayali profile sahip kanat için analiz yapılmıştır. Çeşitli hücum açılarında kamburluksuz profil sonuçlarına göre, taşıyıcı çizgi teorisi ve girdap kafes yöntemi elde edilen taşıma katsayısı değerleri birbirine çok yakındır. Bu değerlerin analitik olarak hesaplanan düz levha ve hayali profil verileri ile karşılaştırılması sonucunda analitik çözüm ile aynı hücum açılarında daha yüksek taşıma katsayısı değerleri hesaplandığı görülmektedir. Çünkü üç boyutlu kanat iki boyutlu kanattan farklı olarak taşımanın yanında indüklenmiş sürükleme bileşenini de içerir ve üç boyutlu kanatta iki boyutlu kanat ile aynı hücum açısında taşıma daha azdır. İkinci uygulama olarak, sayısal yöntemlerin kamburluğa ne tepki vereceğini göstermek için aynı hayali profile bir kamburluk verilmiştir. Taşıyıcı Çizgi Teorisi ile elde edilen taşıma katsayıları, profil kamburluğunu dikkate alacak biçimde bir değişime uğrayarak, Girdap Kafes Yöntemi ile elde edilen taşıma katsayılarından daha büyük değerler almaktadır. Üçüncü uygulama olarak gerçek bir kanat profiline ait doğrusal olmayan bölgeyi de içeren veri kullanılmıştır. Kesit profili olarak NACA 4412 profili seçilmiştir. Her iki yöntemle (Taşıyıcı Çizgi Teorisi ve Girdap Kafes Yöntemi) elde edilen sonuçlarda doğrusal bölge içinde taşıma katsayıları literatürden alınan taşıma katsayısı değerlerine benzer şekilde lineer bir artış göstermektedir. Her iki yöntemde de viskoz etkiler hesaba katılmadığı için, kanat profilinin doğrusal olmayan bölgesinde taşıma eğrisi eğiminde azalma olmamıştır. Pala test çalışmalarının en sonuncusunda deneysel çalışmalarla bir karşılaştırma yapılarak yöntemin genel başarısı hakkında bilgi edinilmesi amaçlanmıştır. Bu çerçevede dikdörtgensel üst görünümlü kanadın kesit profili olarak NACA 4415 profili seçilmiştir. Bu uygulamada dikdörtgensel üst görünümlü kanadın açıklık oranının 6, 9 ve 12 değerleri için hesaplamalar yapılmış olup, bu sonuçlar Ostowari ve Naik tarafından verilen üç-boyutlu deneysel sonuçlarla karşılaşmalı olarak sunulmuştur. Açıklık oranı değeri küçük olduğunda (AR=6) Girdap Kafes Yöntemi sonuçları 0°
Özet (Çeviri)
Propellers with different shaped blades are used in many applications like fans, ships, wind turbines, aircrafts and helicopters. In this study, aerodynamic performance of a propeller is analysed numerically by using vortex lattice method. A FORTRAN based computer program is developed for this project. Propeller geometry information, propeller advance ratio and propeller blade section characteristics is used as program input for numerical calculation. There are many studies in literature which used computational fluid dynamics (CFD) methods, experimental methods and vortex lattice method (VLM) in order to calculate propeller performance in fans, ships, wind turbines, aircrafts and helicopters. There are numerical methods like momentum method, blade element method and Prandtl lifting line method in propeller theory. These methods are also used in literature for propeller calculations. According to studies in literature, thrust and tork values calculated by blade element method are generally %10 greater than values calculated by experimental methods.This difference is resulted from constraint of this method on considering effects of trailing vortices and blade vibration effects on propeller performance. There are also studies in literature comparing blade element method and CFD method. Blade element method and CFD thrust, tork and efficiency results are very similiar to each other and distribution graphics behave in a similiar way. CFD method shows tip loss caused by tip vortices better than blade element method. Because of this, there is a significant reduction in CFD profile's tip region. CFD methods can offer better opportunities for fan performance development rather than old design and analysis technics. CFD programs solve Navier-Stokes Equations for whole flow region. Programs written by VLM method can only calculate velocity and pressure on required location. In slender blades, successful results can be obtained by VLM method by using less computer force. Vortex lattice method helps selecting appropriate blade geometry from different blade geometries and helps rapid pre-design. There are many studies in literature which used vortex lattice method (VLM), in order to calculate propeller performance in ship, aircraft and helicopter propellers. However, there are few studies in literature which are analysing fluid flow around axial fans. Vortex lattice method is selected by researchers because it helps getting rapid and accurate results for propeller design. The weakness of the method is that lifting element is only considered as a camber curve and neglects thickness. Because of this assumption, VLM method can't consider viscose effets. In vortex lattice model which is build up for propeller, blades are divided into panels along span and chord. Flow around propeller can be modelled by horse shoe vortices on every panel. Vortex strenghts are problem's unknown parameters. In this method, horse shoe vortex strengths on blades opposite to each other are same. So, calculating vortex strenghts on one of the blades is enough, it is not necessary to calculate for other blades. Solution of equation system can be done by using surface boundary condition. Surface boundary condition assumes that flow's component vertical to boundary is equal to zero. Aerodynamic forces can be calculated from vortex strenghts calculated by Joukowsky lifting law. Trailing vortices behind 3D blade's wake is propagating in blade's moving frame. In propellers, trailing vortices are propagating behind propeller by helical trajectories on wake at the effect of revolving and advancing movements of propeller. Root vortex is extending along rotor axis lineerly. Going far from propeller, wake vortex straitens and similarity in geometry breaks. In this case it is supposed to be damped in far regions from propeller. Helix wake behind a propeller generated by vortex sheet which is built up by trailing vortices can be modelled by several wake methods. Wake methods can be divided into two basic methods, steady and transient methods. Prescribed wake model is a steady method which describes wake as a function. Free wake model is a transient method which includes rolling on blade tip and calculates wake vortex coordinates iteratively. In this study, free wake model isn't preferred because it needs long computational time and there is a little role on accuracy of the results. Prescribed wake model is selected in order to obtain wake vortex coordinates. Total thrust and moment on a propeller can be calculated if propeller blade geometry, blade section characteristics, propeller advance velocity, rotation velocity and density of air is given. By using these values propeller performance parameters (thrust coefficient, moment coefficient and efficiency) can be calcuated. Program results are validated by using different blade profils. Propeller results are validated by using a test propeller with known geometry and performance. Prandtl lifting line and blade element-momentum based programs which are developed before are used for validation of VLM based study results. Validation studies for propeller blades are performed by using lifting line method program (MAYwingPLL). In first application, uniform sheet shaped blade and blade which has a imaginary profile without camber with a smaller lifting curve slope is analysed. According to results of a profile without camber, lifting coefficient values calculated in different angles of attack with lifting line method and vortex lattice method are similiar to each other. Numerically calculated values are compared with analytically calculated values taken from literature. Lifting coefficient values calculated analytically for sheet shaped blade and a blade which has a imaginary profile, are greater than numerically calculated values because lift in 3D blade is smaller than lift in 3D blade in same angle of attack. Because 3D blade has an induced drag component added to lift which doesn't exist in 2D blade. In second application, camber is added to imaginary profile in order to see effect of camber in numerical method results.Lifting coefficients that are calculated by lifting line method has changed under the effect of profile camber. Lifting coefficients that are calculated by lifting line method are greater than lifting coefficients that are calculated by vortex lattice method. In third application, real blade profile data which includes non-lineer region data of an airfoil is used. NACA 4412 is chosen as a blade section profile. In results which are calculated by two methods (lifting line theory and vortex lattice method), lifting coefficients in lineer region is increasing lineerly similiar to lifting coefficients in literature. There isn't any reduction in slope of lifting curves of two methods in blade profile's non-lineer region because viscose effects are not included. In last blade test study, numerical results are compared with experimental study results in order to examine success of the method. NACA 4415 profile is chosen as rectangular blade section profile. In this application, calculations are made for rectangular blade with aspect ratio values 6, 9 and 12. Results are given compared with 3D experimental results given by Ostowari and Naik. When aspect ratio value is small (AR=6), vortex lattice method results are converging to experimental results between angle of attack interval 0°
Benzer Tezler
- Pervanelerin girdap halkası yöntemi ile aerodinamik analizi
Aerodynamic analysis of propellers by vortex ring method
SADRETTİN KAYNAR
Yüksek Lisans
Türkçe
2017
Havacılık Mühendisliğiİstanbul Teknik ÜniversitesiUçak ve Uzay Mühendisliği Ana Bilim Dalı
PROF. DR. MAHMUT ADİL YÜKSELEN
- Düzensiz akımdaki nozullu pervanelerin dizaynı ve analizi
Design and analysis of ducted propellers in non-uniform velocity field
FAHRİ ÇELİK
Doktora
Türkçe
2005
Gemi MühendisliğiYıldız Teknik ÜniversitesiGemi İnşaat Mühendisliği Ana Bilim Dalı
PROF.DR. MESUT GÜNER
- Bir insansız hava aracı pervanesinin hesaplamalı akışkanlar dinamiği yardımıyla sayısal analizi ve aeroakustik optimizasyonu
Numerical analysis and aeroacoustic optimization of a uav propeller using computational fluid dynamics
VEYSEL EROL
Yüksek Lisans
Türkçe
2024
Havacılık ve Uzay Mühendisliğiİskenderun Teknik ÜniversitesiHavacılık Bilimi ve Teknolojileri Ana Bilim Dalı
DR. ÖĞR. ÜYESİ AHMET ŞUMNU
DR. ÖĞR. ÜYESİ YÜKSEL ERASLAN
- Alüminyum matriksli metal kompozit malzemelrin üretilmesi
The Manufacturing aluminyum metal matrix composite materials
İBRAHİM MUTLU