Campo eléctrico

Vector campo eléctrico

En electrostática se estudia la fuerza que se ejercen mutuamente dos cargas eléctricas puntuales Q, Q' entre ellas; si fijamos sólo en una de las cargas, digamos la Q, nos damos cuenta de que si colocamos cualquier carga Q' en cualquier posición del espacio que rodea a Q, se ejercerá una fuerza sobre Q', es como si el espacio que rodea a Q tuviera una nueva propiedad, la de ejercer una fuerza sobre cualquier carga Q' en cualquier posición que la situemos.

Situemos una pequeña carga de prueba q de valor conocido en cualquier posición del espacio dada por el vector de posición v; si vemos que se ejerce una fuerza F (un vector) sobre la carga por el simple hecho de estar situada en ese punto, inferimos que el espacio está afectado por alguna otra carga Q desconocida que crea ese efecto, y entonces se cumplirá la ley de Coulomb,

F = r·k·qQ/r³ [1]

donde r es el vector de posición de q relativo a Q. La figura 1 muestra las cargas q, Q, sus vectores de posición v, w, el vector de posición relativo de q a Q, r = v - w, y la fuerza resultante F, en la dirección de r, que se supone repulsiva (las cargas son del mismo signo y se repelen entre sí).

Fig.1: fuerza electrostática entre dos cargas

Asumiendo el punto de vista de que el espacio alrededor de Q está afectado por esa carga, pues se ejerce una fuerza a distancia sobre cualquier carga q en cualquier posición v, podemos definir un vector E tal que la fuerza ejercida sobre cualquier carga q sea simplemente F = qE a partir de la expresión [1]:

E = k(r/)·Q [2]

Este vector se denomina vector campo eléctrico, y no depende de la carga de prueba q, sólo de la carga causante Q. Si tomamos dos cargas de prueba q, q', y medimos los vectores fuerza F, F', buscando la intersección de las rectas soporte de las fuerzas encontraremos el origen del campo: la posición de la la carga Q que lo genera (figura 2). Además, la magnitud de la fuerza F, conocida la carga q, nos determina también la carga Q a través de la ecuación [1], y por tanto nos determina el campo E.

Fig.2: dos cargas de prueba q, q' determinan el origen del campo eléctrico, la carga Q, y el vector campo E

La propiedad de que todas las rectas soportes de las fuerzas electrostáticas se corten en un punto se puede expresar diciendo que el campo eléctrico es un campo central pues todas las fuerzas parten de un punto central del espacio.

En Física se utiliza el concepto campo para describir cómo se "reparte" una magnitud física medible por el espacio; así, podemos hablar de campos eléctricos, campos magnéticos, e incluso de campos de velocidades en un fluido. Si la magnitud es escalar, el campo lo será, si es vectorial, el campo es vectorial, y si la magnitud es un tensor, el campo será tensorial. La matemática específica para describir campos se ha estudiado en profundidad, dando lugar a la rama de la Física Matemática conocida como Teoría de Campos.

Campo eléctrico producido por varias cargas puntuales

Fig.3: las líneas del campo eléctrico de una única carga puntual son rectas que se cortan en la carga

En el caso de una única carga puntual hemos visto que todas las las rectas soportes de las fuerzas, rectas que llamaremos líneas de fuerza del campo o simplemente líneas del campo  parten de un punto central donde se sitúa la única carga (fig. 3).

¿Qué pasa con las líneas si el campo es generado por dos cargas? Dado que la fuerza electrostática es acumulativa (se suman las contribuciones de todas las cargas) el campo eléctrico E también lo será, y en cada punto del espacio se sumaran los vectores de campo correspondientes a cada carga. Además, la fuerza F y el campo E decrecen con el cuadrado de la distancia, por ello, las líneas del campo se curvan; por cada punto del espacio pasa una línea de campo, de tal modo que el campo eléctrico en ese punto es tangente a la línea. En la figura 4 se representan dos cargas positivas iguales, y una pequeña carga de prueba en la que se suman los dos vectores E, E' generados por las fuentes del campo. Los vectores resultantes no apuntan a ningún centro, el campo resultante deja de ser central, y las líneas de campo (en color negro en la figura) se curvan.

Fig. 4: campo y lineas de campo generadas por dos cargas puntuales iguales

Hay que tener en cuenta que toda carga generará su propio campo eléctrico que se superpondrá a los ya existentes; por ello en lo expuesto hasta ahora hablamos de colocar "cargas de prueba" que se suponen mucho menores que las cargas generadoras del campo, de forma que se puede despreciar su contribución. Además, las cargas fuente, incluso siendo de mucha mayor magnitud que las de prueba, se supone que son de dimensiones puntuales para evitar complicaciones matemáticas.

Campo del dipolo eléctrico

Un caso particular importante es el del denominado dipolo eléctrico, que son dos cargas de igual magnitud q y de signo contrario separadas por una distancia d pequeña (fig. 5).

Fig. 5: líneas de campo producidas por un dipolo eléctrico, fuente: By Geek3 [GFDL (http://www.gnu.org/copyleft/fdl.html) or CC BY-SA 3.0 (https://creativecommons.org/licenses/by-sa/3.0)], from Wikimedia Commons

Se define el momento eléctrico del dipolo, p, por el vector

p = qd, [3]

siendo d el vector que parte de la carga negativa y acaba en la carga positiva. La distancia d se supone que es mucho menor que las distancias a las cuales colocaremos las cargas de prueba donde mediremos el campo E, ello permite simplificar su expresión matemática, que se obtiene de sumar las contribuciones de cada carga al campo total, resultando ser:

\boldsymbol E=qk\left(3\frac{d\cos\left(\theta\right)}{r^4}\boldsymbol r-\frac{\boldsymbol d}{r^3}\right) [4]

donde r es el vector que parte del punto central del dipolo (entre sus dos cargas); alternativamente, usando [3]:

\boldsymbol E=k\left(3p\frac{d\cos\left(\theta\right)}{r^4}\boldsymbol r-\frac{\boldsymbol p}{r^3}\right) [5]


Ejemplo 1: Situamos una carga de +10⁻³C en el origen de coordenadas, y otra carga de -10⁻³C en el punto (1, 0, 0). Calcular el campo eléctrico en el punto P(0.5, 1, 1). Si colocamos en ese punto una pequeña carga de +10⁻⁶C, ¿qué fuerza ejercerá el campo sobre ella?

El campo E' debido a la primera carga será, aplicando [2]:

\begin{array}{l}E'=k\cdot\lbrack{(0.5,1,1)-(0,0,0)}/\sqrt{(0.5²+1²+1²)\rbrack}^3\cdot(10⁻³)=\\\;(0.5,\;1,\;1)\cdot\frac{10⁻³k}{\left(3/2\right)^3}\end{array}

Para la segunda carga:

\begin{array}{l}E''=k\cdot\lbrack{(0.5,1,1)-(1,0,0)}/\sqrt{(0.5²+1²+1²)\rbrack}^3\cdot(-10⁻³)=\\-\;(-0.5,\;1,\;1)\cdot\frac{10⁻³k}{\left(3/2\right)^3}\end{array}

El campo total será la suma de los anteriores:

E=10^{-3}k\left(\frac23\right)^3\left(1,0,0\right).

La fuerza ejercida sobre la carga de prueba q viene dada por  F = qE:

\boldsymbol F=q\boldsymbol E=10^{-6}\cdot10^{-3}\cdot k\left(\frac23\right)^3\left(1,0,0\right)=9\cdot\cancel{10^9}\cdot\bcancel{10^{-9}}\frac8{27}\left(1,0,0\right)=\left(1,0,0\right)

una fuerza de 1 Newton en la dirección del eje X.


Campo eléctrico creado por una distribución de cargas

Cuando en vez de cargas puntuales tenemos cuerpos materiales cargados los modelamos como si contuvieran cargas puntuales distribuidas por el cuerpo, de forma que el campo eléctrico producido por el cuerpo se obtiene sumando las contribuciones de las cargas puntuales; dependiendo de la forma matemática que demos a la distribución, la suma puede ser más o menos directa, sencilla, o complicada. En el caso límite, que es de hecho el habitual, en el que consideremos que hay una infinidad de cargas puntuales será necesario usar el cálculo diferencial y el integral.

Campo eléctrico creado por una distribución de cargas plana y homogénea

El caso más sencillo de distribución de infinidad de cargas puntuales es el de una barra delgada de longitud 2L que tiene cargas sólo en una cara y además están distribuidas de forma uniforme; llamemos ρ a la densidad de carga eléctrica por unidad de longitud, que evidentemente será ρ = Q/2L siendo Q la carga total de la barra. y limitémonos a calcular al campo eléctrico en un punto P situado sobre la bisectriz de la barra, a una altura h:

Fig. 6: geometría para el cálculo del campo E producido por una línea homogénea de carga en un punto situado sobre la bisectriz de la línea

Tomando un elemento diferencial de longitud dx y situado a una distancia x del centro de la barra, por la geometría del problema vemos que la distancia r² será igual a x² + h², y la carga diferencial de ese elemento será dQ =  ρ·dx, por ello el campo diferencial creado por ese elemento en el punto P (punto azul en la figura) tendrá un módulo

\operatorname dE_1=k\frac{dQ}{r^2}=k\frac{\rho\cdot dx}{x^2+h^2}

Dada la situación simétrica el punto P con respecto a la barra, existirá otro elemento dx situado en -x que producirá un campo diferencial dE_2 tal que al sumar los vectores dE_1 con dE_2 se anularan las componentes horizontales y la resultante será vertical, dE con un módulo:

\operatorname dE=\operatorname dE_1\sin\left(\theta\right)+\operatorname dE_2\sin\left(\theta\right)=2\cdot k\frac{\rho\cdot dx}{x^2+h^2}\cdot\sin\left(\theta\right)

Esto sucederá a lo largo de toda la barra, por ello concluimos que el campo resultante ha de ser vertical. Tenemos ahora que sumar todas las contribuciones diferenciales a lo largo de la barra:

E=\int_0^L\operatorname dE=2k\int_0^L\frac{\rho\cdot dx}{x^2+h^2}\cdot\sin\left(\theta\right)

Notemos que los límites de integración son [0, L] y no [-L, L] pues para cada elemento dx a la derecha de la barra (en x) ya hemos sumado la contribución del elemento simétrico a la izquierda (en -x). Teniendo en cuenta que \sin\left(\theta\right)=\frac hr=\frac h{\sqrt{x^2+h^2}} llegamos a la integral

E=2k\int_0^L\frac{\rho\cdot dx}{x^2+h^2}\cdot\frac h{\sqrt{x^2+h^2}}=2kh\rho\int_0^L\frac{dx}{\left(x^2+h^2\right)^{3/2}}

Sin entrar en detalles del cálculo de la integral (podemos calcularla por ejemplo usando WolframAlpha), nos quedará:

E=2kh\rho\frac L{h^2\sqrt{h^2+L^2}}=\frac{2k\rho L}{h\sqrt{h^2+L^2}}

En el caso de que la barra sea muy larga comparado con la distancia h, o sea L >> h, podemos simplificar el valor del campo:

\lim_{L\rightarrow\infty}E=\frac{2k\rho}h\lim_{L\rightarrow\infty}\frac L{\sqrt{h^2+L^2}}=\frac{2k\rho}h [6]

El valor del campo variará según 1/h.


Flujo del campo vectorial. Ańgulo sólido

En el apartado anterior hemos visto que incluso en un caso simple (línea unidimensional de carga, homogénea, campo en un punto de la bisectriz...) cuando calculamos campos debidos a distribuciones continuas de cargas en seguida aparecen integrales complicadas, o muy complicadas. En este apartado y el siguiente vemos un punto de vista de la cuestión basado en la geometria que a menudo simplifica mucho los cálculos,

Consideremos un campo vectorial central, que es aquel que hace corresponder a cada punto P del espacio un vector V que sigue la dirección OP, siendo O el centro del campo, un punto fijo (figura 7).

Fig.7: Campo central, la masa agente situada en O

El caso especial de la forma V = r·c/r² donde r es el módulo de OP (la distancia al centro del campo), r el vector unitario que indica la dirección de OP, y c es una constante que denominamos masa agente del campo; se denomina campo newtoniano. Son campos newtonianos el campo electrico, el campo magnético y el campo gravitatorio. Vamos a definir el flujo del campo vectorial a través de una superfície infinitesimal dS, que denominamos d\phi, como el producto escalar V·dS, donde dS es el vector perpendicular a la superfície; siendo dS infinitesimal (muy pequeña), consideraremos que su curvatura es despreciable, y por tanto es plana (figura 8):

Fig. 8: elemento diferencial de superficie, vector dS, y campo vectorial V que suponemos pasa por el centro de la superficie

\operatorname d\phi=V\cdot dS=V\cdot dS\cdot\cos\left(\theta\right) [7]

Definamos a continuación el ángulo sólido de la superficie dS con respecto al origen del campo central O. Unamos los extremos de la superficie dS con el punto O usando líneas OP, y definamos una esfera C de radio 1 con centro en O; las lineas OP cortaran a la esfera definiendo sobre ella una pequeña superficie dΩ a la que llamamos ángulo sólido de dS sobre C (figura 9).

Fig. 9: ángulo sólido subtendido sobre la esfera C por el elemento de superficie dS

Veamos ahora una propiedad geometrica de los campos newtonianos: el flujo del campo a través de dS es:

\operatorname d\phi=V\cdot dS=V\cdot dS\cdot\cos\left(\theta\right)=\frac c{r^2}dS\cdot\cos\left(\theta\right)=\frac c{r^2}dS' [8]

donde hemos igualado dS\cdot\cos\left(\theta\right)=dS', que es la proyección de la superficie dS sobre la perpendicular al vector de campo V. Las superficies dS' y dΩ (Fig. 10) son paralelas, y están unidas por las mismas rectas al punto central O, se cumple entonces que la razón de sus áreas es igual a la razón de sus distancias al centro O al cuadrado:

Fig. 10: dS' y dΩ son paralelas, y están unidas por las mismas rectas al punto central O

\frac{\operatorname d\phi}{dS'}=\frac1{r^2}\Leftrightarrow dS'=r^2\operatorname d\phi

Si no ve el por qué el lector, piense que el área de la esfera es 4πr², el área de la esfera unitaria es 4π, y el área de una esfera que pase por dS' es 4πr², luego la razón de áreas es 4πr² : 4π = r². Usando esta proporción en la ecuación [8] obtenemos:

\operatorname d\phi=\frac c{r^2}r^2\cdot d\Omega=c\cdot d\Omega [9]

que nos dice que el flujo del vector campo newtoniano a través de una superficie diferencial cualquiera no depende de la distancia r, y es directamente igual al producto de la masa agente del campo por el ángulo sólido subtendido por la superficie sobre la esfera unidad centrada en el origen del campo.

Flujo del campo a traves de superficies cerradas

El producto escalar definido en [7] puede ser positivo o negativo dependiendo de la orientación relativa de los vectores campo V y dS; consideremos una superficie cerrada S, pdemos calcular el flujo total del campo a través de S integrando para cada elemento dS:

\phi=\int_S\operatorname d\phi=\int_SV\cdot dS  [10]

Si el flujo total es positivo, diremos que es un flujo entrante en S, y si es negativo, será un flujo saliente de S. Si el centro del campo O está en el exterior de S, las lineas OP desde el centro que pasen por S cortaran a S en un número par de puntos; en cambio si O está en el interior de S, las lineas cortaran a S en un número impar de puntos (figura 11).

Fig.11: flujos provenientes del centro O a través de superficies cerradas

Por tanto en el caso exterior cada linea creará una sucesión de flujos entrante-saliente-entrante-saliente-etc en número par, o sea, de flujos positivos y negativos alternados; hemos visto que el flujo no depende de la distancia al centro (ecuación [9]) sino sólo del ángulo sólido subtendido, que será el mismo para cada elemento de superficie sobre las lineas que parten de O. Por ello, para un punto O exterior, los flujos entrantes y salientes tienen el mismo valor y se anulan entre sí, resultando un flujo total nulo:

En un campo newtoniano, el flujo total a través de una superficie cerrada que no contiene al centro del campo es nulo.

En cambio si la superficie contiene al centro del campo tendremos un número impar de flujos entrantes-salientes, y su suma no se anulará; de hecho su valor viene dado por el teorema de Gauss, que se deriva de las ecuaciones [9] y [10]:

\phi=\int_S\operatorname d\phi=\int_Sc\cdot\operatorname d\Omega=c\cdot\int_S\operatorname d\Omega=4\pi c  [11]

Expresado en palabras:

En un campo newtoniano, el flujo total a través de una superficie cerrada que contiene al centro del campo es igual a 4\pi multiplicado por el valor de la masa agente.

En el caso del campo eléctrico, la masa agente vale k·Q, siendo Q la carga eléctrica.

Si tenemos un conjunto de cargas, cada carga creará su flujo de campo, y el flujo total será la suma de todos ellos.

Campo creado por una distribución infinita de cargas, plana y homogenea, usando el teorema de Gauss

Como ejemplo de la utilidad del concepto de flujo de campo y del teoriema de Gauss calcularemos el campo electrico creado por una placa cargada uniformemente, que supondremos de extensión muy grande, sobre un punto P situado a una altura h del plano.  Imaginemos otro punto P' situado al otro lado del plano, simétrico a P, y pensemos en un cilindro con eje PP' y de àrea superior dS (en la figura 12 vemos una vista lateral del plano cargado y de la situación).

Fig. 12: usando el teorema de Gauss para calcular el campo eléctrico E

Por simetría el campo E ha de ser vertical y en las direcciones indicadas en la figura; los vectores dS normales a las bases del cilindro tendrán la misma dirección que E, luego el producto escalar resultante, teniendo en cuenta que no habrá flujo de E en las paredes laterales del cilindro por ser paralelas al campo (luego el vector normal a las paredes es perpendicular a E y su producto escalar, nulo) es \operatorname d\phi=E\cdot dS+E\cdot dS=2E\cdot dS.  Si llamamos \sigma a la densidad de carga por unidad de superfície de la placa, la carga encerrada dentro del cilindro será \operatorname dQ=\sigma\cdot\operatorname dS. Por el teorema de Gauss, el flujo ha de ser entonces \operatorname d\phi=4\mathrm\pi\cdot\mathrm k\cdot\mathrm\sigma\cdot\operatorname d\mathrm S; igualando las dos expresiones para el flujo obtenemos:

d\phi=4\pi k\sigma\cdot\cancel{dS}=2E\cdot\cancel{dS}\Leftrightarrow\boxed{E=2\pi k\sigma}

Vemos que la ingensidad de campo E creada por una placa infinita cargada homogéneamente no depende de la distancia h a la placa, un resultado notable.

Problemas

  1. Calcular el campo eléctrico producido por una esfera cargada con densidad de carga homogenea, tanto en el interior de la esfera como en el exterior.
  2. Una esfera  cargada con densidad de carga homogenea tiene una cavidad esférica en su interior, los centros de la esfera cargada y la cavidad estan a una distancia d. Calcular el campo electrico en la cavidad.

 


Soluciones

Problema 1 - En la figura vemos la geometría del problema: representamos una superficie esférica S interior y concéntrica  a la esfera cargada (en azul) y sobre S un punto cualquiera P, por la que trazamos una línea que pasa por el centro O y divide a las esferas en dos mitades simétricas; por simetría, el vector campo E(P) en el punto P no puede estar dirigido hacia ninguno de las dos mitades en particular, así que ha de ser radial. Además, el punto P podría ser cualquier punto situado en S pues tenemos simetría esférica, luego el módulo del campo E será el mismo en todo S, es decir, el valor de E depende sólo del radio de S.

Por simetría del problema respecto cualquier recta que pase por O, el campo E ha de tener el mismo módulo en toda esfera S interior, y ha de ser radial

Vamos a aplicar el teorema de Gauss a la superficie S: primero de todo damos forma matemática a las consideraciones anteriores sobre el campo E:

\overrightarrow E(P)=\frac1rE(r)\cdot\overrightarrow{OP},

o sea el vector campo E es igual al módulo E(r) por el vector radial OP dividido por el módulo de OP, que es r. Equivalentemente podemos definir el vector unitario radial \widehat r=\frac{\overrightarrow{OP}}r y la expresión del campo en todo punto P de S queda más compacto: \overrightarrow E(P)=E(r)\cdot\widehat r.

Ahora calculamos el flujo de E a través de la superfície S, aplicando [10]:

\phi=\int_S\operatorname d\phi=\int_S\overrightarrow E\cdot d\overrightarrow S=\int_SE\left(r\right)\widehat r\cdot d\overrightarrow S=E\left(r\right)\int_SdS=E\left(r\right)\cdot4\pi r^2

donde hemos aplicado que el elemento diferencial vectorial de superfície dS es un vector radial de módulo dS, y por tanto su producto escalar con el vector radial unitario es simplemente dS, además E(r) es constante sobre S, luego puede salir fuera de la integral, y esta integral sobre S del elemento dS es simplemente la superfície de la esfera S. Este flujo que hemos calculado, según el  teorema de Gauss [11], ha de ser igual a \phi=4\pi kQ, donde Q es la carga contenida en el volumen interior a S; llamando \rho a la densidad de carga por unidad de volumen, tenemos que:

Q=\int_V\rho\cdot\operatorname dV=\frac43\rho\pi r^3\Rightarrow\phi=4\pi k\cdot\frac43\rho\pi r^3=\frac{16}3k\rho\pi^2r^3 [12]

Igualando este flujo dado por el teorema de Gauss con el que hemos calculado antes, hallamos el módulo del campo E:

E(r)\cdot4\pi r^2=\frac{16}3k\rho\pi^2r^3\Rightarrow\boxed{E(r)=\frac43k\rho\pi r} [13]

Vemos que la dependencia E(r) es lineal: aumenta linealmente con r. En el sistema internacional de unidades la constante k se expresa en funcion de la denominada permitividad eléctrica del vacío \varepsilon_0, y el campo se reduce a

E(r)=\frac43\frac1{4\pi\varepsilon_0}\rho\pi r=\frac1{3\varepsilon_0}\rho r. [14]

Esta expresión vale para r\leq a, siendo a el radio de la esfera cargada. Para distancias r al centro de la esfera que sean mayores que el radio a el cálculo es muy parecido, sólo que el valor de la carga Q es constante, siendo la carga total de la esfera, Q=\frac43\pi a^3\rho, y al sustituirla en el teorema de Gauss, el flujo total a través de una superficie de radio r > a vale\phi=\frac43\pi a^3\rho\cdot4\pi k, igualando este flujo con el calculado vectorialmente:

\phi=\frac43\pi a^3\rho\cdot\cancel{4\pi}k=E\left(r\right)\cdot\cancel{4\pi}r^2\Leftrightarrow\boxed{E\left(r\right)=\frac43\pi\frac{a^3}{r^2}\rho k} [14b]

que en función de \varepsilon_0 valdrá

E\left(r\right)=\frac43\pi\frac{a^3}{r^2}\rho\frac1{4\pi\varepsilon_0}=\frac{\rho a^3}{3\varepsilon_0r^2}. [15]

O sea que en el interior de la esfera el campo E crece linealmente con la distancia a centro, y en el exterior decrece cuadráticamente:

Variación del campo E con la distancia r al centro en el caso de una esfera de radio a cargada uniformemente

Problema 2 - La complicación de este problema está en ver como tratar la cavidad dentro de la esfera; la forma más fácil consiste en darse cuenta de que la situación es equivalente, en cuanto al cálculo del campo E, a suponer que la esfera, de radio R y centro en O, está cargada uniformemente en su totalidad con una densidad de carga \rho, y que la cavidad es otra superficie esférica, de radio r < R y centro O', a la que añadimos, superponiéndola, otra densidad de carga de igual valor pero signo contrario, -\rho, con ello, la carga neta en la cavidad será cero. Entonces, debido a que el campo electrico es acumulativo, podemos calcular el campo E por superposición del campo creado por toda la esfera de radio R, al que llamamos E_O, más el campo creado por la carga negativa en la esfera de radio r, al que llamamos E'_O. La geometría de este esquema lo vemos en la imagen, donde hemos dibujado un punto P cualquiera dentro de la cavidad, los dos campos generados en ese punto y el campo total.

Vectores E y geometria del problema

Para calcular cada campo aplicamos [13] pues el punto P es interior a las dos esferas consideradas:

\overrightarrow E\left(P\right)={\overrightarrow E}_O\left(P\right)+\overrightarrow E'_O\left(P\right)=\frac43k\pi\rho r\cdot\widehat r-\frac43k\pi\rho r'\cdot\widehat r'=\frac43k\pi\rho\left(r\cdot\widehat r-r'\cdot\widehat r'\right)

Tenemos que r\cdot\widehat r=\overrightarrow r,\;r'\cdot\widehat r'=\overrightarrow r' y observando la figura deducimos que r\cdot\widehat r-r'\cdot\widehat r=\overrightarrow r-\overrightarrow r'=\overrightarrow{OP}-\overrightarrow{O'P}=\overrightarrow{OO}'.  Por tanto nos queda la siguiente expresión para el campo en el interior de la cavidad:

\overrightarrow E\left(P\right)=\frac43k\pi\rho\cdot\overrightarrow{OO}'=\frac43k\pi\rho\cdot\overrightarrow d

Observemos que E es un vector constante dirigido según la recta OO' (el vector d de módulo igual a la distancia d entre centros).

El campo en el interior de una cavidad dentro de una esfera uniformemente cargada es constante

Publicado en Electricidad y Magnetismo, Fí­sica | Etiquetado , , , , , , | Deja un comentario

Electrostática

Tales de Mileto (600 aC) observó que una barra de ámbar (una resina vegetal) frotada con fuerza atraía y levantaba objetos pequeños; en griego el ámbar se llama "elektron", y de aquí vino que se llamara a estos fenómenos como "eléctricos".  William Gilbert (siglo XVI) descubre que sucede lo mismo con otros materiales, a los que llama materiales eléctricos, que hoy en dia llamamos  materiales aislantes, los cuales presentan este fenómeno de electrización por frotamiento.

Dos tipos de electrización

Si frotamos dos varillas de vidrio y las acercamos nortaremos que se repelen mútuamente; lo mismo ocurre si lo hacemos con dos varillas de ámbar. En cambio si acercamos una varilla electrizada de vidrio y una de ámbar, se atraen.  Experimentalmente se encuentra que, dados dos materiales electrizados, algunos se repelen entre sí y otros se atraen; más aún, se pueden agrupar sólo en dos categorias: los que son atraidos por la varilla de vidrio y repelidos por la varilla de ámbar, y los que reaccionan al contrario. B. Franklin (siglo XVIII) llamó a estas categorias, de forma arbitraria, "positiva" (para las cargas similares a la del vidrio) y "negativa" (para las similares al ámbar). La experiencia probaba que electrificaciones del mismo signo se repelen y de distinto signo se atraen.

Creando péndulos con pequeñas masas electrificadas (fig. 1) se puede deducir la fuerza con que se atraen o repelen, midiendo la desviación de la verticalidad y sabiendo las masas.

Fig. 1: Dos péndulos con material electrificados del mismo signo se repelen entre sí

Supongamos que hemos realizado las mediciones de la fuerza F_A que un cuerpo A ejerce sobre otro C de distinto signo, y también la fuerzaF_B de otro cuerpo B sobre C, también de distinto signo (por tanto A y B deben de ser del mismo signo). Podemos conjeturar que si las fuerzas anteriores no son iguales en magnitud, debe de ser por que los materiales A y B no tienen la misma "cantidad de electrificación", cantidad que llamaremos carga eléctrica, y simbolizaremos por Q. Entonces, la fuerza deberá ser proporcional a la carga, y al dividir las fuerzas entre sí se cumplirá la proporción

\frac{F_A}{F_B}=\frac{Q_A}{Q_B} [1]

Si ahora medimos la fuerza que ejercen A y B (unidos por ejemplo con cola) sobre C, encontraremos que es F_{A+B}=F_A+F_B (figura 2).

Fig. 2: Al añadir cuerpos cargados, la fuerza electrica aumenta en el mismo grado

Usando  la proporción [1] con la fuerza resultante obtenemos:

\frac{F_{A+B}}{F_B}=\frac{F_A+F_B}{F_B}=\frac{F_A}{F_B}+1=\frac{Q_A}{Q_B}+1=\frac{Q_A+Q_B}{Q_B} [2]

expresión que nos dice que el efecto de acumular el objeto con carga Q_A más el objeto con carga Q_B es el mismo que produciría un objeto con la carga Q_A+Q_B, es decir, que la carga es una magnitud física escalar, se pueden comparar cantidad de cargas, sumarlas, etc.

Ley de conservación de la carga eléctrica

Supongamos que henos frotado una varilla de vidrio con un paño suave para electrificarla con una carga Q_A que será positiva. Con instrumentos más precisos que el péndulo se puede medir que también el paño se ha electrificado, y lo ha hecho con una carga de igual magnitud pero de signo contrario a la de la varilla, -Q_A: antes de electrificar la varilla no había carga electrica, y después si consideramos el conjunto varilla + paño tampoco, pues la suma algebraica de cargas sigue siendo cero. Es un caso particular de la ley física de la conservación de la carga eléctrica:

En un sistema aislado la carga eléctrica permanece constante.

En general la materia es eléctricamente neutra; si un objeto adquiere una carga positiva será porque otro cuerpo se la ha cedido, adquiriendo carga negativa, y siendo la suma total de cargas cero.

Electrones ligados y electrones libres

Sabemos que la materia está constituida por átomos, con un núcleo cargado positivamente rodeado de una nube de electrones cargados negativamente; el núcleo es muy estable, y sólo puede verse afectado por reacciones nucleares o por desintegración natural, si es un elemento radiactivo, en cambio la nube de electrones es mucho más inestable, y es relativamente fácil extraer electrones de esa nube. Eso es lo que sucede al frotar un material aislante con un paño: algunos electrones, por efecto del frotamiento, "cambian de bando", por ejemplo pasan de la varilla al paño, con el efecto de alterar la neutralidad de cargas, el material que pierde elecrones queda con carga neta positiva, y el que los recibe con carga negativa.

También sabemos que la carga eléctrica del electrón Q_e es la unidad básica de carga, no existen cargas inferiores. Por ello, se podrá expresar cualquier carga Q como un múltiplo de Q_e.

En los metales los átomos estan en un estado tal que algunos de sus electrones puede "saltar" con facilidad de un átomo a otro; la carga total sigue siendo neutra, pero esos electrones no estan fijados en un átomo en particular sino que se van moviendo por el material, por ejemplo debido a la agitación térmica de las moléculas, que las hace vibrar.  Son los electrones libres, o electrones de conducción.  Los demás electrones que están fijados a un átomo son los electrones ligados. Los materiales se consideran conductores si disponen de abundantes electrones libres y por el contrario aislantes si sus electrones son mayoritariamente ligados.

Al frotar el vidrio, la energía cinética que suministramos es suficientemente fuerte para arrancar algunos de sus electrones ligados, que pasan al paño, y el vidrio queda con déficit de electrones; en cambio al frotar ámbar, siendo sus electrones ligados más fijos que el vidrio, sucede lo contrario, algunos electrones del paño saltan al ámbar, cargándolo negativamente.  Si lo intentamos con una varilla de madera no conseguiremos que atraiga pequeños objetos, pues los electrones de la madera estan más ligados que los del vidrio y además no admite electrones externos, como sí hace el ámbar: decimos que la madera es muy aislante pues no se electrifica.

Electrificación por inducción

Fig.3: separación de cargas en el volumen de un conductor

Si acercamos un cuerpo A cargado, por ejemplo positivamente, a un conductor neutro, sus electrones libres seran atraídos cerca de A, formando una región cargada negativamente (figura 3). Habrá una región central de la que habrán salido electrones hacia la izquierda pero que también recibirán electrones de la derecha, quedando neutra, y habrá una región a la derecha que, como no hay más material a su derecha, habrá perdido electrones y quedará cargada positivamente. Si acercamos un tercer cuerpo C por ejemplo con carga positiva, veremos que es repelido por el metal. Al retirar el cuerpo A, los electrones del metal vuelven a dispersarse por todo su volumen, y el cuerpo C dejaría de ser repelido.

Pero si mantenemos A cerca del metal y conectamos el metal con la tierra usando un cable conductor, los electrones libres de la tierra seran atraídos por la carga positiva del metal B y se incorporaran a éste. Si después retiramos tanto el cable como el cuerpo A, el metal quedará con una carga neta negativa distribuida por todo su volumen.

Fig. 4: electrificación de un conductor por inducción

Ley de Coulomb

Coulomb (siglo XVIII) usó una balanza de torsión para averiguar el efecto de la separación entre cuerpos cargados sobre la fuerza que se ejercen (figura 5).

Fig.5: Esquema de la balanza de torsión de Coulomb

Dispuso dos bolitas de saúco electrizadas A, B, la primera fijada a la tapa del cilindro por una varilla aislante rígida y la segunda en el extremo de otra varilla, balanceada por un peso, ámbos colgando de un hilo muy fino de plata. La repulsión entre A y B hace que éste se aleje, retorciendo el hilo. Midiendo la torsión del hilo puede deducirse la fuerza de respulsión. Colocando A en diversas posiciones, pudo estudiar el efecto que tenia la separación r entre A y B. Su conclusión se conoce como ley o fórmula de Coulomb de la electrostática:

F=k\frac{QQ'}{r^2} [3]

La k es una constante de proporcionalidad que depende del medio (aire, vacío, agua, ...), Q y Q' son las cargas y r la distancia que las separa. El hecho de que k dependa del medio indica que hay una transmisión de fuerza electrostática a través del medio, el cual puede ser más o menos permeable eléctricamente.

Inicialmente se propuso como unidad de carga aquella que, siendo Q = Q' = 1, y situadas a una distancia de 1cm produjeran una fuerza atractiva de 1 dina (que a su vez es la fuerza que acelera una masa de un gramo a un cm/s²). En este sistema de medidas cegesimal la constante k vale la unidad, y a relación entre carga eléctrica y las unidades básicas gr, cm, s son poco intuitivas, pues aparecen exponentes fraccionarios:

\left[gr\frac{cm}{s^2}\right]=\left[\frac{Q^2}{cm^2}\right]\Rightarrow\left[Q\right]=\left[\sqrt{gr\frac{cm^3}{s^2}}\right]=\left[\frac{gr^{1/2}cm^{3/2}}s\right]

Se propuso otra unidad de medida, el Coulomb C, derivada de la corriente eléctrica, evitando exponentes racionales: la unidad de carga se define como aquella que es transportada por una corriente de un Ampere en un segundo, [C] = [A/s]. Además, se dió un valor a k adecuado para simplificar expresiones:

El 4\pi evita que salga ese término en muchas formulas del electromagnetismo. Este sistema es el internacional de unidades.

Distribución de cargas

Como hemos mostrado en [2] las cargas son aditivas; entonces si colocamos diversas cargas en el espacio, para cada par de ellas se aplicará la ley de Coulomb [3]; ¿cuál será la fuerza resultante sobre una de las cargas? Siendo la fuerza un vector, necesitamos una versión vectorial de [3]. Denominando \overrightarrow r_{21} al vector que parte de la carga Q_2 y llega a la carga Q_1, podemos escribir:

\overrightarrow F(Q_1)=k\frac{Q_1Q_2}{r_{21}^3}{\overrightarrow r}_{21} [4]

en donde tenemos que especificar sobre qué carga estamos calculando la fuerza (el sentido sobre Q_1 es el contrario que sobre Q_2); observemos que la fuerza sobre Q_1 utiliza el vector r_{21} que va de la carga 2 a la 1. Para cargas del mismo signo este convenio proporciona una fuerza de repulsión entre cargas. Otra expresión equivalente es:

{\overrightarrow F}_{A\rightarrow B}=k\frac{Q_AQ_B}{\overrightarrow{\left\|AB\right\|}^3}\overrightarrow{AB} [4b]

que expresa la fuerza que ejerce la carga situada en A sobre la carga situada en B, siendo Q_A,Q_B las cargas y \overrightarrow{AB} el vector que va del punto A hasta el B.

Alternativamente, utilizando el vector unitario {\widehat r} que tiene la dirección de\overrightarrow r_{12} pero módulo 1, será:

\overrightarrow F(Q_1)=k\frac{Q_1Q_2}{r^2}{\widehat r}_{21} [5]

Para un conjunto Q_1,Q_2,...,Q_n de cargas, la fuerza ejercida por todas ellas sobre la i-ésima carga Q_i será la suma de las fuerzas de cada pareja posible, o sea sumando las contribuciones dadas por [5]:

\overrightarrow F(Q_i)=k{\textstyle\sum_{j\neq i}}\frac{Q_iQ_j}{r_{ji}^2}{\widehat r}_{ji} [6]

donde el subíndice j recorre todos los valores 1 ... n excepto el propio i (no evaluamos la fuerza de Q_i sobre ella misma).

Ejemplo:  Tenemos tres cargas A, B, C todas de valor 10⁻⁴Coulomb;  A y B son positivas, C es negativa, y estan situadas en un plano, en las coordenadas A(0, 0), B(2,3), C(-1, 2). Calcular la fuerza resultante sobre la carga situada en C.

Aplicamos [4b] a la fuerza de la carga A sobre C, y después a la de B sobre C:

De A hasta C: AC = C(-1, 2) - A(0, 0) = (-1, 2);

{\overrightarrow F}_{A\rightarrow C}=9\cdot10^9\frac{-10^{-4}10^{-4}}{5^{3/2}}\overrightarrow{\left(-1,2\right)}=\frac{90}{5^{3/2}}\overrightarrow{\left(-1,2\right)}

De B hasta C: C(-1, 2) - B(2,3) = (-3, -1);

{\overrightarrow F}_{A\rightarrow C}=9\cdot10^9\frac{-10^{-4}10^{-4}}{10^{3/2}}\overrightarrow{\left(-3,-1\right)}=\frac{90}{10^{3/2}}\overrightarrow{\left(-3,1\right)}

Sumamos las fuerzas:

\begin{array}{l}\overrightarrow F\left(C\right)={\overrightarrow F}_{B\rightarrow C}+{\overrightarrow F}_{A\rightarrow C}=-\frac{90}{5^{3/2}}\overrightarrow{\left(-1,2\right)}-\frac{90}{10^{3/2}}\overrightarrow{\left(-3,-1\right)}=\\\frac{90}{5^{3/2}}\left[\overrightarrow{\left(-1,2\right)}+\frac1{2^{3/2}}\overrightarrow{\left(3,1\right)}\right]\cong\overrightarrow{\left(0.48,20\right)}\end{array}

Cuestiones y ejercicios resueltos

  1. ¿La formula de Coulomb se puede utilizar para calcular la fuerza electrostática entre dos cuerpos de dimensiones arbitrarias?
  2. Si acercamos una varilla de ámbar a otra varilla electrificada de vidrio, y vemos que se atraen, ¿debemos suponer que la de ámbar también está electrificada? ¿Y si en vez de atraerse se repelen?
  3. Calcular la proporción que guardan la fuerza de repulsión electrostática y la fuerza gravitatoria entre losprotones  constituyentes de una partícula alfa (idéntica a un núleo de helio; dos protones y dos neutrones). Comentar el valor obtenido, ¿tiene alguna implicación práctica? Datos: masa cada partícula = 1,5·10⁻²⁷kg, carga del protón = 1,6·10⁻¹⁹C, constante de gravitación universal G = 6,7·1⁻¹¹ Nm²Kg⁻¹.
  4. Dos pequeñas esferas electrificadas de masa 0.1gr  cuelgan de un mismo punto por hilos de longitud 1m. Sabiendo que los hilos forman un ángulo entre sí de 15⁰, calcular la carga eléctrica que poseen.

Soluciones

  1. No, sólo tiene sentido para cargas puntuales o de dimensiones despreciables. Para cuerpos cargados ha de desarrollarse una formula más compleja que usa diferenciales e integración.

2. No necesariamente, incluso si el ámbar es neutro, al acercarlo a la varilla cargada puede producirse el fenómeno de la inducción electrostática, que acumule cargas negativas en la parte del ámbar más cercana al vidrio. En cambio si observamos repulsión sí es seguro que el ámbar está cargado con el mismo signo que el vidrio.

3.  La fuerza de repulsión electrostática viene dada por la ley de Coulomb, en la forma no vectorial, donde d es la distancia entre protones:

F_e=9\cdot10^9\frac{\left(1,6\cdot10^{-19}\right)^2}{d^2}

La fuerza de atracción gravitatoria viene dada por la ley de Newton de la gravitación universal:

F_g=6,7\cdot10^{-11}\frac{\left(1,5\cdot10^{-27}\right)^2}{d^2}

Dividiendo la primera por la segunda obtenemos la proporción de ámbas fuerzas naturales:

\frac{F_e}{F_g}=\frac{9\cdot10^9}{6,7\cdot10^{-11}}\frac{\left(1,6\cdot10^{-19}\right)^2}{\left(1,5\cdot10^{-27}\right)^2}\approx0.2\cdot10^{9+11-38+54}\approx10^{35}

Vemos que la fuerza electrostática es prodigiosamente mayor que la gravitatoria. Imaginemos que tenemos dos máquinas de tren de 100 toneladas, que supondremos son de hierro,  cada una separadas a la distancia de 1 metro. La atracción gravitatoria entre ellas es despreciable comparada con la atracción de la Tierra sobre ellas. Supongamos que, por algún procedimiento, pudiéramos conseguir que uno de cada 100 millones de átomos del hierro (Fe) de las máquinas perdiera un sólo electrón, de los 55 que tiene cada uno. Sabiendo que un mol de Fe son 55gr, y que cada mol contiene el número de Avogradro de átomos, 6·10²³, obtenemos que cada locomotora contiene alrededor de 10³⁰ átomos de Fe, y por tanto 55·10³⁰ electrones. Al quitar uno de cada 100·10⁶ = 10⁸, quitariamos un total del 55·10²² electrones; como cada electrón tiene una carga de 1,6·10⁻¹⁹C, habría más carga positiva que negativa, siendo la carga total neta = 55·10²² · 1,6·10⁻¹⁹ = 8800 Coulomb ... esta es una carga enorme. En efecto, las dos locomotoras electrificadas se repelerian con una fuerza de 9·10⁹·(8800)² / 1 = 7·10¹⁷Newton, y serian aceleradas una respecto la otra con una aceleración a = F/m = 7·10¹⁷/10⁵ = 7·10¹², aproximadamente un billón de veces la aceleración de la gravedad: las locomotoras saldrian despedidas y en millonesimas de segundo alcanzarian una valocidad próxima a la de la luz.

4. Ejercicio para el lector. Si quereis la solución, podeis contactar por la página den Facebook.

 

 

Publicado en Electricidad y Magnetismo, Fí­sica | 1 comentario

Problemas del 2º principio

1. En un ciclo de Carnot ¿en qué isoterma se producirá mayor variación de entropía del sistema?

Respuesta: En un clclo cualquiera la variación de entropía del sistema es nula; un ciclo de Carnot incluye dos transformaciones adiabáticas y dos isotermas, el sistema no sufre variación de entropía en las adiabáticas pues Q = 0 y el ciclo es reversible (luego \begin{array}{l}\triangle S=\int\frac{dQ}T\\\end{array}) así que toda la variación de entropía ocurre en las isotermas, pero siendo nula la global, por fuerza las dos isotermas han de producir la misma variación de entropía del sistema, y con signo contrario, para que una anule a la otra.

2. Dos de las siguientes afirmaciones son falsas, y una es incompleta; identificarlas y explicarlas. A) La entropía es constante en toda transformación isotérmica. B) La entropía es constante en toda transformación adiabática. C) La entropía es constante en toda transformación a volumen V constante.

Respuesta: A) es falsa, pues en una transformación isotérmica la temperatura se mantiene constante, pero puede haber transferencia de calor y trabajo, como sucede en los ciclos de Carnot, en los cuales precisamente en las transformaciones isotérmicas es en donde se producen las variaciones de entropía (ver problema 1). B) es parcialmente cierta, pues sólo se cumple en las transformaciones adiabáticas que además son reversibles, pero no en las irreversibles. C) es falsa, pues V = cte implica que no se produce un trabajo, W = 0, pero por el primer principio Q = ΔE, se puede producir una transferencia de calor y puede variar la entropía.

3. Un sistema recorre un ciclo una de cuyas transformaciones es irreversible. ¿De qué signo será la variación de entropía del sistema?

Respuesta: Siendo un ciclo, la variación de entropía S del sistema ha de ser nula, pues S es una función del estado del sistema, y en un ciclo siempre se retorna al estado inicial.

4. Un gramo de agua helada a 0ºC se calienta poniéndola en contacto con una fuente térmica a 200ºC y 1º) se funde, 2º) se calienta hasta 100ºC y 3º) se evapora, siempre a presión constante de 1 atm. De los tres procesos, ¿cuál aumenta más la entropía? Datos: calor de fusión = 80 cal/gr,  calor de vaporización = 540 cal/gr, calor específico: 1 cal/gr·K

Respuesta: la primera transformación se realiza a T = cte, y P = cte, y podemos considerarla reversible, pues en los cambios de fase el sistema siempre está en equilibrio. Entonces, usando los datos y convirtiendo a sistema mks:

A) \triangle S=\int\frac{\operatorname dQ}T=\frac1T\int\operatorname dQ=\frac1T\triangle Q=\frac1{273}1\cdot80\cdot4.1868=1.23\;J/K

La transformación B en cambio es irreversible pues la diferencia de temperaturas entre el agua y la fuente teŕmica aleja a la primera del equilibrio. Imaginemos una transformación reversible que lleve del estado A (T= 0ºC, P = 1 atm) hasta B (T= 100ºC, P = 1 atm), supondremos que el volumen se mantiene constante (una aproximación bastante buena) y que el agua se va calentando poniéndose en contacto con una fuente sólo ligeramente más caliente; si el agua está a T, la fuente estará a T + dT. Entonces,

B) \triangle S=\int\frac{\operatorname dQ}T=\int_{273}^{373}\frac{mc_p\operatorname dT}T=1\cdot1\cdot4.1868\cdot\ln\left(\frac{373}{273}\right)=1.31\;J/K

La transformación C es reversible pues vuelve a ser un cambio de fase, luego:

\triangle S=\int\frac{\operatorname dQ}T=\frac1T\int\operatorname dQ=\frac1T\triangle Q=\frac1{373}1\cdot540\cdot4.1868=6.06\;J/K

Claramente la evaporación es la que aumenta más la entropía.

5. Sean 3·10⁻² moles de un gas ideal en el estado A definido por P=1 atm, V=R m³, (R: constante de los gases perfectos), T = 300⁰K que evoluciona irreversiblemente hasta el estado B definido por P=2 atm, V=4R/3 m³, T = 400⁰K según la transformación 1 que puede verse en el diagrama PV (figura 1). Supondremos el valor 1 para el coeficiente térmico a volumen constante, c_v=1. Calcular la variación de entropía.

Fig.1: diagrama PV de un sistema que evoluciona irreversiblemente entre los estados A, B

Solución: Para calcular la variación de entropía, buscamos una transformación reversible, o sea, que respete en cada instante los estados de equilibrio dados por la ecuación de estado PV = nRT. Imaginemos el camino 2 de A->C que mantiene el volumen constante, luego P/T = nR/V = cte, llevando el sistema hasta la temperatura del estado B; todos los estados en equilibrio a temperatura T_B=400 están sobre la isoterma en rojo, la isoterma en azul corresponde a T = 300. A continuación, imaginamos que el sistema evoluciona reversiblemente desde C hasta A por la isoterma T =400. La variación de entropía entre A y B será la suma de las entropías A->C más la de C->B. Vamos por ello:

S_{AC}=\int_{AC_{REV}}\frac{\operatorname dQ}T=nc_v\int_{300}^{400}\frac{\operatorname dT}T=3\cdot10^{-2}\ln\left(\frac43\right),

para la entropía del paso A->C; sobre la isotérmica, PV = cte, por el 1r principio, Q = W pues el gas perfecto no cambia su energía interna si no lo hace la temperatura; el trabajo W realizado es:

W_{CA}=\int_{CA_{REV}}P\operatorname dV=\int_{CA_{REV}}\frac{nRT}V\operatorname dV=nR\ln\left(\frac{V_2}{V_1}\right)=12R\ln\left(\frac43\right),

la variación de entropía a T = cte es:

S_{CA}=\frac{Q_{CA}}T=\frac{12R\ln\left({\displaystyle\frac43}\right)}{400}=3\cdot10^{-2}R\ln\left(\frac43\right)

y la variación de entropía total:

S_{AB}=S_{AC}+S_{CA}=3\cdot10^{-2}R\ln\left(\frac43\right)+3\cdot10^{-2}\ln\left(\frac43\right)=3\cdot10^{-2}\left(R+1\right)\ln\left(\frac43\right).

6. Un gramo de agua calentado a 100ºC al convertirse en vapor a la presión de 1atm ocupará un volumen de 1671 cm³, absorbiendo 2257J de energía calorífica. Si evaporamos 1Kg de agua, ¿cuál será su aumento de energía interna y de entropía?

Solución: Los cambios de fase (evaporación, solidificación, sublimación ...) son transformaciones de estado especiales pues siendo reales son reversibles (el sistema permanece en equilibrio durante el cambio de fase). En la figura 2 vemos el diagrama PV de la evaporación del agua, a presión y temperatura constante, el volumen específico aumenta linealmente.

Fig.2: diagrama PV de un cambio de fase líquido a gas, caso del agua ( por cm³)

Por el 1r principio: Q = W + ΔE; el trabajo realizado por el vapor correspondiente a 1gr de agua al expandirse bajo una presión constante de  P = 1atm (101.325 Pascales) es W = PΔV = 101.325Pa · (1671 - 1)·10⁻⁶m³ = 169.21 J. El calor absorbido sabemos que vale 2257J, luego ΔE = Q - W = 2257 - 169.21 = 2087.79 J.

El incremento de entropía, al ser un proceso reversible y suceder a temperatura constante, es simplemente ΔS = Q/T = 2257 J / (100+273)K = 6.05 J·K⁻¹ . Siendo la energía interna y la entropía propiedades extensivas (o sea que son propiedades conservativas), para un Kg de agua simplemente multiplicamos por 1000 los resultados anteriores:

  • ΔE total = 1000 · 2087.79 = 2,088 · 10⁶ J
  • ΔS total = 1000 · 6.05 = 6050 J·K⁻¹.

7. Ponemos en contacto 2kg de agua a T=20ºC con una fuente térmica a 100ºC, dejando que al agua se caliente hasta los 100ºC. Calcular las variaciones de entropía del agua, de la fuente y del Universo.

Solución: cuando hay una transferencia de calor debido a una diferencia de temperaturas el proceso es irreversible pues el sistema (el agua) deja de estar en estado de equilibrio mientras se calienta, para que sea reversible ha de calentarse lentamente, cuasiestáticamente, esto es, con una fuente a T = 20C + dT (una diferencia infinitesimal de temperatura), seguida de otra fuente a T = 20C + 2·dT, ... etc.  Como la entropía S es función de estado, imaginamos un proceso así reversible entre los estados inicial y final para calcular la variación de S, y el resultado valdrá para todo proceso entre esos estados.

Sea pues una fuente a T = 20C + dT; como la diferencia de T es infinitesimal, suponemos el proceso isotérmico (T=cte). tendremos una variación dS = dQ/T = cm·dT / T donde c es el coeficiente de absorción de calor del agua, c = 1 cal/gr, m es la masa de agua. Integrando (o sea sumando todos los pasos dT desde 20C hasta 100C):

\triangle S_{agua}=\int_{293}^{373}cm\frac{dT}T4180\cdot1\cdot\ln\left(\frac{373}{293}\right)=1009\;JK^{-1}.

La fuente térmica permanece a T=cte pues se supone que tiene una capacidad calorífica muy grande, entregando una cantidad de calor total Q (por el convenio de signos, será negativa). Si suponemos que la fuente, por tener una capacidad calorífica muy grande, se mantiene en un estado casi estacionario, su proceso se puede considerar reversible, y entonces:

\triangle S_{fuente}=\frac Q{T_{fuente}}=\frac{-mc\triangle T_{agua}}{T_{fuente}}=\frac{-1\cdot4187\cdot80}{373}=-898\;JK^{-1}, [1]

que es negativo pues el calor es cedido. El incremento de entropía de todo el conjunto será la suma: 1009 - 898 = 111 JK⁻¹, que es positivo de acuerdo al 2º principio de la termodinámica.

Si en cambio suponemos que la fuente también está saliendo del equilibrio (capacidad calorífica grande pero no infinita) el proceso será irreversible, y también deberemos aproximarlo por uno de reversible: supongamos que la fuente cede un dQ estando a temperatura T = cte, y siendo dQ el dado por [1]; entonces:

\triangle S=\int\frac{\operatorname dQ}T=\frac{-mc}T\int_{293}^{373}\operatorname dT=-\frac{1\cdot4187}{373}\left(373-293\right)=898\;J/K,

resultado que coincide con el anterior cálculo, donde considerábamos el proceso de la fuente como reversible, era pues una aproximación muy buena.

8. Demostrar que en una transferencia de calor irreversible entre un sistema y una fuente de calor la entropía del universo siempre aumenta.

Solución: Supongamos que el sistema pasa de la temperatura T a la T', siendo T'> T, y que la fuente es isotérmica. Imaginando un calentamiento reversible del sistema, la variación de entropía es, suponiendo que el coeficiente calorífico c no varía con la temperatura:

\triangle S_s=\int\frac{\operatorname dQ}T=\int_T^{T'}\frac{mc\operatorname dT}T=mc\ln\left(\frac{T'}T\right)

La fuente térmica se mantiene a una temperatura T'' > T', y cede de forma reversible un salor Q  (pues permanece en equilibrio al tener una capacidad calorífica virtualmente infinita), luego:

\triangle S_f=\int\frac{-\operatorname dQ}T=\frac1{T''}\int_{}^{}-\operatorname dQ=-\frac{mc\left(T'-T\right)}{T''}

La variación total de S será la suma:

\triangle S_s+\triangle S_f=mc\ln\left(\frac{T'}T\right)-\frac{mc\left(T'-T\right)}{T''}=mc\left[\ln\left(\frac{T'}T\right)-\frac{T'-T}{T''}\right] [1]

Veamos que esta expresión será siempre mayor que cero; sabiendo que T < T' > T'' se sigue que

\frac{T'-T}{T''}<\frac{T'-T}{T'}=1-\frac T{T'} [2]

Comparemos ahora 1-\frac T{T'} con \ln\left(\frac{T'}T\right): llamemos x=\frac{T'}T; f(x)=\ln\left(x\right); g(x) = 1 - 1/x; derivando las funciones f, g vemos que sus pendientes son f'(x) = 1/x, g'(x) = 1/x², ambas funciones son crecientes, pero la segunda  crece a menor tasa pues 1/x² < 1/x. Los valores mínimos de f, g se dan para x = 1 (que se corresponde con T = T') en donde coinciden, f(1) = g(1) = 0. Así pues, f(x) es siempre mayor que g(x) y por tanto [1] será siempre positiva: la entropía del universo siempre aumenta. En la figura se representan las funciones f(x) (azul) y g(x) (rojo), la diferencia en vertical entre las dos líneas es la variación de entropía del universo en función de x = T' / T. En el caso de T'' > T' la diferencia será mayor.

Comparación de la funciones de variación de entropía del sistema y de la fuente térmica

 

 

 

 

 

Publicado en Fí­sica, Termodinámica | Deja un comentario

Entropía y 2º principio de la termodinámica

Se observa que el calor Q fluye de forma natural desde las fuentes térmicas de más temperatura a las de menos, pero al revés no sucede, excepto si lo forzamos gastando un trabajo W. El 1r principio de la termodinámica establece un balance entre el calor entregado a un sistema, el trabajo realizado por el sistema y su incremento de energía interna, Q = W + ΔE, pero no dice nada del sentido de la transferencia de calor, nada prohíbe que fluya de menos a más temperatura. Se necesita un criterio adicional al 1r principio que regule el sentido de las transferencias de calor, y ese será el 2º principio.

Lo que se explica a continuación sigue el siguiente esquema: 1) exploramos la relación calor transferido/temperatura, Q/T,  en un ciclo reversible, 2) visitaremos brevemente el análisis matemático recordando las nociones de campos vectoriales conservativos, circulación de un vector en el campo y potencial del campo, 3) usaremos esos conceptos matemáticos y la relación Q/T de un ciclo para definir la función de estado entropía S, y a continuación 4) veremos las propiedades más importantes de S y sus aplicaciones prácticas en los cambios de estado termodinámicos.

Relación entre calor transferido y temperatura en un ciclo de Carnot

En Trabajo, calor y primer principio de la Termodinámica se estudian los ciclos de Carnot, que son sucesiones continuas de cambios de estado reversibles (en cada instante el sistema está siempre en equilibrio termodinámico) y  cíclicos (estado final igual al estado inicial) formados por dos transformaciones isotermas (temperatura T = cte) y dos adiabáticas (calor transferido Q = 0). En la figura 1 vemos en un diagrama Presión-Volumen las isotermas (transformaciones desde 1 a 2 y de 3 a 4) y las adiabáticas 2->3 y 4->1.

Fig.1: ciclo de Carnot representado en el diagrama Presión-Volumen

Aplicando la relación PV = cte para las isotermas y PV^\gamma para las adiabáticas del ciclo de Carnot:

\begin{array}{l}P_1V_1=P_2V_2\\P_2V_2^\gamma=P_3V_3^\gamma\\P_3V_3=P_4V_4\\P_4V_4^\gamma=P_1V_1^\gamma\end{array}

si las multiplicamos todas entre sí y eliminamos términos repetidos, nos queda

{V_2}V_4=V_3V_1 [1]

Supongamos ahora que el sistema usado en el ciclo es un gas ideal, entonces la ecuación de estado es PV = nRT, y además el trabajo generado en las isotermas será (para detalles ver Trabajo, calor y primer principio de la Termodinámica) W_{12}=nRT\ln\left(\frac{V_2}{V_1}\right),\;W_{34}=nRT'\ln\left(\frac{V_4}{V_3}\right). Dividiendo la primera por la segunda y aplicando [1]:

\frac{W_{12}}{\;W_{34}}=\frac{T\ln\left({\displaystyle\frac{V_2}{V_1}}\right)}{T'\ln\left(\frac{V_4}{V_3}\right)}=\frac{T\ln\left({\displaystyle\frac{V_2}{V_1}}\right)}{-T'\ln\left(\frac{V_3}{V_4}\right)}=-\frac T{T'}.

Tengamos ahora en cuenta que, pr el 1r principio, Q = W_{12},\;Q'=W_{34}, así que nos queda la relación

\frac Q{Q'}=-\frac T{T'}\Leftrightarrow\boxed{\frac QT+\frac{Q'}{T'}=0} [2]

que relaciona de forma simple los calores transferidos con las temperaturas a las que se transfieren. En particular, si prescindimos de los signos (valores absolutos) implica que si T > T' entonces |Q| > |Q'|: se transfiere más calor del foco térmico que está a más temperatura; la diferencia |Q| - |Q'| es igual al trabajo realizado por el ciclo, y si fuera T = T', entonces |Q| = |Q'| y el trabajo W seria nulo. Pero el trabajo es igual al área encerrada en el diagrama PV por el ciclo, y un trabajo nulo implica que no hay área, luego no se puede recorrer un ciclo manteniendo la temperatura constante en todo él, ha de haber una diferencia de temperatura.

Relación entre calor transferido y temperatura en un ciclo reversible

Generalicemos ahora el resultado del apartado anterior al caso de ciclos reversibles cualesquiera, que representamos en la figura 2.

Fig.2: ciclo reversible cualquiera, parcialmente recubierto por ciclos de Carnot

Hemos añadido algunos ciclos de Carnot (recordemos, isoterma seguida de adiabática seguida de isoterma seguida de adiabática) tales que recubren parcialmente el área del ciclo C, representado en el diagrama PV por una elipse, pero puede ser cualquier figura (con alguna restricción geométrica que no enunciaremos); los ciclos de Carnot tienen las curvas adiabáticas muy pequeñas, de forma que se van ajustando más o menos a la curva C, mientras que las isotermas son más largas y recorren de lado a lado el interior de C, Podemos imaginarnos que dibujamos más ciclos Carnot hasta recubrir en su totalidad toda el área C. Además, observando la figura 1, como cada ciclo de Carnot está adyacente a otro ciclo, podemos pensar que el calor cedido Q' por cada ciclo es igual al absorbido Q por el siguiente, pero no son exactamente iguales, pues cada ciclo de Carnot tampoco coincide exactamente con su contiguo (fig. 3), de forma que el calor cedido Q' por el i-èsimo ciclo de Carnot no es exactamente igual al absorbido por el (i+1): \textstyle Q'_i=Q_{i+1}+\triangle Q_{i,i+1}, hay una diferencia.

Fig.3: detalle de dos ciclos contiguos de Carnot  del recubrimiento de un ciclo C

Además, para cada ciclo de Carnot se cumplirá [2]: \frac{Q_i}{T'_i}+\frac{Q_i'}{T_i'}=0, y la suma anterior extendida a todos los n ciclos también dará cero:

\sum_i\left(\frac{Q_i}T+\frac{Q'_i}{T_i'}\right)=0 [3]

Llevando al límite para n\rightarrow\infty ciclos de Carnot, C quedará recubierto de forma que \triangle Q_i\rightarrow\operatorname dQ, las adiabáticas de cada ciclo de Carnot tenderán a coincidir con el elemento diferencial de longitud dl de la curva C, y la suma [3] pasará a ser una integral extendida a todo el contorno C:

\oint_{C,rev}\frac{\operatorname dQ}T=0 [4]

donde hemos dejado de distinguir calor absorbido Q del emitido Q' y simplemente tenemos un dQ transferido (en más o en menos), y la temperatura T es una variable que toma valores distintos en cada punto de C; también hemos puesto el subíndice rev para indicar que el ciclo se recorre de forma reversible. A continuación hacemos una breve incursión por el análisis matemático, que nos da herramientas potentes para manejar las integrales [4].

Circulación y potencial de un campo conservativo

Un campo vectorial es una función que asigna a cada punto P un vector v(P), como por ejemplo el campo eléctrico, o la asignación a cada punto del espacio (x, y, z) del interior de un gas de una magnitudes (P, V, T) entendidas como un vector en un espacio vectorial. Imaginemos un camino cualquiera (una trayectoria entre dos puntos A, B)  en el espacio a través de una curva C; para cada uno de sus puntos P, existirá un vector campo v(P), y un elemento diferencial dl que será tangente a la curva en P:

Fig.4: curva C, elemento diferencial dl tangente a C en P, y vector campo v(P) en P

Denominamos circulación del campo v(P) entre los puntos A, B a lo largo de la curva C a la suma de los productos escalares v(P)·dl extendida a todos los puntos P entre A y B:

Circ_C\left(\overrightarrow v,A,B\right)={\int_A^B}_C\overrightarrow v\operatorname d\overrightarrow l [5]

En el caso de que C sea un camino cerrado, escribimos la integral así:

Circ_C\left(\overrightarrow v\right)=\oint_C\overrightarrow v\operatorname d\overrightarrow l[6]

Hay un tipo especial de campos vectoriales que en Física son de gran interés, son los denominados campos conservativos, en los que se cumple la condición siguiente para todo camino cerrado C:

Circ_C\left(\overrightarrow v\right)=\oint_C\overrightarrow v\operatorname d\overrightarrow l=0 [7]

Se llaman "conservativos" pues, por las propiedades de las integrales (ver figura 5):

\oint_C\overrightarrow v\operatorname d\overrightarrow l=\int_{PQR}\overrightarrow v\operatorname d\overrightarrow l+\int_{RSP}\overrightarrow v\operatorname d\overrightarrow l=0\Rightarrow\int_{PQR}\overrightarrow v\operatorname d\overrightarrow l=-\int_{RSP}\overrightarrow v\operatorname d\overrightarrow l,

además, invirtiendo el sentido de integración,

\int_{RSP}\overrightarrow v\operatorname d\overrightarrow l=\int_{PSR}-\overrightarrow v\operatorname d\overrightarrow l

y por tanto:

\int_{PQR}\overrightarrow v\operatorname d\overrightarrow l=\int_{PSR}\overrightarrow v\operatorname d\overrightarrow l,

que nos dice que el valor de la circulación de un campo conservativo sólo depende de los puntos inicial P y final R, y no del camino concreto recorrido entre esos dos puntos. [7']

Fig.5: un camino cerrado C y algunos de sus puntos, P, Q, R, S

Debido a esta propiedad de conservar el valor de la circulación entre dos puntos independientemente del camino recorrido, en los campos conservativos se puede la magnitud diferencia de potencial V_{PR} como

V_{PR}=\int_{PR}\overrightarrow v\operatorname d\overrightarrow l [8]

para cualquier trayectoria entre los puntos P y R. Tomando un punto cualquiera O como referencia, el potencial en otro punto P se define como V_P=V_{PO}=V_P-V_O.

Los potenciales tiene gran interés en física; así por ejemplo, un campo vectorial eléctrico es conservativo, y por tanto se puede definir un potencial eléctrico en cada punto del campo, además, la corriente eléctrica solo fluye en el sentido de más potencial a menos potencial.

Entropía

Por analogía, siguiendo el apartado anterior y comparando con el resultado [4] obtenido para la relación Q/T en un ciclo reversible, podemos igualar dQ/T al producto escalar de un cierto campo vectorial v por el elemento diferencial de longitud del camino C; siendo nula la integral para todo C reversible y cerrado, concluimos que ese campo v es conservativo, luego admite un potencial que se calculará según [8], que en el caso termodinámico será

S_{AB}=S_A-S_B=\int_{AB_{REV}}\frac{\operatorname dQ}T [9]

A la cantidad S_A la llamamos entropía del estado A, y seria el potencial del campo en ese punto; de hecho no nos interesa el campo en sí, que sólo es instrumental, sino el potencial, la entropía. Para cada estado (P, V, T) queda definido una entropía (resultado [7']) , de forma que podemos decir que la entropía es una función de estado: dos estados distintos tienen entropía distinta.

Observar que en la integral [9] entre los estados A, B hemos incluido el subíndice REV para indicar que la transformación A->B que se use para el cálculo ha de ser reversible. [9a]

Propiedades inmediatas de la entropía S

  1. Como S es función de estado, la variación de entropía de un sistema que evoluciona siguiendo un ciclo cerrado cualquiera, tanto si es reversible como si es irreversible, es nula.
  2. Como S es función de estado, en una transformación de un sistema desde el estado A hasta el B, la variación de entropía será la misma tanto si se hace de forma reversible como si se hace de forma irreversible.

Transformaciones internas reversibles de sistemas aislados

A continuación aplicaremos la entropía como herramienta para comprender mejor qué sucede en las transformaciones de los denominados  sistemas aislados, que son aquellos que no reciben ni entregan calor ni trabajo, siendo todas sus transformaciones internas. Imaginemos un sistema así, como por ejemplo un recipiente con aislamiento térmico en el que hemos introducido un líquido A a temperatura T_1 superior a la T_2  del aire que contiene el resto del volumen B del recipiente, y acto seguido lo cerramos. Como hay una diferencia de temperatura interna, habrá un flujo de calor Q (fig. 6)

Fig.6: Sistema aislado del exterior, con una transferencia interna de calor Q

Llamando Q_A al calor absorbido por A, y Q_B al absorbido por B, en general en un sistema aislado en el que no se intercambia calor con el exterior, Q_A+Q_B=0. Por el 1r principio de la termodinámica, Q = W + ΔE, siendo W = 0, tendremos también ΔE = 0, o sea \triangle E_A+\triangle E_B=0.

Supongamos que el intercambio de calor es reversible: implica que la diferencia de temperatura entre A y B es tan pequeña que en ningún momento ni A ni B salen de los estados de equilibrio, supondremos T = T'; entonces los incrementos diferenciales de entropía (reversible) de cada subsistema A, B son, aplicando [9]:

\operatorname dS_A=\frac{\operatorname dQ_A}T,\;\operatorname dS_B=\frac{\operatorname dQ_B}T,

pero ha de ser \operatorname dQ_A=-\operatorname dQ_B al ser reversible, y la variación de entropía total es:

\operatorname dS=\operatorname dS_A+\operatorname dS_B=\frac{\operatorname dQ_A}T+\frac{\operatorname dQ_B}T=0,

que puede enunciarse así:

en un sistema aislado que evoluciona reversiblemente, la variación de entropía interna es nula.

En los sistemas reales no podremos en general hacer una suposición tan fuerte, veamos que sucede.

Segundo principio de la termodinámica

En toda transformación irreversible de un sistema aislado, la entropía aumentará.

En las transformaciones reales, que no podemos suponer que mantienen constantemente el equilibrio, incluso si ocurren en un sistema aislado que no transfiere ni recibe calor ni trabajo de su entorno, producen un aumento de la entropía del sistema. En particular este principio se relaciona con el hecho observado que mencionamos al principio, que el calor Q fluye de forma natural desde las fuentes térmicas de más temperatura a las de menos, pero al revés no.

Podemos ver una justificación (que  no una demostración) considerando el caso de un gas ideal que evoluciona  de forma irreversible (ver el diagrama PV en la fig. 7) desde el estado A, definido por (P_A, V_A, T_A) hasta el estado B (P_B, V_B, T_B), de más presión temperatura y volumen; para saber el cambio de entropía entre esos dos estados consideramos la transformación alternativa A->B'->B: la primera parte A->B' es reversible adiabática, la segunda parte B'->B es reversible isotérmica. La variación de entropía será

S_{AB}=S_{AB'}+S_{B'A}=0+\frac Q{T_B}

pues en una adiabática reversible no varía S. El signo de la ΔS depende pues del calor Q absorbido o emitido en la transformación isoterma B'B; por el 1r principio, Q = W + ΔE = Q, pues ΔE en el gas ideal depende exclusivamente de la temperatura, así pues el signo de Q será igual al del trabajo: positivo si el gas realiza trabajo,

Fig. 7: transformación irreversible y camino reversible alternativo

El recorrido completo A-B'-B produce un trabajo W > 0. Por el primer principio, Q = W + ΔE, siendo Q el calor absorbido total, pero sólo en la isoterma, ya que en la adiabática Q = 0,  y ΔE la variación de energía interna total, pero sólo en la adiabática, ya que en los gases perfectos ΔE depende exclusivamente de la temperatura). La variación de entropía total será

S_{AB}=S_{AB'}+S_{B'A}=0+\frac Q{T_B},

ya que A-B' es adiabática reversible, luego Q_{AB'} = 0 y no hay variación de entropía.  Para calcular Q_{B'B} aplicamos el 1r principio, Q = W + ΔE, siendo:

W_{B'A}=\int_{B'}^BP\operatorname dV=nrT\ln\left(\frac{V_B}{V_{B'}}\right)>0

y ΔE = 0 pues T es constante, luego Q_{B'B}>0 y la variación de entropía total S es positiva: aumenta.

El 2º principio afirma que en una transformación irreversible S aumentará incluso si no hay transferencia de calor Q, seria el caso de un proceso adiabático irreversible; llamemos ΔS' a este aumento sin transferencia de calor, si a este proceso la añadiéramos además una transferencia de calor, tendríamos además un ΔS'' debido ese calor Q, por ello, en general podemos decir que el ΔS de toda transformación irreversible puede descomponerse en dos factores:

ΔS = ΔS' + ΔS'' (entropía intrínseca irreversible + debida a Q) [10]

Procesos isoentrópicos son aquellos procesos en los que la entropía S no varia. S es función de estado, y eso significa que cada estado ha de tener un valor de S fijado, pero al mismo tiempo podemos tener varios estados que tengan la misma entropía; un proceso que se mueva entre esos estados será isoentrópico. Un sistema aislado que se transforma irreversiblemente no puede ser isoentrópico, pero si es reversible seguro que lo será.

Diagramas TS

Representan los estados y su evolución usando como coordenadas la temperatura y la entropía. En la figura 8 vemos un ciclo de Carnot para un gas perfecto en el diagrama TS, se ven como rectángulos: las líneas horizontales son las isotermas, en las que se transfiere calor y por tanto hay variación de entropía, las líneas verticales son las adiabáticas con Q = 0, la entropía no varía pero la temperatura sí. Como en el diagrama PV, el área interior al ciclo en el TS equivale al valor de trabajo realizado por el ciclo; en efecto, \operatorname dW=P\operatorname dV, y además \operatorname dS=\frac{\operatorname dQ}T\Rightarrow T\operatorname dS=\operatorname dQ=\operatorname dE+\operatorname dW. Integrando y recordando que la energía interna no varia en un ciclo: \oint T\operatorname dS=\cancel{\oint\operatorname dE}+\oint\operatorname dW\Rightarrow W=\oint T\operatorname dS

Fig.8: un ciclo de Carnot visto en el diagrama TS

 

Máquina térmica perpetua

Se llama así a una máquina (térmica en estos apuntes) teórica con una eficiencia del 100% de forma que se puede invertir todo el trabajo que genera en producir calor para retroalimentarla y que siga funcionando de forma autónoma y perpetua, o sea que no necesita una fuente térmica externa (fig. 8), tal máquina tendría dos partes, una que se encarga de transformar todo el calor Q en trabajo w, y otra que transforma todo el trabajo W de nuevo en calor Q, sin pérdidas).  Se desarrollaron diversos intentos de construirla, todos fallidos.

Fig.8: esquema de una máquina térmica perpetua, que constituye un sistema aislado que se retroalimenta

Veamos que según el 2º principio es imposible tal máquina; en efecto, por ser S una función de estado, su variación es cero en todo ciclo tanto reversible como no reversible, pero una máquina autónoma y perpetua que no interacciona con nada (ni tiene fuentes térmicas ni intercambia trabajo con el exterior) es por definición un sistema aislado, y el 2º principio afirma que todo sistema aislado que evoluciona de forma irreversible ha de aumentar su entropía: llegamos a una contradicción, luego no existe tal máquina térmica perpetua.

Teorema de Carnot

Fig.9: esquema para demostrar el teorema de Carnot con dos máquinas de Carnot conectadas entre sí

Dice así:

el rendimiento de una máquina térmica cualquiera sólo depende de las temperaturas extremas T, T' con las que trabaja.

Equivalentemente, dos máquinas térmicas que trabajen entre las mismas temperaturas extremas tendrán el mismo rendimiento.

Consideremos dos máquinas A y B acopladas que trabajan entre los mismos focos térmicos T > T', en el esquema de la fig. 9 vemos las transferencias de calor y trabajo; el acoplamiento consiste en que el trabajo suministrado por la máquina A se utiliza para hacer funcionar la máquina B, que actúa en modo inverso, o sea como bomba de calor, extrayéndolo del foco frío e inyectándolo en el foco caliente. Ajustemos ambas máquinas de tal forma que Q_A=Q_B y Q'_A=Q'_B; esto equivale a modificar sus ciclos hasta que las áreas interiores sean iguales, pues por el 1r principio para cada máquina Q - Q' = W (fig. 10).

Fig.10: dos ciclos entre distintas isotermas y adiabáticas, ajustados para tener la misma área W

 

Entonces será W = W', además, W = Q_A-Q'_A y W'=Q_B-Q'_B luego Q_A-Q'_A=Q_B-Q'_B. Recordemos que los rendimientos de las máquinas son r_A=W/Q_A,r_B=W'/Q_B o sea  r_A=W/Q_A,r_B=W/Q_B , y vamos a suponer que A tiene más rendimiento que B, r_A > r_B, que implica W/Q_A>W/Q_B, que equivale a Q_A<Q_B: la máquina A toma menos calor de la fuente caliente que el calor devuelt por B.

Fig.11: combinación de máquinas térmicas equivalente a una que viola el 2º principio

Entonces, podemos prescindir de la fuente térmica T acoplando las máquinas como muestra la figura 11; el conjunto A+B toma un calor Q'_B-Q'_A  de un foco T' y devuelve un calor neto Q_B-Q_A>0, obviamente este calor podría pasarse a una tercera máquina C que produciria un trabajo neto  W'' con lo cual el sistema completo A+B+C (linea de trazos en la figura 11) produciria un trabajo neto W''' con un único foco térmico: tendríamos una máquina térmica perpetua en contradicción con el 2º principio.

Fig.12: con rendimentos iguales, la combinación propuesta equivale a una máquina nula que no transfiere calor ni trabajo y por tanto no vulnera ningún principio

 

La contradicción se deshace si consideramos que los rendimientos de A y B son iguales, en este caso la máquina A + B simplemente no hace nada (figura 12) pues extrae y devuelve la misma cantidad de calor de los focos T y T' sin aportar nada de trabajo.

 

Rendimiento de máquinas reales irreversibles

Consideremos de nuevo una máquina térmica M operando entre dos fuentes T > T' pero de forma irreversible; el conjunto máquina + fuentes térmicas será un sistema aislado, con transformaciones internas irreversibles, luego por el 2º principio su entropía S_{total} debe aumentar, y será igual al incremento de entropía de la máquina M más los incrementos de entropía de las fuentes térmicas: \triangle S_{total}=\triangle S_M+\triangle S_T+\triangle S_{T'}. Pero la entropía S_M se mantiene constante en toda transformación cíclica, por ser S función de estado, luego \triangle S_{total}=\triangle S_T+\triangle S_{T'}>0. Para calcular estas entropías consideremos una cesión de calor muy lenta, infinitesimal, dQ, de cada una de las fuentes; teniendo en cuenta que peramenecen a temperatura constante, y tomando la definición de entropía, obtenemos:

\operatorname dS_T=\frac{\operatorname dQ}T,\;\operatorname dS_T'=\frac{\operatorname dQ'}{T'}\Rightarrow\triangle S_T=\int\frac{\operatorname dQ}T=\frac QT,\;\triangle S_{T'}=\int\frac{\operatorname dQ'}{T'}=\frac{Q'}{T'}

luego deducimos que

\frac QT+\frac{Q'}{T'}>0\Leftrightarrow\left|\frac{\displaystyle Q}{\displaystyle T}\right|\neq\left|\frac{\displaystyle Q'}{\displaystyle T'}\right|

Comparando con la expresión [3] vemos que en general,

\frac QT+\frac{Q'}{T'}\geqslant0

siendo válida la igualdad para procesos reversibles y la desigualdad estricta para los irreversibles. Recordemos ahora la expresión [10] y la aplicamos a un ciclo irreversible; para todo ciclo es ΔS = 0, luegoΔS' + ΔS'' = 0, siendo ΔS' debida a la irreversibilidad y ΔS'' debida al intercambio de calor Q, esta última será igual a \frac QT+\frac{Q'}{T'}, luego:

\frac QT+\frac{Q'}{T'}+\triangle S''=0,

despejando Q':

Q'=-T'\left(\triangle S''+\frac QT\right),

como en un ciclo ha de ser W = Q + Q' nos queda:

W=Q-T'\left(\triangle S''+\frac QT\right)=Q\left(1-\frac{T'}T\right)-T'\triangle S'',

entonces el rendimiento del ciclo irreversible será R = W/Q:

R=\left(\frac{T-T'}T\right)-\frac{T'\triangle S''}Q=r-\frac{T'\triangle S''}Q

donde r es el rendimiento de un ciclo de Carnot (reversible por tanto), vemos que el rendimiento irreversible R es siempre menor que el reversible r: los rendimientos de las máquinas reales son inferiores a los de las máquinas ideales reversibles.

 

 

 

Publicado en Fí­sica, Termodinámica | Deja un comentario

Problemas de trabajo, calor y 1r principio de la termodinámica

Problemas resueltos

1 - Determinar el trabajo de expansión realizado por un mol de una sustancia al pasar del estado sólido al estado líquido bajo la presión de 1 atm, sabiendo que su volumen específico en el estado líquido supera al del estado sólido en 0,16 cm³/gr. Otros datos: masa de un mol de la sustancia: 60gr/mol.

Solución: el trabajo de expansión viene dado por (ver Trabajo, calor y 1r principio de la termodinámica. ecuación [1]):

W=\int_{V_1}^{V_2}P\operatorname dV[1]

La fusión del sólido ocurre a presión constante, luego

W=\int_{V_1}^{V_2}P\operatorname dV=P\int_{V_1}^{V_2}\operatorname dV=P\left(V_2-V_1\right)=P\triangle V[a]

¿Cuál es el incremento de volumen de un mol cuando se funde? Tenemos 60gr por mol, y el incremento de volumen por gramo es 0,16 cm³/gr, luego el incremento de volumen total es 60gr/mol · 0,16 cm³/gr = 9.6 cm³/mol. Sólo queda pasar todas la unidades a sistema internacional y aplicar [a]:

1 atm = 1. 013·10⁵ N/m²; 9.6 cm³ = 9.6 · 10⁻⁶m³;

W = 1. 013·10⁵ N/m² · 9.6 · 10⁻⁶m³ = 0.972 Nm = 0.972 J (en Joules).


2 - Se comprime de forma reversible e isotérmica un mol de un gas ideal a T = 18⁰C desde la presión inicial de 0.8 atm hasta la presión final de 3.5 atm. ¿Qué trabajo se ha realizado sobre el gas?

Solución: El trabajo de expansión viene dado por (ver Trabajo, calor y 1r principio de la termodinámica. ecuación [1]):

W=\int_{V_1}^{V_2}P\operatorname dV[1]

En nuestro caso es una compresión, luego el trabajo será negativo (producido sobre el gas, y no producido por el gas). Siendo una gas ideal, se cumple la ecuación de estado PV = nRT, y como T es constante, se cumple PV = cte, o sea que P_1V_1=P_2V_2. Además, P=nRT/V, sustituyendo en [1]:

W=\int_{V_1}^{V_2}\frac{nRT}V\operatorname dV=nRT\int_{V_1}^{V_2}\frac{\operatorname dV}V=nRT\left(\ln\left(V_2\right)-\ln\left(V_1\right)\right)=nRT\ln\left(\frac{V_2}{V_1}\right)

Expresamos el volumen en función de la presión usando P_1V_1=P_2V_2\Leftrightarrow\frac{V_2}{V_1}=\frac{P_1}{P_2}, y el trabajo en función de las presiones queda

W=nRT\ln\left(\frac{P_1}{P_2}\right)

Sustituyendo los valores numéricos:

0.8 atm = 0.8 · 1. 013 · 10⁵ N/m² = 0.8104 · 10⁵ N/m² ,

3.5 atm = 3.5 · 1. 013 · 10⁵ N/m² = 3.5455 · 10⁵ N/m² ,

18⁰C = 18 + 273⁰K = 291⁰K,

W=1\cdot8314\cdot291\cdot\ln\left(\frac{0.8104}{3.5455}\right)=-3496\;J.


3. Un gas ideal encerrado en un émbolo a 25⁰C y presión inicial P = 2 atm (estado A) se calienta suministrándole un calor Q sin variar el volumen (émbolo fijado). A continuación, una vez absorbido el calor por el gas (estado B) se libera el émbolo permitiendo que el gas se expanda adiabáticamente, de forma que la presión final es 1,1 atm y la temperatura final la misma que la inicial, 25⁰C (estado C). Calcular la presión y la temperatura en el estado intermedio B, y el calor suministrado. Se supone que todos los cambios son reversibles. Datos adicionales: masa del gas m = 1Kg, calor específico molar a volumen constante c = 5R/2, siendo R la constante universal del los gases perfectos.

Solución. En la imagen 1 vemos los estados del sistema y su evolución en un diagrama PV. La línea vertical AB es a volumen constante, y el sistema absorbe un calor Q, aumentando su presión; por el punto A se ha dibujado en discontinuo la curva isoterma correspondiente a la T = 25⁰C, observemos que el estado final C, al estar a la misma temperatura, está también sobre la isoterma, debe de ser así pues al ser todo el proceso reversible, el estado C es de equilibrio, y al estar A, C  a la misma temperatura, han de estar unidos por una isoterma. La curva BC es adiabática de temperatura variable.

Fig.1: las transformaciones del problema en el diagrama PV

Resumimos lo que sabemos y lo que desconocemos:

  • Estado A: PA = 2 atm, VA = ? , TA = 25⁰C
  • Estado B: PB = ?, VB = VA, TB = ?
  • Estado C: PC = 1.1 atm, VC = ?, TC = 25⁰C
  • Transformación A->B a volumen constante
  • Transformación B->C adiabática Q = 0

Para el estado A es immediato calcular VA, usando PV = nRT -> VA = nRTA / PA = (100·R·295) / (2·1.013·10⁵) = 295R/2026 m³. Luego también tenemos VB = VA = 295R/2026 m³.

Lo mismo sucede con el estado C: VC = nRTC / PC = (100·R·295) / (1.1·1.013·10⁵) = 2950R/11143 m³. Los estados A y C quedan totalmente determinados, y del B nos falta la presión y el volumen.

Transformación B->C. Es una adiabática, aplicaremos lo visto en Trabajo, calor, 1r principio de la Termodinámica apartado Transformaciones adiabáticas; aplicación a los gases perfectos: \mathrm{PV}^\mathrm\gamma=\mathrm{cte} siendo \gamma=c_p/c_v y además se cumple la relación de Mayer, c_p-c_v=R. Con  los datos del enunciado, encontramos el valor de la constante gamma:

\mathrm\gamma=\frac{{\mathrm c}_\mathrm p}{{\mathrm c}_\mathrm v}=\frac{\mathrm R+{\mathrm c}_\mathrm v}{{\mathrm c}_\mathrm v}=\frac{\mathrm R+{\displaystyle\frac{5\mathrm R}2}}{\displaystyle\frac{5\mathrm R}2}=\frac75=1.4

Por tanto, en la transformación B->C, PV^{1.4}=cte, y podemos escribir:

P_BV_B^{1.4}=P_CV_C^{1.4}\Rightarrow P_B=P_C\frac{V_C^{1.4}}{V_B^{1.4}}=1.1\frac{2950\cdot2026}{295\cdot11143}=2.31\;atm

Nos falta TB . Veamos la transformación A->B.

Transformación A-B. Siendo un gas ideal se cumple PV = nRT, y como el volumen es fijo, se cumple la ley de Gay-Lussac (que se deduce directamente de la ecuación PV=nRT considerando el volumen constante):

\frac{P_A}{T_A}=\frac{P_B}{T_B}\Leftrightarrow T_B=T_A\frac{P_B}{P_A}=298\frac{2.31}2=344,

Las presiones no las hemos convertido a sistema internacional pues el cociente \frac{P_B}{P_A} vale lo mismo en atmósferas que en Pascales. Ahora podemos obtener el calor absorbido por el sistema, que es absorbido a volumen constante; el incremento de temperatura viene dado por:

Q=nc_V\triangle T=100\cdot\frac528,314\cdot\left(344-298\right)=95611J\;\simeq\;96KJ.,

o en kilocalorias, aproximadamente 23 Kcal.


4. Una máquina térmica que sigue un ciclo de Carnot trabaja entre las temperaturas de 150⁰C y 50⁰C. ¿Cuál es la cantidad de calor que absorbe si la máquina tiene una potencia de 75 KW?

Solución: En la teoría Trabajo, calor, 1r principio de la Termodinámica apartado Máquina térmica de Carnot para gases perfectos se muestra que "el ciclo de Carnot opera entre dos fuentes térmicas a temperaturas distintas T, T', absorbiendo calor Q de la fuente más caliente y cediendo calor Q' a la fuente más fría, y realizando un trabajo neto W = Q - Q'."; una potencia de 75 KW equivale a desarrollar un trabajo de 75 KJ por segundo, luego Q - Q' = 75000 J. Además, sabemos que el rendimiento de los ciclos de Carnot cumplen:

r=\frac WQ=\frac{T-T'}T

En nuestro caso,

r=\frac{75000}Q=\frac{\left(150+270\right)-\left(50+270\right)}{\left(150+270\right)}=\frac5{21}\Rightarrow Q=\frac{21}575000=315000J

Hay que suministrar 315 KJ por segundo, o sea 315 KW en forma de energía térmica a la máquina para obtener 75 KW de potencia; la máquina tiene un rendimiento de 5/21, en tanto por ciento es de solo el 23.8%: de cada 100 J suministrados en forma de calor, 23.8 J se transforman en calor, y el resto, 76.2 J, se devuelven en forma de calor no aprovechable. El rendimiento de las máquinas térmicas es pues bajo, y mucho menor que el de los motores eléctricos, que rondan el 80% de eficiencia.


Preguntas cortas

1 - Un gas perfecto se dilata reversiblemente en dos fases, la primera es isoterma, la segunda adiabática, recorriendo tres estados, el inicial A, el intermedio B y el final C. Tenemos conocimiento completo de les estados A y C, pero no sabemos nada del estado intermedio B. ¿Es posible calcular el trabajo de expansión realizado por el sistema?

Respuesta. No es posible, y no lo es porque el trabajo realizado en una transformación de diversos estados depende de cada uno de los estados intermedios. En la figura 2 vemos en un diagrama PV dos posibles caminos ABC y AB'C entre los mismos estados A y C, recordando que el trabajo realizado coincide con el área bajo la curva, vemos claramente que las áreas delimitadas por los dos caminos son distintas, luego conocer el estado intermedio B es necesario para calcular el trabajo.

Fig. 2: el trabajo realizado en una transformación de estados depende de todos los estados, no sólo del inicial y el final


2 - Un gas perfecto se dilata reversible y adiabáticamente desde el estado inicial A al final B. Conocemos las temperaturas en A y B pero  no los volúmenes. ¿Es posible calcular el trabajo de expansión realizado por el sistema?

Respuesta. Sí es posible, pues en un proceso adiabático Q = 0, y por el 1r principio Q = W + E -> 0 = W + E -> W = -E, el trabajo realizado será igual al cambio de energía interna cambiado de signo (pues en la expansión se pierde energía interna), además sabemos que en los gases perfectos el cambio de energía interna sólo depende del cambio en la temperatura \triangle E=nc_V\triangle T.


3 - Es conocido el hecho de que, al inflar a mano un neumático con una bomba de aire, ésta se calienta, el efecto es muy notable con los neumáticos de bicicleta de carretera, que han de hincharse a presiones elevadas (de unos 7kg/cm² para un ciclista de 70kg), y no puede explicarse por simple rozamiento, ¿a que es debido? Considerar que el aire se comporta como un gas perfecto.

Respuesta. La bomba comprime el aire y lo inyecta en la cámara del neumático; para hacerlo realizamos trabajo sobre el aire (la fuerza que aplicamos a la bomba), lo que no suministramos es calor: nosotros no calentamos el aire poniéndolo en contacto con una fuente térmica cálida, nos limitamos a realizar un trabajo. Por el 1r principio, Q = W + E, el incremento de energía interna es igual al trabajo suministrado (signo negativo, recordemos el convenio de signos, el trabajo realizado por el sistema es positivo), -W = E, pero en los gases perfectos el cambio de energía interna sólo depende del cambio en la temperatura \triangle E=nc_V\triangle T, así pues, el gas se ha de calentar si aumenta su energía interna.


4. Un ciclo de un gas ideal tiene forma rectangular en el diagrama Presión-Volumen. ¿Qué aspecto tendrá en un diagrama Presión-Temperatura?

Respuesta. En la figura 2 vemos el diagrama PV, indicando el carácter de cada transformación entre los cuatro estados ABCD. Siendo un gas ideal, tenemos la ecuación de estado PV = nRT. Veamos que no dice para cada transformación respecto a las variables P, T:

Fig.2: diagrama PV en forma de rectángulo

A->B: dilatación a presión constante, se produce trabajo, V=\left(\frac{nR}P\right)T, la temperatura depende linealmente del volumen, a más volumen, más temperatura; esto sólo es posible si hay una aportación de calor Q al gas, que se expandirá y calentará a presión constante. P = cte implica que será una linea horizontal en el plano PT.

B->C: Descompresión a volumen constante, P=\left(\frac{nR}V\right)T, dependencia lineal, la temperatura disminuye con la presión, manteniendo el volumen constante, implica que no se produce trabajo, y por el 1r principio la disminución de energía interna (debido al descenso de temperatura) es igual al calor emitido: se debe quitar calor al gas para disminuir la presión a volumen constante. En el diagrama PT será una línea con pendiente positiva \frac{nR}V

C->D: reducción de volumen a presión constante, se recibe trabajo (trabajo negativo), V=\left(\frac{nR}P\right)T, la temperatura depende linealmente del volumen, desciende con éste, y se desprende calor (1r principio, Q = W + E, ambos negativos). En el diagrama PT será una línea horizontal.

D->A: compresión a volumen constante, P=\left(\frac{nR}V\right)T, dependencia lineal, la temperatura aumenta con la presión, manteniendo el volumen constante, implica que no se produce trabajo, luego Q = E, el gas absorbe calor. En el diagrama PT será una línea con la misma pendiente que B->C.

Con todo ello nos resulta un ciclo con forma de paralelogramo:

Fig. 3: diagrama PT equivalente


 

Publicado en Fí­sica, Termodinámica | Etiquetado , , , , , | Deja un comentario

Trabajo, calor, 1r principio de la Termodinámica

La Termo-dinámica se ocupa de la dinámica (movimiento, transferencia) térmica, y utiliza las magnitudes temperatura, volumen, presión, energía, calor y trabajo para relacionarlas entre sí. El intercambio calor-energía-trabajo supondremos que se realiza entre un "sistema" - una parte del Universo delimitada por un volumen cerrado - y su entorno - el resto del Universo - , aunque a menudo restringiremos el entorno a las denominadas "fuentes térmicas" que definiremos. Entonces, hay magnitudes que utilizaremos para definir el estado termodinámico del sistema, y otras para definir su evolución o cambio de estado. Si tenemos alguna ecuación matemática que relacione entre sí las magnitudes que definen el estado, la llamaremos ecuación de estado; un ejemplo conocido de ecuación de estado es la de los gases perfectos, PV = nRT.

Diagramas PV de Clapeyron

Una herramienta útil para entender y trabajar con cambios de estado es el diagrama de presión y volumen, que muestra como un punto cada valor de esas magnitudes, y como una línea su evolución estado a estado.  Por ejemplo, consideremos un gas encerrado en un cilindro con un émbolo encima del cual hemos colocado un conjunto de pequeños pesos, que vamos quitando uno a uno (figura 1).

Fig. 1: un gas comprimido en un pistón con una presión P que variamos quitando unos pequeños pesos

El  gas, cuando no tocamos los pesos, y si lo hemos dejado en reposo, estará en equilibrio termodinámico: no hay flujo de calor y la presión y volumen están estables; tendrá una presión P y un volumen V representados por el punto A del diagrama PV (figura 2).

Fig.2: diagrama PV de Clapeyron, mostrando cambios irreversibles AB, BA', A'B', ...

Entonces levantamos uno de los pesos y la presión disminuye bruscamente: el sistema está ahora en el punto B, que no es un estado de equilibrio (no es estable), entonces el gas se expande espontáneamente hasta alcanzar el punto A' que vuelve a ser de equilibrio.  Vamos repitiendo el procedimiento con los otros pesos, y el gas va pasando por los estados  A",  A''', ...etc hasta llegar al estado final C. La línea continua azul que une todos los estados de equilibrio A, A', ..., C representa una transformación suave desde el estado inicial A hasta el final C sin que el sistema deje de estar en equilibrio.

En general, un cambio brusco en una magnitud saca del equilibrio al sistema: decimos que es un cambio de estado irreversible; sólo si los cambios son tan progresivos y suaves que, en cada momento, el sistema sólo se desvía del equilibrio una cantidad infinitesimal, de forma que podemos considerar que el sistema evoluciona estando siempre en equilibrio, entonces diremos que es un cambio de estado reversible. Los cambios que se ven como una función escalonada en el diagrama de Clapeyron son todos irreversibles, mientras que la línea azul continua que une estados de equilibrio A y C representa una transformación reversible.

Para estar siempre en equilibrio y no obstante cambiar de estado, las variaciones han de ser continuas: imaginemos que en lugar de los pesos de la figura 1 usamos arena, y que vamos quitando la arena grano a grano, entonces tendremos una buena aproximación a una transformación reversible; se llama así porque en cualquier estado intermedio entre A y C podríamos volver a poner los granos de arena y el sistema recorrería el camino inverso de C hacia A de forma suave: se puede revertir con una modificación infinitesimal. En cambio si realizamos una transformación brusca irreversible como la A->B, si queremos "deshacerla" no basta con una modificación infinitesimal de la presión, pues la presión exterior ha variado bruscamente y una modificación infinitesimal de la presión no detendrá la expansión del gas que seguirá evolucionando hacia el estado B. En definitiva:

Una transformación reversible mantiene al sistema siempre en equilibrio, y puede revertirse cambiando infinitesimalmente las condiciones. Una transformación irreversible aleja al sistema del equilibrio y no  puede revertirse cambiando infinitesimalmente las condiciones.

El estudio de las tranformaciones irreversibles, la termodinámica de los procesos irreversibles, es considerablemente más complicada que la de los procesos reversibles, a las cuales dedicaremos el resto de este artículo.

Cambios isobáricos, isotérmicos, isócoros

En el ejemplo del gas en un cilindro de la figura 1, en el cambio brusco de estado AB varia la presión pero no el volumen, que aumenta: decimos que es una transformación isócora. Cuando el gas se expande en el cambio BA', lo hace a presión constante: es una transformación isobárica.

Fig. 3: sistema de la figura 1 envuelta en un baño térmico a temperatura constante

Si mantenemos la temperatura constante, por ejemplo fabricando el cilindro de un material muy conductor del calor y rodeándolo de un baño térmico (un entorno que no varia de temperatura por su gran capacidad calorífica, por ejemplo, una habitación acondicionada o una piscina, figura 3) entonces es una transformación isotérmica; supongamos que el sistema es un gas perfecto, entonces se cumplirá, para una transformación isotérmica, que PV = cte, y la representación gráfica de la curva de estados de equilibrio del gas ideal en el diagrama PV será una hipérbola con asíntotas en los ejes, tal como la línea azul de la figura 2.

Trabajo realizado por dilatación

Supongamos que tenemos un sistema delimitado por una superficie S, sujeto a una presión exterior P, que se está dilatando. El trabajo W realpizado al desplazarse una distancia x contra una fuerza F es W = S·x, y la fuerza total ejercida por la presión F sobre la superfície S es F = P·S. Pensemos en una dilatación infinitesimal del sistema, de forma que sus límites se muevan una distancia dx y su superfície aumente en dS; el trabajo elemental realizado contra la presión exterior será dW = dF·dx = P·dS·dx, pero el producto dS·dx es igual al incremento de volumen dV = dS·dx. luego dW = P·dV, integrando entre el volumen inicial y el volumen final nos queda

W=\int_{V_1}^{V_2}P\operatorname dV [1]

En el diagrama PV, esta integral coincide con el área bajo la curva de cambio de estado (figura 4)

Fig. 4: El trabajo realizado por el sistema al expandirse de A a B es el área rayada

Si en vez de una expansión tenemos una contracción, el desplazamiento cambia de signo, y el trabajo resultante es negativo: es trabajo realizado sobre el sistema. Tenemos pues: si el sistema realiza un trabajo W, éste tendrá un valor positivo, y si se realiza un trabajo sobre el sistema, tendrá un valor negativo.

Transformaciones cíclicas

Fig.5: una transformación termodinámica PV cíclica, o brevemente, un ciclo

Son transformaciones que empiezan y acaban en el mismo estado. En la figura 5 vemos representado un ciclo en el diagrama PV recorrido en el sentido de las agujas del reloj; al recorrer la rama superior el trabajo de dilatación realizado será igual al área bajo la curva roja, en cambio al recorrer la rama inferior el sistema se contrae, luego el trabajo es negativo, y será igual al área bajo la curva azul cambiada de signo. El trabajo total del ciclo será la suma de los dos anteriores, y se ve claramente que será igual al área contenida dentro del ciclo y tendrá un valor positivo (trabajo neto realizado por el sistema).

Si recorremos el ciclo en sentido inverso a las agujas de reloj, los signos de los trabajos cambian, y el trabajo del ciclo será negativo, el ciclo recibirá un trabajo neto.

Primer principio de la Termodinámica

El principio de conservación de la energía enuncia que ésta nunca se pierde sino que se transforma de unos tipos a otros; en termodinámica hemos visto que los tipos de energías que estudiamos son la energía térmica o calor Q, y el resto de energías que puedan tener un papel en las transformaciones las agrupamos todas bajo la denominación genérica "energía E", con mención especial a la energía genérica interna del sistema que estudiamos, que denominaremos energía interna del sistema. La energía se puede invertir en realizar un trabajo, y aquí es dónde llegamos a enunciar el primer principio de la termodinámica, que relaciona calor y trabajo:

1r principio: Si un sistema absorbe una cantidad de calor Q i una cantidad de trabajo W, y emite una cantidad de calor Q' y realiza una cantidad de trabajo W, se cumple siempre que:

Q + Q' = W + W' [2]

Fig. 6: balance entre calor absorbido y emitido, y trabajo absorbido y realizado

El 1r principio no es más que un balance energético donde usamos el trabajo para explicar lo que sucede con el calor que "desaparece"; observemos los convenios de signo: el trabajo realizado por el sistema y el calor absorbido por él se consideran de signo positivo, mientras que el calor emitido por el sistema y el trabajo recibido sobre él son negativos. Pero no significa que el trabajo sea una forma de energía; parecería pues que el 1r principio violaría la conservación de la energía, veamos que esto no es así.

Energía interna y magnitudes conservativas

Consideremos la transformación A->B de la figura 5, que es parte de un clclo, una expansión del sistema que generará un trabajo W > 0 dado por la expresión 1. Supongamos que el sistema absorbe un calor Q > 0  en esa transformación; la diferencia entre calor absorbido y trabajo realizado, Q - W, debe de "ir a algún sitio"; diremos que se utiliza para aumentar la energía interna del sistema:

\triangle E_{A\rightarrow B}=Q-W

Cerrando el ciclo, la transformación B->A es una compresión, luego se genera trabajo sobre el sistema, W' < 0, y quizá el sistema libere calor Q' < 0, de nuevo, la diferencia Q' - W' se empleará en variar la energía interna del sistema:

\triangle E_{B\rightarrow A}=Q'-W'

Si consideramos ahora todo el ciclo A->B->A, como el sistema queda en el mismo estado inicial, la energía interna también ha de ser la misma; sumando las dos igualdades anteriores:

\begin{array}{l}\begin{array}{l}\triangle E_{AB}=Q-W\\\triangle E_{BA}=Q'-W'\end{array}\\---------\\\triangle E_{ABA}=\left(Q+Q'\right)-\left(W-W'\right)\end{array}

Luego

\left(Q+Q'\right)-\left(W-W'\right)=0\Leftrightarrow\boxed{Q+Q'=W-W'}

recuperamos el 1r principio; vemos que nos faltaba para entender la conservación de la energía la energía interna del sistema: la energía térmica, el trabajo, y la energía interna se van equilibrando en cada transformación termodinámica. Enunciamos ahora unas propiedades importante:

P1) La energía interna E del sistema es una función del estado del sistema: a cada estado le corresponde una única energía definida.

Por ejemplo si un sistema queda totalmente definido por las magnitudes presión P y Volumen V, entonces la energía interna es función de P, V: E = f(P, V).

P2) El trabajo W no es una función de estado del sistema, para cada estado no es posible asignar un único valor de W.

En efecto, sabemos que el valor del trabajo es igual al área definida por la curva de cambio de estado PV y el eje V; en la figura 7 vemos tres trayectorias distintas entre los mismos estados A, B, y cada una representa un trabajo distinto. Este hecho también se enuncia diciendo que el trabajo no es una magnitud conservativa, significando que si hacemos un ciclo ABA, el trabajo total no tiene porque ser nulo, pues dependerá de las trayectorias seguidas. En cambio la energía interna es una magnitud conservativa, pues siempre la recuperaremos en un ciclo (variación total nula) independientemente del camino seguido.

Fig. 7: El trabajo ralizado en una tranasformación no depende sólo de los estados inicial A y final B

Con esto podemos reformular el 1r principio en una forma equivalente:

1r principio: El calor absorbido por un sistema se invierte íntegramente en el trabajo realizado por él más su incremento de energía interna, Q=W+\triangle E. [3]

Energía interna y temperatura

Diferenciando la expresión 3, \operatorname dQ=\operatorname dW+\operatorname dE, y recordando la expresión 1 para el trabajo: \operatorname dQ=P\operatorname dV+\operatorname dE. [4]

Por otra parte, el calor absorbido y el incremento de temperatura están relacionados con el coeficiente calorífico c, \triangle Q=nc\triangle T, siendo una medida de la masa del sistema (como por ejemplo el número de moles); sucede que c toma valores distintos en transformaciones isóbaras que en isócoras:

a V = cte, \triangle Q=nc_V\triangle T,

a P = cte, \triangle Q=nc_P\triangle T.

Tomemos un proceso a V = cte y diferenciemos: \operatorname dQ=nc_V\operatorname dT. Comparando con 4:

P\operatorname dV+\operatorname dE=nc_V\operatorname dT\Leftrightarrow\operatorname dE=nc_V\operatorname dT

pues V es constante y dV = 0. integrando:

\triangle E=nc_V\triangle T [5]

que nos da la variación de energía interna respecto a la variación de temperatura en un proceso a volumen constante. Pero la energía es una función de estado, y su variación sólo depende de los estados, no de los procesos, por ello, la expresión 5 es válida en general, y nos da la variación de la energía interna en función del cambio de temperatura para cualquier transformación.

Transformaciones adiabáticas; aplicación a los gases perfectos.

Si no hay intercambio de calor, el 1r principio dice que Q=0\Rightarrow W+\triangle E=0. Hemos visto que podemos relacionar el incremento de energía interna en un sistema termodinámico con un incremento de temperatura (ecuación 5):

\operatorname dW+nc_v\operatorname dT=0\Leftrightarrow P\operatorname dV+nc_v\operatorname dT=0 [6]

ecuación que nos da la relación entre diferenciales de volumen y temperatura cuando Q = 0.

Aplicación a los gases perfectos

La ecuación de estado es PV = nRT, aislando P y sustituyendo en 6:

\frac{nRT}V\operatorname dV+nc_v\operatorname dT=0

dividimos todo por nCvT e integramos:

\int\frac R{c_v}\frac{\operatorname dV}V+\int\frac{\operatorname dT}T=0\Leftrightarrow\frac R{c_v}\ln\left(V\right)+\ln\left(T\right)=\mathrm C

donde C es una constante de integracion; operamos con los logaritmos, y además definimos la constante \gamma=c_p/c_v, y usamos la relación de Mayer, c_p-c_v=R para obtener:

\left.\begin{array}{r}\frac R{c_v}\ln\left(V\right)+\ln\left(T\right)=\mathrm C\\c_p-c_v=R;\;\gamma=\frac{c_p}{c_v}\end{array}\right\}\Rightarrow\frac{c_p-c_v}{c_v}\ln\left(V\right)+\ln\left(T\right)=\mathrm C\Rightarrow\ln\left(\mathrm{TV}^{\mathrm\gamma-1}\right)=\mathrm C

Volvemos a utilizar la ecuación de estado ideal PV = nRT para eliminar T de la expresión anterior:

\ln\left(\frac{PV}{nR}V^{\mathrm\gamma-1}\right)=\mathrm C\Leftrightarrow\ln\left(\frac1{\mathrm{nR}}\mathrm{PV}^\mathrm\gamma\right)=\mathrm C

Como 1/nR es una constante, lo incorporamos a la constante de integración, y concluimos:

\ln\left(\frac1{\mathrm{nR}}\mathrm{PV}^\mathrm\gamma\right)=\mathrm C\Rightarrow\ln\left(\mathrm{PV}^\mathrm\gamma\right)=\mathrm{cte}\Rightarrow\boxed{\mathrm{PV}^\mathrm\gamma=\mathrm{cte}}

que nos da la ecuación de los estados de las transformaciones adiabáticas de gases perfectos;  sabiendo que c_p-c_v=R>0, deducimos que c_p=R+c_v>c_v y por tanto que \gamma>1.

La figura 8 compara gráficamente las transformaciones adiabáticas con las isotérmicas de gases perfectos, en las que PV = cte:

Fig.8: transformaciones isoterma y adiabáticas de gases perfectos en el diagrama PV

Observamos que la isotérmica tiene una pendiente menor que todas las adiabáticas, pues \gamma>1. Asintóticamente todas coinciden en los ejes P y V.

Máquina térmica de Carnot para gases perfectos

Pensemos ahora en un ciclo de transformaciones reversibles como el mostrado en la figura 9, recorido en el sentido de las agujas del reloj, 1-2-3-4, y con las siguientes condiciones:

  • 1-2: isotérmica, T=cte, absorbe Q > 0
  • 2-3: adiabática, Q = 0, T variable pasa a ser T'
  • 3-4: isotérmica, T'=cte, emite Q' < 0
  • 4-1: adiabática, Q = 0, T'=cte

El diagrama tiene la geometría correcta, pues hemos visto que las curvas adiabáticas tienen mayor pendiente que las isotérmicas, al menos lo hemos visto para gases perfectos, y nos restringimos a este caso en lo que sigue.

En las isotérmicas la energía interna debe de ser constante según la ecuación 5, y en este caso por el 1r principio, el trabajo ha de ser igual al calor, Q = W(12) y Q'=W(34). Además, el valor del área 112'1' coincide con el valor de Q, y el del área 44'33' coincide con el valor de Q',

En cambio en las adiabáticas hay cambios de temperatura y de energía interna, y por el 1r principio W(23) = incremento E(23) y W(41) = incremento E(41). Además, los valores de las áreas 11'44' y 22'3'3 coinciden con los cambios de energía interna. Como el cambio de energía interna del ciclo completo ha de ser cero, se deduce que las áreas 11'44' y 22'3'3 han de ser iguales (en el diagrama no lo son).

Las variaciones de temperatura, en las transformaciones adiabáticas, viene dadas por (5) y son:

\triangle T=\frac{nc_v}{\triangle E}=\frac{nc_v}{-W}\left\{\begin{array}{l}<0\text{ en expansión adiabática}\\>0\text{ en compresión adiabática}\end{array}\right.

Como en la adiabática 23 la variación de temperatura es negativa, implica que T > T': el ciclo de Carnot opera entre dos fuentes térmicas a temperaturas distintas, absorbiendo calor Q de la fuente más caliente y cediendo calor Q' a la fuente más fría, y realizando un trabajo neto W = Q - Q'.

Ciclo inverso de Carnot: bomba de calor

Siendo un ciclo reversible puede invertirse en cualquier momento y hacerse en sentido contrario al del reloj; en el sentido 4321 la 32 es una compresión adiabática, que implica un aumento de temperatura, luego T > T' como en el ciclo directo; en la 21 tenemos compresión isoterma que libera calor Q, y en cambio en 43 se absorbe: el ciclo inverso absorbe calor de la fuente fría, cede calor a la fuente caliente y consume un trabajo W, por ello decimos que bombea calor de la fuente fría a la cálida, es decir, que la refrigera.

Rendimiento de un ciclo

Nos preguntamos, ¿que trabajo neto W obtenemos a cambio de inyectar un calor neto Q + Q' en el ciclo? (recordemos que Q' tendrá signo negativo). Este último calor suele provenir de la combustión (fuente térmica a alta temperatura), se genera un trabajo y se libera un calor residual al medio ambiente Q'.  Definimos el índice de rendimiento del ciclo r = W / Q, que idealmente seria 1 (todo el calor se convierte en trabajo). Vemos que r = (Q  + Q') / Q. Se puede demostrar que, para un ciclo de Carnot, se cumple

r=\frac{T-T'}T=1-\frac{T'}T [6]

que sólo depende de las temperaturas de los focos térmicos, y será siempre menor que 1 excepto en el caso extremo de T' = 0 (medio exterior en el cero absoluto de temperaturas). Se puede definir un índice análogo para bombas de calor.

 

 

Publicado en Fí­sica, Termodinámica | Etiquetado , , , , , , , , , , , , , | 2 comentarios

Física del sonido

En anteriores artículos de la categoría "Acústica" hemos sobre todo realizado una descripción de las características del sonido y de la música, un tipo especial de sonido. En este  artículo tratamos los detalles desde un punto de vista más matemático.

Intensidad del sonido

El sonido es una onda elástica de presión, la intensidad que percibimos en el sonido es proporcional al cuadrado de de la amplitud de la vibración (esto es así para el movimiento ondulatorio en general). También podemos expresar la intensidad en función de la variación máxima de presión del aire originada por el sonido al propagarse, como vamos a ver a continuación.

La presión del aire, cuando propaga el sonido, es variable en cada punto, y por tanto el volumen que ocupa una cantidad fija de aire también varia con la presión.

Fig. 1: Onda longitudinal de presión: desplazamientos instantáneos

Imaginemos (figura 1) un volumen de aire dentro de un pequeño cilindro de área S y longitud dx, situado inicialmente entre las posiciones x, x + dx, con volumen V = S·dx,  que es atravesado por la onda de expansión (el sonido) en la dirección de el eje X; decimos que sea pequeño para aproximar el frente de onda a una onda plana, cuando en realidad es esférica. El cilindro se desplazará y ensanchará sólo en la dirección X (la base S es transversal a X, y la onda es longitudinal según X) de forma que las posiciones de los extremos serán (x + dy) y (x + dx + dy), donde llamamos y a la elongación de la onda, que recordemos es longitudinal y por tanto paralela al eje X.

El el nuevo volumen será S(dx + dy), luego el incremento de volumen es dV = S·dy. Luego:

\frac{\operatorname dV}V=\frac{S\operatorname dy}{Sdx}=\frac{\operatorname dy}{\operatorname dx} [1]

No confundamos la variable y con la dimensión vertical del eje de ordenadas: aquí y representa la elongación de la onda longitudinal en la dirección horizontal.

Por otro lado, la velocidad de propagación de una onda longitudinal en un medio elástico sólo depende de su densidad \rho y del denominado módulo de compresibilidad \chi=-V\frac{\triangle p}{\triangle V} [2], donde p es la presión, según la igualdad:

v=\sqrt{\frac\chi\rho} [3]

Combinando las expresiones [1] y [2], e identificando las variaciones finitas \triangle con las diferenciales:

\chi=-\cancel{\operatorname dV}\frac{\operatorname dx}{\operatorname dy}\frac{\operatorname dp}{\cancel{\operatorname dV}}=-\operatorname dp\frac{\operatorname dx}{\operatorname dy}

Sustituimos en la expresión de la velocidad [3]:

v=\sqrt{-\operatorname dp\frac{\operatorname dx}{\rho\operatorname dy}}\Rightarrow-\operatorname dp=v^2\rho\frac{\operatorname dy}{\operatorname dx} [4]

Supongamos ahora que el sonido es un movimiento ondulatorio armónico (por tanto, es un sonido musical simple, puro) con elongación dada por

y=A\sin\omega\left(t-\frac xv\right) [5]

Derivamos respecto x manteniendo t constante:

\frac{\operatorname dy}{\operatorname dx}=-\frac Av\omega\cos\left(t-\frac xv\right)

Sustituimos el valor de dy/dx en [4]:

\operatorname dp=Av\rho\omega\cos\omega\left(t-\frac xv\right)=\widehat A\cos\omega\left(t-\frac xv\right) [6]

que nos dice que la variación de presión debida al sonido se propaga como una onda de presión, con amplitud \widehat A y que, comparada con la onda de elongación [5], está adelantada en fase 90⁰ (ya que depende de la función coseno en vez de la seno).

El valor máximo de la sobrepresión es \operatorname dp_{MAX}=Av\rho\omega; en general, la intensidad de una onda (energía que atraviesa la unidad de superficie normalmente a ella en la unidad de tiempo) es I=\frac12A^2\rho\omega^2v, de donde deducimos, usando [6]:

\frac{\operatorname d{p^2}_{MAX}}I=\frac{A^2v^2\rho^2\omega^2}{\frac12A^2\rho\omega^2v}=2v\rho\Rightarrow\boxed{I=\frac{\operatorname d{p^2}_{MAX}}{2v\rho}}

La intensidad del sonido es proporcional al cuadrado de la sobrepresión  máxima asociada a la onda de presión.

Cálculo de la velocidad del sonido en el medio ambiente

En el aire y en condiciones normales, las variaciones de presión son muy rápidas, y se puede considerar que no da tiempo al intercambio de calor (un gas al ser comprimido se calienta y al ser expandido se enfría), en términos termodinámicos, son variaciones de presión adiabáticas. Supongamos que el aire se comporta aproximadamente como un gas perfecto, cumpliendo la ecuación de los gases ideales en expansión adiabática, pV^\gamma=cte [7], donde \gamma es el coeficiente adiabático. Diferenciando: \operatorname dp\cdot V^\gamma+\gamma V^{\gamma-1}\cdot\operatorname dV=0\Rightarrow\operatorname dp\cdot V+\gamma\operatorname dV=0, de donde \operatorname dp=-\gamma p\frac{\operatorname dV}V. Sustituyendo en la expresión [2]:

\operatorname dp=-\gamma p\frac{\operatorname dV}V,\;\chi=-V\frac{-\gamma p\frac{\operatorname dV}V}{\operatorname dV}=\gamma p

Sustituyendo en la velocidad [3]:

v=\sqrt{\frac\chi\rho}=\sqrt{\frac{\gamma p}\rho} [8]

Para el aire es \gamma=1.4; además, a 0⁰ grados Celsius y presión p = 1atm = 1.013·10⁵ N/m², la densidad del aire es ρ = 1.293  kg/m³, sustituyendo en [8] obtenemos

v=\sqrt{\frac{1.4\cdot1.013\cdot10^5}{1.293}}=331

Vemos que con consideraciones teóricas obtenemos un valor muy de acuerdo con las mediciones experimentales. Nos damos cuenta, además, de que la velocidad del sonido en el aire dependerá de la temperatura, pues con ella variará la presión atmosférica. Usando la ecuación de los gases perfectos pV = nRT, para eliminar la presión p de la ecuación [8] se obtiene (T en grados Kelvin):

v=\sqrt{\frac{\gamma RT}M} [9]

En la figura 2 vemos la variación de v según la temperatura obtenida aplicando [9] para un rango de [0, 30] grados centígrados, es prácticamente lineal, variando unos 0,6 m/s por cada grado.

Fig. 2: Variación de la velocidad del sonido con la temperatura

Efecto Doppler

Cuando oímos un tren que se acerca con el silbato actuando, nos parece que el sonido es más agudo que si el tren se aleja de nosotros. El mismo efecto se observó con la luz, su frecuencia parece más alta cuando el foco emisor se acerca a nosotros que cuando se aleja, si la luz es visible, en el primer caso vemos la luz "más azul" (corrimiento al azul) y en el segundo, "más roja" (corrimiento al rojo). Veamos la expresión exacta que nos da el corrimiento en frecuencia.

Fig. 3: foco de sonido y receptor en movimiento mutuo. Efecto Doppler.

En la figura 3 representamos la situación: un observador O se mueve hacia la derecha con velocidad v y un foco de luz se mueve en dirección contraria con velocidad u (todo con respecto a un sistema de referencia que suponemos fijo); los frentes de onda se suponen planos (realmente son esféricos). La velocidad del sonido relativa al foco es w, y la absoluta (respecto al suelo) será u + w, mientras que la velocidad relativa al observador O será u + w - v (estamos suponiendo que las velocidades son mucho menores que la de la luz, y por ellos despreciando los efectos relativistas).

La frecuencia percibida del sonido es el número de ondas que llegan al receptor por segundo; en el caso de observador y foco en reposo, esa frecuencia será ν, y está relacionada con la velocidad w de la onda y su longitud λ por la relación w = λ·ν, de donde ν = w/ λ.  Cuando el foco de sonido se mueve, el número de ondas por segundo que  llegan a un receptor no será el mismo, y se percibirá una frecuencia ν' distinta. La onda de sonido en sí no cambia: su longitud de onda, su amplitud, su forma, no cambian, sólo que su velocidad relativa al observador sí lo hace (figura 4).

Fig. 4: onda sinusoidal, se desplaza a velocidad w relativa a una referencia fija

Desde el punto de vista del foco, que se mueve en la misma dirección que los frentes de onda y con velocidad u, se emiten ν ondas por segundo, y se "ven" los frentes de onda alejándose del foco a velocidad (w - u), por ello el foco verá una longitud de onda de λ = (w - u) / ν.

Veamos ahora el punto de vista del observador O, que se mueve a velocidad v acercándose al foco, y por tanto verá moverse a los frentes de onda a velocidad (v + w). Razonando de la misma forma que para el foco, verá una longitud de onda λ = (v + w ) / ν'. Igualando las dos expresiones para λ:

\lambda=\frac{w-u}\nu=\frac{w+v}{\nu'}

que nos proporciona la relación buscada entre frecuencia emitida f y frecuencia observada f'. Por ejemplo, si el silbato de un tren emite un sonido de frecuencia f = 440 Hz (440 ondas por segundo), estamos parados respecto al suelo (luego v = 0) y el tren se acerca hacia nosotros a velocidad u = 30 m/s, la frecuencia percibida será:

\frac{331-30}{440}=\frac{331+0}{f'}\Rightarrow f'=484 m/s

donde hemos tomado para la velocidad del sonido w = 331 m/s. Cuando el tren llega a nuestra posición para empezar a alejarse de nosotros, tendremos una velocidad  u = -30 m/s, luego:

\frac{331-(-30)}{440}=\frac{331+0}{f'}\Rightarrow f'=403 m/s

Vemos que se aprecian claramente variaciones de frecuencia de aproximadamente un 9% del valor real, en más o en menos.

Acústica de salas

Para terminar con este breve paseo por la física del sonido, vemos la física del acondicionamiento de salas para escuchar música.

Cuando en una sala un emisor emite un sonido breve (una nota musical o una sílaba si es un orador), la onda sonora se expande por la sala, reflejándose en las paredes, suelo y techo. A un oyente que esté en el otro extremo le llegará el sonido primero por el camino directo más corto, pero un instante más tarde recibirá las ondas reflejadas, con menor intensidad. Las reflejadas siguen expandiéndose por la sala y se vuelven a reflejar una y otra vez, llegando cada vez más atenuadas al oyente, hasta que son demasiado débiles para ser percibidas. Decimos que hay reverberación en la sala para resumir este fenómeno de estar oyendo el sonido directo y reflejado. La reverberación refuerza el sonido: si salimos al aire libre se oirá menos. Tiene también un grave inconveniente: si se emiten más notas o sílabas, puede pasar que se mezclen en nuestro oído varias notas: las que llegan de forma directa y las anteriores que todavía están reverberando, las notas "se pisan" unas a otras.

En las iglesias antiguas se tenia en cuenta esto: la reverberación es grande allí, ayudando a que se oiga bien en el otro extremo, pero obliga a hablar lentamente, para dar tiempo a que se disipen las palabras anteriores. El tiempo de reverberación de una sala se define como el tiempo que tarda en hacerse inaudible el sonido reflejado. más exactamente, se define como es tiempo después del cual la intensidad del sonido es la millonésima parte. Se intuye que el tiempo de reverberación óptimo depende de la finalidad de la sala: para conferencias estará cerca de un segundo, para un concierto de rock será inferior y para un concierto de órgano será mayor.

Para ajustar la reverberación de una sala podemos usar materiales más o menos absorbentes del sonido en las paredes; en cada reflexión se pierde una parte de la energía: el coeficiente de absorción sonora α es la fracción de energía absorbida: si I es la energía incidente,  I' es la reflejada, entonces se cumple I' = I·(1-α). Para α cercano a 1, el sonido es absorbido totalmente por la superficie. En la página del portal de Acústica y Sonido encontraremos más detalles del coeficiente α y una tabla de valores para distintos materiales. Sucede que α no es constante para cada material, sino que depende de la frecuencia del sonido; algunos valores típicos son 0.02 para paredes de hormigón, 0.08 para el yeso, o 0.80 para el poliuretano, a una frecuencia de 1 KHz. Para ajustar la reverberación podemos poner o quitar superficies con α cercano a 1: alfombras, cortinas gruesas, etc. Incluso la ropa que llevan puesta los espectadores influye en la reverberación: una sala repleta de gente en verano con poca ropa reverbera más que en invierno con más ropa.

Puede demostrarse, haciendo un estudio estadístico de todas las múltiples reflexiones en las diversas superficies, que la intensidad del sonido decrece exponencialmente con el tiempo, y que el tiempo de reverberación de la sala viene dado aproximadamente por la fórmula de Sabine:

T=0.16\frac V{S\overline\alpha}

donde S es la superficie total de las paredes, V el volumen de la sala, y \overline\alpha es el coeficiente de absorción medio de la sala, calculado promediando todos las superficies de la sala:

\overline\alpha=\frac{{\overline\alpha}_1S_1+{\overline\alpha}_2S_2+\dots}{S_1+S_1+\dots}

Se puede tener en cuenta también la absorción del sonido por el aire: en salas grandes el recorrido del sonido entre reflexiones será también mayor, y el efecto será más notable, especialmente a frecuencias altas. Se puede demostrar que el recorrido libre medio L, que se define como la media de los tramos recorridos entre reflexiones sucesivas, vale

L=\frac{4V}S

Podemos pues complementar el coeficiente de absorción α con un término adicional que tenga en cuenta la absorción por el aire, que a su vez depende del recorrido libre medio:

\alpha'=\overline\alpha+kL

donde k es una constante que depende de las condiciones físicas del aire (temperatura, humedad, presión ...).

Ejemplo práctico

Una habitación destinada a audición musical tiene un tiempo de reverberación de 1s. Se pone un tabique hecho del mismo material que las paredes, de forma que la superficie total S de la habitación aumenta un 20%, evidentemente el volumen V no varia. Calcular el nuevo tiempo de reverberación.

Antes de poner el tabique teníamos, por la formula de Sabine,

T=0.16\frac V{S\overline\alpha}=1

Con el tabique será

T'=0.16\frac V{1.2S\overline\alpha}

dividiendo la 1a por la 2a:

\frac T{T'}=1.2\Rightarrow T'=\frac T{1.2}=\frac{\displaystyle1}{\displaystyle1.2}=\frac56\approx0.8s

Publicado en Acústica, Fí­sica | Etiquetado , , , , , , , , | Deja un comentario

Momentos de vectores

Introducción

En el artículo Vectores en Física se habló de algunas propiedades geométricas de los vectores (la invariancia respecto de transformaciones de coordenadas) que son importantes para representar magnitudes físicas como la fuerza o la velocidad angular. En este artículo vemos el concepto teórico de momento de un vector respecto a un punto o a una recta, que físicamente tiene importancia para calcular el efecto que un vector ejerce respecto a ese punto o recta; un ejemplo práctico es la antigua ley de la palanca de Arquímedes, que nos muestra cómo varía el efecto de una fuerza con el punto de aplicación: la fuerza P, más alejada del punto de apoyo, ejerce una acción de giro igual a la fuerza mayor R que está más cerca; esa "acción de giro" se materializa usando el concepto de momento de cada fuerza respecto al punto de apoyo.

Fig. 1: Ley de la palanca (By Dnu72 (Own work) [GFDL (http://www.gnu.org/copyleft/fdl.html) or CC BY-SA 4.0-3.0-2.5-2.0-1.0 (https://creativecommons.org/licenses/by-sa/4.0-3.0-2.5-2.0-1.0)], via Wikimedia Commons)

Otro ejemplo lo vemos en la siguiente figura, en la que se representa la sección de un rodillo situado en un plano horizontal, al que hemos atado una cuerda r también horizontal en su parte superior.

Fig. 2: un cilindro sujeto a una fuerza tangente efectuada tirando de una cuerda, el efecto de giro no depende del punto exacto de aplicación de la fuerza

Si ahora aplicamos una fuerza a lo largo de la cuerda, digamos F_1 o F_2, el efecto será que el rodillo adquirirá una velocidad angular w, que no dependerá más de la magnitud de la fuerza: fuerzas de magnitud igual producirán la misma velocidad angular, independientemente del punto de aplicación de la fuerza (círculos en azul). Físicamente, diremos que la rotación del cilindro es causada por el momento de la fuerza aplicada, y ese momento depende de la magnitud y de la dirección de la fuerza, pero no de su posición a lo largo de la recta r.

Cuando en una situación física encontramos vectores que se comportan de este modo, causando el mismo efecto independientemente de si movemos su punto de aplicación a lo largo de una recta, decimos que los vectores son deslizantes. Tiene pues sentido estudiar los momentos ya no de fuerzas, sino de vectores deslizantes en general.

Momento de un vector respecto a un punto

Si un vector v tiene su origen (o está aplicado en) en punto P, el momento m de v con respecto a otro punto O se define por:

\boldsymbol m=\left(P-O\right)\times\boldsymbol v [1]

donde la expresión (P-O) simboliza la diferencia de las coordenadas de los dos puntos, y el producto \times representa el producto vectorial.

Fig. 3: momento m del vector v, aplicado en P, respecto del punto O

Recordemos que (P-O) puede verse como el vector OP con orígen en O y extremo en P; además, recordando las propiedades del producto vectorial, el vector momento m será perpendicular al plano formado por los vectores v y OP, y será un pseudovector, o vector polar (ver  Vectores en Física), su módulo valdrá el doble del área del triángulo formado por los vectores v y OP, o analíticamente,

\left|\mathbf m\right|=\left|\left(\mathbf P\boldsymbol-\mathbf O\right)\boldsymbol\times\mathbf v\right|=\left|\left(\mathbf P\boldsymbol-\mathbf O\right)\right|\cdot\left|\boldsymbol v\right|\cdot\sin\left(\alpha\right)

siendo \alpha el ángulo formado por los vectores v y OP.

Ejemplo 1: en el caso de la palanca (figura 1) llamamos O al punto de apoyo, y los puntos P, R de aplicación de las fuerzas y R estan situados a una distancia Bp y Br respectivamente, además, los vectores OP y OR son perpendiculares a las fuerzas y R y todos estos vectores están en un mismo plano que llamamos XY (son coplanarios). En estas condiciones, si aplicamos la definición para calcular el moment total respecto a O:

\begin{array}{l}\boldsymbol m(\boldsymbol P)=\boldsymbol O\boldsymbol P\times\boldsymbol P=\begin{vmatrix}\widehat i&\widehat j&\widehat k\\B_P&0&0\\0&-P&0\end{vmatrix}=-\widehat kPB_P;\\\boldsymbol m(R)=\boldsymbol O\boldsymbol R\times\boldsymbol R=\begin{vmatrix}\widehat i&\widehat j&\widehat k\\-B_R&0&0\\0&-R&0\end{vmatrix}=\widehat kRB_R\end{array}

Luego la suma de momentos es -\widehat k\left(RB_R-PB_P\right), un vector en la dirección del versor \widehat k, que "sale" de la pantalla. En la condición de equilibrio el momento debe de ser cero, para ello debe de cumplirse que PB_P=RB_R, expresión que coincide con la ley de la palanca.

Propiedades del vector momento

Una primera propiedad es que el vector momento es un pseudovector (o vector axial) debido a que se obtiene del producto vectorial de dos vectores polares (ver por ejemplo Vectores en Física).

Si el vector v lo desplazamos a lo largo de su recta soporte r hasta aplicarlo en otro punto P' (figura 2), geométricamente el nuevo triángulo formado por OP' - v tendrá la misma área que el OP - v , pues la base v es la misma y la altura h (obtenida trazando la perpendicular a la recta r pasando por O) de los dos triángulos es la misma. Además, el plano que contiene al vector v y a OP es el mismo que contiene a v y a OP' (es el plano definido por el punto O y la recta r), por tanto vemos que:

Propiedad 1: El momento m de un vector v respecto a un punto O no depende del punto de aplicación del vector, siempre que esté sobre la recta r que contiene al vector.

Fig. 4: el momento m de un vector v respecto un punto fijo O no depende del punto de aplicación P, P', etc

Por tanto los vectores en lo que respecta a su momento respecto a un punto fijo O se comportan como vectores deslizantes. En cambio si variamos del punto O al O', entonces si obtenemos un cambio en el vector momento m: podemos verlo planteando P - O' = (P - O) + (O - O') y sustituyendo en la expresión del momento m':

\begin{array}{l}\boldsymbol m\boldsymbol'=\left(P-O'\right)\times\boldsymbol v=\left[\left(P-O\right)+\left(O-O'\right)\right]\times\boldsymbol v\Rightarrow\\\boldsymbol m\boldsymbol'=\left(P-O\right)\times\boldsymbol v+\left(O-O'\right)\times\boldsymbol v\Rightarrow\\\boxed{\mathbf m\boldsymbol'\boldsymbol=\mathbf m\boldsymbol+\left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\times\mathbf v}\;\lbrack2\rbrack\end{array}

Viendo al producto \left(\mathrm O-\mathrm O'\right)\boldsymbol\times\mathbf v como el momento de v estando aplicado en O con respecto a O' (comparar con la definición [1]), podemos expresarlo en palabras:

Propiedad 2: El momento de un vector aplicado en P respecto a un punto O' es igual al momento de ese vector aplicado en P respecto a otro punto O más el momento respecto a O' del vector aplicado en el punto O.

Otra propiedad del vector momento es el denominado teorema de Varignon:

Propiedad 3 (Varignon): El momento total respecto a un punto O cualquiera de un conjunto de vectores concurrentes en un punto P es igual al momento de la suma del conjunto de vectores respecto a ese punto O.

Cálculo vectorial con vectores deslizantes

Las operaciones con vectores las realizamos usando un sistema de coordenadas y los componentes de los vectores respecto a ese sistema; así, el vector PQ con origen en el punto P y extremo en el punto Q, al restar las coordenadas Q - P obtenemos un vector v con origen en el origen de coordenadas, no en el punto P. Si deslizamos el vector PQ a lo largo de su recta soporte, pasando a estar en los puntos P', Q', el nuevo vector P'Q' tendrá las mismas coordenadas Q'-P' coincidentes con v. Dos vectores de la misma magnitud y direcciones paralelas se dice que son equipolentes; así pues, cuando trabajemos con las coordenadas de vectores deslizantes, realmente estaremos trabajando con las coordenadas de vectores equipolentes con origen en el origen de coordenadas, pero los  momentos que calculemos los supondremos aplicados en los puntos dados.

Fig. 5: vector equipolente a un vector deslizante

Ejemplo 1: obtener el momento respecto al punto O(0,0) de los vectores v, w, ambos aplicados en el punto P(1,1), y con los extremos en Q(2,2) y R(2,3) respectivamente.

Fig. 6: momento de dos vectores concurrentes respecto a un punto O

Las coordenadas son en dos dimensiones, pero el momento, calculado según la definición [1] es un vector perpendicular al plano que contiene los vectores v, w y el punto O; "ampliamos" pues nuestra referencia con una tercera coordenada que "saldrá" del plano de la pantalla hacia nuestro rostro (regla de la mano derecha):

O(0, 0, 0), P(1, 1, 0), OP = P - O = (1, 1, 0);

El vector PQ = (2, 2, 0) - (1, 1, 0) = (1, 1, 0) es paralelo a v con el mismo módulo y origen en (0,0, 0) (v y PQ son vectores equipolentes), por tanto el producto \boldsymbol OP\times\boldsymbol v  es el mismo que el producto \boldsymbol OP\timesPQ, que calculamos usando determinantes:

\boldsymbol OP\times \boldsymbol PQ=\begin{vmatrix}i&j&k\\1&1&0\\1&1&0\end{vmatrix}=\left(0,0,0\right)

obviamente el momento es cero, pues OP y v estan sobre la misma recta; procedemos igual con w: tomamos en su lugar el vector PR  = R-P = (0, 1, 0) y obtenemos el producto vectorial \boldsymbol OP \times \boldsymbol PR:

\boldsymbol OP\times \boldsymbol PR=\begin{vmatrix}i&j&k\\1&1&0\\0&1&0\end{vmatrix}=\left(0,0,-1\right)

este vector momento "entra" verticalmente en la pantalla. El momento total de los dos vectores es la suma de sus momentos. Como v, w son concurrentes, también podemos llegar al mismo resultado aplicando la propiedad 3: sumamos los dos vectores y calculamos el momento de la suma: sustituimos v + w por los vectores que parten de O: (1, 1, 0) + (0, 1, 0) = (1, 2, 0), calculamos el momento:

OP\times PQ=\begin{vmatrix}i&j&k\\1&2&0\\0&1&0\end{vmatrix}=\left(0,0,-1\right)

Momento de un vector respecto a un eje

En los apartados que siguen, se estudian las propiedades de los momentos y de los sistemas de momentos respecto a rectas (ejes); los resultados nos serviran para reducir conjuntos de vectores a un conjunto mínimo equivalente a efectos del momento resultante. También nos interesará encontrar los puntos respecto a los cuales el momento resultante de un sistema de vectores resulta ser mínimo, la cual cosa permitirá resolver problemas de equilibrio.

Fig. 7: momentos de un vector v respecto a los puntos situados en un eje E

Imaginemos un eje E, esto es, una recta sobre la que hemos definido un vector unitario u para definir un sentido, y sobre el eje dos puntos distintos O, O'. Nos preguntamos por los momentos m, m' de un vector cualquiera v respecto a esos puntos. Recordemos la igualdad [2]:

\mathbf m\boldsymbol'\boldsymbol=\mathbf m\boldsymbol+\left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\times\mathbf v

Multiplicando escalarmente los dos lados de esta igualdad por el vector unitario u:

\mathbf m\boldsymbol'\boldsymbol\cdot\boldsymbol u\boldsymbol=\mathbf m\boldsymbol\cdot\mathbf u\boldsymbol+\left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\times\mathbf v\boldsymbol\cdot\mathbf u

El vector \left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\times\mathbf v es perpendicular al plano que contiene a u y v, por tanto el producto escalar es cero, nos queda pues

\mathbf m\boldsymbol'\boldsymbol\cdot\boldsymbol u\boldsymbol=\mathbf m\boldsymbol\cdot\mathbf u

Recordando la interpretación geométrica del producto escalar, tenemos que la proyección sobre el eje E del momento del vector v respecto a cualquier punto del eje E es constante; a este valor, que es un escalar, se le llama momento del vector v respecto al eje E.

En particular, el momento de un vector respecto a un eje paralelo al vector será nulo, pues el vector y el eje serán coplanarios, y el momento será perpendicular a ese plano, luego el producto escalar del vector u con el momento será cero, siendo perpendiculares.

Entonces para un vector v cualquiera, podemos descomponerlo en suma de dos vectores, uno paralelo al eje E y otro perpendicular; el primero tendrá momento respecto al eje nulo, por tanto vemos que:

Propiedad 4: el momento respecto al eje E de un vector v será igual al momento respecto al eje E del vector proyección ortogonal de v respecto a un plano perpendicular a E.

Ejemplo 2: obtener el momento respecto a los ejes Y, Z del vector w, aplicado en el punto P(1,1), y con extremo en R(2,3).

El vector u para el eje Y es u(0, 1, 0), y para calcular el momento de w respecto al eje Y nos vale cualquier punto O sobre el eje, por ejemplo, el (0, 0, 0) que ya hemos obtenido en el ejemplo 1, siendo m = (0, 0, -1); el producto escalar es m·u = (0, 0, -1)·(0, 1, 0) = 0, resultado esperado, pues el momento m es perpendicular al plano que contiene el eje Y. Para el momento respecto del eje Z tomamos u(0, 0, 1), entonces m·u = (0, 0, -1)·(0, 0, 1) = -1.

Momento resultante de un sistema de fuerzas

Si tenemos un conjunto de vectores v_1,v_2,...,v_n podemos calcular sus momentos m_1,m_2,..._m_n respecto a un único punto fijo O; si consideramos otro punto O', aplicando [2] a cada vector y sumando:

\begin{array}{l}m'_i=m_i+\left(O-O'\right)\times v_i\Rightarrow\\\sum_{}m'_i=\sum_{}m_i+\left(O-O'\right)\times\sum_{}v_i\end{array}

Llamando M y M' al momento suma, y R al vector suma (llamado resultante del sistema de vectores), nos queda:

\boldsymbol M'=\boldsymbol M+\left(O-O'\right)\times\sum_{}{\boldsymbol v}_i=\boldsymbol M+\left(O-O'\right)\times\boldsymbol R [3]

Propiedad 5: El momento resultante respecto a O' es la suma del momento resultante respecto a O y el momento respecto a O' del vector  resultante del sistema aplicado en O.

En el caso especial de que la resultante sea nula, R = 0, vemos que debe de ser M = M': para sistemas de vectores con resultante nula el momento del sistema es independiente del punto que se tome. Un sistema de vectores con resultante nula pero con momento total distinto de cero se llama un par vectorial. En el caso particular de vectores fuerza, será un par de fuerzas.

Otro caso especial es cuando el vector OO' es paralelo a la resultante R, pues entonces el producto vectorial de [3] es cero, y el  momento respecto a O es igual al momento respecto a O':

Propiedad 6: El momento resultante de un sistema de vectores es el mismo para cualquier punto situado en una recta paralela a la resultante del sistema.

Trinomio invariante

Si multiplicamos la expresión [3] por la resultante R (producto escalar) obtenemos:

\boldsymbol M'\cdot\boldsymbol R=\boldsymbol M\boldsymbol\cdot\boldsymbol R+\cancel{\left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\times\mathbf R\boldsymbol\cdot\mathbf R}=\boldsymbol M\boldsymbol\cdot\boldsymbol R

luego el producto escalar del momento resultante por la resultante \boldsymbolM\cdot\boldsymbolR es invariante, no depende del punto escogido para el cálculo de los momentos. Cuando se expresa en función de las componentes x, y, z en la forma T=M_xR_x+M_yR_y+M_zR_z se llama trinomio invariante. Para un sistema de un sólo vector T siempre vale cero, pues el vector es perpendicular a su momento.

Eje central de momentos

Momento mínimo

El que la expresión \boldsymbol M'\cdot\boldsymbol R=\boldsymbol M\boldsymbol\cdot\boldsymbol R=\left|M\right|\cdot\left|R\right|\cdot\cos\left(\overset{\boldsymbol\hat{}}{\mathbf M\mathbf R}\right) sea invariante para cualquier punto de cálculo del momento, implica que si en un cierto punto O el momento resultante es paralelo a la resultante R, entonces para ese punto \cos\left(\overset{\boldsymbol\hat{}}{\mathbf M\mathbf R}\right)=1, y comparando con otro punto O' con momento no paralelo a R, será \left|M'\right|\cdot\left|R\right|\cdot cos\left(\widehat{\mathrm{MR}}\right)=\left|M\right|\cdot\left|R\right|,  como cos\left(\widehat{\mathrm{MR}}\right) es en módulo menor que uno (estamos suponiendo que R y M no estan en la misma dirección) , se deduce que \left|M'\right|\cdot\left|R\right|>\left|M\right|\cdot\left|R\right|para M paralelo a R, que significa que el momento resultante M respecto a puntos O en los cuales M es paralelo a la resultante R es el momento mínimo posible.

Lugar geométrico de puntos que generan el mínimo momento

La pregunta siguiente que nos hacemos es: ¿cuáles son esos puntos en los cuales el momento resultante M del sistema es paralelo a la resultante R, y por tanto el momento M es el mínimo posible?

Fijémonos en que si existieran dos puntos O, O' con esa propiedad, tendríamos que sus momentos M, M' son iguales, y por [3], que

\boldsymbol M'=\boldsymbol M+\left(O-O'\right)\times\boldsymbol R=\boldsymbol M\Rightarrow\left(O-O'\right)\times\boldsymbol R=0

por tanto OO' sería un vector paralelo a la resultante R; esto valdría para cualquier punto que proporcionara un momento paralelo a R:

El lugar geométrico de los puntos respecto a los cuales los momentos de un sistema de vectores C son paralelos a la resultante R, es una recta también paralela a R, denominada eje central de momentos del sistema.

Caracterización del eje central de momentos

Encontremos ahora las condiciones que ha de cumplir un punto O' del eje central.  Sea O el origen de coordenadas, y O' el punto de intersección de la perpendicular al eje central que pasa por O (figura 8).

Fig. 8: situación para la deducción de las ecuaciones del eje central

El  momento M' serà paralelo a la resultante R y al eje central; si aplicamos [3], multiplicando vectorialmente toda la expresión por la izquierda por la resultante R obtenemos:

\boldsymbol R\times\boldsymbol M'=\boldsymbol R\times\boldsymbol M+R\times\left(O-O'\right)\times\boldsymbol R

Como R y M' son paralelos, el producto del primer miembro es nulo. Para el doble producto vectorial aplicamos la propiedad A x (B x C) = B · (A·C) - C · (A·B), y resulta:

\begin{array}{l}\mathbf0=\boldsymbol R\times\boldsymbol M+\left(O-O'\right)\cdot\left(\mathbf R\boldsymbol\cdot\mathbf R\right)-\boldsymbol R\boldsymbol\cdot\left(\cancel{\left(\mathbf O\boldsymbol-\mathbf O\boldsymbol'\right)\boldsymbol\cdot\mathbf R}\right)\Rightarrow\\0=\boldsymbol R\times\boldsymbol M+\left(O-O'\right)\cdot R^2\Leftrightarrow\\\left(O'-O\right)=\boldsymbol O\boldsymbol O\boldsymbol'=\frac1{R^2}\boldsymbol R\times\boldsymbol M\end{array} [4]

El producto escalar (O-O')·R es nulo porque OO' es perpendicular a R. La expresión [4] es el vector de posición de un punto O' del eje central, en función de la resultante R y del momento respecto al orígen de coordenadas M.  Ele eje central es la recta paralela a R que pasa por el punto O', quedando así totalmente determinada.

Ejemplo 3: dado el sistema de  vectores deslizantes v(1, 2, 3) aplicado en P(1, 1, 1) y w(0, 1, -1) aplicado en Q(0, 0, 1),  obtener la ecuación de su eje central de momentos, así como el momento mínimo del sistema.

Utilizamos vectores los equipolentes  v' = v - P = (0, 1, 2) y w' = w - Q = (0, 1, -2), obtenemos la resultante R = v' + w' = (0, 2, 0) y seguidamente el momento resultante respecto al origen de coordenadas:

M(v') = (P - O) x v' = [(1, 1, 1) - (0, 0, 0)] x (0, 1, 2) = (1, 1, 1) x (0, 1, 2) = (1, -2, 1);

M(w') = (Q - O) x v' = [(0, 0, 1) - (0, 0, 0)] x (0, 1, -2) = (0, 0, 1) x (0, 1, -2) = (-1, 0, 0);

MM(v') + M(w') = (0, -2, 1), el módulo de M es √5.

Calculamos ahora el punto O' del eje central:

R = |(0, 2, 0)| = 2; R x M = (0, 2, 0) x (0, -2, 1) = (2, 0, 0), por tanto las coordenadas de O' son (1, 0, 0); el eje central pasa por O' y es paralelo a R(0, 2, 0), sus ecuaciones paramétricas son (x, y, z) = (1, 0, 0) + t·(0, 2, 0) -> x = 1, y = t, z = 0. El momento del sistema respecto del eje central, M',  puede obtenerse usando el valor del trinomio invariante y el hecho de que es paralelo a R:

condición trinomio invariante: M'·R = M·R = (0, -2, 1)·(0, 2, 0)= -4

es paralelo a R: M' = t·(0, 2, 0)

Luego M'·R = t·(0, 2, 0)·(0, 2, 0) = 4t, igualando al trinomio invariante, 4t = -4 -> t = -1, el momento M' es (0, -2, 0). Observemos que el módulo de M' es 2, que es menor que el módulo de M, como esperábamos, ya que M' es el momento mínimo.

Reducción de sistemas de vectores deslizantes

Sistemas de vectores equivalentes

Consideramos que dos sistemas de vectores deslizantes son equivalentes si ambos tiene la misma resultante R y el mismo momento resultante M respecto a un punto fijo O. De hecho, aplicando la identidad [3], el momento resultante será también idéntico para cualquier punto del espacio; en efecto, llamando M_1, R_1, M_2,R_2 a los momentos y resultantes del sistema 1 y 2, se cumple:

\begin{array}{l}M_1'=M_1+\left(O-O'\right)\times R_1,\\M_2'=M_2+\left(O-O'\right)\times R_{1;}\\M_1=M_2=M,\;R_1=R_2=R\Rightarrow\\M_1'=M_{}+\left(O-O'\right)\times R_{}=M_2'\end{array}

Operaciones para obtener un sistema equivalente

Las siguientes transformaciones, aplicadas a un sistema de vectores deslizantes, proporcionan otro sistema equivalente:

  1. Decomposición de un vector en varios concurrentes
  2. Composición de vectores concurrentes en uno solo
  3. Sumar o restar dos vectores situados en una misma línea

Respecto a la obtención de un sistema equivalente que sea más simple, tenemos el siguiente teorema:

Teorema: todo sistema de vectores deslizantes puede reducirse a otro que contenga como máximo sólo dos vectores.

Veamos como: si R es la resultante y M el momento resultante respecto a O, un sistema equivalente será el constituido por un vector equipolente a R que pase por O (por tanto de momento nulo) más un par (dos vectores opuestos v, w con resultante nula, uno de ellos, digamos el v, pasando por O) cuyo momento sea M; si sumamos el vector v con R, obtenemos un sistema de sólo dos vectores, (v+R) y w, equivalente al sistema inicial.

 Ejemplo 4:  Sea el sistema de vectores  deslizantes a(1,2,3), b(-1,0,1), c(0, 2, 0), con puntos de aplicación p(0,0,0), q(0, 1,0) y r(0,0,1) respectivamente, y sea el punto O(-1, 2, 1). Encontrar un sistema equivalente de sólo dos vectores.

Primero hallamos las coordenadas de los vectores equipolentes aplicados en el origen:

a' = (1,2,3) - (0,0,0) = (1,2,3); b' = (-1,0,1) - (0, 1,0) = (-1,-1,1); c' = (0, 2, 0) - (0,0,1) = (0, 2, 0).

La resultante será la suma R = a' + b' + c' = (0, 3, 4); el momento resultante respecto a O lo calculamos para cada vector según la definición. recordando de utilizar los vectores equipolentes a', b', c' para los vectores (P-O) y los vectores originales a, b, c para el producto vectorial:

\begin{array}{l}M_a=\left[(0,0,0)-\left(-1,2,1\right)\right]\times\left(1,2,3\right)=\left(-4,-4,4\right);\\M_b=\left[(0,1,0)-\left(-1,2,1\right)\right]\times\left(-1,-1,1\right)=\left(-2,0,-2\right);\\M_c=\left[(0,0,1)-\left(-1,2,1\right)\right]\times\left(0,2,0\right)=\left(0,0,2\right);\\M=M_a+M_b+M_c=\left(-6,-4,4\right).\end{array}

Un sistema equivalente a {a, b, c} aplicados a los puntos p, q, r respecto al punto O será el formado por el vector R aplicado en O más un par de vectores v, w con momento M, estando v aplicado en O.

Como la resultante del par (v,w) ha de ser nula, los vectores equipolentes v', w' situados en el origen de coordenadas han de cumplir v' + w' = 0. Llamando x, y, z a las componentes de v:

v'=\left(x,y,z\right)-\left(-1,2,1\right)=\left(x-1,y-2,z-1\right)\Rightarrow w'=\left(1-x,2-y,1-z\right)

Como el momento de v respecto a O es nulo, ha de ser que el momento de w respecto a O sea igual a M. Lo planteamos así:

M=\left(-6,-4,4\right)=(S-O)\times w'\left(0,0,1\right)

Detallando las coordenadas, y llamando al punto S(s,t,u):

\begin{array}{l}\begin{array}{l}M=\left(-6,-4,4\right)=(s+1,t-2,u-1)\times\left(1-x,2-y,1-z\right)\Rightarrow\\\begin{vmatrix}i&j&k\\s+1&t-2&u-1\\1-x&2-y&1-z\end{vmatrix}=\left(-6,-4,4\right)\end{array}\\\end{array}

Resulta un sistema  de ecuaciones:

\left.\begin{array}{r}\left(t-2\right)\left(1-z\right)-\left(u-1\right)\left(2-y\right)=-6\\-\left(s+1\right)\left(1-z\right)+\left(u-1\right)\left(1-x\right)=-4\\\left(s+1\right)\left(2-y\right)-\left(t-2\right)\left(1-x\right)=4\end{array}\right\}

Tenemos 6 incógnitas, las coordenadas del punto S y las del vector w, y sólo tres ecuaciones; esto significa que tenemos tres grados de libertad al escoger el par, aunque no todas las combinaciones son posibles, pues algunas de ellas pueden conducirnos a sistemas incompatibles. Una posible elección es la que sigue (se definen  algunos parámetros y de ellos se deducen los otros):

\begin{array}{l}t=2;\;s=2;\Rightarrow y=\frac23,\;u=\frac{11}2;\\\;x=z\Rightarrow x=-\frac{11}6=z\\\\\end{array}

Los vectores del par resultan:

v'=\left(-\frac{11}6,\frac23,-\frac{11}6\right),\;w'=\left(1+\frac{11}6,2-\frac23,1+\frac{11}6\right)=\left(\frac{17}6,\frac43,\frac{17}6\right)

y estan aplicados en los puntos O(-1, 2, 1) y S(2,2,11/2). Junto con el vector R'= (0, 3, 4) aplicado también en O, constituyen un sistema equivalente al inicial. Sumamos ahora R' y v' (ámbos aplicados al mismo punto O) para obtener R''= (0, 3, 4) + (-11/6, 2/3, -11/6) = (-11/6, 11/6, -11/6) = (11/6)·(-1, 1, -1), vector aplicado en O. Hemos reducido el sistema original a sólo dos vectores:

  1. R''=\frac{11}6\left(-1,1,-1\right) aplicado en O(-1, 2, 1)
  2. vector equipolente w'=\left(\frac{17}6,\frac43,\frac{17}6\right) aplicado en S(2,2,11/2)

 

 

Publicado en Fí­sica, Vectores | Etiquetado , , , , , , | Deja un comentario

Movimiento ondulatorio

Movimiento ondulatorio: ejemplos

Cuando arrojamos una piedra en un estanque vemos que se forman en el agua ondulaciones en forma de círculos concéntricos, que se extienden a partir del punto en que se produjo la perturbación.

Fig. 1: ondas en el agua. Fuente: http://hyperphysics.phy-astr.gsu.edu/hbasees/Sound/wavplt.html

Si a  una cuerda sujeta por un extremo le damos un movimiento de vaivén, veremos que se forman ondulaciones que se trasladan a lo largo de la cuerda.

Fig.2: Al mover la mano arriba y abajo la cuerda oscila. Los puntos de la cuerda se mueven transversalmente a la dirección de propagación de la oscilación.

Si ejercemos una fuerza bruscamente sobre un extremo de un resorte que tenemos sujeto por el otro extremo, y después lo liberamos, el resorte oscilará; observando su compresión, podremos ver que hay una zona de compresión que avanza a lo largo del muelle. Cuanto más largo sea el resorte y menor su rigidez (constante elástica), más fácil será observar la onda de compresión.

Fig. 3: al comprimir rápidamente un muelle podemos observar una onda de compresión que avanza por él.

En todos los movimientos ondulatorios el modelo físico es el mismo: una serie de partículas que ejercen un movimiento oscilante idéntico, pero cada una de ellas con cierto retraso (diferencia de fase) entre ellas, además las partículas no se trasladan en la dirección de propagación de la onda, como se muestra en el ejemplo de la figura 4.

Fig. 4: las olas no transportan agua, el agua sube y baja al paso de las olas, como indica el movimiento de la barca, que no avanza.

Otra propiedad importante de las ondas es su independencia entre sí, si han estado generadas por sucesos independientes: si en las ondas de la figura 1 tiramos otra piedra en un punto cualquiera, se forma otro sistema de ondas con centro en ese punto, que en todo caso puede interferir con el original, sumándose o restándose en cada punto. Es por este motivo que, cuando oímos sonidos procedentes de varias fuentes, podemos distinguirlas por separado, ya que cada sonido es independiente del otro.

Ondas transversales y ondas longitudinales

En las ondas en el agua o en una cuerda, la oscilación de cada partícula se produce en una dirección transversal a la dirección de propagación de la perturbación, hay unas crestas y unos valles, por ello se llaman ondas transversales. En cambio en un resorte las partículas vibran en la misma dirección de propagación de la perturbación, son una sucesión de compresiones y dilataciones, por lo que se llaman ondas longitudinales. (figura 5).

Fig. 5: Ejemplos de ondas longitudinales y transversales

En el caso de las ondas longitudinales las partículas, al desplazarse de su posición de equilibrio, presionan a sus vecinas, poniéndolas también en movimiento siempre que el medio sea elástico, ya que si fuera totalmente rígido no podría oscilar. En las ondas transversales para que el movimiento de una partícula afecte (arrastre) a las vecinas, es necesario que entre ellas exista cierta cohesión, que sólo existe en el caso de medios sólidos, pues en los fluidos cada partícula se mueve con fluidez respecto a las otras, por ello estas ondas en general solo existen en los sólidos; la excepción son las ondas transversales que se propagan por la superficie de un líquido, como las olas del mar, en ellas la cohesión la aporta la tensión superficial del líquido.

Matemática del movimiento ondulatorio

En la figura 5 vemos una instantánea de un movimiento ondulatorio que se propaga según el eje horizontal X; cada punto, como por ejemplo los puntos A, B, oscilan verticalmente según un eje vertical Y, con una amplitud A. A medida que la onda se propaga, alcanza más puntos que estaban en reposo que pasan también a oscilar. La distancia entre "crestas" es la longitud de onda \lambda. Los puntos que están en el mismo estado oscilatorio, o sea tienen la misma ordenada Y, y la misma velocidad v, se dicen que están en fase, y la distancia que los separa también es igual a la longitud de onda \lambda.

Fig. 5: longitud de onda

El movimiento oscilatorio vertical de amplitud A de cada punto (figura 6) se puede describir imaginando un movimiento circular asociado de forma que al punto oscilante P le sigue como su sombra el punto P', que localizamos por su posición angular \theta. Si tomamos Y = 0 en el punto medio de la oscilación, que corresponde a \theta=0, entonces la posición Y = +A/2 corresponde a \theta=\pi/2, y en general, Y = R\sin(\theta), siendo R el radio del movimiento circular, R = A/2.

Fig. 6: descripción del movimiento oscilatorio

Para introducir el tiempo t en la ecuación, definimos la frecuencia \nu de la oscilación como el número de oscilaciones completas realizadas por unidad de tiempo; en el movimiento circular asociado, cada oscilación corresponde a una rotación completa \theta=2\pi, por tanto la velocidad angular será \nu\cdot2\pi, y el ángulo girado en un tiempo t es \theta=t·\nu\cdot2\pi\cdot. Sustituyendo:

Y(t) = (A/2)\sin(t·\nu\cdot2\pi)

Hemos obtenido el desplazamiento vertical Y de un punto cualquiera en función del tiempo, sin especificar su posición horizontal X observamos que dos puntos separados una distancia igual a la longitud de onda, \triangle X=\lambda, deben oscilar en fase (el ángulo theta del movimiento circular asociado ha de ser idéntico), o sea, Y\left(x,t\right)=Y\left(x+\triangle X,t\right)=Y\left(x+\lambda,t\right); por ello, el desplazamiento Y(x,t) ha se ser una función periódica en la variable x con período \lambda. La función \sin(x) tiene periodo 2\pi, si queremos que tenga periodo \lambda hay que hacer un cambio de escala del eje X, transformándolo en X',  de forma que la unidad original pasa a valer \lambda, esto es, \frac{2\pi}\lambda=\frac x{x'}\Rightarrow x=\frac{\displaystyle2\pi}{\displaystyle\lambda}x'. Nos queda el detalle de la dirección de transmisión de la onda: si suponemos que va de izquierda a derecha, entonces el ángulo de fase decrecerá conforme x aumente, y suficientemente lejos del origen valdrá cero (la onda todavía no habrá llegado). Así pues, ensayemos la función:

Y\left(x,t\right)=\frac A2\sin(2\pi\nu t-\frac{2\pi}\lambda x)=\frac A2\sin(2\pi\left(\nu t-\frac x\lambda\right))

Observemos que el término en x se resta del ángulo total, tal como hemos comentado. Comprobemos si es periódica en X con periodo \lambda :

\begin{array}{l}Y\left(x+\lambda,t\right)=\frac A2\sin(2\pi\left(\nu t-\frac{x+\lambda}\lambda\right))=\frac A2\sin(2\pi\left(\nu t-\frac x\lambda+1\right))=\\\frac A2\sin(2\pi\left(\nu t-\frac x\lambda\right)+2\pi)=\frac A2\sin(2\pi\left(\nu t-\frac x\lambda\right))=\\Y\left(x,t\right)\end{array}

y cumple con la condición exigida. Llamando Y_0 a A/2 y definiendo el período de oscilación, T=1/\nu, tiempo transcurrido en una oscilación completa de cada punto, la ecuación de la onda nos queda en la forma:

Y\left(x,t\right)=Y_0\sin(2\pi\left(\frac tT-\frac x\lambda\right))

Para el instante inicial t = 0, y en el punto inicial x = 0, resulta Y(0, 0) = 0, que no siempre será cierto, por ejemplo en la figura 2 movemos inicialmente la cuerda para que empiece a oscilar, por tanto Y no será cero; para tener en cuenta este hecho, consideramos que la fase puede ser distinta de cero en el instante inicial, y definimos esa fase inicial por  \theta_0, con ello resulta la ecuación de onda (unidimensional, en la dirección X, propagación según la dirección X > 0):

Y\left(x,t\right)=Y_0\sin(2\pi\left(\frac tT-\frac x\lambda\right)+\theta_0) [1]

Ejemplo 1: en la figura 2, en el tiempo inicial t = 0, el extremo de la cuerda está desplazado verticalmente 20cm hacia arriba, y a continuación se hace oscilar con una frecuencia de oscilación es de 2 ciclos por segundo.  La longitud de la onda que se forma en la cuerda es de 40cm. Cuál será la posición Y del punto x = 100cm en el instante t = 5s? Supondremos que la perturbación se mueve con una velocidad constante.

Identificamos parámetros de la ecuación [1]: Y(0, 0) = 0.2,\; Y_0=0.2,\;\lambda=0.4,\;\nu=2, sustituyendo:

Y\left(x,t\right)=0.2\sin(2\pi\left(\frac t{1/2}+\frac x{0.4}\right)+\theta_0)

donde hemos usado T=1/\nu; para encontrar la fase inicial sustituimos el valor de Y(0,0):

\begin{array}{l}Y\left(0,0\right)=0.2=0.2\sin(2\pi\left(\frac0{1/2}+\frac0{0.4}\right)+\theta_0)=0.2\sin(\theta_0)\Rightarrow\\\sin(\theta_0)=1\Rightarrow\theta_0=\frac{\mathrm\pi}2\end{array}

Ya tenemos la ecuación de la onda en la cuerda:

Y\left(x,t\right)=0.2\sin(2\pi\left(2t+\frac x{0.4}\right)+\frac{\mathrm\pi}2)

Contestemos la pregunta: cuál será la posición Y del punto x = 100cm en el instante t = 5s?

Y\left(1,5\right)=0.2\sin(2\pi\left(10+\frac1{0.4}\right)+\frac{\mathrm\pi}2)=-0.2

Pero, un momento, esa posición sólo será correcta si, en el tiempo t = 5 la onda ya ha alcanzado la posición x = 100cm, en otro caso, la cuerda estaría todavía en reposo; tenemos que comprobarlo: la velocidad de transmisión, también conocida por velocidad de fase, en el sentido en que es la velocidad a la que la fase de la oscilación se propaga a lo largo de la cuerda, es v=\lambda\cdot\nu=0.4/2=0.2m/s ya que velocidad = espacio/tiempo, la longitud de onda es el espacio recorrido, y la frecuencia es el número de veces que se recorre esa distancia por segundo, a esa velocidad, en 5 segundos la onda se hallará en x = v·t = 0.2·5 = 1m, justo acabará de llegar a ese punto, nuestra respuesta era correcta.

Ondas de propagación

Las ondas que hemos visto hasta ahora son generadas desde una fuente (el punto donde lanzamos la piedra en el estanque, el extremo de la cuerda que agitamos...) y se propagan por un medio que, en teoría, es ilimitado; si no lo fuera, la onda llegaría a un "final de trayecto" donde se reflejaría, y al reflejarse interferiría consigo misma, creando patrones de ondas que no estudiaremos en este artículo. Además, suponemos que la fuente de la perturbación no se detiene, sino que sigue impulsando la creación de ondas (la mano que hace oscilar la cuerda no se detiene), si no fuera así, la onda se amortiguaría rápidamente debido a la resistencia del medio (ondas amortiguadas). Cuando tenemos una fuerza motriz actuando que genera las oscilaciones, diremos que formamos oscilaciones forzadas. Las oscilaciones forzadas aplicadas a un medio ilimitado generan las denominadas ondas de propagación.

La fuerza motriz puede actuar sólo en un instante (la piedra que cae en el agua) o continuamente, en este último caso si se aplica una oscilación armónica (sinusoidal) obtendremos ondas armónicas de propagación.

Publicado en Fí­sica, Mecánica, Oscilaciones y ondas | Etiquetado , , , , , , , , , | Deja un comentario

Acústica: música

La sucesión de sonidos con ciertas características especiales la denominamos música. El oído es capaz de percibir, dados dos sonidos emitidos en sucesión, la relación entre sus frecuencias (o tonos), esa relación se denomina intervalo musical. Por ejemplo, dados dos sonidos a 300Hz y 150Hz, el oído notará que entre ellos existe la misma relación que entre otros dos de 800Hz y 400Hz, ya que el intervalo es el mismo para las dos secuencias, 300:150 = 800:400 = 2. En la siguiente grabación puede escucharse dos secuencias de 11 notas cada una formando las denominadas escalas musicales (concretamente la escala pentatónica);  las relaciones de frecuencias entre las notas son las mismas en ambas escalas, pero la primera escala tiene una frecuencia superior.

Si se producen dos sonidos simultáneamente, o muy cercanos en el tiempo, nos producirá una sensación agradable (acorde musical) cuanto más simple sea su intervalo musical. Los más simples son 2:1 (intervalo llamado octava), 3:2 (quinta), 4:3 (cuarta) y 5:4 (tercera mayor), que eran los más utilizados en la música clásica antigua.

Escalas musicales

Como la musicalidad de una secuencia de sonidos no depende de sus frecuencias sino de la relación entre sus frecuencias, podemos escoger libremente una frecuencia base cualquiera. Tomemos un sonido de frecuencia f como base (tono fundamental, o tónico), y lo denominamos do. A continuación definimos los intervalos musicales, por ejemplo, si tomamos los intervalos (5:4)f, (4:3)f, (3:2)f, (5:3)f y (2:1)f obtendremos los sonidos que, relativos al base do, se denominan mi, fa, sol, la, do'.  Habremos obtenido una escala musical. ¿Que relación guardan entre si las frecuencias de esta escala? En la siguiente tabla vemos las relaciones f:f' de cada nota respecto a la fundamental do, una columna donde, multiplicando las relaciones f:f' por el factor 12 se convierten a enteros, y una última columna que descompone en producto de dos factores esos enteros.

 Nota f:f’ 12·f:f’ factores
do 1 12 4x3
mi 1 1/4 15 5x3
fa 1 1/3 16 4x4
sol 1 1/2 18 6x3
la 1 2/3 20 5x4
do’ 2 24 6x4

Hay subgrupos de tres acordes con relaciones especialmente simples: do-mi-sol y fa-la-do', sus tablas de relaciones son:

 Nota f:f’ 4·f:f’
do 1 4
mi 1 1/4 5
sol 1 1/2 6
 Nota f:f’ 3·f:f’
fa 1 1/3 4
la 1 2/3 5
do’ 2 6

Vemos que sus frecuencias tienen relaciones simples iguales entre sí, de 4:5:6. La denominada escala diatónica natural añade otras dos notas, siendo 8 en total: re, con relación (9:8)f, y si, con relación (15/8)f. La tabla de relaciones queda:

Nota f:f’ 24·f:f’ factores
do (base)
1 24 8x3
re (segunda mayor)
1 1/8 27 9x3
mi (tercera mayor)
1 1/4 30 6x5
fa (cuarta)
1 1/3 32 8x4
sol (quinta)
1 1/2 36 9x4
la (sexta)
1 2/3 40 8x5
si (séptima)
1 7/8 45 9x5
do’ (octava)
2 48 8x6

Con estas dos notas, encontramos un nuevo acorde de tres notas especialmente simple: re-sol-si, con una relación mutua de frecuencias de  3:4:5. En Música, además de la nomenclatura do, re, mi, etc. se usa también la abreviada con letras A, B, C, ...

Instrumentos musicales de cuerda

Una cuerda tensa, pulsada o bien frotada, oscilará con un movimiento armónico simple amortiguado; al hacerlo, presionará el aire que la rodea de forma que generará una onda de presión de frecuencia constante: un sonido de tono dado. La frecuencia de vibración es más alta como menor sea la longitud de la cuerda, y como más tensa esté; el grosor de la cuerda también interviene, cuanto más grosor, menos frecuencia.

Instrumentos como el piano o el arpa tienen una cuerda dedicada a cada nota (o frecuencia); en la guitarra acústica y en el violín el músico varía la frecuencia modificando la longitud de la cuerda vibrante. Si pulsamos la cuerda no por el centro exacto sino por diferentes puntos se obtendrá la misma nota pero con diferente timbre (vibraciones adicionales, o armónicos, añadidas a la vibración fundamental de la cuerda), y un oído entrenado podrá detectarlo.

Fig. 1: oscilaciones del tono fundamental y del primer armónico (octava del fundamental)

Cuerdas vibrantes en música

Una cuerda es un hilo elástico con sección pequeña comparada con su longitud. Cuando una cuerda tensada por sus extremos se golpea, pulsa o frota en un punto, la perturbación se propaga por la cuerda hasta los extremos, donde se reflejan e interactúan entre sí formando vibraciones en forma de ondas estacionarias con los nodos en los extremos, y cada punto de la cuerda describe un movimiento armónico simple; para ello, la cuerda ha de contener o bien media onda (dos nodos), o una onda completa (tres nodos, dos en los extremos y uno en el centro), o una onda y media (cinco nodos), y en general un número impar de nodos y un múltiplo entero de "medias ondas completas".  Por ello, las vibraciones de la cuerda cumplen que la longitud l de la cuerda ha de ser un múltiplo entero de la semilongitud de la onda: l=n\frac\lambda2, y la frecuencia sonora será f=\fracv\lambda=n\fracv{2l} para n = 1, 2, 3, .... , donde v es la velocidad de transmisión de la perturbación en la cuerda, que depende del material y de la tensión de la cuerda. A estas frecuencias se les llama frecuencias propias de la cuerda. Si intentamos hacer vibrar a una cuerda con frecuencias distintas de las propias, la vibración se amortiguará muy rápidamente debido a que la interferencia de la onda transmitida y la reflejada será destructiva: una anula a la otra.

Escala atemperada

En las guitarras, también en los violines, el músico puede variar la frecuencia del sonido a voluntad de forma continúa, variando la posición del dedo que pulsa la cuerda.

Fig. 2: Parte del teclado de un piano mostrando la escala atemperada de 13 notas, las teclas negras son los semitonos.

En cambio en el piano, instrumento en el que la cuerda se golpea con un pequeño martillo siempre en la misma posición, las frecuencias son fijas. Cuando se quiere acompañar con el instrumento musical a la voz de una cantante, con las frecuencias variables se puede ajustar la escala musical a la frecuencia de la voz

demos hacerlo, pero para solventar en parte este problema se recurre a la denominada escala atemperada: se parte de la escala natural que se divide en 12 intervalos iguales (por tanto contiene 13 notas) añadiendo nuevas notas llamadas semitonos, o notas sostenidas, marcadas con el símbolo #. Las nuevas notas sostenidas se destacan en negro.

La escala atemperada completa queda así:

do, do#, re, re#, mi, fa, fa#, sol, sol#, la, la#, si, do'

Al dividir la escala natural en 12 intervalos en vez de los 7 originales sucede que las frecuencias de las notas de la escala atemperada no coinciden con la natural, excepto en las notas inicial y final (do y do'). Esta escala facilita la acomodación de frecuencias a la voz, pero incluyendo una pequeña variación de las relaciones de frecuencias, perceptible al oído entrenado musicalmente.

Más sobre escalas musicales.

Instrumentos de viento

Cuando se vacía una vasija de agua, se percibe el borboteo como una sucesión de notas que van disminuyendo de tono a medida que se vacía, o sea, cuando la cavidad de aire se hace mayor. Se produce una vibración del aire (onda longitudinal) dentro de la vasija, un sonido, con una frecuencia inversamente proporcional a la longitud de la columna de aire vibrante, igual que sucede con las cuerdas vibrantes.

En los tubos sonoros, cilindros llenos de aire, se generan sonidos al hacer vibrar el aire que contienen. También se presentan modos fundamentales de vibración y modos armónicos en las columnas de aire, sólo que en este caso hay que diferenciar los casos en que el tubo sea semi-abierto al aire (fig. 3), o sea cerrado sólo en un extremo, o totalmente abierto; en el primer caso, en la pared de cierre para todos los modos vibratorios siempre habrá un nodo (punto de vibración nula), mientras que para los tubos abiertos los nodos están en el interior, y en los extremos abiertos siempre encontramos vientres (puntos de vibración máxima).

Fig. 3: tercer armónico y modo fundamental de vibración en un tubo sonoro semiabierto

En los tubos sonoros la onda longitudinal de aire generada en un extremo se refleja en el otro extremo, la onda interfiere consigo misma, formando un sistema de ondas estacionarias. Si el tubo es semiabierto, la onda al chocar con la pared del extremo cerrado cambia 90⁰ de fase, y las dos ondas, transmitida y reflejada, se anulan en ese punto, formando un nodo; en el extremo abierto, que es donde generamos la vibración, tenemos un vientre. En el modo fundamental (fig. 3 a la derecha) en el tubo hay un sólo nodo y un sólo vientre, como la onda completa comprende dos nodos y dos vientres, llamando \lambda a la longitud de onda y l a la longitud del tubo, tendremos la relación, para el modo fundamental, l=\frac\lambda4, y la frecuencia f de ese modo será f=\fracv\lambda=\fracv{4l}, donde v es la velocidad del sonido. En general si se presentan n nodos en el tubo, estando uno de ellos siempre en el extremo cerrado, teniendo en cuenta que cada dos nodos corresponden a una longitud de onda de distancia, tendremos que los n nodos ocupan dentro del tubo una distancia n\lambda/2, desde el último nodo hasta el extremo abierto del tubo hay una distancia \lambda/4 (pues hay un vientre en el extremo), y por tanto la longitud total l del tubo se distribuye así:

l=n\frac\lambda2+\frac\lambda4=\frac{2n\lambda+\lambda}4=\left(2n+1\right)\frac\lambda4

y por tanto la frecuencia del modo con n nodos es:

f=\frac v\lambda=\left(2n+1\right)\frac v{4l}.

El factor (2n + 1) nos indica que en los tubos semiabiertos sólo se emiten los armónicos impares (las frecuencias de los armónicos sin múltiplos impares de la frecuencia fundamental). Queda como ejercicio demostrar que esto no sucede en los tubos abiertos.

En los órganos hay un tubo de longitud fija para cada tono, de forma semejante a los pianos (fig. 4), en cambio en la flauta la longitud efectiva del tubo la puede variar el músico con los dedos.

Fig. 4: tubos del órgano de la basílica de la Sagrada Familia, en Barcelona

Resonancia acústica

Dos sistemas vibratorios con la misma frecuencia pueden interactuar y producir fenómenos de resonancia, en el que las energías de los dos sistemas se suman, o bien se transmite la energía de un sistema al otro.

Fig. 5: resonancia entre un diapasón y una columna de aire cercana

En la figura 5 se presenta un experimento de resonancia acústica: tomamos un tubo semiabierto con agua dispuesto de forma que podamos variar la altura del líquido (vaso comunicante), cerca del extremo abierto hacemos sonar un diapasón, veremos que el tubo también emite un sonido, pues el aire que contiene vibra por efecto del diapasón. Si vamos probando a variar la altura del líquido hasta que una de las frecuencias del tubo coincida con la del diapasón, más exactamente, cuando la longitud de la columna de aire sea un múltiplo de un cuarto de la longitud de onda del diapasón, en ese momento entrarán en resonancia y el tubo emitirá un sonido intenso, ya que se estará transmitiendo eficientemente la energía vibratoria del diapasón al aire del tubo.  Esta transmisión eficiente provoca que el sonido del diapasón se disipe rápidamente, pues transfiere toda su energía al aire del tubo. Las cajas de madera de algunos instrumentos (guitarra, violín) se utilizan para reforzar los sonidos por resonancia, y por ello se denominan cajas de resonancia.

Bibliografía usada

Ambas obras están descatalogadas, siendo textos excelentes, he intentado aprovechar la información que daban para que no se pierda en el olvido. La figura 5 proceden de la primera obra:

  • M. Catalán, Andres León: Física y Química – Madrid, 1945
  • J. Fernandez, M. Pujal: Iniciación a la Física – Barcelona, 1975

 

Publicado en Acústica, Fí­sica | Etiquetado , , , , , , , , , | Deja un comentario