Processing math: 0%
User: No user logged in.


Ecuación de Fokker Planck

Storyboard

Para describir procesos irreversibles se puede estudiar como kas distribuciones de probabilidades evolucionan en en tiempo afectando la posibilidad de que un sistema vuelva a su estado inicial.

>Model

ID:(1140, 0)



Probabilidad de transición

Definition

Si se supone que una partícula tiene en un momento t_0 una velocidad v_0 existe una probabilidad\\n\\n

P(t,v|t_0,v_0)

que en un tiempo t posterior tenga una velocidad v.

Nota: Un proceso en que el nuevo estado depende del estado anterior se denomina un proceso de Markoff.

ID:(9129, 0)



Ecuación de Fokker Planck

Storyboard

Para describir procesos irreversibles se puede estudiar como kas distribuciones de probabilidades evolucionan en en tiempo afectando la posibilidad de que un sistema vuelva a su estado inicial.

Variables

Symbol
Text
Variable
Value
Units
Calculate
MKS Value
MKS Units
k_B
k_B
Constante de Boltzmann
J/K
m
m
Masa de la partícula
kg
M_n
M_n
Momento n-esimo de la distribución
-
n
n
Orden del momento de la distribución
-
M_1
M_1
Primer momento de la distribución
-
s
s
Primer tiempo relativo
s
P(v,s|v_0)
P_vsv0
Probabilidad de transición
-
M_2
M_2
Segundo momento de la distribución
-
\tau
tau
Segundo tiempo relativo
s
T
T
Temperatura
K
t
t
Tiempo
s
t_0
t_0
Tiempo inicial
s
v(\tau)
v_tau
Velocidad al final del segundo intervalo
m/s
v(0)
v0
Velocidad al principio del segundo intervalo
m/s
v
v
Velocidad final de la evolución
m/s
v_0
v_0
Velocidad inicial de la evolución
m/s
v_2
v_2
Velocidad intermedia de la evolución
m/s
\xi
xi
Velocidad relativa
m/s

Calculations


First, select the equation:   to ,  then, select the variable:   to 
s=t-t_0\displaystyle\frac{\partial P}{\partial s}dv\tau=-\displaystyle\int_{-\infty}^{\infty}du P(v,s|v_0)P(u,\tau|v)+\displaystyle\int_{-\infty}^{\infty}du P(u,s|v_0)P(v,\tau|u)\displaystyle\int_{-\infty}^{\infty}du P(u,s|v)=1\displaystyle\frac{\partial P}{\partial s}\tau=- P(v,s|v_0)+\displaystyle\int_{-\infty}^{\infty}d\xi P(v-\xi,s|v_0)P(v,\tau|v-\xi)\displaystyle\frac{\partial}{\partial s}P(v,s|v_0)=\sum_{n=1}^{\infty}\displaystyle\frac{(-1)^n}{n!}\displaystyle\frac{\partial^n}{\partial v^n}[M_nP(v,s|v_0)]M_n=\displaystyle\frac{1}{\tau}\displaystyle\int_{-\infty}^{\infty}d\xi P(v+\xi,\tau|v)\xi^n=\displaystyle\frac{1}{\tau}\langle (v(\tau)-v(0))^n\rangle\displaystyle\frac{\partial P}{\partial s} = \displaystyle\frac{\partial}{\partial v}(M_1P)+\displaystyle\frac{1}{2}\displaystyle\frac{\partial^2}{\partial v^2}(M_2P)M_1=\displaystyle\frac{1}{\tau}(v(\tau)-v(0)) = -\displaystyle\frac{1}{\tau}vM_2=\displaystyle\frac{1}{\tau}\langle (v(\tau)-v(0))^2\rangle = \displaystyle\frac{1}{\tau}\displaystyle\frac{6k_BT}{m}\tau\displaystyle\frac{\partial P}{\partial s} = \displaystyle\frac{\partial}{\partial v}(vP)+\displaystyle\frac{3k_BT}{m}\displaystyle\frac{\partial^2}{\partial v^2}PP(v,s|v_0)=\sqrt{\displaystyle\frac{m}{2\pi k_BT(1-e^{-2s/\tau})}}e^{-m(v-v_0e^{-s/\tau})^2/2k_BT(1-e^{-2s/\tau})}k_BmM_nnM_1sP_vsv0M_2tauTtt_0v_tauv0vv_0v_2xi

Symbol
Equation
Solved
Translated

Calculations

Symbol
Equation
Solved
Translated

 Variable   Given   Calculate   Target :   Equation   To be used
s=t-t_0\displaystyle\frac{\partial P}{\partial s}dv\tau=-\displaystyle\int_{-\infty}^{\infty}du P(v,s|v_0)P(u,\tau|v)+\displaystyle\int_{-\infty}^{\infty}du P(u,s|v_0)P(v,\tau|u)\displaystyle\int_{-\infty}^{\infty}du P(u,s|v)=1\displaystyle\frac{\partial P}{\partial s}\tau=- P(v,s|v_0)+\displaystyle\int_{-\infty}^{\infty}d\xi P(v-\xi,s|v_0)P(v,\tau|v-\xi)\displaystyle\frac{\partial}{\partial s}P(v,s|v_0)=\sum_{n=1}^{\infty}\displaystyle\frac{(-1)^n}{n!}\displaystyle\frac{\partial^n}{\partial v^n}[M_nP(v,s|v_0)]M_n=\displaystyle\frac{1}{\tau}\displaystyle\int_{-\infty}^{\infty}d\xi P(v+\xi,\tau|v)\xi^n=\displaystyle\frac{1}{\tau}\langle (v(\tau)-v(0))^n\rangle\displaystyle\frac{\partial P}{\partial s} = \displaystyle\frac{\partial}{\partial v}(M_1P)+\displaystyle\frac{1}{2}\displaystyle\frac{\partial^2}{\partial v^2}(M_2P)M_1=\displaystyle\frac{1}{\tau}(v(\tau)-v(0)) = -\displaystyle\frac{1}{\tau}vM_2=\displaystyle\frac{1}{\tau}\langle (v(\tau)-v(0))^2\rangle = \displaystyle\frac{1}{\tau}\displaystyle\frac{6k_BT}{m}\tau\displaystyle\frac{\partial P}{\partial s} = \displaystyle\frac{\partial}{\partial v}(vP)+\displaystyle\frac{3k_BT}{m}\displaystyle\frac{\partial^2}{\partial v^2}PP(v,s|v_0)=\sqrt{\displaystyle\frac{m}{2\pi k_BT(1-e^{-2s/\tau})}}e^{-m(v-v_0e^{-s/\tau})^2/2k_BT(1-e^{-2s/\tau})}k_BmM_nnM_1sP_vsv0M_2tauTtt_0v_tauv0vv_0v_2xi



Equations


Examples

Si se supone que una part cula tiene en un momento t_0 una velocidad v_0 existe una probabilidad\\n\\n

P(t,v|t_0,v_0)

que en un tiempo t posterior tenga una velocidad v.

Nota: Un proceso en que el nuevo estado depende del estado anterior se denomina un proceso de Markoff.

Para simplificar la formulaci n se considera el tiempo desde el tiempo inicial t_0 al que se le define como s con list

equation

Para comprender como cambia la probabilidad en el tiempo de que una part cula tenga una velocidad entre v y v+dv \\n\\n

\displaystyle\frac{\partial P}{\partial s}dv

\\n\\ndebemos considerar aquellas part culas que una vez llegaron a tener una velocidad en dicho rango y ahora presentan una velocidad v' distinta. Estas representan una reducci n (signo negativo) y deben ser sumadas sobre todas las velocidades v' posibles:\\n\\n

-\displaystyle\int_{-\infty}^{\infty}du P(v,s|v_0)P(u,\tau|v)

\\n\\nEn forma analog se sumaran part culas a aquellas en el rango de la velocidad si antes presentaban una distinta y son transferidas al nuevo rango:\\n\\n

+\displaystyle\int_{-\infty}^{\infty}du P(u,s|v_0)P(v,\tau|u)



Con ello la ecuaci n para la evoluci n temporal de la probabilidad es con list

equation

Como la probabilidad P(v,s|v_0) cubre todos los cambios posibles, una integraci n en v debe llevar a la unidad con list

equation

Con la ecuaci n diferencial de la probabilidad con list=9131

equation=9131



con la condici n de normalizaci n con list=9132

equation=9132



y el cambio de variable u=v-\xi se tiene la ecuaci n maestra de la probabilidad con list

equation

El producto P(v-\xi,s|v_0)P(v,\tau|v-\xi) se puede desarrollar en \xi lo que nos arroja\\n\\n

P(v-\xi,s)v_0)P(v,\tau|v-\xi)=\displaystyle\sum_{n=0}^{\infty}\displaystyle\frac{(-\xi)^n}{n!}\displaystyle\frac{\partial^n}{\partial v^n}[P(v,s)v_0)P(v+\xi,\tau|v)]



Con ello la ecuaci n con list=9133

equation=9133



se puede escribir con list como

equation



en donde el factor M_n se define con list=9135 como

equation=9135

El coeficiente M_n de la serie se define como\\n\\n

M_n=\displaystyle\frac{1}{\tau}\displaystyle\int_{-\infty}^{\infty}d\xi P(v+\xi,\tau|v)\xi^n



Estas expresiones corresponde a la promediaci n en el tiempo de las potencias de las diferencias de las velocidades v en el tiempo \tau con la inicial se puede escribir con list como

equation

Si se toma la ecuaci n maestra con la expansi n de Taylor con list=9134

equation=9134



hasta segundo orden se obtiene una ecuaci n para la distribuci n de probabilidad P(v,s|v_0) con list

equation

que se denomina la ecuaci n de Fokker-Planck

El factor M_1 que se obtiene con list=9134

equation=9134



corresponde a la velocidad con list

equation

El factor M_2 que se obtiene con list=9134 de

equation=9134



corresponde al factor de dispersi n estimado para la ecuaci n de Langevin con list

equation

Si se toma la ecuaci n de Fokker Planck con list=3136

equation=9136



y se empelan para los M_n los valores con list=3137

equation=9137



y con list=3138

equation=9138



se obtiene la ecuaci n con list

equation

Si se considera la ecuaci n de Fokker Planck con list=9139

equation=9139\\n\\ny se asume una soluci n del tipo\\n\\n

P(v,s)=e^{s/\tau}Q(u,s)

\\n\\nse obtiene la ecuaci n para Q(u,s) de la forma\\n\\n

\displaystyle\frac{\partial Q}{\partial s}=\displaystyle\frac{3k_BT}{m\tau}\displaystyle\frac{\partial^2Q}{\partial u^2}

\\n\\nque con el cambio de variable\\n\\n

\xi=\displaystyle\frac{1}{2}\tau(e^{2s/\tau}-1)

\\n\\nse obtiene la funci n de difusi n\\n\\n

\displaystyle\frac{\partial Q}{\partial s}=C\displaystyle\frac{\partial^2Q}{\partial\xi^2}

\\n\\ncon la constante\\n\\n

C=\displaystyle\frac{k_BT}{m\tau}

\\n\\nResolviendo esta ecuaci n se obtiene la soluci n\\n\\n

Q=\displaystyle\frac{1}{\sqrt{4\pi C\xi}}e^{-(u-u_0)^2/4C\xi}



con u_0 el valor inicial. Si se integran las distintas funciones se obtiene finalmente list

equation


>Model

ID:(1140, 0)