computational fluid dynamic
در نشریات گروه مکانیک-
کار حاضر شبیه سازی کمک فنر مورد استفاده در تجهیزات فرود یک هواپیما می باشد. تمامی معادلات به روش المان محدود حل میگردند و نوع شبکه بندی بی سازمان برای همگرایی سریعتر بهره گرفته میشود. در گام نخست برای اعتبارسنجی مدل مورد استفاده مقایسه ای با دادههای تجربی موجود صورت میگیرد که خطایی نزدیک 5 درصد را نشان میدهد که خطای مطلوبی میباشد. نتایج نشان میدهند که با فشار آمدن به ارابه فرود جابجایی آغاز میشود که در زمان 05/0 ثانیه این امر مشاهده میگردد. در این حالت بیشترین تنش به پایه قلابی شکل ارابه فرود وراد میآید. با مقایسه این حالت با حالت قبلی دیده میشود که حداقل ده برابر شده است. با ورود سامانه جاذب شوک به سیستم این تنش کاهش مییابد. کاهش تنش بیشینه از 108*3 از به 108*6/1 پاسکال بیانگر همین مطلب میباشد. با جابجایی ارابه فرود کاهش تنش ادامه مییابد. همچنین دیده میشود که دامنه جابجایی نسبی در ابتدا بسیار زیاد است و حتی تا 35 سانتی متر هم میرسد اما به دلیل اتلاف انرژی در کمک فنر و لاستیک ها، به سرعت از بین میرود. به طوریکه بعد 45/0 ثانیه تقریبا زیر 5 سانتی متر است و بعد از 7/0 ثانیه میرا میشود. سرعت نسبی نیز در ابتدا بسیار زیاد است و تا 10 متر بر ثانیه نیز پیش میرود ولی به دلیل اتلاف انرژی در کمک فنر و لاستیکها در طول فرود این رقم به 5 متر بر ثانیه رسیده و نهایتا بعد از 45/0 ثانیه به زیر 1 متر بر ثانیه کاهش یافته و در 7/0 ثانیه به صفر میرسد.
کلید واژگان: تحلیل تنش، انتقال حرارت، کمک فنر چرخ هواپیما، دینامیک سیالات محاسباتیPresent work is a simulation of shock absorbers of aircraft landing gear. Equations are solved by finite element method and grid type is unorganized. Comparison with experimental data shows an error of about 5%. The results show that the displacement starts when the landing gear is pushed, which is observed in 0.05 seconds. In this case, the greatest stress comes to the base of the landing gear. By comparing this state with the previous state, it can be seen that it has increased at least tenfold. With the introduction of the shock absorber system, this tension is reduced. The reduction of the maximum stress from 3*108 to 1.6*108 pascal indicates the same thing. As the landing gear moves, the tension continues to decrease. It can also be seen that the range of relative displacement is very large at the beginning and even reaches 35 cm, but due to energy loss in the shock absorber and tires, it quickly disappears. So that after 0.45 seconds it is almost below 5 cm and after 0.7 seconds it is damped.
Keywords: Stress Analysis, Heat Transfer, Airplane Wheel Shock Absorber, Computational Fluid Dynamic -
In this study, the effects of the structural geometries of the channels used for cooling blades in the heated regions of the gas turbine exposed to high temperature were investigated. It is aimed to cool the gas turbine blade using 10 different types of ribbed channels by simulation method. Roof, inverted roof, slope and wedge ribs with standard and stepped arrangements developed for improving thermal performance in a rectangular cooling duct with a 4:1 ratio aspect were studied. Square type rib, the basic geometry is designed to make thermal comparisons. Details of turbulent flow and numerical calculations were made using of the standard k-ε turbulence model by means of the Computational Fluid Dynamics (CFD) method. The heat transfer in ribbed walls was examined in four different Reynolds numbers, 10000, 20000, 40000 and 80000 according to the channel input cross section. As a result of the calculations, the temperature changes of the turbine blade depending on Re number, the heat transfer improvement occurring in the internal channels inside the blade and the overall thermal performance were compared. The different types of rib evaluated in the study were compared with the standard-square rib; higher Nu number was obtained in stepped-wedge rib, stepped-roof rib, stepped-slope rib and stepped-inverted roof rib, respectively. It has been observed that a stepped-wedge rib can improve overall thermal performance and is promising for internal turbine blade cooling applications. The best operating range of all models was found to be between Re=40000 and Re=80000. The highest PEC results were obtained in the stepped-wedge rib model. This is 3.37% higher than the closest performing stepped-inverted roof rib model with 5.04 at Re=8000.Keywords: Gas Turbines, Blade Cooling, Computational Fluid Dynamic, Rib, Optimization
-
To assess the impact of different heights of the viaduct on air pollutant dispersion and the prediction accuracy of pollutant concentration in urban street canyons, simulation results based on LES and RANS models are analyzed. The presence of a viaduct generated a poorly ventilated region underneath it, and RANS significantly underestimated the wind speed and grossly overestimated the pollutant concentration. LES gives better results for the flow pattern, distribution of turbulent kinetic energy and mean pollutant concentration. With a fluctuation of less than 15% of the pollutant concentration, both RANS and LES cases show that an increase in the viaduct height has a weak impact on the concentration of pollutants in most areas of the canyon except windward, and cases with a viaduct height of 0.75H had the lowest predicted pollutant concentration relative to other cases with a viaduct as a result of better ventilation. In addition, LES found a subregion of pollutant accumulation above the ground, but RANS did not.Keywords: Street Canyon, Viaduct, Pollution Distribution, Large Eddy Simulation, Computational Fluid Dynamic
-
بیشینه سازی راندمان آیرودینامیکی نسبت به دامنه امواج دریا با استفاده از تلفیق نتایج عددی و شبکه عصبیهدف از این مقاله با توجه به گسترش روز افزون استفاده از هواپیماهای اثرسطحی و پرواز این نوع وسیله بر روی آب های آزاد به واسطه استفاده از پدیده اثر سطح، بررسی تاثیر امواج دریا بر روی ضرایب آیرودینامیک هواپیما به منظور بیشینه سازی راندمان آیرودینامیکی می باشد. جهت انجام این کار، شبیه سازی سه بعدی هواپیمای کامل در شرایط پروازی عبور از امواج دریا توسط دینامیک سیالات محاسباتی انجام گرفت و بر اساس نتایج حاصل از آن مدلسازی شبکه عصبی برای ماکزیمم راندمان آیرودینامیکی نسبت به دامنه امواج دریا صورت پذیرفته است. تلفیق دینامیک سیالات محاسباتی و شبکه عصبی نشان می دهد که هواپیما جهت پرواز در مجاورت سطح دریا با دامنه امواج مختلف بایستی در سرعت خاصی پرواز کرده تا به ماکزیمم راندمان آیرودینامیکی برسد. این نتایج می تواند به عنوان یکی از ورودی های سیستم کنترلی جهت استفاده در خلبان خودکار هواپیما مورد استفاده قرار گیرد.کلید واژگان: امواج دریا، اثر زمین، دینامیک سیالات محاسباتی، هواپیمای اثرسطحی، ضرایب آیرودینامیکی، شبکه عصبی، الگوریتم ژنتیک، سیستم کنترل پروازThe purpose of this article is to maximize the aerodynamic efficiency of an airplane during flight in the vicinity of the sea waves. A whole 3D aircraft flight at the closed proximity to the wavy sea water has been numerically simulated by the Fluent commercial software and the results were used to train a neural network. This network was then employed to determine the maximum aerodynamic efficiency in the vicinity of sea waves. The combination of the computational fluid dynamics and the neural network has approved that the flight must be performed at a certain speed to attain the maximum aerodynamic efficiency with a respect to a varied range of sea waves amplitude. The results can be considered as the inputs for autopilot system to provide the most efficient flight over the sea waves.Keywords: Sea Waves, Ground Effect, Computational Fluid Dynamic, WIG, Aerodynamic Coefficients, Neural Network, Genetic Algorithm, FLIGHT CONTROL SYSTEM
-
بررسی عددی و تجربی روش برآمدگی های موجی شکل هم زمان در بهبود عملکرد آیرودینامیکی ملخبهبود عملکرد ملخ به عنوان یکی از اجزاء اصلی پیشران پرنده ها همواره درحال توسعه قرار دارد. در این میان استفاده از راهکاری باقابلیت استفاده در گستره بالایی از انواع ملخ و با امکان حفظ پارامترهای طراحی اولیه مفید است. در این تحقیق بررسی روش بهبود بازدهی لبه های موجی شکل هم زمان با ایده اولیه الهام از طبیعت و مبنای ایجاد تغییر گام ذاتی در طول پره ملخ چوبی با قطر 7/45 سانتی متر مبتنی بر روش های عددی و تجربی انجام شده است. در روش عددی از راهکار شبیه سازی قاب مرجع چرخان و حل معادلات ناویر - استوکس تراکم پذیر پایا با محاسبه عبارات تنش در معادله مومنتوم از روش میانگین عدد رینولدز و روابط K-ω-SST استفاده شده است. این تحلیل ضمن بررسی استقلال نتایج از ابعاد و تعداد المان محاسباتی میدان های با قاب مرجع چرخان و ثابت اطراف ملخ نسبت به نتایج، دقت مناسبی در تخمین ضرایب عملکردی در مقایسه با آزمایش تجربی داشته است. آزمایش تجربی در تونل باد زیر صوت مرکز تحقیقات آیرودینامیک قدر با نسبت سطح ملخ به مقطع آزمون 2/0 برای محدوده ای از نسبت پیشروی عملکردی ملخ انجام شد. داده برداری ضریب گشتاور و تراست با استفاده از بالانس متصل به موتور براشلس نسبت به تغییرات سرعت چرخش و سرعت جریان در تونل باد انجام شده است. نتایج افزایش بازدهی ملخ با مقادیر طول موج 10 میلی متر و دامنه 5 درجه در مقایسه با ملخ اولیه معادل 5/5 درصد در نسبت پیشروی مطابق با محدوده پروازی را نشان داد.کلید واژگان: بهبود عملکرد ملخ، برآمدگی های موجی شکل هم زمان، دینامیک سیالات محاسباتی، آیرودینامیک تجربیNumerical and experimental investigation of synchronic wavy tubercles method in improving propeller aerodynamic performanceThis research aims to enhance propeller performance, a critical aspect of avian propulsion systems, without altering the original design of various propeller types. The study explores a method to improve the efficiency of wavy edges, inspired by nature, by creating an inherent pitch change along a 45.7 cm diameter wooden blade. Both numerical and experimental methods were utilized. The numerical approach involved simulating a rotating reference frame and solving the Navier-Stokes equations using the K-ω-SST method to calculate stress expressions in the momentum equation. The analysis demonstrated good accuracy in estimating performance coefficients compared to experimental tests. The experimental tests were conducted in a subsonic wind tunnel, with a propeller surface to test section ratio of 0.2, covering a range of the propeller’s functional advance ratio. Data on torque and thrust coefficients were acquired using a balance connected to a brushless motor. Results indicated an increase in propeller efficiency with specific wavelength values and amplitudes, showing a 5.5% increase in the advance ratio according to the flight range compared to the original propeller.Keywords: Propeller Performance Improvement, Synchronic Wavy Tubercles, Computational Fluid Dynamic, Experimental Aerodynamic
-
از دیرباز عدم قطعیت به عنوان عامل غیرقابل انکار در مطالعه پدیده های فیزیکی موردتوجه بوده است. ازاین رو جهت حصول نتایج قابل اطمینان در کاربردهای مهندسی از طریق پیش بینی های محاسباتی، تمامی منابع عدم قطعیت موجود در سیستم باید در نظر گرفته شود. پیشرفت های اخیر در دینامیک سیالات محاسباتی این امکان را فراهم کرده است که اثرات عدم قطعیت بر میدان های جریان و انتقال حرارت پیچیده در کاربردهای مهندسی موردمطالعه قرار گیرند. در مقاله حاضر آنالیز عدم قطعیت میدان جریان گذر صوتی اطراف ایرفویل RAE 2822 به عنوان یک مسئله چالشی در حوزه دینامیک سیالات محاسباتی تحت اثر عدم قطعیت هندسه ایرفویل ناشی از تلرانس های ساخت ارزیابی شده است. در ابتدا کد توسعه یافته بسط چندجمله ای آشوب بر روی چند تابع چالشی و غیرخطی اعتبارسنجی گردیده است. سپس برای ایجاد عدم قطعیت هندسی در ایرفویل RAE 2822، روش کارهونن لوو (KL) با بهره گیری از تعداد 18 متغیر تصادفی و انجام 2660 شبیه سازی مختلف به کار گرفته شده است. نتایج به دست آمده از ضریب فشار غیرقطعی اطراف ایرفویل نشان می دهند که عدم قطعیت های هندسی با شدت بیشتری مکان وقوع و قدرت شوک نرمال را تحت تاثیر قرار می دهند.کلید واژگان: کمی سازی عدم قطعیت، بسط چندجمله ای آشوب، روش های کارا، تحلیل حساسیت، دینامیک سیالات محاسباتیUncertainty has been known as an unavoidable parameter since early steps of investigation on physical phenomena. Therefore, to ensure validity of numerical simulations in the engineering applications, all of the related uncertainty sources must be considered. The examination of uncertainty effects on the flow field and heat transfer in the complex applications have been possible by recent advances in the computational fluid dynamics (CFD) methods. In this paper, the analysis of the uncertainty quantification (UQ) of the transonic flow field around the RAE 2822 airfoil is evaluated as a challenging problem in the CFD field under the effect of airfoil geometry uncertainties caused by manufacturing tolerances. At first, the developed code of polynomial chaos expansion (PCE) is validated on the several challenging and nonlinear test functions. Then, to construct geometrical uncertainties in the RAE 2822 airfoil, the Karhunen-Loeve (KL) method is employed by using 18 random variables and performing 2660 different CFD simulations. The obtained results from the non-deterministic pressure coefficient around the airfoil show that the geometric uncertainties more strongly affect the place of occurrence and the strength of the normal shock.Keywords: Uncertainty quantification, Polynomial chaos expansion, Efficient methods, Sensitivity analysis, computational fluid dynamic
-
A vital process for converting heavy petroleum productions is Fluid Catalytic Cracking (FCC). As a major source of CO2 emissions, the regenerator reactor in the FCC unit accounts for about 20-35% of the refinery's total emissions. A common method for reducing CO2 emissions from the FCC regenerator is oxy-combustion, which has different advantages with regard to reducing energy penalties and associated costs. In this study, a computational fluid dynamic (CFD) study was used to examine the hydrodynamic characteristics of solid particles and gas inside the FCC regenerator, allowing CO2 to be captured more efficiently. Utilizing Ansys Fluent platform, the Eulerian-Eulerian model was applied with granular flow kinetic theory. In the simulations, different mesh sizes were tested, and the hydrodynamics of the oxy-combustion regenerator were evaluated by adjusting CO2 flow rates to achieve similar fluidization behaviors. The CFD results indicated that the conventional drag model accurately predicted the density phases within the bed. In oxy-combustion, CO2, due to its density, naturally creates a smaller dense phase compared to air-combustion. Moreover, optimizing the fluidizing gas velocities resulted in enhanced particle mixing, resulting in a distributed flow with vortices within the dense phases due to a reduction in gas velocity. To improve the environmental performance of the FCC unit, this research provides valuable insight into the hydrodynamics of solid catalysts used in the oxy-combustion process.Keywords: Oxy-combustion, Computational fluid Dynamic, Hydrodynamic, Multiphase flow, Regenerator
-
رفتار جابجایی نفت و میزان به دام افتادن آن در ناحیه اطراف چاه توسط سیال پیش شوی فرآیند اسیدکاری حایز اهمیت است. در این مطالعه، رفتار ویسکوز و مویینگی جریان دو فازی در مقیاس منفذی با استفاده از روش دینامیک سیالات محاسباتی مورد تجزیه و تحلیل قرار گرفته است. در این کار یک مدل دو بعدی، بر اساس معادلات میدان فازی کان هیلیارد و ناویر استوکس، ایجاد و با استفاده از روش اجزای محدود حل شد. با توجه به اهمیت نیروهای حاکم در این جابجایی، نمودار فازی براساس عدد مویینگی و نسبت ویسکوزیته تهیه و سپس با کارهای تجربی گزارش شده قبلی مقایسه گردید. پس از شناسایی حداکثر و حداقل دامنه عدد مویینگی (Nc) و نسبت ویسکوزیته (M)، پایدارترین منطقه جابجایی در Log M ≈ 0.5 و Log Nc ≈ -2 قرار گرفت. علاوه بر این، تاثیر چهار متغیر مستقل شامل حجم تزریق (5>1>PV)، عدد مویینگی (0>Log Nc<-6)، نسبت ویسکوزیته (20.95) در Log M> 0 و 1 π/6 و(2->Log Nc<-4.5) اتفاق می افتد. شایان ذکر است که برای Log M<0، شرایط بهینه در Log M ≈ 0، Log Nc ≈ -3.5، PV ≈ 4 و θ ≈ π/6 رخ می دهد.کلید واژگان: اسید کاری، سیال پیش تزریق، دینامیک سیالات محاسباتی، آنالیز حساسیت، عدد موئینگیOptimum oil displacement is the most important goal of pre-flush stage during matrix acidizing. In this study, the visco-capillary behavior of the two-phase flow in the pore-scale is analyzed using computational fluid dynamics. A two-dimensional model, based on Cahn–Hilliard phase-field and Navier–Stokes equations, was established and solved using the finite element method. To recognize the effective forces of two-phase flow displacement, a stability phase diagram based on the Logarithm of capillary number (Nc) versus the Logarithm of viscosity ratio (M) was constructed and then compared with the reported experimental data. After identifying both viscous fingering and capillary fingering regions, the most stable displacement region was found to be located at Log M ≈ 0.5 and Log Nc ≈ -2. Furthermore, the impact of four independent variables that critically affect the efficiency of the pre-flush stage, including pore volume of injection (1<PV<5), capillary number (-60.95) occurred at Log M > 0, -4.5 1, and θ> π/6 conditions. It is worth mentioning that for Log M< 0, the optimum condition occurred at Log M ≈ 0, Log Nc ≈ -3.5, PV ≈ 4 and θ ≈ π/6.Keywords: Acidizing, Pre-flush fluid, Computational fluid dynamic, Sensitivity analysis, Capillary number
-
استفاده از پنجره های بزرگ در ساختمان های مدرن یکی از طرح های محبوب در بین معماران و طراحان ساختمان است ولی افزایش مصرف انرژی در ساختمان به دلیل تشدید درفت ناشی از پنجره های بزرگ از جمله محدودیت ها در مقابل اجرای این طرح می باشد. هم چنین عدم امکان استفاده از سیستم های گرمایشی سنتی در محل پنجره ها از مشکلات دیگر فرا روی اجرای این طرح می باشد. در این مطالعه رویکرد استفاده از سیستم گرمایش جابه جایی اجباری با هدف کاهش اثرات کوران در قالب سه چیدمان مختلف از دمنده هوای گرم که دارای جانمایی الف) در زیر پنجره، ب) در دیوارهای جانبی و ج) به صورت ترکیبی در زیر و بالای پنجره قرار گرفته است بااستفاده از تکنیک دینامیک سیالات محاسباتی شبیه سازی شده است. بررسی نتایج بااستفاده از مقایسه توزیع دما و توزیع سرعت هوا برای هر سه چیدمان انجام شد و هم چنین آنالیز آسایش حرارتی بااستفاده از مدل فنگر برای هر سه حالت انجام گرفت. نتایج نشان دادند که حالت سوم با داشتن دمنده ترکیبی که در آن از دمنده هایی در پایین و بالای پنجره استفاده شده است دارای مطابقت قابل قبولی با استانداردهای آسایش حرارتی اشری و ایزو 7730 است، به طوری که در ارتفاع حساس مچ پا و ناحیه تنفسی در حالت نشسته و ایستاده، شاخص های مورد تحلیل در آسایش حرارتی در محدوده تعریف شده استاندارد قرار گرفته است.کلید واژگان: آسایش حرارتی، کوران، پوسته شیشه ای، دینامیک سیالات محاسباتیEmploying large windows in the modern buildings is one the most popular approach among architects and buildings designers but the increase of energy consumption of the building due to down draught effects caused by large windows is one of the challenges against this type of design. Also, preventing of utilizing traditional heating system adjacent to the window is another concern of the approach. In this study, the use of forced convection heating system to mitigate the down-draught effects is simulated by computational fluid dynamic method for three different configurations of the warm air registers mounted (a) under the window; (b) alongside walls and (c) in mixed form under and above the window. The results were investigated by comparing air temperature and velocity distribution in the space and also the thermal comfort conditions were analyzed through Fenger model. The results showed that the third configuration of the inlet registers which they are installed beneath and above the window, is in satisfactory concord to ASHRAE standards and ISO 7730 such that at ankle and breathing height of sitting and standing position, the comfort condition indexes are in regular range of the standards.Keywords: Thermal comfort, down-draught, glazing envelope, computational fluid dynamic
-
The sound generated by axial flow fans has become increasingly important in several industrial areas. These devices represent considerable noise sources that must be considered from the design stage. The scope of this research work is to present an axial fan design methodology based on the theory of the wing lift and to compare the aerodynamic and aeroacoustics behavior with the free vortex and the non-free vortex conditions. The effect of forward circumferential sweep on the rotor blades is analyzed in both conditions, which consists of a displacement in the tangential direction in the direction of rotor rotation. A cubic polynomial function is proposed for the construction of the blade geometry with circumferential sweep. It defines the center line of the blade and it is taken at 50% of the leading edge of each blade profile. The methodology for analyzing aeroacoustics behavior in fans is based on the integration of Computational Aeroacoustics and Computational Fluid Dynamics in order to analyze aerodynamic behavior, sound power level and sound pressure level at the preliminary design stage. It was verified that although both rotors present a similar aerodynamic behavior (the NFV rotor has a slightly lower performance), according to the analysis of local sound sources, the NFV rotor presents a slight decrease in the sound sources. It was also verified that the NFV rotor presents 2.2 dB less in sound pressure level than the FV rotor.Keywords: CFD analysis, Design, Axial flow fan, Circumferential sweep, Aeroacoustics, Computational fluid Dynamic, free vortex, Non-Free Vortex
-
با ظهور ایرفویل های تغییرشکل پذیر، آیرودینامیک ایرفویل توربین های بادی و بال ها دچار تغییرات زیادی شد. در این پژوهش ضرایب آیرودینامیکی ایرفویل تغییر شکل پذیر بر مبنای ناکا0015 در محدوده عدد رینولدز 105 تا 106 و زاویه حمله ی 0 تا 12 درجه به کمک شبکه عصبی مصنوعی و الگوریتم ژنتیک بهینه سازی شده است. ابتدا ایرفویل ها به وسیله نقاط کنترل تصادفی در نرم افزار متلب تولید و در نرم افزار گمبیت شبکه بندی شدند، سپس در نرم افزار انسیس به صورت دو بعدی شبیه سازی شدند. نتایج شبیه سازی شامل ضرایب برآ و پسآ، نقطه جدایش و مرکز فشار به همراه نقاط کنترل برای آموزش شبکه عصبی مورد استفاده قرار گرفت. تابع آموزش دیده شبکه عصبی به عنوان تابع ورودی به الگوریتم ژنتیک داده می شود تا ضرایب مورد نظر بهینه سازی شوند. ضریب برآ، مرکز فشار، نقطه ی جدایش و نسبت ضریب برآ به پسآ به صورت تک هدفه بهینه سازی شدند، در بهینه سازی تک هدفه ضریب برآ با استفاده از ایرفویل تغییرشکل پذیر مقدار ضریب برآ 18% افزایش یافت. همچنین ضریب برآ و مرکز فشار، ضریب برآ و پسآ به صورت دو هدفه بهینه سازی شدند. در بهینه سازی دو هدفه ضریب برآ و پسآ، با تغییر شکل ایرفویل مقادیر آنها به ترتیب در 0/8 و 0/03 کنترل شد.
کلید واژگان: ایرفویل تغییرشکل پذیر، دینامیک سیالات محاسباتی، شبکه عصبی مصنوعی، الگوریتم ژنتیکWith the advent of morphing airfoils, the aerodynamics of wind turbines and wings underwent many changes. In this study, the aerodynamic coefficients of morphing airfoil based on NACA 0015 are optimized in the range of Reynolds number 105 to 106 and the angle of attack 0 to 12 degrees using Artificial Neural Network (ANN) and Genetic Algorithm(GA). First, the airfoils were created in MATLAB software by random control points and mesh generated in Gambit software, then in two-dimensional Ansys software were simulated. The simulation results, including lift and drag coefficients, separation point and pressure center, with control points were used to train the Artificial Neural Network (ANN). The trained function is given as an input function to the Genetic Algorithm(GA) to optimize the desired coefficients.Lift coefficient, the center of pressure, separation point and lift to drag ratio were optimized as a single objective, In single-objective optimization, the lift coefficient was increased by 18% using the morphing airfoil. Also, the lift coefficient and the center of pressure, the lift coefficient and the drag coefficient were optimized as the dual-objectives optimization. In the optimization of the dual objectives, lift and drag coefficients were controlled by 0.8 and 0.03, respectively, by the morphing airfoils.
Keywords: Morphing Airfoil, Computational Fluid Dynamic, Artifical Neural Networks, Genetic Algorithm -
تحقیقات در مورد استفاده از مواد هوشمند در بال هواپیماها برای افزایش عملکرد آنها و به دنبال آن استفاده از مواد هوشمند در ایرفویل توربین های باد سال هاست که آغاز شده است و نتایج این پژوهش ها نیز مطلوب بوده است. در این پژوهش، از دینامیک سیالات محاسباتی و معادلات URANS برای توربین بادی داریوس مجهز به ایرفویل تغییر شکل پذیر به منظور تعیین مقطع بهینه پره استفاده شده است، 250 ایرفویل توسط نقاط کنترل تصادفی تولید، درنرم افزار گمبیت، به طور غیرسازمان یافته و به صورت شبکه لغزان شبکه بندی شدند، سپس در انسیس به صورت دو بعدی با استفاده از مدل آشفتگی k-ω SST و الگوریتم PISOشبیه سازی شدند، نقاط کنترل و ضریب توان برای آموزش شبکه عصبی مورد استفاده قرار گرفتند و از الگوریتم ژنتیک برای بهینه سازی ضریب توان بهره گرفته شد. در این پژوهش از ایرفویل ناکا0015 به عنوان ایرفویل مبنا استفاده شده است، بهینه سازی برای دوحالت مختلف انجام شد. برای حالت اول(تعیین مقطع بهینه توربین در یک دور) میزان ضریب توان توربین بادی داریوس با مقطع بهینه واحد به میزان 42% افزایش یافت که مقطع پره (ایرفویل) نیز ترسیم شد. برای حالت دوم (تعیین مقطع بهینه در هریک از چهار ناحیه چرخش روتور) نیز بهینه ترین مقاطع (ایرفویل) بدست آمدند، افزایش 60% ضریب توان توربین نتیجه ی بهینه سازی برای حالت دوم بود
کلید واژگان: ضریب توان، دینامیک سیالات محاسباتی، توربین بادی داریوس، شبکه عصبی مصنوعی، الگوریتم ژنتیکWith the advent of smart materials in recent years, the aviation industry and airfoils have undergone many changes. Research into the use of smart materials in aircraft wings to increase their performance and then the use of smart materials in wind turbine airfoils has begun. In this study, computational fluid dynamics and URANS equations for a three-bladed Darrieus wind turbine equipped with a morphing airfoil was used to determine the optimum blade cross-section. 250 airfoils were generated by random control points, in Gambit software, they were unstructured and generated as a sliding mesh then they were simulated in 2D Ansys using the k-ω SST turbulence model and PISO algorithm. Control points and power coefficient were used for Artificial Neural Network (ANN)training and the Genetic Algorithm(GA) was used to optimize the power coefficient. In this study, the base airfoil is NACA0015. The results of that have been very effective. For the first case (Determine the optimal cross-section of the turbine at a full round) the power coefficient of Darrieus wind turbine with the optimal cross-section increased by 42%, and the blade section(airfoil) was also drawn. For the second case(Determine the optimal cross-section in each of the four zones of the rotor), the most efficient sections(airfoils) in four-zone were obtained, increasing the turbine power coefficient by 60% was the result of this optimization.
Keywords: Power Coefficient, Computational Fluid Dynamic, Darrieus Wind Turbine, ANN, GA -
تغییر فاز بخار به مایع که کاربردهای گسترده ای در صنعت دارد را چگالش می نامند. چگالش قطره ای بر روی سطوح آبگریز و ابر آبگریز اتفاق می افتد. تحقیقات نشان داده اند که ایجاد ساختار های میکرو-نانو ساختار بر روی سطوح،یکی از زوش های ابگریز نمودن آنها می باشد.اختلاف دمای بخار و مایع میتواند منجر به ایجاد گرادیان دما و بالطبع همرفت مارانگونی شود. در این مقاله به شبیه سازی تک قطره بر روی سطوح با ساختار های کیسی، ونزل در چهار حالت مختلف و سطح صاف پرداخته شده است و نرخ انتقال حرارت برای تک قطره چگالیده شده بر روی این سطوح و در دو حالت با و بدون مارانگونی محاسبه گشته است. نتایج نشان داده اند نرخ انتقال حرارت عبوری از قطره ساکن در حالتی که اثر مارانگونی احتساب گشته است از حالتی که از آن صرف نظر شده است، بیشتر می باشد.میزان انتقال حرارت عبوری از قطره در حالت بدون مارانگونی بر روی سطح با زبری کیسی از سطح صاف350% و سطح صاف نیز از سطح با زبری ونزل128% بیشتر است.در حالت با مارانگونی نحوه تغییر نرخ انتقال حرارت از این سطوح متاثر از زاویه تماس قطره بر روی سطح می باشد.
کلید واژگان: چگالش قطره ای، دینامیک سیالات محاسباتی، ساختار ونزل، ساختار کیسی، زاویه تماسThe change of the vapor to liquid phase, which is widely used in industry, is called condensation. Dropwise condensation(DWC) occurs on the hydrophobic and the super hydrophobic surfaces. Research has shown that the creation of micro-nanostructured structures on surfaces is one of the ways . The difference in temperature between the vapor and the liquid can lead to temperature gradient, which causes the convection Marangoni convection. This paper simulates a single drop on surfaces with Cassie and Wenzel structures in four different modes and a smooth surface, and the heat transfer rate for a single drop condensed on these surfaces is calculated in two modes with and without marangoni.The results show that the rate of heat transfer passing through the static drop in the case of the presence of marangoni is higher than the non-Maragoni mode. The rate of heat transfer from the droplet in the without-Marangoni state on the surface with the roughness of the cassie from the smooth surface is 350% ,And the smooth surface is 128% higher than the surface with the roughness of the the Wenzel. In the case of Marangoni, the mode of change of heat transfer rate from these surfaces is affected by the contact angle of the drop on the surface.
Keywords: Dropwise Condensation, computational fluid dynamic, Wenzel Structure, Cassie Structure, Contact Angle -
تبدیل انرژی شیمیایی به انرژی الکتریسیته به واسطه دستگاه الکترو-شیمیایی پیل سوختی انجام می گیرد. به وسیله دمای عملیاتی بالای پیل سوختی اکسید جامد (بین 700 و 1000)، تنش گرمایی به وجود آمده که عامل اصلی شروع و انتشار ترک است. این پدیده ممکن است باعث نشت گاز، ناپایداری ساختار و توقف عملیات SOFC قبل از عمر مفید آن باشد. هدف این پژوهش ارایه روشی برای پیش بینی شروع ترک ها در یک SOFC صفحه ای متخلخل ناهمسانگرد است. دما و توزیع تنش محاسبه می شود. در این روش از ضریب شدت تنش و انتگرال-جی مواد برای پیش بینی شروع ترک در داخل آند و کاتد متخلخل استفاده شده است. نتایج تحقیق نشان داد که بالاترین تنش گرمایی در گوشه های بالای کاتد و در گوشه های پایین آند رخ می دهد. به علاوه، ضخامت الکترود کاتد در سمت چپ، به میزان %5/1 افزایش می یابد. درنهایت، شروع ترک در سمت چپ بین گوشه های بالایی و پایینی کاتد رخ می دهد.کلید واژگان: پیل سوختی اکسید جامد، دینامیک سیالات محاسباتی، المان های محدود، تنش گرمایی، ضریب شدت تنش، شروع ترکA fuel cell is an electro-chemical tool capable of converting chemical energy into electrical energy. The high operating temperature of the solid oxide fuel cell (SOFC), (between 700oC to 1000oC), causes thermal stress which is the origin of crack initiation and propagation. Thermal stress causes gas escape, structure variability and SOFC operation cessation before its lifetime. The purpose of the current paper is to present a method that predicts the thermal stress distribution and forecasts the beginning of fissure or crack occurrences in an anisotropic porous electrode of the planar SOFC. The governing coupled non-linear differential equations of heat transfer, fluid flow, mass transfer, mass continuity, and momentum are solved numerically. A code based on computational fluid dynamics (CFD), computational structural mechanics and finite element method (FEM) is developed and utilized. The results show that the highest thermal stress occurs at the lower corners of anode and the upper corners of cathode. The cathode’s thickness at the left side increases by 1.5% and the concentrated temperature and thus the fissure occurs between the top and bottom left corners of the cathode.Keywords: solid oxide fuel cell, computational fluid dynamic, Finite Element, thermal stress, Stress Intensity Factor, crack initiation
-
استفاده از سوخت های متنوع در موتورهای دیزل باعث شده تا این دسته از موتورها به صورت دوگانه سوز مورداستفاده قرار گیرند. کاهش آلایندگی و مصرف سوخت کمتر از ویژگی های موتور دیزل دوگانه سوزی است که از ترکیب گاز- دیزل بهره می گیرند. در این تحقیق سعی شد تا به کمک روش دینامیک سیالات محاسباتی اثر تغییرات زمان پاشش و درصد اختلاط سوخت در موتور دیزل دوگانه سوز OM355 EU2 در دو سرعت مختلف بررسی گردد. اثرات متقابل متغیرهای ورودی و عدم ارایه مدل یکنواخت برای پارامترهای عملکردی موتور موجب شد مدل سازی شبکه عصبی مصنوعی جهت مطالعه اثرات توام متغیرهای زمان پاشش و درصد اختلاط سوخت ها استفاده شود. بهینه سازی به کمک الگوریتم ژنتیک چندهدفه نیز برای تعیین سطوح بهینه متغیرها بکار گرفته شد. نتایج این تحقیق نشان داد که شبکه عصبی مصنوعی RBF پیش بینی های قابل قبولی (RMSE = 0.01 و R2=0.99) از روند تغییرات خروجی ها داشته و به کمک سطح پاسخ شبکه عصبی محدوده بهینه کارکرد موتور در سرعت های مشخصه به دست آمد. همچنین قابلیت های بهینه سازی الگوریتم ژنتیک چندهدفه این امکان را فراهم می سازد تا در سرعت های مختلف سطوح بهینه ای از متغیرهای ورودی و خروجی به دست آید.
کلید واژگان: دینامیک سیالات محاسباتی، شبکه عصبی مصنوعی RBF، الگوریتم ژنتیک چندهدفه، موتور دیزل دوگانه سوز OM55 EU2The dual fuel diesel engine (DDF) is a one of the various IC engine that can used alternative fuel for power generation. Emission and fuel consumption reduction are some of the properties that use a combination of fuel mixture in dual fuel diesel engine (gas-diesel). We focused in this research to study the effect of start of injection (SOI) and blend of fuel variation in OM355 EU2 dual fuel diesel engine at two various speeds with the help of the computational fluids dynamic. The modeling of the artificial neural network (ANN) has been used to study the interaction effects of SOI and blend of fuel on the operational parameters. The non-sorted genetic algorithm (NSGA II) has been used for determining the optimized levels of the variables. The results of this study show that the RBFNN has acceptable predictions about the outputs’ variables (R^2=0.99 ,RMSE=0.01), and the optimized range of the engine function in the specific speeds has been attained with the help of the responding surface of the neural network. Besides, the optimizing capabilities of the NSGA II have provided the optimized levels of the input and output variables in the various speeds.
Keywords: OM355 EU2 Dual fuel diesel engine, computational fluid dynamic, ANN, RBFNN, NSGA II -
امروزه به دلیل کاربرد وسیع موتورهای اشتعال تراکمی در حمل و نقل شهری و نیز لزوم کاهش انتشار آلاینده های خروجی و همچنین مهار مصرف سوخت باعث شده تا استفاده از سوخت های جایگزین در موتور دیزل متداول گردد. به طوریکه یکی از رایج ترین سوختهای جایگزین در موتورهای دیزلی، سوخت های گازی هستند. استفاده از سوخت های جایگزین در موتورهای اشتعال تراکمی مستلزم مطالعه اثرات احتراقی و نحوه تشکیل مخلوط قابل اشتعال در محفظه احتراق این نوع از موتورهاست. متغیرهای پاشش سوخت در کیفیت خود اشتعالی موتورهای دوگانه سوز دیزل/گاز نقش بسزایی را ایفا می کنند. از جمله متغیرهای مرتبط با پاشش سوخت، زمان پاشش و نحوه اختلاط سوخت ها با یکدیگر است. به منظور بررسی اثر زمان پاشش سوخت بر کیفیت احتراق موتورهای دوگانه سوز دیزل/گاز مطالعه ای انجام شد. در این مطالعه به کمک دینامیک سیالات محاسباتی، کیفیت احتراق موتوردوگانه سوز دیزل/گاز ارزیابی شد. برای دستیابی به این هدف، شروع زمان پاشش سوخت در سه حالت 22، 18و 14 درجه قبل از نقطه مکث بالا و درصد اختلاط سوخت متان - ان هپتان با درصدهای10-90، 15-85 و 25-75 مد نظر قرار گرفت. نتایج مطالعات عددی نشان داد که با افزایش زمان پاشش سوخت در موتور دوگانه سوز دیزل/گاز، دمای محفظه احتراق و قله فشار استوانه افزایش یافت و همچنین موجب افزایش طول تاخیر احتراق شد. افزایش قله فشار استوانه موجب افزایش متغیر فشار متوسط نامی شد و متعاقب آن گشتاور خروجی به مقدار 16 درصد افزایش یافت. با بیشتر شدن نرخ جرمی سوخت مایع در مخلوط قابل اشتعال کاهش تاخیر احتراق و افت قله فشار متوسط نامی رخ داد به طوری که گشتاور خروجی به مقدار 17درصد و توان نامی به مقدار 70 درصد کاهش پیدا کرد. با کاهش مدت زمان پاشش سوخت در دورهای تند توان نامی بهبود یافت و نتایج انتشار آلایندگی نیز نشان داد که پیش رسی پاشش موجب افزایش 31 درصد اکسید ازت و 16درصد کربن دی اکسید میشود که نشان دهنده بهبود کیفیت خود اشتعالی در موتور دوگانه سوز دیزل/گاز است.
کلید واژگان: موتور دوگانه سوز دیزل، گاز، متغیرهای عملکردی، دینامیک سیالات محاسباتی، انتشار آلایندگی، زمان پاشش سوختToday, due to the various usage of compression ignition engines in urban transportation, as well as the need to reduce exhaust emissions and control fuel consumption, the use of alternative fuels has become common in diesel engines. Gaseous fuel is one of the most common alternative fuels that can be used in diesel engines. The utilization of alternative fuels in compression ignition engines requires the study of the combustion quality and the type of flammable mixture formation in the combustion chamber. Fuel injection parameter play a significant role in the combustion quality of dual-fuel diesel/gas engines. The injection timing and gaseous fuel mixing percentage together are among the parameters associated with fuel injection. In order to investigate the effect of fuel injection timing on the combustion quality of dual-fuel diesel / gas engines, a study was conducted. In this research, dual-fuel diesel engine is evaluated at two different speeds using the CFD method. To this research, three modes were selected for the fuel IT: BTDC 22°, BTDC 18° and BTDC 14° and three modes for gaseous fuel mixing percentage: 90%-10%, 85%-15% and 75%-25%. The results of numerical studies showed that with increasing fuel injection timing in a diesel / gas dual-fuel engine, the combustion chamber temperature and cylinder pressure peak (PCP) was increased. Increasing the IT also increased the ignition delay, diesel Knocking and distanced the PCP curve from TDC. The IMEP and output torque increased by 16% as the PCP increased. As the mass flow rate of pilot fuel increased, the ignition delay decreased, causing the PCP to fall. The output torque and indicator power decreased by 17% and 70% respectively. Reducing the fuel IT at high speed, the indicator power improved and the emission results showed that the IT advancing increased the NOx and CO2 by 31% and 16% respectively, indicating improved combustion quality in the dual-fuel diesel engine.
Keywords: Dual-fuel diesel engine, Parametric performance, Computational fluid dynamic, Emission of pollutants, Injection timing -
در این مطالعه، اثر اندازه و تعداد حفره های بلوک بتنی سبک وزن بر میزان اتلاف حرارت با استفاده از شبیه سازی عددی مورد بررسی قرار گرفت. به منظور به حداقل رساندن ساده سازی ها، انتقال حرارت ترکیبی جابجایی و تشعشع در حفره های بلوک که در آن ها هوا محبوس شده است بکار گرفته شده است. همچنین در جداره بلوک انتقال حرارت هدایت در نظر گرفته شد. جهت بررسی نحوه قرارگیری حفره های هوا در بلوک، دو قید هندسی در طراحی آن ها اعمال شد که درنتیجه آن 22 پیکربندی ممکن به دست آمد که در آن ها می بایستی حداقل دو و حداکثر نه حفره هوا وجود داشته باشد. به منظور شبیه سازی جریان و انتقال حرارت از روش دینامیک سیالات محاسباتی مبتنی بر روش حجم کنترل محدود استفاده گردید. انتقال حرارت طبیعی سه بعدی به صورت پایا با جریان آرام و تراکم ناپذیر در محفظه ها در نظر گرفته شد. به منظور بررسی و مقایسه عملکرد حرارتی بلوک ها، ضریب هدایت معادل تعیین و نتایج به صورت کانتورهای دما و سرعت ارایه گردیدند و درنهایت بهینه ترین پیکربندی های حفره های هوا از دیدگاه انتقال حرارت معرفی شدند. نتایج نشان می دهد که نحوه قرارگیری و اندازه حفره ها تاثیر قابل توجهی در نرخ انتقال حرارت دارد. همچنین لحاظ کردن انتقال حرارت تشعشعی در حفره ها بسیار با اهمیت می باشد.کلید واژگان: دینامیک سیالات محاسباتی، انتقال حرارت، ضریب هدایت معادل، بلوک سبک بتنی، حفرهIn this study, heat loss through light weight concrete bricks as a function of the size and number of bricks’ cavities was examined via numerical simulation. To minimize the simplifications, conjugate convection and radiation heat transfer within the brick’s cavities filled with air was considered. Also, in the solid of the brick, conduction heat transfer was taken into account. To assess the orientation of the cavities in the brick, two geometrical constrains were employed in the design of the bricks leading to 22 different layouts for which there should be at least two cavities and at most nine enclosures inside a brick. Computational fluid dynamic approaches based on finite volume method were used for the simulations. 3D natural heat transfer with incompressible laminar flow was assumed in the cavities under steady state conditions. To analyze and compare the results in terms of thermal characteristics, the equivalent coefficient of conduction heat transfer was defined. Also, the results were presented with isotherms and velocity contours. Finally, the best configurations of air cavities in the view point of heat transfer were found. The results demonstrate that the configuration and size of the cavities have profound impact on the rate of heat transfer. Also, considering radiation heat transfer in the cavities is rather important.Keywords: computational fluid dynamic, Heat Transfer, equivalent conductivity, lightweight concrete brick, Cavity
-
A fuel cell is an electro-chemical tool capable of converting chemical energy into electricity. High operating temperature of solid oxide fuel cell, between 700oC to 1000oC, causes thermal stress. Thermal stress causes gas escape, structure variability and cease operation of the SOFC before its lifetime.The purpose of the current paper is to present a method that predicts the thermal stress distribution in an anisotropic porous anode of planar SOFC. The coupled governing non-linear differential equations, heat transfer, fluid flow, mass transfer, mass continuity, and momentum are solved numerically. A code based on computational fluid dynamics (CFD), computational structural mechanics and finite element method (FEM) is developed and utilized. The code uses the generated data inside the porous anode in order to detect the temperature and the stress distribution using the Darcy’s law and the Navier-Stokes equations. The numerical results used to govern the areas of high values of stresses were higher than the yield strength of materials. The results show that a highest thermal stress occurs at lower corners of the anode. The concentrated temperature occurs at the middle of the electrolyte-anode whereas the maximum pressure occurs at the middle of the upper and lower section of the anode.
Keywords: Solid oxide fuel cell (SOFC), Computational Fluid Dynamic, Finite elements method, Thermal Stress, anode -
یک مبادله کن گرمایی پوسته- لوله با بافل میله ای با استفاده از نرم افزار HTRI و روش دینامیک سیالات محاسباتی (CFD) شبیه سازی گردید. نتایج به دست آمده از هر دو مدلسازی با داده های تجربی در دسترس مقایسه شد. همچنین اثر فاصله بافل ها و سرعت جریان سیال سرد بر کارایی مبادله کن گرمایی مورد بررسی قرار گرفت. نتایج حاصل نشان داد که کاهش فاصله بافل ها سبب افزایش انتقال گرما در مبادله کن می شود. با این حال، کم کردن بیش از حد فاصله بافل ها می تواند موجب کاهش قابل توجه انتقال گرما گردد. در مقایسه دو روش به کار گرفته شده برای شبیه سازی مبادله کن گرمایی پوسته- لوله با بافل میله ای، می توان گفت استفاده از نرم افزار HTRI آسانتر و سریعتر از روش CFD است و پاسخ های قابل قبولی را به دست می دهد. از طرف دیگر دقت نتایج حاصل از CFD بیشتر است، همچنین جزئیات پروفایل های دما و سرعت سیال در داخل مبادله کن، قابل محاسبه و تجزیه و تحلیل است. نتایج این تحقیق می تواند گامی موثر در توسعه کاربرد این نوع مبادله کن و بازطراحی آن باشد.
کلید واژگان: مبادله کن گرمایی، بافل میله ای، دینامیک سیالات محاسباتی، HTRIA Rod Baffle shell-tube heat exchanger was simulated using HTRI software and Computational Fluid Dynamics (CFD). The results from both models were compared with available experimental data. In addition, the effects of distance of baffles and cold fluid’s velocity on the heat exchanger’s performance were analyzed. The obtained results showed that decreasing the distance between baffles can improve the heat transfer. However, diminution of this distance more than a threshold can lead to a significant reduction in the heat transfer. By comparing the two employed methods for simulation of shell-tube heat exchangers with bar baffles, one can say that the use of HTRI software is easier and faster than the CFD method and achieves acceptable responses, on the other hand, the CFD results are more accurate. Furthermore, the details of the temperature and velocity profiles of the fluid inside the heat exchanger can be calculated and analyzed. The results of this research can be an effective step in developing the application of this type of heat exchanger and redesigning of it.
Keywords: Heat Exchanger, Rod Baffle, Computational Fluid Dynamic, CFD, HTRI -
مبادله کن گرمایی که از چرخه میسوتسنکو بهره می برند علیرغم اینکه از سرمایش تبخیری استفاده می کنند، بازدهی بالایی نیز دارند. لذا در این پژوهش به کمک دینامیک سیالات محاسباتی مبادله کن چرخه میسوتسنکو، شبیه سازی و مورد بررسی قرار گرفته است. در این مقاله تاثیر کلیه پارامترهای عملیاتی از قبیل دما، رطوبت هوای ورودی، سرعت هوا در کانال های تر و خشک و همچنین پارامترهای هندسی مثل طول مبادله کن و ارتفاع کانال ها بر عملکرد مبادله کن مورد بررسی قرار گرفته است. نتایج شبیه سازی نشان داد که این مبادله کن بازدهی مطلوبی در شرایط آب و هوایی مختلف دارد به طوری که بازده حباب تر برای این مبادله کن از 103 تا 115 تغییر می کند اما افزایش رطوبت هوا باعث کاهش بازده و کارایی این مبادله کن می گردد. همچنین نتایج شبیه سازی نشان داد که افزایش سرعت در کانال های خشک باعث کاهش بازده و افزایش سرعت در کانال تر باعث افزایش بازده می شود. نهایتا نسبت بهینه ی دبی جرمی کانال تر به خشک برابر با 3/1 به دست آمده است.
کلید واژگان: مبادله کن گرمایی، چرخه میسوتسنکو، دینامیک سیالات محاسباتیHeat exchangers that use the Maisotsenko Cycle have high efficiency despite using evaporative cooling. Therefore, in this study, based on computational fluid dynamics, the simulation of the Maisotsenko Cycle has been investigated. In this paper, the effects of all operational parameters such as temperature, inlet air huidity, air velocity in wet and dry channels, as well as geometric parameters such as the length of the heat exchanger and the height of the channels on the efficiency of heat exchanger have been studied. The simulation results showed that the heat exchanger has a favorable efficiency in different climatic conditions, so that wet efficiency for this heat exchanger varies from 103 to 115, but the increase in humidity will reduce the efficiency and the effectiveness of the heat exchanger.Also, the results showed that an increase in the velocity inside the dry and wet channels, respectively, decreases and increases the efficiency. Finally, the optimum ratio of the mass flow rate of wet channel to dry one is determined to be 1.3.
Keywords: Heat exchanger, Maisotsenko cycle, computational fluid dynamic
- نتایج بر اساس تاریخ انتشار مرتب شدهاند.
- کلیدواژه مورد نظر شما تنها در فیلد کلیدواژگان مقالات جستجو شدهاست. به منظور حذف نتایج غیر مرتبط، جستجو تنها در مقالات مجلاتی انجام شده که با مجله ماخذ هم موضوع هستند.
- در صورتی که میخواهید جستجو را در همه موضوعات و با شرایط دیگر تکرار کنید به صفحه جستجوی پیشرفته مجلات مراجعه کنید.