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

Resumo Técnico: Avaliação de Sensibilidade Paramétrica do Escoamento com Deslizamento e Efeitos Termoforéticos sobre uma Superfície de Estiramento Vertical Usando a Metodologia de Superfície de Resposta

Enunciado do Problema
Este estudo investiga o escoamento de camada limite magnetohidrodinâmico (MHD) de um nanofluido eletricamente condutor sobre uma chapa de estiramento vertical. A pesquisa foca nos efeitos acoplados do deslizamento de velocidade (velocity slip) e do deslizamento de energia (deslizamento térmico/thermal slip) e da termoforese nas características de transferência de momento, calor e massa. Embora a literatura anterior tenha explorado diversas condições de deslizamento e comportamentos de nanofluidos, este trabalho aborda especificamente a interação sistemática dos parâmetros de deslizamento de velocidade ($L_1$) e de deslizamento de energia ($L_2$), juntamente com o parâmetro de termoforese ($N_t$), para otimizar o desempenho térmico. Os fenômenos físicos governantes incluem campos magnéticos induzidos, radiação térmica, movimento browniano e reações químicas dentro de um regime de escoamento laminar bidimensional e estacionário.

Metodologia
O estudo emprega uma abordagem numérica e estatística híbrida:

1. Modelagem Matemática: O problema físico é descrito por um sistema de equações diferenciais parciais (PDEs) não lineares acopladas que regem a continuidade, o momento, a indução magnética, a energia e a concentração de espécies. Estas equações incorporam efeitos como o efeito Dufour, o efeito Soret, dissipação viscosa e energia de ativação de Arrhenius.
2. Transformação de Similaridade: As PDEs governantes são reduzidas a um sistema de equações diferenciais ordinárias (ODEs) não lineares utilizando variáveis de similaridade apropriadas.
3. Solução Numérica: O problema de valor de contorno resultante é resolvido numericamente usando o solver bvp4c do MATLAB, que utiliza um método de colocação com uma tolerância relativa de $10^{-6}$.
4. Análise Estatística e Otimização:
   • Metodologia de Superfície de Resposta (RSM): Um Design Composto Central (CCD) é utilizado para gerar execuções de simulação. Modelos de regressão são desenvolvidos para prever o coeficiente de atrito superficial ($C_f$) e o número de Nusselt ($Nu$) com base nos parâmetros de entrada ($L_1, L_2, N_t$). A ANOVA é realizada para validar a significância dos termos lineares, quadráticos e de interação.
   • Redes Neurais Artificiais (ANN): Um modelo de ANN é desenvolvido para prever as respostas de saída, e seu desempenho é comparado ao RSM para avaliar a precisão preditiva.
   • Análise de Sensibilidade: Índices de sensibilidade baseados em derivadas são calculados para determinar a influência de cada parâmetro nas respostas de saída.

Principais Contribuições e Resultados
O artigo apresenta os seguintes achados específicos em relação à dinâmica do escoamento e às características de transferência de calor:

• Deslizamento de Velocidade ($L_1$): O aumento do parâmetro de deslizamento de velocidade de 0,7 para 1,9 altera o perfil de velocidade, causando uma diminuição na velocidade próxima à parede devido à redução da interação fluido-parede, seguida por um aumento na camada limite externa. Isso leva a uma redução no atrito superficial de aproximadamente 19,3%. No entanto, o número de Nusselt diminui aproximadamente 3,3%, indicando uma redução na eficiência da transferência de calor por convecção.
• Deslizamento de Energia ($L_2$): O aumento do parâmetro de deslizamento de energia de 0,3 para 0,7 enfraquece o acoplamento térmico entre a parede e o fluido. Isso resulta em uma redução do perfil de temperatura e em uma diminuição do número de Nusselt em 19,3%. Os perfis de concentração também diminuem conforme as forças motrizes de difusão térmica são mitigadas.
• Termoforese ($N_t$): Um aumento no parâmetro de termoforese (de 0,01 para de 0,05) potencializa a migração de partículas para longe da superfície aquecida. Isso leva a perfis elevados de temperatura, dipolo e concentração. O estudo observa uma ligeira redução no atrito superficial (0,2%) e uma diminuição no número de Nusselt (2,8%) sob estas variações específicas.
• Validação do Modelo: Os resultados numéricos mostram excelente concordância com os dados de referência de Mustafa et al. [22].
• Desempenho Comparativo: O estudo conclui que as Redes Neurais Artificiais (ANN) demonstram uma precisão preditiva superior em comparação com a Metodologia de Superfície de Resposta (RSM) para esta configuração de escoamento específica.
• Otimização: A otimização multiobjetivo usando funções de desejabilidade identificou conjuntos de parâmetros ideais. Por exemplo, uma combinação de $L_1 \approx 0,60$, $L_2 \approx 0,25$ e $N_t \approx 0,02$ foi encontrada para minimizar o atrito superficial enquanto mantém ou maximiza as taxas de transferência de calor, alcançando um escore de desejabilidade de 1,000.

Significância e Alegações
Os autores afirmam que este estudo é o primeiro a investigar sistematicamente os efeitos combinados dos parâmetros de deslizamento de velocidade e de deslizamento de energia no escoamento de camada limite sobre uma chapa vertical para melhorar o desempenho térmico. A principal significância reside na avaliação de sensibilidade paramétrica, que revela que o deslizamento de velocidade tem um efeito mais potente no aumento do desempenho de transferência de calor em sistemas avançados de gestão térmica em comparação com o deslizamento de energia e a termoforese isoladamente.

O estudo enfatiza que, embora as condições de deslizamento geralmente reduzam o atrito superficial (benéfico para a redução de arrasto), elas também tendem a reduzir o número de Nusselt (prejudicial para a transferência de calor). O uso de RSM e ANN permite a identificação de janelas operacionais específicas onde esses objetivos conflitantes podem ser equilibrados. Os achados são apresentados como uma base para a otimização de sistemas de engenharia envolvendo escoamentos de nanofluidos MHD, tais como dispositivos de gestão térmica, embora os autores afirmem explicitamente que o modelo atual é restrito ao escoamento laminar bidimensional e estacionário, com propriedades constantes, deixando regimes não estacionários e turbulentos para investigações futuras.