Efficient 2.5-D FEM-Based Scattering Analysis of the Human Body for RF Sensing

Título: Análise Eficiente de Espalhamento do Corpo Humano para Sensoriamento RF Baseada em FEM 2.5-D

1. Problema e Motivação

O sensoriamento de Radiofrequência (RF) e a Localização sem Dispositivos (DFL) dependem fortemente de grandes conjuntos de dados para o treinamento de modelos de aprendizado de máquina. No entanto, a coleta de dados reais é cara, demorada e difícil de escalar, exigindo medições extensas em múltiplos cenários, poses corporais e condições ambientais.

• Limitações das abordagens existentes:
  • Métodos de Onda Completa (3-D FEM, FDTD, MoM): Oferecem alta precisão, mas exigem tempo computacional e recursos excessivos, tornando-os inviáveis para a geração de grandes bancos de dados sintéticos.
  • Modelos Analíticos/Semi-analíticos: Usam geometrias simples (esferas, cilindros infinitos), oferecendo eficiência, mas com realismo severamente limitado na representação das características de espalhamento humano.
  • Aproximações de Alta Frequência: Eficientes para propagação, mas falham na precisão do espalhamento de campo próximo e em excitações fora do eixo.
    O desafio é encontrar um equilíbrio entre eficiência computacional e precisão física para simular a interação complexa entre o corpo humano e campos eletromagnéticos em cenários internos.

2. Metodologia Proposta

O artigo propõe um método rápido de Elementos Finitos 2.5-Dimensional (2.5-D FEM) para calcular os campos espalhados por um modelo humano aproximado como um Corpo de Revolução (BoR - Body of Revolution) sob excitação de um dipolo orientado no eixo z.

• Aproximação BoR: O corpo humano é modelado como um corpo de revolução, preservando características essenciais de espalhamento enquanto explora a simetria rotacional. Isso permite desacoplar as contribuições de diferentes harmônicos, reduzindo drasticamente o tempo de cálculo e o uso de memória em comparação com solvers 3-D completos.
• Formulação Matemática:
  • O campo elétrico espalhado é expandido em uma soma de harmônicos circulares no ângulo azimutal ($\phi$) e resolvido no plano $(\rho, z)$.
  • Utiliza-se o método de Galerkin com uma base mista de elementos finitos (funções de aresta para componentes no plano e funções nodais de Lagrange para a componente fora do plano).
  • O sistema de equações resultante possui uma estrutura bloqueio-diagonal, permitindo a solução independente de cada harmônico, o que acelera a computação.
• Tratamento de Fontes e Fronteiras:
  • A fonte é um dipolo de Hertz orientado em z, posicionado fora do domínio de simulação.
  • Camadas Perfeitamente Casadas (PML) são usadas para simular o espaço livre.
  • Para calcular campos fora do domínio de malha FEM (onde estão os receptores), utiliza-se o Teorema da Equivalência: correntes elétricas e magnéticas equivalentes são calculadas em uma superfície cilíndrica ao redor do corpo e usadas para reconstruir o campo em qualquer ponto externo.
  • O efeito do solo é modelado aproximadamente através de uma fonte imagem com um coeficiente de reflexão ($\Gamma$), lidando com a incerteza dos materiais do chão em ambientes internos.

3. Principais Contribuições

1. Formulação 2.5-D FEM para Dipolo Excitado: Desenvolvimento de uma formulação específica para excitação por dipolo em modelos BoR humanos. Reduz a complexidade computacional mantendo alta precisão numérica, sendo ideal para simulações em larga escala.
2. Validação para Previsão de RSSI: Validação da eficácia do método na previsão de distribuições espaciais de Intensidade do Sinal Recebido (RSSI) através da comparação com simulações de onda completa (Feko) e medições experimentais em um cenário DFL real.
3. Ferramenta de Geração de Dados Sintéticos: Demonstração de que o método pode gerar dados sintéticos rotulados com propriedades estatísticas consistentes com medições reais, incluindo a análise de micro-movimentos humanos (respiração, pequenos deslocamentos) e seu impacto nas estatísticas do RSSI.

4. Resultados e Validação

• Validação Numérica: A comparação com o software comercial Feko (usando Método dos Momentos) em um modelo de cilindro dielétrico mostrou excelente concordância nos padrões de diretividade, validando a precisão do método 2.5-D.
• Eficiência Computacional: O método 2.5-D FEM foi implementado em MATLAB. Para um cenário complexo, o tempo de solução foi de aproximadamente 60 segundos em um processador Intel i9, enquanto a simulação equivalente em Feko levaria várias horas. A aceleração é estimada em pelo menos um fator de 100 em comparação com a discretização 3-D completa.
• Comparação Experimental: Em um cenário DFL interno com 20 nós de rádio (2.43 GHz), os valores de RSSI simulados (após calibração escalar) apresentaram boa concordância com as medições reais em termos de tendências espaciais e atenuação relativa.
• Análise de Micro-movimentos: O estudo mostrou que, embora o RSSI absoluto seja pouco sensível a pequenos deslocamentos, a diferença de RSSI ($\Delta$RSSI) e as distribuições de probabilidade (PDF) capturam bem as flutuações estatísticas causadas por micro-movimentos. O método conseguiu reproduzir as tendências centrais e a dispersão das medições experimentais.

5. Significado e Conclusão

O método proposto oferece uma solução prática e eficiente para a geração de grandes bancos de dados de RF rotulados, um gargalo crítico para o desenvolvimento de sistemas de localização e sensoriamento baseados em aprendizado de máquina.

• Impacto: Permite a criação de modelos físicos-prioritários (physics-prior) mais robustos para DFL, acelerando a calibração e o ajuste de métodos de reconstrução tomográfica.
• Futuro: Os autores planejam estender o método para antenas direcionais e ambientes com múltiplos alvos (interação entre vários corpos humanos), visando simulações ainda mais realistas.

Em resumo, o trabalho preenche a lacuna entre a precisão física e a eficiência computacional, fornecendo uma ferramenta viável para simular cenários complexos de sensoriamento humano sem a necessidade de campanhas de medição exaustivas.