جریان آرام
در نشریات گروه مکانیک-
پاسخ دینامیکی یک ایرفویل جریان آرام طبیعی که تحت نوسانات هماهنگ پیچشی کم دامنه قرار می گیرد، با استفاده از شبیه سازی گردابه های بزرگ با عدد رینولدز مبتنی بر وتر 750.000 انجام شده است. در طی چرخه های پیچشی، تغییرات زیاد در موقعیت گذرش و همچنین جدایش در لبه ی فرار مشاهده می شود. این منجر به پاسخ غیرخطی نیروهای آیرودینامیکی می شود. اگرچه ماهیت جریان بسیار غیرخطی است، تکامل لایه مرزی روی سطح ایرفویل را می توان با استفاده از یک مفهوم ساده ی تاخیر فاز مدل سازی کرد که تکامل شبه پایای لایه ی مرزی را پیشنهاد می دهد. بر اساس این فرض تاخیر فاز، یک مدل مهندسی ساده ی جدید توسعه داده شده که به خوبی با داده های تجربی انداز ه گیری شده مطابقت دارد. با کمک این مدل، جمله ی شبه پایا به عنوان منبع اولیه ی رفتار غیرخطی در نیروهای آیرودینامیکی ناپایا شناسایی می شود که می تواند از مشخصه های ایرفویل استاتیک استخراج شود. قدرت این جمله با تغییر فرکانس کاهش یافته بدون تغییر باقی می ماند، اما قدرت عبارت نوسانی که اثرات ناپایا را مدل می کند، با افزایش فرکانس افزایش می یابد.کلید واژگان: جریان آرام، شبیه سازی گردابه های بزرگ، نوسان پیچشی، پاسخ ناپایاThe dynamic response of a natural laminar flow airfoil undergoing harmonic small-amplitude pitching oscillations is investigated using Large-eddy simulations with a chord-based Reynolds number of Re = 750;000. Throughout the pitch cycles, large changes in the transition point are seen as well as trailing-edge separation. This leads to a nonlinear response of the aerodynamic forces. Although the nature of the flow is highly nonlinear, the development of the boundary layer over the airfoil surface can be modeled using a simple phase-lag concept which suggests a quasi-steady evolution of the boundary layer. Based on this phase-lag assumption, a simple new empirical model is developed which agrees very well with the measured experimental data. With the aid of this model, the primary source of non-linearities in the unsteady aerodynamic forces is identified to be the quasi-steady term, which can be evaluated from the static airfoil characteristics. The strength of this remains unchanged with the variation of reduced frequency, but the harmonic term's strength, which models the unsteady effects, increases with increasing frequency.Keywords: Laminar Flow, Large Eddy Simulation, Pitching Oscillation, Unsteady Response
-
هدف از انجام پژوهش حاضر، بررسی تجربی مشخصه های انتقال حرارت نانوسیالات پایهروغن درون میکروکانال های موجی با چینش سری و موازی می باشد. آزمایش ها بر روی سیال پایه روغن حاوی نانوذرات فلزی و غیرفلزی TiO2و SiO2 در غلطت های حجمی 0.05 و 0.1 درصد در سه دبی مختلف 0.5 ، 1 و 1.5 لیتر بر دقیقه و چهار دمای ورودی 40، 45، 50، 55 درجه سانتی گراد انجام گرفت. نتایج نشان دهنده افزایش عدد ناسلت سیال پایه تا 41.8 درصد در چینش سری و کاهش دمای سطح در چینش سری نسبت به موازی می باشد. همچنین نانوسیالات TiO2 و SiO2 با غلظت حجمی 0.1 درصد به ترتیب بیشترین میزان افزایش انتقال حرارت را تا 56 و 52.7 درصد در چینش موازی و 45.8 و 42 درصد در چینش سری نسبت به سیال پایه داشته اند. در نهایت، افت فشار مقطع تست در چینش سری 82.1 درصد بیشتر از چینش موازی ثبت گردید.
کلید واژگان: افت فشار، انتقال حرارت جابه جایی، جریان آرام، عدد ناسلت، میکروکانال، نانوسیال پایه روغنThe development of microchannel manufacturing technology has led to a growing interest in using them as heat exchangers. Microchannels are used to control the temperature of equipment and components that generate a high amount of heat flux. Heat transfer can be further increased by dispersing particles (nano-sized particles) with a low volume fraction into the base fluid (nanofluid). A nanofluid can change the thermophysical properties of a base fluid and improve its thermal performance. The purpose of this study was to examine the heat transfer characteristics of oil-based nanofluid within wavy microchannels in series and parallel arrangements of microchannels. In order to examine the performance of each microchannel separately and to make it easier to draw their diagrams, they are named 1 and 2. Experiments were performed on TiO2 and SiO2 oil-based nanofluids in volume fractions of 0.05 and 0.1, flow rates of 0.5, 1.0, and 1.5 lit/min, and inlet temperatures of 40°C, 45°C, 50°C, 55°C. The results show an increase in the Nusselt number of the base fluid up to 41.8% in the series arrangement and also a decrease in the surface temperature in the series arrangement compared to the parallel arrangement. Also, TiO2 and SiO2 nanofluids with a volume fraction of 0.1 caused the highest increase in heat transfer, up to 56% and 52.7% in parallel arrangement and up to 45.8% and 42% in series arrangement compared to the base fluid, respectively. The pressure drop of the test section in series arrangement was up to 82.1% higher than parallel.
Keywords: Microchannel, Oil-Based Nanofluid, Convective Heat Transfer, Nusselt Number, Laminar Flow, Pressure Drop -
ارتعاش القاشده توسط جریان از عوامل موثر تخریب مکانیکی سازه هایی است که در معرض جریان سیال قرار دارند. در این مطالعه با روش گردابه تصادفی-المان مرزی، جریان آرام (رینولدز200) دوبعدی، حول چهار سیلندر استوانه ای یک و دو درجه آزادی در چیدمان مستطیلی، شبیه سازی و ضرایب نیروهای هیدرودینامیکی، خطوط جریان و جابجایی سیلندرها رسم شده است. میدان چرخشی به تعدادی گردابه عددی منفصل و تغییرات آن در هر گام زمانی، با تعقیب گردابه ها در دیدگاه لاگرانژی، تحت اثر دو مکانیزم جابجایی و نفوذ، همچنین ارضای شرط مرزی سرعت صفر روی جدار هندسه، با خلق گردابه ی صفحه ای صورت می گیرد. سیلندر متحرک با سیستم جرم، فنر و دمپر مدل سازی شده است. نتایج نشان داد تغییرات متوسط ضریب پسای سیلندر یک و دو درجه آزادی نسبت به سیلندر ساکن به ترتیب 0/84 و 0/97 برابر می باشد. دامنه ی ارتعاش سیلندرهای عقبی از سیلندرهای جلویی کمتر بود. دامنه ی عرضی نوسان، سه برابر دامنه طولی مشاهده شد. بیشینه دامنه ی عرضی ارتعاش سیلندرهای یک درجه آزادی 1/51 برابر سیلندرهای دو درجه آزادی بدست آمد. مقایسه حل جریان حول تک سیلندر با بکارگیری روش المان مرزی (و عدم احتیاج به نگاشت همدیس یا درنظر گرفتن تصاویر گردابه ها) با نرم افزار انسیس فلوئنت % 25 کاهش زمانی و % 2/3 افزایش دقت محاسبات را نشان داد.
کلید واژگان: ارتعاش القاشده توسط جریان، گردابه تصادفی، المان مرزی، چهار سیلندر استوانه ای، جریان آرامFlow-induced vibration is the effective factor in mechanical destruction of structures that are exposed to fluid flow. In this study, random vortex- boundary element methods (RVM-BEM), is used to simulate two-dimensional laminar fluid flow around four one/two-degrees-of-freedom cylindrical cylinders in a rectangular arrangemen. Hydrodynamic force coefficients, streamlines and cylinder displacements were plotted. The vorticity distribution is separated into blob-vortexes and its changes are studied by tracking these vortexes in the Lagrangian approach by considering two mechanisms of convection and diffusion in each time step. Satisfying no-slip boundary condition, vortex sheets were created on boundary. The cylinder vibrations were modeled as a system of mass, spring, and damper. The results showed that for 1 and 2DoF compared to the stationary cylinder, the average drag coefficient changes are 0.84 and 0.97, respectively. The rear cylinders vibrations amplitude were less than in the front cylinders. The y-amplitude was three times larger than x-amplitude. The maximum x-amplitude vibration of 1DoF cylinders was 1.51 times larger than 2DoF ones. Solving flow over 2DoF single cylinder by BEM(no need for homodis mapping or considering the vortexe images) with a similar solution in Ansys-Fluent software, showed 25% reduction in runtime and 2.3% increase in calculations accuracy.
Keywords: Flow-induced Vibration, Random vortex, boundary element, four cylindrical cylinders, laminar flow -
نشریه یافته های نوین کاربردی و محاسباتی در سیستم های مکانیکی، سال سوم شماره 3 (پاییز 1402)، صص 35 -43
خنک کاری به عنوان یکی از مهمترین چالش های موجود در صرفه جویی انرژی و افزایش بهره وری بسیاری از صنایع مطرح می-باشد. اولین مانع جدی در فشرده سازی و کارآمد کردن دستگاه های انتقال حرارت، خواص ضعیف سیالات انتقال حرارت است. راندمان انرژی تجهیزات انتقال حرارت به تغییرات فلاکس حرارتی ایجاد شده در آن ها بستگی دارد. استفاده از سیالات متداولی همچون آب و اتیلن گلیکول به عنوان سیالات خنک کننده در سیستم خنک کاری دستگاه های مختلف انتقال حرارت جوابگوی دفع شار حرارتی بسیار بالا (ده ها مگا وات بر متر مربع) نیست. امروزه نانو سیالات به عنوان سیالات جدید انتقال حرارت و با اضافه کردن نانو ذرات به سیالاتی که انتقال حرارت پایینی دارند و با ایجاد تغییرات در چگالی، گرمای ویژه و ویسکوزیته این سیالات، به منظور افزایش هدایت حرارتی و بهبود عملکرد انتقال حرارت ساخته می شوند. در این پژوهش نتایج تحقیقات آزمایشگاهی در مورد انتقال حرارت در نانوسیالات تحت جریان های آرام، مغشوش و همچنین جوشش استخری بحث و بررسی می شود.
کلید واژگان: نانو سیال، ضریب هدایت حرارتی، جریان آرام، جریان مغشوش، جوشش استخریJournal of New Applied and Computational Findings in Mechanical Systems, Volume:3 Issue: 3, 2023, PP 35 -43Cooling is one of the most important challenges in saving energy and increasing the productivity of many industries. The first serious obstacle in compressing and making heat transfer devices more efficient is the weak properties of heat transfer fluids. The energy efficiency of heat transfer equipment depends on the heat flux changes created in them. The use of common fluids such as water and ethylene glycol as cooling fluids in the cooling system of various heat transfer devices is not the answer to the removal of very high heat flux (tens of megawatts per square meter). Today, nanofluids are made as new heat transfer fluids by adding nanoparticles to fluids that have low heat transfer and by creating changes in the density, specific heat and viscosity of these fluids in order to increase thermal conductivity and improve heat transfer performance. In this research, the results of laboratory research on heat transfer in nanofluids under laminar, turbulent flows and pool boiling are discussed and investigated.
Keywords: Nanofluid, Thermal conductivity coefficient, Laminar flow, Turbulent flow, Pool Boiling -
حل بسیاری از مسایل مهم صنعتی نیاز به دانستن مقادیر ضریب انتقال حرارت جابجایی یک یا چند جریان سیال در تجهیزات، سیستم ها یا لوله های مختلف دارند. در مطالعه حاضر، یک مدل عددی برای شبیه سازی جریان سیال تراکم پذیر در قسمت ورودی یک لوله داغ با زوایای مختلف نسبت به افق ایجاد شده است. در این ناحیه، لایه های مرزی هیدرودینامیکی و حرارتی جریان سیال در حال توسعه هستند. با توجه به آشفتگی جریان سیال در اثر برهمکنش جریان حرارت و سیال در داخل این لوله، از مدل آشفته سه بعدی برای این شبیه سازی استفاده شد. برای این منظور، معادلات پیوستگی، ناویر-استوکس تراکم پذیر، مدل استرس رینولدز، و معادلات انرژی آشفته و تراکم پذیر بصورت همزمان حل شده اند. سپس با استفاده از مجموعه ای از اجراهای عددی بوسیله مفاهیم طراحی آزمایش و روش های بهینه سازی، یک فرمول پیش بینی برای عدد ناسلت برای این جریان ها به دست آمده است. در نهایت توانایی این فرمول با استفاده از مجموعه ای از داده های آزمایشگاهی مورد بررسی قرار گرفته است.
کلید واژگان: عدد ناسلت، جریان تراکم پذیر، جریان آرام، منطقه ورودی، لوله داغSolving many important industrial problems requires knowing the values of the heat transfer coefficient of passing of one or more fluid streams in different equipment, systems or pipes. In the present study, a numerical model has been developed to simulate compressible fluid flow at the inlet of a hot pipe with different angles to the horizon. In this zone, the hydrodynamic and thermal boundary layers of the fluid flow are developing. Due to the turbulence of the fluid flow due to the interaction of heat and fluid flow inside this tube, a three-dimensional turbulent model was used for this simulation. For this purpose, continuity and compressible Navier-Stokes equations, Reynolds stress model, and turbulent and compressible energy equation have been solved simultaneously. Then, using a set of numerical runs by the concepts of experimental design and optimization methods, a predictive formula for the Nusselt number for these flows has been obtained. Finally, the ability of this formula has been investigated using a set of laboratory data.
Keywords: Nusselt number, Compressible flow, Laminar flow, Entrance zone, Hot tube -
در این مقاله به بررسی انتقال حرارت جابجایی نانوسیال آب/ Fe3O4 در کانال با سطح مقطع مربعی با ابعاد cm1⨯cm1⨯cm80 تحت تاثیر شار حرارتی یکنواخت در جریان آرام فروسیال در حضور میدان مغناطیسی پرداخته می شود. ابتدا به تولید فروسیال با درصدهای حجمی %0/5 و %1 ، بررسی کیفیت فروسیال و روش تولیدی پرداخته می شود. نتایج آزمون های پتانسیل زتا و اشباع مغناطیسی فروسیال نشان از کیفیت و پایداری فروسیال تولیدی دارد. خواص ترمو فیزیکی فروسیال ساخته شده با روابط تجربی موجود مقایسه و ارزیابی می شوند. انتقال حرارت جابجایی فروسیال تولیدی تحت تاثیر شارهای حرارتی کل 546 – 134 وات، بدون حضور میدان مغناطیسی خارجی بررسی و سپس تاثیر میدان مغناطیسی خارجی بر انتقال حرارت فروسیال با درصد حجمی%0/5، تحت تاثیر شار حرارتی کل 1258/2 وات مطالعه می شود. میزان افزایش ضریب انتقال حرارت جابجایی نسبت به آب خالص، بدون وجود میدان خارجی، برای فروسیال با درصد حجمی %1، تحت شار حرارتی کل 134 وات، معادل %30، تحت شار حرارتی کل 545 وات، معادل %48 و تحت شار حرارتی کل 321/3 وات، معادل %38 می باشند. در شار حرارتی 1258/2 وات و برای فروسیال با غلظت %0/5، ضریب انتقال حرارت جابجایی در حضور میدان مغناطیسی خارجی نسبت به حالت عدم وجود میدان %3/16 رشد دارد.
کلید واژگان: میدان مغناطیسی، جریان آرام، انتقال حرارت، عدد ناسلت، فروسیالThis study aimed to investigate the heat transfer of water/Fe3O4 nanofluid in a square cross-sectional channel with dimensions of 80 cm ⨯1 cm ⨯1 cm under the influence of a uniform heat flux perpendicular to the laminar flow of ferrofluid in the presence of a magnetic field. Firstly, the production of ferrofluid with concentrations of 0.5 vol.% and 1vol.%, their quality, and the quality of the production method was investigated. The results of zeta potential and vibrating-sample magnetometer tests show the good quality and stability of the produced ferrofluid. The thermophysical properties of the made ferrofluid are compared and evaluated with existing experimental correlations. The heat transfer of the produced ferrofluids under the influence of heat fluxes of 134-546 Watts is investigated in the absence of an external magnetic field. Then, the effect of the external magnetic field on the heat transfer at 0.5 vol.%, under the influence of a total heat flux of 1258.2 Watts is investigated. The magnitude of increase of heat transfer coefficient compared to pure water, without external field, for ferrofluid with 1 vol.%, under total heat fluxes of 134, 545, and 321.3 Watts, are 30%, 48%, and 38% respectively. At a heat flux of 1258.2 Watts, the heat transfer coefficient in the presence of an external magnetic field increases by 3.16% at 0.5 vol.% compared to the absence of a magnetic field.
Keywords: Magnetic Field, Laminar Flow, Heat Transfer, Nusselt Number, Ferrofluid -
در بسیاری از جریان های سیالات بخش عمده ای از میدان جریان دارای پیچیدگیهای زیادی برای مدل سازی نبوده و می توان از شبکه بندی درشت و نیز معادلات حاکم ساده تر جهت تحلیل عددی استفاده نمود. به منظور کاهش زمان محاسباتی و حافظه کامپیوتری، با بلوک بندی میدان جریان از شبکه بندی مجزای با ابعاد و اندازه مختلف در هر بلوک استفاده شده است. هندسه موردنظر در این تحقیق یک پله عقب گرد دو بعدی آرام است که از سه بلوک بندی مختلف در تحلیل عددی آن استفاده شده است. حافظه مصرفی کامپیوتر و زمان محاسباتی به همراه روند همگرایی و نیز دقت نتایج عددی حاصل شده از سه نوع بلوک بندی مختلف به کار گرفته شده بررسی شد. نتایج نشان داد که در صورت انتخاب ترکیب مناسبی از بلوک بندی و شبکه بندی هر بلوک با صرف زمان و حافظه کمتر می توان پیش بینی هایی در حد دقت نزدیک به نتایج تجربی و سایر روش ها بدست آورد.
کلید واژگان: روش چندبلوکی، جریان آرام، حافظه کامپیوتری، زمان محاسباتیA major part of the flow field has no complicated behavior in many laminar & turbulent flows which can use large mesh and simplified governing equation for numerical simulation. In order to reduce memory and CPU time, the flow field was decomposed to several blocks with different size and dimension in each block. A two dimensional backward facing step was considered in this research that used three different type of blocking in numerical modeling. Computer memory and CPU time consumption in addition to numerical result of main fluid parameters. Consideration showed that, suitable combination of used blocks led to the results with the same accuracy as model was employed only one block and experimental results, in addition to the reduction of memory and CPU time consumption.
Keywords: Multi-Bolck Method, Laminar Flow, Computer memory, CPU time -
شبیه سازی عددی دوبعدی شعله اسپری تشکیل شده در یک پیکربندی جریان متقابل آرام، به منظور بررسی رژیم های شعله و تاثیر آن بر تغییرات گونه ها و رادیکال های مهم، انجام شده است. از اتانول به عنوان یک سوخت اسپری مایع استفاده شده و مکانیسم اسکلتی اتانول/هوا با 40 گونه و 180 واکنش مقدماتی برای واکنش های احتراق به کار گرفته شده است. قطرات سوخت با قطر یکسان و به صورت تصادفی با استفاده از یک کد UDF در ورودی هوا تزریق می شوند و حرکت قطرات براساس دیدگاه لاگرانژی محاسبه شده است. با توجه به نتایج در نسبت های هم ارزی بالا، نرخ های کرنش بالا و قطرهای بزرگ قطرات، رژیم غالب شعله به صورت غیرپیش آمیخته است. همچنین، در نواحی مرکز شعله، به دلیل صفرشدن مقدار کسر جرمی اکسیدکننده و کاهش مقدار کسر جرمی رادیکال OH در این ناحیه، واکنش های احتراقی فرونشانده شده و همچنین مقدار عدد شاخص شعله صفر می شود که نشان دهنده احتراق گروه داخلی قطرات در این ناحیه است.
کلید واژگان: پیکربندی جریان متقابل، احتراق اسپری، جریان آرام، رژیم پیش آمیخته، رژیم غیرپیش آمیختهThe two-dimensional numerical simulation of the spray flame formed in a laminar counterflow configuration has been performed to investigate flame regimes and its effect on changes in important species and radicals. Ethanol is used as a liquid spray fuel and the skeletal mechanism of ethanol/air with 40 species and 180 preliminary reactions is used for combustion reactions. The Sauter mean diameter of fuel droplets and randomly injected into the air inlet are considered using a UDF code, and the motion of the droplets is calculated using the Lagrangian approach. The results show that the predominant flame regime is Non-Premixed in the high equivalences ratio, high strain rates and large droplet diameters. Also, the combustion reactions are suppressed in the center areas of the flame, due to the reduction of the oxidant fraction and the OH mass fraction, and the flame index is zero, which indicates the internal group of droplets in this area.
Keywords: Configuration of Counter-flow, Spray Combustion, Laminar Flow, premixed, non-premixed -
در مقاله حاضر، هندسه بهینه یک کانال چندراهی دوبعدی با یک ورودی و سه خروجی در عدد رینولدز 10 Re=مورد ارزیابی قرار گرفته است. در این ارزیابی، بهینه سازی توپولوژی بر اساس روش تخلخل و به کمک شبیه سازی شبکه بولتزمن صورت گرفته و درنهایت طرح بهینه با استفاده از آنالیز حساسیت یک تابع هدف ارایه شده است. برخلاف تحقیقات پیشین به جای ثابت نگاه داشتن عرض مجرا در خروجی ها، سرعت متوسط در این مقاطع یکسان در نظر گرفته شده است درحالیکه مقایسه نتایج نشان می دهد مقدار اتلافات انرژی در طرح بهینه چندراهی به میزان 04/26 درصد کاهش خواهد یافت که این خود موید مزیت تغییر در عرض مجرا در مقایسه با تغییر در سرعت متوسط جریان است. در این حالت، شرایط هندسه مجرای ورودی شامل عرض و موقعیت مجرا نقش به سزایی در طرح نهایی خواهد داشت که پرداختن به آن هدف اصلی این تحقیق به شمار می رود. نتایج عددی حاکی از آن است که با فرض تغییرات خطی عرض مجرا در خروجی ها، کمترین افت توان جریان زمانی حاصل خواهد شد که عرض ورودی برابر با بزرگترین خروجی باشد. قراردادن ورودی در مقابل عریض ترین خروجی، هندسی نهایی را هموارتر از سایر حالات خواهد کرد اما این حالت از نظر اتلافات کاملا بهینه نیست. برای رسیدن به کمترین اتلافات انرژی، ورودی جریان می بایست هم راستا با فضای میانی دو خروجی با عرض های متوسط و بزرگ قرار گیرد. در این حالت اتلافات به میزان 18/23% کمتر از حالتی است که ورودی در مقابل باریک ترین مجرا قرار گرفته باشد.
کلید واژگان: کانال چندراهی، جریان آرام، بهینه سازی، بهینه سازی توپولوژی، روش لتیس-بولتزمنIn the present paper, the optimal geometry of a two-dimensional multi-channel with one inlet and three outlets is investigated at Re=10. In this study, the topology optimization based on a porosity method using lattice Boltzmann simulation is adopted to find the optimal layout by computing the sensitivity analysis of an objective function. In contrary to previous studies of the channel with constant width of the outlets, the average velocity of the outlets is maintained as the same while the energy dissipation is reduced by 26.04% in comparison of the results, showing the adavntages of changing in width of ducts rather than the change in their average velocities. In this case, the geometry conditions for the inlet duct including width and position of the duct, play an important role in the final design as the main goal of this research. Assuming that the duct width in the outlets changes linearly, the numerical results showed that the equality of inlet width with the largest outlet caused the lowest power loss for the flow. The final geometry was getting much smooth when the inlet was located in front of the outlet with the largest width however it cannot be definitely optimum. In order to obtain the lowest energy loss, the inlet should be parallel to the space between the two outlets with medium and large widths. A reduction of 23.18% of loss was found in this case rather than the inlet was set in front of the duct having the smallest size.
Keywords: Multi-channel, laminar flow, Optimization, Topology Optimization, Lattice-Boltzmann Method -
در این مطالعه، تاثیر یک استوانه یا مفتول کمکی در بالادست جریان سیال حول استوانه به قطر D در یک کانال دوبعدی به ابعاد D50 × D20 و در عدد رینولدز 100 بر الگوهای گرما-شاره استوانه در حالت های ساکن و نوسانی مفتول بررسی شده است. فاصله افقی مفتول از D5/1 تا D8، قطر آن از D1/0 تا D9/0 و همچنین فاصله عمودی مفتول از محور افقی از D1/0 تا D9/0 متغیر است. معادلات حاکم در این مطالعه، معادلات پایستاری با فرض جریان دوبعدی، ناپایا، تراکم ناپذیر، آرام و لزج می باشند. شرط مرزی عدم لغزش در دیواره استوانه ها و شرط دما ثابت در استوانه اصلی و شرط عایق دیواره برای مفتول برقرار است. معادلات حاکم به روش پسرو و با نرم افزار کامسول حل شده اند. نتایج نشان می دهند که وجود مفتول یا استوانه کمکی در بالادست جریان در حالت ساکن باعث کاهش ضریب پسا و عدد نوسلت استوانه اصلی می شود ولی نوسان عرضی آن با قطر D3/0، فاصله D3 و دامنه D25/0، باعث افزایش این مقادیر نسبت به حالت ساکن می گردند که این افزایش در حالت قفل شدگی یعنی 1=F بیش ترین مقدار است.
کلید واژگان: استوانه دایره ای، جریان ارام، مفتول کنترل گر، انتقال حرارتIn this paper, the effects of tripwire in upstream of flow over a circular cylinder in a 2D channel and Re=100 has been investigated. The dimensions of the channel are 20D×50D. The distance between cylinder and wire is considered from 1.5D to 6D, and the relative size of wire is considered from 0.1D to 0.9D. The governing equations are continuity, momentum, and energy for 2D, incompressible, unsteady and vicious conditions. The no-slip condition for both cylinder wall, constant temperature for cylinder, and adiabatic wall for wire have been considered in this paper. The followed equations have been solved by the back-ward differentiate method and the CFD software - COMSOL MULTYPHISIC. The results show that the higher diameter of the wire and horizontal and vertical distances of the wire causes a higher decrease in values of cylinder’s drag and heat transfer. At all, It is found that existing of wire in upstream, decreases the drag and heat transfer of the cylinder but oscillating wire with the amplitude of 0.25D and relative distance and diameters of 3D and.3D causes increasing in discussed parameters relative to the stationary state of them.
Keywords: Circular Cylinder, controller wire, Laminar flow, Heat transfer -
یکی از عوامل موثر در کاهش راندمان دیگ های بخار در طول دوره کارکرد رسوب است. رسوب به صورت یک لایه بر روی سطوح انتقال حرارت دیگ بخار می نشیند و در اثر مقاومت حرارتی که ایجاد می کند راندمان حرارتی دیگ بخار را کاهش می دهد. به همین دلیل تشخیص میزان رسوب و اندازه گیری دقیق آن در لوله های دیگ بخار از اهمیت بسیار بالایی برخوردار است ، با اندازه گیری دقیق رسوب می توان زمان عملیات رسوب زدایی راتعیین کرد. در این مقاله فرض براین است که هیچ اطلاعی در رابطه با نحوه توزیع رسوب در لوله های دیگ بخار موجود نیست، لذا پروفیل رسوب با استفاده از شبیه سازی دما در دیواره داخلی لوله های دیگ بخار و همچنین استفاده از روش گرادیان مزدوج تخمین زده شده است. در این روش میزان خطای حاصل از تخمین تابع پروفیل رسوب چند ضابطه ای محاسبه شد، این بیانگر دقت بالای روش گرادیان مزدوج در تخمین رسوب می باشد. بنابراین می توان از این روش به عنوان یک سیستم هشداری دهنده برای تمیزکاری لوله های دیگ بخار زمانی که میزان رسوب از یک مقدار تعیین شده بیشتر شده باشد، استفاده کرد
کلید واژگان: دیگبخار، رسوب، روش انتقال حرارت معکوس، گرادیان مزدوج، جریان آرامOne of the effective factors in decreasing the efficiency of boilers during the operation period is sediment. The sediment sits as a layer on the boiler's heat transfer surfaces and reduces the boiler's thermal efficiency due to the thermal resistance it creates. For this reason, it is very important to determine the amount of sediment and to accurately measure it in boiler pipes.In this paper, it is assumed that no information is available on the distribution of sediment in boiler pipes, so the sediment profile is estimated using temperature simulation on the inner wall of boiler pipes as well as conjugate gradient method. In this method, the error rate was calculated by estimating the function of multi-criteria sediment profile , which indicates the high accuracy of the conjugate gradient method in sediment estimation. The amount of sediment exceeded a specified value.
-
این پژوهش رفتارهای حرارتی و هیدرودینامیک جریان آرام سیال عامل آب درون لوله خمیده افقی وایپرتکس، در محدوده اعداد گراشهف پایین و در محدوده انتقال حرارت جابه جایی ترکیبی را به روش عددی مورد بررسی قرار داده است. این لوله دارای خم 180درجه، شعاع انحنا بی بعد 6.62=R/D2 و زبری های کروی شکل با ارتفاع بی بعد 0.1=e/D و شار حرارتی ثابت روی دیواره خود است. معادلات پیوستگی، مومنتوم و انرژی با روش حجم محدود به صورت سه بعدی حل شده اند. برای حل مساله از کد دینامیک سیالات محاسباتی نرم افزار انسیس فلوئنت استفاده شده است. نتایج حاصل بیانگر تاثیر عدد گراشهف و نیروهای شناوری ناشی از آن، همچنین اثر تقابلی نیروهای گریز از مرکز و شناوری در کنار زبری سطح، روی پارامترهای هیدرودینامیک و حرارتی از جمله سرعت محوری بی بعد، جریان های ثانویه، دمای بی بعد، ضریب انتقال حرارت جابه جایی و ضریب اصطکاک سطحی است. به گونه ای که در یک رینولدز ثابت با افزایش گراشهف به دلیل تقابل شدیدتر نیروهای شناوری و گریز از مرکز، ضریب انتقال حرارت جابه جایی کاهش می یابد حال آن که در لوله خمیده صاف (بدون زبری) افزایش گراشهف منجر به افزایش ضریب انتقال حرارت جابه جایی می شود. با این وجود لوله خمیده وایپرتکس در محدوده گراشهف و رینولدزهای بسیار پایین، در هر گراشهف و رینولدز یکسان، انتقال حرارت بالاتری نسبت به لوله خمیده صاف داشته است. همچنین نتایج نشان داد که ضریب اصطکاک سطحی در این نوع لوله ها با افزایش گراشهف می تواند تا 3/5برابر بیشتر از مقدار آن در لوله خمیده صاف باشد.کلید واژگان: جریان آرام، لوله خمیده افقی، وایپرتکس، انتقال حرارت جابه جایی ترکیبیThermal and hydrodynamic behavior of a laminar flow of water within a horizontal curved Vipertex tube with mixed convection heat transfer, in the range of low Grashof numbers, has been numerically studied. The curved horizontal Vipertex tube has geometry of 180o, fixed radius of centerline curvature of 2R/D=6.62, roughness height e/D=0.1, and a constant heat flux is exerted on the walls. The three-dimensional governing equations were using a finite volume method. To solve the problem, the computational fluid dynamics of ANSYS Fluent The results reveal that not only Grashof number and the buoyancy forces arising from it, but the mutual effects of the centrifugal and the buoyancy forces affect the thermal and hydrodynamic characteristics such as axial velocity contours, secondary flow vectors, temperature contours, heat transfer coefficient, and skin friction coefficient. So that, for a given Reynolds number, increasing due to more interaction between buoyancy and centrifugal forces, results in the Vipertex tube. Therefore, the buoyancy forces decrease and lead to the lower heat transfer coefficient, but in smooth curved Grashof number leads to the higher heat transfer coefficient. Nevertheless, the Vipertex curved tube in of Grashof and Reynolds, in each Grashof and Reynolds equally, has a higher heat transfer than a smooth curved pipe. The results also indicated that the skin friction coefficient in these types of tubes can be up to 3.5 times higher than that of smooth one with a Grashof increase.Keywords: Laminar Flow, Horizontal Curved Tube, Vipertex, Mixed Convection Heat Transfer
-
در این پژوهش، یک سیستم فتوولتاییک/ حرارتی با متمرکزکننده سهموی- خطی شبیه سازی شده و اثر استفاده از نانوسیال اکسید آلومینیوم/ اتیلن گلیکول: آب 50:50 دارای نانوذرات مختلف پلاکتی، استوانه ای، تیغه ای و آجری شکل بر کارایی این سیستم از دیدگاه انرژی و اگزرژی در دو جریان آرام و مغشوش بررسی شده است. مدل پیشنهاد شده با استفاده از نتایج آزمایشگاهی موجود اعتبارسنجی شده است و تطابق مناسبی بین نتایج مشاهده گردیده است. باتوجه به نتایج، استفاده از نانوسیال با نانوذرات استوانه ای شکل در جریان آرام و نانوذرات آجری شکل در جریان مغشوش منجر به حداکثر بازده سیستم می گردد. به علاوه، استفاده از نانوسیال در جریان آرام نسبت به جریان مغشوش برای بهبود عملکرد سیستم فتوولتاییک/ حرارتی موثرتر است.کلید واژگان: سیستم فتوولتاییک مغشوش حرارتی، شکل نانوذرات، نانوسیال، جریان آرام، جریانIn this study, a linear parabolic trough concentrating photovoltaic/thermal system has been simulated and the effects of using Al2O3/ethylene glycol-water 50:50 nanofluid with different nanoparticle shapes including platelet, cylindrical, blade and brick shapes from energy and exergy standpoints in the laminar and turbulent regimes have been numerically investigated. The proposed model has been validated using existing experimental results where good agreement was observed. The results indicated that using nanoparticles of cylindrical shape in laminar regime and brick-shaped one in turbulent regime lead to the best system performance compared to others. In addition, applying nanofluid in laminar regime is more effective compared to turbulent regimeKeywords: Photovoltaic, Thermal system, Nanoparticles shape, Nanofluid, Laminar flow, Turbulent flow
-
در مطالعه حاضر، اثر جابجایی آشوبناک بر اختلاط درون یک همزن جریان آرام به صورت تجربی و عددی بررسی شده است. همزن شامل دو روتور و یک استاتور استوانه ای با سطح مقطع دایروی است به صورتی که می توان سرعت چرخش روتورها را با زمان کنترل نمود. این نوع همزن در صنایع غذایی و دارویی کاربرد فراوان دارد. عملکرد همزن به صورت تجربی از طریق تزریق رنگ و ثبت تصاویر مورد بررسی قرار گرفت و به صورت عددی نیز اثر جریان آشوبناک بر اختلاط با ردگیری لاگرانژی ذرات، محاسبه میزان کشیدگی المان های سیال و ایجاد مقاطع پوآنکاره به صورت کمی و کیفی نتایج تجربی نشان می دهد با اغتشاش سینوسی اعمال شده به سرعت دورانی روتورها نواحی ضعیف اختلاط در طول زمان از بین رفته و ذرات سیال به خوبی در سطح همزن پخش می شوند. با استفاده از شبیه سازی عددی نیز نشان می دهند هنگامی که سرعت دورانی روتورها متغیر است جریان درون همزن حساس به شرایط اولیه می باشد که این امر یکی از شاخصه های مهم جریان آشوبناک است. همچنین با محاسبه مقدار کشیدگی المان های سیال نشان داده شد که مقدار متوسط نمایی کشیدگی در جریان آشوبناک نسبت به جریان غیر آشوبناک افزایش 2 برابری دارد. در این تحقیق بر اساس مفاهیم علمی و کاربردی تئوری آشوب در مکانیک سیالات، راندمان همزن های دور پایین در اختلاط سیالات بسیار لزج بدون افزایش انرژی مصرفی، افزایش چشمگیری یافته است.کلید واژگان: همزن آشوبناک، اختلاط، سیال بسیار لزج، جریان آرامIn this study, the effect of chaotic advection on laminar mixing is analyzed experimentally and numerically. Mixer includes two circular rotors and a circular stator in the way that rotational speed of each rotor could be controlled by time. This kind of mixer is widely used in the food and pharmaceutical industry. The performance of the mixer was investigated experimentally with dye injection and image recording. Moreover, the effect of chaotic flow on mixing numerically is studied qualitatively and quantitatively with Lagrangian fluid particle tracing, calculating the stretching rate of fluid elements and Poincare sections. The experimental results indicate that as a result of applying sinusoidal perturbation on rotational velocity of rotors, weak mixing zones will disappear with the passage of time, and fluid particles will be distributed uniformly in the surface of the mixer. Numerical simulation shows that the mixer flow is sensitive to initial condition when rotational velocity of rotors is variable, and this is one of important factors in chaotic flow. Further, calculation of stretching rate of fluid elements indicates that average exponential value of fluid elements stretching in the chaotic flow is two times more than non-chaotic one. This research, based on scientific concepts and theoretical application of chaos in fluid mechanics, denotes that the efficiency of low-speed mixers in highly viscous fluid mixing is increased without any increase in the consuming energy.Keywords: Chaotic mixer, mixing, highly viscous fluid, laminar flow
-
در پژوهش حاضر تاثیر حضور مفتول دایروی در مقایسه با استوانه بدون مفتول روی پارامترهای ضریب پسا و انتقال گرما برآورد شده و سپس تاثیر قطر و زاویه قرارگیری مفتول بررسی می شود. معادلات پیوستگی، مومنتوم و انرژی به همراه شرایط مرزی با روش حجم کنترل حل شده است. برای گسسته سازی عبارتهای جابجایی از طرح بالا دست مرتبه دوم و برای پیوند سرعت و فشار الگوریتم PISO به کار رفته است. موقعیت زاویه ای قرارگیری مفتول، گستره 0 تا 180 درجه از نقطه سکون جلویی و نسبت قطر مفتول به قطر استوانه مقادیر 05/0 تا 4/0 در عددهای رینولدز جریان آرام کار شده است. افزودن مفتول تا موقعیت زاویه ای 90o درجه موجب کاهش انتقال گرما نسبت به استوانه ساده می شود اما بعد از این زاویه، انتقال گرما نسبت به حالت ساده افزایش می یابد. همچنین با کاهش عدد رینولدز تاثیر افزودن مفتول بیشتر می گردد.کلید واژگان: استوانه دایروی، مفتول کنترلگر، جریان آرام، معادلات ناویر، استوکس، انتقال گرماIn this work, the effect of a controller wire on the hydro-thermal flow characteristics around a circular cylinder in a channel was studied. The two-dimensional unsteady and incompressible flow with heat transfer was investigated by solving the Navier-Stokes and energy equations using control volume methodology. The second-order upwind scheme is used to discretize convective terms. The discretized equations are solved by PISO algorithm. The angular position for the wire investigated from 0° to 180° from the forward stagnation point and the diameter ratio for wire and the cylinder is set at 0.05 to 0.4, in laminar Reynolds numbers. The surface Nusselt number of cylinder is increased by adding a wire on the cylinder surface between 90ᵒ and 180ᵒfrom the forward stagnation point. In addition, as the Reynolds number decreases, the effects of adding a wire on the main cylinder are becoming stronger.Keywords: Circular cylinder, Controller wire, laminar flow, Navier-Stokes equations, Heat transfer
-
در این پژوهش، انتقال حرارت جابجایی اجباری جریان آرام نانوسیال آب - مس درون کانال با صفحات موازی با منابع حرارتی دما ثابت، به صورت عددی بررسی شده است. هدف، بررسی عددی تاثیر عدد رینولدز، درصد حجمی نانوذرات و تعداد منابع حرارتی روی میدان جریان و نرخ انتقال حرارت می باشد. معادلات حاکم به روش حجم کنترل جبری شده و توسط الگوریتم سیمپل به طور هم زمان حل شدند. نتایج نشان می دهند که افزایش درصد حجمی نانوذرات و افزایش عدد رینولدز، هم چنین تقسیم یک منبع حرارتی به قسمت های کوچکتر و توزیع آنها، نرخ انتقال حرارت را افزایش می دهند.کلید واژگان: جابجایی اجباری، نانوسیال، جریان آرام، کانال، منابع حرارتیIn this study, the laminar forced convection heat transfer of water-coper nanofluid is numerically investigated within a parallel plate channel. Fixed temperature heat sources with the specified sizes and distances are embedded on the walls of the channel. The entry and exit sections of the channel as well as the sections between the heat sources are thermally insulated. The fluid flow with uniform velocity and temperature enters the channel. The channel length is considered large enough, so the flow in the channel output is assumed fully developed. The aim of this research is the numerical investigation of the effects of the Reynolds number, the solid volume fraction and the number of the heat sources on the flow field and heat transfer rate. For this purpose, the governing equations are discretized by finite difference method based on the control volume formulation and are solved using the SIMPLE algorithm. In order to validate the computer program, the results of this study have been compared with the results of the previous numerical studies. This comparison has confirmed the accuracy of the performance of the computer program. The results show that the rate of heat transfer increases by increasing the solid volume fraction and Reynolds number. The results also show that, heat transfer rate increases when the heat sourse is divided into smaller sections and these sections are distributed on the channel wall.Keywords: Forced convection, nanouid, Laminar Flow, channel, discrete heat sources
-
نشریه شبیه سازی و تحلیل تکنولوژی های نوین در مهندسی مکانیک، سال نهم شماره 2 (تابستان 1395)، صص 269 -280با توجه به نیاز افزایش راندمان انتقال حرارت در سیستم های حرارتی بخصوص در سیستم هایی با ابعاد میکرو و نانو وجود ایده استفاده از نانوسیالات در تجهیزات انتقال حرارت برای بهبود عملکرد آن ها بسیار مناسب می باشد.
دراین مقاله جریان و انتقال حرارت جابجایی اجباری نانوسیال آب-مس در یک میکرو لوله در رژیم لغزشی تحت شار حرارتی ثابت در مقادیر رینولدز پایین به صورت عددی موردبررسی قرار می گیرد. شرط مرزی لغزش و پرش دمایی در دیواره برای اولین بار مورداستفاده قرار گرفت و نتایج به صورت پروفیل های سرعت و دما و محاسبه مقادیر ناسلت و افت فشار موضعی در ناحیه ورودی و توسعه یافته کانال ارائه شده است. اثر لغزش و استفاده از نانو سیالات بر انتقال حرارت در میکرو لوله بررسی گردیده است.
نتایج نشان می دهد که عدد ناسلت در ناحیه توسعه یافته با افزایش لغزش در دیواره، افزایش و میزان افت فشار در این ناحیه کاهش می یابد و نتایج در این ناحیه مستقل از عدد رینولدز جریان هست. در ناحیه ورودی افزایش عدد رینولدز جریان نیز باعث افزایش این مقادیر در این ناحیه می-گردد. همچنین در این پژوهش مشاهده گردید افزودن نانو ذرات به سیال باعث افزایش انتقال حرارت جریان گردیده اما این مقدار ناچیز است.کلید واژگان: نانوسیال آب مس، میکرو لوله، جریان لغزشی، پرش دمایی، جریان آرامJournal of Simulation and Analysis of Novel Technologies in Mechanical Engineering, Volume:9 Issue: 2, 2016, PP 269 -280Given the need to increase the efficiency of heat transfer in thermal systems, especially systems using nanofluids in microscale and nanoscale heat transfer equipment ideas to improve their performance is very good.In present study, the flow and heat transfer of Water-Cu nanofluid in micro-tube with slip regime with constant wall heat flux numerically simulated with low Reynolds numbers. Slip velocity and temperature jump boundary conditions are also considered along the microtube walls, for first time. The results are presented as the profiles of temperature and velocity. Nusselt number and pressure drop coefficient calculated in interance and full developed region. The effect of slip and using nano particle considerd.
It is observed that Nusselt number increases with slip velocity coefficient and pressure drop coefficient decreases; att intrance region the Raynolds of flow has effect on Nusselt and pressure drop coefficient,too.
Likewise observed nano particle adding to water has low effect to increases Nusselt number and pressure drop coefficientt.Keywords: Water-Cu nanofluid, Microtube, Slip flow, Temperature jump, Laminar flow -
مساله بررسی عددی انتقال حرارت در کانال با مقطع مثلث متساوی الاضلاع با قطر هیدرولیکی 8 میلیمتر، طول 1 متر با استفاده از نانوسیال آب-اکسیدمس، در 2 قطر 20 و80 نانومتری و کسرحجمی 1، 2 و 4 درصد، ابتدا در حالت شار ثابت در هر 3 وجه کانال و سپس حالت شارثابت بر روی 2 وجه و دمای ثابت بر روی وجه کف(صفحه داغ) می باشد. نانوسیال بصورت تکفازی در نظر گرفته می شود. در تکنولوژی نانو، اولین اثر کاهش اندازه ذرات، افزایش سطح است. افزایش نسبت سطح به حجم نانوذرات باعث می شود که اتم های واقع در سطح، اثر بسیار بیشتری نسبت به اتم های درون حجم ذرات، بر خواص فیزیکی ذرات داشته باشند. همچنین ترکیب نانوذرات با سیال، بواسطه حرکت براونی نانوذرات و اثر خوشه ایآنها، موجب افزایش ضریب هدایت گرمایی و ضریب انتقال حرارت جابجایی سیال می شود. از آنجا که در بحث انتقال حرارت، افزایش ضریب انتقال حرارت جابجایی به تنهایی کافی نبوده و عدد ناسلت نیز باید بررسی شود. در این پایان نامه سعی کردیم هر دو پارامتر را در کنار هم بررسی کنیم. رژیم جریان آرام در نظر گرفته شده است. ضریب انتقال حرارت جابجایی و عدد ناسلت، تاثیر قطر نانوذرات بر میزان انتقال حرارت، تاثیر کسرحجمی نانوذرات بر ماکزیمم سرعت در مقطع عرضی بررسی شده است.
کلید واژگان: ضریب انتقال حرارت جابجایی، نانوسیال، عدد ناسلت، جریان آرام -
در این پژوهش معادلات ناویر-استوکس برای تحلیل حرکت شتابدار کره صلب داخل سیال تراکم ناپذیر به کار رفته اند. ایده حرکت شتاب لحظه ای کره برای اولین بار مطرح شده است. تکانش در یک لحظه کوچک به وجود آمده و سرعت قبل و بعد از تکانش ثابت در نظر گرفته می شود. جریان آرام و ناپایا در رینولدزهای مختلف حل شده است. در یک لحظه عدد رینولدز از یک مقدار ثابت تغییر پیدا کرده و موجب تغییرات اساسی پسا می شود. برای حل سایر حرکتهای شتابدار (مانند حرکت نوسانی) نوع مدل دیگر (مانند زیردریایی های کره مانند) این روش قابل اعمال می باشد.
کلید واژگان: شتاب لحظه ای، تطبیق شبکه، معادلات ناویر، استوکس، جریان آرام، جرم افزوده، ضریب پسا، جریان ناپایا، جریان متقارن محوریIn this study, accelerating motion of a sphere in incompressible fluid has been nvestigated by the Navier-Stokes equations. The motion is set in a small time interval, having the constant velocities out of it.The laminar and unsteady flow has been solved in various Reynolds numbers. In a short moment the Reynolds number is rised which cause large variations in the drag force. This method is applicable to other accelerating motions such as oscillation and also spherical submarines.Keywords: Acceleration, Adaptive grid, Navier, Stokes equations, Laminar flow, Added mass, Drag coefficient -
در این تحقیق، انتقال گرمای همرفت آزاد از روی یک صفحه سینوسی عمودی دما ثابت با حل معادلات ناویر- استوکس، در رژیم جریان آرام و به روش حجم محدود تحلیل شده است. همچنین، مدل جدیدی بر پایهی فشار برای محاسبه شارهای جابه جایی ارائه شده است. نسخه های مرتبه اول و دوم این مدل پیشنهادی قادرند تا شارهای جابه جایی را با دقت خوبی بدون نیاز به اتلاف مصنوعی محاسبه نمایند. نتایج عددی به خوبی با نتایج سایر پژوهشگران مطابقت دارد. در این بررسی، پارامترهای تغییرکننده شامل عدد ریله و نسبت دامنه به طول موج صفحه سینوسی میباشند. نتایج نشان میدهد که فرکانس تغییر ضرایب انتقال گرمای محلی برای صفحه در فضای آزاد برابر فرکانس خود صفحه میباشد. در ضمن، با افزایش عدد ریله ضریب انتقال گرمای میانگین افزایش مییابد.
کلید واژگان: همرفت آزاد، معادلات ناویر، استوکس، فرض بوسینسک، صفحه موج دارسینوسی، روش حجم محدود، جریان آرامIn this research, free convection heat transfer over vertical sinusoidal wall in constant temperature has been performed, using Navier-Stokes equations in laminar flow regime with finite volume method. Also, a new pressurebased model for the convective fluxes is offered, which minimizes the need for artificial dissipation. The first and second-order versions of the suggested model are unable to calculate the convective fluxes with a good accuracy.Numerical results are well adapted with the numerical or experimental results of other researchers. In this study, the variable parameters are Reyleigh number and the ratio of domain length to the sinusoidal wave length. The results show that the local heat transfer coefficient variation frequency for space is equal to the wall frequency. Meanwhile,with increasing the Reyleigh number, the average heat transfer coefficient increases.Keywords: Free Convection, Navier, Stokes Equations, Boussinesq Approximation, Sinusoidal Wavy Wall, Finite
- نتایج بر اساس تاریخ انتشار مرتب شدهاند.
- کلیدواژه مورد نظر شما تنها در فیلد کلیدواژگان مقالات جستجو شدهاست. به منظور حذف نتایج غیر مرتبط، جستجو تنها در مقالات مجلاتی انجام شده که با مجله ماخذ هم موضوع هستند.
- در صورتی که میخواهید جستجو را در همه موضوعات و با شرایط دیگر تکرار کنید به صفحه جستجوی پیشرفته مجلات مراجعه کنید.