Parametric Sensitivity Assessment of Slip Flow and Thermophoretic Effects on a Vertical Stretching Surface Using Response Surface Methodology

Technische Samenvatting: Parametrische Gevoeligheidsanalyse van Slipstroom en Thermoforetische Effecten op een Verticaal Rekendervlak met Gebruik van de Response Surface Methodologie

Probleemstelling
Deze studie onderzoekt de magnetohydrodynamische (MHD) grenslaagstroom van een elektrisch geleidende nanofluïde over een verticaal rekend vlak. Het onderzoek richt zich op de gekoppelde effecten van snelheidsslip, energyslip (thermische slip) en thermoforese op de momentum-, warmte- en massaoverdrachtskarakteristieken. Terwijl eerdere literatuur diverse slipcondities en nanofluïde gedragingen heeft verkend, richt dit werk zich specifiek op de systematische interactie van snelheidsslip-parameters ($L_1$) en energyslip-parameters ($L_2$) naast de thermoforese-parameter ($N_t$) om de thermische prestaties te optimaliseren. De onderliggende fysische verschijnselen omvatten geïnduceerde magnetische velden, thermische straling, Brownse beweging en chemische reacties binnen een stationair, tweedimensionaal laminair stromingsregime.

Methodologie
De studie maakt gebruik van een hybride numerieke en statistische benadering:

1. Wiskundige Modellering: Het fysische probleem wordt beschreven door een systeem van gekoppelde nietlineaire partiële differentiaalvergelijkingen (PDE's) die continuïteit, momentum, magnetische inductie, energie en concentratie bepalen. Deze vergelijkingen incorporeren effecten zoals het Dufour-effect, het Soret-effect, viskeuze dissipatie en de Arrhenius-activeringsenergie.
2. Gelijkenistransformatie: De onderliggende PDE's worden gereduceerd tot een systeem van nietlineaire gewone differentiaalvergelijkingen (ODE's) met behulp van geschikte gelijkenisvariabelen.
3. Numerieke Oplossing: Het resulterende randwaardeprobleem wordt numeriek opgelost met de MATLAB bvp4c solver, die gebruikmaakt van een collocatiemethode met een relatieve tolerantie van $10^{-6}$.
4. Statistische Analyse en Optimalisatie:
   • Response Surface Methodologie (RSM): Een Central Composite Design (CCD) wordt gebruikt om simulatieruns te genereren. Regressiemodellen worden ontwikkeld om de output-responsen, de wrijvingscoëfficiënt ($C_f$) en het Nusselt-getal ($Nu$), te voorspellen op basis van de inputparameters ($L_1, L_2, N_t$). ANOVA wordt uitgevoerd om de significantie van lineaire, kwadratische en interactietermen te valideren.
   • Artificiële Neurale Netwerken (ANN): Een ANN-model wordt ontwikkeld om de output-responsen te voorspellen, en de prestaties ervan worden vergeleken met RSM om de voorspellende nauwkeurigheid te beoordelen.
   • Gevoeligheidsanalyse: Op afgeleiden gebaseerde gevoeligheidsindices worden berekend om de invloed van elke parameter op de output-responsen te bepalen.

Belangrijkste Bijdragen en Resultaten
Het artikel presenteert de volgende specifieke bevindingen met betrekking tot de stromingsdynamica en warmteoverdrachtskarakteristieken:

• Snelheidsslip ($L_1$): Het verhogen van de snelheidsslipparameter van 0,7 naar 1,9 verandert het snelheidsprofiel, wat leidt tot een afname van de snelheid nabij de wand door verminderde interactie tussen wand en vloeistof, gevolgd door een toename in de buitenste grenslaag. Dit leidt tot een reductie in de wrijvingscoëfficiënt van ongeveer 19,3%. Echter, het Nusselt-getal neemt met ongeveer 3,3% af, wat duidt op een verminderde convectieve warmteoverdrachtsefficiëntie.
• Energyslip ($L_2$): Het verhogen van de energyslipparameter van 0,3 naar 0,7 verzwakt de thermische koppeling tussen de wand en de vloeistof. Dit resulteert in een reductie van het temperatuurprofiel en een afname van het Nusselt-getal met 19,3%. De concentratieprofielen nemen eveneens af naarmate de drijvende krachten voor thermische diffusie worden verminderd.
• Thermoforese ($N_t$): Een toename van de thermoforeseparameter (van 0,01 naar 0,05) bevordert de migratie van deeltjes weg van het verhitte oppervlak. Dit leidt tot verhoogde temperatuur-, dipool- en concentratieprofielen. De studie merkt een lichte reductie in de wrijvingscoëfficiënt op (0,2%) en een afname van het Nusselt-getal (2,8%) op onder deze specifieke variaties.
• Modelvalidatie: De numerieke resultaten vertonen een uitstekende overeenstemming met de benchmarkgegevens van Mustafa et al. [22].
• Vergelijkende Prestaties: De studie concludeert dat Artificiële Neurale Netwerken (ANN) een superieure voorspellende nauwkeurigheid vertonen vergeleken met de Response Surface Methodologie (RSM) voor deze specifieke stromingsconfiguratie.
• Optimalisatie: Multi-objective optimalisatie met behulp van wenselijkheidstermen identificeerde optimale parametersets. Bijvoorbeeld, een combinatie van $L_1 \approx 0,60$, $L_2 \approx 0,25$ en $N_t \approx 0,02$ werd gevonden om de wrijvingscoëfficiënt te minimaliseren terwijl de warmteoverdrachtssnelheden worden gehandhaafd of gemaximaliseerd, waarbij een wenselijkheidsscore van 1,000 werd bereikt.

Significantie en Claims
De auteurs claimen dat deze studie de eerste is die systematisch de gecombineerde effecten van snelheidsslip en energyslip parameters op de grenslaagstroom over een verticaal vlak onderzoekt om de thermische prestaties te verbeteren. De primaire significantie ligt in de parametrische gevoeligheidsanalyse, die onthult dat snelheidsslip een krachtiger effect heeft op het verbeteren van de warmteoverdrachtprestaties in geavanceerde thermische beheersystemen vergeleken met energyslip en thermoforese alleen.

De studie benadrukt dat hoewel slipcondities over het algemeen de wrijvingscoëfficiënt verminderen (gunstig voor drag-reductie), ze ook de neiging hebben het Nusselt-getal te verlagen (nadelig voor warmteoverdracht). Het gebruik van RSM en ANN maakt het mogelijk om specifieke operationele vensters te identificeren waar deze concurrerende doelstellingen in balans kunnen worden gebracht. De bevindingen worden gepresenteerd als een basis voor het optimaliseren van technische systemen met MHD-nanofluïde stromingen, zoals in thermische beheersapparaten, hoewel de auteurs expliciet vermelden dat hun huidige model beperkt is tot stationaire, tweedimensionale laminaire stroming met constante eigenschappen, waardoor onstabiele en turbulente regimes aan toekomstig onderzoek worden overgelaten.