%46تخفیف

شبیه­سازی سه­بعدی تأثیر جریان عرضی القایی بر عملکرد پیل سوختی غشاءپلیمری با آرایش میدان جریان موازی

تعداد 89 صفحه  درword


گروه آموزشی مهندسی مکانیک

 

کارشناسی ارشد در رشته مهندسی مکانیک گرایش تبدیل انرژی

 

 

شبیه­سازی عددی حباب در جریان ترکیبی برشی و پوآسوی با عدد رینولدز محدود در حالت سه بعدی

 

 

چكيده:

در پژوهش حاضر به شبیه­سازی عددی سه­بعدی مهاجرت عرضی یک حباب تغییر شکل­پذیر تحت جریان ترکیبی برشی ساده و پوآسوی در اعداد رینولدز محدود پرداخته شده است. معادلات ناویر-استوکس با رویکرد بقائی برای سیالات تراکم­ناپذیر با استفاده از روش اختلاف محدود روی یک شبکۀ منظم، ساکن و جابجا شده حل شده است. وجه مشترک بصورت صریح توسط اتصال نقاط نشانگر از طریق روش ردیابی جبهه روی یک شبکۀ نامنظم، مثلثی و متحرک ردیابی شده است. اثرات کشش سطحی نیز از طریق اضافه کردن یک جملۀ منبع مناسب به معادلات حاکم به­حساب آمده است. برای تزریق نیروهای کشش سطحی به شبکۀ ساکن از روشي به­نام مرز غوطه­ور استفاده شده است.  هدف از انجام شبیه­سازی­ها بررسی تأثیر اعداد بدون بعد حاکم، نظیر عدد کاپیلاری، عدد رینولدز، نسبت هندسی  بر نحوۀ مهاجرت جانبی است. نتایج نشان می­دهند که یک حباب صرفنظر از موقعیت اولیۀ آن در یک موقعیت تعادلی میان دیواره و خط مرکزی تثبیت خواهد شد. همچنین شبیه­سازی­ها نشان می­دهند که تغییر شکل حباب به­شدت وابسته به عدد کاپیلاری است. بنابراین، عدد بدون بعد مناسب برای کشش سطحی، عدد کاپیلاری خواهد بود. مشاهده شده است که با افزایش عدد رینولدز اثرات چسبندگی ضعیف­تر شده و نيروي روانسازي كاهش مي­يابد و درنتيجهحباب در یک موقعیت نزدیکتری نسبت به دیواره قرار خواهد گرفت. همچنین تأثیر شعاع حباب بر روی مهاجرت بررسی شده است. مشاهده می­گردد که افزایش شعاع حباب باعث خواهد شد كه مرکز حباب در یک وضعیت تعادلی نزدیکتری نسبت به خط مرکزی قرار بگيرد. گرادیان فشار منفی موجب خواهد شد که تغییر شکل حباب بیشتر شود و در نتیجه  حباب در یک وضعیت تعادلی نزدیکتری نسبت به خط مرکزی قرار بگیرد. تنش برشي بر روي ديواره­ها در حضور حباب، داراي تغييراتي وابسته به موقعيت اولية حباب مي­باشد. نتايج نشان داد كه با افزايش شعاع حباب و يا افزايش عدد كاپيلاري تنش برشي بر روي ديواره پايين كاهش مي­يابد. در مجموعۀ دوم از شبیه­سازی­ها، جملۀ نیروی گرانشی به معادلات ناویر-استوکس اضافه شده است. در اینقسمت نیز به تأثیر اعداد بدون بعد حاکم در حضور شتاب گرانشی، نظیر عدد فرود پرداخته می­شود. در اين مورد تغيير شكل حباب به­شدت وابسته به عدد اتوش مي­باشد.

واژه­های کلیدی:حباب، جریان ترکیبی برشی ساده و پوآسوی، روش اختلاف محدودردیابی جبهه، عدد کاپیلاری، عدد اتوش

 

 

فهرست مطالب

عنوان                                           صفحه

فهرست مطالب. الف

فهرست جداول. د

فهرست شکل­ها. ه

فهرست علائم اختصاری. ط

فصل اول: مقدمه. 1

1-1-                   طبقهبنديجريان. 4

1-1-1-   براساسرژيمجريان. 4

1-1-2-   برحسبعددرينولدز. 5

1-1-3-   سایرتقسیمبندی­ها. 6

1-1-4-   مکانیزم­هایمهاجرتذرات. 6

1-2-                   جریان­هایبنیادین. 7

1-2-1-   جریانبرشیسادهیاجریانکوئت. 7

1-2-2-   جریانپوآسوی. 7

1-3-                                     جریانترکیبیبرشیوپوآسوی. 7

فصل دوم: بررسی منابع. 9

2-1-    مقدمه. 11

2-2-                   پیشینۀپژوهش. 11

2-3-  فرضیاتپژوهش. 24

2-4-               اهدافپژوهش. 24

فصل سوم: معادلات حاکم. 25

3-1-………………………………………………………………………………… مقدمه. 27

3-2-                   دینامیکحباب. 28

3-3-               معادلاتحاكم. 29

3-4-            شرايطمرزي. 32

3-5-                      پارامترهايبدونبعد. 32

3-6-                                                                ماهيتكششسطحيدرمعادلاتناوير- استوكس  34

3-7-                                  معادلاتحاكمبرميدانفشار. 35

3-8-                                                             معادلاتحاكمبرميدانهايچگاليوچسبندگي  35

فصل چهارم: روش عددی. 39

4-1-………………………………………………………………………………… مقدمه. 41

4-2-                                مختصريدربارهروشهايVOFوLevel Set 44

4-3-                   روشردیابیجبهه. 46

4-4-                        انتگرال­گیریدرزمان. 49

4-5-               ساختارجبهه. 50

4-6-               بازسازيجبهه. 52

4-7-                                                                 انتقالوپوشاندنخواصجبههبررويشبکهساکن  53

4-8-                          بهروزکردنخواصمواد. 55

4-9-                          محاسبۀکششسطحي. 57

4-9-1-وجهمشترکدوبعدی. 58

4-9-2-وجهمشترکسهبعدی. 60

4-10-حلمعادلهفشار. 61

4-11-پيشرويجبهه. 62

فصل پنجم: بحث و نتایج. 65

5-1-………………………………………………………………………………… مقدمه. 67

5-2-                                                                            شبیه­سازیمهاجرتجانبییکحباببدوننیرویشناوری. 68

5-2-1-   استقلالشبکه. 69

5-2-2-   نتایجکلی. 70

5-2-3-   تأثیرعددکاپیلاری. 78

5-2-4-   تأثیرعددرینولدز. 84

5-2-5-   تأثیراندازۀحباب. 86

5-2-6-   تأثیرگرادیانفشار. 90

5-2-7-   تأثیرنسبتچسبندگی. 94

5-2-8-   تأثیرنسبتچگالی. 95

5-2-9-   تأثیرحضورحباببررویتنشبرشیدیواره. 96

5-3-                                                                                           شبیه­سازیمهاجرتجانبییکحباببادرنظرگرفتننیرویشناوری. 100

5-3-1-   مطالعۀشبکه. 101

5-3-2-   نتایجکلی. 101

فصل ششم: نتیجه­گیری و پیشنهادات. 111

6-1-            نتیجه­گیری. 113

6-2-            پیشنهادات. 115

6-2-1-شبیه­سازیانتقالحرارت. 116

6-2-2-شبیه­سازیدویاچندحباب. 116

6-2-3-شبیه­سازیلایهمرزییکحباب. 116

6-2-3-شبیه­سازیصعودیکحبابدریککانالعمودی. 117

6-2-5-استفادهازروشهایدیگرشبیه­سازیجریانهایچندفازی. 117

پیوست (الف): گسسته سازی معادلات حاکم . 119

1-مقدمه. 121

2-گسسته­سازیمکانییافضایی. 121

3- گسسته­سازیجملۀجابجایی. 125

4-گسسته­سازیجملۀنفوذ. 126

5- چسبندگی­هایمختلف. 127

6- معادلۀفشار. 128

7-  دامنۀمحاسباتی. 130

8-  شرایطمرزی. 131

9- الگوریتمعددی. 134

10- انتگرال­گیریمرتبۀدومدرزمان. 136

پیوست (ب): کد عددی. 137

1- مقدمه. 139

2-  زیرروالمحاسبۀجملاتجابجاییوچسبندگی. 139

3- زیرروالبرایحلمعادلۀفشار. 142

4- زیرروالشرطپریودیکدرجهاتxوy. 145

منابع و مراجع. 147

 

فهرست جداول

 

عنوان جدول                                      صفحه

جدول‏1‑1 مکانیزم­هایمهاجرتجانبیذراتصلبوشکل­پذیردرجریان­هایمختلف   6

جدول‏4‑1 توابعهمواربرایانتقالاطلاعاتمیانجبههوشبکۀثابت [51].   55

 

 

 

فهرست اشكال

 

عنوان شکل                                       صفحه

شکل‏1‑1 جریانترکیبیکوئتوپوآسویبینصفحاتموازی [6]… 8

شکل‏2‑1جزئیاتمیدانجریاندرزمان 4/64 [16] … 14

شکل‏2‑2 حالتپایدارصعودحباببرایاعداداتوشومورتونمختلف. سطربالاEo=1(M=10-4,10-5,10-6,10-7)،سطروسطEo=10(M=10-1,10-2,10-3,10-4)،وسطرپایینEo=104(M=102,10,1,10-1) [18].   15

شکل‏2‑3 خطوطجریانبرایشکل­هایردیفبالاوردیفوسطشکل‏2‑2 [18]… 16

شکل‏2‑4 میدانفشاربرایEo=1: (الف) M=10-8و (ب) M=10-5 . فشاربعلتاثرهیدرواستاتیکبهسمتپایینافزایشمی­یابد [18]… 17

شکل‏2‑5 تنشبرشیمتوسطبررویدیوارۀبالاوپاییندرمقابلزمانبرایجریانبا،وبدونحباب [28]… 20

شکل‏2‑6 تأثیرنسبتچگالیبررویسرعتصعودحباب [36]… 22

شکل‏3‑1 جریانچندفازی. سیالاتمختلفمی­توانندیکمنطقۀمتصلویاجداشدهازهمرادرسیالدیگربعنوانحبابویاقطرهتشکیلدهند… 27

شکل‏3‑2 هندسةحلبرايشبيه­سازيحركتيكحبابدريککانالپريوديک… 32

شکل‏3‑3 یکحجممعیارکوچکیازجبهه.. 34

شکل‏4‑1 موقعيتشبکهلاگرانژيرويشبکهساکناويلريبرايبيشازيكفاز [39].   48

شکل‏4‑2 ساختارجبهه،الف) جبههدوبعدیوب) جبههسه­بعدي [39]… 50

شکل‏4‑3 انتقالاطلاعاتهرگرهازجبههبهچهارگرهمجاورآندرشبکهساکن [51].   51

شکل‏4‑4 بازسازيجبههدرحالتسهبعدي [49]… 53

شکل‏4‑5 تفسیرضرایبوزنی… 54

شکل‏4‑6 الف) تاخوردنمرزجبههرويخودشوب) نزديکشدنبيشازحددوجبههبههمديگر.   56

شکل‏4‑7 ایجادتابعنشانگرازگرادیانآنبررویشبکۀجابجاشده… 57

شکل‏4‑8 محاسباتانحناءدوبعدی… 59

شکل‏4‑9 توزیعنیروبراییکوجهمشترکردیابیشدهتوسطیکحلقۀ(Chain)نقاطنشانگر[51].   61

شکل‏5‑1 الف) هندسةجريانموردمطالعهوب) شبکه­بندیمورداستفادهبرایشبیه­سازیحرکتیکحبابدریککانالباشرطمرزیپریودیک… 68

شکل‏5‑2 اثرکیفیتشبکهبررویمهاجرتجانبییکحبابدرجریانترکیبیبرشیسادهوپوآسویبراییکحبابدرچهارشبکۀمختلف. الف) موقعیتعرضیبرحسبموقعیتمحوری،ب) سرعتمحوریبرحسبزمانبی­بعد،ج) تغییرشکلحباببرحسبزمانبدونبعدود) نرخمهاجرتبرحسبزمانبدونبعد. پارامترهایجریانعبارتنداز: 10Red=، 9/0Ca=، 8/0α=λ=… 70

شکل‏5‑3 خطوطجریاندرمقطعمیانیدرجهتyشرایطجریانعبارتنداز: (10Red=، 5/2We=، 9/0α=λ=، 125/0ξ=)… 71

شکل‏5‑4 الف) موقعیتجانبی،ب) سرعتمحوری،ج) نرخمهاجرتود) سرعتلغزشیبرحسبزمانبدونبعد… 72

شکل‏5‑5 مقایسهباشبیه­سازیانجامشدهتوسطفنگوهمکارانش [20] برایجریانکوئت.   74

شکل‏5‑6 مقایسهباشبیه­سازیسه­بعدیبیارهومرتضوی [22] برایقطرهدرجریانکوئت.   74

شکل‏5‑7 مقایسهباشبیه­سازیانجامشدهتوسطفنگوهمکارانش [20] برایجریانپوآسوی… 75

شکل‏5‑8 مقایسهباشبیه­سازیسه­بعدینوربخشومرتضوی [26] برایقطرهدرجریانپوآسوی… 76

شکل‏5‑9 تغییرشکلیکحبابتحتجریانترکیبیبرشیسادهوپوآسويدرکاپیلاریثابت 3/0Ca=… 77

شکل‏5‑10 تغییرشکلیکحبابتحتجریانترکیبیبرشیسادهوپوآسويدروبرثابت 5/2We=.   78

شکل‏5‑11 تأثیرعددکاپیلاریبرالف) موقعیتعرضیوب) سرعتمحوریبرحسبزمانبدونبعددرعددرینولدز 10Red=… 79

شکل‏5‑12 تأثیرعددکاپیلاریبرروینرخمهاجرتحباب… 79

شکل‏5‑13 تحولتدریجیحبابدرزمان­هایمختلف… 80

شکل‏5‑14 الف) سرعتلغزشی،وب) تغییرشکلاسکالرتیلوربرحسبزمانبدون­بعد.   81

شکل‏5‑15 خطوطجریاندرمقطعمیانیدرجهتyبراییکحباببااعدادکاپیلاری: الف) 3/0Ca=وب) 7/0Ca= . پارامترهایجریان 10Red=، 8/0α=λ=و 125/0ξ=هستند.   83

شکل‏5‑16 فرآیندشکستهشدنحبابدراعدادبدونبعد: 10Red=، 45/0α=λ=، 2Ca=و 125/0ξ=… 84

شکل‏5‑17 تأثیرعددرینولدزبرالف) موقعیتعرضیبرحسبموقعیتمحوری،ب) سرعتمحوری،ج) سرعتلغزشیود) تغییرشکلتیلوربرحسبزمانبدونبعددرسهعددرینولدزمختلف… 86

شکل‏5‑18 تأثیراندازۀحباببررویالف)موقعیتعرضیبرحسبموقعیتطولیوب)نرخمهاجرتبرحسبزمانبدون­بعدبرایحبابباسهاندازۀمختلفدر 10Red=، 2/0Ca=و 8/0α=λ=.   87

شکل‏5‑19 تأثیراندازۀحباببررویتغییرشکلتیلور… 88

شکل‏5‑20 الف) موقعیتعرضیدرمقابلزمانبدونبعدبرایمقایسهباجریانکوئتوب) موقعیتعرضیدرمقابلموقعیتمحوریبرایمقایسهباجریانپوآسوی.. 88

شکل‏5‑21 خطوطجریانبرایحباببانسبتهندسی 18/0ξ=درزمانبدون­بعد 20τ=.   89

شکل‏5‑22 خطوطجریانبرایحباببانسبتهندسی 22/0ξ=درزمانبدون­بعد 20τ=.   90

شکل‏5‑23 پروفیلسرعتبرحسبارتفاعکانال… 91

شکل‏5‑24 تأثیرگرادیانفشاربررویالف) موقعیتعرضی،وب) تغییرشکلتیلوربرحسبزمانبدونبعد… 91

شکل‏5‑25 تأثیرگرادیانفشاربررویالف) نرخمهاجرت،ب) سرعتلغزشیبرحسبزمانبدونبعد… 92

شکل‏5‑26 پروفیلسرعتبرایجریانمختلنشده (خطممتد) وجریانباحضورحباب (خطمربع)،الف) گرادیانفشارمنفیوب) گرادیانمثبتفشار… 93

شکل‏5‑27 خطوطجریانبرایالف) گرادیانفشارمنفی،وب) گرادیانفشارمثبت.   94

شکل‏5‑28 تأثیرنسبتچسبندگیبررویالف) مهاجرتجانبیوب) تغییرشکلتیلور   95

شکل‏5‑29 تأثیرنسبتچگالیبررویالف) مهاجرتجانبیوب) سرعتمحوری.   95

شکل‏5‑30 تأثیرنسبتچگالیبرروینرختغییرشکل… 96

شکل‏5‑31 تأثیرموقعیتاولیۀحباببررویتنشبرشیبررویدیوارۀپایین.   97

شکل‏5‑32 تأثیرموقعیتاولیۀحباببررویتنشبرشیبررویدیوارۀبالایی.   97

شکل‏5‑33 تأثیرشعاعحباببررویتنشبرشیرویدیوارۀپایین… 98

شکل‏5‑34 تأثیرعددکاپیلاریبررویتنشبرشیرویدیوارۀپایین… 98

شکل‏5‑35 الف) تاثيرشعاعحبابوب) عددكاپيلاريبررويتنشبرشيبررويديوارةبالا.   99

شکل‏5‑36 شماتیکجریانبرشیروبهبالاوجریانبرشیروبهپایین… 100

شکل‏5‑37 تأثیرکیفیتشبکهبررویالف) مهاجرتجانبی،ب) تغییرشکلیکحباببرایچهارشبکۀمختلف. پارامترهایحاکم: 8/0α=λ=، 10Red=، 5/0Ca=و 10Fr=… 101

شکل‏5‑38 بررسیتاثیرجهتجاذبهبررویمهاجرتعرضیحباب… 102

شکل‏5‑39 تأثیرجهتنیرویگرانشبررویسرعتلغزشیبرحسبزمانبدون­بعد.   103

شکل‏5‑40 تأثیرعددفرودبرروی (الف) مهاجرتجانبی،و (ب) سرعتلغزشی   103

شکل‏5‑41 تأثیرعددفرودبرنرخمهاجرتحباب… 104

شکل‏5‑42 تأثیرعددفرودبرتغییرشکلحباب… 104

شکل‏5‑43 خطوطجریاندرعددفرود: الف) 10Fr=وب) 50Fr=… 105

شکل‏5‑44 تأثیرعدداتوشبررویشکلحباب. الف) 1Eo=،ب) 10Eo=وج) 100Eo=.   106

شکل‏5‑45 تأثیرعدداتوشبررویخطوطجریان. الف) 1Eo=،ب) 10Eo=وج) 100Eo=.   108

شکل‏6‑1 شماتیکلایهمرزییکحباب. 116

شکل (الف)- 1 طرحوارۀشبکۀساختاریافتۀمنظم… 122

شکل (الف)- 2 نمادگذاریمورداستفادهبراییکشبکۀجابجاشدۀاستاندارد. فرضشدهکهفشاردرمرکزحجمکنترلیکهپیرامونشیکخطصلبضخیمکشیدهشدهاست،شناختهشدهمی­باشد… 123

شکل (الف)- 3 نمادگذاریمورداستفادهبراییکشبکۀجابجاشدۀاستاندارد. مؤلفه­هایسرعتافقیدروسطلبه­هایسمتچپوراستذخیرهشده­اند،ونیزمؤلفه­هایسرعتعمودیدروسطلبه­هایبالاوپایینذخیرهشده­اند… 124

شکل (الف)- 4 نمودارباقیماندۀفشار.. 130

شکل (الف)- 5 نمادگذاریاستانداردبراییکشبکۀجابجاشده. فرضشدهاستکهفشاردرمرکزحجممعیارهاییکهاطرافآنیکخطصلبضخیمکشیدهشدهاست،معلومهستند. مؤلفه­هایسرعتافقیدروسطلبه­هایسمتچپوراستذخیرهشدهاند،ونیزمؤلفه­هایسرعتعمودیدروسطلبه­هایبالاوپایینذخیرهشده­اند… 131

شکل (الف)- 6 یکحجممعیاردرکنارمرزیکهسرعتنرمالآنشناختهشدهاست.   132

شکل (الف)- 7 مرزجریانخروجی… 134

 

 

فهرست علائم اختصاری

a

شعاعحباب

b

کوچکترین فاصلۀ سطح تا مرکز حباب

Ca

عدد کاپیلاری

d

قطر حباب

D

تغییر شکل تیلور

مشتق مادی

Eo

عدد اتوش

f

نیروی حجمی به غیر از جاذبه

نیروی ناشی از کشش سطحی

Fr

عدد فرود

g

شتاب جاذبه

h

اندازۀ هر المان

H

ارتفاع کانال

تابع نشان گذار

L

بزرگترین فاصلۀ سطح تا مرکز حباب

M

عدد مورتون

n

بردار عمود بر سطح المان

p

فشار

گرادیان فشار

عدد رینولدز حجمی

عدد رینولدز بر مبنای قطر حباب

t

بردار مماس بر لبۀ هر المان

u

میدان سرعت

سرعت ماکزیمم در کانال

We

عدد وبر

موقعیت اویلری

موقعیت لاگرانژی

فهرست علائم يوناني

نسبت چگالی

تعیین کننده دو یا سه بعدی بودن جریان

تابع دلتا

المان کوچک سطح مشترک

مساحت المان

کمیت سطح مشترک

کمیت شبکه

k

انحناء

نسبت چسبندگی

چسبندگی

چسبندگی سیال داخلی (حباب)

چسبندگی سیال بیرونی (محیط)

چگالی

چگالی سیال داخلی (حباب)

چگالی سیال بیرونی (محیط)

ضریب کشش سطحی

زمان بدون بعد شده

اندازۀ بی­بعد شعاع حباب یا نسبت هندسی

نقد و بررسی‌ها

هنوز بررسی‌ای ثبت نشده است.

اولین کسی باشید که دیدگاهی می نویسد “شبیه­سازی سه­بعدی تأثیر جریان عرضی القایی بر عملکرد پیل سوختی غشاءپلیمری با آرایش میدان جریان موازی”

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *

قبلا حساب کاربری ایجاد کرده اید؟
گذرواژه خود را فراموش کرده اید؟
Loading...
enemad-logo