Notas sobre Transferência de Calor e Mecânica dos Fluidos Computacional

 Durante a implementação de um código para modelar a transferência de calor por radiação em um domínio bidimensional, me deparei com o seguinte problema: mesmo levando em conta as simplificações, a grandeza fundamental do problema, a intensidade da radiação, precisaria ser descrita por quatro dimensões, ou seja, é uma variável quadridimensional. As dimensões são: duas coordenadas espaciais: r e z (o problema é em coordenadas cilíndricas) e duas coordenadas que definem a direção de propagação da radiação: o ângulo polar theta e o ângulo azimutal phi . Não é incomum que o modelo numérico de problemas de CFD e CHT produzam matrizes esparsas, ou seja, matrizes onde apenas a diagonal principal e as diagonais adjacentes são formadas por elementos não nulos. Neste caso, o uso de matrizes para computar a relação entre um ente discreto (ponto, volume ou elemento) e os vizinhos faz com que seja necessário alocar muita memória que não será usada. Além de limitar o tamanho do problema que o...