بررسی عددی اختلاط دو سیال با لزجت متفاوت در یک میکرو کانال در نسبت‌های منظری مختلف پره به روش شبکه بولتزمن

نوع مقاله : مقاله مستقل

نویسندگان

1 دانشجوی دکتری، دانشکده مکانیک، دانشگاه یزد، یزد، ایران

2 دانشیار، دانشکده مکانیک، دانشگاه یزد، یزد، ایران

3 استادیار، دانشکده مهندسی، دانشگاه میبد، میبد، ایران

چکیده

در مطالعه حاضر اختلاط دو سیال با لزجت‌های مختلف در یک میکرو کانال مجهز به پره نوسانی به روش شبکه بولتزمن با زمان آسایش چندگانه شبیه‌سازی‌شده و اثرات نسبت منظری پره (AR) بر بازده اختلاط بررسی‌شده است. شبیه‌سازی در عدد رینولدز 50، عدد اشمیت 10 انجام‌شده است. نتایج نشان می‌دهد در دامنه نوسان (K) کم، در اعداد استروهال (St) پایین، افزایش AR باعث بهبود بازده اختلاط ولی در سایر اعداد استروهال، ابتدا باعث کاهش و سپس افزایش بازده شده است. در K متوسط، افزایش AR در St کم و میانی، باعث کاهش بازده و در St بالاتر، ابتدا باعث کاهش و سپس افزایش بازده شده است. در K بالا، در St کم و زیاد، افزایش AR موجب بهبود در بازده نمی‌شود و در St میانی، بازده ابتدا کاهش و سپس افزایش می‌یابد. در محدوده مطالعه حاضر در نسبت لگاریتمی لزجت 2 می‌توان مقدار بهینه را در حالت 7/0AR=، 5/0K= و 1St= در نظر گرفت.
در بررسی بازده اختلاط در نسبت‌های مختلف لگاریتمی لزجت (R)، نتایج نشان می‌دهد که در دو سیال با لزجت متفاوت، در اعداد استروهال خیلی کوچک، با افزایش R، بازده افزایش و در بقیه اعداد استروهال، بازده کاهش چشمگیری می‌یابد. به‌طورکلی در دو سیال با لزجت متفاوت، بازده اختلاط در میکرو کانال با پره در نسبت‌های منظری کوچک بهتر است.

کلیدواژه‌ها

موضوعات


[1]  Nguyen NT (1970). Micromixers: Fundamentals, design and fabrication. William Andrew.
[2] Hessel V, Löwe H, Schönfeld F (2005) Micromixers—a review on passive and active mixing principles. Chem Eng Sci 60(8-9): 2479-2501.
[3] Capretto L, Cheng W, Hill M, Zhang X (2011) Micromixing within microfluidic devices. In Microfluidics (pp. 27-68). Springer, Berlin, Heidelberg.
[4] Park, JY, Kim, YD, Kim SR., Han, SY, Maeng JS (2008) Robust design of an active micro-mixer based on the Taguchi method. Sensor Actuat B-Chem 129(2): 790-798.
[5] Koch M, Witt H, Evans AGR, Brunnschweiler A (1999) Improved characterization technique for micromixers. J Micromech Microeng 9(2), 156.
[6] Liu YZ, Kim BJ, Sung HJ (2004) Two-fluid mixing in a microchannel. Int J Heat Fluid Fl 25(6): 986-995.
[7] Liu RH, Stremler MA, Sharp KV, Olsen MG, Santiago JG, Adrian RJ, Beebe DJ (2000) Passive mixing in a three-dimensional serpentine microchannel. J MIcroelectromech S 9(2): 190-197.
[8] An SJ, Kim YD, Heu S, Maeng JS (2006) Numerical study of the mixing characteristics for rotating and oscillating stirrers in a microchannel. J Korean Phys Soc 49(2): 651-659.
[9] Jin SY, Liu YZ, Wang WZ, Cao ZM, Koyama HS (2006) Numerical evaluation of two-fluid mixing in a swirl micro-mixer. J Hydrodyn 18(5): 542-546.
[10] Kim YD, An SJ, Maeng JS (2007) Numerical analysis of the fluid mixing behaviors in a microchannel with a circular cylinder and an oscillating Stirrer. J Korean Phys Soc 50: 505-513.
[11] Celik, B., Akdag, U., Gunes, S., & Beskok, A. (2008). Flow past an oscillating circular cylinder in a channel with an upstream splitter plate. Phys Fluids, 20(10), 103603.
[12] Celik B, Beskok A (2009) Mixing induced by a transversely oscillating circular cylinder in a straight channel. Phys Fluids 21(7), 073601.
[13] Im M, Park JY, Oh YK, Kim YD, Maeng JS, Han SY (2009) Microfluidic analysis of a micro-mixer with an oscillating stirrer.
[14] Shamsoddini R, Sefid M, Fatehi R (2014) ISPH modelling and analysis of fluid mixing in a microchannel with an oscillating or a rotating stirrer. Eng Appl Comp Fluid 8(2): 289-298.
[15] Ortega-Casanova J (2016) Enhancing mixing at a very low Reynolds number by a heaving square cylinder. J Fluid Struct 65: 1-20.
[16] Ortega-Casanova J (2017) CFD study on mixing enhancement in a channel at a low Reynolds number by pitching a square cylinder. Comput Fluids 145: 141-152.
[17] Ghanbari S, Sefid M, Shamsoddini R (2016) Numerical Analysis of two-fluid mixing with various Density and Viscosity in a microchannel with forced oscillating stirrer. Modares Mechanical Engineering 16(8): 109-119.
[18] Khozeymeh-Nezhad H, Niazmand H (2017) A numerical analysis of an active micromixer with the oscillating stirrer at the different aspect ratios by LBM. Modares Mechanical Engineering 17(9): 417-426.
[19] Khozeymeh-Nezhad H, Niazmand H (2018) A double MRT-LBM for simulation of mixing in an active micromixer with rotationally oscillating stirrer in high Peclet number flows. Int J Heat Mass Tran 122: 913-921.
[20] Talon L, Meiburg E (2011) Plane Poiseuille flow of miscible layers with different viscosities: instabilities in the Stokes flow regime. J Fluid Mech 686: 484-506.
[21] Lallemand P, Luo LS (2003) Theory of the lattice Boltzmann method: Acoustic and thermal properties in two and three dimensions. Phys Rev E 68(3): 036706.
[22] Guo Z, Shu C (2013) Lattice Boltzmann method and its applications in engineering (Vol. 3). World Scientific.
[23] Li L, Mei R, Klausner JF (2017) Lattice Boltzmann models for the convection-diffusion equation: D2Q5 vs D2Q9. Int J Heat Mass Tran 108: 41-62.
[24] Zou Q, He X (1997) On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys Fluids 9(6): 1591-1598.
[25] Guo Z, Zheng C, Shi B (2002) An extrapolation method for boundary conditions in lattice Boltzmann method. Phys Fluids 14(6): 2007-2010.
[26] Li L, Mei R, Klausner JF (2013) Boundary conditions for thermal lattice Boltzmann equation method. J Comput Phys 237: 366-395.
[27] Lallemand P, Luo LS (2003) Lattice Boltzmann method for moving boundaries. J Comput Phys 184(2): 406-421.