Search

In the present paper, two-dimensional unsteady flows over circular cross-section cylinders are analyzed numerically. The effects of placement of the cylinders are investigated for two different arrangements: tandem and side-by-side. Several turbulence models are tested, and it is found that Spalart-Allmaras turbulence model is the best among one- and two-equation turbulence models. The most appropriate time step, which is one of the important parameters in unsteady simulations, is found as 0.002 seconds. After successful validations, the cylinders are positioned as side-by-side and tandem. The effects of the arrangement on flow regime, drag coefficient, lift coefficient and Strouhal number are presented for various distances between the cylinders. It is found that the flow is almost steady without any vortex in the gap when cylinders are in tandem and the gap between them is low. In contrast, the interactions are strong in case of side-by-side arrangement at the lowest gap. When the gap increases, the flow is affected that results in change on the global parameters.

In the present three-dimensional study, laminar flow over a backward-facing step was investigated in the framework of the finite volume method for various Reynolds numbers ranging from 100 to 1000. Attempts were made to clarify the effects of the slope of the downstream heated bottom wall behind the step on the reattachment length and heat transfer. It was found that the length of the recirculation region changes slightly with the slope angles at lower Reynolds numbers while it changes dramatically as Reynolds number increases. The longest length of the recirculation region was obtained at the highest slope angle in the negative vertical-direction while the shortest one was determined at the maximum slope angle in the positive vertical direction. Increasing slope angles in both directions also provide a higher local Nusselt number in some distance behind the step.

Bilindiği üzere doğadaki hemen hemen tüm akışlar türbülanslıdır. Türbülansın görüldüğü akışlardan biri de çapraz akıştaki jet akışları ya da diğer adıyla çapraz jetlerdir. Çapraz akıştaki jet akışları jet akış yönünün içerisine püskürtüldüğü serbest akış yönünden tamamen farklı olduğu bir akış türünü ifade eder. Serbest akış, sahip olduğu momentumun etkisiyle lüleden püskürtülerek sağlanan jet akışını kendi doğrultusunda akmaya zorlayarak yönünü değiştirir. Çapraz akıştaki jetler, düşük hızlı baca dumanlarından, yüksek hızlara sahip jet motorlarının yanma odalarındaki akışa ve hipersonik roket ve füze kontrol sistemlerine, düşük kütle akışlı sınır tabaka kontrol sistemlerinden, dikey olarak inip havalanabilen V/STOL uçaklara kadar pek çok mühendislik olayında kendini göstermektedir.Bu çalışmada çapraz akıştaki jet akışları problemi deneysel ve hesaplamalı akışkanlar dinamiği (HAD) kullanılarak ele alınmıştır. HAD ile çapraz akıştaki jet akışlarının incelenmesi amacıyla, gerek endüstriyel uygulamalarda mühendislerce, gerekse de üniversite ve araştırma enstitülerinde bilim adamları ve araştırmacılar tarafından yaygın olarak kullanılan FLUENT 6.3.26 programından yararlanılmıştır. Sayısal kısımda lüle kesit geometrisi, lüle sayısı, hız oranlarının etkisinin yanında ayrıca lüle ile çapraz akış arasındaki eğim açısı da dikkate alınarak, bu gibi değişkenlerin akış yapısı üzerindeki etkileri kapsamlı olarak ele alınmıştır. Bunun için dairesel ve kare kesitli olmak üzere iki farklı lüle kesit geometrisi, bir lüle ve üç tane lüle yan yana olmak üzere iki farklı lüle sayısı, R=0.2, 0.5, 1.0 ve 2.0 olmak üzere dört farklı hız oranı ve a=30°, 45°, 60°, 75° ve 90° olmak üzere beş farklı eğim açısı incelenmiştir. Çapraz akıştaki baskın girdap yapılarından en önemli olan zıt-dönüşlü girdap çiftinin hız oranı, eğim açısı ve lüle kesit geometrisi ile nasıl değiştiği araştırılmıştır. Eğim açısı ve hız oranının artmasıyla zıt-dönüşlü girdap çiftinin ölçek olarak büyüdüğü ve yan yana üç tane lüle kullanılması durumunda oluşan zıt-dönüşlü girdap çiftinin, tek bir lüle kullanılmasına göre daha farklı bir yapıda olduğu görülmüştür. Bunun yanında özellikle gaz türbinlerinin film soğutulmasında büyük öneme sahip film soğutma verimi de en uygun lüle kesit geometrisi, eğim açısı ve hız oranı için araştırılarak sonuçlar irdelenmiştir. Film soğutma veriminin en yüksek değeri, eğim açısı 30° için elde edilmiş ve dairesel kesitli lüle kullanılmasının daha iyi olacağı sonucuna varılmıştır.Deneysel ölçümler 304 mm´304 mm´914 mm test kesidine sahip bir rüzgâr tünelinde değişik hız oranları için farklı yatay ve düşey mesafelerde gerçekleştirilmiştir. Lüleden sağlanan jet akışı ve rüzgâr tünelinde gerçekleştirilen serbest akış hızı sabit sıcaklık hızölçeri (sıcak tel hızölçeri) kullanılarak ölçülmüştür. Jet akışının sağlandığı lüle dairesel kesitli olup, ölçümler 90° eğim açısında R=0.5, 1.0 ve 1.5 için gerçekleştirilmiştir. Elde edilen deneysel veriler Fourier, dalgacık (wavelet) ve uygun ortogonal ayrıklaştırma (POD) metotları kullanılarak incelenmiş, akışa ait güç tayfları (power sprectrum), dalgacık haritaları ve frekans-zaman alanındaki değişimleri sunulmuştur. Fourier analizine göre, bir serbest jette lüle çıkışındaki güç tayfı ilerleyen diğer mesafelere göre laminar karakter göstermiş ve lüle çıkış hızı artırıldığında frekans enerjisi düşmüştür. Üç farklı hız oranı (R=0.5, 1.0 ve 1.5) için yapılan mukayesede daha yüksek hız oranında daha düşük frekans enerjisi saptanmıştır. Bu durum lülenin merkezinde jetin çapraz akış önünde bir engel gibi davranarak çapraz akış etkisini sönümlemesi ve ilerleyen mesafelerde etkisini kaybetmesine bağlanmıştır. Verilerin dalgacık (wavelet) analizinde en düşük hız oranı olan R=0.5'te baskın frekansın 45 Hz'te kendini gösterdiği, ancak R=1.0 olduğunda bu değerin 85 Hz ile 105 Hz arasında olduğu gözlenmiştir. Enerjinin modlara göre dağılımlarının gösterildiği uygun ortogonal ayrıklaştırma (POD) metodunda ise düşük hız oranındaki modların yüksek hız oranındaki modlara nazaran daha fazla enerji taşıdıkları görülmüştür. Lüle çıkışından uzaklaşıldıkça toplam enerjinin daha fazla sayıda moda dağıldığı gösterilmiştir.Akım görüntüleme tekniği kullanılarak çapraz akıştaki jet akışlarında görülen büyük ölçekli yapılar ve akış yapısı gösterilmiştir. Bunun için bir duman üreteci ve dijital video kayıt cihazı kullanılarak rüzgâr tüneli ve jet hızları değiştirilerek akış alanı incelenmiştir. Yüksek jet hızlarında lüleden püskürtülen jet, çapraz akış içerisine daha dik bir şekilde nüfuz ederek akmaya devam etmiştir. Lüle hızı kademeli olarak azaltıldığında ise rüzgâr tünelinden sağlanan çapraz akış daha baskın hale gelmiş ve jet akışını yenerek, kendi doğrultusunda akmaya zorlamıştır. Rüzgâr tüneli yüksek hızda çalıştırıldığında bu baskın özellik daha çok görülmüş ve jet akışı gittikçe düzleşerek test kesidinden akmaya başlamıştır. Rüzgâr tüneli ve lülenin düşük hızlarda aktığında saptanan kayma tabakası girdapları, özellikle çapraz akışın yüksek hızlarında tespit edilememiştir.Anahtar Kelimeler: Hesaplamalı akışkanlar dinamiği (HAD), türbülans modelleri, çapraz jet, sabit sıcaklık hızölçeri, dalgacık ayrıklaştırması, uygun ortogonal ayrıklaştırma, Fourier ayrıklaştırması, güç tayfı. As it is well known that, most of the flows in the nature are turbulent. One of the turbulent flows seen in the nature is jet in cross flow or transverse jet. Jet in cross flow expresses the flow at which its direction is completely different from the cross flow. Free stream forces the jet flow to change its direction due to its momentum. Jet in crossflow is encountered from the stack smoke to flow in combustion room and hypersonic misille and rocket control systems, low mass flux boundary layer control systems to V/STOL aircrafts.In this study, jet in cross flow is taken into account by using experimental and computational fluid dynamics (CFD). In computational fluid dynamics calculations a CFD code of FLUENT 6.3.26 which has been used commonly in industry by engineers and in institutions and universities by academicians and researchers are considered. In the numerical section, nozzle geometry, number of nozzles, velocity ratio and inclination angle between the crossflow and nozzle are investigated to explore the effect of these parameters on the flow structure. For this reason two different nozzle section geometries; two different nozzle arrangements; circular and square nozzles, one and three-side by side nozzles, four different velocity ratios; R=0.2, 0.5, 1.0 and 2.0, five different inclination angle; a=0°, 45°, 60°, 75° ve 90° are considered. It was realized that, the high velocity ratio and inclination angle between the channel and nozzle, the more penetration of jet and the enlarge the counter-rotating vortex pair (CRVP) in size. It was also reported that using three nozzles side by side lead to high CRVP in size than the use of one nozzle. Film cooling which is very important for cooling technique was also investigated for different jet inclination angle and blowing ratios. It can be concluded that, the high thermal film cooling efficiency was obtained for the lowest inclination angle of 30°.Experimental measurements are done in a wind tunnel which has 305 mm´305 mm´914 mm test section for different jet-to-crossflow velocity ratios at different stations. Jet in crossflow are measured by using constant temperature anemometer (CTA) which is traditional but most commonly used in the experiments. Nozzles used in the experiments are considered as circular and the inclination angle between the nozzle and crossflow are designed as 90°. All measurements are conducted for R=0.5, 1.0 and 1.5. The results are then decomposed by using Fourier, wavelet and proper orthogonal decomposition (POD) techniques. Besides power spectrums and histograms dealing with the flow, velocity distributions are showed. According to Fourier analysis, the power spectrum at the center of the nozzle are more similar to the laminar character that the streamwise stations. When jet exit velocity was increased, it was seen that the energy of the flow decreased ironically. Comparisons were made for three velocity ratios, (R=0.5, 1.0 and 1.5) and it was seen that, the highest velociy ratio has the lowest energy. At the center of the nozzle, jet acts as an obstacle to the coming crossflow and damps the effect of the crossflow. The most dominant frequency of the flow was detected as 45 Hz at the ratio of R=0.5 by using wavelet technique. At R=1.0, dominant frequencies were seen at 85 Hz to 105 Hz. Energy contribution of the flow was calculated by using proper orthogonal decomposition (POD) technique and it was seen that the lower velocity ratio leads to carry more energy of the flow than the higher velocity ratio.Flow visualization technique was also performed to reveal the flow structure. A fog generator and digital video camera were used for this purpose. The flow field was investigated for different jet to crossflow velocity ratio by changing the velocity of the jet and wind tunnel, respectively. At the highest velocity ratio, jet penetrate directly into the crossflow due to its high momentum and acts as obstacle on the coming crossflow. It was also showed that, when the jet velocity was increased, crossflow become more dominant and sweep the jet flow, bending it to the crossflow direction. Shear layer vortices were detected especially at lower velocity ratio while it could not be seen at high crossflow.Keywords: Computational fluid dynamics (CFD), turbulence models, transverse jet, constant temperature aneometer (CTA), wavelet transform, proper ortogonal decomposition (POD), Fourier decomposition, power sprectrum. 210

ÖZET Bu çalışmada, içerisinde düşey bir engel bulunan borudaki laminar, sürekli ve sıkıştırılamaz akım ve ısı transferi problemleri silindirik koordinatlarda iki boyutlu ve nümerik olarak incelenmiştir. Boru içerisindeki akışı modellemek için süreklilik, enerji, momentum denklemlerinin yatay ve düşey yöndeki bileşenleri ve girdap-akım fonksiyonu metodu kullanılmıştır. Sözü edilen denklemler sonlu faiklar metodu kullanılarak ayrıklaştırılıp eş zamanlı çözüme hazır hale getirildikten sonra, Fortran dilinde yazılmış bir bilgisayar programı ile çözdürülmüşlerdir. Çalışmada, engelin değişik geometrik oranlan ve boru girişine uzaklığının akım alanı ve ısı transferi üzerindeki etkileri gösterilmiştir. Borudan Prandtl sayıları sırasıyla 0.7, 7.0 ve 546 olan hava, su ve kullanılmamış motor yağı geçtiği kabulüyle, Reynolds sayısının 50 ve 500 değerleri için akım hatları ve ısı transferi durumu incelenmiştir. Sonuçlar, literatürden elde edilen çalışmalarla karşılaştırılmış ve benzer sonuçların elde edildiği görülmüştür. Anahtar Kelimeler: Laminar akış, silindirik koordinatlar, engelli boru, Navier-Stokes vııı ABSTRACT In this study, considering the flow is laminar, steady and incompressible, the flow and heat transfer through a pipe employing a vertical obstacle are investigated in cylindrical - coordinates with 2 dimensional numerical analysis. To form a model vertical and horizontal components of the continuity, energy, momentum equations and vorticity-stream function method are included. These equations to be solved simultaneously are discretized employing the finite difference method and they are solved using a computer that has been writtenin Fortran. Through this study, the effects of different ratios relating to the obstacle and the position of the obstacle from pipe entrance on the streamlines and heat transfer are investigated. Air, water and clean engine oil with Prantl numbers of 0.7, 7.0, 546 are let flow through the pipe and the streamlines and heat transfer for the Reynolds numbers of 50 and 500 are investigated. Finally, it is shown that the results obtained are consistent with the reported results in the literature. Keywords : Laminar flow, cylindrical coordinates, pipe with obstacle, Navier-Stokes. IX 46

In this study the experimental and computational results for a fast ship model is presented. The Reynolds Averaged Navier Stokes (RANS) equations and the nonlinear free surface boundary conditions are discretized by means of an overset grid finite volume scheme. The experiments are performed at Istanbul Technical University Towing Tank basin. In the numerical turbulent flow calculations, the relationship between the Boussinesq's hypothesis of turbulence viscosity and the velocities are obtained through the standard k-ε turbulence model. Simulations of turbulent free surface flows around the model are performed by using Star CCM+ solver and Volume of Fluid (VOF) model to capture the free surface between air and water. The total resistance of the ship model is compared with the experimental results. Bow and aft wave form developments are also investigated qualitatively. For Froude (Fn) numbers less than 0.25, the computations are found to be well satisfactory, giving efficient and accurate tool to predict curves of resistance. For relatively higher speeds (Fn>0.25) a low Reynolds number turbulence model may be more suitable to predict the resistance.

Bu çalışmada dikdörtgen kesitli bir düzlemde düzenli akım kullanılarak oluşturulmuş bir model içerisine, akıma dik olacak şekilde yerleştirilmiş iki farklı yüzey geometrisine sahip (yivli ve pürüzsüz) silindirler etrafındaki akışın sayısal olarak incelenmesi yapılmıştır. Bunun için sonlu hacimler metodunu kullanan ANSYS-FLUENT 19 paket programı kullanılmıştır. Çalışmada ilk olarak D=0.049m ölçülerinde pürüzsüz yüzeyli 2 silindir modellenerek düzenli akış içerisine tandem düzende yerleştirilmiştir. Re=65×103'de elde edilen değerlerin (Cd ve St) deneysel verilerle karşılaştırması yapılmıştır. Daha sonra model geliştirilerek yine aynı hesap alanında olacak şekilde yivli yüzeyli silindirler yan yana yerleştirilmiştir. Silindirler arasındaki uzaklık (L) değiştirilerek yeni analizler yapılmıştır. Yapılan yüzey geometrisi değişikliğinin ilk modelde elde edilen değerler üzerindeki etkileri incelenmiştir. Elde edilen yeni değerlerin deneysel verilere olan yakınlığı gösterilmiştir. In this study, the flow around the cylinders having two different surface (grooved and smooth) placed perpendicular to the flow into a model formed by using regular flow in a rectangular cross-sectional plane was investigated numerically. For this purpose, ANSYS-FLUENT 19 package program was used which uses finite volume method. In the study, firstly, two cylinders with smooth surface D = 0.049m were modeled and the values obtained in Re=65×103 in regular flow (Cd and St) were compared with experimental values. Then the model was developed and the rollers with rough surfaces were placed side by side in the same calculation area. New analyzes were performed by changing the distance (L) between the cylinders. The effects of surface geometry change on the values obtained in the first model were investigated. The proximity and usability of the obtained new values to the experimental values were determined. 66

Günümüzde mühendislerin karşılaştığı en zorlu problemlerden biri girdap kaynaklı titreşim problemidir. Girdap kaynaklı titreşimler, bir yapının maruz kaldığı akış sonucunda oluşan girdaplardan dolayı meydana gelir. Bu girdaplar periyodik hale gelerek yapının salınım yapmasına sebep olur. Bu olguya, girdap kaynaklı titreşim denmektedir. Bu salınım, yapıların bozulmasına neden ve ciddi hasarlar almasına sebep olabilir. Günümüzde bilgisayar teknolojilerindeki son gelişmeler sayesinde, girdap kaynaklı titreşimler, birçok araştırmacı tarafından incelenmeye başlanmıştır. Şimdiye kadar kapsamlı araştırmalar yapılmış olmasına rağmen, hala bu konuda geniş çapta araştırılacak yönler bulunmaktadır.Bu çalışmada, girdap kaynaklı titreşimlerin en önemli parametresinden birisi olan sönümleme parametresi incelenmiş ve sonuçlar önceki daha önceden yapılmış deneysel bir çalışma ile karşılaştırılmıştır. Simülasyonlar, tek serbestlik derecesinden akış yönüne dik hareket edebilen, silindirik bir yapıda gerçekleştirilmiştir. Yapılan simülasyonlarda; genel olarak, deneysel sonuçlarla aynı davranış yakalansada, simülasyonlar iki boyutlu alanda gerçekleştirildiğinden dolayı bazı sonuçlarda tutarsızlık gözlemlenmiştir. Birimsiz deplasman, frekans oranı, girdap modları ve faz açısı sonuçları gösterilmiş, genel olarak şu sonuca ulaşılmıştır; sönümleme artırıldığında genliğin azaldığı ve senkronize bölgeye geçişte bir ertelenme olduğu gözlemlenmiştir. Aynı sonuçlar karşılaştırılan deneysel çalışmada da gözlemlenmiştir. One of the difficult challenges that engineers encounter nowadays is the Vortex Induced Vibration (VIV). VIV occurs when a bluff body is subjected to external fluid flow. Vortices shed from the structures cause bluff bodies to oscillate. This oscillation can lead to failures of the structures. With the recent development of computer technologies, VIV have been caught attention by many researchers. Although, extensive research has been done until now, there are still research areas to be investigated on VIV. In this study, one of the most important parameters of VIV, which is damping of the structure, is investigated and the results are compared with the experimental study of the previous research. The simulations are done with 2D cylindrical cross-section structure that is free to move in transverse direction. Overall, the general trend is well captured, however due to 2D limitations some results are not consistent with experimental data. The results are shown and overall, it is concluded that the response of the amplitude decreases with the increased damping of the structure and a shift occurs in the beginning of the lock-in region with increasing damping. 88

Günümüzde sayısal yöntemler mühendislik problemlerinin çözümünde yaygın olarak kullanılırlar. Bunun nedeni, sayısal çalışmaların deneysel çalışmalara göre hem ekonomik hem de zaman açısından daha avantajlı olmasıdir. Çarpan jet akışları, son yıllarda üzerinde önemli ölçüde çalışmalar yapılan bir konudur. Yüksek ısı transfer oranına sahip olduklarından dolayı endüstride yaygın olarak kullanılırlar. Çarpan jetlerin ısı transfrer karakteristiklerine olan bu ilgi, çarpan jetlerin gaz türbinlerinde ve elektronik parçalarda ısınan parçaların soğutulmasında, tekstil ve kağıt ürünlerinin kurutulması gibi bir çok mühendislik uygulamalarında ki kullanımından kaynaklanır. Ayrıca, kütle transferinde de etkili olduklarından maden havalandırılmasında, tünel operasyonlarında, boya spreylerinde de tercih edilirler.Bu çalışmada, jetlerin farkli düzlemler üzerine çarptırılmasıyla oluşan akış yapıları ve bunların ısı transferi üzerindeki etkileri farklı Reynolds sayıları ve jet-çarpma düzlemi mesafeleri için incelenmiştir. Akışkan yoğunluğundaki değişimleri gözlemleyebilmek için üç farklı metod uyarlanmıştır. Hesaplama ortamı ICEM programında modellenmiş ve daha sonra Ansys Fluent programında analizler gerçekleştirilmiştir. Oluşturulan ağ örgüleri ve türbülans modelleri arasından en uygun olanı belirlenmiştir. Veriler yapılmış olan bir çalışmadan elde edilen sonuçlarla kıyaslanmıştır.Sonuç olarak çarpan jet akışlarında Re sayısıyla Nu sayısı arasında doğru bir orantı olduğu görülmüştür. Maksimum ısı transfer oranı H/D parametresinin 6 olduğu ve Re sayısının 20000 olduğu durumda elde edilmiştir. Sürtünme katsayısı ile Re sayısı arasında ters bir korelasyon olduğu gözlemlenmiştir. Nowadays numerical methods have been generally used to solve engineering problems. This is because numerical methods more economical than the experimental study and more beneficial in terms of time. In the last decades significant studies have been carried on the impinging jets. This impinging jets widely used in industry because of they have high heat transfer rate. To interest of studying this heat transfer mechanism through impinging jets is owing to its many engineering applications such as the drying of textile and paper products, cooling of electronic parts, cooling of gas turbine blades etc.In addition impinging jets are also effective on the mass transfer, they are use in main ventilation, tunnel operation and also in paint spray.In this study impinging jets have been investigated numerically on the different types of planes to observe heat transfer mechanism and flow structure around them for different Re numbers and several jet-target surface distances. In order to detect the alteration on fluid density, three different approaches have been adjusted. The experimental setup drawing and meshing in ICEM CFD programme and then analyzis made by FLUENT. In generated meshs and turbulence models the best have been selected. Parameters according to the plane made by comparing the results obtained from the study will be proposed.Consequently, in impinging jets, the correlation between the Re number and Nu number has been detected. Maximum heat transfer rate has been obtained when Re number equals to 20000 and H/D commensurates to 6. The correlation has been observed between friction coefficient and Re number. 117