مقایسه کامسول (COMSOL) با نرم افزارهای شبیه‌ساز
انتخاب ابزار شبیه‌سازی چندفیزیکی به معادلات حاکم، روابط سازنده، شدت بازخورد میدان‌ها و شواهد اعتبار مدل وابسته است. این مقاله مروری، کامسول را در کنار حلگرهای صنعتی، بسته‌های متن‌باز، چارچوب‌های توسعه عددی و کتابخانه‌های کوپل بررسی می‌کند. مقایسه بر پایه دامنه کاربرد و معماری مدل‌سازی انجام می‌شود و شامل آزمون کنترل‌شده سرعت یا دقت محصولات نیست. انتقال حرارت مزدوج، اندرکنش سیال و سازه، رئولوژی و الکتروشیمی به‌عنوان مسائل نمونه بررسی می‌شوند. تحلیل نشان می‌دهد که آزادی تعریف معادله و ترکیب میدان‌ها، دلیل مهمی برای بررسی کامسول در پژوهش است؛ در مقابل، حلگرهای تخصصی و تجربه تیم می‌توانند در پروژه‌های صنعتی تعیین‌کننده باشند. انتخاب نهایی باید با یک مسئله نماینده، مطالعه شبکه و زمان و مقایسه با شواهد مستقل پشتیبانی شود. یک پروژه آموزشی کانال گرم‌شونده نیز برای اتصال مفاهیم به اجرای عملی، مقایسه با پاسخ تحلیلی و کنترل موازنه انرژی ارائه می‌شود.
1- مقدمه و دامنه مطالعه
این راهنما چه پرسشی را پاسخ می‌دهد
برای انتخاب نرم‌افزار، ابتدا باید بدانیم کدام پدیده‌ها تعیین‌کننده‌اند، چه کمیتی را می‌خواهیم پیش‌بینی کنیم و چه دقتی برای تصمیم مهندسی لازم است. پرسش «کامسول بهتر است یا انسیس» بدون معرفی پروژه پاسخ علمی یکتایی ندارد. همین ابزار ممکن است برای یک حسگر چندفیزیکی انتخاب مناسبی باشد و برای تحلیل برخورد یک سازه، مسیر کم‌دردسرتری در نرم‌افزار دیگری وجود داشته باشد. این متن ابتدا مبانی را توضیح می‌دهد، سپس ابزارها را بر اساس معماری و دامنه کاربرد مقایسه می‌کند و در پایان روش انتخاب و آزمون مدل را ارائه می‌دهد. پیشنهادهای انتخاب ابزار، تحلیل مهندسی نویسنده بر پایه قابلیت‌های مستند هستند؛ حاصل یک آزمون سرعت یا دقت کنترل‌شده بین محصولات نیستند. بنابراین از امتیازدهی ساختگی و معرفی یک برنده همیشگی پرهیز شده است.
ابزارهایی که بررسی می‌شوند
COMSOL Multiphysics ؛ اکوسیستم ANSYS شامل Fluent ، Mechanical ، Maxwell و System Coupling ؛ Simcenter STAR-CCM+ ؛ Abaqus ؛ CST Studio Suite ؛ مجموعه Altair شامل AcuSolve ، OptiStruct ، Flux و SimLab ؛ OpenFOAM ؛ Elmer ؛ MOOSE ؛ FEniCSx ؛ deal.II ؛ SU2 ؛ SimScale و کتابخانه کوپل preCICE بررسی می‌شوند. ابزارهای تخصصی الکترومغناطیس و فوتونیک نیز در زمینه کاربرد خود معرفی می‌شوند. این انتخاب نماینده خانواده‌های اصلی است، نه فهرست تمام محصولات بازار. نرم‌افزار آماده، مجموعه حلگر، کتابخانه عددی و ابزار کوپل در جدول‌ها از یکدیگر تفکیک شده‌اند. حضور یک نام در دسته «چندفیزیک» به معنی پشتیبانی آماده از هر ترکیب دلخواه نیست.
مسیر پیشنهادی برای خواننده
مبتدی: بخش‌های مبانی، جریان کار و سناریوهای انتخاب. کاربر متوسط: معرفی ابزارها، جدول‌های مقایسه و اعتبارسنجی. پژوهشگر پیشرفته: معماری کوپل، خطاهای انتقال، رئولوژی، انتخاب پیش‌شرط‌گذار و طراحی آزمون مقایسه. قابلیت‌ها بر پایه مستندات رسمی بررسی‌شده در ۱۳ اکتبر ۲۰۲۶بیان شده‌اند. انتخاب نهایی باید با نسخه نصب‌شده، مجوز ماژول‌ها، محدودیت محاسباتی و مستندات همان نسخه تطبیق داده شود[۲۵–۱]. بخش پروژه با مستندات رسمی نسخه ۶٫۴ در ۱۵ اکتبر ۲۰۲۶ تطبیق داده شده است.
2- مبانی شبیه‌سازی چندفیزیکی
فیزیک مدل با نام نرم‌افزار فرق دارد
در مدل انتقال حرارت، مجهول اصلی می‌تواند دما باشد؛ در مکانیک جامدات، جابه‌جایی؛ و در جریان سیال، سرعت و فشار. هر یک به معادلات حاکم، خواص ماده، شرایط مرزی و شرایط اولیه نیاز دارند. مسئله چندفیزیکی زمانی ایجاد می‌شود که پاسخ یک میدان بر معادله، ضرایب، هندسه یا شرایط میدان دیگر اثر بگذارد. مثلاً جریان هوا گرما را از یک قطعه می‌گیرد. اگر تغییر دما چگالی هوا را تغییر دهد و شناوری مهم شود، دما نیز جریان را تغییر می‌دهد. در این حالت رابطه دوطرفه است. اگر اثر دما بر جریان ناچیز باشد، شاید حل جریان و سپس انتقال حرارت کافی باشد. اضافه‌کردن فیزیک دوم بدون وجود رابطه اثرگذاری، به‌تنهایی مدل را علمی‌تر نمی‌کند.
چندفیزیک با چندفازی یکی نیست
جریان آب و هوا یک مسئله چندفازی است. اگر مدل فقط حرکت دو فاز را حل کند، ممکن است هنوز تنها در حوزه مکانیک سیالات باشد. برعکس، گرمایش الکتریکی یک قطعه تک‌فاز، با کوپل میدان الکتریکی و دما، چندفیزیکی است. مدل چندمقیاسی نیز موضوع دیگری دارد: ارتباط توصیف میکروسکوپی با پاسخ در مقیاس بزرگ‌تر. یک پروژه می‌تواند هر سه ویژگی را داشته باشد.
نخستین تصمیم علمی
پیش از انتخاب ابزار، خروجی موردنظر را مشخص کنید: بیشینه دما، افت فشار، تنش، نرخ تبدیل، امپدانس، فرکانس تشدید یا ویسکوزیته مؤثر. سپس بپرسید حذف کدام اثر می‌تواند این خروجی را به‌طور معنادار تغییر دهد«. معنادار» باید نسبت به نیاز پروژه، عدم‌قطعیت داده و خطای مجاز تعریف شود، نه با یک آستانه عمومی برای همه مسائل. برای نمونه، در خنک‌کاری یک مدار، تحلیل تغییر شکل ممکن است ضروری نباشد؛ اما اگر تغییر شکل فاصله یک کانال را عوض کند، همان میدان مکانیکی می‌تواند مسیر انتقال حرارت را تغییر دهد. کیفیت مدل به همین تشخیص وابسته است[۱، ۲، ۳].
3-معادلات حاکم و روابط سازنده
میدان‌ها چگونه به هم متصل می‌شوند
برای یک سیال نیوتنی تراکم‌ناپذیر، معادله پیوستگی بیان می‌کند که واگرایی سرعت صفر است. معادله مومنتوم تغییرات سرعت را به فشار، تنش ویسکوز و نیروهای حجمی مرتبط می‌کند. معادله انرژی نیز دما را به رسانش، همرفت و تولید گرما متصل می‌کند. در یک مدل غیرنیوتنی، رابطه تنش و نرخ تغییرشکل باید متناسب با ماده جایگزین شود.
در روابط بالا u سرعت p ، فشار tau ، تنش انحرافی rho ، چگالی Cp ، ظرفیت گرمایی ویژه k ، رسانندگی و Q تولید گرمای حجمی است. این روابط، یک چارچوب آموزشی هستند؛ انتخاب فرم دقیق برای جریان تراکم‌پذیر، چندفازی یا ماده دارای حافظه به فرضیات دیگری نیاز دارد.
یک مثال روشن از کوپل
در گرمایش ژولی، جریان الکتریکی گرما تولید می‌کند. اگر رسانندگی الکتریکی تابع دما باشد، حل میدان الکتریکی نیز به دما وابسته می‌شود. در حالت رسانندگی ایزوتروپ، چگالی تولید گرما از ضرب داخلی چگالی جریان و میدان الکتریکی به‌دست می‌آید. برای ماده ناهمسانگرد، جایگزین‌کردن یک رسانندگی اسکالر بدون بررسی، ممکن است مدل را تغییر دهد.
رابطه سازنده بخش اصلی انتخاب مدل است
قانون ماده از معادلات بقا جداست. حل «Navier–Stokes» به‌تنهایی به معنی مدل‌کردن ویسکوالاستیسیته یا تیکسوتروپی نیست. همین‌طور وجود یک رابط سازه‌ای، مشخص نمی‌کند پلاستیسیته، آسیب یا رفتار وابسته به زمان با چه فرض‌هایی تعریف شده‌اند. یک نرم‌افزار تنها معادلاتی را حل می‌کند که مدل‌ساز به آن داده یا در رابط انتخاب کرده است. در نتیجه، پیشرفته‌ترین بخش کار اغلب انتخاب رابطه سازنده و شرایط مرزی مناسب است. توسعه یک رابطه جدید در محیط معادله‌محور ممکن است ارزش بیشتری از استفاده از ده‌ها ماژول آماده داشته باشد[۲، ۸، ۱۶].
4- معماری کوپل و الگوریتم حل
یکطرفه و دوطرفه توصیف رابطه فیزیکی هستند
در کوپل یکطرفه، نتیجه مسئله اول ورودی مسئله دوم است و پاسخ دوم به اول بازنمی‌گردد. در کوپل دوطرفه، این بازخورد لحاظ می‌شود. مثلاً فشار سیال می‌تواند سازه را خم کند و تغییر شکل سازه، هندسه جریان را تغییر دهد. اگر اثر دوم حذف شود، مدل یکطرفه است؛ حتی اگر هر دو حلگر بسیار پیشرفته باشند.
یکپارچه و تفکیکی توصیف الگوریتم هستند
در رویکرد یکپارچه یا، monolithic مجهولات میدان‌ها در یک دستگاه جبری مشترک حل می‌شوند. در رویکرد تفکیکی یا، partitioned حلگرهای جدا داده تبادل می‌کنند. داخل یک نرم‌افزار واحد نیز می‌توان گروه‌های مجهول را به‌صورت segregated حل کرد. بنابراین «یک رابط کاربری» با «یک دستگاه معادلات» یکسان نیست. تفکیک حل، معمولاً استفاده از حلگرهای تخصصی را آسان‌تر می‌کند و می‌تواند نیاز حافظه هر زیرمسئله را کاهش دهد. در مقابل، همگرایی تکرارهای بین میدان‌ها و خطای انتقال داده باید کنترل شود. حل یکپارچه می‌تواند کوپل‌های قوی را مستقیم‌تر بازنمایی کند، اما دستگاه بزرگ‌تر و انتخاب پیش‌شرط‌گذار مناسب را به همراه دارد. هیچ‌کدام تضمین عمومی برای سریع‌تر یا دقیق‌تر بودن نیستند.
هم‌شبیه‌سازی و چند زمان‌بندی
در، co-simulation ابزارها در یک جریان کار هماهنگ اجرا می‌شوند. تبادل فقط یک فایل نتیجه در پایان، با تبادل تکراری در هر گام زمانی تفاوت دارد. اگر گرما طی ثانیه‌ها تغییر کند ولی ارتعاش در میلی‌ثانیه‌ها رخ دهد، انتخاب زمان تبادل اهمیت پیدا می‌کند. کوچک‌کردن گام زمانی حلگر سیال، خطای ناشی از تبادل دیرهنگام داده بین دو حلگر را لزومًا حذف نمی‌کند. در پروژه باید مشخص باشد چه داده‌ای، روی چه سطحی، با چه فاصله زمانی و با چه معیار همگرایی منتقل می‌شود. مستندات COMSOL ، System Coupling و preCICE نمونه‌های این معماری‌ها را توضیح می‌دهند. [۲، ۴، ۱۵[
5- فرایند تدوین مدل عددی
جریان کار مشترک همه ابزارها
هدف و خروجی را تعریف کنید، هندسه مؤثر را انتخاب کنید، خواص ماده و فرضیات را بنویسید، شرایط مرزی را تعیین کنید و سپس گسسته‌سازی و حلگر را انتخاب کنید. پس‌پردازش پایان کار نیست؛ نتیجه باید از نظر عددی بررسی و با شواهد مستقل مقایسه شود. تفاوت نرم‌افزارها در میزان خودکارسازی و آزادی هر مرحله است. ساده‌سازی هندسه باید نسبت به خروجی توجیه شود. حذف یک لبه کوچک ممکن است برای افت فشار بی‌اثر باشد، اما همان لبه در بیشینه تنش اهمیت داشته باشد. تقارن نیز زمانی مجاز است که هندسه، بارگذاری، خواص و پاسخ موردانتظار با آن سازگار باشند. یک مسئله متقارن می‌تواند ناپایداری نامتقارن داشته باشد؛ تحمیل تقارن ممکن است آن پاسخ را حذف کند.
داده ماده باید در محدوده درست باشد
رسانندگی در یک دمای ثابت، برای بازه دمایی گسترده کافی نیست. پارامتر ویسکوزیته استخراج‌شده از یک نرخ برش، قانون عمومی جریان نیست. برای هر داده باید واحد، روش آزمون، بازه دما و نرخ تغییرشکل و عدم‌قطعیت ثبت شود. داده‌های واردشده از منابع مختلف باید پیش از استفاده سازگار شوند.
مدل ساده یک نقطه کنترل است
پیش از اجرای هندسه پیچیده، یک مسئله ساده با حل تحلیلی یا پاسخ مرجع بسازید: رسانش یک‌بعدی، جریان پوازی، خمش تیر یا یک مسئله انتشار. این آزمون نمی‌تواند مدل نهایی را به‌تنهایی اعتبارسنجی کند، اما اشتباه در واحدها، قرارداد علامت، شرایط مرزی و پیاده‌سازی معادله را آشکار می‌کند.
6- گسسته‌سازی و انتخاب شبکه
FEM و FVM چه تفاوتی دارند
روش اجزای محدود معمولاً از فرم ضعیف معادلات و توابع تقریب استفاده می‌کند؛ روش حجم محدود معادلات بقا را روی حجم‌های کنترل گسسته می‌کند FEM. در بسیاری از مسائل سازه، میدان و معادلات سفارشی کاربرد دارد و FVM در بسیاری از حلگرهای جریان صنعتی رایج است. این توصیف به‌معنی انحصار نیست FEM: نیز می‌تواند جریان سیال را حل کند و FVM نیز برای بعضی مسائل جامد و چندفیزیک استفاده می‌شود. انتخاب روش تنها یکی از اجزای مسئله است. مرتبه تقریب، کیفیت شبکه، پایدارسازی، گسسته‌سازی زمان و نوع حلگر نیز مهم‌اند. نمی‌توان از نام FEM یا FVM نتیجه گرفت که کدام محصول برای یک مسئله همیشه دقیق‌تر است. برای ادعای دقت، به آزمونی با فیزیک، داده، شرایط مرزی و معیار خطای هم‌تراز نیاز داریم[۳، ۱۰، ۱۳].
مقیاس بدون بعد چه کمکی می‌کند
عدد رینولدز نسبت اثر اینرسی به ویسکوزیته را نشان می‌دهد. عدد پکله حرارتی نسبت انتقال همرفتی به پخش حرارتی است. عدد وایسنبرگ نسبت زمان آسایش ماده به زمان تغییرشکل است. این اعداد مشخص می‌کنند کدام اثر اهمیت دارد و کجا حل ممکن است سخت شود؛ اما عددی مانند Wi به‌تنهایی مشخص نمی‌کند کدام رابطه سازنده درست است.
در این روابط U و L مقیاس سرعت و طول alpha ، پخشندگی حرارتی و lambda زمان آسایش‌اند. ویسکوزیته یک سیال غیرنیوتنی ثابت نیست؛ برای تعریف Re باید مقیاس ویسکوزیته و شرایط ارزیابی آن اعلام شود. همچنین معتبر بودن پیوستار در مقیاس‌های بسیار کوچک باید مستقل بررسی شود.
شبکه مناسب برای هر میدان
لایه مرزی سرعت، گرادیان دما، شکاف تماس و طول موج صوتی نیاز یکسانی ندارند. شبکه‌ای که دما را خوب حل می‌کند ممکن است برای تنش یا موج کافی نباشد. در مدل چندفیزیکی، همگرایی باید برای هر خروجی مهم و برای کمیت تبادل‌شده در مرز بررسی شود.
نویسنده:
مهندس سارینا داورزنی
کارشناسی ارشد مهندسی پلیمر دانشگاه تهران
(برای مطالعه بیشتر روی نام یا تصویر ایشان کلیک کنید)