A robust consistent hybrid finite-volume/particle method for solving the pdf model equations of turbulent reactive flows
Reaksiyon içeren türbülanslı akışların modellenmesi için birleşik olasılık yoğunluk fonksiyonu (boyf) model denklemlerinin çözümünde kullanılan yeni bir gürbüz hibrit sonlu-hacimler/parçacık-tabanlı monte carlo metodu
- Tez No: 367649
- Danışmanlar: DOÇ. DR. METİN MURADOĞLU
- Tez Türü: Doktora
- Konular: Makine Mühendisliği, Mechanical Engineering
- Anahtar Kelimeler: Belirtilmemiş.
- Yıl: 2014
- Dil: İngilizce
- Üniversite: Koç Üniversitesi
- Enstitü: Fen Bilimleri Enstitüsü
- Ana Bilim Dalı: Makine Mühendisliği Ana Bilim Dalı
- Bilim Dalı: Belirtilmemiş.
- Sayfa Sayısı: Belirtilmemiş.
Özet
Reaksiyon içeren türbülanslı akışların modellenmesinde yaygın olarak kullanılan birleşik olasılık yoğunluk fonksiyonu (BOYF) model denklemlerinin istatistiki olarak zamandan bağımsız çözümü için yeni bir gürbüz hibrit sonlu-hacimler/parçacık-tabanlı Monte Carlo metodu geliştirilmiştir. Bu yöntem Muradoğlu vd. [40, 39] tarafından geliştirilen hibrit yöntemin bazı önemli yetersizliklerini gidermek için dizayn edilmiştir. Yapılan birçok alev simülasyonunda, orijinal hibrit algoritmasındaki yoğunluk-tabanlı sonlu-hacimler çözücüsünün aşırı disipatif olduğu ve yeterince gürbüz olmadığı gözlenmiştir. Bu yetersizlikleri gidermek amacıyla, açık kaynak kodlu, OpenFOAM, paketinde mevcut olan basınç-tabanlı PISO algoritması, Favre-ortalamalı kütle ve momentum denklemlerinin çözümü için kullanılmıştır. Buna paralel olarak çalkantı hızları-türbülans frekansı-kimyasal kompozisyon BOYF'nun taşınım denklemi iseparçacık-tabanlı Monte Carlo algoritması kullanılarak çözülmüştür.Bu sonlu-hacimler ve partikül yöntemleri daha sonra tutarlı hibrit metodu çerçevesinde birleştirilmiştir. Sonlu- hacimler ve partikül yöntemleriyle çözülen bütün denklemler doğrudan modellenmiş BOYF taşınım denkleminden elde edildiğinden, mevcut hibrit metot çözülen denklemler açısından tamamıyla tutarlıdır. Buna ilave olarak, eski hibrit metodundan farklı olarak, mevcut hibrit yöntemde ortalama yoğunluk parçacık kodu tarafından hesaplanmakta ve sonlu-hacim algoritmasına gönderilmektedir. Bu sayede orijinal hibrit metodunda gereğinden fazla olan yoğunluk alanı elimine edilmiş ve dolayısı ile yeni hibrit yönteminin sayısal çözüm seviyesinde çok daha tutarlı olmasını sağlamıştır. Ayrıca basınç tabanlı sonlu-hacimler yöntemi sayesinde aşırı sayısal difüzyonun oluşması engellenmiş ve her akış için kalibre edilmesi gereken hiçbir parametreye ihtiyaç kalmamıştır. Yeni hibrit yöntemi ilk olarak reaksiyon içermeyen (soğuk) ve reaksiyon içeren küt-cisim akışlarına uygulanmıştır. Yöntemin istatistiki olarak zamandan bağımsız çözüme başarı ile ulaştığı ve çözüm ağı yakınsaması gösterilmiştir. Orijinal hibrit yönteminin aksine, yeni hibrit algoritmasının çözüm ağının sıklaştırmasına karşı son derece gürbüz olduğu ve hiçbir şekilde fiziksel olmayan akım ayrılmasına yol açmadığı tespit edilmiştir. Ayrıca deterministtik bias hatasının hibrit yöntemde esas itibarı ile elimine edildiği ve dolayısı ile eski hibrit yönteminin en önemli avantajının muhafaza edildiği gösterilmiştir. Buna ek olarak mevcut yöntemle elde edilen sonuçların, daha önceki BOYF hesaplamalarıyla ve deneysel verilerle uyumlu olduğu görülmüştür. Daha sonra yeni hibrit yöntemi çok daha kompleks olan Sydney girdaplı küt-cisim alevi (SM1)'e başarıyla uygulanmıştır. Yöntemin son derece zor olan bu alev için de çok gürbüz olduğu ve hesaplama sonuçlarının deneysel verilerle uyumlu olduğu gösterilmiştir. Son olarak mevcut hibrit yöntemi ISAT algoritmasıyla birleştirilmiş ve bu sayede türbülanslı alevlerin simülasyonlarında detaylı kimyasal kinetik modellerinin kullanılabilmesine imkân sağlanmıştır. ISAT ile birleştirilen yeni hibrit yöntemi ARM2 detaylı kimyasal kinetik modeli kullanılarak girdapsız ve girdaplı küt-cisimle stabilize edilmiş türbülanslı alevlere başarı ile uygulanmış ve sayısal sonuçların daha önceki BOYF simülasyonları ve deneysel verilerle uyumlu olduğu gösterilmiştir. Simülasyonu yapılan bütün akışlarda, partikül kodunda hesaplanan ortalama yoğunluk alanının çok fazla istatistiki gürültü hatası içerse bile basınç-tabanlı sonlu hacimler yönteminde kararsızlığa yol açmadığı ve sonlu- hacimler yönteminin son derece gürbüz olduğu tespit edilmiştir.
Özet (Çeviri)
A new robust hybrid finite-volume (FV)/particle method is developed for solving joint probability density function (JPDF) model equations of statistically stationary turbulent reacting flows. The method is designed to remedy the deficiencies of the hybrid algorithm developed by Muradoglu et al. [J. Comp. Phys.,154, 342(1999); J. Comp. Phys.,172, 842(2001)]. The density-based FV solver in the original hybrid algorithm has been found to be excessively dissipative and yet not very robust. To remedy these deficiencies, a pressure-based PISO algorithm in the open source FV package, OpenFOAM, is used to solve the Favre-averaged mean mass and momentum equations while a particle-based Monte Carlo algorithm is employed to solve the fluctuating velocity-turbulence frequency-composition JPDF transport equation. The mean density is computed as a particle field and passed to the FV method. Thus the redundancy of the density fields in the original hybrid method is removed making the new hybrid algorithm more consistent at the numerical solu- tion level. The new hybrid algorithm is first applied to simulate non-swirling cold and reacting bluff-body flows. The convergence of the method is demonstrated. In contrast with the original hybrid method, the new hybrid algorithm is very robust with respect to grid refinement and achieves grid convergence without any unphysical vortex shedding in the cold bluff-body flow case. In addition, the results are found to be in good agreement with the earlier PDF calculations and also with the available experimental data. Then the new hybrid algorithm is successfully applied to simulate the more complicated Sydney swirling bluff-body flame 'SM1'. The method is also very robust for this difficult test case and the results are in good agreement with the available experimental data. In all the cases, the PISO-FV solver is found to be highly resilient to the noise in the mean density field extracted from the particles. Fi- nally the method is coupled with the ISAT algorithm [53] for efficient treatment of detailed chemistry. Extensive calculations have been performed for the swirling and non-swirling bluff-body flames using ISAT.
Benzer Tezler
- Bazı tekstil boyalarının penicillium funiculosum küfü ile biyosorpsiyon ve biyoakümülasyon özelliklerinin incelenmesi
Investigation of biosorption and bioaccumulation properties of some textile dyes with penicillium funiculosum
NURSENA DEMİR
Yüksek Lisans
Türkçe
2025
BiyokimyaSakarya ÜniversitesiKimya Ana Bilim Dalı
DOÇ. DR. SEMRA YILMAZER KESKİN
- Peridinamik yaklaşım ile eş yönlü fiber takviyeli karma malzemelerin mekanik davranışlarının incelenmesi
Investigation of the mechanical behavior of unidirectional fiber-reinforced composites by using peridynamic approach
ENES OCAK
Yüksek Lisans
Türkçe
2025
Savunma ve Savunma Teknolojileriİstanbul Teknik ÜniversitesiSavunma Teknolojileri Ana Bilim Dalı
PROF. DR. ŞAFAK YILMAZ
- Frekans seçici yüzeyler (FSY) ile duyarlılığı iyileştirilmiş mikroşerit yama antenli hibrit yakın alan glikoz tespit sensörü tasarımı
Design of a hybrid near-field glucose sensing sensor based on a microstrip patch antenna with sensitivity enhancement using frequency selective surfaces (FSS)
UMUT KÖSE
Doktora
Türkçe
2026
Elektrik ve Elektronik Mühendisliğiİstanbul Teknik ÜniversitesiElektronik ve Haberleşme Mühendisliği Ana Bilim Dalı
PROF. DR. MESUT KARTAL
- Karma sistemlerin tümleyen değişkenli modelleri
Complementarity modeling of hybrid system
SELİM TÜRKYILMAZ
Yüksek Lisans
Türkçe
1999
Elektrik ve Elektronik Mühendisliğiİstanbul Teknik ÜniversitesiElektronik ve Haberleşme Mühendisliği Ana Bilim Dalı
DOÇ. DR. KÜLMİZ ÇEVİK
- Dört rotorlu insansız hava araçları için fizik rehberli artık geri beslemeli sinir ağı ile çok ajanlı dağıtılmış lineer olmayan model öngörülü formasyon kontrolü
Multi-agent distributed nonlinear model predictive formation control for quadrotor UAVs using a physics-guided residual feedback neural network
DENİZ MUHACİR
Yüksek Lisans
Türkçe
2026
Elektrik ve Elektronik Mühendisliğiİstanbul Teknik ÜniversitesiKontrol ve Otomasyon Mühendisliği Ana Bilim Dalı
PROF. DR. YAPRAK YALÇIN