El método de Monte Carlo y la solución de la ecuación de transferencia radiativa

Abstract

El paso de la radiación que atraviesa un medio donde ocurren los fenómenos de absorción, dispersión y emisión, es modelado por la ecuación de transferencia radiativa (ETR). Para medios puramente absorbentes, el fenómeno de dispersión es despreciado reduciendo la ETR es una ecuación diferencial no lineal, sin solución analítica y su solución es aproximada por métodos numéricos determinísticos. Monte Carlo es un método estocástico, basado en probabilidades y utiliza un generador de números aleatorios para simular la ocurrencia de un evento, en este caso, el paso de la radiación por un medio. En este trabajo se utilizó el método de diferencias finitas junto con el de ordenadas discretas para resolver la ETR en una dimensión, los resultados fueron comparados con los obtenidos por el método de Monte Carlo. Algunos ejemplos de prueba son presentados, así también los programas utilizados para la solución de ambos métodos en el SciLab.

Keywords

• Ecuación de trasferencia radiativa
• método de diferencias finitas
• método de ordenadas discretas
• método de Monte Carlo
• SciLab