Topology optimization of passively moving rigid bodies in unsteady flows

Este artigo propõe um método de otimização de topologia para o projeto de corpos rígidos de movimento passivo em fluxos fluidos não estacionários, acoplando a dinâmica de corpo rígido com equações de fluidos cinéticos de rede e aplicando uma análise de sensibilidade baseada em adjuntos para otimizar formas para movimentos translacionais e rotacionais.

O essencial

Imagine que você esteja tentando projetar a vela perfeita para um barco ou a pá perfeita para uma turbina eólica. Normalmente, os engenheiros começam com uma forma, colocam-na no vento, veem como ela se move, ajustam a forma e tentam novamente. Este artigo apresenta uma maneira automatizada e superinteligente de fazer isso, mas com um toque: em vez de a forma apenas ficar lá parada enquanto o vento sopra sobre ela, a forma tem permissão para se mover porque o vento a está empurrando.

Aqui está uma análise de como essa "mágica" funciona, usando analogias simples:

1. A Grande Ideia: O "Dançarino Passivo"

A maioria dos programas de computador que projetam formas assume que o objeto está colado ao chão (como uma ponte ou um tubo estacionário). Se você quiser projetar uma parte móvel, como uma pá de ventilador, você geralmente diz ao computador: "Gire esta pá a 100 rotações por minuto", e o computador calcula o fluxo de ar.

Este artigo inverte o roteiro. Ele trata o objeto como um dançarino em uma pista de dança.

• O Jeito Antigo (Ativo): Você diz ao dançarino exatamente como se mover e observa como o ar se move ao redor dele.
• O Novo Jeito (Passivo): Você não diz ao dançarino como se mover. Você apenas define a música (o vento) tocando e pede ao computador para projetar o corpo do dançarino de modo que a música naturalmente o empurre para girar ou deslizar o máximo possível. O movimento do dançarino é um resultado do vento, não um comando.

2. O Truque das Duas Grades: O Mapa vs. O Terreno

Para fazer isso funcionar, o computador usa um truque inteligente chamado "grades separadas". Imagine que você está desenhando o mapa de uma ilha móvel em uma folha de papel quadriculado.

• A Grade de Design (O Mapa): É onde a forma é desenhada. É como um bloco de notas. O computador decide onde está o material "sólido" (a ilha) e o espaço "vazio" (a água) aqui.
• A Grade de Análise (O Terreno): É onde a física acontece. É uma grade fixa de água e vento.

A cada fração minúscula de segundo, o "Mapa" (a forma) fisicamente se move e rotaciona. Então, o computador projeta esse mapa móvel sobre a grade de "Terreno" fixa para calcular como o vento empurra contra ele. Após o vento empurrar, o computador calcula como o objeto deve se mover a seguir, atualiza o Mapa e repete o ciclo. É como tirar uma foto de um objeto em movimento, calcular a força do vento, mover o objeto e tirar a próxima foto instantaneamente.

3. A Força "Fantasma" (Força de Brinkman)

Como o computador sabe onde está o objeto sólido? Ele usa um conceito chamado força de Brinkman.
Pense na área de design como uma sala cheia de mel invisível e pegajoso.

• Onde há material sólido, o mel é super espesso e pegajoso. O vento não pode passar por ele; ele apenas empurra contra a superfície.
• Onde há espaço vazio, o mel é fino ou inexistente, então o vento flui livremente.
  O computador não precisa desenhar uma linha dura; ele apenas ajusta a "pegajosidade" do mel em cada ponto. Se a pegajosidade for alta, é uma parede; se for baixa, é ar. Isso permite que a forma se transforme suavemente de uma forma para outra.

4. A Matemática da "Viagem no Tempo" (Método Adjoint)

Para encontrar a forma perfeita, o computador tem que saber: "Se eu mudar este minúsculo ponto de material aqui, o quanto o objeto se moverá melhor?"
Calcular isso para cada um desses pontos levaria uma eternidade. Por isso, os autores usam um método chamado Método da Variável Adjoint.

• A Analogia: Imagine que você está tentando encontrar o melhor caminho para subir uma montanha no escuro. Em vez de percorrer todos os caminhos possíveis para frente para ver qual é o melhor, você aponta uma lanterna para trás a partir do cume. A luz mostra exatamente quais passos levam para cima de forma mais eficiente.
• Neste artigo, a "lanterna" corre para trás através do tempo, calculando como as forças do vento e o movimento do objeto teriam reagido a cada pequena mudança na forma. Isso fornece ao computador um "mapa de sensibilidade", dizendo exatamente onde adicionar ou remover material para obter o melhor resultado.

5. Os Resultados: O Que Eles Construíram?

A equipe testou isso em três cenários:

1. A Vela 2D: Eles projetaram uma forma que começa parada e é empurrada pelo vento para deslizar horizontalmente. O resultado parece uma asa de avião curva (aerofólio). O vento empurra com mais força na parte superior do que na inferior, criando sustentação que puxa o objeto para frente.
2. A Turbina 2D: Eles projetaram uma forma que gira. O resultado parece uma hélice de quatro pás. O vento atinge as pás curvas, criando uma torção que a faz girar.
3. A Turbina 3D: Eles fizeram o mesmo em 3D. O resultado parece uma turbina eólica do mundo real.

6. A Questão do "Escala de Cinza"

Nesses designs de computador, as bordas das formas nem sempre são linhas perfeitas de preto e branco. Às vezes, elas são "escala de cinza" — um pouco de ambos, sólido e ar.

• Nos exemplos 2D, os autores descobriram que mesmo que tornassem a forma perfeitamente nítida (preto e branco), o desempenho era quase o mesmo. As bordas "esfumaçadas" não prejudicaram o resultado.
• No exemplo 3D, as bordas "esfumaçadas" importaram mais. Como a grade do computador era um pouco "robusta" (baixa resolução), as bordas esfumaçadas alteraram a forma como o vento atingia as pás. Isso sugere que, para formas 3D complexas, precisamos de um "mapa" mais fino para obter um resultado perfeito.

Resumo

Este artigo apresenta uma nova maneira de projetar máquinas móveis (como velas ou turbinas) onde o computador descobre a forma e o movimento simultaneamente. Ele trata o objeto como um dançarino passivo empurrado pelo vento, usa o truque do "mel pegajoso" para definir a forma e executa uma simulação matemática de tempo reversível para encontrar a forma mais eficiente. O resultado são formas que naturalmente parecem asas e hélices, otimizadas para se mover o máximo possível ou girar o mais rápido possível sob forças fluídas.

Resumo técnico

Resumo Técnico: Otimização Topológica de Corpos Rígidos de Movimento Passivo em Fluxos Não Estacionários

Definição do Problema
Este estudo aborda a otimização topológica de corpos rígidos que são impulsionados por forças induzidas pelo fluxo de fluido, como velas e turbinas. Diferente da otimização topológica tradicional em dinâmica de fluidos, onde o movimento do corpo rígido é prescrito (ativo) para impulsionar o fluxo, ou onde partículas são rastreadas em um acoplamento de via única onde elas não afetam o fluido, este trabalho foca em um sistema de acoplamento de via dupla. Aqui, o movimento do corpo rígido é passivo, o que significa que sua trajetória e rotação são inteiramente determinadas pelas forças fluídas que atuam sobre ele. O desafio governante envolve resolver as equações de fluxo de fluido não estacionário simultaneamente com as equações de movimento do corpo rígido, onde as forças do fluido servem como entradas para as equações de movimento, e a posição do corpo serve como uma condição de contorno para o fluido.

Metodologia
O framework proposto integra diversas técnicas numéricas para lidar com a complexa interação entre o sólido móvel e o fluido não estacionário:

• Equações Governantes: O fluxo de fluido é modelado usando as equações de Navier-Stokes incomprimíveis, resolvidas via Esquema Cinético de Rede (LKS), uma versão estendida do Método de Lattice Boltzmann (LBM) adequada para fluxos não estacionários. A dinâmica do corpo rígido é governada pelas equações de Newton-Euler, onde as forças e torques são computados como a reação da força de Brinkman dentro do domínio do fluido.
• Grades de Design-Análise Separadas: Um componente central da metodologia é a separação da grade de design (onde a forma é definida) da grade de análise (onde as variáveis de estado do fluido são computadas). A grade de design sofre movimento de corpo rígido em cada passo de tempo e é mapeada na grade de análise usando uma função de peso. Isso permite a representação direta do movimento do corpo rígido sem a necessidade de remalha do domínio do fluido.
• Análise de Sensibilidade Adjoint: Para computar as sensibilidades de design de forma eficiente, o estudo emprega o Esquema Cinético de Rede Adjoint (ALKS). Os autores derivam equações adjoint contínuas que acoplam as equações de fluxo de fluido com as equações de movimento do corpo rígido. Esta formulação trata as variáveis de movimento do corpo rígido (posição, velocidade, ângulo de rotação e velocidade angular) como campos de estado adicionais dentro do sistema adjoint.
• Algoritmo de Otimização: A forma é atualizada usando o Método das Assíntotas Móveis (MMA). Para garantir designs manufaturáveis e binários, utiliza-se um filtro de densidade com projeção de Heaviside. Um esquema de continuação é aplicado ao parâmetro de convexidade e ao parâmetro de rigidez da projeção para guiar a otimização de um estado inicial difuso para uma distribuição sólido-fluido clara.
• Aproximação Periódica: Para problemas envolvendo movimento cíclico (ex: turbinas), é empregada uma aproximação de problema periódico. Os valores iniciais para os campos de estado e adjoint no início de uma iteração de otimização são definidos como os valores terminais da iteração anterior, reduzindo o custo computacional ao evitar a simulação de fases de partida transientes para cada iteração.

Principais Contribuições

1. Formulação de Movimento Passivo: O artigo apresenta um novo framework de otimização topológica especificamente para corpos rígidos impulsionados por forças fluídas, distinguindo-se do movimento ativo ou do rastreamento de partículas de acoplamento de via única.
2. Derivação de Adjoint Acoplado: Fornece uma derivação rigorosa de equações adjoint e sensibilidades de design que acoplam totalmente a dinâmica de fluidos com as equações de movimento do corpo rígido, tratando as variáveis cinemáticas do corpo como campos de estado.
3. Implementação de Separação de Grades: O estudo demonstra a eficácia da abordagem de separação de grades para problemas de corpo rígido móvel e não estacionário, facilitando a transferência de quantidades como o coeficiente de Brinkman e a velocidade do sólido entre grades móveis e estacionárias.
4. Validação Numérica: O método é verificado contra aproximações de diferenças finitas para os casos translacional e rotacional, confirmando a precisão das sensibilidades derivadas.

Resultados
O método foi aplicado a três exemplos numéricos:

• Vela 2D: A otimização maximizou o deslocamento horizontal de um corpo rígido. A forma resultante assemelhou-se a um aerofólio, gerando forças de sustentação que impulsionaram o corpo para frente. A forma otimizada superou significativamente uma forma de referência com uma superfície superior plana.
• Turbina 2D: O objetivo foi maximizar o ângulo de rotação. Empregando uma condição periódica cíclica, a otimização produziu uma forma com quatro lâminas semelhantes a aerofólios que rotacionam efetivamente sob forças fluídas. Os resultados indicaram que a presença de regiões de tons de cinza na forma otimizada teve um efeito desprezível no desempenho quando binarizada.
• Turbina 3D: Uma forma rotativa tridimensional foi otimizada, resultando em uma estrutura do tipo turbina. No entanto, devido a limitações computacionais resultantes de uma resolução de malha mais grosseira, a versão binarizada apresentou diferenças não desprezíveis em fase e período em comparação com a forma otimizada, destacando a sensibilidade de estruturas finas 3D à resolução da malha.

Significância e Alegações
Os autores afirmam que o método proposto otimiza efetivamente as formas de dispositivos como velas e turbinas, onde o desempenho é ditado pela interação entre forças fluídas e o movimento de corpo rígido. O estudo demonstra que a análise não estacionária é necessária para tais problemas, pois o ângulo de ataque efetivo e as forças mudam dinamicamente com o movimento do corpo.

O artigo nota modestamente que, embora o método seja eficaz, existem limitações, particularmente em casos tridimensionais, onde restrições de resolução de malha podem levar a discrepâncias entre as formas otimizadas e binarizadas. Os autores sugerem que trabalhos futuros devem focar na melhoria da resolução da malha (ex: via métodos adaptativos) e na extensão do framework para fluxos turbulentos. O estudo conclui que a abordagem de otimização topológica acoplada proposta fornece uma ferramenta robusta para o design de componentes móveis impulsionados por fluidos, oferecendo insights sobre os mecanismos físicos de geração de sustentação e arrasto em regimes não estacionários.