دانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Enhancing Defect Detection in Additively Manufactured PLA and ABS Polymer Parts Using PCA and PPT Image Processing in Active Thermographyبهبود آشکارسازی عیوب در قطعات پلیمری PLA و ABS تولیدشده به روش افزایشی با ترکیب پردازش تصویر PCA و PPT در ترموگرافی فعال7597712787810.48311/mme.2025.117176.82871FAزهراحسنیگروه مهندسی مکانیک، دانشکدگان فنی، دانشگاه تهران، تهران، ایرانمحمدرضافراهانیگروه مهندسی مکانیک، دانشکدگان فنی، دانشگاه تهران، تهران، ایران0000-0001-5515-4283Journal Article20251029<span style="font-size: 9.0pt; mso-bidi-font-size: 9.5pt; line-height: 97%; font-family: 'Calibri',sans-serif; mso-fareast-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">This study investigates the performance of active infrared thermography for detecting surface and subsurface defects in two commonly used polymers, PLA and ABS. To enhance detection accuracy and image contrast, two advanced image-processing techniques—Principal Component Analysis (PCA) and Pulsed Phase Thermography (PPT)—were applied individually and in combination to the thermal image sequences. Results show that PCA efficiently reduces noise and emphasizes dominant thermal patterns, providing superior identification of shallow defects. In contrast, PPT, through frequency-domain analysis, exhibits lower sensitivity to non-uniform heating and emissivity variations, enabling more reliable detection of deeper defects. Comparative thermal behavior analysis revealed that PLA, due to its higher thermal conductivity, identifies surface defects more rapidly and with greater contrast, whereas ABS, with its higher specific heat capacity, offers improved stability and precision in subsurface defect detection. Moreover, the hybrid application of PCA and PPT proved sequence-dependent: in PLA, the PCA→PPT approach yielded a balanced visualization of both surface and subsurface anomalies, while in ABS, the PPT→PCA sequence significantly enhanced subsurface defect contrast. These findings demonstrate that combining complementary thermographic processing methods can overcome the limitations of single techniques and, when tailored to material-specific thermal properties, provide a robust framework for non-destructive evaluation of additively manufactured polymer components</span>در این پژوهش کارایی ترموگرافی فعال مادونقرمز در شناسایی عیوب سطحی و نیمهسطحی در دو پلیمر پرکاربرد PLA و ABS مورد بررسی قرار گرفت. بهمنظور ارتقای وضوح و افزایش دقت آشکارسازی، دو روش پردازش تصویر تحلیل مؤلفههای اصلی (PCA) و ترموگرافی پالسی فازی (PPT) بهطور جداگانه و ترکیبی بر دادههای حرارتی اعمال شدند. نتایج نشان داد که PCA در حذف نویز و برجستهسازی مؤلفههای غالب عملکرد مؤثری دارد و عیوب کمعمق را با وضوح بالایی آشکار میکند، در حالی که PPT با تحلیل دادهها در حوزه فرکانس، نسبت به ناهمگنی گرمایش و تغییرات ضریب گسیل حساسیت کمتری داشته و عیوب نیمهسطحی را با دقت بیشتری نمایان میسازد. مقایسه رفتار حرارتی دو ماده نشان داد که PLA به دلیل هدایت حرارتی بالاتر، عیوب سطحی را سریعتر و با کنتراست بالاتر آشکار میکند، اما ABS به دلیل ظرفیت گرمایی ویژهی بیشتر، در آشکارسازی عیوب نیمهسطحی پایدارتر و دقیقتر عمل میکند. تحلیل ترکیبی دو روش نیز نشان داد که ترتیب اعمال آنها نقش تعیینکنندهای دارد؛ بهطوریکه در PLA رویکرد پردازش تصویر PCA و بعد PPT متوازنتر بوده و امکان آشکارسازی همزمان عیوب سطحی و نیمهسطحی را فراهم ساخته است، در حالی که در ABS رویکرد PPT و بعد PCA با تقویت دادههای فرکانسی، وضوح عیوب نیمهسطحی را به شکل قابلملاحظهای افزایش داده است. این یافتهها بیانگر آن است که ترکیب روشهای پردازش تصویر میتواند محدودیتهای روشهای منفرد را جبران کرده و با تطبیق ترتیب اعمال آنها با ویژگیهای حرارتی ماده، ابزاری کارآمد برای بازرسی غیرمخرب قطعات تولید افزایشی فراهم سازد.https://mme.modares.ac.ir/article_27878_4a2d84630e650155dd3f0d107c9794a9.pdfدانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Numerical Simulation of Flexible Composite Skin Bending in Morphing Wings Using the Radial Point Interpolation Mesh free Method: Comparison with Experimental Resultsحل عددی خمش پوسته کامپوزیتی انعطاف پذیر مورد استفاده در بال مورفینگ به روش بدون شبکه ی درونیابی نقاط شعاعی و مقایسه با نتایج تجربی7737832787110.48311/mme.2025.96923.0FAمرتضیشکیباسرشتدانشگاه صنعتی مالک اشتر، مجتمع دانشگاهی هوافضا، تهران، ایرانمحمودذبیح پوردانشگاه صنعتی مالک اشتر، مجتمع دانشگاهی هوافضا، تهران، ایران0000-0002-2648-8070کرامتملک زادهدانشگاه صنعتی مالک اشتر، مجتمع دانشگاهی هوافضا، تهران، ایرانجعفراسکندری جمدانشگاه صنعتی مالک اشتر، مجتمع دانشگاهی هوافضا، تهران، ایرانJournal Article20250923<span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Calibri',sans-serif; mso-fareast-font-family: Calibri; color: black; mso-themecolor: text1; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">In this study, the bending of a flexible composite shell (with an elastomeric matrix) used in a morphing wing skin is analyzed numerically, and the numerical results are compared with those of a reputable experimental study conducted on the bending behavior of this shell. First, the governing equations are derived using the total potential energy minimization method and the Euler–Lagrange equations, resulting in highly nonlinear differential equations. To discretize these equations, the strong-form meshfree Radial Point Interpolation Method(RPIM) is employed, and the Newton–Raphson method is used to obtain the numerical solution of the nonlinear equations. The boundary conditions are considered clamped on all four edges. In this analysis, the effects of pretension forces as well as out-of-plane loads on the maximum displacement of the composite shell are investigated. Comparison of the meshfree numerical results with the experimental data shows that the proposed numerical solution has good accuracy and exhibits only a small deviation from the experimental results</span>در این مقاله خمش یک پوسته کامپوزتی انعطافپذیر(ماتریس الاستومر) مورد استفاده در پوسته بال مورفینگ به روش عددی مورد آنالیز قرار گرفته است و نتایج حاصل از حل عددی با یک مقاله معتبر که بصورت تجربی برروی خمش این پوسته تحقیقی انجام شده، مقایسه شدهاست. در ابتدا معادلات حاکم با استفاده از روش حداقل انرژی پتانسیل کل و معادلات اویلر-لاگرانژ استخراج شده است که شامل معادلات دیفرانسیل بشدت غیرخطی میباشند، جهت گسسته سازی این معادلات از روش بدون شبکه درونیابی نقاط شعاعی به فرم قوی استفاده شده است و همچنین از روش نیوتن-رافسون جهت حل عددی معادلات غیر خطی بهره بردهشده است. شرایط مرزی بصورت گیردار در چهار طرف درنظر گرفته شده است. در این تحلیل تاثیر نیروهای پیش کشش و همچنین بارهای خارج ازصفحه در میزان جابجایی حداکثر پوسته کامپوزیتی بررسی شده است. مقایسه نتایج حل به روش عددی بدون شبکه با نتایج تجربی نشان داد که حل به روش مذکور از دقت خوبی برخوردار است واختلاف اندکی با نتایج تجربی دارد.https://mme.modares.ac.ir/article_27871_7cefdcbff0ddf29c3bbaeec362015c30.pdfدانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Parametric Study on the Influence of Moving Surface Length and Speed on Performance Enhancement of Offshore Horizontal-Axis Wind Turbineمطالعه پارامتریک اثر طول و سرعت صفحه متحرک بر بهبود عملکرد توربین بادی محور افقی فراساحلی7858002784810.48311/mme.2025.96919.0FAمجتبیطحانیدانشکدگان علوم و فناوریهای بینرشتهای، دانشگاه تهران، تهران، ایرانعلیرضاشهرابی فراهانیدانشکدگان علوم و فناوریهای بینرشتهای، دانشگاه تهران، تهران، ایراندارینرضویفردانشکدگان علوم و فناوریهای بینرشتهای، دانشگاه تهران، تهران، ایرانفاطمهفرامرزلودانشکدگان علوم و فناوریهای بینرشتهای، دانشگاه تهران، تهران، ایرانامیرحسینحسیندانشکدگان علوم و فناوریهای بینرشتهای، دانشگاه تهران، تهران، ایرانJournal Article20250820<span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Calibri',sans-serif; mso-fareast-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">The implementation of a moving surface in the aerodynamic profile of an airfoil is recognized as an effective active flow control mechanism. In this study, the effect of this mechanism on the aerodynamic and mechanical performance of an offshore horizontal-axis wind turbine blade at a Reynolds number of 7.5×10<sup>5 </sup>was investigated. The S809 airfoil was selected as the blade section, and part of its suction surface was replaced with a moving surface. Simulations were performed using the commercial software ANSYS Fluent with the k-ω SST turbulence model, and the results were validated against available experimental data. To evaluate the influence of the moving surface length and speed, relative lengths (l) ranging from 0.05 to 0.55 chord and speed ratios of (k) 1,3 and 5 times the free-stream velocity were examined at a fixed installation location of s = 0.01 chord. The moving surface delayed stall angle by up to 80%, such that the drop in the lift curve shifted from an AoA of 8</span><span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Cambria Math',serif; mso-fareast-font-family: Calibri; mso-bidi-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">°</span><span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Calibri',sans-serif; mso-fareast-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;"> to 12</span><span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Cambria Math',serif; mso-fareast-font-family: Calibri; mso-bidi-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">°</span><span style="font-size: 12.0pt; mso-bidi-font-size: 11.0pt; line-height: 97%; font-family: 'Calibri',sans-serif; mso-fareast-font-family: Calibri; color: black; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">, and the mechanical performance improved by a factor of five to seven depending on the AoA. At higher angles of attack, shorter lengths maintained flow control while reducing power consumption. In the aerodynamic analysis using the Blade Element Momentum method, the maximum power coefficient at a TSR of 4 increased by up to 4.6 times, while <span style="mso-spacerun: yes;"> </span>the torque coefficient at the same point increased by 364%. The moving surface showed the greatest effect at low TSRs, improving both aerodynamic and mechanical performance</span>ایجاد صفحه متحرک در پروفیل آیرودینامیکی ایرفویل بهعنوان مکانیز مؤثر کنترل جریان فعال شناخته شده است. در این مطالعه اثر این مکانیزم بر عملکرد آیرودینامیکی و مکانیکی پرهٔ توربین بادی محور افقی فراساحلی در عدد رینولدز 750 هزار بررسی شده است. ایرفویل اس809 بهعنوان مقطع پره انتخاب و بخشی از سطح مکش آن با صفحهای متحرک جایگزین شده است. شبیهسازیها با استفاده از نرمافزار تجاری انسیس فلوئنت و مدل آشفتگی کی - امگا اساستی انجام شد و نتایج با دادههای تجربی موجود اعتبارسنجی گردید. بهمنظور ارزیابی تأثیر طول و سرعت صفحه متحرک، طولهای نسبی بین 0/05 تا 0/55 طول وتر ایرفویل و سرعتهای نسبی 1، 3 و 5 برابر سرعت جریان آزاد در موقعیت نصب 0/01 موردبررسی قرار گرفت. صفحه متحرک موجب تأخیر زاویه واماندگی تا حدود 80 درصد شد، بهطوری که افت شیب منحنی برآ از زاویه حمله 8 درجه به 12 درجه منتقل و عملکرد مکانیکی بسته به زاویه حمله بین پنج تا هفت برابر بهبود یافت. در زوایای حمله بالا، طولهای کوتاهتر ضمن حفظ خاصیت کنترل جریان، توان مصرفی را کاهش داد. در تحلیل آیرودینامیکی با روش مومنتوم المان پره بیشینه ضریب توان در نسبت سرعت نوک معادل 4 تا مقدار 4/6 برابر افزایش یافت و ضریب گشتاور نیز در همان نقطه 364 درصد بیشتر شد. صفحه متحرک بیشترین اثر را در مقادیر پایین نسبت سرعت نوک پره داشته و هر دو عملکرد آیرودینامیکی و مکانیکی را بهبود بخشیده است.https://mme.modares.ac.ir/article_27848_5cc9bbd6c9c018e4642e8bf7292d65ce.pdfدانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Stress Field Prediction Using Conditional Generative Adversarial Networks and Image Processing in a Perforated Plate under Static Loadingپیشبینی میدان تنش با استفاده از شبکه مولد تخاصمی شرطی و پردازش تصویر در یک صفحه سوراخدار تحت بارگذاری استاتیکی8018182791310.48311/mme.2025.96917.0FAبهنامانبارلوییدانشگاه شهید بهشتی، دانشکده فناوری های نوین و مهندسی هوافضاملیکاملکی راددانشگاه شهید بهشتی، دانشکده فناوری های نوین و مهندسی هوافضاJournal Article20250818<span style="font-size: 8.0pt; mso-bidi-font-size: 8.5pt; line-height: 97%; mso-fareast-font-family: Calibri; mso-bidi-font-family: Calibri; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">The study of the mechanical properties and behavior of materials, as well as stress and strain fields, has been carried out using methods such as experiments, numerical methods, and precise mathematical solutions over the decades. In recent years, machine learning, and especially deep learning, have become one of the most commonly used methods in various engineering fields. One of its important applications is the prediction of material behavior in numerous structures. These methods have drawn significant attention due to their rapid execution, apposite accuracy, and implementation convenience, and are used as an alternative or supplementary tool for traditional analysis methods. Using the machine learning method, in case the problems are properly characterized, they can provide a much more powerful tool in a machine learning process compared to other tools. The purpose of this paper is to predict the stress field and maximum stress on a perforated plate under static loading using a deep learning method based on a conditional adversarial generative network (CGANs) and to quantify the results using an image processing method. Also, at the end, the numerical results obtained from this model are extracted and compared with the results attained from finite element analysis to evaluate the accuracy of the proposed model</span>بررسی خواص و رفتار مکانیکی مواد و میدان های تنش و کرنش از روشهایی مانند آزمایشهای تجربی، شبیهسازیهای عددی و حل دقیق ریاضی صورت می پذیرد. در سالهای اخیر، یادگیری ماشین و بهویژه یادگیری عمیق به یکی از روشهای پرکاربرد در حوزههای مختلف مهندسی تبدیل شدهاند. یکی از کاربردهای مهم آن، پیشبینی رفتار مواد در سازه های مختلف میباشد. این روشها به دلیل سرعت بالا، دقت مناسب و سهولت در پیادهسازی، توجه ویژهای را به خود جلب کردهاند و بهعنوان جایگزین یا ابزار کمکی برای روشهای سنتی تحلیل مورد استفاده قرار میگیرند و با استفاده از فرآیند یادگیری ماشین در صورت مشخصه سازی صحیح مسائل می تواندد در یک فرآیند یادگیری ماشین ابزاری به مراتب قوی تر از ابزارهای پیشین فراهم سازند. هدف از این مقاله، پیشبینی میدان تنش و حداکثر تنش وارده بر یک صفحه سوراخدار تحت بارگذاری استاتیکی با استفاده از روش یادگیری عمیق مبتنی بر شبکه مولد تخاصمی شرطی و کمی سازی نتایج با استفاده از روش پردازش تصویر میباشد. همچنین در پایان، نتایج عددی بهدستآمده از این مدل استخراج شده و با نتایج حاصل از تحلیل اجزاء محدود مقایسه میگردد تا میزان دقت مدل پیشنهادی مورد ارزیابی قرار گیرد.<br><br> https://mme.modares.ac.ir/article_27913_6f8c673720168eb65f33ef15f9802a0b.pdfدانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Stability of Rotating Beam Made of Ultra-light Metamaterials Subjected to Distributed Follower Forceپایداری تیر دوار ساخته شده از فرامواد فوق سبک تحت نیروی پیرو گسترده8198282787310.48311/mme.2025.96869.0FAمحمدتشکریاندانشکده مهندسی مکانیک، دانشگاه شیراز، شیراز، ایرانسید احمدفاضل زاده حقیقیدانشکده مهندسی مکانیک، دانشگاه شیراز، شیراز، ایراناسمعیلقوانلودانشکده مهندسی مکانیک، دانشگاه شیراز، شیراز، ایرانJournal Article20250607<span style="font-size: 9.0pt; line-height: 97%; mso-fareast-font-family: Calibri; mso-bidi-font-family: Calibri; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">Today, the implementation of new materials in aero-structures such as propellers is crucial in the aeronautical industry. Also, ultralight metamaterials with capabilities to reduce the weight of the structure have received much interest. Therefore, it is essential to investigate the stability of aero-structures under non-conservative forces such as the follower force, which is generally considered as the fluid-structure interaction effect. In this paper, the stability of the rotating Timoshenko beam subjected to distributed follower force with clamped-free boundary condition is investigated. The beam consists of a sandwich structure and the core is made of ultralight metamaterials. The governing equations of the beam are derived using Hamilton's principle. For numerical analysis, the differential quadrature method (DQM) is used and the results are compared with those obtained from previous studies. Four core structures including square, hexagonal, triangular, and mixed are studied and their effects on the critical force of the beam are studied. The results showed that among the studied metamaterial structures, hexagonal structure has a better performance than other structures. Also, based on the numerical results, the use of metamaterials as the core of the sandwich beam increases the critical force per unit weight. In addition, the increase in shear effects in short beams and the effect of beam rotation reduces the effects of the core structure in increasing the critical force compared to the non-rotating beam; therefore, at high rotational speeds and long beams, the critical force of the sandwich beam and the uniform beam are quite the same</span>امروزه بهکارگیری مواد نوین در سازههای هوایی مانند ملخها، امری حیاتی در صنعت هوافضا بهشمار میرود. از سوی دیگر، فرامواد فوقسبک با قابلیت کاهش وزن سازه، توجه زیادی را به خود جلب کردهاند. بنابراین، بررسی پایداری سازههای هوایی تحت نیروهای ناپایستار مانند نیروی پیرو که معمولاً بهعنوان مدل سادهسازی شده اثر اندرکنش سیال و سازه شناخته میشود ضروری است. در این مقاله، پایداری تیر دوار تیموشنکو تحت نیروی پیرو گسترده یکنواخت و شرط مرزی یکسرگیردار مورد بررسی قرار گرفته است. تیر دارای ساختار ساندویچی با هستهای از فرامواد فوقسبک است و معادلات حاکم بر تیر با استفاده از اصل همیلتون استخراج شدهاند. برای تحلیل عددی از روش تربیع دیفرانسیل (DQM) استفاده شده و نتایج با مطالعات پیشین مقایسه شدهاند. چهار نوع ساختار هسته شامل مربعی، ششضلعی، مثلثی و ترکیبی مورد بررسی قرار گرفته و تأثیر آنها بر نیروی بحرانی تیر تحلیل شده است. نتایج نشان داد که در میان ساختارهای فرامواد بررسیشده، ساختار ششضلعی عملکرد بهتری نسبت به سایر ساختارها دارد. همچنین، بر اساس نتایج بدست آمده، کاربرد فرامواد در هسته تیر ساندویچی باعث افزایش نیروی بحرانی به ازای وزن واحد میشود. علاوه بر آن، افزایش اثرات برشی در تیرهای کوتاه و نیز اثر دوران تیر، تأثیر ساختار هسته را در افزایش نیروی بحرانی نسبت به تیر غیر دوار کاهش میدهد؛ در نتیجه، در هر دو حالت سرعتهای دورانی بالا و طول تیر زیاد، نیروی بحرانی تیر ساندویچی و تیر یکنواخت تقریباً یکسان میشود.https://mme.modares.ac.ir/article_27873_f2f25c1b2ea4636e7cf2190a36677f98.pdfدانشگاه تربیت مدرسمهندسی مکانیک مدرس2476-6909251220251122Investigating Transport and Deposition of Nanoparticles in Voronoi Porous Foams with Various Pore Densities using Euler–Lagrange Method with Pore-Scale Perspectiveمطالعه انتقال و رسوب نانوذرات در جریان نانوسیال در فومهای متخلخل ورونویی با چگالی حفرات متفاوت به روش اویلر-لاگرانژ با رویکرد مقیاس حفره8298412793610.48311/mme.2025.116918.82866FAمحمددهقاندانشکده مهندسی مکانیک، دانشگاه علم و صنعت ایران، تهران، ایرانمجیدسیاوشیدانشکده مهندسی مکانیک، دانشگاه علم و صنعت ایران، تهران، ایرانJournal Article20251019<span style="font-size: 8.0pt; line-height: 115%; font-family: 'Times New Roman',serif; mso-fareast-font-family: Calibri; mso-ansi-language: EN-US; mso-fareast-language: EN-US; mso-bidi-language: AR-SA;">Pore-scale simulations, by explicitly incorporating the geometrical features of the porous medium, enable accurate examination of fluid flow and particle transport within the foam. The present study numerically investigates nanoparticle transport and deposition in open-cell metallic foams with Voronoi-based geometries and a fixed porosity of 80%. Foam structures were reconstructed using the Laguerre–Voronoi algorithm, and simulations were performed with an Eulerian–Lagrangian framework implemented in the open-source software OpenFOAM. The model accounts for drag, Brownian motion, gravity, buoyancy, and Saffman lift forces, as well as van der Waals and electrostatic double-layer interactions to capture nanoparticle–wall adhesion. Results indicate that increasing the pore density from 30 to 60 pores per inch (PPI) enlarges the interfacial area while reducing flow velocity, thereby extending particle residence time and significantly enhancing deposition. Quantitatively, this increase in pore density leads to an approximate 26.5% rise in nanoparticle deposition ratio. Among the cases examined, the Voronoi foam with 80% porosity and 30 PPI exhibited the lowest deposition rate, whereas the foam with 80% porosity and 60 PPI showed the highest</span>شبیهسازی جریان در مقیاس حفره با وارد کردن ساختار هندسی دقیق محیط متخلخل، این امکان را فراهم میآورد تا بتوان جریان سیال در فوم متخلخل را با دقتی در مقیاس حفرات محیط متخلخل بررسی نمود. هدف از این مطالعه، بررسی تاثیر پارامترهای مختلف نانوسیال و محیط متخلخل بر روی میزان رسوب نانوذرات بر روی سطوح محیط متخلخل است. انتقال و رسوب نانوذرات در فومهای فلزی سلول باز با ساختار هندسی ورونوی و تخلخل ثابت %80 بهصورت عددی بررسی شده است. هندسه فومها با استفاده از الگوریتم لاگر- ورونوی تولید شده و شبیهسازیها با رویکرد اویلری–لاگرانژی در نرمافزار متنباز اپنفوم انجام شده است. در شبیهسازی، نیروهای پسا، براونی، گرانش، شناوری و برآ سافمن درنظر گرفته شده است. همچنین، نیروهای واندروالس و دولایهی الکترواستاتیکی برای شبیهسازی برهمکنش نانوذرات با دیواره محیط متخلخل لحاظ گردیده است. نتایج نشان میدهند که با افزایش چگالی حفره از 30 به 60 حفره در اینچ (PPI) سطح تماس افزایش یافته، سرعت جریان کاهش پیدا میکند و در نتیجه نرخ رسوبگذاری افزایش مییابد. به طور کمی، این افزایش چگالی حفره باعث رشد حدود 26.5% در نسبت رسوب نانوذرات شده است. در میان فومهای بررسیشده، فوم فلزی ورونوی با تخلخل 80% و چگالی حفره PPI 30 کمترین نرخ رسوبگذاری و فوم فلزی ورونوی با تخلخل 80% و چگالی حفره PPI 60 بیشترین نرخ رسوب را نشان دادهاند.https://mme.modares.ac.ir/article_27936_a73d731082fabe880672f123d4b3dcfe.pdf