Ecuación de onda - Wave equation

Un pulso que viaja a través de una cuerda con puntos finales fijos según el modelo de la ecuación de onda.
Ondas esféricas provenientes de una fuente puntual.
Una solución a la ecuación de onda 2D

La ecuación de onda es una ecuación diferencial parcial lineal de segundo orden para la descripción de ondas, como ocurren en la física clásica , como las ondas mecánicas (por ejemplo , ondas de agua , ondas de sonido y ondas sísmicas ) u ondas de luz . Surge en campos como la acústica , electromagnética y dinámica de fluidos . Debido al hecho de que la segunda ecuación de onda orden describe la superposición de una onda entrante y una onda de salida (es decir, en lugar de un campo de onda estacionaria) que también se llama "ecuación de onda de dos vías" (en contraste, la primera orden de una vía La ecuación de onda describe una sola onda con una dirección de propagación de onda predefinida y es mucho más fácil de resolver debido a las derivadas de primer orden).

Históricamente, el problema de una cuerda vibrante como la de un instrumento musical fue estudiado por Jean le Rond d'Alembert , Leonhard Euler , Daniel Bernoulli y Joseph-Louis Lagrange . En 1746, d'Alembert descubrió la ecuación de onda unidimensional y, en diez años, Euler descubrió la ecuación de onda tridimensional.

Introducción

La ecuación de onda (bidireccional) es una ecuación diferencial parcial de segundo orden que describe las ondas. Este artículo se centra principalmente en la ecuación de onda escalar que describe ondas en escalares mediante funciones escalares u = u ( x 1 , x 2 ,…, x n ; t ) de una variable de tiempo t (una variable que representa el tiempo) y una o más variables espaciales x 1 , x 2 ,…, x n (variables que representan una posición en un espacio en discusión) mientras que hay ecuaciones de ondas vectoriales que describen ondas en vectores como ondas para campo eléctrico, campo magnético y potencial vectorial magnético y ondas elásticas . En comparación con las ecuaciones de ondas vectoriales, la ecuación de ondas escalares puede verse como un caso especial de las ecuaciones de ondas vectoriales; En el sistema de coordenadas cartesianas , la ecuación de onda escalar es la ecuación que debe satisfacer cada componente (para cada eje de coordenadas, como el componente x del eje x) de una onda vectorial sin fuentes de ondas en el dominio considerado ( es decir, un espacio y un tiempo). Por ejemplo, en el sistema de coordenadas cartesiano, ya que como representación de una onda de campo de vector eléctrico en ausencia de fuentes de onda, cada componente del eje de coordenadas ( i = x , y , o z ) debe satisfacer la ecuación de onda escalar. Otras soluciones de ecuaciones de ondas escalares u son para cantidades físicas en escalares como la presión en un líquido o gas, o el desplazamiento , a lo largo de alguna dirección específica, de partículas de un sólido en vibración alejándose de sus posiciones de reposo (equilibrio).

La ecuación de onda escalar es

donde c es un coeficiente real fijo no negativo .

Usando las notaciones de la mecánica newtoniana y el cálculo vectorial , la ecuación de onda se puede escribir de manera más compacta como

donde el punto doble denota la derivada de tiempo doble de u , es el operador nabla , y 2 = ∇ · ∇ es el operador laplaciano (espacial) (no vector laplaciano):

Una notación aún más compacta que a veces se usa en física se lee simplemente

donde todos los operadores se combinan en el operador d'Alembert :

Una solución de esta ecuación puede ser bastante complicada, pero se puede analizar como una combinación lineal de soluciones simples que son ondas planas sinusoidales con varias direcciones de propagación y longitudes de onda, pero todas con la misma velocidad de propagación c . Este análisis es posible porque la ecuación de onda es lineal ; de modo que cualquier múltiplo de una solución también es una solución, y la suma de dos soluciones cualesquiera es nuevamente una solución. Esta propiedad se llama principio de superposición en física.

La ecuación de onda por sí sola no especifica una solución física; Por lo general, se obtiene una solución única planteando un problema con condiciones adicionales, como las condiciones iniciales , que prescriben la amplitud y la fase de la onda. Otra clase importante de problemas ocurre en espacios cerrados especificados por condiciones de contorno , para los cuales las soluciones representan ondas estacionarias , o armónicos , análogos a los armónicos de los instrumentos musicales.

La ecuación de onda es el ejemplo más simple de una ecuación diferencial hiperbólica . Este, y sus modificaciones, juegan un papel fundamental en la mecánica del continuo , la mecánica cuántica , la física del plasma , la relatividad general , la geofísica y muchas otras disciplinas científicas y técnicas.

Ecuación de onda en una dimensión espacial

El científico francés Jean-Baptiste le Rond d'Alembert descubrió la ecuación de onda en una dimensión espacial.

La ecuación de onda en una dimensión espacial se puede escribir de la siguiente manera:

Esta ecuación se describe típicamente como teniendo solo una dimensión espacial x , porque la única otra variable independiente es el tiempo t . No obstante, la variable dependiente u puede representar una segunda dimensión espacial, si, por ejemplo, el desplazamiento u se produce en la dirección y , como en el caso de una cuerda que se ubica en el plano xy .

Derivación de la ecuación de onda

La ecuación de onda en una dimensión espacial se puede derivar en una variedad de entornos físicos diferentes. Lo más famoso es que se puede derivar para el caso de una cuerda que vibra en un plano bidimensional, con cada uno de sus elementos tirados en direcciones opuestas por la fuerza de la tensión .

Otro escenario físico para la derivación de la ecuación de onda en una dimensión espacial utiliza la Ley de Hooke . En la teoría de la elasticidad , la ley de Hooke es una aproximación para ciertos materiales, que establece que la cantidad por la cual se deforma un cuerpo de material (la deformación ) está relacionada linealmente con la fuerza que causa la deformación (la tensión ).

De la ley de Hooke

La ecuación de onda en el caso unidimensional se puede derivar de la ley de Hooke de la siguiente manera: imagine una matriz de pequeños pesos de masa m interconectados con resortes sin masa de longitud h . Los resortes tienen una constante de resorte de k :

Matriz de masas.svg

Aquí la variable dependiente u ( x ) mide la distancia desde el equilibrio de la masa situada en x , de modo que u ( x ) mide esencialmente la magnitud de una perturbación (es decir, deformación) que viaja en un material elástico. Las fuerzas ejercidas sobre la masa m en el lugar x + h son:

La ecuación de movimiento para el peso en la ubicación x + h se obtiene al igualar estas dos fuerzas:

donde la dependencia temporal de u ( x ) se ha hecho explícita.

Si la matriz de pesos consta de N pesos espaciados uniformemente sobre la longitud L = Nh de la masa total M = Nm , y la constante de resorte total de la matriz K = k / N , podemos escribir la ecuación anterior como:

Tomando el límite N → ∞, h → 0 y asumiendo suavidad se obtiene:

que es de la definición de una segunda derivada . KL 2 / M es el cuadrado de la velocidad de propagación en este caso particular.

Onda estacionaria 1-d como superposición de dos ondas que viajan en direcciones opuestas

Pulso de estrés en un bar

En el caso de un impulso de tensión que se propaga longitudinalmente a través de una barra, la barra actúa como un número infinito de resortes en serie y puede tomarse como una extensión de la ecuación derivada de la ley de Hooke. Una barra uniforme, es decir, de sección transversal constante, hecha de un material elástico lineal tiene una rigidez K dada por

Donde A es el área de la sección transversal y E es el módulo de Young del material. La ecuación de onda se convierte en

AL es igual al volumen de la barra y por lo tanto

donde ρ es la densidad del material. La ecuación de onda se reduce a

Por tanto, la velocidad de una onda de tensión en una barra es E / ρ .

Solución general

Enfoque algebraico

La ecuación de onda unidimensional es inusual para una ecuación diferencial parcial en el sentido de que se puede encontrar una solución general relativamente simple. Definición de nuevas variables:

cambia la ecuación de onda en

que conduce a la solución general

o equivalente:

En otras palabras, las soluciones de la ecuación de onda 1D son sumas de una función de desplazamiento a la derecha F y una función de desplazamiento izquierda G . "Viajar" significa que la forma de estas funciones arbitrarias individuales con respecto ax permanece constante, sin embargo, las funciones se trasladan de izquierda a derecha con el tiempo a la velocidad c . Esto fue derivado por Jean le Rond d'Alembert .

Otra forma de llegar a este resultado es observar que la ecuación de onda se puede "factorizar" en dos ecuaciones de onda unidireccionales :

Como resultado, si definimos v así,

luego

De esto, v debe tener la forma G ( x + ct ) , y de esto se puede deducir la forma correcta de la solución completa u . Además de la descomposición matemática de la ecuación de onda de segundo orden, la ecuación de onda unidireccional también se puede derivar directamente de la impedancia.

Para un problema de valor inicial, las funciones arbitrarias F y G se pueden determinar para satisfacer las condiciones iniciales:

El resultado es la fórmula de d'Alembert :

En el sentido clásico, si f ( x ) ∈ C k y g ( x ) ∈ C k −1 entonces u ( t , x ) ∈ C k . Sin embargo, las formas de onda F y G también pueden ser funciones generalizadas, como la función delta. En ese caso, la solución puede interpretarse como un impulso que viaja hacia la derecha o hacia la izquierda.

La ecuación de onda básica es una ecuación diferencial lineal y, por lo tanto, se adherirá al principio de superposición . Esto significa que el desplazamiento neto causado por dos o más ondas es la suma de los desplazamientos que habría causado cada onda individualmente. Además, el comportamiento de una onda se puede analizar dividiendo la onda en componentes, por ejemplo, la transformada de Fourier rompe una onda en componentes sinusoidales.

Modos propios de onda plana

Otra forma de resolver la ecuación de onda unidimensional es analizar primero sus modos propios de frecuencia . Un llamado modo propio es una solución que oscila en el tiempo con una frecuencia angular constante bien definida ω , de modo que la parte temporal de la función de onda toma la forma e - iωt = cos ( ωt ) - i sin ( ωt ) , y la amplitud es una función f ( x ) de la variable espacial x , dando una separación de variables para la función de onda:

Esto produce una ecuación diferencial ordinaria para la parte espacial f ( x ) :

Por lo tanto:

que es precisamente una ecuación de valor propio para f ( x ) , de ahí el nombre modo propio. Tiene las conocidas soluciones de onda plana .

con número de onda k = ω / c .

La función de onda total para este modo propio es entonces la combinación lineal

donde los números complejos A, B dependen en general de cualquier condición inicial y de frontera del problema.

Los modos propios son útiles para construir una solución completa de la ecuación de onda, porque cada uno de ellos evoluciona trivialmente en el tiempo con el factor de fase . para que una solución completa se pueda descomponer en una expansión de modo propio

o en términos de ondas planas,

que está exactamente en la misma forma que en el enfoque algebraico. Las funciones s ± ( ω ) se conocen como componente de Fourier y están determinadas por las condiciones iniciales y de contorno. Este es un método denominado en el dominio de la frecuencia , alternativo a las propagaciones directas en el dominio del tiempo , como el método FDTD , del paquete de ondas u ( x , t ) , que es completo para representar ondas en ausencia de dilataciones en el tiempo. La completitud de la expansión de Fourier para representar ondas en presencia de dilataciones de tiempo ha sido desafiada por soluciones de onda de chirrido que permiten una variación de tiempo de ω . Las soluciones de onda de chirrido parecen estar particularmente implicadas por residuos de radar muy grandes pero previamente inexplicables en la anomalía de sobrevuelo , y difieren de las soluciones sinusoidales en que pueden recibirse a cualquier distancia solo en frecuencias desplazadas proporcionalmente y dilataciones de tiempo, correspondientes a estados de chirp pasados ​​de la fuente.

Ecuación de onda escalar en tres dimensiones espaciales

El matemático y físico suizo Leonhard Euler (n. 1707) descubrió la ecuación de onda en tres dimensiones espaciales.

Se puede obtener una solución del problema de valor inicial para la ecuación de onda en tres dimensiones espaciales a partir de la solución correspondiente para una onda esférica. El resultado también se puede utilizar para obtener la misma solución en dos dimensiones espaciales.

Ondas esféricas

La ecuación de onda se puede resolver mediante la técnica de separación de variables . Para obtener una solución con frecuencias constantes, primero transformemos de Fourier la ecuación de onda en el tiempo como

Entonces obtenemos,

Esta es la ecuación de Helmholtz y se puede resolver mediante la separación de variables. Si se utilizan coordenadas esféricas para describir un problema, entonces la solución a la parte angular de la ecuación de Helmholtz viene dada por armónicos esféricos y la ecuación radial ahora se convierte en

Aquí kω / cy la solución completa ahora está dada por

donde h(1)
l
( KR )
y h(2)
l
( kr )
son las funciones esféricas de Hankel .

Ejemplo

Para comprender mejor la naturaleza de estas ondas esféricas, retrocedamos y observemos el caso en el que l = 0 . En este caso, no hay dependencia angular y la amplitud depende solo de la distancia radial, es decir, Ψ ( r , t ) → u ( r , t ) . En este caso, la ecuación de onda se reduce a

Esta ecuación se puede reescribir como

donde la cantidad ru satisface la ecuación de onda unidimensional. Por tanto, existen soluciones en la forma

donde F y G son soluciones generales a la ecuación de onda unidimensional y pueden interpretarse como una onda esférica entrante o saliente, respectivamente. Dichas ondas son generadas por una fuente puntual y hacen posible señales nítidas cuya forma se altera solo por una disminución en la amplitud a medida que aumenta r (ver una ilustración de una onda esférica en la parte superior derecha). Tales ondas existen solo en casos de espacio con dimensiones impares.

Para ver ejemplos físicos de soluciones de ondas no esféricas para la ecuación de ondas 3D que poseen dependencia angular, consulte radiación dipolo .

Onda esférica monocromática

Corte de frentes de onda esféricos, con una longitud de onda de 10 unidades, que se propagan desde una fuente puntual.

Aunque la palabra "monocromática" no es exactamente exacta ya que se refiere a la luz o radiación electromagnética con una frecuencia bien definida, el espíritu es descubrir el modo propio de la ecuación de onda en tres dimensiones. Siguiendo la derivación en la sección anterior sobre modos propios de onda plana , si nuevamente restringimos nuestras soluciones a ondas esféricas que oscilan en el tiempo con una frecuencia angular constante bien definida ω , entonces la función transformada ru ( r , t ) tiene simplemente soluciones de onda plana,

o

De esto podemos observar que la intensidad máxima de la oscilación de la onda esférica, caracterizada como la amplitud de la onda cuadrada

cae a una tasa proporcional a 1 / r 2 , un ejemplo de la ley del cuadrado inverso .

Solución de un problema de valor inicial general

La ecuación de onda es lineal en u y no se ve alterada por las traslaciones en el espacio y el tiempo. Por tanto, podemos generar una gran variedad de soluciones traduciendo y sumando ondas esféricas. Sea φ ( ξ , η , ζ ) una función arbitraria de tres variables independientes, y sea la forma de onda esférica F una función delta: es decir, sea F un límite débil de funciones continuas cuya integral es la unidad, pero cuyo soporte (la región donde la función no es cero) se reduce al origen. Sea una familia de ondas esféricas con centro en ( ξ , η , ζ ) , y sea r la distancia radial desde ese punto. Por lo tanto

Si u es una superposición de tales ondas con función de ponderación φ , entonces

el denominador 4 πc es una conveniencia.

De la definición de la función delta, u también se puede escribir como

donde α , β , y γ son coordenadas en la unidad de esfera S , y ω es el elemento de área en S . Este resultado tiene la interpretación de que u ( t , x ) es t multiplicado por el valor medio de φ en una esfera de radio ct centrada en x :

Resulta que

El valor medio es una función par de t , y por tanto si

luego

Estas fórmulas proporcionan la solución al problema del valor inicial de la ecuación de onda. Ellos muestran que la solución en un punto dado P , dado ( t , x , y , z ) sólo depende de los datos de la esfera de radio ct que está intersectado por las cono de luz hacia atrás extraídas de P . No , no depende de los datos en el interior de esta esfera. Así, el interior de la esfera es una laguna para la solución. Este fenómeno se llama principio de Huygens . Es cierto para números impares de dimensión espacial, donde para una dimensión la integración se realiza sobre el límite de un intervalo con respecto a la medida de Dirac. No se satisface en las dimensiones espaciales uniformes. El fenómeno de las lagunas se ha investigado ampliamente en Atiyah , Bott y Gårding (1970, 1973).

Ecuación de onda escalar en dos dimensiones espaciales

En dos dimensiones espaciales, la ecuación de onda es

Podemos usar la teoría tridimensional para resolver este problema si consideramos u como una función en tres dimensiones que es independiente de la tercera dimensión. Si

entonces la fórmula de solución tridimensional se convierte en

donde α y β son las dos primeras coordenadas en la esfera unitaria, y d ω es el elemento de área en la esfera. Esta integral se puede reescribir como una integral doble sobre el disco D con centro ( x , y ) y radio ct :

Es evidente que la solución en ( t , x , y ) depende no solo de los datos del cono de luz donde

pero también en datos que son internos a ese cono.

Ecuación de onda escalar en dimensión general y fórmulas de Kirchhoff

Queremos encontrar soluciones para u tt - Δ u = 0 para u  : R n × (0, ∞) → R con u ( x , 0) = g ( x ) y u t ( x , 0) = h ( x ) . Consulte a Evans para obtener más detalles.

Dimensiones impares

Suponga que n ≥ 3 es un número entero impar y gC m +1 ( R n ) , hC m ( R n ) para m = ( n + 1) / 2 . Sea γ n = 1 × 3 × 5 × ⋯ × ( n - 2) y sea

luego

uC 2 ( R n × [0, ∞))
u tt - Δ u = 0 en R n × (0, ∞)

Incluso dimensiones

Suponga que n ≥ 2 es un número entero par y gC m +1 ( R n ) , hC m ( R n ) , para m = ( n + 2) / 2 . Sea γ n = 2 × 4 × ⋯ × ny sea

luego

uC 2 ( R n × [0, ∞))
u tt - Δ u = 0 en R n × (0, ∞)

Problemas con los límites

Una dimensión espacial

La formulación de Sturm-Liouville

Una cuerda flexible que se estira entre dos puntos x = 0 y x = L satisface la ecuación de onda para t > 0 y 0 < x < L . En los puntos limítrofes, u puede satisfacer una variedad de condiciones limítrofes. Un formulario general que es apropiado para las aplicaciones es

donde un y b son no negativo. El caso en el que se requiere que u desaparezca en un punto final es el límite de esta condición cuando el respectivo a o b se acerca al infinito. El método de separación de variables consiste en buscar soluciones a este problema en la forma especial

Una consecuencia es que

El valor propio λ debe determinarse de modo que no exista una solución trivial del problema del valor en la frontera

Este es un caso especial del problema general de la teoría de Sturm-Liouville . Si un y b son positivos, los valores propios son todos positivos, y las soluciones son funciones trigonométricas. Se puede obtener una solución que satisfaga las condiciones iniciales integrables al cuadrado para u y u t a partir de la expansión de estas funciones en la serie trigonométrica apropiada.

Investigación por métodos numéricos

Aproximando la cuerda continua con un número finito de puntos de masa equidistantes, se obtiene el siguiente modelo físico:

Figura 1: Tres puntos de masa consecutivos del modelo discreto para una cuerda

Si cada punto de masa tiene la masa m , la tensión de la cuerda es f , la separación entre los puntos de masa es Δ x y u i , i = 1,…, n son el desplazamiento de estos n puntos de sus puntos de equilibrio (es decir, su posición en una línea recta entre los dos puntos de unión de la cuerda) la componente vertical de la fuerza hacia el punto i + 1 es

 

 

 

 

( 1 )

y la componente vertical de la fuerza hacia el punto i - 1 es

 

 

 

 

( 2 )

Tomando la suma de estas dos fuerzas y dividiendo por la masa m se obtiene el movimiento vertical:

 

 

 

 

( 3 )

Como la densidad de masa es

esto se puede escribir

 

 

 

 

( 4 )

La ecuación de onda se obtiene dejando Δ x → 0, en cuyo caso u i ( t ) toma la forma u ( x , t ) donde u ( x , t ) es una función continua de dos variables,··tu yotoma la forma 2 u / ∂ t 2 y

Pero la formulación discreta ( 3 ) de la ecuación de estado con un número finito de puntos de masa es la adecuada para una propagación numérica del movimiento de la cuerda. La condición de frontera

donde L es la longitud de la cuerda toma en la formulación discreta la forma que para los puntos más externos u 1 y u n las ecuaciones de movimiento son

 

 

 

 

( 5 )

y

 

 

 

 

( 6 )

mientras que para 1 < i < n

 

 

 

 

( 7 )

donde c = f / ρ .

Si la cuerda se aproxima con 100 puntos de masa discretos, se obtienen las 100 ecuaciones diferenciales de segundo orden acopladas ( 5 ), ( 6 ) y ( 7 ) o, de manera equivalente, 200 ecuaciones diferenciales de primer orden acopladas.

Propagando estos hasta los tiempos

Utilizando un método multipaso de octavo orden , se encuentran los 6 estados que se muestran en la figura 2:

Figura 2: La cuerda en 6 épocas consecutivas, la primera (roja) correspondiente al tiempo inicial con la cuerda en reposo
Figura 3: La cuerda en 6 épocas consecutivas
Figura 4: La cuerda en 6 épocas consecutivas
Figura 5: La cuerda en 6 épocas consecutivas
Figura 6: La cuerda en 6 épocas consecutivas
Figura 7: La cuerda en 6 épocas consecutivas

La curva roja es el estado inicial en el tiempo cero en el que la cuerda se "suelta" en una forma predefinida con todos . La curva azul es el estado en el tiempo, es decir, después de un tiempo que corresponde al tiempo que necesitaría una onda que se mueve con la velocidad de onda nominal c = f / ρ para un cuarto de la longitud de la cuerda.

La Figura 3 muestra la forma de la cuerda en esos momentos . La onda viaja en la dirección derecha con la velocidad c = f / ρ sin estar activamente restringida por las condiciones de contorno en los dos extremos de la cuerda. La forma de la onda es constante, es decir, la curva tiene la forma f ( x - ct ) .

La Figura 4 muestra la forma de la cuerda en esos momentos . La restricción en el extremo derecho comienza a interferir con el movimiento evitando que la ola eleve el extremo de la cuerda.

La Figura 5 muestra la forma de la cuerda en los momentos en que se invierte la dirección del movimiento. Las curvas roja, verde y azul son los estados en los momentos, mientras que las 3 curvas negras corresponden a los estados en los momentos en que la onda comienza a moverse hacia la izquierda.

La figura 6 y la figura 7 finalmente muestran la forma de la cuerda en los tiempos y . La ola ahora viaja hacia la izquierda y las restricciones en los puntos finales ya no están activas. Cuando finalmente se encuentre en el otro extremo de la cuerda, la dirección se invertirá nuevamente de una manera similar a lo que se muestra en la figura 6.

Varias dimensiones de espacio

Una solución de la ecuación de onda en dos dimensiones con una condición de límite de desplazamiento cero a lo largo de todo el borde exterior.

La teoría del valor de frontera inicial unidimensional puede extenderse a un número arbitrario de dimensiones espaciales. Considere un dominio D en m -dimensional x espacio, con frontera B . Entonces, la ecuación de onda debe satisfacerse si x está en D y t > 0 . En el límite de D , la solución u satisfará

donde n es la unidad normal exterior a B , y una es una función no negativa definida en B . El caso en el que u desaparece en B es un caso límite para un infinito que se aproxima. Las condiciones iniciales son

donde f y g se definen en D . Este problema puede ser resuelto mediante la ampliación de f y g en las funciones propias de la Laplaciano en D , que satisfacen las condiciones de contorno. Por tanto, la función propia v satisface

en D , y

en B .

En el caso de dos dimensiones espaciales, las funciones propias pueden interpretarse como los modos de vibración de una piel de tambor estirada sobre el límite B . Si B es un círculo, entonces estas funciones propias tienen un componente angular que es una función trigonométrica del ángulo polar θ , multiplicado por una función de Bessel (de orden entero) del componente radial. Más detalles se encuentran en la ecuación de Helmholtz .

Si el límite es una esfera en tres dimensiones espaciales, los componentes angulares de las funciones propias son armónicos esféricos y los componentes radiales son funciones de Bessel de orden medio entero.

Ecuación de onda no homogénea en una dimensión

La ecuación de onda no homogénea en una dimensión es la siguiente:

con condiciones iniciales dadas por

La función s ( x , t ) a menudo se denomina función fuente porque en la práctica describe los efectos de las fuentes de ondas en el medio que las transporta. Los ejemplos físicos de funciones de fuente incluyen la fuerza que impulsa una onda en una cuerda, o la carga o la densidad de corriente en el medidor de electromagnetismo de Lorenz .

Un método para resolver el problema del valor inicial (con los valores iniciales como se planteó anteriormente) es aprovechar una propiedad especial de la ecuación de onda en un número impar de dimensiones espaciales, a saber, que sus soluciones respetan la causalidad. Es decir, para cualquier punto ( x i , t i ) , el valor de u ( x i , t i ) sólo depende de los valores de f ( x i + ct i ) y f ( x i - ct i ) y la valores de la función g ( x ) entre ( x i - ct i ) y ( x i + ct i ) . Esto se puede ver en la fórmula de d'Alembert , enunciada anteriormente, donde estas cantidades son las únicas que aparecen en ella. Físicamente, si la velocidad de propagación máxima es c , entonces ninguna parte de la onda que no pueda propagarse a un punto dado en un tiempo dado puede afectar la amplitud en el mismo punto y tiempo.

En términos de encontrar una solución, esta propiedad de causalidad significa que para cualquier punto dado de la línea que se está considerando, la única área que debe considerarse es el área que abarca todos los puntos que podrían afectar causalmente al punto que se está considerando. Denotan el área que casualmente afecta a punto ( x i , t i ) como R C . Suponga que integramos la ecuación de onda no homogénea en esta región.

Para simplificar esto en gran medida, podemos usar el teorema de Green para simplificar el lado izquierdo y obtener lo siguiente:

El lado izquierdo es ahora la suma de tres integrales de línea a lo largo de los límites de la región de causalidad. Estos resultan ser bastante fáciles de calcular

En lo anterior, el término a integrar con respecto al tiempo desaparece porque el intervalo de tiempo involucrado es cero, por lo que d t = 0 .

Para los otros dos lados de la región, vale la pena señalar que x ± ct es una constante, es decir, x i ± ct i , donde el signo se elige apropiadamente. Usando esto, podemos obtener la relación d x ± c  d t = 0 , nuevamente eligiendo el signo correcto:

Y de manera similar para el segmento de límite final:

Sumando los tres resultados y volviéndolos a poner en la integral original:

Resolviendo para u ( x i , t i ) llegamos a

En la última ecuación de la secuencia, los límites de la integral sobre la función fuente se han hecho explícitos. Al observar esta solución, que es válida para todas las opciones ( x i , t i ) compatibles con la ecuación de onda, está claro que los dos primeros términos son simplemente la fórmula de d'Alembert, como se indicó anteriormente como la solución de la ecuación de onda homogénea en una dimensión. La diferencia está en el tercer término, la integral sobre la fuente.

Otros sistemas de coordenadas

En tres dimensiones, la ecuación de onda, cuando se escribe en coordenadas cilíndricas elípticas , puede resolverse mediante la separación de variables, lo que lleva a la ecuación diferencial de Mathieu .

Más generalizaciones

Ondas elásticas

La ecuación de onda elástica (también conocida como ecuación de Navier-Cauchy ) en tres dimensiones describe la propagación de ondas en un medio elástico homogéneo isótropo . La mayoría de los materiales sólidos son elásticos, por lo que esta ecuación describe fenómenos como las ondas sísmicas en la Tierra y las ondas ultrasónicas que se utilizan para detectar defectos en los materiales. Si bien es lineal, esta ecuación tiene una forma más compleja que las ecuaciones dadas anteriormente, ya que debe tener en cuenta tanto el movimiento longitudinal como el transversal:

dónde:

  • λ y μ son los llamados parámetros de Lamé que describen las propiedades elásticas del medio,
  • ρ es la densidad,
  • f es la función de la fuente (fuerza motriz),
  • y u es el vector de desplazamiento.

Al usar ∇ × (∇ × u ) = ∇ (∇ ⋅ u ) - ∇ ⋅ ∇ u = ∇ (∇ ⋅ u ) - ∆ u, la ecuación de onda elástica se puede reescribir en la forma más común de la ecuación de Navier-Cauchy.

Tenga en cuenta que en la ecuación de onda elástica, tanto la fuerza como el desplazamiento son cantidades vectoriales . Por lo tanto, esta ecuación a veces se conoce como ecuación de onda vectorial. Como ayuda para la comprensión, el lector observará que si f y ∇ ⋅ u se ponen a cero, esto se convierte (efectivamente) en la ecuación de Maxwell para la propagación del campo eléctrico E , que solo tiene ondas transversales.

Relación de dispersión

En los fenómenos de ondas dispersivas , la velocidad de propagación de las ondas varía con la longitud de onda de la onda, lo que se refleja en una relación de dispersión.

donde ω es la frecuencia angular y k es el vector de onda que describe las soluciones de onda plana . Para ondas de luz, la relación de dispersión es ω = ± c | k | , pero en general, la velocidad constante c se reemplaza por una velocidad de fase variable :

Ver también

Notas

Referencias

enlaces externos