Limite continuo
Storyboard
En base a las distintas propiedades que se pueden calcular con las distribuciones de velocidades se pueden también establecer ecuaciones diferenciales continuas que equivalen a las conocidads ecuación de continuidad, Navier Stokes y transporte de calor.
ID:(1222, 0)
Introducción de la aceleración
Ecuación
Como la masa es constante, la fuerza con es
$\displaystyle\frac{\partial f}{\partial t}+\vec{v}\cdot\nabla_{\vec{x}} f+\displaystyle\frac{1}{m}\vec{F}\cdot\nabla_{\vec{v}} f=0$ |
puede ser reemplazada por la masa por la aceleración quedando con la ecuación de Boltzmann de la forma
$\displaystyle\frac{\partial f}{\partial t}+\displaystyle\sum_iv_i\displaystyle\frac{\partial f}{\partial x_i}+\displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=\displaystyle\frac{df}{dt}\mid_c$ |
en que se incluyo un termino que representa las colisiones.
ID:(9709, 0)
Calculo de la densidad
Ecuación
Si se multiplica la ecuación con
$\displaystyle\frac{\partial f}{\partial t}+\displaystyle\sum_iv_i\displaystyle\frac{\partial f}{\partial x_i}+\displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=\displaystyle\frac{df}{dt}\mid_c$ |
por la masa y se integra sobre las velocidades se obtiene con que
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
ID:(9710, 0)
Calculo de la densidad, primer termino
Ecuación
El primer termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
puede reescribirse con la expresión para la densidad con
$\rho(\vec{x},t) = m\displaystyle\int f(\vec{x},\vec{v},t)d\vec{v}$ |
con como
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v f=\displaystyle\frac{\partial\rho}{\partial t}$ |
ID:(9711, 0)
Calculo de la densidad, segundo termino
Ecuación
El primer termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
puede reescribirse con la expresión para la velocidad media con
$\vec{u}(\vec{x},t) = \displaystyle\frac{m}{\rho}\int \vec{v}f(\vec{x},\vec{v},t)d\vec{v}$ |
con como
$m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if=\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho u_i)=\vec{\nabla}\cdot(\rho\vec{u})$ |
ID:(9712, 0)
Calculo de la densidad, tercer termino
Ecuación
El tercer termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
\\n\\npuede reescribirse como\\n\\n
$m\displaystyle\int_V d^3v \displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int_V d^3v\vec{\nabla}\cdot(\vec{a}f)$
\\n\\nlo que se puede reducir, mediante el teorema de la divergencia, a una integral sobre la superficie \\n\\n
$m\displaystyle\int_V d^3v\vec{\nabla}\cdot(\vec{a}f)=m\displaystyle\int_S dS\hat{n}\cdot(\vec{a}f)$
Como se puede asumir que la velocidad en el infinito es cero se tendrá con que el tercer termino es nulo:
$m\displaystyle\int d^3v \displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=0$ |
ID:(9713, 0)
Calculo de la densidad, termino colisiones
Ecuación
Como la masa se conserva, el termino de colisiones en la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
debe ser nulo, o sea con
$m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c=0$ |
ID:(9714, 0)
Ecuación de continuidad
Ecuación
Partiendo de las expresiones, con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v f=\displaystyle\frac{\partial\rho}{\partial t}$ |
con
$m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if=\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho u_i)=\vec{\nabla}\cdot(\rho\vec{u})$ |
con
$m\displaystyle\int d^3v \displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=0$ |
y con
$m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c=0$ |
la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3vf+m\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_if+m\displaystyle\sum_i\displaystyle\int d^3v a_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v \displaystyle\frac{df}{dt}\mid_c$ |
se reduce a la ecuación de continuidad con :
$\displaystyle\frac{\partial\rho}{\partial t}+\vec{\nabla}\cdot(\rho\vec{u})=0$ |
ID:(9715, 0)
Calculo del momento
Ecuación
Si se multiplica la ecuación con
$\displaystyle\frac{\partial f}{\partial t}+\displaystyle\sum_iv_i\displaystyle\frac{\partial f}{\partial x_i}+\displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=\displaystyle\frac{df}{dt}\mid_c$ |
por la masa y componente
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
ID:(9716, 0)
Calculo del momento, primer termino
Ecuación
El primer termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
puede reescribirse con la expresión para la velocidad media con
$\vec{u}(\vec{x},t) = \displaystyle\frac{m}{\rho}\int \vec{v}f(\vec{x},\vec{v},t)d\vec{v}$ |
con como
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf=\displaystyle\frac{\partial}{\partial t}(\rho u_j)$ |
ID:(9717, 0)
Calculo del momento, segundo termino
Ecuación
El segundo termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
involucra un tensor del promedio del producto de las velocidades con :
$\rho\langle v_jv_i\rangle\equiv m\displaystyle\int d^3v v_jv_if$ |
ID:(9718, 0)
Velocidad de las partículas en el sistema local
Ecuación
La velocidad
$\vec{w}=\vec{v}-\vec{u}$ |
ID:(9720, 0)
Valor medio de la velocidad de las partículas en el sistema local
Ecuación
Por definición de la velocidad media de las partículas es igual a la velocidad del flujo:\\n\\n
$\langle\vec{v}\rangle=\vec{u}$
y por ello la velocidad de las partículas en el sistema local con es
$\langle\vec{w}\rangle=\langle\vec{v}\rangle-\langle\vec{u}\rangle=\langle\vec{v}\rangle-\vec{u}=0$ |
ID:(9721, 0)
Tensor del promedio del producto de las velocidades
Ecuación
Con la definición de la velocidad en el sistema de coordenadas locales
$\vec{w}=\vec{v}-\vec{u}$ |
se puede reescribir el tensor del promedio del producto de las velocidades con como
$m\displaystyle\int d^3v v_jv_if=\rho\langle v_jv_i\rangle=\rho u_ju_i + \rho\langle w_jw_i\rangle$ |
ID:(9722, 0)
Calculo del momento, tercer termino
Ecuación
El tercer termino de la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
\\n\\ncorresponde a la derivada temporal de la densidad por el producto de la velocidad y la derivada de la función de distribución\\n\\n
$m\displaystyle\int d^3v \displaystyle\sum_i v_ja_i\displaystyle\frac{\partial f}{\partial v_i}$
\\n\\nen que se puede realizar una integración parcial quedando\\n\\n
$m\displaystyle\sum_i a_i\displaystyle\int d^3v v_j\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\sum_i a_i\displaystyle\int d^3v\left[\displaystyle\frac{\partial}{\partial v_i}(fv_i)-\delta_{ij}f\right]$
\\n\\nEn este caso se puede realizar la integración en que el primer termino es nulo en la medida que la distribución lo sea en el infinito, por lo que se tiene\\n\\n
$m\displaystyle\sum_i a_i\displaystyle\int d^3v\left[\displaystyle\frac{\partial}{\partial v_i}(fv_i)-\delta_{ij}f\right]=-a_jm\displaystyle\int d^3v f$
que con la densidad con
$\rho(\vec{x},t) = m\displaystyle\int f(\vec{x},\vec{v},t)d\vec{v}$ |
resulta con
$m\displaystyle\int d^3v \displaystyle\sum_i v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=-\rho a_j$ |
ID:(9719, 0)
Calculo del momento, termino colisiones
Ecuación
Como el momento se conserva en las colisiones, el termino de colisiones en la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
debe ser nulo, o sea con
$m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c=0$ |
ID:(9723, 0)
Ecuación de momento
Ecuación
De las expresiones, con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf=\displaystyle\frac{\partial}{\partial t}(\rho u_j)$ |
con
$m\displaystyle\int d^3v v_jv_if=\rho\langle v_jv_i\rangle=\rho u_ju_i + \rho\langle w_jw_i\rangle$ |
con
$m\displaystyle\int d^3v \displaystyle\sum_i v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=-\rho a_j$ |
y con
$m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c=0$ |
la ecuación con
$m\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v v_jf+m\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v v_jv_if+m\displaystyle\sum_i\displaystyle\int d^3v v_ja_i\displaystyle\frac{\partial f}{\partial v_i}=m\displaystyle\int d^3v v_j\displaystyle\frac{df}{dt}\mid_c$ |
se reduce a la ecuación de continuidad con :
$\displaystyle\frac{\partial}{\partial t}(\rho u_j)+\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho u_ju_i + \rho\langle w_jw_i\rangle)-\rho a_j=0$ |
ID:(9724, 0)
Ecuación de flujo
Ecuación
La ecuación con
$\displaystyle\frac{\partial}{\partial t}(\rho u_j)+\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho u_ju_i + \rho\langle w_jw_i\rangle)-\rho a_j=0$ |
se puede reescribir con la ecuación de continuidad con
$\displaystyle\frac{\partial\rho}{\partial t}+\vec{\nabla}\cdot(\rho\vec{u})=0$ |
con como:
$\rho\displaystyle\frac{\partial}{\partial t}\vec{u}+\vec{u}\cdot\nabla\vec{u}+\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho\langle w_jw_i\rangle)=\rho\vec{a}$ |
ID:(9725, 0)
Tensor de tensión viscosa
Ecuación
Para escribir la ecuación del flujo en la forma tradicional se puede introducir con el tensor de tensión viscosa definido con la presión:
$p=\displaystyle\frac{1}{3}\rho\displaystyle\sum_i\langle w_i^2\rangle$ |
con por
$\pi_{ij}=p\delta_{ij}-\rho\langle w_iw_j\rangle$ |
ID:(9727, 0)
Ecuación de flujo en forma tensorial
Ecuación
Con el tensor de tensión viscosa con :
$\pi_{ij}=p\delta_{ij}-\rho\langle w_iw_j\rangle$ |
la ecuación de flujo con
$\rho\displaystyle\frac{\partial}{\partial t}\vec{u}+\vec{u}\cdot\nabla\vec{u}+\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}(\rho\langle w_jw_i\rangle)=\rho\vec{a}$ |
se escribe con como
$\rho\displaystyle\frac{\partial}{\partial t}\vec{u}+\vec{u}\cdot\nabla\vec{u}+\vec{\nabla} p-\displaystyle\sum_i\displaystyle\frac{\partial}{\partial x_i}\pi_{ij}=\rho\vec{a}$ |
ID:(9728, 0)
Ecuación de energía
Ecuación
Si se multiplica la ecuación con
$\displaystyle\frac{\partial f}{\partial t}+\displaystyle\sum_iv_i\displaystyle\frac{\partial f}{\partial x_i}+\displaystyle\sum_ia_i\displaystyle\frac{\partial f}{\partial v_i}=\displaystyle\frac{df}{dt}\mid_c$ |
por la energía cinética y se integra sobre las velocidades se obtiene que con
$\displaystyle\frac{\partial}{\partial t}\displaystyle\int d^3v \displaystyle\frac{1}{2}mv^2f+\displaystyle\sum_i \displaystyle\frac{\partial}{\partial x_i}\displaystyle\int d^3v \displaystyle\frac{1}{2}mv^2v_if+\displaystyle\sum_i\displaystyle\int d^3v \displaystyle\frac{1}{2}mv^2a_i\displaystyle\frac{\partial f}{\partial v_i}=\displaystyle\int d^3v \displaystyle\frac{1}{2}mv^2\displaystyle\frac{df}{dt}\mid_c$ |
ID:(9729, 0)
Flujo de energía por conducción
Ecuación
El flujo de energía por conducción se puede calcular con de
$\vec{F}=\displaystyle\frac{1}{2}\langle\vec{w}(\vec{w}\cdot\vec{w})\rangle$ |
ID:(9731, 0)
Fracción de disipación viscosa
Ecuación
La fracción de disipación viscosa se puede calcular con de
$\Psi=\displaystyle\sum_{i,j}\pi_{ij}\displaystyle\frac{\partial u_i}{\partial x_j}$ |
ID:(9732, 0)
Ecuación de energía
Ecuación
La energía se puede calcular con
$\rho\displaystyle\frac{\partial\epsilon}{\partial t}+\rho\vec{u}\cdot\vec{\nabla}\epsilon=-p\vec{\nabla}\cdot\vec{u}-\vec{\nabla}\cdot\vec{F}+\Psi$ |
ID:(9733, 0)