Método de Celdas de Boltzmann (LBM)
Storyboard
El método de celdas de Boltzmann o lattice Boltzmann Model emplea un sistema de ecuaciones basados en la teoría de transporte de Boltzmann para calcular la velocidad de un fluido.
ID:(1030, 0)
Densidad
Ecuación
Si los parámetros se calculan con\\n\\n
$\chi = m c(\vec{x},t)$
y se promedia sobre la velocidad mediante
$ \chi_k(\vec{x},t) =\displaystyle\frac{1}{c(\vec{x},t)}\displaystyle\int d\vec{v} f(\vec{x},\vec{v},t) \chi_k(\vec{x},\vec{v},t)$ |
se obtiene mediante la masa la estimación de la densidad mediante:
$\rho(\vec{x},t) = m\displaystyle\int f(\vec{x},\vec{v},t)d\vec{v}$ |
ID:(8458, 0)
Velocidad de flujo
Ecuación
Si los parámetros se calculan con\\n\\n
$\chi_k = v_k$
promediando sobre la velocidad mediante
$ \chi_k(\vec{x},t) =\displaystyle\frac{1}{c(\vec{x},t)}\displaystyle\int d\vec{v} f(\vec{x},\vec{v},t) \chi_k(\vec{x},\vec{v},t)$ |
\\n\\ny con\\n\\n
$c(\vec{x},t)=\displaystyle\frac{1}{m}\rho(\vec{x},t)$
la velocidad del flujo se calcula integrando la función distribución de velocidad sobre todas las velocidades ponderando sobre las velocidades:
$\vec{u}(\vec{x},t) = \displaystyle\frac{m}{\rho}\int \vec{v}f(\vec{x},\vec{v},t)d\vec{v}$ |
ID:(8459, 0)
Temperatura
Ecuación
Con el teorema de equipartición en que\\n\\n
$\displaystyle\frac{1}{2}m\vec{v}\cdot\vec{v}=\displaystyle\frac{3}{2}k_B T$
\\n\\ncon el parámetro se calculan con\\n\\n
$\chi = T = \displaystyle\frac{m\vec{v}\cdot\vec{v}}{3k_B}=\displaystyle\frac{\vec{v}\cdot\vec{v}}{3R}\displaystyle\frac{c(\vec{x},t)}{\rho(\vec{x},t)}$
y se promedia promediando sobre la velocidad mediante
$ \chi_k(\vec{x},t) =\displaystyle\frac{1}{c(\vec{x},t)}\displaystyle\int d\vec{v} f(\vec{x},\vec{v},t) \chi_k(\vec{x},\vec{v},t)$ |
y se considera el teorema de equipartición, la temperatura se podrá estimar integrando la energía cinética ponderada por la distribución de velocidad dividida por la constante de los gases:
$T(\vec{x},t) = \displaystyle\frac{m}{3R\rho}\displaystyle\int (\vec{v}\cdot\vec{v})f(\vec{x},\vec{v},t)d\vec{v}$ |
ID:(8460, 0)
Tensor de tensión
Ecuación
Si los parámetros se calculan con\\n\\n
$\chi = m c(\vec{x},t)(v_i-u_i)(v_j-u_j)$
y se promedia sobre la velocidad mediante
$ \chi_k(\vec{x},t) =\displaystyle\frac{1}{c(\vec{x},t)}\displaystyle\int d\vec{v} f(\vec{x},\vec{v},t) \chi_k(\vec{x},\vec{v},t)$ |
el tensor del flujo se calcula integrando la función distribución de velocidad sobre todas las velocidades ponderando sobre las diferencias de velocidades:
$\sigma_{ij} = m\displaystyle\int (v_i-u_i)(v_j-u_j)f(\vec{x},\vec{v},t)d\vec{v}$ |
ID:(8461, 0)
Función de discretización
Ecuación
En el caso de la discretización en los modelos LBM se trabaja no con funciones de la velocidad si no que con componentes discretas. De esta forma se define la componente
$f_i(\vec{x},t)=w_if(\vec{x},\vec{v}_i,t)$ |
en donde
ID:(8466, 0)
Densidad del Gas
Ecuación
Con la descritización
$f_i(\vec{x},t)=w_if(\vec{x},\vec{v}_i,t)$ |
la ecuación
$\rho(\vec{x},t) = m\displaystyle\int f(\vec{x},\vec{v},t)d\vec{v}$ |
pasa a ser
$\rho(\vec{x},t)=m\sum_if_i(\vec{x},t)$ |
ID:(8492, 0)
Densidad de momento del gas
Ecuación
Con la descritización
$f_i(\vec{x},t)=w_if(\vec{x},\vec{v}_i,t)$ |
la ecuación
$\vec{u}(\vec{x},t) = \displaystyle\frac{m}{\rho}\int \vec{v}f(\vec{x},\vec{v},t)d\vec{v}$ |
pasa a ser
$\rho(\vec{x},t)\vec{u}(\vec{x},t)=m\sum_i\vec{e}_if_i(\vec{x},t)$ |
ID:(8493, 0)
Temperatura del Gas
Ecuación
Con la descritización
$f_i(\vec{x},t)=w_if(\vec{x},\vec{v}_i,t)$ |
la ecuación
$T(\vec{x},t) = \displaystyle\frac{m}{3R\rho}\displaystyle\int (\vec{v}\cdot\vec{v})f(\vec{x},\vec{v},t)d\vec{v}$ |
pasa a ser
$T(\vec{x},t)=\displaystyle\frac{m}{3R\rho}\sum_i(\vec{e}_i\cdot\vec{e}_i)f_i(\vec{x},t)$ |
ID:(8897, 0)
Ecuación de Boltzmann
Ecuación
La función de Boltzmann describe el transporte de un sistema de partículas descrito por la función de distribución de velocidades:
$\displaystyle\frac{\partial f}{\partial t}+v_i\displaystyle\frac{\partial f}{\partial x_i}=C(f)$ |
donde el termino
ID:(8462, 0)
Teoría de Grad de los 13 momentos
Ecuación
La distribución de velocidades se puede representar como un polinomio ortogoan de Hermite
$f^N(\vec{x},\vec{v},t)=\omega(\vec{v})\displaystyle\sum_{n=0}^N\displaystyle\frac{1}{n!}a^{(n)}(\vec{x},t)\mathcal{H}(\vec{v})$ |
ID:(8463, 0)
Coeficiente de Orden 0
Ecuación
El coeficiente de orden cero es
$a^{(0)}=\displaystyle\int f(\vec{x},\vec{v},t)d\vec{v}=\rho$ |
ID:(8464, 0)
Coeficiente de Orden 1
Ecuación
El coeficiente de orden cero es
$a^{(1)}=\displaystyle\int \vec{v}f(\vec{x},\vec{v},t)d\vec{v}=\rho\vec{u}$ |
ID:(8465, 0)
Coeficientes
Ecuación
Los coeficientes son:
$a^{(n)}(\vec{x},\vec{v},t)=\displaystyle\int f(\vec{x},\vec{v},t)\mathcal{H}^{(n)}(\vec{v})d\vec{v}=\displaystyle\sum_if_i(\vec{x},t)\mathcal{H}^{(n)}(\vec{v})$ |
ID:(8467, 0)