عبور از یک سیلندر

View Categories

عبور از یک سیلندر

9 min read

PDF

عبور از یک سیلندر
معرفی
محاسبه جریان سیال در پشت بدنه ای مانند خودرو به دلیل جریان های ناپایدار دشوار است. دنباله پشت چنین بدنه ای از گرداب های نامرتب در هر اندازه تشکیل شده است که کشش زیادی روی بدنه ایجاد می کند. در مقابل، آشفتگی در لایه‌های مرزی نازک در کنار بدنه‌های ساده هواپیما و ماهی تنها اختلالات ضعیفی در جریان ایجاد می‌کند.
یک استثنا در این مورد زمانی اتفاق می‌افتد که یک بدن باریک را در زوایای قائم با جریان آهسته قرار دهید زیرا گرداب‌ها سازماندهی می‌شوند. یک خیابان گرداب فون کارمان با فرکانس قابل پیش بینی ظاهر می شود و شامل ریختن گرداب ها از طرف های متناوب است. نمونه‌های روزمره این پدیده شامل آواز خواندن سیم‌های تلفن و ارتعاش آنتن رادیویی خودرو در جریان هوا است.
از نقطه نظر مهندسی، پیش‌بینی فرکانس ارتعاشات در سرعت‌های مختلف سیال و در نتیجه اجتناب از تشدید نامطلوب بین ارتعاشات سازه‌های جامد و ریزش گرداب مهم است. برای کمک به کاهش چنین اثراتی، مهندسان کارخانه یک مارپیچ را در قسمت بالایی دودکش های بلند قرار می دهند. تغییر شکل حاصل از تداخل سازنده عناصر گردابی که سازه از موقعیت های مختلف می ریزد، جلوگیری می کند.
تعریف مدل
برای نشان دادن اینکه چگونه می‌توانید چنین اثراتی را مطالعه کنید، مدل زیر جریان ناپایدار و تراکم‌ناپذیری را از کنار یک استوانه بلند که در کانالی با زاویه قائم به سیال روبرو قرار گرفته است، بررسی می‌کند. با مشخصات سرعت ورودی متقارن، جریان به نوعی عدم تقارن نیاز دارد تا تولید گرداب را آغاز کند. این را می توان با قرار دادن سیلندر با فاصله کمی از مرکز جریان به دست آورد.
پیش‌بینی زمان شبیه‌سازی لازم برای ظاهر شدن یک الگوی جریان دوره‌ای دشوار است. یک پیش بینی کننده کلیدی عدد رینولدز است که بر اساس قطر سیلندر است. برای مقادیر کم (زیر 100) جریان ثابت است. در این شبیه‌سازی، عدد رینولدز برابر با 100 است که یک خیابان گرداب فون کارمان توسعه‌یافته را نشان می‌دهد، اما جریان هنوز کاملاً آشفته نیست.
فرکانس و دامنه نوسانات ویژگی های پایداری هستند، اما جزئیات جریان به شدت به آشفتگی ها حساس هستند. برای به دست آوردن قدردانی از این حساسیت، می‌توانید تصاویر جریانی را که همزمان گرفته شده‌اند، اما با تفاوت‌های جزئی که توسط تلورانس‌های مختلف برای مرحله زمانی ایجاد می‌شوند، مقایسه کنید. توجه به این نکته مهم است که این حساسیت یک واقعیت فیزیکی است و صرفاً یک مصنوع عددی نیست.
قبل از محاسبه نیروهای متغیر زمان بر روی استوانه، می توانید روش محاسبه را با عدد رینولدز پایین تر با استفاده از حل کننده غیرخطی مستقیم تأیید کنید. این باعث صرفه جویی در زمان می شود زیرا می توانید خطاها و اشتباهات ساده را قبل از شبیه سازی نهایی وابسته به زمان پیدا کرده و تصحیح کنید که به زمان قابل توجهی نیاز دارد.
نیروهای چسبناک روی سیلندر با گرادیان میدان سرعت در سطح سیلندر متناسب است. ارزیابی گرادیان سرعت روی مرز با تمایز مستقیم راه‌حل FEM ممکن است اما چندان دقیق نیست. این تمایز چند جمله ای های مرتبه 1 را هنگامی که عناصر درجه دوم برای میدان سرعت استفاده می شود، تولید می کند. یک روش به مراتب بهتر استفاده از یک جفت عملگر نیروی واکنش برای محاسبه انتگرال نیروهای ویسکوز است که با انتگرال های دقیق مرتبه دوم نیروهای ویسکوز قابل مقایسه است. یک رویکرد جایگزین استفاده از یک جفت متغیر با محدودیت ضعیف برای اعمال شرایط بدون لغزش است. هنگام محاسبه انتگرال نیروهای واکنش یا شار در پس پردازش، ترجیحاً از عملگر نیروی واکنش به جای قیود ضعیف استفاده کنید.
خود نیروهای درگ و لیفت به اندازه ضرایب درگ و لیفت بدون بعد جالب نیستند. اینها فقط به عدد رینولدز و شکل یک جسم بستگی دارند، نه اندازه آن. ضرایب به صورت تعریف شده است
با استفاده از پارامترهای زیر:
F D و F L نیروهای درگ و بالابر هستند
ρ چگالی سیال است
میانگین U میانگین سرعت است
A مساحت پیش بینی شده (محصول ضخامت و قطر سیلندر) است.
نتایج و بحث
شکل 1 الگوی جریان حاصل از هندسه را نشان می دهد.
شکل 1: نمودار آخرین مرحله زمانی به وضوح مسیر فون کارمان را نشان می دهد.
جریان در اطراف یک استوانه یک آزمایش معیار رایج برای الگوریتم‌های CFD است. تیم‌های تحقیقاتی مختلف با استفاده از تکنیک‌های مختلف نقاط قوت خود را روی این مشکل امتحان کرده‌اند. نتایج برخی از این آزمایش‌ها توسط Schäfer و Turek جمع‌آوری شده است .
شکل 2 نشان می دهد که چگونه ضریب بالابر با تشکیل ساختار گرداب فون کارمان، یک تغییر دوره ای ایجاد می کند.
شکل 2: ضریب لیفت، L ، به عنوان تابعی از زمان.
شکل 3: ضریب کشیدن، D ، به عنوان تابعی از زمان.
ارجاع
. _ _ _ ، 1996.
مسیر کتابخانه برنامه: COMSOL_Multiphysics/Fluid_Dynamics/cylinder_flow
دستورالعمل مدلسازی
از منوی File ، New را انتخاب کنید .
جدید
در پنجره جدید ، روی  Model  Wizard کلیک کنید .
مدل جادوگر
1
در پنجره Model  Wizard روی  2D کلیک کنید .
2
در درخت Select  Physics ، Fluid  Flow>Single-Phase  Flow>Laminar  Flow  (spf) را انتخاب کنید .
3
روی افزودن کلیک کنید .
4
 روی مطالعه کلیک کنید .
5
در درخت انتخاب  مطالعه ، General  Studies>Time  Dependent را انتخاب کنید .
6
 روی Done کلیک کنید .
تعاریف جهانی
پارامترهای 1
1
در پنجره Model  Builder ، در قسمت Global  Definitions روی Parameters  1 کلیک کنید .
2
در پنجره تنظیمات برای پارامترها ، بخش پارامترها را پیدا کنید .
3
در جدول تنظیمات زیر را وارد کنید:
 
نام
اصطلاح
ارزش
شرح
U_mean
1[m/s]
1 متر بر ثانیه
میانگین سرعت ورودی
اچ
0.41 [m]
0.41 متر
ارتفاع
دبلیو
2.2[m]
2.2 متر
عرض
آر
0.05 [m]
0.05 متر
شعاع سیلندر
در مرحله بعد، یک ویژگی تابع مرحله هموار ایجاد کنید که برای افزایش سرعت جریان از آن استفاده خواهید کرد.
مرحله 1 (مرحله 1)
1
در نوار ابزار Home ، روی  Functions کلیک کنید و Global>Step را انتخاب کنید .
2
در پنجره تنظیمات برای مرحله ، قسمت پارامترها را پیدا کنید .
3
در قسمت متن مکان ، 0.05 را تایپ کنید .
هندسه 1
مستطیل 1 (r1)
1
در نوار ابزار Geometry ، روی  Rectangle کلیک کنید .
2
در پنجره تنظیمات برای Rectangle ، بخش Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، W را تایپ کنید .
4
در قسمت متن ارتفاع ، H را تایپ کنید .
5
 روی Build  Selected کلیک کنید .
دایره 1 (c1)
1
در نوار ابزار Geometry ، روی  Circle کلیک کنید .
2
در پنجره تنظیمات برای دایره ، قسمت موقعیت را پیدا کنید .
3
در قسمت متن x ، 0.2 را تایپ کنید .
4
در قسمت متن y ، 0.2 را تایپ کنید .
5
قسمت Size  and  Shape را پیدا کنید . در قسمت متن Radius ، R را تایپ کنید .
6
 روی Build  Selected کلیک کنید .
تفاوت 1 (dif1)
1
در نوار ابزار Geometry ، روی  Booleans  and  Partitions کلیک کنید و Difference را انتخاب کنید .
2
فقط شی r1 را انتخاب کنید.
3
در پنجره تنظیمات برای تفاوت ، بخش تفاوت را پیدا کنید .
4
زیربخش اشیاء را  برای  تفریق پیدا کنید . برای انتخاب دکمه ضامن  فعال کردن  انتخاب کلیک کنید .
5
فقط شی c1 را انتخاب کنید.
6
در نوار ابزار Geometry ، روی  ساختن  همه کلیک کنید .
7
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
مواد
مواد 1 (mat1)
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Materials راست کلیک کرده و Blank  Material را انتخاب کنید .
2
در پنجره تنظیمات برای Material ، قسمت Material  Contents را پیدا کنید .
3
در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
تراکم
rho
1
کیلوگرم بر متر مکعب
پایه ای
ویسکوزیته دینامیکی
که در
1e-3
پس
پایه ای
جریان آرام (SPF)
ورودی 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Laminar  Flow  (spf) کلیک راست کرده و Inlet را انتخاب کنید .
2
فقط مرز 1 را انتخاب کنید.
ما با استفاده از پارامترهای از پیش تعریف شده، یک پروفیل سرعت سهموی تعریف می کنیم. با استفاده از تابع گامی که قبلاً تعریف شده است، سرعت را افزایش دهید. براکت واحد معکوس [1/s] را به متغیر زمان t اضافه کنید زیرا تابع step انتظار یک آرگومان بدون بعد را دارد.
3
در پنجره تنظیمات برای ورودی ، بخش Velocity را پیدا کنید .
4
در قسمت متنی 0 ، 6*U_mean*y*(Hy)/H^2*step1(t[1/s]) را تایپ کنید .
خروجی 1
1
در نوار ابزار Physics ، روی  Boundaries کلیک کنید و Outlet را انتخاب کنید .
2
فقط مرز 4 را انتخاب کنید.
مش 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Mesh  1 کلیک کنید .
2
در پنجره تنظیمات مش ، قسمت Physics-Controlled  Mesh را پیدا کنید .
3
از فهرست اندازه عنصر  ، Finer را انتخاب کنید .
4
 روی ساخت  همه کلیک کنید .
اگر روی ورودی و سیلندر بزرگنمایی کنید، می توانید لایه های مرزی را که مش کنترل شده فیزیک ایجاد می کند، ببینید.
مطالعه 1
مرحله 1: وابسته به زمان
1
در پنجره Model  Builder ، در بخش مطالعه  1 ، روی Step  1:  Time  Dependent کلیک کنید .
2
در پنجره تنظیمات مربوط به زمان  وابسته ، قسمت تنظیمات مطالعه  را پیدا کنید .
3
در قسمت متن زمان خروجی  ، range(0,0.2,3.4) range(3.5,0.02,7) را تایپ کنید .
راه حل 1 (sol1)
1
در نوار ابزار مطالعه ، روی  Show  Default  Solver کلیک کنید .
2
در پنجره Model  Builder ، گره Solution   (sol1) را گسترش دهید ، سپس روی Time-Dependent  Solver  1 کلیک کنید .
3
در پنجره تنظیمات برای حل وابسته به زمان  ، برای گسترش بخش Time Steping کلیک کنید .
4
از لیست مراحل  انجام شده  توسط  حل کننده ، Intermediate را انتخاب کنید .
5
در نوار ابزار مطالعه ،  روی محاسبه کلیک کنید .
نتایج
سرعت (spf)
برای بازتولید نمودار در شکل 1 ، یک ردیابی ذرات با گره جرمی را به اولین گروه نمودار پیش فرض اضافه کنید .
ردیابی ذرات با جرم 1
1
در نوار ابزار Velocity (spf) ، روی  More  Plots کلیک کنید و Particle  Tracing  with  Mass را انتخاب کنید .
2
در پنجره تنظیمات برای ردیابی ذرات  با جرم ، برای گسترش بخش Mass and Velocity کلیک کنید .
3
در قسمت Mass text، 4*pi/3*1e-9 را تایپ کنید .
4
زیربخش سرعت اولیه را  پیدا کنید . در قسمت متن x-component ، u را تایپ کنید .
5
در قسمت متن y-component ، v را تایپ کنید .
6
قسمت Particle  Positioning را پیدا کنید . در قسمت متن y ، range(0.1,0.05,0.3) را تایپ کنید .
7
در قسمت متن x ، 0 را تایپ کنید .
8
برای گسترش بخش Release کلیک کنید . از لیست Release  particles ، در  فواصل زمانی را انتخاب کنید .
9
تیک Start  time را انتخاب کنید . در فیلد متن مرتبط، 3.6 را تایپ کنید .
10
در قسمت متن زمان  بین  انتشار ، 0.4 را تایپ کنید .
11
برای گسترش بخش Coloring  and  Style کلیک کنید . زیربخش Line  style را پیدا کنید . از لیست Type ، None را انتخاب کنید .
12
زیربخش Point  style را پیدا کنید . از لیست نوع ، نقطه را انتخاب کنید .
13
زیربخش Point  motion را پیدا کنید . از لیست وقتی  ذرات  دامنه را ترک می کند  ، ناپدید شدن را انتخاب کنید .
14
در نوار ابزار Velocity (spf) ، روی  Plot کلیک کنید .
برای بازتولید شکل 2 و شکل 3 از ضرایب بالابر و درگ، ابتدا یک مجموعه داده انتگرال برای محاسبه نیروی واکنش کل روی سیلندر اضافه کنید.
انتگرال 1
در نوار ابزار نتایج ، روی  More  Datasets کلیک کنید و Evaluation>Integral را انتخاب کنید .
انتخاب
1
در نوار ابزار نتایج ، روی  ویژگی ها کلیک کنید و Selection را انتخاب کنید .
2
در پنجره تنظیمات برای انتخاب ، قسمت انتخاب موجودیت هندسی  را پیدا کنید .
3
از لیست سطح نهاد هندسی  ، Boundary را انتخاب کنید .
4
فقط مرزهای 5-8 را انتخاب کنید.
گروه طرح 1 بعدی 3
1
در نوار ابزار نتایج ، روی  1D  Plot  Group کلیک کنید .
2
در پنجره Settings for 1D  Plot  Group ، بخش Data را پیدا کنید .
3
از لیست Dataset ، Integral  1 را انتخاب کنید .
نمودار نقطه 1
1
روی 1D  Plot  Group  کلیک راست کرده و Point  Graph را انتخاب کنید .
2
در پنجره تنظیمات نمودار نقطه‌ای  ، بخش y-Axis Data را پیدا کنید .
3
در قسمت متن Expression ، (-reacf(v)[N]*2/(spf.rho*U_mean^2*(2*R)[m^2]))[1/m] را تایپ کنید .
4
چک باکس Description را انتخاب کنید . در فیلد نوشتاری مرتبط، Lift Coefficient را تایپ کنید .
5
در نوار ابزار 1D Plot Group 3 ، روی  Plot کلیک کنید .
نمودار را با نمودار نشان داده شده در شکل 2 مقایسه کنید .
گروه طرح 1 بعدی 4
در نهایت با استفاده از مراحل زیر ضریب درگ را تجسم کنید:
1
در نوار ابزار صفحه اصلی ، روی  Add  Plot  Group کلیک کنید و 1D  Plot  Group را انتخاب کنید .
2
در پنجره Settings for 1D  Plot  Group ، بخش Data را پیدا کنید .
3
از لیست Dataset ، Integral  1 را انتخاب کنید .
نمودار نقطه 1
1
روی 1D  Plot  Group  کلیک راست کرده و Point  Graph را انتخاب کنید .
2
در پنجره تنظیمات نمودار نقطه‌ای  ، بخش y-Axis Data را پیدا کنید .
3
در قسمت متن Expression ، (-reacf(u)[N]*2/(spf.rho*U_mean^2*(2*R)[m^2]))[1/m] را تایپ کنید .
4
چک باکس Description را انتخاب کنید . در قسمت متن مرتبط، ضریب Drag را تایپ کنید .
5
در نوار ابزار 1D Plot Group 4 ، روی  Plot کلیک کنید .
با شکل 3 مقایسه کنید .