research centers


Search results: Found 182

Listing 1 - 10 of 182 << page
of 19
>>
Sort by

Article
NUMERICAL SIMULATION OF A NATURALLY FRACTURED RESERVOIR.
المحاكاة العدديه للمكامن النفطية المتشققة

Authors: Mohammad S. Al-Jawad محمد صالح الجواد --- Ahmed J. Mahmood احمد جبير محمود
Journal: Journal of university of Anbar for Pure science مجلة جامعة الانبار للعلوم الصرفة ISSN: ISSN: 19918941 Year: 2010 Volume: 4 Issue: 2 Pages: 82-91
Publisher: University of Anbar جامعة الانبار

Loading...
Loading...
Abstract

This research is concerned with the development of a numerical model for stratified normally fractured reservoirs. Three dimensional three phase flow black oil simulation model is adopted. The dual porosity-dual permeability model is used. The IMPES (Implicit Pressure Explicit Saturation) method is used to solve the difference equations. The Tertiary trap in K oil field (an Iraqi oil field) was simulated by the numerical model. The trap consists of six layers having different properties. Equally spaced Cartesian grids were used to divide each layer into 1600 cells in the x-y plane with the thickness as the dimension of each grid block in z direction. Applying the two IMPES pressure equations to each grid of the simulated domain resulting in a block seven diagonal coefficient matrix.

في هذه الد ا رسة تم بناء نموذج محاكاة عددي یتناول حالات الجریان الثلاثي الأبعاد والثلاثي الطور في الأوساط المسامیةلحل معادلات الفروق IMPES المتشققة. اس تُخدم نموذج ثنائي المسامیة-ثنائي النفاذیه لوصف الوسط المسامي المتشقق،واستخدمت طریقةاحد الحقول النفطیة الع ا رقیة)والذي یتألف من ستة طبقات حاویه على النفط ) K الجزئیة . تم تطبیق النموذج على المكمن الثلاثي في حقلبكمیات متباینة إضافة إلى الغاز والماء . تختلف هذه الطبقات بالصفات الفزیاویه وشدة التشققات فیها.المكمن له قبة غازیه كبیره نسبیاوعمود نفطي متوسط إضافة إلى الماء في جوانب التكوین.بلغ عدد الآبار المحفورة في المكمن خمسة عشر بئ ا رً ولكل بئر تاریخ إنتاجمختلف. تم إدخال المعلومات الخاصة بالحقل للمودیل الریاضي و قُسِمت كل طبقه من طبقات المكمن الست إلى ١٦٠٠ خلیه في المستويالأفقي باستخدام شبكه ذات شكل غیر منتظم ومتساویة الأبعاد في الاتجاهین السیني والصادي مع اعتماد سمك ألطبقه كبعد بالاتجاهالعمودي.كانت . Gauss-Seidel على كل خلیه في الشبكة نتجت عنه مصفوفة تم حلها باستخدام طریقة التك ا رر (IMPES) إن تطبیق معادلتيوسیلة التحكم في الخطوة الزمنیة المستخدمة في الحسابات هما مقدار التغیر في نسب التشبع ومقدار الخطأ في المی ا زن المادي.تم استخدام تاریخ الإنتاج لخمسة عشر بئ ا ر المحفورة في الحقل في المودیل الریاضي للحصول على تصرف المكمن لفترة الإنتاج الممتدة منعام 1990 إلى عام 2005 أظهرت مقارنه معدل الضغط ألمكمني خلال مدة الإنتاج المحسوب بواسطة المودیل الریاضي والمقاس حقلیافي عدد من الآبار تطابقا جیدا. كما إن هناك تطابقا جیدا بین نسبة الغار إلى النفط المنتج المحسوب ریاضیا والمقاس حقلیا .


Article
Theoretical and Numerical Analysis of Fuel Droplet Combustion Process with Conventional and Transient Models

Authors: Ali Sabri --- Basma Abbas Abdul Majeed
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2008 Volume: 14 Issue: 4 Pages: 2940-2958
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

The study introduces a theoretical analysis and numerical solution for the combustion of a hydrocarbon fuel droplet inside the combustion chamber. The study employs two mathematical models to analyze the combustion process. Conventional (Classical) model and Transient model. The combustion process of a stagnant droplet in the steady state was analyzed in the classical model, while, In the transient model, it was assumed that there is a period of time in the stages of the droplet combustion in which the droplet is heated before combustion. The effect of change in temperature on the thermo-physical properties of the fuel was adopted through the two models. For the classical model, a convenient approximation was adopted for the heat transferred inside the droplet.

Keywords

droplet --- combustion --- numerical --- modeling --- transient.


Article
2D SEDIMENT TRANSPORT MODELLING TO THE SHATT AL-ARAB

Author: A.N.A. Hamdan
Journal: Journal of Basrah Researches (Sciences) مجلة ابحاث البصرة ( العلميات) ISSN: 18172695 Year: 2008 Volume: 34 Issue: 4A Pages: 50-59
Publisher: Basrah University جامعة البصرة

Loading...
Loading...
Abstract

This study is an attempt to determine the suspended sediment transport of Shatt Al-Arab estuary (North west Arabian Gulf).Two dimensional time dependent hydrodynamics model coupled with sediment transport model were applied and solved numerically by using the explicit finite difference method, computer program used to simulate the flow and suspended sediment transport which shows a good agreement with the satellite image .The two cases of tide (spring and neap tides) were applied to the estuary ,which shows that the spring tide case was more verified than the neap tide case regarding the situation at the time the satellite image was taken .

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


Article
Investigation into the Propagation Characteristics of Photonic Crystal Fiber
خصائص الانتشار في الليف الكرستالي الفوتوني

Authors: Raad H. Thaher رعد حمدان ظاهر --- Oras A.Shareef أوراس احمد شريف --- Sondus D.Yousif سندس داود يوسف
Journal: Al-Ma'mon College Journal مجلة كلية المامون ISSN: 19924453 Year: 2012 Issue: 20 Pages: 285-300
Publisher: AlMamon University College كلية المامون الجامعة

Loading...
Loading...
Abstract

Photonic crystal fiber (PCF) is a new class of optical fiber based on the properties of photonic crystals. It has the ability to confine light in hollow cores. In this paper, the propagation characteristics of PCF described, using different materials, was investigated such as effective area, effective refractive index , numerical aperture and material dispersion for three different materials . It is concluded that zero dispersion occurs at less wavelength when there is no silica in the constriction of PCF (i.e. 0.83µm for sapphire (extra ordinary wave). Also, the three dimension relation between dispersion ,bit error rate ,and fiber length for the three different materials was presented.

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

Keywords

BER --- Dispersion --- Numerical aperture --- PCF


Article
Numerical Simulations of Cloud Condensation Processes

Author: Enas A. Habeeb
Journal: Journal of Education for Pure Science مجلة التربية للعلوم الصرفة ISSN: 20736592 Year: 2019 Volume: 9 Issue: 2 Pages: 172-179
Publisher: Thi-Qar University جامعة ذي قار

Loading...
Loading...
Abstract

Abstract:Clouds are collections of tiny particles of water and/or ice that are large enough to be visible. Thetwo required ingredients for the formation of cloud are water vapor and aerosols. The spontaneous growthof Cloud Condensation Nuclei (CCN) into cloud droplets under supersaturated water vapor conditionsis described by classic Köhler theory. The aim of this work is to model Köhler theory and investigateeffects of mass, solubility, and composition of the CCN in the formation and growth of cloud droplets.Sodium Chloride (NaCl) and Ammonium Sulfate ((NH4)2SO4) were used as a solute material. The resultsshowed that both mass and solute material affect the nucleation of cloud droplet. Increasing mass of solutewould increase the critical radius and decrease the supersaturation required for activating the droplet. Itwas found that the use of NaCl produce critical radius larger than that of (NH4)2SO4 and the requiredsupersaturation for NaCl solute is lower than that of (NH4)2SO4 solute. This is due the fact that NaCl haslower molecular weight than (NH4)2SO4. The results also indicate that cloud temperature has minor effecton the cloud droplet activation. The model can be used to investigate more suitable solute materials todetermine which is best for cloud seeding projects.


Article
Developing Laminar Mixed Convection Heat Transfer Through Vertical Concentric Annuli With Adiabatic Inner Cylinder

Authors: Ghada A. Sadiq --- Akeel Abdullah Mohammed --- Ihsan Y. Hussain
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2011 Volume: 17 Issue: 5 Pages: 1289-1307
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

ABSTRACTTheoretical and experimental investigations have been carried out on developing laminar combined free and forced convection heat transfer in a vertical concentric annulus with uniformly heated outer cylinder (constant heat flux) and adiabatic inner cylinder for both aiding and opposing flows. The theoretical investigation involved a mathematical modeling and numerical solution for two dimensional, symmetric, simultaneously developing laminar air flows was achieved. The governing equations of motion (continuity, momentum and energy) are solved by using implicit finite difference method and the Gauss elimination technique. The theoretical work covers heat flux range from (200 to 1500) W/m2, Re range from 400 to 2000 and (1.36×105 ≤ Ra ≤ 1.1×107) with radius ratio of 0.555 which is the same radius ratio used in the experimental part of this study and Pr=0.7. The experimental work includes construct a rig consists essentially of an annulus with uniformly heated outer cylinder and adiabatic inner cylinder to give clear insight into heat transfer process and compare its results with that obtained in theoretical part, the range of the study are (Re= 383, 724, 1000, 1500) and heat flux equal to (q =370, 422, 588, 980) W/m2. Numerical results were represented by the temperature profile, axial velocity profile, outer surface temperature and the distribution of local Nusselt number along the dimensionless axial distance. The velocity and temperature profile results have revealed that the secondary flow created by natural convection have significant effects on the heat transfer process. Results reveal also that the experimental local Nusselt number along the annulus follows the same trend as present theoretical results with mean difference 10.23 %.

أُجريتْ دراسة نظرية وعملية لانتقال الحرارة بالحمل المختلط لجريان الهواء المتطور الطباقي خلال تجويف حلقي عمودي ذو أسطوانتين متمركزتين؛ الخارجية مسخنة تسخين منتظم ( فيض حراري ثابت) و الداخلية معزولة (لكلا الجريانين المساعد والمعاكس). البحث النظري تضمن عمل نموذج رياضي و دراسة عددية لجريان الهواء الطباقي المشكّـل تراكبيا ثنائي البعد. تم حل المعادلات الحاكمة للحركة (الاستمرارية , الزخم , الطاقة) باستخدام طريقة الفروقات المحددة و طريقة معكوس المصفوفة (طريقة كاوس). غطت الدراسة النظرية فيضا حراريا مداه من 200 واطم2 إلى1500 واطم2 , معدل رقم رينولدز Re من 400 إلى 2000 و (1.36×105 ≤ Ra ≤ 1.1×107) مع نسبة نصف قطر تساوي 0.555 وهي نفسها المستخدمة في الجزء العملي من هذا البحث ورقم براندتل =0.7 Pr. اشتمل الجزء العملي بناء جهاز تجريبي يتكون من تجويف حلقي متمركز بأسطوانة خارجية مسخنة تسخين منتظم و أسطوانة داخلية معزولة لاعطاء نظرة واضحة لعملية انتقال الحرارة ومقارنة نتائجها مع تلك التي تم الحصول عليها في الجزء النظري مدى الدراسة (Re= 383, 724, 1000, 1500) والفيض الحراري مساوي الى (q =370, 422, 588, 980) W/m2 . تم تمثيل النتائج النظرية برسم توزيع درجة الحرارة، توزيع السرعة المحورية على طول التجويف الحلقي، درجة حرارة الاسطوانة الخارجية و توزيع رقم نسلت الموقعي مع المسافة المحورية اللابعدية. اظهرت نتائج توزيع السرعة و درجات الحرارة ان الجريان الثانوي الناتج من الحمل الحر(الطبيعي) له تاثير مهم على عملية انتقال الحراة. النتائج بينت بان قيم رقم نسلت الموقعي العملية تتخذ نفس سلوك قيم رقم نسلت الموقعي النظرية لكن أعلى تقريباً مع معدل فرق 10.23 %.


Article
Stress Analysis of Composite Plates with Different Types of Cutouts
تحليل الاجهادات للألواح المركبة ذات الأنواع المختلفة من القطوعات

Authors: Ahmed N.E. أحمد نوري عويد --- Riyah N.K. رياح نجم كطر
Journal: Anbar Journal of Engineering Sciences مجلة الأنبار للعلوم الهندسية ISSN: 19979428 Year: 2009 Volume: 2 Issue: 1 Pages: 11-29
Publisher: University of Anbar جامعة الانبار

Loading...
Loading...
Abstract

This research presents an experimental and theoretical investigation of the effect of cutouts on the stress and strain of composite laminate plates subjected to static loads. The experimental program covers measurement of the normal strain at the edges of circular and square holes with different number of layers and types of composite materials by using strain gages technique under constant tensile loads. A numerical investigation has been achieved by using the software package (ANSYS), involving static analysis of symmetric square plates with different types of cutouts. The numerical results include the parametric effects of lamination angle, hole dimensions, types of hole and the number of layers of a symmetric square plate. The experimental results show good agreement compared with numerical results. It is found that increasing the number of layers reduces the value of normal strain at the edges of circular and square holes of a symmetric plate and the maximum value of stress occurs at a lamination angle of (30o) and the maximum value of strain occurs at a lamination angle of (50o) for the symmetric square plates subjected to uni-axial applied load. The hole dimensions to width of plates ratio is found to increase the maximum value of stress and strain of a symmetric square plate subjected to uniaxial applied load. Moreover, the value of maximum stress increases with the order of type of circular, square, triangular and hexagonal cutout, whereas the value of maximum strain increases with the order of type of circular, square, hexagonal and triangular cutout.

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


Article
NUMERICAL SIMULATIONS OF PARTICLE IN A DOUBLE OSCILLATORS
المحاكاة العددية لحركة جسيم في جهد توافقي ثنائي

Author: S. N. T. AL-RASHID سعيد نايف تركي
Journal: Journal of university of Anbar for Pure science مجلة جامعة الانبار للعلوم الصرفة ISSN: ISSN: 19918941 Year: 2007 Volume: 1 Issue: 3 Pages: 86-95
Publisher: University of Anbar جامعة الانبار

Loading...
Loading...
Abstract

In the present work, we extend the Hermann and Al-Rashid works to the problem of particle in a double oscillators potential. In this problem, one can take a special case when oscillation quantum number ( ) is none negative integer. Computer programming is built to make numerical simulations to this problem. The probability density of finding particle in a double oscillators potential is calculated without using Schrِdinger equation or any conventional quantum mechanics. This probability is compared with probability of conventional quantum mechanics.

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


Article
Numerical Study of Natural Convection in A Cavity with Wavy Vertical Walls

Author: SATTAR J. HABEEB
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2008 Volume: 14 Issue: 3 Pages: 2826 -2846
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

This paper describes a numerical study of natural convection heat transfer and fluid flow characteristics inside a cavity with wavy vertical walls. The bottom wall is heated by spatially varying temperature and other three walls are kept at cooled temperature. Governing equation was discretized using the finite volume-method with staggered variables arrangement in curvilinear coordinates. Two geometrical configurations were used in this study for symmetrical and unsymmetrical wavy vertical walls (total of 132 cases) for range of Ra=100 to 106 and fixed Prandtl number (0.71). The effects of the wave geometry, wave amplitude, number of undulation, and Rayliegh number on flow behavior, thermal field, local Nusselt number and Nusselt number ratio (NNR) factor have been studied. Streamline, velocity vector, and isothermal contour are used to present the corresponding flow and thermal field inside the cavity. The Results show that the enhanced of heat transfer rate seems to depend on geometrical configuration.

البحث الحالي يصف دراسة عددية لأنتقال الحرارة بالحمل الطبيعي وسمات جريان المائع داخل تجويف ذا جدران عمودية متموجة. الجدار السفلي مسخن بتغير حيزي لدرجة الحرارة ( T=F(x)) بينما بقية الجدران الثلاثة حفظت بدرجة حرارة ثابتة. المعادلات الحاكمة حٌلت بطريقة الحجوم الثابتة مع ترتيب متخالف للمتغيرات المدروسة بأحداثيات مطابقة الجسم. تم أخذ حالتان لتركيب الشكل الهندسي للتجويف, حالة التماثل, وحالة عدم التماثل للجدران العمودية المتموجة (كلياً 132 حالة) لمدى من Ra=100 الى 106 وثبوت Pr=0.71 , تم دراسة تأثير كل من الشكل الهندسي المتموج, مدى الموجة, عدد الموجات, Ra على تصرف الجريان والمجال الحراري داخل التجويف, Nu الموقعي, العامل NNR , تم الاعتماد على خطوط الانسياب, متجه السرع, خطوط ثبوت درجة الحرارة لأظهار الشكل المتكون للجريان والمجال الحراري داخل التجويف. النتائج أظهرت هناك تحسس لمعدل انتقال الحرارة يعتمد على تركيب الشكل الهندسي للتجويف.


Article
Numerical Simulation of Two Dimensional Transient Natural Convection Heat Transfer from Isothermal Horizontal Cylindrical Annuli

Authors: Jasim M. A. Al-lateef --- Ayad K. Hassan
Journal: Journal of Engineering مجلة الهندسة ISSN: 17264073 25203339 Year: 2007 Volume: 13 Issue: 2 Pages: 1429-1444
Publisher: Baghdad University جامعة بغداد

Loading...
Loading...
Abstract

Numerical solutions are presented for the transient natural convection heat transfer problem in horizontal isothermal cylindrical annuli, enclosed in heated inner and cooled outer cylinders. Solutions for laminar case were obtained within Grashof number based on the inner diameter which varied from 1x10^2 to 1x10^5 in air. Both vorticity and energy equations were solved using alternating direction implicit (ADI) method and stream function equation by successive over relaxation (SOR) method. The structure of fluid flow such as a velocity vector and temperature distribution as well as Nusselt number were obtained and the effect of diameter ratio on them is examined. In addition, the Grashof number was changed with the influence of variation Prandtle number and diameter ratio. Our numerical calculation are summarized by Nussult number vs. Grashof number curves with diameter ratios and prandtl as a parameter, which serves as a guide to natural convection heat transfer calculated from annulus. Good agreement with previous data were obtained.

يقَدَّم هذا البحث حَلّ عددي لمسالةِ انتقال الحرارةِ العابرةِ بالحمل ِ الطبيعيِ في العمود الحلقي الإسطوانيِ ثابت درجة الحرارةِ بوضع الأفقيِ، حيث تكونَ الاسطوانة الداخليةِ مسخنِة وَ الإسطوانة َ الخارجيةَ مبردة .تم حل المسالة لحالةِ الجريان الطباقي ولرقم كراشوف مستندة على القطرِ الداخليِ يتراوح بين 1x10^2 الى 1x10^5 للهواء. ِمعادلة الدوامية ومعادلة الطاقةِ حلتا بإستعمال بطريقة الإتّجاهَ الضمنيَ المتناوبُ ومعادلة الانسياب بطريقة فوق التراخي المتعاقبة.تم الحصول علىّ تركيبَ جريان المائعِ كمتجه السرعةِ وتوزيعِ درجةِ حرارة بالإضافة إلى رقم نسلت، تم اختبار تأثيرِ نسبةِ القطرِ عليها. بالإضافة الى دراسة تأثيرِ إختلافِ رقم كراشوف مع إختلافِ رقم برانتل و نسبة القطرَ عليهاَ. تم اجمال الحسابات العددية بواسطة منحنيات رقم نسلت مع رقم كراشوف و نسبة القطرَ ورقم برانتل. التي يمكن الاستفادة منها كدليل لحسابات انتقال الحرارة بالحمل الطبيعي من الشكل الحلقي.تم الحصول على توافق جيد بين النتائج الحالية وبياناتِ البحوث السابقةِ.

Listing 1 - 10 of 182 << page
of 19
>>
Sort by
Narrow your search

Resource type

article (181)

journal (1)


Language

English (142)

Arabic (19)

Arabic and English (17)


Year
From To Submit

2019 (14)

2018 (30)

2017 (18)

2016 (11)

2015 (9)

More...