9.1: Derivación de la Ecuación de Difusión
- Page ID
- 116971
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)
( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\id}{\mathrm{id}}\)
\( \newcommand{\Span}{\mathrm{span}}\)
\( \newcommand{\kernel}{\mathrm{null}\,}\)
\( \newcommand{\range}{\mathrm{range}\,}\)
\( \newcommand{\RealPart}{\mathrm{Re}}\)
\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)
\( \newcommand{\Argument}{\mathrm{Arg}}\)
\( \newcommand{\norm}[1]{\| #1 \|}\)
\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)
\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)
\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)
\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)
\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vectorC}[1]{\textbf{#1}} \)
\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)
\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)
\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)
\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)
\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)
\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)Para derivar la ecuación de difusión en una dimensión espacial, imaginamos un quieto líquido en una tubería larga de área transversal constante. Se coloca una pequeña cantidad de tinte en una sección transversal de la tubería y se deja que se difumine hacia arriba y hacia abajo de la tubería. El colorante se difunde de regiones de mayor concentración a regiones de menor concentración.
Definimos\(u(x, t)\) que es la concentración del tinte en posición a\(x\) lo largo de la tubería, y deseamos encontrar el pde satisfecho por\(u\). Es útil realizar un seguimiento de las unidades de las diversas cantidades involucradas en la derivación e introducimos la notación de corchete\([X]\) para significar las unidades de\(X\). Las unidades dimensionales relevantes utilizadas en la derivación de la ecuación de difusión son masa\(m\)\(l\), longitud y tiempo\(t\). Suponiendo que la concentración de tinte es uniforme en cada sección transversal de la tubería, las dimensiones de concentración utilizadas aquí son\([u] = m/l\).
La masa de tinte en el volumen del tubo infinitesimal ubicado entre la posición\(x_1\) y la posición\(x_2\) en el tiempo\(t\), con\(x_1 < x < x_2\), se da al orden\(\Delta x = x_2 − x_1\) por\[M=u(x,t)\Delta x.\nonumber\]
La masa de tinte en este volumen de tubo infinitesimal cambia por difusión dentro y fuera de los extremos de sección transversal situados en posición\(x_1\) y\(x_2\) (Figura\(\PageIndex{1}\)). Suponemos que la velocidad de difusión es proporcional al gradiente de concentración, una relación conocida como ley de difusión de Fick. La ley de difusión de Fick asume el flujo másico\(J\), con unidades\([J] = m/t\) a través de una sección transversal de la tubería viene dada por \[\label{eq:1}J=-Du_x,\]donde la constante de difusión\(D > 0\) tiene unidades\([D] = l^2/t\), y hemos utilizado la notación\(u_x = \partial u/\partial x\). El flujo másico es opuesto en signo al gradiente de concentración de manera que el flujo es de alta concentración a baja concentración. La tasa de tiempo de cambio en la masa de tinte entre\(x_1\) y\(x_2\) viene dada por la diferencia entre el flujo másico hacia adentro y el flujo másico fuera de la sección transversal infinitesimal
volumen. Un flujo másico positivo significa difusión de izquierda a derecha. Por lo tanto, la tasa de tiempo de cambio de la masa de tinte viene dada por\[\frac{dM}{dt}=J(x_1, t)-J(x_2, t),\nonumber\] o reescritura en términos de\(u(x,t)\):
\[u_t(x,t)\Delta x=D(u_x(x_2,t)-u_x(x_1,t)).\nonumber\]
Dividir por\(\Delta x\) y tomar el límite\(\Delta x\to 0\) da como resultado la ecuación de difusión:
\[u_t=Du_{xx}.\nonumber\]
Observamos que la ecuación de difusión es idéntica a la ecuación de conducción de calor, donde\(u\) está la temperatura, y la constante\(D\) (comúnmente escrita como\(\kappa\)) es la conductividad térmica.