Aplicación al calculo de Dosis
Storyboard
Como la ecuación de Boltzmann incluye las colisiones entre partículas, es posible estimar la energía entregada al tejido y con ello la dosis en radioterapia. La mayor eficiencia se logra al trabajar con distribuciones de partículas (fotones, electrones y positrones) en vez de las partículas individuales como en el caso de Monte Carlo. La introducción de celdas mediante el método de celdas de Boltzmann debiese ademas permitir la posibilidad de elegir la precisión en función de los recursos computacionales que se ponen a disposición. Esto significa que el algoritmo debiese poder pasar en forma continua desde un sistema equivalente a Pencil Beam, Convolución hasta Monte Carlo.
ID:(1164, 0)
Caso Photones
Ecuación
Para el caso en que se consideran fotones térmicos uniformemente distribuidos su número por celda será según la distribución de Bose-Einstein
$f_i^{eq}=\displaystyle\frac{1}{e^{\hbar\omega/kT}-1}$ |
donde
ID:(8561, 0)
Caso Electrones
Ecuación
En el caso de los electrones se tiene que la distribución en equilibrio es la de Fermi-Direca por lo que en la situación de equilibrio la distribución tendrá a ser de la forma
$f^{eq}_i=\displaystyle\frac{1}{e^{\beta (m_ev_i^2/2-\mu)}+1}$ |
Fuera de ello los posibles scatterings corresponden a
- Absorción
- Choques elásticos
- Colisión electrón-electrón
- Excitación y deexitación
ID:(9165, 0)
Scattering
Imagen
Los scattering que contribuyen (in) o describen el abandono de partículas (out) se pueden graficar de la siguiente forma:
Gráfica Scattering entre dos partículas
Hay que notar que el termino colisión:
- integra sobre todas las velocidades externas a las del volumen
- incluye la probabilidad de que existan ambas velocidades que llevan al scattering simultaneamente
- la velocidad relativa multiplicado por la sección eficaz total representa el flujo de partículas hacia el target
Esto ultimo se puede mostrar en forma simple mediante
ID:(9177, 0)
Colisiones
Ecuación
En caso de que las partículas colisionan la función distribución
$\displaystyle\frac{df}{dt}\neq 0$
Las colisiones lleva a que partículas de celdas vecinas sufran un colisión que las lleva a la celda en consideración y partículas dentro de la celda ser expulsadas. Lo primero lleva a un incremento de partículas
$\displaystyle\frac{df}{dt}=\displaystyle\frac{1}{\tau}(f_{in}-f_{out})$ |
ID:(9077, 0)
Colisiones que abandonan la celda
Ecuación
En el caso que abandonan la celda se considera
$f(\vec{x},\vec{v}_1,t)f(\vec{x},\vec{v}_2,t)|\vec{v}_2-\vec{v}_1|\sigma(\vec{v}_1,\vec{v}_2\rightarrow\vec{v}_12,\vec{v}_22)d\vec{v}_12d\vec{v}_22$ |
integrando sobre una de las velocidades que inician la colisión y ambas resultantes ya que la otra es la contribución a la función distribución local
$\displaystyle\frac{1}{\tau}f_{out}(\vec{v})=\displaystyle\int d\vec{v}_1d\vec{v}_12d\vec{v}_22f(\vec{x},\vec{v}_1,t)f(\vec{x},\vec{v},t)|\vec{v}-\vec{v}_1|\sigma(\vec{v},\vec{v}_1\rightarrow\vec{v}_12,\vec{v}_22)$ |
ID:(9080, 0)
Colisiones que contribuyen
Ecuación
En el caso de contribuciones a la celda se considerar
$f(\vec{x},\vec{v}_1,t)f(\vec{x},\vec{v}_2,t)|\vec{v}_2-\vec{v}_1|\sigma(\vec{v}_1,\vec{v}_2\rightarrow\vec{v}_12,\vec{v}_22)d\vec{v}_12d\vec{v}_22$ |
integrando sobre las velocidades que inician la colisión y una de las resultantes ya que la otra es la contribución a la función distribución local
$\displaystyle\frac{1}{\tau}f_{in}(\vec{v})=\displaystyle\int d\vec{v}_1d\vec{v}_2d\vec{v}_12f(\vec{x},\vec{v}_1,t)f(\vec{x},\vec{v}_2,t)|\vec{v}_2-\vec{v}_1|\sigma(\vec{v}_1,\vec{v}_2\rightarrow\vec{v}_12,\vec{v})$ |
ID:(9079, 0)