Distributed Roughness-Induced Transition on a Blunt Body at Mach 6: a Numerical Investigation

Este estudo apresenta a primeira simulação numérica direta de rugosidade distribuída em um cilindro rombo Mach 6, revelando que o arranjo da rugosidade dita o mecanismo de transição ao promover ou modos de estrias sinuosas ou ondas T-S 2D, sendo estas últimas mantidas unicamente por um loop de feedback acústico interno que elimina a necessidade de forçamento externo.

O essencial

Imagine uma espaçonave cortando a atmosfera a seis vezes a velocidade do som. Para sobreviver ao calor intenso, sua superfície é coberta por um material especial que se queima lentamente (ablação) para proteger a nave. No entanto, à medida que esse material queima, ele não deixa uma superfície perfeitamente lisa; ele deixa para trás uma textura áspera e irregular, parecida com uma lixa.

Este artigo é uma simulação computacional de alta velocidade que faz uma pergunta simples, mas crítica: Como essas pequenas saliências na superfície transformam um fluxo de ar suave e ordenado em um fluxo de ar caótico e turbulento?

Aqui está o detalhamento de suas descobertas usando analogias do cotidiano:

1. A Configuração: O Cilindro de "Lixa"

Os pesquisadores construíram um modelo digital de um cilindro rombo (como o nariz de um foguete) voando a Mach 6. Em vez de uma superfície lisa, eles cobriram o objeto com pequenas "protuberâncias" artificiais (rugosidade) para imitar a textura arenosa deixada pelo material em combustão.

Eles testaram três maneiras diferentes de organizar essas protuberâncias:

• Alinhadas: Como soldados em fileiras e colunas perfeitas.
• Escalonadas: Como uma parede de tijolos, onde as protuberâncias de uma fileira são deslocadas em relação à fileira de trás.
• Aleatórias: Como pedras espalhadas em uma calçada, sem nenhum padrão.

2. A Teoria Antiga vs. A Nova Descoberta

Por muito tempo, os cientistas pensaram que a transição para a turbulência era causada por um "acúmulo lento" de energia, semelhante a um balanço que ganha altura se você o empurrar no ritmo certo ao longo do tempo. Isso é chamado de "crescimento transiente".

A Descoberta do Artigo:
A simulação mostrou que essa teoria do "acúmulo lento" não explica realmente o que está acontecendo aqui. As protuberâncias na superfície não apenas amplificaram a energia lentamente; elas atuaram como desestabilizadores. Elas pegaram o fluxo de ar e o tornaram instável imediatamente, transformando-o em um tipo específico de onda que cresce muito rápido.

Pense desta forma: a teoria antiga achava que as protuberâncias estavam apenas dando um empurrãozinho suave em um dominó para que ele caísse. A nova descoberta mostra que as protuberâncias estão, na verdade, chutando o dominó, fazendo com que ele colida com o próximo imediatamente.

3. Os Dois Tipos de "Ondas"

Dependendo de como as protuberâncias foram organizadas, o fluxo de ar reagiu de duas maneiras diferentes:

• A "Serpente" (Modo Sinuoso): Quando as protuberâncias estavam alinhadas (fileiras perfeitas), o fluxo de ar começou a oscilar de um lado para o outro como uma serpente. Este é um balanço muito específico e organizado.
• A "Onda Plana" (Ondas de Tollmien-Schlichting ou T-S): Quando as protuberâncias estavam escalonadas ou aleatórias, o fluxo de ar começou a ondular para cima e para baixo em um padrão de onda plana 2D. Este é um tipo clássico de onda geralmente encontrado em fluxos de ar muito mais lentos e de baixa velocidade, o que foi surpreendente de encontrar neste ambiente de alta velocidade.

O Insight Principal: A disposição das protuberâncias ditou qual "dança" o ar faria. A dança da "serpente" aconteceu com as protuberâncias alinhadas, enquanto a dança da "onda plana" aconteceu com as outras.

4. O Final do "Grampo de Cabelo" (Hairpin)

Uma vez que essas ondas cresceram o suficiente, elas desencadearam o estágio final do desastre. As "estrias" constantes de ar criadas pelas protuberâncias (que são como longas fitas invisíveis de ar lento) subitamente torceram e se quebraram em vórtices de grampo de cabelo (hairpin vortices).

Imagine um elástico que foi esticado. De repente, ele torce e forma um laço que se parece com um grampo de cabelo. Esses laços são o nascimento da turbulência. Uma vez que esses grampos se formam, o ar suave se quebra completamente em caos, e o calor na superfície da espaçonave aumenta drasticamente.

5. A Surpresa da "Câmara de Eco"

Uma das descobertas mais fascinantes foi como a turbulência começou de fato nos casos escalonados e aleatórios.

Normalmente, os cientistas pensam que você precisa de um "empurrão" externo (como uma rajada de vento ou uma vibração) para iniciar essas ondas. Mas a simulação mostrou algo autossustentável:

1. A turbulência começa em um ponto no cilindro.
2. Como o ar atrás da onda de choque está se movendo mais devagar que o som (subsônico), o ruído dessa turbulência viaja para trás, rio acima, como um eco.
3. Esse "eco" atinge a parte lisa da superfície à frente da turbulência e excita o ar ali, causando o início de uma nova turbulência.
4. Isso cria um ciclo de feedback: a turbulência gera ruído, o ruído viaja para trás e o ruído cria mais turbulência.

É como um microfone captando sua própria saída de som e criando um loop de feedback estridente, mas, neste caso, o "guicho" é o ar tornando-se turbulento.

Resumo

Este artigo usou um supercomputador para observar o fluxo de ar sobre um cilindro irregular de alta velocidade. Eles descobriram que:

• O padrão das protuberâncias decide exatamente como o ar se torna turbulento.
• A ideia antiga de "acúmulo lento de energia" não é a principal culpada; em vez disso, as protuberâncias desestabilizam diretamente ondas específicas.
• Essas ondas crescem até se torcerem em formas de "grampo de cabelo", fazendo com que o ar se torne caótico.
• Em alguns casos, a turbulência cria seu próprio "eco" que viaja para trás para iniciar o processo novamente, sem precisar de qualquer ajuda externa.

Isso ajuda os engenheiros a entender que as pequenas protuberâncias aleatórias deixadas pelos escudos térmicos em combustão não são apenas imperfeições menores; elas são as arquitetas primárias de como e quando a superfície de uma espaçonave atinge temperaturas perigosamente altas.

Resumo técnico

Resumo Técnico: Transição Induzida por Rugosidade Distribuída em um Corpo Rombo a Mach 6

Definição do Problema
A previsão da transição laminar-turbulenta (LTT) sobre geometrias rombas de alta velocidade permanece um desafio crítico para o design de engenharia, particularmente para veículos de entrada atmosférica onde a ablação da superfície cria rugosidade distribuída. Embora existam correlações empíricas (ex: Hollis 2017, 2019, 2025), estas carecem de uma base física clara. Esforços teóricos anteriores, como a correlação de Reshotko & Tumin (2004), baseiam-se em mecanismos de crescimento transiente, mas enfrentam limitações: assumem a escala linear de perturbações com a altura da rugosidade (válida apenas para $k/\delta$ pequeno), requerem constantes de ajuste para dados experimentais e têm dificuldade em explicar a transição em casos onde as alturas de rugosidade excedem a espessura da camada limite ($k/\delta > 1$). Além disso, investigações numéricas prévias frequentemente simulavam manchas de rugosidade em vez de rugosidade totalmente distribuída, ou falharam em observar a amplificação de instabilidade modal, deixando os mecanismos específicos que impulsionam a transição em fluxos de corpo rombo hipersônico com rugosidade totalmente distribuída pouco claros.

Metodologia
Este estudo apresenta a primeira Simulação Numérica Direta (DNS) de um cilindro rombo em escoamento transversal de Mach 6 com elementos de rugosidade distribuídos sobre toda a superfície. As simulações utilizam o solver CHAMPS para integrar as equações de Navier-Stokes compressíveis (CNSE) para um gás ideal. A geometria é um cilindro de 0,1524 m de diâmetro, modelado com base nas condições experimentais de Hollis (2017) (Run 38, Test 6975), com um número de Mach de fluxo livre de 5,98.

A rugosidade é modelada como uma distribuição sinusoidal com uma amplitude fixa ($A \approx 35\%$ da espessura da camada limite) e comprimento de onda correspondente à rugosidade nominal de grão de areia. Três configurações distintas de faseamento transversal foram investigadas:

1. Deslocada (Staggered): As linhas estão $180^\circ$ fora de fase.
2. Alinhada (Aligned): As linhas estão em fase ($0^\circ$ de deslocamento).
3. Aleatória (Random): As linhas possuem deslocamentos de fase distribuídos aleatoriamente.

Para analisar o processo de transição, os autores empregaram:

• DNS: Para capturar a evolução não linear completa do escoamento.
• Equações de Perturbação Linear (LDE): Resolvidas sobre fluxos médios temporais do DNS para isolar mecanismos de instabilidade usando forçamento de pulso.
  im
• Teoria de Estabilidade Linear (LST): Aplicada a perfis com média transversal para comparação.
• Análise Espectral: Transformadas Rápidas de Fourier (FFT) dos campos de velocidade e pressão para identificar frequências e números de onda dominantes.

Principais Resultados

• Fluxo Médio e Listras (Streaks): Todas as configurações geraram listras de velocidade longitudinal via um mecanismo de lift-up impulsionado por pares de vórtices contra-rotativos (CVPs). A configuração alinhada produziu as listras estáticas mais fortes devido ao faseamento consistente que eleva o fluido de baixo fluxo de massa para partes mais altas da camada limite. As configurações deslocada e aleatória produziram padrões de listras mais fracos e complexos. Contrário às previsões da teoria de crescimento transiente ótimo, as amplitudes das listras exibiram baixo crescimento relativo (menos de um fator de 2) e permaneceram não lineares por todo o domínio, desafiando a utilidade do crescimento transiente linear como ferramenta preditiva para este regime.

• Mecanismos de Instabilidade: O processo de transição foi impulsionado por instabilidades convectivas, mas o tipo de modo específico dependeu do faseamento da rugosidade:

  • Caso Alinhado: Dominado por um modo de listra sinuoso fundamental. As ondas de Tollmien-Schlichting (T-S) foram fracamente desestabilizadas.
  • Casos Deslocado e Aleatório: Dominados por ondas T-S 2D (modo varicoso). Estas configurações exibiram uma desestabilização significativa das ondas T-S, levando a um crescimento exponencial rápido.
  • Ruptura (Breakdown): Em todos os casos, uma vez que as ondas de instabilidade dominante cresceram até amplitude suficiente, elas desencadearam a formação de vórtices de cabelo (hairpin vortices) a partir das listras estáticas, levando à ruptura e à turbulência.

• Locais de Transição: A DNS previu locais de transição ($\Theta_{tr}$) significativamente mais a jusante (aprox. $11^\circ$–$17^\circ$) do que as correlações empíricas de última geração (que previam $6^\circ$–$13^\circ$). O caso alinhado transitou mais cedo entre os casos de rugosidade devido ao forçamento necessário para iniciar a LTT, enquanto os casos deslocado e aleatório transitaram mais tarde.

• Mecanismo de Feedback: Um achado inovador nos casos deslocado e aleatório foi a identificação de um ciclo de feedback. A turbulência gerada na região subsonica do cilindro dispersou ondas acústicas que propagaram para montante. Essas ondas acústicas lentas excitaram ondas T-S no lado oposto do cilindro (através da linha de estagnação) em seu ponto neutro. Este mecanismo, que não requer forçamento externo, sustenta o processo de transição.

Significância e Alegações
Os autores alegam que este trabalho fornece uma nova compreensão física da transição induzida por rugosidade em corpos rombos hipersônicos, desafiando a dependência predominante da teoria de crescimento transiente. As principais contribuições são:

1. Rejeição do Crescimento Transiente como Principal Impulsionador: O estudo demonstra que para rugosidade totalmente distribuída com $k/\delta \approx 0,35$, o desenvolvimento não linear das listras e as baixas razões de amplificação observadas tornam a teoria de crescimento transiente ótimo um preditor insuficiente do processo de transição real.
2. Desestabilização de Instabilidades Modais: O artigo estabelece que a rugosidade distribuída pode efetivamente desestabilizar instabilidades modais (especificamente ondas T-S) mesmo em corpos rombos onde o escoamento de parede lisa é tipicamente estável em termos de eigenmodo. A configuração da rugosidade dita se a transição é impulsionada por modos de listra ou ondas T-S.
3. Descoberta de um Ciclo de Feedback Acústico: A identificação de um mecanismo de feedback auto-sustentável, onde a acústica gerada pela turbulência excita ondas T-S no lado oposto do corpo, oferece uma explicação potencial para a transição em regiões subsonicas de escoamento hipersônico, um fenômeno anteriormente não explicado por mecanismos lineares.
4. Sensibilidade ao Faseamento: Os resultados destacam a extrema sensibilidade da LTT ao faseamento relativo dos elementos de rugosidade, sugerindo que topologias simples podem delimitar o comportamento de padrões de rugosidade mais complexos e aleatórios.

Os autores concluem que, embora o mecanismo de feedback acústico proposto exija uma prova mais rigorosa, a exposição desses mecanismos de instabilidade via DNS de alta fidelidade sugere uma mudança fundamental na compreensão das perturbações responsáveis pela LTT em corpos rombos hipersônicos.