A general purpose large-eddy simulation/probability density function simulator on block structured grids

Abstract

Bu doktora tezinde, türbülanslı alevlerin sayısal modellenmesi için genel amaçlı büyük çalkantı simülasyonu (Large Eddy Simulation-LES)/olasılık yoğunluk fonksiyonu (Probability Density Function-PDF) simülatörü geliştirilmiştir. LES/PDF çözücüsü hibrit bir çözüm yöntemidir. Bu yöntemde, filtre edilmiş kütle ve momentum korunumu denklemlerinin çözümü için sonlu hacimler yöntemini (LES çözücüsü) kullanılmıştır. Türbülans-yanma etkileşimleri kimyasal bileşenlerin birleşik olasılık yoğunluk fonksiyonu yöntemiyle modellenmiş ve model denklemlerinin sayısal çözümü için Lagrangian Monte Carlo tabanlı bir partikül algoritması (PDF çözücüsü) kullanılmıştır. LES ve PDF çözücülerinin her ikisi de açık kaynaklı OpenFOAM yazılım paketi içerisinde geliştirilmiş ve pakete tümüyle entegre edilmiştir. Detaylı kimyasal yanma hesaplamalarının etkin bir şekilde yapılabilmesi için tablolama (the in situ adaptive tabulation (ISAT)) metodu yeni LES/PDF çözücüsü ile birleştirilmiştir. Bu çalışmada geliştirilen hibrit yöntem yapısal olmayan çözüm ağları esas alınarak geliştirilmiş ama bütün hesaplamalar blok-yapılı (block-structured) çözüm ağlarında yapılmıştır. Sayısal yöntem yapısal olmayan çözüm ağlarına kolaylıkla adapte edilebilir formdadır. Zhang ve Haworth (2004) tarafından geliştirilen üç aşamalı hız düzeltme algoritması mevcut LES/PDF yöntemine eklenmiş ve performansı gerçekçi türbülanslı alevlere uygulanarak gösterilmiştir. Yeni LES/PDF çözücüsü ayrıca geleneksel alt bölgeye ayırma (domain decomposition) yöntemi kullanılarak paralelleştirilmiştir. LES ve PDF çözücüleri arasındaki tutarlılık Wang ve Pope (2011) tarafından geliştirilen tek-yönlü birleştirme yöntemi kullanılarak test edilmiştir. Mevcut LES ve PDF çözücülerinin sayısal çözüm seviyesinde çok tutarlı sonuçlar verdiği gösterilmiştir. Bu LES ve PDF çözücülerinin doğru olarak birleştirildiklerine delalet etmektedir. Daha sonrasında ise, üç aşamalı hız düzeltme algoritmasının her bir aşaması, LES ve PDF çözücülerinin tutarlılıkları üzerindeki etkilerini belirlemek ve hesaplama maliyetini tespit etmek için incelenmiştir. Düzeltme metodu sayısal çözüm seviyesinde tutarlılığın sağlanması açısından çok faydalı bulunmuştur. İlk iki aşama basit türbülans alevleri için yeterli olmakla birlikte özellikle geri çevrimli ve karmaşık yapılı akış durumlarında üçüncü aşama düzeltmeye ihtiyaç duyulduğu tespit edilmiştir. Daha sonrasında ise, yeni hibrit LES/PDF çözücüsü bir detaylı kimyasal kinetik model kullanılarak pilotlu bir jet difüzyon türbülanslı alevine (Sandia Flame-D) uygulanmıştır. Metan/hava yanması indirgenmiş bir kimyasal mekanizma (ARM1) ile temsil edilmiştir. Sonuçların deneysel veriler ve daha önceki LES/PDF simülasyon sonuçları ile çok uyumlu oldukları gösterilmiştir. Son olarak, yeni LES/PDF çözücüsü Cambridge seyrekleştirilmiş alev serisindeki ön-karışımlı türbülanslı bir aleve uygulanmıştır. Sonuçların deneysel sonuçlar ile makul seviyede uyumlu oldukları görülmüştür.

A general-purpose large-eddy simulation (LES)/probability density function (PDF) methodology is developed for simulations of turbulent reacting flows. The LES/PDF solver is a hybrid solution methodology consisting of (i) a finite volume (FV) method for solving the filtered mass and momentum equations (LES solver), and (ii) the Lagrangian Monte Carlo based particle algorithm (PDF solver) for solving the modeled transport equation of the filtered joint PDF of compositions. Both the LES and PDF methods are developed and combined to form a hybrid LES/PDF simulator within the OpenFOAM framework. The in situ adaptive tabulation (ISAT) method developed by Pope /cite{Pope-1997} is incorporated into the new LES/PDF solver for efficient computations of combustion chemistry with detailed reaction kinetics. The method is designed to utilize a block structured mesh and can be readily extendible for the block unstructured grids. The three-stage velocity correction method of Zhang and Haworth /cite{ZhangHaworth} is also incorporated into the hybrid algorithm to interpolate the LES velocity field onto particle locations accurately and to enforce the consistency between the LES and PDF solvers at the numerical solution level. The hybrid algorithm is also fully parallelized using the conventional domain decomposition approach.First, the consistency of the FV-LES solver and the Lagrangian-PDF solver is examined by using the one-way coupling methodology developed by Wang and Pope /cite{Wang-2011}. It is found that the LES and the PDF solvers exhibit very good consistency at the numerical solution level demonstrating accurate coupling of the LES and PDF algorithms. Then the performance of the three-stage velocity correction algorithm is investigated to determine effects of each stage on the consistency of the LES and PDF solvers and the computational cost required by each stage. The correction algorithm is found to be very effective in enforcing the consistency conditions at the numerical solution level and first two stages are sufficient for simple reacting flows while the third stage is needed in the case of more complex flows involving recirculation regions. The predictive capability of the LES/PDF solver with the detailed chemistry representation is then examined by simulating a turbulent piloted methane/air jet diffusion flame (Sandia Flame-D). An augmented reduced mechanism (ARM1) /cite{ARM1} is used for the description of methane/air combustion. The results are found to be in very good agreement with the experimental measurements as well as with the earlier LES/PDF simulations. Finally, the new LES/PDF solver is applied to study a turbulent premixed flame from the Cambridge turbulent stratified flame series /cite{StratifiedFlame1}. The results are found to be in reasonably good agreement with the experimental data.