Monte Carlo, espalhamento da radiação, método dos volumes finitos, transferência de calor tridimensional transiente
O câncer é alvo de intensa pesquisa para aprimorar procedimentos terapêuticos. Figuram nesse cenário terapias de hipertermia, como a Terapia a Laser Intersticial (ILT) e a Terapia Fototérmica Plasmônica (PPTT). Modelos são empregados na simulação computacional de tratamentos hipertérmicos, sendo crucial validá-los para assegurar confiabilidade. Este estudo valida modelos de hipertermia ao comparar distribuições de temperatura de simulações às obtidas em um protocolo real de ILT para câncer de mama. A irradiação por um laser é administrada através de uma fibra óptica no centro de um tumor aproximadamente esférico. O protocolo inclui dispositivos para o controle da energia entregue ao tecido, como uma bomba de fluido, evitando danos à ponta da fibra óptica, e sondas com termopares para medir temperaturas. A simulação abrange a resolução numérica da Equação do Biocalor (PBHTM) e da Equação de Transferência Radiativa (RTE), sendo o divergente do fluxo de calor, obtido pela RTE, o termo fonte na PBHTM. A solução numérica fornece distribuições de temperatura na região irradiada, abrangendo o tumor e os tecidos circundantes. Propriedades termofísicas e ópticas dos tecidos foram obtidas na literatura, enquanto as do laser foram obtidas de um manual técnico de operação do laser DIOMED DELTA15. As propriedades do meio foram assumidas como as de um tecido mamário saudável, em um domínio de cálculo cilíndrico de 4 cm de comprimento e 2 cm de raio. A ponta da fibra óptica foi localizada no centro do domínio. Baseando-se no protocolo, o tempo total de irradiação foi de 2100 s, considerando um laser semicondutor de diodo com comprimento de onda de 805 nm, operando a 5 W e diâmetro de feixe de luz de 600 $\mu$m. Resultados foram obtidos para averiguar a sensibilidade do campo de temperaturas em relação à parâmetros associados a incertezas, incluindo regime de operação do laser (pulsado ou contínuo), posicionamento dos termopares e abrangência da região de arrefecimento na ponta da fibra ótica. Os resultados obtidos através das simulações apresentaram excelente concordância com os experimentais provenientes do protocolo, conferindo confiabilidade à modelagem, hipóteses adotadas e códigos utilizados.