research centers


Search results: Found 108

Listing 1 - 10 of 108 << page
of 11
>>
Sort by

Article
Theoretical Investigation of Starved Porous Journal Bearings Anisotropic Permeability with
دراسة نظرية لاداء المساند المقعدية ذاتية التزيت تحت ظروف عدم اكتمال شاء الزيت والنفاذية الغير منتظمة

Authors: Basim A. Abass باسم عجيل عباس --- Liqaa H. Abd-Alsheed لقاء حميد عبد الشهيد
Journal: AL-NAHRAIN JOURNAL FOR ENGINEERING SCIENCES مجلة النهرين للعلوم الهندسية ISSN: 25219154 / eISSN 25219162 Year: 2009 Volume: 12 Issue: 1 Pages: 13-25
Publisher: Al-Nahrain University جامعة النهرين

Loading...
Loading...
Abstract

The present study is a theoretical analysis to evaluate the static characteristics of oil – filled porous journal bearing with non – uniform permeability using an improved boundary condition. The analysis was carried out for a porous bearing with a permeability of the bearing ends is lower than that of the middle. The circumferential boundary condition for the oil film pressure is obtained by an integral momentum equation to the oil film region in the bearing clearance. It was found that a porous bearing with anisotropic permeability enhances the load carrying capacity of the bearing, also the bearing oil film is extended along specified extent

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

Keywords

Design --- Boundary layer


Article
Using boundary integral equation method to solve Plane strain problems of asymmetric elasticity
استخدام طريقة المعادلة التكاملية الحدودية لحل مسائل الإجهاد المستوي للمرونة غير المتماثلة

Authors: Auni M. Gaftan عوني محمد كفطان --- Muayad Sultan Wahyeb مؤيد سلطان وهيب
Journal: Tikrit Journal of Pure Science مجلة تكريت للعلوم الصرفة ISSN: 18131662 Year: 2013 Volume: 18 Issue: 3 Pages: 196-198
Publisher: Tikrit University جامعة تكريت

Loading...
Loading...
Abstract

In this paper we try to solve plane strain problems of asymmetric elasticity by using boundary integral equation method and we compare the results which we obtained from this method with the results of other methods we see that accuracy of our result is more than other results.

حاولنا في هذا البحث حل ( مسائل الإجهاد المستوي للمرونة غير المتماثلة ) باستخدام طريقة المعادلة التكاملية الحدودية , وتمت مقارنة النتائج التي حصلنا عليها عند الحل باستخدام هذه الطريقة مع نتائج الحل باستخدام الطرق الأخرى حيث لاحظنا إن النتائج الجديدة كانت أدق من النتائج الأخرى .


Article
Removal of a Single-Component NAPL under Nonequilibrium Condition: a Zero Mass Flux Moving Boundary
ازالة السائل ذو الطور اللامائي والمركب الواحد تحت ظروف عدم الأتزان: ظروف حدودية متحركة ذات تدفق كتلي صفري

Authors: Talib R. Abbas --- Mohammed A. Abdul-Majeed --- Inmar N.Ghazi
Journal: Engineering and Technology Journal مجلة الهندسة والتكنولوجيا ISSN: 16816900 24120758 Year: 2014 Volume: 32 Issue: 8 Part (A) Engineering Pages: 2096-2105
Publisher: University of Technology الجامعة التكنولوجية

Loading...
Loading...
Abstract

A new zero flux moving boundary conditions to solve the non-equilibrium single-component dispersion-advection transport problem containing NAPL have been introduced. The resultant model can overcome many approximations in the model which is done by previous study and give more precise analytical solution to the aqueous (or gas) phase concentration profile, NAPL saturation and the moving front speed. From the results, it is clear that the number of pore volumes (P), which are required for removing all NAPL under equilibrium conditions, plays a role in the differences between the two models' predictions. The NAPL saturation profile predicted by the two models is the same. At very low P the predicted value of the front speed by the present model is half of the value predicted by Brouwers and Augustijn model and increases gradually with P. It reaches 0.99 when P is 100.

تم أعتماد ظروف حدودية متحركة ذات تدفق كتلي صفري للحل الرياضي لمسألة الحمل-الأنتشار ولمركب واحد وتحت ظروف عدم الأتزان. يتجاوز الموديل الرياضي الناتج عدة تقريبات للحل الذي تم في دراسة سابقة. ويعطي حل تحليلي ادق لتركيز الطور المائي (الغازي) وكذلك لتركيز الطور السائل اللامائي وسرعة تقدم الجبهه المتحركة. تشير النتائج الى ان P ،عدد مرات حجم الفراغات في التربة من المائع النظيف واللازمة لأزالة كل الطور السائل الامائي وتحت ظروف الأتزان، يلعب دورا في الأختلافات بين تنبؤات كل من الموديلين. يعطي كلا الموديلين حلا متشابها لتركيز الطور السائل اللامائي. عند القيم الصغيرة جدا الى P فأن سرعة الجبهه التي يعطيها الحل في هذا البحث هي نصف القيمه التي يعطيها الحل السابق وتزداد هذه النسبة مع زيادة قيمة P. تصل قيمة هذه النسبة 0,99 عندما تكون قيمة P 100.


Article
Legendre Wavelets Method for Solving Boundary Value Problems
طريقة لجندر الموجية لحل مسائل القيم الحدودية

Authors: Dr. Suha N. Shihab --- Hayat A. Ali --- Bushra E. Kasheem
Journal: journal of the college of basic education مجلة كلية التربية الاساسية ISSN: 18157467(print) 27068536(online) Year: 2012 Volume: 18 Issue: 76 / ملحق Pages: 73-86
Publisher: Al-Mustansyriah University الجامعة المستنصرية

Loading...
Loading...
Abstract

Two techniques for solving nth order boundary value problem using continuous Legendre wavelets on the interval [0, 1] are presented. The first algorithm solves the boundary value problem BVP directly use the operational matrix of derivative of Legendre wavelets while the second algorithm converts the BVP into a system of Volterra integral equations then using the operational matrix of integration for Legendre wavelets, the system of integral equations is reduced to solve a set of linear algebraic equations, some examples are presented to illustrate the ability of the algorithms

تم تقديم اسلوبين لحل مسائل القيمة الحدودية من الرتبة النونية باستخدام متعددات لجندر الموجية المستمرة ضمن الفترة [0, 1]. الخوارزمية الاولى تحل المسألة مباشرة باستخدام مصفوفة العمليات المشتقة لمتعددات لجندر الموجية بينما الخوارزمية الثانية تحول المسألة الى منظومة من معادلات فولتيرا التكاملية ومن ثم استخدام مصفوفة العمليات التكاملية لنفس المتعددات, منظومة المعادلات التكاملية تختزل الى حل مجموعة من المعادلات الجبرية الخطية , بعض الامثلة قدمت لبيان كفاءة الخوارزميات


Article
High Resolution Biostratigrahpy of the K/T Boundary in the Higran Section, Shaqlawa Area, Northern Iraq
لطباقية الحياتية الدقيقة للحد الفااصل بين العصر الكريتاسي والتيرشري في مقطع هجران، منطقة شقلاوة، شمال العراق

Author: Rund A. Hammoudi رند حمودي
Journal: Iraqi National Journal of Earth Science المجلة العراقية الوطنية لعلوم الارض ISSN: 16823222 Year: 2011 Volume: 11 Issue: 1 Pages: 23-48
Publisher: Mosul University جامعة الموصل

Loading...
Loading...
Abstract

The Higran section in northern Iraq exposes the contact between two widely recognized formations. They are the Shiranish and Kolosh Formations, considered to be of Cretaceous (Campanian-Maastrichtian) and Tertiary (Paleocene-Eocene) ages respectively. Samples were taken from above and below the physical contact between the two formations at 5 to 7 cm apart. Planktonic foraminiferal biostratigraphic analysis of the Higran section indicates that the Cretaceous
/Tertiary boundary in this section is hitherto the most complete and expanded section in Iraq up to date. Quantitative high resolution foraminiferal biostratigraphic analysis has shown that a major biotic change in the planktonic foraminifera occurred during the end of the Cretaceous to the beginning of the Tertiary. The Maastrichtian Gansserina gansseri and the Abathomphalus mayaroensis Zones were recognized. The Abathomphalus mayaroensis Zone was subdivided into three subzones they are: Racemiguembelina fructicosa Subzone, Pseudoguembelina hariaensis Subzone and part of Pseudoguembelina palpebra Subzone. The defined Tertiary zones are: Guembelitria cretacea Zone, Parvularugoglobigerina eugubina Zone, Parasubbotina pseudobulloides Zone and part of the Parasubbotina varianta Zone. This study has established the absence of the Plummerita hantkeninoides Zone within the Cretaceous succession indicating a hiatus (Diastem) between the Cretaceous and Tertiary successions. The planktonic to benthonic ratio of the Shiranish Formation across the K/T boundary shows the role of the outer shelf to upper slope environments well above and below the K/T boundary in the section, hence it is characterized by lithologic continuity above and below the missing interval. The planktonic foraminiferal extinction in the Higran section occurred gradually and over relatively long period. Seventeen species (40.5% of the Cretaceous species) became extinct in the Late Maastrichtian before the first appearance of the Tertiary species, whereas nineteen species (45.2%) disappeared exactly with the hiatus that defines the K/T boundary in the studied section. Most of the extinct forms are large, complex, tropical and subtropical. Nevertheless, five species of the small cosmopolitan generalists with simple morphologies (14.3%of the Cretaceous species) survived the K/T boundary and the drastic change in the ecosystem. They did however become extinct in the Early Danian. The fine biostratigraphic analysis also revealed that the K/T boundary is located in the same lithology of blue marls within the Shiranish Formation, rather than at the contact between the Shiranish and Kolosh formations.This study has shown that the well-known Shiranish Formation conventionally taken to be of Cretaceous age, actually extends into the Paleocene in the Higran section. In essence, the turnover of the Cretaceous foraminifera which defines the end of the Cretaceous is located one meter below the physical contact between the Shiranish and the Kolosh formations.

يظهر مقطع هجران في شمال العراق الحد الفاصل بين تكوينين واسعي الانتشار في العراق هما تكوين شرانش ( الكامبانيان - الماساترختيان ) وتكوين كولوش ( الباليوسين – الايوسين ). جمعت النماذج لغرض هذه الدراسة من أعلى وأسفل الحد الطبيعي الفاصل بين التكوينين, المسافات بين النماذج تتراوح بين خمسة وسبعة سنتيمترات. الدراسة التحليلية البايوستراتغرافية لمتحجرات الفورامنيفرا الطافية اظهرت ان الحد الفاصل بين الكريتاسي والتيرشري في مقطع الدراسة يمثل اكثر مقطع تكاملا لهذه الترسبات الزمنية في العراق لحد الان. كما اظهرت الدراسة الطباقية الاحصائية لأنواع متحجرات الفورامنيفرا الطافية إن اكبر التغيرات الحياتية في أنواع الفورامنيفرا الطافية قد حصلت عند نهاية عصر الكريتاسي وبداية التيرشري. تم تحديد الانطقة الحياتية التالية من خلال الدراسة الطباقية لمتحجرات الفورامنيفرا الطافية وهي لعصر الكريتاسي:Gansserina gansseri Zone.Abathomphalus mayaroensis ZoneRacemiguembelina fructicosa Subzone.Pseudoguembelina hariaensis Subzone.Pseudoguembelina palpebra Subzone.أما الانطقة التي حددت للعصر التيرشري فهي:Guembelitria cretacea Zone.Parvularugoglobigerina eugubina Zone.Parasubbotina pseudobulloides Zone and part of t.Parasubbotina varianta Zone (Part).وقد أثبتت الدراسة الحالية غياب النطاق Plummerita hantkeninoides ضمن تتابعات الكريتاسي مما يدل على وجود فجوة زمنية صغيرة في تتابع ترسبات العصر الكريتاسي والتيرشري في مقطع الدراسة. إن نسبة الفورامنيفرا الطافية إلى القاعية ضمن تكوين شرانش في أعلى وأسفل حد الكريتاسي التيرشري أثبتت سيادة بيئة الرصيف الخارجي والمنحدر الاعلى في هذا المقطع مع تميزه بوجود استمرارية صخرية في أعلى وأسفل الحد الفاصل. إن انقراض أنواع الفورامنيفرا الطافية في أسفل الحد الفاصل كان تدريجيا وخلال فترة زمنية طويلة نسبيا. أثبتت الدراسة انقراض 17 نوعا من الفورامنيفرا الطافية (40.5%) قبل ظهور أول أنواع الفورامنيفرا العائدة للتيرشري في حين انقرض 19 نوعا تماما عند الحد الفاصل (45.2%) في مقطع الدراسة. معظم الأنواع المنقرضة تمثل الأنواع الكبيرة الاستوائية وشبه الاستوائية. وتم تشخيص خمسة أنواع من الفوارمنيفرا الطافية (14.3%) التي قاومت الظروف البيئية واستمرت خلال الحد الفاصل وبالتالي انقرضت ضمن ترسبات الدانيان المبكر.أثبتت الدراسة كذلك أن الحد الفاصل بين الكريتاسي والتيرشري يقع ضمن ترسبات المارل العائدة لتكوين شرانش, وليس عند الحد الفاصل بين تكوين شرانش وتكوين كولوش كما هو متعارف عليه في الادبيات السابقة. أثبتت الدراسة إن عمر تكوين شرانش المعروف تقليديا بالكريتاسي يمتد إلى الباليوسين في مقطع هجران, اعتمادا على وجود الحد البايوستراتغرافي الفاصل بين الكريتاسي والتيرشري أسفل الحد الطبقي العلوي لتكوين شرانش بمسافة متر واحد.


Article
Simulation Study of WATM Using MATLAB

Loading...
Loading...
Abstract

Several measures for evaluating the performance of wireless Asynchronous Transfer Mode (WATM) have been published in the literature. Among them: Offered Traffic Load, Handover Traffic Load, Call Admission Control and others. To the best of our knowledge, there is no study that collects all these measures into one GUI simulator. In the present work, eight measures have been discussed and then simulated into friendly used GUI using MATLAB. The proposed simulator enables the network administrator to: enter or select the values of the input parameters; plot and analyze the performance of WATM at different conditions. The graphical results, obtained from the simulator, can assist the administrator to take right decisions to improve the performance of the WATM.

تحتوي الادبيات على عدة معايير لتقييم أدائية شبكات ATM اللاسلكية ومن بين هذه المعايير حمل Offered Traffic وحمل Handover Traffic و تحكم Call Admission. حسب أخر معلوماتنا لا توجد دراسة جمعت كل هذه المعايير في برنامج واحد بواجهة مستخدم رسومية. في الدراسة الحالية أستخدم تطبيق MATLAB لدراسة وتمثيل ثمانية معايير في واجهة رسومية واحدة ، التمثيل المقترح يمكن مدير الشبكة من أدخال أو أختيار المتغيرات ورسم وتحليل أدائية الشبكة في ظروف عملية مختلفة. النتائج الرسومية للتمثيل المقترح تساعد مدير الشبكة في أتخاذ القرارات الصحيحة لتطوير أدائية شبكة ATM اللاسلكية.


Article
Natural Convection Heat Transfer from a Plane Wall to Thermally Stratified Environment

Authors: Naseem K. Ali --- Ihsan Y. Hussain
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2012 Volume: 18 Issue: 2 Pages: 223-238
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

The effect of linear thermal stratification in stable stationary ambient fluid on free convective flow of a viscous incompressible fluid along a plane wall is numerically investigated in the present work. The governing equations of continuity, momentum and energy are solved numerically using finite difference method with Alternating Direct implicit Scheme. The velocity, temperature distributions and the Nusselt number are discussed numerically for various values of physical parameters and presented through graphs. ANSYS program also used to solve the problem. The results show that the effect of stratification parameter is marginalized with the increase in Prandtl number, and the increase in Grashof number does not practically vary the effect of stratification parameter.


Article
Determination of Depth of Placement of Tunnels and Cavities by the Boundary Element Method

Authors: Rana A. J. Al-Adhamii --- Mohammed Yousif Fattah --- Omar Al-Farouk Al-Damluji
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2006 Volume: 12 Issue: 1 Pages: 53-70
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

A boundary element numerical algorithm has been developed for the determination of stresses and deformations around cavities and tunnels. A study of the influence of depth below the ground surface on the distribution of stresses and deformations around cavities and tunnels is presented in this paper. The soil is assumed to behave linearly elastic.A computer program has been built to perform the numerical computations. The results show that with increasing the depth of placement of tunnel or opening below the ground surface, the settlements decrease. The maximum stresses occur at the haunches of the tunnel rather than at the crown.For the circular cavity that is considered in this paper, it was found that with increasing the depth below the ground surface (depth/tunnel diameter > 3), the surface settlements do not exceed 6 % from those obtained for the case of no-cavity condition.

لقد طورت طريقة العناصر الحدودية كخوارزمية عددية لايجاد الاجهادات و الازاحات حول الانفاق و الفتحات الارضية. في هذا البحث أجريت دراسة حول تاثير عمق الفتحة أو النفق تحت سطح الارض على توزيع الاجهادات و الازاحات حول الانفاق أو الفتحات الارضية. و قد أعتبر تصرف التربة مرنا خطيا.و قد بني برنامج حاسبة الكترونية لاجراء الحسابات العدديه. وبينت النتائج أنه مع زيادة عمق الفتحه أو النفق تحت سطح الارض تقل الاجهادات و الهبوطات و أن الاجهادات القصوى تحدث عند الحافات الجانبية للنفق بدلا من قمته.و في حالة الفتحة الدائرية التي تم تحليلها في هذا البحث، وجد أنه مع زيادة العمق تحت سطح الأرض، (العمق / قطر النفق > 3)، لا تتجاوز الازاحات السطحية 6% من تلك التي تحصل في حالة عدم وجود فتحة.


Article
Boundary Element Analysis of Capped Pile Groups
التحلیل بطریقة العناصرالمتاخمة لمجامیع الركائز المربوطة بغطاء

Authors: Madhat S. Al-Soud --- Mohammed Y. Fattah --- Kais T. Shlash
Journal: Engineering and Technology Journal مجلة الهندسة والتكنولوجيا ISSN: 16816900 24120758 Year: 2009 Volume: 27 Issue: 11 Pages: 2321-2334
Publisher: University of Technology الجامعة التكنولوجية

Loading...
Loading...
Abstract

The boundary element method (BEM) has become one of the most powerfulnumerical techniques which has already established itself within the scientificcommunity. The most striking feature of this technique is that, in principle, only theboundaries of the region being investigated have to be discretized, which thereforeleads to many fewer discrete elements than any scheme requiring internal subdivisionof the whole body. This means that the number of unknowns is reduced dramatically,especially for 3D problems, as the unknowns occurred only on the boundary of theproblem.This paper is devoted to make use of the boundary element method (BEM) as apractical problem solving tool to analyze a soil - structure interaction problem. Theprogram (MRBEM) is adopted in this study for the analysis process. It is a generalpurpose boundary element method program for solving elasticity and potentialproblems with multiple regions. This program is written by FORTRN-90 language anddeveloped during this study to solve a three dimensional problem represented by agroup of piles. The results were compared with those findings in some experimentaland theoretical researches and good agreements were obtained.It was found that when using the BEM in the analysis, the stresses anddisplacements need only to be calculated where the details of interest occur on theboundary or are localized to a particular part of the domain, and hence an entiredomain solution is not required. Moreover, boundary conditions at infinity can bemodeled exactly without the need to extend the region a long distance away or to applyartificial boundary conditions as a result to the arbitrary truncation of the outer region.

واحدة من أقوى التقنيات العددية والتي أسست لنفسها أن (BEM) أصبحت طريقة العناصر المتاخمةتكون ضمن المجتمع العلمي. الصفة الأكثر تميزا" لهذه التقنية أن عملية تجزئة المنطقة المراد التحريعنها الى عناصر يقتصر على تخوم هذه المنطقة, الأمر الذي يؤدي الى أن عدد العناصر المجزئةيكون أقل من مخطط التجزئة الداخلي المطلوب للجسم كله. هذا يعني أن عدد المجاهيل سيقل بصورةهائلة, خصوصا" للمسائل ثلاثية الأبعاد, بسبب اقتصار ظهور المجاهيل عند تخوم المسألة.هذا البحث مكرس للأستفادة من طريقة العناصر المتاخمة كوسيلة حل عملية لتحليل مسائللأجراء عمليات التحليل في هذه (MRBEM) التداخل بين التربة والمنشآت. تم اعتماد البرنامجالدراسة. و هو برنامج متعدد الأغراض ويعتمد على طريقة العناصر المتاخمة لتحليل مسائل المرونةوالطاقة لمناطق متداخلة متعددة. كتب هذا البرنامج بلغة فورتران- 90 وقد تم تطويره خلال مرحلةالبحث لحل مسألة ثلاثية الأبعاد متمثلة بمجموعة من الركائز. تم مقارنة النتائج المستحصلة مع بعضالأبحاث العملية والنظرية وأثبتت توافقا" جيدا" مع هذه الأبحاث.وجد أنه عند استعمال طريقة العناصر المتاخمة فأن الاجهادات والازاحات سيقتصراحتسابها عند مناطق مهمة محددة على التخم أو أجزاء معينة ضمن مجال المنطقة دون الحاجة لأيجادحل لكامل الأجزاء الداخلية للمسألة. اضافة الى ذلك, فأن تمثيل الشروط التخومية للمسألة في المناطقالبعيدة يمكن تمثيلها بصورة دقيقة دون اللجوء الى بسط المنطقة لمسافات بعيدة وتطبيق شروطتخومية غير حقيقية نتيجة" للقطع العشوائي للمنطقة الخارجية


Article
The Effect of Boundary Air Flow on Premixed Stationary Burning Velocity
تأثير جريان الهواء المحيط باللهب الساكن مسبق الخلط على سرعة الاحتراق

Author: Ali Rasool M. Al-Zubaidi
Journal: Engineering and Technology Journal مجلة الهندسة والتكنولوجيا ISSN: 16816900 24120758 Year: 2009 Volume: 27 Issue: 14 Pages: 2611-2623
Publisher: University of Technology الجامعة التكنولوجية

Loading...
Loading...
Abstract

In this study the effect of boundary condition on premixed methane-airstationary Bunsen flames has been experimentally investigated. Laminar burning velocity is calculated by the concentrated cons method (CCM) and Schlieren photography technique, under the effect of laminar boundary horizontal stream air bulk with a range of air flow speeds. The experimental results have shown that the effect of these boundary conditions in general is small on calculated burning velocity if air flow speed is around between (0-50 cm/s). So, it is suggested thatthis effect be neglected. This effect so characterized increases at/around Stoichiometric ratio flames because of the increase of temperature difference between flame and boundary. The results of the experimental findings were compared with the latest published work and showed a good agreement with it, with a maximum discrepancy of (±2.5% at j =1.1).

تم في هذا البحث دراسة تأثير الظروف المحيطة بلهب بنسن الطباقي لخليط غاز الميثان مع الهواء مسبق الخلط. سرعة الاحتراق تم حسابها باستخدام طريقة المخروطين المتمركزين وتقنية تصوير شليرن تحت تأثير ظرف جريان هواء طباقي لسرع هواء مختلفة. (CCM) التجارب بينت وبشكل عام إن هذا التأثير على سرعة الاحتراق يكون صغيراً بحدود سرعة هواء 0-50 ) واقترحت الدراسة إهمال هذا التأثير. لكن هذا التأثير يزداد عند/حول نسبة الخلط cm/s) المثالية نتيجة زيادة الفرق بين درجة حرارة اللهب والمحيط. النتائج العملية لهذه الدراسة تم مقارنتها مع احدث الاعمال المنشورة واعطت توافق جيد معها. بنسبة خطاً.(j = اقصاها (% 2.5 ± عند 1.1

Listing 1 - 10 of 108 << page
of 11
>>
Sort by
Narrow your search

Resource type

article (108)


Language

English (93)

Arabic and English (7)

Arabic (6)


Year
From To Submit

2019 (13)

2018 (14)

2017 (8)

2016 (10)

2015 (8)

More...