@article { author = {Farahani, Majid and Jalali, Ali and Saghafi, Mohammad Ali}, title = {A second order closure model for calculating the turbulence fluxes in atmospheric boundary layer}, journal = {Journal of the Earth and Space Physics}, volume = {40}, number = {1}, pages = {107-122}, year = {2014}, publisher = {Institute of Geophysics, University of Tehran}, issn = {2538-371X}, eissn = {2538-3906}, doi = {10.22059/jesphys.2014.36700}, abstract = {Turbulent fluxes in the atmospheric boundary layer are calculated using different models. One of the most popular and operational models in this groundwork is the second order turbulent model of Mellor and Yamada that has been widely used to simulate the planetary boundary layer (PBL). Despite its popularity, it has been shown that this model has several weeknesses that affect its accuracy. The main deficiency of this model is its inability in distinguishing between vertical and horizontal components of turbulent kinetic energy. Also, this model can not accurately predict the height of the boundary layer. This deficiency in turn can lead to errors in the parameterization of physical processes. Canuto and his colleagues attempt to overcome these deficiencies by introducing a modified model that is an improved version of Mellor and Yamada model. This later model in theory has proven to have more reliable results in computing turbulent fluxes. To verify the results of this theoretical model, the model was coded and verified using some of the available data sets. The data sets used in this study include, the Sodar data of Institute of Geophysics of University of Tehran and Mehrabad airport sounding data. Also the standard Large Edy Simulation data set (LES) introduced in GABLS (Cuxart, 2006) project was examined. For better certainty of the new model performance, the Nakanishi's LES data that used in his study was used and applied to the model. The data used here present a boundary layer of the atmosphere in neutral or Ekman layer condition. Potential temperature does not exist in the Sodar data set and therefore it was computed using measured laps rate in the same data set. The output of model includes the main flux terms, ,, . Results show that the model produces acceptable values for fluxes and have almost the same vertical profile and trend. There are some differences in the predicted values of parameters. It seems that this model could recognize eddies of boundary layer with more accuracy. The differences between the vertical and horizontal TKE distribution is clear. Some of the differences are explained. The results of model as long as the boundary layer conditions match to those that theory proposed or assumed, are reasonable. The analytical data feed to the model results in better output. Compare to a well-known intercomparison program (GABLS project) the results are quite good. The PBL height is defined as the height at which the turbulent kinetic energy or the magnitude of the momentum flux decreases to a small fraction of the corresponding surface value. Results of experiments show that boundary layer height is estimated well in this model. The critical Richardson number predicted by the model is around unity rather than 0.2 as given by previous models. The larger critical Richardson number is in agreement with LES data. Experiments revealed that the presented model is not suitable for boundary layers in which the Richardson numbers within them approach 1 or are larger than 1.    }, keywords = {Turbulent fluxes,Turbulance in boundary layer,Closure problem,Mellor-Yamamda model}, title_fa = {فرایابی شارهای تلاطمی در لایه‌ مرزی جوّ با استفاده ازیک مدل بستار مرتبه دوم}, abstract_fa = {در این پژوهش، لایه‌ مرزی جوّ بر اساس یک مدل بستار مرتبه دوم که چنگ و همکاران (2002) معرفی کرده‌اند، شبیه‌سازی عددی می‌‌شود و مورد نقد و بررسی قرار می­گیرد. با استفاده از مدل یاد­شده، شارهای تلاطمی در لایه‌ مرزی  فرایابی می‌‌شوند و با منظور کردن آنها در معادلات دینامیکی وابسته به زمان لایه‌ مرزی ، متغیرهای تلاطمی محاسبه می‌شوند. در سیر تحول این مدل ، نظریه مطرح­شده در مدل ملور- یامادا، به­کار گرفته شده و جملات تلاطمی در جوّی نزدیک به شرایط ترمودینامیکی خنثی با استفاده از آن محاسبه شده است. همچنین نتایج حاصل از آن با داده­های پیچک­های بزرگ­مقیاسِ شبیه‌سازی شده مستند می‌شود و به‌صورت فرایابی مورد ارزیابی و مقایسه قرار می‌گیرد. بررسی صورت گرفته در مورد نتایج به‌دست آمده از اجرای این مدل، نشان‌دهنده کاهش و رفع برخی ضعف­های مدل ملور– یامادا در مدل جدید است. در مدل جدید، جمله ناشی از باد افقی در انرژی جنبشی پیچکی، چندین برابر جمله قائم آن محاسبه شده است. نتایج شارهای فرایابی شده،  نشان‌دهنده توانایی مدل در شناسایی پریشیدگی‌های موجود در میدان­های باد و دمای پتانسیل و مطابقت آنها با شارهای خروجی است.    }, keywords_fa = {Turbulent fluxes,Turbulance in boundary layer,Closure problem,Mellor-Yamamda model}, url = {https://jesphys.ut.ac.ir/article_36700.html}, eprint = {https://jesphys.ut.ac.ir/article_36700_89e0b3178e57d9823a0da5f5135e277b.pdf} }