# Energía, fuerzas y métodos de solución en sistemas capacitivos: del campo eléctrico a los algoritmos numéricos

El proceso de carga de un condensador constituye uno de los ejemplos más reveladores sobre la conservación de la energía y la naturaleza del campo eléctrico. Cuando se transfieren portadores de carga de una armadura a la otra, se debe realizar un trabajo para superar la diferencia de potencial que se genera acumulativamente entre ambas. Este trabajo no se pierde, sino que se almacena en el sistema en forma de energía potencial electrostática. La derivación clásica, partiendo de la relación entre carga instantánea, capacidad y diferencia de potencial, permite obtener las expresiones fundamentales que gobiernan el almacenamiento energético en estos dispositivos (Griffiths, 2017, cap. 2; Purcell y Morin, 2013, cap. 3).

## El trabajo de carga y la energía almacenada

Consideremos un condensador de capacidad \(C\) inicialmente descargado. Para transferir un elemento infinitesimal de carga \(dq\) desde la armadura negativa a la positiva, cuando existe una diferencia de potencial instantánea \(v = q/C\) entre ellas, se requiere un trabajo \(dW = v\,dq = (q/C)\,dq\). El trabajo total para llevar la carga desde \(0\) hasta \(Q_f\) es entonces

\[
W = \int_0^{Q_f} \frac{q}{C}\,dq = \frac{Q_f^2}{2C}.
\]

Este trabajo queda almacenado como energía potencial electrostática \(U\) en el sistema. Utilizando la relación \(Q_f = C \varphi_{12}\), donde \(\varphi_{12}\) es la diferencia de potencial final entre las armaduras, obtenemos las expresiones equivalentes

\[
U = \frac{Q_f^2}{2C} = \frac{1}{2} C \varphi_{12}^2 = \frac{1}{2} Q_f \varphi_{12},
\]

que constituyen las ecuaciones fundamentales de la energía almacenada en un condensador (Jackson, 1999, sec. 2.5; Hyperphysics, 2025).

Un aspecto particularmente fascinante de esta formulación es que la energía no reside en las placas del condensador como tales, sino en el campo eléctrico que existe en el espacio entre ellas (Purcell y Morin, 2013, cap. 3). Para el condensador de placas paralelas con área \(A\) y separación \(s\), la capacidad es \(C = \varepsilon_0 A / s\) y la diferencia de potencial es \(\varphi_{12} = E s\), donde \(E\) es la magnitud del campo eléctrico uniforme entre las placas. Sustituyendo en \(U = C\varphi_{12}^2/2\):

\[
U = \frac{1}{2} \left( \frac{\varepsilon_0 A}{s} \right) (E s)^2 = \frac{1}{2} \varepsilon_0 E^2 (A s).
\]

El producto \(A s\) es precisamente el volumen \(V\) donde existe el campo eléctrico uniforme. Se deduce entonces que la **densidad de energía electrostática** en el vacío es

\[
u = \frac{U}{V} = \frac{1}{2} \varepsilon_0 E^2.
\]

Este resultado, confirmado por la teoría y la experimentación, revoluciona la concepción de que la energía reside únicamente en las cargas: el campo eléctrico es un ente físico real capaz de almacenar energía localmente, incluso en el vacío (OpenStax, 2021; Smythe, 1968, cap. 2). La densidad de energía es proporcional al cuadrado del campo, lo que implica que duplicar el campo cuadruplica la energía almacenada por unidad de volumen.

## El método de los trabajos virtuales: fuerzas electrostáticas desde la energía

El método de los trabajos virtuales ofrece una herramienta poderosa para determinar fuerzas mecánicas en sistemas electrostáticos sin recurrir al cálculo directo de interacciones coulombianas, que suele ser sumamente complejo debido a la distribución tridimensional de la carga (Landau y Lifshitz, 1984, cap. 2). La idea central es que la fuerza mecánica asociada a una coordenada generalizada \(x\) puede obtenerse a partir de la variación de la energía electrostática mediante

\[
F_x = -\left( \frac{\partial U}{\partial x} \right)_{\text{variables eléctricas fijas}},
\]

donde las variables que se mantienen constantes durante la derivada dependen de las condiciones de contorno del problema.

### Fuerza a carga constante

Para un condensador de placas paralelas con carga constante \(Q\) y separación variable \(x\), la energía almacenada es

\[
U(x) = \frac{Q^2}{2C(x)} = \frac{Q^2 x}{2\varepsilon_0 A},
\]

pues \(C(x) = \varepsilon_0 A / x\). La fuerza atractiva entre las placas se obtiene derivando:

\[
F = -\frac{dU}{dx} = -\frac{Q^2}{2\varepsilon_0 A},
\]

donde el signo negativo indica que la fuerza es atractiva (tiende a disminuir la separación). En magnitud,

\[
|F| = \frac{Q^2}{2\varepsilon_0 A}.
\]

Expresando esta fuerza en términos del campo eléctrico neto entre las placas, dado que \(E = \sigma / \varepsilon_0 = Q / (\varepsilon_0 A)\), obtenemos

\[
|F| = \frac{1}{2} Q E,
\]

donde el factor \(1/2\) tiene una interpretación física profunda: la fuerza sobre una placa no se debe al campo eléctrico neto total \(E\), sino al campo generado exclusivamente por la otra placa, que tiene una magnitud exactamente igual a la mitad del campo total (Purcell y Morin, 2013, cap. 3). En efecto, para dos placas paralelas infinitas con densidades de carga \(+\sigma\) y \(-\sigma\), cada placa genera un campo \(\sigma/(2\varepsilon_0)\); el campo total entre ellas es la suma \(\sigma/\varepsilon_0\), pero la fuerza sobre una placa proviene solo del campo creado por la otra.

### El caso a potencial constante: matices energéticos

Un análisis detallado del problema a potencial constante revela interesantes matices sobre el balance energético cuando el condensador permanece conectado a una batería (OCW-UC3M, 2025). En este caso, al variar la separación entre placas manteniendo \(V = V_0\) constante, la batería debe retirar carga del condensador para mantener el potencial fijo, realizando un trabajo eléctrico sobre el sistema. La conservación de la energía establece que

\[
dW_{\text{mec}} + dW_{\text{batería}} = dU_{\text{condensador}}.
\]

Para un desplazamiento \(dx\) que aumenta la separación, la capacidad disminuye (\(dC < 0\)), y como \(Q = CV_0\), la carga también disminuye (\(dQ < 0\)). La batería absorbe energía: \(dW_{\text{batería}} = V_0 dQ\). El cambio en la energía del condensador es \[ dU = \frac{1}{2} V_0^2 dC = \frac{1}{2} V_0^2 \left( -\frac{\varepsilon_0 A}{x^2} dx \right) = -\frac{1}{2} \frac{\varepsilon_0 A V_0^2}{x^2} dx. \] El trabajo mecánico realizado por el agente externo es \(dW_{\text{mec}} = F\,dx\). Sustituyendo en la ecuación de conservación: \[ F\,dx + V_0 dQ = dU. \] Dado que \(dQ = V_0 dC = -(\varepsilon_0 A V_0 / x^2) dx\), obtenemos \[ F\,dx - \frac{\varepsilon_0 A V_0^2}{x^2} dx = -\frac{1}{2} \frac{\varepsilon_0 A V_0^2}{x^2} dx, \] de donde \[ F = \frac{\varepsilon_0 A V_0^2}{2 x^2} = \frac{Q E}{2}, \] que es idéntica en magnitud a la fuerza calculada a carga constante (Smythe, 1968, cap. 2). Sorprendentemente, aunque el papel de la batería altera significativamente el balance de energía —absorbiendo el doble de la energía que el condensador pierde—, el valor absoluto de la fuerza de atracción entre las placas resulta ser el mismo. Esto ilustra que la fuerza electrostática es una propiedad del estado del sistema, no del proceso termodinámico seguido para alcanzarlo. ## Solución de la ecuación de Laplace: métodos analíticos y numéricos Cuando los sistemas electrostáticos presentan fronteras sin simetrías perfectas, las soluciones analíticas directas de la ecuación de Laplace \(\nabla^2 \varphi = 0\) no son viables. Dos alternativas metodológicas fundamentales son la **representación conforme** mediante funciones de variable compleja y el **método numérico de relajación**. ### Representación conforme En problemas bidimensionales, las funciones analíticas de variable compleja ofrecen una herramienta elegante. Si \(w = f(z) = u(x,y) + i v(x,y)\) es una función analítica, entonces \(u\) y \(v\) satisfacen las ecuaciones de Cauchy-Riemann y, en consecuencia, ambas cumplen la ecuación de Laplace bidimensional: \[ \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0, \qquad \frac{\partial^2 v}{\partial x^2} + \frac{\partial^2 v}{\partial y^2} = 0. \] Además, las curvas de nivel de \(u\) y \(v\) son ortogonales entre sí, lo que permite identificar una de ellas con las líneas equipotenciales y la otra con las líneas de campo (Smythe, 1968, cap. 4; Jackson, 1999, cap. 2). La transformación conforme \(w = f(z)\) mapea un dominio complicado en el plano \(z\) a un dominio más simple en el plano \(w\) donde la ecuación de Laplace es fácil de resolver. Este método es particularmente útil para geometrías con bordes afilados, esquinas y discontinuidades, como en el caso de condensadores con placas no paralelas o electrodos de geometría compleja. ### Método de relajación de Gauss-Seidel El método de relajación se fundamenta en la discretización de la ecuación de Laplace sobre una malla de puntos. Consideremos una región bidimensional sin cargas libres, donde el potencial satisface \(\nabla^2 \varphi = 0\). Discretizando el espacio con un espaciado uniforme \(h\) en ambas direcciones, el potencial en un punto \((x_i, y_j) = (ih, jh)\) se denota como \(\varphi_{i,j}\). Mediante expansiones en serie de Taylor del potencial alrededor de \((i,j)\): \[ \varphi_{i+1,j} = \varphi_{i,j} + h \frac{\partial \varphi}{\partial x} + \frac{h^2}{2} \frac{\partial^2 \varphi}{\partial x^2} + \mathcal{O}(h^3), \] \[ \varphi_{i-1,j} = \varphi_{i,j} - h \frac{\partial \varphi}{\partial x} + \frac{h^2}{2} \frac{\partial^2 \varphi}{\partial x^2} + \mathcal{O}(h^3). Sumando ambas expresiones y despejando la segunda derivada: \[ \frac{\partial^2 \varphi}{\partial x^2} \approx \frac{\varphi_{i+1,j} + \varphi_{i-1,j} - 2\varphi_{i,j}}{h^2}. \] Procediendo análogamente para \(y\), la ecuación de Laplace discretizada resulta: \[ \frac{\varphi_{i+1,j} + \varphi_{i-1,j} - 2\varphi_{i,j}}{h^2} + \frac{\varphi_{i,j+1} + \varphi_{i,j-1} - 2\varphi_{i,j}}{h^2} = 0, \] de donde se obtiene la condición fundamental: \[ \varphi_{i,j} = \frac{1}{4} \left( \varphi_{i+1,j} + \varphi_{i-1,j} + \varphi_{i,j+1} + \varphi_{i,j-1} \right). \] Es decir, el potencial en cualquier nodo interno libre de carga es aproximadamente igual al promedio aritmético del potencial de sus cuatro nodos vecinos adyacentes (Griffiths, 2017, cap. 3; Smythe, 1968, cap. 3). El algoritmo de relajación (Gauss-Seidel) consiste en iterar repetidamente sobre la malla, actualizando cada nodo con el promedio de sus vecinos, hasta alcanzar la convergencia deseada. Para ilustrarlo, consideremos un problema bidimensional con cuatro nodos internos \((a, b, c, d)\) en una malla \(2 \times 2\), rodeados por fronteras con potenciales fijos: 100 V en el borde superior e izquierdo, y 0 V en el borde inferior y derecho. Inicialmente, asignamos \(\varphi = 0\) a todos los nodos internos. En la primera iteración, aplicando el método de Gauss-Seidel: - Nodo \(a\) (frontera izquierda a 100 V, superior a 100 V, vecinos derecho e inferior a 0 V): \[ \varphi_a = \frac{100 + 100 + 0 + 0}{4} = 50\,\mathrm{V}. \] - Nodo \(b\) (vecino izquierdo \(a\) recién actualizado a 50 V, superior a 100 V, derecho a 0 V, inferior a 0 V): \[ \varphi_b = \frac{50 + 100 + 0 + 0}{4} = 37.5\,\mathrm{V}. \] - Nodo \(c\) (vecino izquierdo a 100 V, superior \(a\) a 50 V, derecho a 0 V, inferior a 0 V): \[ \varphi_c = \frac{100 + 50 + 0 + 0}{4} = 37.5\,\mathrm{V}. \] - Nodo \(d\) (vecino izquierdo \(c\) a 37.5 V, superior \(b\) a 37.5 V, derecho a 0 V, inferior a 0 V): \[ \varphi_d = \frac{37.5 + 37.5 + 0 + 0}{4} = 18.75\,\mathrm{V}. \] Sin embargo, notemos que el texto original del usuario daba valores de 50 V, 62.5 V, 37.5 V y 25 V para los nodos \(a, b, c, d\), lo que corresponde a condiciones de contorno ligeramente diferentes (posiblemente con el nodo \(b\) teniendo un vecino superior a 100 V y vecino derecho a 100 V). En cualquier caso, el principio es el mismo: este proceso iterativo debe repetirse sucesivamente en un código computacional hasta que la variación de potencial en todos los nodos entre dos ciclos consecutivos sea despreciable (inferior a un criterio de tolerancia predefinido, típicamente \(10^{-6}\) o \(10^{-8}\)). La convergencia puede acelerarse mediante la técnica de **sobrerrelajación sucesiva** (SOR), donde se introduce un parámetro \(\omega > 1\) para sobreamplificar la corrección en cada iteración (Jackson, 1999, cap. 2).

## Implicaciones conceptuales y aplicaciones

La densidad de energía del campo eléctrico, \(u = \varepsilon_0 E^2 / 2\), tiene profundas implicaciones conceptuales. Establece que el campo eléctrico no es meramente un artificio matemático para calcular fuerzas entre cargas, sino una entidad física que posee energía, momento y, como muestra la teoría de la relatividad, masa equivalente (\(m = U/c^2\)). En un condensador de placas paralelas con \(A = 0.01\,\mathrm{m}^2\), \(s = 1\,\mathrm{mm}\) y \(V = 100\,\mathrm{V}\), el campo es \(E = 10^5\,\mathrm{V/m}\), la densidad de energía \(u \approx 0.044\,\mathrm{J/m^3}\), y la energía total almacenada es \(U \approx 4.4 \times 10^{-7}\,\mathrm{J}\). Aunque modesta, esta energía es suficiente para aplicaciones en circuitos de temporización y filtrado.

El método de los trabajos virtuales, por su parte, encuentra aplicación directa en el diseño de actuadores electrostáticos, micromáquinas (MEMS) y sistemas de posicionamiento de precisión, donde la fuerza entre conductores cargados debe calcularse con exactitud sin recurrir a complejas integrales de campo (Landau y Lifshitz, 1984, cap. 2). La distinción entre condiciones de carga constante y potencial constante es crucial en el diseño experimental: un electrómetro de cuadrante, por ejemplo, opera a potencial constante, mientras que un generador de Van de Graaff funciona a carga constante.

Finalmente, los métodos numéricos como el de relajación constituyen la base de los modernos simuladores de campos electromagnéticos utilizados en la industria electrónica y de semiconductores. La combinación de la elegancia matemática de la representación conforme con la potencia computacional de los métodos iterativos permite abordar problemas de ingeniería que van desde el diseño de condensadores de precisión hasta el modelado de líneas de transmisión en circuitos integrados de alta frecuencia (Harrington, 1968, cap. 3). La ecuación de Laplace, simple en apariencia, revela así una riqueza de soluciones y métodos que conectan la física teórica con la práctica computacional.

—

## Referencias

– Griffiths, D. J. (2017). *Introduction to Electrodynamics* (4.ª ed.). Cambridge University Press.
– Harrington, R. F. (1968). *Field Computation by Moment Methods*. Macmillan.
– Hyperphysics (2025). *Energía de Campo Eléctrico en Condensador*. Georgia State University. http://hyperphysics.phy-astr.gsu.edu/hbasees/electric/capeng.html
– Jackson, J. D. (1999). *Classical Electrodynamics* (3.ª ed.). Wiley.
– Landau, L. D. y Lifshitz, E. M. (1984). *Electrodynamics of Continuous Media* (2.ª ed.). Pergamon Press.
– OCW-UC3M (2025). *Energía del campo electrostático: Condensadores*. Universidad Carlos III de Madrid. https://ocw.uc3m.es/pluginfile.php/1353/mod_page/content/19/condensadores_tema5.pdf
– OpenStax (2021). *8.3 Energía almacenada en un condensador*. Física Universitaria Volumen 2. https://openstax.org/books/f%C3%ADsica-universitaria-volumen-2/pages/8-3-energia-almacenada-en-un-condensador
– Purcell, E. M. y Morin, D. J. (2013). *Electricity and Magnetism* (3.ª ed.). Cambridge University Press.
– Smythe, W. R. (1968). *Static and Dynamic Electricity* (3.ª ed.). McGraw-Hill.