A carrier-wave factored one-way Navier--Stokes method for boundary-layer instability modelling

Resumo Técnico: M-OWNS para Modelagem de Instabilidade de Camada Limite

Enunciado do Problema
Prever a transição laminar–turbulenta requer o rastreamento da evolução das perturbações através da amplificação em escoamentos de base em desenvolvimento. Enquanto a análise de autovalores locais (LST) é computacionalmente barata, mas limitada por suposições de fluxo paralelo, e os resolvedores de Navier–Stokes harmônicos linearizados globais (LHNS) são gerais, mas proibitivamente caros devido a restrições de memória e inversão de matrizes, os métodos de marcha espacial oferecem um meio-termo. No entanto, a marcha espacial padrão enfrenta uma dificuldade fundamental: o operador de Navier–Stokes linearizado suporta modos acústicos que se propagam para montante, o que pode desestabilizar uma marcha descendente ingênua.

As soluções existentes possuem compensações (trade-offs):

• Equações de Estabilidade Parabolizadas (PSE): Utilizam o fatoramento de onda portadora para deslocar o espectro de jusante em direção à origem, permitindo uma resolução mais grosseira. Contudo, dependem de estratégias de regularização (por exemplo, negligenciar gradientes de pressão no sentido do escoamento ou amortecimento implícito) para suprimir modos de montante, o que introduz erros característicos. Além disso, a PSE sofre de uma restrição de passo mínimo; refinar o passo para capturar perturbações mono-modais ou não-modais frequentemente recaptura os modos instáveis de montante, desestabilizando a marcha.
• Navier–Stokes Unidirecional (OWNS): Remove os modos de montante inteiramente via uma projeção recursiva (OWNS-R), retendo o gradiente de pressão no sentido do escoamento e evitando os erros de regularização da PSE. No entanto, a OWNS não emprega o fatoramento de onda portadora. Consequentemente, deve resolver os autovalores em seus números de onda totais a partir da origem, necessitando de uma resolução streamwise fina. Essa resolução fina amplia a janela espectral do integrador, exigindo uma projeção mais complexa para suprimir uma gama mais ampla de modos de montante, o que é limitado pelo acúmulo de erro numérico nos parâmetros de recursão.

Metodologia: M-OWNS
O artigo apresenta o M-OWNS, um método de marcha espacial que integra o fatoramento de onda portadora da PSE na estrutura de Navier–Stokes unidirecional recursiva (OWNS-R).

1. Fatoramento de Onda Portadora: O estado da perturbação $\phi$ é decomposto em um envelope de variação lenta $\tilde{\phi}$ e uma onda portadora de variação rápida $\exp(i\Theta_0)$, onde $d\Theta_0/dx = \alpha_0$. Isso desloca o espectro do operador espacial por um escalar $\alpha_0$.
2. Projeção Unidirecional: O operador de projeção $P$ da OWNS-R é aplicado ao sistema fatorado. Os autores provam que o mesmo operador de projeção que paraboliza o sistema não fatorado também paraboliza o sistema fatorado por onda portadora, removendo modos de propagação para montante sem a necessidade de modificar os parâmetros de projeção.
3. Seleção da Portadora ($\alpha_0$):
   • Variante Iterativa: $\alpha_0$ é atualizado em cada estação usando o quociente de Rayleigh, convergindo para o autovalor dominante (semelhante à PSE).
   • Variante de Portadora Fixa: $\alpha_0$ é definido como um valor representativo derivado do escoamento de base (por exemplo, o ponto final do ramo de vorticidade/entropia $\alpha_1 = (\omega - \beta W_\infty)/U_\infty$). Esta variante não requer iterações e mantém o mesmo custo por passo da OWNS não fatorada.

Principais Contribuições

• Vantagem de Resolução Espectral: Uma análise de primeira ordem demonstra que, para qualquer autofunção excitada cujo autovalor esteja mais próximo do número de onda da portadora $\alpha_0$ do que da origem, o sistema fatorado resolve o modo em um passo streamwise mais grosseiro do que o sistema não fatorado. O fatoramento efetivamente translada o "disco de resolução" do integrador da origem para o agrupamento (cluster) dos modos excitados.
• Desacoplamento de Resolução e Regularização: Ao combinar o fatoramento com a projeção unidirecional, a M-OWNS alcança os benefícios de resolução da PSE (passos mais grosseiros) sem as restrições de tamanho de passo ou erros de regularização (ex: gradientes de pressão negligenciados) que assolam a PSE.
• Robustez em Regimes Multimodais: O método mostra-se eficaz tanto para perturbações mono-modais (onde a iteração é ótima) quanto para perturbações multimodais ou de banda larga (onde uma portadora fixa é preferível para evitar problemas de convergência do quociente de Rayleigh associados a autofunções não ortogonais).

Resultados e Validação
O método foi validado em camadas limite de placa plana incomprimíveis e subsônicas, perturbações de crossflow tridimensionais sobre um cilindro varrido e uma camada limite hipersônica de Mach 4.5. Quatro configurações de excitação foram testadas: excitação de autofunção na entrada, sucção/sopro na parede, excitação de freestream multimodal e excitação de entrada randomizada.

• Eficiência de Resolução: Em todos os casos de teste, a M-OWNS capturou amplitudes de perturbação, campos de radiação acústica e sequências de sincronização modal em resoluções streamwise mais grosseiras que a OWNS não fatorada.
  • Para o caso do cilindro varrido, a OWNS não fatorada exigiu $N_x \approx 8000$ passos para convergência, enquanto a M-OWNS de portadora fixa convergiu em $N_x \approx 1000$ (uma redução de 8 vezes), e a M-OWNS iterativa convergiu em $N_x \approx 150$ (uma redução de 47 vezes).
  • Para o caso de entrada randomizada hipersônica (o teste multimodal mais rigoroso), a OWNS não fatorada exigiu $N_x \approx 6000$, enquanto a M-OWNS de portadora fixa convergiu em $N_x \approx 3000$ (uma redução de 2 vezes).
• Acurácia: Os resultados da M-OWNS foram validados contra PSE, OWNS não fatorada e resolvedores LHNS totalmente elípticos. No caso de sucção/sopro na parede, os resultados da M-OWNS coincidiram com a referência LHNS, enquanto a OWNS não fatorada, na mesma resolução, subestimou as amplitudes.
• Custo: O custo computacional total (total de resoluções lineares) foi reduzido por fatores de dois a oito nos casos de teste. A variante de portadora fixa alcançou essas reduções com custos por passo idênticos aos da OWNS não fatorada, pois elimina o overhead iterativo.

Significância
O artigo afirma que a M-OWNS oferece uma vantagem distinta de resolução numérica e eficiência na modelagem da evolução do estado de perturbação e instabilidade. Ao deslocar a região espectral resolvida para o agrupamento fisicamente relevante de autovalores, o método permite passos streamwise significativamente mais grosseiros sem sacrificar a precisão ou introduzir os erros de regularização inerentes à PSE. A variante de portadora fixa é destacada como particularmente significativa, pois alcança esses ganhos sem aumentar o custo computacional por passo ou exigir convergência iterativa, tornando-a uma alternativa robusta para problemas de instabilidade mono-modais e multimodais complexos. Os autores observam que o método é aplicável a escoamentos gerais onde o espectro excitado se agrupa próximo a um número de onda específico, fornecendo uma solução prática para as restrições de resolução da OWNS padrão.