Divergencia

Operador de divergencia como tasa de flujo neto saliente por unidad de volumen en campos vectoriales.

Que es

La divergencia es un operador diferencial que actúa sobre un campo vectorial diferenciable $\vec{F}(x,y,z)$ para producir un campo escalar $\vec{\nabla}\cdot\vec{F}$. En coordenadas cartesianas se define como:

$\vec{\nabla}\cdot\vec{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}$

Físicamente, la divergencia mide la densidad volumétrica de flujo neto saliente de un campo vectorial en cada punto del espacio. Se define formalmente mediante el límite del flujo a través de una superficie cerrada $S$ que encierra un volumen $V$, cuando dicho volumen se contrae hacia cero:

$\vec{\nabla}\cdot\vec{F} = \lim_{V \to 0} \frac{1}{V} \oint_S \vec{F} \cdot d\vec{S}$

• Si $\vec{\nabla}\cdot\vec{F}(P) > 0$: El punto $P$ actúa como una fuente neta (el fluido o las líneas de campo se originan o expanden desde allí).
• Si $\vec{\nabla}\cdot\vec{F}(P) < 0$: El punto $P$ actúa como un sumidero neto (el fluido o las líneas de campo convergen y desaparecen allí).
• Si $\vec{\nabla}\cdot\vec{F}(P) = 0$: El campo es solenoidal o incompresible en $P$ (todo lo que entra en un volumen elemental infinitesimal sale exactamente en la misma cantidad).

Por qué aparece en física

La divergencia es el lenguaje natural para expresar las leyes fundamentales de conservación y las ecuaciones de Maxwell:

• Ley de Gauss electrostática: $\vec{\nabla}\cdot\vec{E} = \frac{\rho}{\varepsilon_0}$. Las cargas eléctricas positivas son fuentes del campo eléctrico; las negativas son sumideros.
• Ley de Gauss magnética: $\vec{\nabla}\cdot\vec{B} = 0$. No existen monopolos magnéticos aislados en la naturaleza; las líneas de campo magnético nunca nacen ni mueren, son siempre cerradas.
• Ecuación de continuidad de fluidos: $\frac{\partial \rho}{\partial t} + \vec{\nabla}\cdot(\rho \vec{v}) = 0$. La divergencia del flujo de masa cuantifica la tasa local de acumulación o pérdida de masa fluida.
• Fluidos incompresibles: $\vec{\nabla}\cdot\vec{v} = 0$. El agua y los líquidos a velocidades moderadas conservan su volumen local durante el movimiento.

Idea visual o geométrica

Imagina una esponja sumergida en agua:

• Si abres una llave sumergida que inyecta agua nueva al espacio, los vectores velocidad del agua apuntan hacia afuera en todas direcciones: la divergencia en la boca de la llave es positiva (+).
• Si abres un desagüe que succiona agua, los vectores convergen hacia el agujero: la divergencia en el desagüe es negativa (-).
• Si el agua simplemente fluye de manera uniforme por una tubería sin crearse ni destruirse masa líquida, la divergencia es exactamente cero.

Lenguaje simbólico mínimo

• Operador en cartesianas: $\vec{\nabla}\cdot\vec{F} = \frac{\partial F_x}{\partial x} + \frac{\partial F_y}{\partial y} + \frac{\partial F_z}{\partial z}$.
• En cilíndricas $(r, \theta, z)$: $\vec{\nabla}\cdot\vec{F} = \frac{1}{r} \frac{\partial (r F_r)}{\partial r} + \frac{1}{r} \frac{\partial F_\theta}{\partial \theta} + \frac{\partial F_z}{\partial z}$.
• En esféricas $(r, \theta, \varphi)$: $\vec{\nabla}\cdot\vec{F} = \frac{1}{r^2} \frac{\partial (r^2 F_r)}{\partial r} + \frac{1}{r \sin\theta} \frac{\partial (\sin\theta F_\theta)}{\partial \theta} + \frac{1}{r \sin\theta} \frac{\partial F_\varphi}{\partial \varphi}$.
• Teorema de Gauss (o de la divergencia): $\int_V (\vec{\nabla}\cdot\vec{F}) \, dV = \oint_S \vec{F} \cdot d\vec{S}$.
• Campo solenoidal: $\vec{\nabla}\cdot\vec{F} = 0 \iff \vec{F} = \vec{\nabla}\times\vec{A}$ (existe un potencial vector $\vec{A}$).

Que significa el resultado

El escalar resultante $\vec{\nabla}\cdot\vec{F}$ no posee dirección: representa la tasa neta por unidad de volumen a la que el campo "nace" o "muere" en ese punto. Sus unidades físicas son siempre las del campo vectorial divididas por metros ($[\vec{F}]/[\text{m}]$).

Cuando no debes usarlo

• No intentes calcular la divergencia de un campo escalar $\phi$. La divergencia exige un vector con componentes direccionales.
• No omitas la derivada del producto métrico en coordenadas no cartesianas: por ejemplo, en esféricas es un error mortal calcular $\partial F_r / \partial r$ en lugar de $\frac{1}{r^2} \frac{\partial (r^2 F_r)}{\partial r}$.
• No asumas que porque un campo sea radial su divergencia es no nula: el campo gravitatorio newtoniano $\vec{g} = -\frac{GM}{r^2} \hat{e}_r$ tiene divergencia idénticamente cero en todo el espacio fuera del origen ($\vec{\nabla}\cdot\vec{g} = 0$ para $r > 0$).

Resumen operativo

Para calcular $\vec{\nabla}\cdot\vec{F}$: diferencia cada componente espacial respecto a su coordenada correspondiente (incluyendo los factores de escala métrica en cilíndricas o esféricas) y suma algebraicamente los términos. El resultado es un escalar local que delata fuentes positivas o sumideros negativos de flujo.

Procedimiento

Antes de calcular

Confirma que el campo a evaluar es estrictamente vectorial: $\vec{F} = (F_x, F_y, F_z)$ o $\vec{F} = F_r \hat{e}_r + F_\theta \hat{e}_\theta + F_z \hat{e}_z$. Identifica la geometría del dominio.

Paso 1: identifica el objeto matemático

Desglosa las componentes del campo vectorial en el sistema de coordenadas adoptado. Asegúrate de verificar las dependencias funcionales de cada componente.

Paso 2: elige la operación

Selecciona la fórmula correcta del operador divergencia:

• Cartesianas: $\partial_x F_x + \partial_y F_y + \partial_z F_z$.
• Cilíndricas: $\frac{1}{r} \frac{\partial (r F_r)}{\partial r} + \frac{1}{r} \frac{\partial F_\theta}{\partial \theta} + \frac{\partial F_z}{\partial z}$.
• Esféricas: $\frac{1}{r^2} \frac{\partial (r^2 F_r)}{\partial r} + \frac{1}{r \sin\theta} \frac{\partial (\sin\theta F_\theta)}{\partial \theta} + \frac{1}{r \sin\theta} \frac{\partial F_\varphi}{\partial \varphi}$.

Paso 3: conserva unidades y signos

Diferencia cada término algebraicamente. En coordenadas no cartesianas, multiplica primero la componente por su factor métrico ($r$ o $r^2$) antes de derivar respecto a la coordenada radial.

Paso 4: calcula

Suma todos los términos escalares resultantes. Si se evalúa en un punto numérico $P(x_0, y_0, z_0)$, sustituye las coordenadas para obtener el valor escalar final.

Paso 5: interpreta el resultado

Analiza el signo del escalar obtenido:

• Positivo ($> 0$): densidad de fuentes (expansión neta).
• Negativo ($< 0$): densidad de sumideros (compresión o absorción neta).
• Cero ($= 0$): flujo incompresible o libre de fuentes.

Checklist de autocontrol

• ¿El resultado obtenido es un número escalar (sin vectores unitarios $\hat{i}, \hat{j}, \hat{k}$)?
• ¿Las unidades del resultado equivalen a $[\vec{F}]/[\text{m}]$?
• ¿Se incluyó el término $r$ o $r^2$ dentro de la derivada radial en cilíndricas/esféricas?

Ejemplo guiado

Situación física

Un campo electrostático no uniforme en una región del espacio viene dado por: $\vec{E}(x,y,z) = (2x^2 y) \hat{i} - (y^2 z) \hat{j} + (3z^2 x) \hat{k} \quad [\text{N/C}]$ donde $x,y,z$ están en metros. Se desea calcular la densidad volumétrica de carga eléctrica libre $\rho(x,y,z)$ en el punto $P(1, 2, -1)$ mediante la ley de Gauss diferencial: $\vec{\nabla}\cdot\vec{E} = \frac{\rho}{\varepsilon_0}$ donde $\varepsilon_0 \approx 8.854 \times 10^{-12} \text{ F/m}$.

Traducción matemática

Calculamos la divergencia cartesiana: $\vec{\nabla}\cdot\vec{E} = \frac{\partial E_x}{\partial x} + \frac{\partial E_y}{\partial y} + \frac{\partial E_z}{\partial z}$ y despejamos la densidad de carga como $\rho = \varepsilon_0 (\vec{\nabla}\cdot\vec{E})$.

Calculo paso a paso

1. Derivada parcial de $E_x = 2x^2 y$ respecto a $x$: $\frac{\partial E_x}{\partial x} = \frac{\partial}{\partial x}(2x^2 y) = 4xy$

2. Derivada parcial de $E_y = -y^2 z$ respecto a $y$: $\frac{\partial E_y}{\partial y} = \frac{\partial}{\partial y}(-y^2 z) = -2yz$

3. Derivada parcial de $E_z = 3z^2 x$ respecto a $z$: $\frac{\partial E_z}{\partial z} = \frac{\partial}{\partial z}(3z^2 x) = 6zx$

4. Suma de la divergencia: $\vec{\nabla}\cdot\vec{E} = 4xy - 2yz + 6zx \quad [\text{N/(C}\cdot\text{m)}]$

5. Evaluación en el punto $P(1, 2, -1)$: $\vec{\nabla}\cdot\vec{E}(P) = 4(1)(2) - 2(2)(-1) + 6(-1)(1) = 8 - (-4) - 6 = 8 + 4 - 6 = +6 \text{ N/(C}\cdot\text{m)}$

6. Determinación de la densidad de carga $\rho$: $\rho = \varepsilon_0 (\vec{\nabla}\cdot\vec{E}) = (8.854 \times 10^{-12} \text{ C}^2/(\text{N}\cdot\text{m}^2)) \times (6 \text{ N}/(\text{C}\cdot\text{m}))$ $\rho = 5.3124 \times 10^{-11} \text{ C/m}^3 = +53.12 \text{ pC/m}^3$

Lectura del resultado

La divergencia es positiva ($+6$), indicando que en el punto $P$ las líneas de campo eléctrico experimentan una dispersión neta hacia afuera. Por la ley de Gauss, esto requiere físicamente una densidad de carga positiva neta de $+53.12 \text{ pC/m}^3$.

Error común evitado

Sumar las derivadas conservando vectores unitarios (como $4xy \hat{i} - 2yz \hat{j} + 6zx \hat{k}$). La divergencia es un producto escalar operacional que produce estrictamente un campo escalar, no un vector.

Ejemplo inverso

Resultado dado

El campo de velocidades de un flujo incompresible en dos dimensiones tiene componente horizontal: $v_x(x,y) = 3x^2 - 2y$

Que pregunta física permite responder

¿Cuál debe ser la forma matemática de la componente vertical $v_y(x,y)$ para garantizar la conservación de la masa líquida? La condición de incompresibilidad impone: $\vec{\nabla}\cdot\vec{v} = \frac{\partial v_x}{\partial x} + \frac{\partial v_y}{\partial y} = 0$ Calculamos la derivada de $v_x$: $\frac{\partial v_x}{\partial x} = \frac{\partial}{\partial x}(3x^2 - 2y) = 6x$ Por tanto: $\frac{\partial v_y}{\partial y} = -6x$ Integrando respecto a $y$: $v_y(x,y) = -6xy + g(x)$ donde $g(x)$ es una función arbitraria que depende exclusivamente de $x$. Cualquier campo con esta estructura vertical garantiza divergencia idénticamente nula y conservación perfecta de volumen.

Errores frecuentes

Error 1: Devolver un vector en lugar de un escalar

Escribir la divergencia manteniendo los vectores unitarios cartesianos $\hat{i}, \hat{j}, \hat{k}$. La divergencia es un producto escalar diferencial y da como resultado un número escalar en cada punto.

Error 2: Derivar la componente radial sin el factor $r^2$ en esféricas

Calcular $\partial_r F_r$ en lugar de $\frac{1}{r^2} \partial_r (r^2 F_r)$. Esto llevaría a la conclusión errónea de que el campo de Coulomb tiene divergencia no nula en el vacío.

Error 3: Confundir un campo sin divergencia con un campo constante

Creer que $\vec{\nabla}\cdot\vec{F} = 0$ implica que $\vec{F}$ es constante en el espacio. Campos altamente complejos como el campo magnético de un toroide o la velocidad de un vórtice tienen divergencia nula y son fuertemente variables.

Error 4: Olvidar la singularidad de Dirac en el origen

Para un campo $\vec{F} = \frac{\hat{r}}{r^2}$, afirmar que su divergencia es cero en todo el espacio sin excluir $r=0$. En el origen, $\vec{\nabla}\cdot\left(\frac{\hat{r}}{r^2}\right) = 4\pi \delta^3(\vec{r})$.

Error 5: Aplicar la divergencia a un campo escalar

Intentar calcular $\vec{\nabla}\cdot\phi$. La divergencia solo actúa sobre campos vectoriales; sobre escalares, la operación correspondiente es el gradiente $\vec{\nabla}\phi$.

Regla rápida de detección

Si el resultado de una divergencia contiene componentes direccionales, o si en coordenadas polares se omitió el término métrico $1/r$, el planteamiento contiene un error de formalismo.

Aplicaciones físicas

1. Electromagnetismo y leyes de Maxwell

Las dos primeras ecuaciones de Maxwell gobiernan la divergencia:

• $\vec{\nabla}\cdot\vec{E} = \rho/\varepsilon_0$ (relaciona el campo eléctrico con la densidad de carga fuente).
• $\vec{\nabla}\cdot\vec{B} = 0$ (establece la naturaleza solenoidal del campo magnético y la ausencia de monopolos).

2. Hidrodinámica e incompresibilidad de fluidos

La ecuación de Navier-Stokes para fluidos incompresibles (agua, sangre, petróleo en tuberías) se simplifica enormemente mediante la condición solenoidal $\vec{\nabla}\cdot\vec{v} = 0$.

3. Dinámica de gases y ondas de choque acústicas

En flujo compresible (aerodinámica supersónica), la divergencia de la velocidad $\vec{\nabla}\cdot\vec{v}$ cuantifica la tasa de compresión o dilatación del gas y predice la aparición de ondas de choque.

4. Transferencia de calor y generación interna

En termodinámica, el balance de energía en un sólido con fuentes internas de calor (reacciones exotérmicas, decaimiento radiactivo) se modela mediante $-\vec{\nabla}\cdot\vec{q} + Q = \rho c_p \frac{\partial T}{\partial t}$.

5. Astrofísica y equilibrio estelar

En el interior de una estrella, el flujo radiativo satisface $\vec{\nabla}\cdot\vec{F}_{\text{rad}} = \rho \epsilon$ (donde $\epsilon$ es la tasa de generación nuclear de energía en el núcleo estelar).

Practica

Problema 1: Divergencia de un campo de fuerza radial en coordenadas esféricas

Un campo de fuerzas newtoniano varía como $\vec{F}(\vec{r}) = r^n \hat{e}_r$, donde $r = \|\vec{r}\|$ y $n$ es una constante entera.

1. Utilizando la expresión de la divergencia en coordenadas esféricas, halle $\vec{\nabla}\cdot\vec{F}$ para todo $r > 0$.
2. ¿Para qué valor particular del exponente $n$ el campo tiene divergencia idénticamente cero (campo solenoidal)?

Solución:

1. En coordenadas esféricas con dependencia puramente radial ($F_\theta = F_\varphi = 0$): $\vec{\nabla}\cdot\vec{F} = \frac{1}{r^2} \frac{\partial}{\partial r}(r^2 F_r) = \frac{1}{r^2} \frac{\partial}{\partial r}(r^2 r^n) = \frac{1}{r^2} \frac{\partial}{\partial r}(r^{n+2})$ $\vec{\nabla}\cdot\vec{F} = \frac{1}{r^2} (n+2) r^{n+1} = (n+2) r^{n-1}$
2. Para que $\vec{\nabla}\cdot\vec{F} = 0$, se requiere: $n + 2 = 0 \implies n = -2$ Esto demuestra rigurosamente por qué las leyes de la gravitación universal y de Coulomb deben seguir una ley del inverso del cuadrado $\sim 1/r^2$: es el único exponente para el cual el campo es solenoidal en todo el espacio libre de fuentes.

Problema 2: Compresibilidad en dinámica atmosférica

El viento en una capa atmosférica se modela mediante el campo de velocidades: $\vec{v}(x,y) = (U_0 \cos(k x) \cosh(k y)) \hat{i} - (U_0 \sin(k x) \sinh(k y)) \hat{j}$

1. Calcule la divergencia $\vec{\nabla}\cdot\vec{v}$.
2. Determine si el aire experimenta convergencia barométrica o si se comporta de forma incompresible.

Solución:

1. Calculamos derivadas: $\frac{\partial v_x}{\partial x} = \frac{\partial}{\partial x}[U_0 \cos(k x) \cosh(k y)] = -k U_0 \sin(k x) \cosh(k y)$ $\frac{\partial v_y}{\partial y} = \frac{\partial}{\partial y}[-U_0 \sin(k x) \sinh(k y)] = -k U_0 \sin(k x) \cosh(k y)$
2. Sumando las componentes: $\vec{\nabla}\cdot\vec{v} = -2 k U_0 \sin(k x) \cosh(k y)$ Como la divergencia no es nula, el flujo atmosférico no es incompresible. Las zonas donde $\vec{\nabla}\cdot\vec{v} < 0$ corresponden a regiones de convergencia de masa donde se acumula humedad y se forman borrascas meteorológicas.

Problema 3: Corriente de convección y ley de conservación de carga

La densidad de corriente en un conductor es $\vec{J}(x,y,z,t) = J_0 e^{-\alpha t} (x \hat{i} + y \hat{j} + z \hat{k})$. Utilizando la ecuación de continuidad $\frac{\partial \rho}{\partial t} + \vec{\nabla}\cdot\vec{J} = 0$, determine la tasa de variación de la densidad de carga $\frac{\partial \rho}{\partial t}$.

Solución:

1. Divergencia de $\vec{J}$: $\vec{\nabla}\cdot\vec{J} = J_0 e^{-\alpha t} \left( \frac{\partial x}{\partial x} + \frac{\partial y}{\partial y} + \frac{\partial z}{\partial z} \right) = J_0 e^{-\alpha t} (1 + 1 + 1) = 3 J_0 e^{-\alpha t}$
2. Por la ecuación de continuidad: $\frac{\partial \rho}{\partial t} = -\vec{\nabla}\cdot\vec{J} = -3 J_0 e^{-\alpha t}$ La carga neta disminuye exponencialmente a una tasa proporcional a $3 J_0$.