Canoa Quebrada

O que é e como funciona na prática

A expressão canoa quebrada designa, no contexto em que costumo ouvir, a situação em que uma simulação de engenharia ou um código numérico simplesmente se separa — os valores explodem, o resíduo fica NaN ouinf, e o trabalho de uma tarde vai embora. Isso não é abstrato. Acontece porque algo no modelo foge do regime de estabilidade que você assumiu no início.

Por que a maioria dos casos de canoa quebrada acontece

Na prática, os motivos se repetem. Malha muito distorcida perto de geometrias finas gera elementos com razão de aspecto ruim e o condition number da matriz dispara. Passo de tempo grande demais em esquema explícito viola CFL. Condição de contorno incompatível com o campo inicial gera choque na primeira iteração. Material com rigidez ou viscosidade absurda quebra o equilíbrio de unidades. E, claro, tolerância de convergência folgada demais faz o solver parar antes de chegar em algum lugar utilizável, o que deixa resíduos disfarçados que só aparecem horas depois.

Um problema real que eu tive e como resolvi

Havia um caso em que um modelo de escoamento compressível com malha adaptativa começava bem, mas quebrava na iteração 47, sempre no mesmo local. O sintoma clássico: pressão e temperatura subiam e o resíduo oscilava sem tendência. Eu revisei a malha inteira, aumentei o tempo de setup, mudei esquema de transporte, e nada. A raiz era um detalhe chatíssimo: o contorno de saída estava demasiado próximo de uma região de gradiente forte gerada por uma obstrução interna, e a condição de pressão imposta refletia onda de pressão de volta para o domínio. A solução foi mover o outlet para pelo menos dez comprimentos característicos abaixo da geometria, adicionar uma zona de amortecimento com difusividade aumentada gradualmente, e refinar só aquela camada limite com elementes alongados alinhados ao fluxo. O tempo de setup dobrou, mas a simulação passou de 47 para mais de 200 iterações estáveis, e os resultados deixaram de depender do refinamento local.

Como evitar que a canoa quebre, passo a passo

Comece devagar. Se o modelo é não linear ou transiente, rode uma versão simplificada primeiro com materiais lineares, passo de tempo pequeno e condições suaves. Valide a malha em regiões de alto gradiente antes de confiara no resultado final. Cheque números adimensionais locais — número de Mach, Reynolds, Courant — em cada cédula crítica, não apenas na média do domínio. Use inicição física, não chute aleatório; campos de pressão e velocidade coerentes reduzem drasticamente a chance de divergência nas primeiras iterações. Monitore resíduos e grandezas físicas ao mesmo tempo. Resíduo baixo com energia cinética variando selvagemente quase sempre indica ilusão de convergência. Configuração de contorno merece atenção específica. Evite impor perfis completos em regiões onde o fluxo ainda não se desenvolveu. Prefira inlet de velocidade ou massa com pressão livre no outlet, ou vice-versa, conforme o regime. Se usar paredes sem slip, verifique y+ compatível com o modelo de turbulência adotado. Em casos compressíveis, cuidado com reflexões em boundaries — zonas de relaxamento ou técnicas de sponge ajudam, mas precisam de extensão suficiente para não introduzir artefato numérico visível.

O que mais costuma dar errado

Unidades inconsistentes são o erro mais barato e mais frequente. Misturar mm com m, g com kg, ou segundos com horas gera coeficientes fora da escala esperada e o solver entra em colapso sem aviso. Outro ponto cego é a definição de propriedades dependentes de temperatura ou pressão. Funções mal contínuas, especialmente com descontinuidade numérica em intervalo de iteração, geram saltos na matriz Jacobiana e divergência repentina. Ainda há o problema de malhas móveis ou de reconexão automática que, se mal parametrizada, cria volumes negativos e mata a simulação no instante seguinte.

Ferramentas e onde encontrar referências

Não existe um pacote único chamado canoa quebrada que se possa baixar e resolver tudo. O termo descreve o comportamento, não um software. O que vale é reunir boas práticas de pré-processamento, configuração de Solver, e análise de diagnósticos. Bibliotecas abertas como o OpenFOAM têm exemplos que cobrem a maioria dos regimes comuns, e a comunidade oferece scripts de verificação de malha, como checkMesh, que identificam problemas antes de rodar. Para Python, pacotes como FEniCS, FiPy e deal.II permitem implementar esquemas próprios, mas exigem cuidado extra com estabilização numérica, especialmente em regimes adominados por convecção. Se o foco é CFD estruturado, soluções como SU2 ou codes acadêmicos com esquemas de captura de choque ajudam em fluxos compressíveis. A documentação técnica desses projetos costuma ser mais honesta sobre limitações do que manuais comerciais, o que economiza tempo na fase de diagnóstico.

Quando a abordagem padrão não funciona

Há cenários em que o modelo convencional simplesmente não aguenta. Geometrias com escalas muito separadas, como microcanais acoplados a volumes grandes, geram stiffness numérica que métodos padrão resolvem com dificuldade. Regimes altamente transientes com tempos de resposta distintos exigem esquemas implícitos semi-adiabáticos ou métodos de multiple scaling. Fluidos não newtonianos com shear rate extremo podem precisar de regularização de viscosidade artificial para não travar o solucionador. Nesses casos, a opção mais segura é particionar o domínio, resolver subdomínios com métodos adequados e acoplar por interfaces conservativas, ou recorrer a reformulações matemáticas que suavizam as descontinuidades sem alterar a física relevante.

Checklist rápido antes de ligar o solver

Revise unidades e conversões. Verifique a qualidade da malha e os parâmetros de aspect ratio, skewness e orthogonality. Confirme condições de contorno compatíveis com o regime esperado. Teste passo de tempo com CFL abaixo do limite recomendado pelo esquema. Inicie com campo fisicamente plausível. Monitore Grandezas globais, não apenas resíduos. anote o histórico completo para posterior comparação. E, acima de tudo, não confie em convergência única. Rodadas múltiplas com parâmetros levemente diferentes revelam instabilidades que uma execução isolada esconde. A experiência mostra que a maior parte dos casos de canoa quebrada se resolve com paciência na configuração inicial e diagnóstico sistemático, não com ajuste cego de tolerâncias. O solver muitas vezes está certo ao falhar; o modelo é que está pedindo para ser reconsiderado.