شبیه‌سازی جریان سیال و انتقال حرارت در کامسول (COMSOL)

در این بلاگ به شبیه‌سازی جریان سیال و انتقال حرارت در کامسول (COMSOL) می‌پردازیم. یک کانال دوبعدی دارای جریان آرام آب و دیواره گرم را از ابتدا مدل‌سازی می‌کنیم. هر ورودی، مرز، پارامتر و خروجی را با دلیل انتخابش توضیح می‌دهیم تا خواننده بتواند مدل را مستقل بازسازی کند. تصاویر هندسه، مش و پروفیل سرعت در این مطلب «شماتیک آموزشی» هستند، نه اسکرین‌شات یا خروجی محاسبه‌شده COMSOL؛ اعدادِ مرجع نیز از روابط تحلیلی به دست آمده‌اند و باید پس از اجرای واقعی مدل بررسی شوند.

۱. پیش از شروع؛ قرار است چه چیزی بسازیم؟

آب از چپ وارد کانالی با طول ۵۰ میلی‌متر و ارتفاع ۱ میلی‌متر می‌شود. کانال بین دو دیواره جامد قرار دارد؛ از سطح بیرونی دیواره پایینی گرما وارد می‌کنیم. نرم‌افزار سرعت و فشار جریان، دمای آب و دمای دو دیواره را محاسبه می‌کند. در پایان بررسی می‌کنیم که سرعت روی دیواره صفر است، افت فشار معقول است و گرمای ورودی با افزایش انرژی جریان سازگار است.

پیش‌نیاز: دسترسی به COMSOL Multiphysics 5.6 و مجوز رابط‌های Laminar Flow و Heat Transfer in Solids and Fluids / Conjugate Heat Transfer. اگر گزینه‌ای در Model Wizard دیده نمی‌شود، نخست فهرست ماژول‌های نصب‌شده را بررسی کنید؛ مسیرهای این آموزش بر مستندات رسمی 5.6 تکیه دارند.

۱-۱. راهنمای شناخت پنجره COMSOL

Model Builder در سمت چپ مانند فهرست مراحل مدل است؛ هر گره بخشی از ساخت مدل را کنترل می‌کند. Settings در سمت راستِ انتخاب هر گره ظاهر می‌شود و تنظیمات همان گره را نشان می‌دهد. Graphics محل نمایش هندسه، انتخاب ناحیه‌ها و نتایج است. نوار Home و Physics ابزارهای افزودن فیزیک، حل و رسم نمودار را دارند. هر بار که روی گره‌ای کلیک می‌کنید، پیش از واردکردن عدد عنوان بالای Settings را کنترل کنید.

۱-۲. پارامترهای ثابت پروژه

واحدها را همراه مقدار داخل براکت بنویسید. COMSOL واحدها را بررسی می‌کند؛ مثلاً اگر برای سرعت به‌اشتباه m را وارد کنید، ناسازگاری واحدی رخ می‌دهد. این اعداد یک مجموعه آموزشی با خواص ثابت‌اند، نه داده اندازه‌گیری‌شده یک نمونه واقعی.

نام پارامتر	عبارت دقیق ورود	معنی
L	0.05[m]	طول کانال
H	0.001[m]	ارتفاع سیال
ts	0.0005[m]	ضخامت هر دیواره
Um	0.01[m/s]	سرعت متوسط ورودی
Tin	293.15[K]	دمای ورودی
q0	100[W/m^2]	شار گرمایی کف
rho_f	1000[kg/m^3]	چگالی سیال
mu_f	0.001[Pa*s]	ویسکوزیته دینامیکی سیال
cp_f	4180[J/(kg*K)]	گرمای ویژه سیال
k_f	0.6[W/(m*K)]	رسانندگی سیال
rho_s	7800[kg/m^3]	چگالی جامد
cp_s	500[J/(kg*K)]	گرمای ویژه جامد
k_s	15[W/(m*K)]	رسانندگی جامد

۲. ایجاد فایل مدل در نسخه 5.6

۲-۱. ورود به Model Wizard

۱) نرم‌افزار COMSOL 5.6 را باز کنید.

۲) New → Model Wizard را انتخاب کنید.

۳) بُعد 2D را تعیین کنید.

 ۴) در Add Physics از شاخه Heat Transfer → Conjugate Heat Transfer، گزینه Laminar Flow را انتخاب و Add کنید.

 ۵) در Study، گزینه Stationary و سپس Done را بزنید. اگر شاخه Conjugate Heat Transfer در مجوز شما نبود، می‌توانید Laminar Flow و Heat Transfer in Solids and Fluids را جداگانه اضافه کنید و سپس کوپل Nonisothermal Flow را از Add Multiphysics بیفزایید؛ تطابق رابط‌ها را قبل از ادامه کنترل کنید.

چرا Stationary؟ در این تمرین فرض می‌کنیم با گذشت زمان پاسخ به حالت پایا رسیده است. اگر هدف گرم‌شدن از لحظه روشن‌شدن سامانه باشد، باید مطالعه Time Dependent بسازیم و شرایط اولیه را جدی‌تر مشخص کنیم.

کنترل پایان مرحله: در Model Builder باید بخش‌های Component 1، Geometry، Materials، Laminar Flow (spf)، Heat Transfer in Solids and Fluids (ht)، Multiphysics با Nonisothermal Flow، Mesh و Study قابل مشاهده باشند. اگر افزودن فیزیک به‌صورت جدا انجام شد، وجود کوپل Nonisothermal Flow را مستقلاً کنترل کنید.

۲-۲. ورود پارامترها

در Model Builder روی Global Definitions راست‌کلیک کنید و Parameters را باز کنید. هر سطر جدول را در ستون Name و Expression وارد کنید. در ستون Description می‌توانید توضیح فارسی یا انگلیسی کوتاه بنویسید. پس از پایان، از File → Save As فایلی مانند heated_channel_56.mph بسازید. از این‌جا به بعد مرتب فایل را ذخیره کنید.

۳. ساخت هندسه کانال

۳-۱. سه مستطیلِ متصل بسازید

در Geometry 1 واحد طول را m قرار دهید. سپس با راست‌کلیک Geometry 1 → Rectangle، اولین مستطیل را بسازید: Width=L، Height=H، Position=(0,0). این ناحیه آب است. Rectangle دوم برای دیواره زیرین: Width=L، Height=ts، Position=(0,-ts). Rectangle سوم برای دیواره بالایی: Width=L، Height=ts، Position=(0,H). بعد از هر مستطیل Build Selected را بزنید تا شکل را ببینید.

در پایان روی Form Union کلیک و Build All کنید. در صورت نمایش گزینه Keep interior boundaries آن را فعال نگه دارید تا سه دامنه از هم قابل تشخیص باشند و تماس حرارتی در مرزهای مشترک حفظ شود. Form Assembly برای این تمرین لازم نیست.

۳-۲. انتخاب ناحیه و مرز را یاد بگیریم

در نوار Graphics دکمه Selection را روی Domain بگذارید و روی قسمت آبیِ کانال کلیک کنید. COMSOL شناسه عددی دامنه را نشان می‌دهد؛ شناسه‌ها ممکن است بعد از تغییر هندسه عوض شوند. سپس انتخاب را روی Boundary قرار دهید و به‌ترتیب ضلع چپ سیال، ضلع راست سیال و لبه بیرونی دیواره پایین را انتخاب کنید. انتخاب مرز بیرونی گرم‌شونده با سطح مشترک سیال و جامد متفاوت است؛ این اشتباه یکی از دلایل رایج نتایج نادرست است.

توصیه: در Definitions → Selections از Explicit Selection برای fluid، solid_lower، solid_upper، inlet، outlet و heated_bottom استفاده کنید. اگر ساختن Selection را هنوز بلد نیستید می‌توانید مرزها را دستی انتخاب کنید، ولی هر بار پس از ویرایش هندسه دوباره آن‌ها را کنترل کنید.

۴. مواد و خواص آن‌ها

۴-۱. تعریف ماده سیال

روی Materials راست‌کلیک و Blank Material را انتخاب کنید. نام آن را Educational water (constant properties) بگذارید. فقط ناحیه وسط را در Domain Selection انتخاب کنید. در Material Contents، Density را rho_f، Dynamic viscosity را mu_f، Heat capacity at constant pressure را cp_f و Thermal conductivity را k_f تنظیم کنید. اگر ویژگی به‌صورت پیش‌فرض در جدول نبود، آن را در رابط فیزیکی مربوط وارد کنید و از ناسازگاری میان مقدار ماده و فیزیک جلوگیری کنید.

۴-۲. تعریف دیواره جامد

یک Blank Material دیگر برای دو مستطیل بالا و پایین بسازید. مقدار Density=rho_s، Heat capacity=cp_s و Thermal conductivity=k_s را وارد کنید. مطمئن شوید دامنه سیال در ماده جامد انتخاب نشده و برعکس. در مسئله پایای حاضر چگالی و گرمای ویژه جامد ممکن است مستقیماً در پاسخ پایا نقش نداشته باشند، اما برای توسعه زمانمند ضروری می‌شوند.

۵. تنظیم جریان آرام Laminar Flow

۵-۱. دامنه جریان را محدود کنید

روی Laminar Flow (spf) کلیک کنید و در Domain Selection فقط مستطیل میانی را انتخاب کنید. در Fluid Properties مطمئن شوید چگالی و ویسکوزیته از پارامترهای ثابت تعریف‌شده خوانده می‌شوند. اگر درخت مربوطه Fluid Properties 1 دارد، همان را بررسی کنید. دو دیواره جامد نباید جزو دامنه جریان باشند.

۵-۲. ورودی، خروجی و دیواره

در Model Builder روی Laminar Flow (spf) راست‌کلیک و Inlet را اضافه کنید؛ فقط ضلع چپ دامنه سیال را انتخاب کنید. برای مدل دوبعدی حاضر، پروفیل سهموی ورودی مطابق رابطه (1) تعریف می‌شود. در Settings قسمت Velocity، نوع شرط را روی Velocity قرار دهید. اگر گزینه Velocity field در دسترس است، مؤلفه افقی u را با رابطه (1) و مؤلفه عمودی v را برابر صفر تنظیم کنید. اگر فقط Normal inflow velocity یا U0 نمایش داده می‌شود، مقدار نرمال و علامت آن به جهت نرمال خارجیِ ورودی وابسته است؛ در این حالت از گزینه مناسب برای تعریف مؤلفه‌های سرعت استفاده کنید یا راهنمای Inlet همان پنجره را بررسی کنید، نه اینکه همان عبارت را بدون کنترل علامت وارد کنید. برای بررسی، پیکان سرعت در ورودی باید رو به راست باشد؛ بیشینه سرعت در میانه کانال مطابق رابطه (2) است.

	(1)
	(2)

برای خروجی، Outlet را روی لبه راست سیال تعریف و فشار گیج را 0[Pa] بگذارید. در Wall مطمئن شوید مرزهای بالا و پایین سیال No slip هستند. No slip یعنی سرعت سیال روی دیواره ساکن صفر است. هیچ شرط Inlet یا Outlet روی لبه دیواره جامد اعمال نکنید.

۵-۳. چرا جریان آرام معتبر است؟

برای کانال تختِ پهن، قطر هیدرولیکی تقریبی برابر 2H است. عدد رینولدز از رابطه (3) محاسبه می‌شود و با مقادیر این پروژه برابر 20 است. این عدد برای مدل آموزشی حاضر بسیار کوچک است و انتخاب جریان آرام را توجیه می‌کند. با وجود این، تعریف قطر هیدرولیکی باید متناسب با هندسه انتخاب شود؛ در مدل سه‌بعدی با عرض محدود، تقریب 2H همیشه دقیق نیست.

(3)

۶. تنظیم انتقال حرارت و کوپل دو فیزیک

۶-۱. دامنه سیال و جامد

گره Heat Transfer in Solids and Fluids (ht) باید هر سه دامنه را پوشش دهد: آب و دو دیواره. گره Fluid برای آب و Solid برای دیواره‌ها را بررسی کنید. اگر رابط فقط برخی دامنه‌ها را انتخاب کرده، انتخاب‌ها را اصلاح کنید. در Multiphysics → Nonisothermal Flow مطمئن شوید جفت رابط‌های spf و ht درست متصل‌اند تا سرعت محاسبه‌شده در انتقال حرارت همرفتی به‌کار رود.

۶-۲. دمای ورودی و شار گرمایی

در Heat Transfer in Solids and Fluids، از مسیر راست‌کلیک روی رابط حرارت، شرط Inflow را فقط به ضلع چپ سیال اضافه کنید و دمای سیال ورودی را Tin قرار دهید. برای لبه خارجیِ پایینیِ دیواره جامد، Heat Flux را اضافه کنید؛ نوع شار را General inward heat flux و مقدار q0 را تعیین کنید. در این حالت q0 مثبت به‌معنی ورود گرما به دامنه است. مرز گرم‌شونده، سطح بیرونی جامد است؛ نه فصل مشترک آب و دیواره. سایر سطوح بیرونی جامد را Thermal Insulation در نظر بگیرید، مگر آنکه شرط دیگری صریحاً انتخاب شده باشد. توجه: انتخاب‌ها در گره‌های موجود را نیز کنترل کنید تا شرط‌های هم‌پوشان ناخواسته باقی نمانند.

برای خروجی سیال از Outflow استفاده کنید. دقت کنید دمای ورودی با ورودی سرعت متفاوت است؛ هر کدام در رابط فیزیکی خودش تنظیم می‌شود. روی فصل مشترک آب و جامد در Form Union، پیوستگی دما و شار را حفظ کنید و بدون دلیل مقاومت تماس اضافه نکنید.

۷. ابتدا جریان را جداگانه حل کنیم

برای اولین آزمون، ابتدا در Mesh 1 یک شبکه پیش‌فرض بسازید (Build All) و از نبود خطای مش مطمئن شوید. سپس q0 را موقتاً به 0[W/m^2] تغییر دهید. در Study → Step 1: Stationary و بخش Physics and Variables Selection، در صورت امکان حل Heat Transfer را غیرفعال و Laminar Flow را فعال نگه دارید؛ سپس Compute را انتخاب کنید. اگر کوپل فعال مانع این آزمون جداگانه شد، موقتاً Nonisothermal Flow را نیز غیرفعال کنید. پس از کنترل میدان سرعت و فشار، رابط حرارت و کوپل را برای حل نهایی دوباره فعال کنید.

نتیجه اول باید قابل پیش‌بینی باشد: شرط عدم لغزش روی دیواره‌ها برقرار است؛ نزدیک محور میانی سرعت به بیشینه تحلیلی می‌رسد و در کانال بلند با ورودی توسعه‌یافته، مؤلفه عمودی سرعت باید نزدیک صفر باشد. این رفتارها در رابطه (4) خلاصه شده‌اند. برای دیدن میدان از Results → 2D Plot Group → Surface مقدار spf.U یا مؤلفه مناسب سرعت را انتخاب کنید؛ نام دقیق متغیر را از Replace Expression در نسخه نصب‌شده بردارید.

(4)

برای نمودار مقطعی، از Results → Data Sets → Cut Line 2D در x=L/2 و y=0 تا H استفاده کنید؛ سپس 1D Plot Group → Line Graph را به آن مجموعه داده مرتبط و سرعت افقی را رسم کنید. اگر خروجی خطی پله‌ای یا وارونه است، ترتیب نقاط خط و واحد محورها را بررسی کنید.

۸. شبکه‌بندی و اطمینان از دقت عددی

۸-۱. شبکه مستطیلی در Mesh 1

برای کانال مستطیلی، Mapped Mesh گزینه مناسبی برای کنترل ساختار المان‌هاست. روی Mesh 1 راست‌کلیک کنید و Mapped را به دامنه‌های موردنظر بیفزایید. برای تعیین تعداد تقسیم روی لبه‌ها، از Distribution استفاده کنید؛ لبه‌های روبه‌رو باید تقسیم‌بندی سازگار داشته باشند. برای آزمون اولیه Nx=100 در طول، Ny=10 در ارتفاع سیال و Ns=5 در هر دیواره را هدف بگیرید. اگر ساخت Mapped Mesh با خطای تقسیم‌بندی مواجه شد، ابتدا مطمئن شوید سه مستطیل به‌صورت Union با مرزهای داخلی محفوظ ساخته شده‌اند؛ در غیر این صورت برای شروع Free Triangular با اندازه المان مناسب را انتخاب کنید. در پایان Build All را بزنید. این اعداد آموزشی‌اند و باید با مطالعه استقلال از شبکه بررسی شوند.

شبکه	Nx طولی	Ny سیال	Ns هر دیواره	تعداد چهارضلعی برنامه‌ریزی‌شده
درشت	100	10	5	2000
متوسط	200	20	10	8000
ریز	400	40	20	32000

۸-۲. استقلال از شبکه یعنی چه؟

مدل را روی سه شبکه اجرا و افت فشار، بیشینه سرعت، دمای توده خروجی و بیشینه دمای جامد را ثبت کنید. اگر تغییر نسبی خروجی بین دو شبکه متوالی کمتر از حدود ۱ درصد شد، می‌توان آن را یک معیار آموزشی اولیه دانست؛ برای نتیجه پژوهشی باید تلرانس با حساسیت مسئله، عدم‌قطعیت و هزینه حل تنظیم شود. درصد تغییر را با رابطه (5) حساب کنید و در صورت نزدیک‌صفر بودن مخرج، معیار خطای مطلق به‌کار برید.

(5)

۹. حل کوپل‌شده و مشاهده دما

q0 را دوباره 100[W/m^2] بگذارید. هر دو رابط جریان و حرارت و Nonisothermal Flow را در Study فعال کنید. Study 1 → Compute را اجرا کنید. اگر حل ناموفق بود، ابتدا Selectionها، واحدها، خواص ماده، وجود مرز فشار مرجع، جهت شار حرارت و شکل شبکه را بررسی کنید. تغییر بی‌دلیل خواص ماده برای رسیدن به همگرایی، حل علمی مشکل نیست.

از Results → 2D Plot Group → Surface برای رسم دما استفاده کنید. دمای سیال باید از Tin بالاتر برود، اما شکل دقیق میدان فقط پس از اجرای واقعی مشخص می‌شود. یک Plot جدا برای سرعت و یک Plot برای دما بسازید تا واحد و راهنمای رنگ هر شکل روشن باشد. دما را در نمایش Celsius می‌توان نشان داد، ولی مراقب تفاوت K و °C در محاسبات باشید.

۱۰. کنترل کمّی نتیجه؛ فقط به تصویر رنگی اعتماد نکنیم

۱۰-۱. افت فشار مرجع

برای جریان پوازی کاملاً توسعه‌یافته بین دو صفحه تخت، گرادیان فشار و افت فشار مرجع از رابطه (6) به دست می‌آیند. با اعداد آموزشی، افت فشار در طول 0.8L حدود 4.8 Pa است. مقاطع x1=0.1L و x2=0.9L را تعریف و فشار متوسط آن‌ها را مقایسه کنید. فشار دو نقطه منفرد نزدیک ورودی یا گوشه خروجی الزاماً معیار مقاومی نیست.

(6)

۱۰-۲. موازنه انرژی

توان گرمای ورودی به‌ازای عمق واحد در مدل دوبعدی مطابق رابطه (7) است. دبی جرمی به‌ازای عمق واحد از رابطه (8) به دست می‌آید. در تقریب صرف‌نظر از رسانش محوری مرزی، افزایش دمای توده از رابطه (9) برآورد می‌شود. این فقط پاسخ مرجع تقریبی است و نتیجه محاسبه COMSOL نیست؛ برای ارزیابی دقیق باید شار رسانشی ورودی/خروجی هم در موازنه انرژی منظور شود.

	(7)
	(8)
	(9)

دمای توده خروجی باید با وزن‌دهی بر اساس دبی جرمی محاسبه شود، نه با میانگین ساده دما روی لبه. در Component 1 → Definitions → Component Couplings یک Integration به نام intout روی ضلع خروجی سیال تعریف کنید. برای خروجی عمود بر محور x و جریان رو به راست، از رابطه (10)  استفاده   کنید؛ با چگالی ثابت، rho_f در صورت و مخرج حذف می‌شود. در Results → Derived Values → Global Evaluation آن را محاسبه کنید. از وجود متغیرهای u و T در همان Dataset، علامت سرعت خروجی و واحد K مطمئن شوید. برای محاسبه موازنه انرژی دقیق، شارهای رسانشی و همرفتی همه مرزهای باز را نیز در نظر بگیرید.

(10)

۱۱. مطالعه پارامتری و تفسیر نتایج

در Study راست‌کلیک و Parametric Sweep را اضافه کنید. برای Um سه مقدار 0.005[m/s]، 0.01[m/s] و 0.02[m/s] قرار دهید؛ سایر پارامترها و شبکه را یکسان نگه دارید. Compute را اجرا کنید و برای هر سرعت افت فشار و Tb_out را یادداشت کنید. انتظار تحلیلی مدل خواص ثابت این است که افت فشار با Um افزایش یابد و در تقریب موازنه انرژی ساده، ΔTb با Um کاهش پیدا کند؛ اما بیشینه دمای جامد باید از خروجی واقعی مدل ارزیابی شود.

۱۲. اگر نتیجه اشتباه شد، از کجا شروع کنیم؟

علامت مشکل	اولین بررسی عملی
سرعت در ناحیه جامد دیده می‌شود	Domain Selection در Laminar Flow
سرعت ورودی علامت منفی دارد	جهت محور x و مؤلفه ورودی
دما با گرم‌کردن کف کاهش یافته	علامت و مکان Heat Flux
افت فشار بسیار متفاوت است	واحد H، ویسکوزیته، وضعیت جریان توسعه‌یافته
حل متوقف می‌شود	واحدها، مرز فشار، ناسازگاری Selection، شبکه
افزایش دمای خروجی غیرمنطقی است	وزن‌دهی Tb_out و موازنه انرژی
پلات سفید یا خالی است	Dataset و Expression انتخاب‌شده

۱۳. گسترش مدل به سه‌بعد

پس از موفقیت مدل دوبعدی، پروژه جدیدی در 3D بسازید و هندسه سه‌ناحیه‌ای را با عرض معین W به‌صورت Block ایجاد کنید. اکنون دیواره‌های جانبی نیز وجود دارند و شرط No slip روی آن‌ها پروفیل سرعت را عوض می‌کند. بنابراین رابطه سهموی دوبعدی و تقریب قطر هیدرولیکی از رابطه (11) تنها زمانی مناسب‌اند که عرض نسبت به ارتفاع بسیار بزرگ باشد. ابتدا نتایج مرکزی را با مدل دوبعدی مقایسه کنید و سپس حساسیت به عرض و تعداد المان در بعد سوم را بسنجید؛ هیچ مقدار سه‌بعدی بدون اجرای واقعی گزارش نشود.

(11)

۱۴. چک‌لیست پایان پروژه

پیش از انتشار یا استفاده از نتایج، مطمئن شوید مدل دارای هندسه و واحدهای مشخص، جدول پارامترها، انتخاب‌های صحیح مرز و دامنه، خواص دو ماده، رابط‌های Laminar Flow و Heat Transfer، کوپل Nonisothermal Flow، مشِ ساخته‌شده، نمودار سرعت و دما، بررسی جواب تحلیلی فشار و موازنه انرژی است. تمام نمودارهای شبیه‌سازی را فقط پس از اجرای واقعی COMSOL ذخیره کنید و در زیرنویس شکل، وضعیت «خروجی نرم‌افزار» یا «شماتیک آموزشی» را روشن بنویسید. مدل را با پسوند .mph ذخیره کنید.

پانویس اصطلاحات تخصصی (معادل انگلیسی)

شبیه‌سازی چندفیزیکی: Multiphysics Simulation؛

 انتقال حرارت مزدوج: Conjugate Heat Transfer؛

جریان آرام: Laminar Flow؛

 شرط مرزی: Boundary Condition؛

 شبکه‌بندی: Meshing؛

 استقلال از شبکه: Mesh Independence؛

حل پایا: Stationary Study؛

 شار گرمایی: Heat Flux؛

دمای توده: Bulk Temperature؛

 سرعت کاملاً توسعه‌یافته: Fully Developed Flow؛

انتخاب دامنه: Domain Selection؛

 انتگرال‌گیری مرزی: Boundary Integration؛

حلگر: Solver؛

 موازنه انرژی: Energy Balance.

منابع و مستندات راهنما

[1] COMSOL Multiphysics 5.6, Introduction to the Heat Transfer Module, official documentation: https://doc.comsol.com/5.6/doc/com.comsol.help.heat/IntroductionToHeatTransferModule.pdf

[2] COMSOL Multiphysics 5.6, Thermal Modeling of a Microchannel Heat Sink, official example: https://doc.comsol.com/5.6/doc/com.comsol.help.models.heat.microchannel_heat_sink/models.heat.microchannel_heat_sink.pdf

[3] COMSOL Multiphysics 5.6 Reference Manual، بخش Heat Flux و General inward heat flux: https://doc.comsol.com/5.6/doc/com.comsol.help.comsol/COMSOL_ReferenceManual.pdf

[4] COMSOL Multiphysics 5.6، Fluid–Structure Interaction example؛ نمونه رسمیِ ورود پروفیل سهموی سرعت: https://doc.comsol.com/5.6/doc/com.comsol.help.models.mems.fluid_structure_interaction/models.mems.fluid_structure_interaction.pdf

[5] روابط مرجع جریان پوازی بین صفحات موازی و موازنه انرژی یک‌بعدی؛ در این مقاله به‌صورت روابط تحلیلی آموزشی استفاده شده‌اند و جایگزین اعتبارسنجی عددی نیستند.