research centers


Search results: Found 29

Listing 1 - 10 of 29 << page
of 3
>>
Sort by

Article
Free Convection in an Inclined Concentric Annular Square Cavities Filled With Porous Medium and Heated By Non-Uniform
انتقال الحرارة بالحمل الحر في تجويف مائل بين فجوتين مربعتين متحدتين المراكز ممتلئ بوسط مسامي و مسخن بدرجة حرارة غير منتظمة

Author: Nabeel Mohammed Jasim نبيل محمد جاسم
Journal: Al-Khwarizmi Engineering Journal مجلة الخوارزمي الهندسية ISSN: 18181171 23120789 Year: 2010 Volume: 6 Issue: 3 Pages: 45-60
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

A numerical study of the two-dimensional steady free convection flow in an inclined annulus between two concentric square cavities filled with a porous medium is presented in this paper for the case when the side outer walls are kept with differentially heated temperature while the horizontal outer walls and the inner walls are insulated. The heated wall is assumed to have spatial sinusoidal temperature variation about a constant mean value. The Darcy model is used and the fluid is assumed to be a standard Boussinesq fluid. For the Cartesian coordinate system, the governing equations which were used in stream function form are discretized by using the finite difference method with successive under – relaxation method (SUR) and are solved by Gauss-Siedel iterative method. The upwind scheme was used for the transport terms in the energy conservation equation. The results are presented to demonstrate the streamlines, the isotherms, and the Nusselt number depending on the Rayleigh number ranging from (Ra =10 to 1000), dimension ratio from (Dr = 0.15 to 0.45), and the inclination angle from ( = 0o to 45o). Also the effects of the amplitude ( =0 to 1) and the wave number (f =0 to 5) of the heated side wall temperature variation on the free convection are investigated. The results show the effect of previous parameters (Ra, Dr, , , and f) on the flow fields and temperature profiles. It also show that the average Nusselt number is a strong function of the Rayleigh number, inclination angle, dimension ratio, and temperature variation. The peak value of the average Nusselt number based on the hot wall temperature is observed to occur at dimension ratio of (0.15), inclination angle of (40.1°), amplitude and wave number (1 & 0.75) for Rayleigh number of (1000).

نعرض في هذا البحث دراسة عددية لتدفق الحمل الحر المستقر في بعدين داخل تجويف مائل بين فجوتين مربعتين متحدتين المراكز ممتلئ بوسط مسامي حيث تكون الجدران الخارجية الجانبية مسخنة بدرجة حرارة مختلفة بينما تكون جميع الجدران الأخرى معزولة. الجدار المسخن له درجة حرارة تتغير وفقا لدالة جيبيه حول قيمة متوسط درجة الحرارة. تم الاعتماد على موديل دارسي في التمثيل الرياضي لانتقال الحرارة في الوسط المسامي وقد افترض إن المائع هو مائع بوسين المعياري. الحل العددي للمعادلات أنجز بواسطة حل معادلات كل من دالة الجريان و الطاقة باستخدام طريقة الفر وقات المحددة و تطبيق طريقة (كاوس – سيدل) مع استخدام طريقة تحت التراخي. تم تمثيل النتائج العددية بدلالة خطوط الانسياب و درجات الحرارة بالاعتماد على قيم مختلفة لعدد رايلي (من 10 إلى 1000) , نسبة البعد (من 0.15 إلى 0.45), و زاوية الميل (من o0 إلى o45). كما تم دراسة تأثير كل من سعة الموجة (من 0 إلى 1) و عدد الموجات (من 0 إلى 5) لدرجة الحرارة المتغيرة للجدار الجانبي الخارجي المسخن على انتقال الحرارة. النتائج بينت تأثير المعاملات السابقة, and f) ,φ (Ra, Dr, على خطوط الجريان و درجة الحرارة. كما بينت النتائج إن عدد نسلت هو دالة قوية من عدد رايلي , زاوية الميل, نسبة البعد, و التغير في درجة الحرارة. وتتمركز القيمة القصوى لمعدل عدد نسلت عندما تكون قيمة نسبة البعد مساوية إلى (0.15), و زاوية الميل (40.1o) و عندما تكون سعة و عدد الموجات مساوية إلى (1 و 0.75) لقيمة عدد رايلي (1000).


Article
AIDING AND OPPOSING MIXED CONVECTION FROM AN ARRAY OF CIRCULAR CYLINDERS IN A SATURATED POROUS MEDIUM
الحمل المختلط المساعد والمعاكس من صف اسطوانات دائرية في وســط مسامي مشبع

Author: Mohammed A. Ahmed محمد عبد احمد
Journal: Anbar Journal of Engineering Sciences مجلة الأنبار للعلوم الهندسية ISSN: 19979428 Year: 2008 Volume: 1 Issue: 2 Pages: 15-26
Publisher: University of Anbar جامعة الانبار

Loading...
Loading...
Abstract

Numerical solutions are presented for mixed convection from an array of circular cylinders embed in a saturated porous medium. The cylinders are at constant temperature(isothermal) and arranged in a staggered tube bank. Both aiding and opposing flow conditions are considered. Numerical calculations using finite difference method with body-fitted coordinates have covered a wide range of governing parameters(i.e.,10 ،ـ Re ،ـ 100, 0 ،ـ Gr ،ـ 400 and Pr = 0.7). Results are presented for streamline, isotherms and the local and the average Nusselt number at different values of the governing parameters. The present results are compared with previous theoretical results and show good agreement

أجريت دراسة عددية للحمل المختلط من صف اسطوانات دائرية موضوعة في وسط مسامي مشبع. الاسطوانات عند درجة حرارة ثابتة وحزمة الأنابيب مرتبة بشكل متخالف. تم دراسة كل من شرط الجريان المساعد والمعاكس.استخدمت الحلول العددية طريقة الفروق المحددة مع نظام مطابقة إحداثياتالجسم وغطت مدى واسع من المتغيرات الحاكمة (≤100 عدد رينولدز ≤ 10 , 400 ≤ عدد كرشوف ≤ 0 و عدد براندتل= 0.7 ). مثلت النتائج لخطوط الانسياب, درجة الحرارة وعدد نسلت الموضعي و المعدل عند قيم مختلفة للمتغيرات الحاكمة. قورنت النتائج الحالية مع نتائج عددية متوفرة في البحوث السابقة و أوضحت توافق جيد


Article
On the Effect of Flow Direction on Mixed Convection from a Horizontal Cylinder in a Saturated Porous Medium
تأثير اتجــاه الجريان على الحمل المختلط من اسطوانة أفقيــة في وســط مسامي مشبع

Author: Mohammed A. Ahmed محمد عبد احمد
Journal: Anbar Journal of Engineering Sciences مجلة الأنبار للعلوم الهندسية ISSN: 19979428 Year: 2007 Volume: 1 Issue: 1 Pages: 77-87
Publisher: University of Anbar جامعة الانبار

Loading...
Loading...
Abstract

A numerical investigation of mixed convection from a horizontal cylinder in a saturated porous medium is presented. The governing equations based on Darcy’s law are expressed in a body- fitted coordinate system and solved numerically by explicit method. The direction of the flow varies between the vertically up ward(assisting flow) and vertically downward(opposing flow). Results are presented for Reynolds number Re from 10 to 100 with Grashof numbers up to Gr =5Re. The Prandtl number was kept at a constant value of 0.7. results are presented for the streamlines and isotherms as well as the local and average Nusselt number at different values of governing parameters. Comparison with previous theoretical results show good agreement.

أجريت دراسة عددية للحمل المختلط من اسطوانة أفقية في وسط مسامي مشبع. المعادلات الحاكمة المعتمدة على قانون دارسي تم التعبير عنها بنظام مطابقة إحداثيات الجسم وحلت عدديا بالطريقة البينية. اتجاه الجريان يتغير بين شاقوليا باتجاه الأعلى إلى شاقوليا باتجاه الأسفل. النتائج مثلت لعدد رينولدز من 10 إلى 100 مع عدد كر شوف إلى (5 *عدد رينولدز).مع بقاء قيمة عدد براندتل ثابتة عند القيمة 0.7 . مثلت النتائج لخطوط الانسياب , درجة الحرارة و عدد نسلت الموضعي والمعدل عند قيم مختلفة للمتغيرات الحاكمة. المقارنة مع نتائج نظرية سابقة توضح توافق جيد.


Article
Numerical Investigation of Laminar Natural Convection in Rectangular Enclosures of Porous Media

Authors: Selah M. Salih --- Majid H. Majeed --- Ala'a A. Mahdi
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2007 Volume: 13 Issue: 3 Pages: 1666-1677
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

In this investigation, steady tow-dimensional natural convective heat transfer in rectangular porous cavity, (heated from below) with horizontal walls heated to uniform but different temperatures and adiabatic sides has been studied numerically. The numerical results of heat transfer rates are presented for porous Rayleigh numbers (Ra*), based on width of cavity, in the range (Ra*≤500), with layer aspect ratios (Ar), (height / width) ranging between (0.5≤Ar≤5). Plots of streamlines and isotherms to show the behavior of the flow and temperature distribution are presented. The current study shows that the Nusselt number is a strong function of the porous Rayleigh number, and the geometry of the cavity is represented by aspect ratios. Porous Rayleigh number has a large effect on the flow field, whereas any increase in (Ra*) results in changing the flow pattern from unicellular to multicellular flow. Correlation equation has been obtained to show the dependence of Nusselt number on the porous Rayleigh number, and aspect ratio (Ar), as this correlation will be beneficial in design of systems of thermal insulators in the energy storage engineering applications.

في هذا البحث، تم اجراء دراسة عددية لانتقال الحرارة بالحمل الطبيعي المستقر ثنائي البعد في تجويف مسامي مستطيل الشكل، (مسخن من الأأسفل) له جدران أفقية مسخنة الى درجات حرارة منتظمة ولكن مختلفة والجوانب معزولة,. تم تمثيل النتائج العددية لمعدلات انتقال الحرارة بعدد ر الي المسامي (Ra*) المحسوب على أساس عرض التجويف، في المدى (Ra*≤500) ولنسب ابعاد (Ar)، (الارتفاع / العرض) تمتد من(0.5≤Ar≤5) لقد تم عرض رسومات لخطوط الانسياب ولخطوط درجات الحرارة لوصف سلوك الجريان وتوزيع درجات الحرارة. بالاضافة، بينت ان عدد نسلت هو دالة قوية من عدد ر الي المسامي ونسبة الابعاد. ان عدد ر الي المسامي يمتلك تأثير كبير عللى مجال الجريان، حيث ان اي زيادة بعدد رالي (Ra*)المسامي يسبب تغيير لشك الجريان من جريان احادي الخلايا الى جريان متعدد الخلايا. تم ايجاد معادلة تقريبية تبين اعتمادية عدد نسلت على عدد ر الي المسامي ونسبة الابعاد، وان هذه العلاقة يمكن الاستفادى منها في تصميم العوازل الحرارية من اجل حفظ الطاقة في التطبيقات الهندسية.


Article
Numerical Study of Natural Convection Heat Transfer of a Square Eccentric Body Buried in a Porous Media
دراسة عددية لانتقال الحرارة بالحمل الطبيعي لجسم مربع مغمور لا مركزيا في وسط مسامي

Authors: Ramiz Ibraheem saeed رامز إبراهيم سعيد --- Ra'ed Ahmed Ali رائد احمد علي
Journal: AL Rafdain Engineering Journal مجلة هندسة الرافدين ISSN: 18130526 Year: 2009 Volume: 17 Issue: 4 Pages: 51-62
Publisher: Mosul University جامعة الموصل

Loading...
Loading...
Abstract

AbstractNatural convection in a saturated porous medium bounded by two horizontal eccentric ducts is studied numerically by solving the governing Darcy-Boussinsq equations using finite differences on very fine grid. This paper reports the effect of the position of a cooled square body with respect to the centre of the square cavity vertically and horizontally. The numerical results of heat transfer are presented for modified Rayleigh number Ra* (based on width of cavity) in the range (50 to 400) and for different size of the body with respect to the width of the cavity (0.25x0.25,0.3x0.3,0.5x0.5) and six positions of the body respect to the center of cavity. The results are represented by streamlines and isotherms to show the behavior of the flow and temperatures distribution. Also some graphics had been obtained which represent the relationship between Nusselt and Rayleigh numbers with eccentricity. (Key words: natural convection, porous medium, cavity)Key Words : Porous medium , Thermal nonequilibrium , Heat transfer , Natural convection.

تم إجراء دراسة عددية للحمل الطبيعي لوسط مسامي مشبع محصور بين مجريين مربعين غير متمركزين وذلك بواسطة حل المعادلات المتحكمة لدارسي – بوسنسك باستخدام طريقة الفروق المحددة على شبكة ناعمة جدا. في هذا البحث تم تناول تأثير موقع الجسم المربع البارد نسبة إلى مركز التجويف عموديا وأفقيا. تمثلت النتائج العددية لانتقال الحرارة بعدد رالي المطور (نسبة إلى عرض التجويف ) في المجال (50-400) وإحجام متعددة للجسم ( (0.25x0.25,0.3x0.3,0.5x0.5 ولستة مواقع للجسم نسبة الى مركز التجويف. تمثلت النتائج بخطوط الانسياب وخطوط درجات الحرارة التي تظهر صيغة الانسياب وتوزيع درجات الحرارة. كذلك تم الحصول على بعض الرسومات البيانية التي مثلت العلاقة بين عدد نسلت وعدد رالي مع اللا تمركزية.


Article
Numerical Study Of Natural Convective Heat Transfer In An Inclined Square Porous Layer
دراسة عددية لانتقال الحرارة بالحمل الطبيعي في طبقة مسامية مربعة مائلة

Authors: Muyassar E. Ismaeel ميسر ادريس اسماعيل --- Amir S. Dawood أمير سلطان داؤد
Journal: AL Rafdain Engineering Journal مجلة هندسة الرافدين ISSN: 18130526 Year: 2008 Volume: 16 Issue: 3 Pages: 81-95
Publisher: Mosul University جامعة الموصل

Loading...
Loading...
Abstract

AbstractIn this paper, steady two-dimensional natural convective heat transfer in an inclined square porous cavity with two parallel walls kept at constant different temperatures, while the other two parallel walls were well insulated, has been studied numerically. The governing equations have been solved using finite difference method. Results have been obtained for modified Rayleigh numbers between 0 and 300 and inclination angle between 0° (heated from below) and 90° (heated from side). The rate of heat transfer was found a strong function of modified Rayleigh number and inclination angle. The maximum heat transfer was occurred at about ( 50  ≤  ≤ 57  ) degrees of inclination. Keywords:Natural Convective Heat Transfer, Porous Medium, Inclined Enclosure.

الخلاصةفي هذا البحث ، تم إجراء دراسة عددية لانتقال الحرارة بالحمل الطبيعي المستقر ثنائي البعد في تجويف مسامي مربع الشكل له جداران متوازيان ذو درجات حرارة ثابتة ولكن مختلفة والجداران المتوازيان الآخران معزولان.حصلت الحلول للمعادلات الحاكمة بطريقة الفروق المحددة. ان معدلات انتقال الحرارة المحسوبة قدمت لعدد رالي المطور Ra* في المدى ( 0 ) الى (300 ) ولزاوية ميلان  تتغير من (0 ) (مسخن من الاسفل) الى (90) ( مسخن من الاعلى). ان الدراسة الحالية بينت ان معدل انتقال الحرارة هو دالة قوية من عدد رالي المطور و زاوية الميلان. كذلك بينت ان اقصى انتقال للحرارة يحدث عند زوايا ميلان حوالي (50 ≤  ≤ 57 ).


Article
The effect of slip condition on heat transfer of MHD oscillatory third order fluid flow in a channel filled with porous medium
التأثير المشترك لحقل مغناطيسي وانتقال الحرارة بالاشعاع في الجريان اللامستقر لمائع من الرتبة الثالثة في قناة مملوءة بوسط مثقب

Authors: Hala KahtanHamdi هالة قحطان حمدي --- Ahmed M. Abdulhadi احمد مولود عبد الهادي
Journal: Iraqi Journal of Science المجلة العراقية للعلوم ISSN: 00672904/23121637 Year: 2012 Volume: 53 Issue: 4 Pages: 887-893
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

The combined effect of a transverse magnetic field and radiative heat transfer to unsteady flow of a conducting optically third order fluid through a channel filled with saturated porous medium and non-uniform walls temperature is investigated.It is assumed that the no-slip condition between the wall and the fluid remains no longer valid. The third order fluid equations of continuity momentum and energy are obtained. Analytical solutions for problem are established.The effect of wall slip on velocity field is presented by figures. The basic properties of the flow are studied.

في هذا البحث تمت دراسة التأثير المشترك لحقل مغناطيسي وانتقال الحرارة بالاشعاع في الجريان اللامستقر لمائع من الرتبة الثالثة في قناة مملوءة بوسط بوريس وذات جدران مختلفة الحرارة وقد اعتبر أن شرط عدم الانزلاق بين الجدران والمائع غير متحقق .أوجدنا معادلات الاستمرارية والحركة والطاقة للمائع من الرتبة الثالثة وقد حصلنا على حلول تحليلية للمعادلات أعلاه وكذلك وضحنا الصفات الرئيسية للمسألة من خلال الرسوم .


Article
Conjugate Heat transfer in a Differentially Heated Porous Cavity Filled with Nanofluid
انتقال الحرارة المقترن خلال فجوة مسامية مشبعة بمائع نانوي و مسخنة تباينيا

Authors: Muneer A. Ismael منير عبد الجليل اسماعيل --- Ahmed Abdulkareem Mahdi احمد عبدالكريم مهدي
Journal: Basrah Journal for Engineering Science مجلة البصرة للعلوم الهندسية ISSN: Print: 18146120; Online: 23118385 Year: 2013 Volume: 13 Issue: 1 Pages: 123-139
Publisher: Basrah University جامعة البصرة

Loading...
Loading...
Abstract

The conjugate natural convection-conduction heat transfer in a domain composed of nanofluids filled porous cavity heated by a vertical solid wall is studied under steady-state conditions. The vertical left wall of the solid is kept isothermal at hot temperature Th. The vertical right wall of the solid is in contact with the nanofluid saturated porous medium contained in the cavity. The right vertical wall of the cavity is kept isothermally at the lower temperature Tc. The upper and lower horizontal walls are kept adiabatic. The governing equations of the heat transfer in the solid wall and heat and nanofluid flow, based on the Darcy model, in the nanofluid-saturated porous medium together with the derived relation of the interface temperature are solved numerically using the over-successive relaxation finite-difference method. A temperature independent nanofluids properties model is adopted. The investigated parameters are the nanoparticles volume fraction (0-0.2), Rayleigh number Ra (10-1000), solid wall to base-fluid saturated porous medium thermal conductivity ratio kwf (0.1, 1, 10), and the solid wall thickness D (0.05-0.5). The results are presented in the conventional form; contours of streamlines and isotherms and the average Nusselt number. At a very low Rayleigh number Ra=10, an enhancement in heat transfer within the porous cavity with is observed. Otherwise, the heat transfer may be unchanged or deteriorated with depending on the wall thickness D and the conductivity ratio kwf.

تم في هذا البحث دراسة انتقال الحرارة المقترن (حمل-توصيل) خلال حيز مكون من فجوة مسامية مشبعة بمائع نانوي و مسخنة بجدار عمودي تحت ظرف الأستقرار. سطح الجدار العمودي الأيسر مثبت عند درجة حرارة ساخنة Th و سطحه الأيمن في حالة تلامس مع المائع النانوي المشبع داخل الفجوة المسامية. جدار الفجوة الأيمن مثبت عند درجة حرارة باردة Tc . جدارا الفجوة الأعلى و الأسفل معزولان حرارياً. تم حل المعادلات الحاكمة لأنتقال الحرارة خلال الجدار الأيمن و معادلات انتقال الحرارة و جريان المائع النانوي (الذي بني على أساس نموذج دارسي) خلال الفجوة وكذلك المعادلة المشتقة للسطح المشترك Interface بأستخدام طريقة الفروقات المحددة. خواص المائع النانوي اعتبرت ثابتة مع درجة الحرارة. المتغيرات التي تمت دراستها هي النسبة الحجمية للجسيمات النانوية , رقم رايليRa=(10-1000) , نسبة الموصلية الحرارية للجدار العمودي الى المائع الأساس kwf=(0.1, 1, 10) وسمك الجدار العمودي D=(0.05-0.5). تم تمثيل النتائج بالطريقة المعروفة, أي خطوط دالة الجريان و التحارر و متوسط رقم نسلت. بينت النتائج بأنه عند قيمة Ra=10 هنالك تحسن في عملية انتقال الحرارة داخل الفجوة المسامية بزيادة النسبة الحجمية للجسيمات النانوية. فيما عدا ذلك, فأن انتقال الحرارة قد لا يتغيراو يتناقص مع اعتمادا على سمك الجدار D و نسبة الموصلية الحرارية kwf


Article
Study Natural Convection in a Porous Trapezoidal Cavity with a Square Body at the Center of the Enclosure
دراسة الحمل الحر داخل فجوة على شكل شبھ منحرف مملؤة بالمادة المتسامیة مع جسم مربع في مركز الغلاف

Author: Zainab K. Radhi
Journal: Engineering and Technology Journal مجلة الهندسة والتكنولوجيا ISSN: 16816900 24120758 Year: 2013 Volume: 31 Issue: 14 Part (A) Engineering Pages: 2632-2649
Publisher: University of Technology الجامعة التكنولوجية

Loading...
Loading...
Abstract

Steady state free convection heat transfer in a two dimensional trapezoidalenclosure filled with a fluid-saturated porous medium with a square solid bodylocated at the center of cavity, is performed in this study. The bottom wall of thecavity was heated with a sinusoidal temperature distributionq = l(1 - cos( 2pX )) ,the inclined walls are insulated (adiabatic) and the top wall is maintained at q = 0.To obtain the effects of the presence of a square body on heat transfer and fluidflow inside the enclosure, three different temperature boundary conditions wereapplied for the body as cold ( θbody=0) , heated (θbody=1) and adiabatic ( = 0)¶¶nq atdifferent Ra numbers. In this study, the governing equations were solvednumerically using finite element software package (FLEXPDE). Results for themean Nusselt number, Num, the contour maps of the streamlines and isotherms arepresented. It is observed that fluid flow and temperature fields strongly depend onthermal boundary conditions of the body. A comparison of the flow field, isothermfield and averaged Nusselt number, Num with previous work, which revealed agood agreement.


Article
Slip Effect on the Peristaltic Transport of MHD Fluid through a Porous Medium with Variable Viscosity
تأثير الانزلاق على الانتقال التموجي لمائع مغناطيسي هايدروديناميكي خلال وسط مسامي مع لزوجة متغيره

Authors: Ahmad M. Abdulhadi احمد مولود عبد الهادي --- Aya H. Al-hadad آية هادي الحداد
Journal: Iraqi Journal of Science المجلة العراقية للعلوم ISSN: 00672904/23121637 Year: 2015 Volume: 56 Issue: 3B Pages: 2346-2363
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

The present paper concerns with peristaltic analysis of MHD viscous fluid in a two dimensional channel with variable viscosity through a porous medium under the effect of slip condition. Along wave length and low Reynolds number assumption is used in the problem formulation. An analytic solution is presented for the case of hydrodynamic fluid while for magneto hydrodynamic fluid a series solution is obtained in the small power of viscosity parameter. The salient features of pumping and trapping phenomena are discussed in detail through a numerical integration. The features of the flow characteristics are analyzed by plotting graphs and discussed in detail. When .

هذا البحث يهتم بتحليل التموج لمائع مغناطيسي هايدروديناميكي لزج في قناة ثنائية مع لزوجة متغيره خلال وسط مسامي تحت تأثير شرط الانزلاق. فرضية طول الطول الموجي وعدد رينولد صغير قد استخدمت في هذه المسالة. قدم الحل التحليلي في حالة المائع الهايدروديناميكي بينما للمائع المغناطيسي الهايدروديناميكي قد حصلنا على سلسلة حلول لقوة صغيرة لعامل اللزوجة. السمات البارزه للضخ ومحاصرة الضواهر قد نوقشت بالتفاصيل من خلال تكامل عددي. حللنا ميزات خصائص الجريان بالرسوم ونوقشت بالتفاصيل. عندما

Listing 1 - 10 of 29 << page
of 3
>>
Sort by
Narrow your search

Resource type

article (29)


Language

English (24)

Arabic (3)

Arabic and English (2)


Year
From To Submit

2019 (1)

2018 (2)

2017 (2)

2016 (2)

2015 (2)

More...