Entendendo ondas mesmo
A maioria dos materiais didáticos trata ondas mecânicas e ondas eletromagnéticas como se fossem categorias estanques, mas na prática elas aparecem juntas em situações que ninguém avisa antes. Se você está mexendo com simulação de antenas, acústica de salas ou até medições de campo EM em laboratório, acabar enfrentando os dois tipos ao mesmo tempo é padrão. Vou explicar como eu lido com isso.
Ondas mecanicas eletromagneticas: o que separa uma da outra
Ondas mecânicas precisam de um meio material para se propagar. Sonoridade em ar, vibrações em barras metálicas, ondas em cordas tensionadas. Sem meio, simplesmente não existe. Ondas eletromagnéticas, por outro lado, se propagam no vácuo porque são oscilações acopladas dos campos elétrico e magnético. Isso já explica metade da diferença. A outra metade tá nos números. No ar, o som viaja a aproximadamente 343 metros por segundo a 20 graus Celsius. Luz no vácuo viaja a 299.792.458 metros por segundo. A ordem de grandeza entre elas é de cerca de nove casas decimais. Isso muda tudo: o comprimento de onda para uma mesma frequência é completamente diferente, e o tipo de simulação ou medição que você vai precisar também.
Uma onda sonora de 1 kHz tem comprimento de onda de cerca de 34 cm no ar. Uma onda eletromagnética de 1 kHz, se existisse como radiação livre, teria comprimento de onda de 300 km. Frequências diferentes geram comportamentos completamente distintos com a mesma geometria. Isso é o erro mais comum que eu vejo em quem tá começando: calcular o comprimento de onda errado porque confundiu qual velocidade usar.
Como escolher o método certo
Se o seu problema envolve propagação em materiais homogêneos e contínuos, método dos elementos finitos (MEF) resolve bem. Para ondas eletromagnéticas em geometrias complexas com condutores, método dos momentos (MoM) costuma ser mais eficiente. Quando o meio é irregular — tipo tecidos biológicos em imageamento médico ou solo com camadas heterogêneas — FDTD (diferenças finitas no domínio do tempo) se torna a opção padrão. Não é questão de preferência, é questão de escala do problema. Para ondas mecânicas, MEF domina porque lida naturalmente com condições de contorno estruturais. Já pensei em usar FDTD para vibração em vigas, mas a dispersão numérica do método distorce demais o resultado em frequências altas. Acabei abandonando essa abordagem há anos.
Se você quer uma implementação prática e rápida, o código aberto GetFEM++ roda simulações 3D de acústica com MEF de forma relativamente direta. Para eletromagnetismo, o Meep (do MIT) usa FDTD e funciona bem para problemas de difração e guiamento. Ambos exigem que você configure malha adequadamente, mas os exemplos de uso são acessíveis. Eu costumava usar esses dois pacotes juntos em projetos de transdutores ultrassônicos, onde o domínio elétrico e o domínio mecânico precisavam ser acoplados numericamente.
👉 Clique no botão abaixo para saber mais sobre o assunto!
O problema que ninguém conta sobre interface entre os dois domínios
Eu passei três semanas tentando fazer uma simulação de um transdutor piezoelétrico que convertesse sinal elétrico em onda mecânica (ultrassom) e depois recuperasse o sinal de eco. O modelo parecia correto nas equações, mas os resultados numéricos ficavam instáveis sempre que a frequência passava de 2 MHz. A solução não era nenhuma equação nova. Era um problema de escala de malha. A malha eletromagnética precisava de elementos na casa de micrômetros porque o comprimento de onda do campo EM no cerâmico piezoelétrico era da ordem de milímetros divididos pela constante dielétrica elevada. Já a malha mecânica precisava resolver ondas acústicas no material de acoplamento, que tinha velocidade bem menor. Quando eu usava uma única malha para tudo, o condition number da matriz travava o solver. A solução foi particionar: malha fina no domínio piezoelétrico, malha mais grosseira no meio de acoplamento, e acoplar as fronteiras com condições de transmitância. Isso reduziu o tempo de simulação de cerca de 6 horas para 40 minutos em minha estação de trabalho.
Se você está fazendo algo similar, use solver iterativo com pré-condicionador algébrico (ILU ou AMG) em vez de solver direto. A diferença em memória é brutal. Eu tive que matar processos porque o solver direto consumia mais de 128 GB de RAM em um problema que o solver iterativo resolveu com 16 GB.
Pegadinhas que custam tempo
Condições de contorno absorventes (PML — perfectly matched layer) funcionam bem para eletromagnetismo, mas no domínio mecânico elas são muito mais sensíveis à direção de incidência. Se a onda sonora vem em ângulo oblíquo forte, a PML pode refletir parte do sinal de volta para o domínio computacional. A correção é usar PML anisotrópica ou, se o custo permitir, aumentar a espessura da camada absorvente para pelo menos um comprimento de onda na frequência mais baixa do interesse. Outro ponto: a Impedância característica do meio determina quanta energia é refletida na interface entre dois materiais. No caso de ondas mecânicas, Z = c (densidade vezes velocidade do som). No caso eletromagnético, Z = (/). Quando há grande descontinuidade de impedância — como ar para aço — a maior parte da energia é refletida. Isso parece óbvio, mas na prática leva a erros de modelagem porque as pessoas esquecem de incluir camadas de adaptação de impedância na simulação e acabam interpremando reflexões numéricas como se fossem fenômeno físico real.
Quando não adianta simular
Se o meio tem heterogeneidades na escala do comprimento de onda, métodos determinísticos perdem eficiência rapidamente. Espalhamento múltiplo em materiais granulares, concreto com agregados variados, tecidos biológicos com microestrutura complexa — nesses casos, abordagem estocástica ou homogeneização efetiva costuma ser mais viável. Simular tudo explicitamente pode exigir milhões de elementos e ainda assim não capturar a variabilidade real. Nesses cenários, modelos fenomenológicos de atenuação e velocidade efetiva, calibrados com medições, dão resultado útil com uma fração do esforço computacional. Também não adianta usar eletromagnetismo de alta frequência (aproximação óptica) quando as dimensões do problema são menores que o comprimento de onda. A aproximação vale quando ka >> 1, onde k é o número de onda e a é a dimensão característica. Abaixo disso, você precisa da solução completa de Maxwell. Já vi gente aplicar ray tracing para antenas de tamanho comparável ao comprimento de onda e se perguntar por que os resultados não batiam com a medida. O problema não era o software, era o regime físico inadequado.
Referências práticas
Para ondas mecânicas, o livro do Kinsler e Frey ainda é a base mais confiável, apesar de antigo. A parte sobre impedância acústica e reflexão em interfaces segue sendo referenciada em papers recentes sem alterações. Para eletromagnetismo, o Jin do IEEE Press tem a formulação variacional mais clara que eu já vi para quem precisa implementar. Se você quer começar com algo que rode hoje, procure por tutoriais do GetFMM++ para acústica e do Meep para óptica/eletromagnetismo. Ambos têm examples compiláveis e a documentação cobre os casos mais comuns. Eu configurei meu primeiro modelo de guided wave em placas usando o Meep seguindo o tutorial de grating coupler, e em duas horas tinha algo queproduzia resultados compatíveis com a literatura.
