@article { author = {Ashrafi, Khosro and Ahmadi Orkomi, Ali and Shafiepour Motlagh, Majid}, title = {Plume rise simulation via hybrid RANS-LES CFD method in neutral atmospheric conditions}, journal = {Iranian Journal of Geophysics}, volume = {11}, number = {1}, pages = {66-83}, year = {2017}, publisher = {Iranian Geophysical Society}, issn = {2008-0336}, eissn = {2783-168X}, doi = {}, abstract = {Most physical systems in nature have nonlinear dynamics and the system linearization of these physical systems is only a simplified assumption. The atmospheric motions, with the same philosophy, have a turbulent structure and ommiting the turbulent motions in the atmosphere is only a problem simplification assumption. When a buoyant jet of a chimney enters the atmosphere, it behaves like a turbulent flow, in which the atmospheric turbulence and self-generated turbulence of the plume play major parts. The present survey aimed at demonstrating the effects of turbulence on plume dynamics through computer simulation. The well-known turbulent flow simulation methods commonly used to simulate plume dynamics and atmospheric processes are: Direct numerical simulation (DNS), Reynolds averaged Navier-Stokes (RANS) and large eddy simulation (LES), which most distinguishing feature is their way of parameterizing the turbulence. As far as computational requirements, accuracy and turbulence simulation, the former two models are the two extremes while LES occupies an intermediate position between them, directly simulating the large-scale eddies and parameterizing the less important sub-grid scale (SGS) dissipative processes using sub-grid models (SGM). Most often, LES can predict the unsteadiness and intermittency of the turbulence structure, which is the most important feature of a buoyancy-driven jet. It should be noted that it is not efficient to employ full LES method when tackling an issue with certain unimportant zones. Furthermore, in the case of strong turbulent motions, the scale of flow structures near the rigid bodies are small and LES method method requires very fine grids that can increase its computational cost as large as DNS. To surmount this drawback, a hybrid RANS-LES method with a new mixed scale sub-grid parameterization model was applied to simulate the turbulent plume dynamics in ANSYS Fluent 14.5 software. The effectiveness and the accuracy of the mentioned turbulence simulation method was demonstrated through simulation study and experimental data in the neutral atmospheric conditions. Comparing the simulation results of the RANS method, the default hybrid RANS-LES method with static sub-grid scale parameterization and the new RANS-LES method with dynamic mixed scale parameterization indicated that the mean temperature profile at stack downstream was more accurately predicted by the new hybrid method. The root mean square error of plume rise estimation of the RANS method, the default RANS-LES method and the new hybrid RANS-LES method were 0.0437, 0.054 and 0.0323, respectively. It was further demonstrated that the Briggs integral plume rise model could not properly predict the plume rise in the presence of turbulence, because it does not consider the updraft and downdraft turbulence-induced motions in the atmosphere. Ultimately, we checked the capability of the new hybrid method to resolve the substantial parts of the turbulent motions. The turbulent energy transfer from the energy containing scales to inertial sub-range followed the well-known  law. The capability of the new hybrid method in predicting the mean profile and the turbulent structure can be employed in the study of the effects of turbulent parameters on plume rise in different atmospheric stability classes.}, keywords = {Plume,numerical simulation,Turbulence,hybrid RANS-LES method}, title_fa = {شبیه‌سازی صعود پَره دودکش با استفاده از روش ترکیبی دینامیک شاره‌های محاسباتی RANS-LES در شرایط خنثای جوّی}, abstract_fa = {اغلب سامانه­ها در طبیعت دارای دینامیک غیرخطی هستند و خطی­سازی آنها تنها یک فرض ساده­کننده است. تلاطم در جریان­های جوّی نیز اغلب حاکم است و آرام بودن یا آرام فرض کردن آنها، به‌ عنوان ساده­سازی مسأله تلقی می­شود. در صعود و پراکنش پَره دودکش Plume) (Stack، بویژه در فواصل نزدیک به آن، تلاطم جوّی و تلاطم ناشی از پَره دود خروجی از دودکش تأثیر قابل توجهی دارند. هدف اصلی این پژوهش بررسی تأثیر تلاطم‌ جوّی بر پره دودکش است. رسیدن به این مهم، نیازمند پیش‌بینی مناسب نقش تلاطم در رفتار پَره دودکش و صرف هزینه محاسباتی کمتر در مقایسه با دقت شبیه‌سازی است؛ بنابراین در این پژوهش برای نخستین بار از ترکیب روش‌های شبیه‌سازی پیچک­های بزرگ و میانگین­گیری رینولدز به منظور شبیه‌سازی رفتار پَره دودکش و یک روش ترکیبی دینامیک برای پارامترسازی نقش پیچک­های ریز استفاده شده است. شبیه‌سازی عددی رفتار پره دود با استفاده از روش مذکور و متداول انجام گرفت. توزیع دما در فواصل مختلف پایین­دست دودکش با داده‌های تجربی مقایسه و صحت نتایج بررسی شد. نتایج نشان داد که روش ترکیبی پیشنهادی نسبت به روش‌های میانگین­گیری رینولدز و روش ترکیبی موجود در نرم­افزار فلوئنت، توزیع دما در پایین­دست را با دقت بیشتری پیش‌بینی می‌کند. همچنین خطای تخمین صعود پَره دود محاسبه شده با روش‌های میانگین­گیری رینولدز، روش ترکیبی موجود در نرم­افزار فلوئنت و روش ترکیبی پیشنهادی در حالت خنثای جوّی به ترتیب برابر با 0437/0، 054/0 و 0323/0 است. با مقایسه صعود پَره دودکش با معادله انتگرالی بریگز، مشاهده شد که روش انتگرالی به دلیل عدم در نظر گرفتن اختلاط قائم ناشی از تلاطم، میزان صعود را بیش از مقدار واقعی تخمین می‌زند. صحت شبیه‌سازی تلاطم با استفاده از نمودار چگالی طیفی انرژی جنبشی متلاطم و پارامتر نسبت تلاطم شبیه‌سازی شده به تلاطم مدلسازی شده، بررسی شد. از نتایج این پژوهش می‌توان در بهینه‌سازی و اصلاح روابط موجود برای تخمین صعود پَره دودکش در شرایط مختلف پایداری، با در نظر گرفتن نقش تلاطم مکانیکی و گرمایی جوّ و تلاطم پَره دودکش استفاده کرد.}, keywords_fa = {پَره دودکش,شبیه‌سازی عددی,تلاطم,روش ترکیبی LES - RANS}, url = {https://www.ijgeophysics.ir/article_46716.html}, eprint = {https://www.ijgeophysics.ir/article_46716_15c83982f0253ebfa4f2c73db169070c.pdf} }