La ecuación de onda es una ecuación diferencial parcial lineal de segundo orden importante 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 ecuación de onda de segundo orden describe la superposición de una onda entrante y saliente (es decir, más bien un campo de onda estacionaria), también se denomina "ecuación de onda bidireccional" (en contraste, el primer ordenLa ecuación de onda unidireccional 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 . [1] [2] [3] [4] [5] En 1746, d'Alembert descubrió la ecuación de onda unidimensional, y en diez años Euler descubrió la ecuación de onda tridimensional. [6]
Introducción
La ecuación de onda (bidireccional) es una ecuación diferencial parcial de segundo orden que puede restringir alguna función escalar u = u ( x 1 , x 2 ,…, x n ; t ) de una variable de tiempo ty una o más variables espaciales x 1 , x 2 ,…, x n . La cantidad u puede ser, por ejemplo, la presión en un líquido o gas, o el desplazamiento , a lo largo de alguna dirección específica, de las partículas de un sólido en vibración alejándose de sus posiciones de reposo. La ecuación 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) :
Una notación aún más compacta que se usa a veces en física es
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
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 . [7]
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 :
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.
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: [8]
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 . [9]
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 . [10] 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 puede derivarse directamente de la impedancia. [11]
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 a 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 ω . [12] Las soluciones de onda de chirp 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 se pueden recibir 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
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 usan 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 [13]
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. Estas 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. [ cita requerida ]
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
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 permanece inalterada por 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 g ∈ C m +1 ( R n ) , h ∈ C m ( R n ) para m = ( n + 1) / 2 . Sea γ n = 1 × 3 × 5 × ⋯ × ( n - 2) y sea
luego
- u ∈ C 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 g ∈ C m +1 ( R n ) , h ∈ C m ( R n ) , para m = ( n + 2) / 2 . Sea γ n = 2 × 4 × ⋯ × ny sea
luego
- u ∈ C 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 donde u se requiere para desaparecer en un punto final es el límite de esta condición cuando el respectivo una 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:
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,toma 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:
La curva roja es el estado inicial en el tiempo cero en el que la cuerda se "suelta" en una forma predefinida [14] con todos. La curva azul es el estado en el momentoes decir, después de un tiempo que corresponde al tiempo que necesitaría una onda que se mueve con la velocidad nominal de onda c = √ f / ρ para un cuarto de la longitud de la cuerda.
La figura 3 muestra la forma de la cuerda en los 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 los momentos . La restricción en el extremo derecho comienza a interferir con el movimiento evitando que la ola suba el extremo de la cuerda.
La figura 5 muestra la forma de la cuerda en los momentos cuando se invierte la dirección del movimiento. Las curvas roja, verde y azul son los estados en el momento mientras que las 3 curvas negras corresponden a los estados a veces con la ola comenzando a moverse hacia la izquierda.
La figura 6 y la figura 7 finalmente muestran la forma de la cuerda en los momentos 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
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 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 conduce 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.
Usando ∇ × (∇ × 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
- Atenuación acústica
- Ecuación de ondas acústicas
- Transformación de Bateman
- Ecuación de ondas electromagnéticas
- Ecuación de Helmholtz
- Ecuación de onda electromagnética no homogénea
- Operador de Laplace
- Matemáticas de oscilación
- Ecuaciones de Maxwell
- Ecuación de onda unidireccional
- Ecuación de Schrödinger
- Onda estacionaria
- Vibraciones de una membrana circular
- Teoría del absorbedor de Wheeler-Feynman
Notas
- ^ Cannon, John T .; Dostrovsky, Sigalia (1981). La evolución de la dinámica, teoría de la vibración desde 1687 hasta 1742 . Estudios de Historia de las Matemáticas y Ciencias Físicas. 6 . Nueva York: Springer-Verlag. págs. ix + 184 págs. ISBN 978-0-3879-0626-3.
- ^ GREY, JW (julio de 1983). "RESEÑAS DE LIBROS". Boletín de la American Mathematical Society . Series nuevas. 9 (1). (consultado el 13 de noviembre de 2012).
- ^ Gerard F Wheeler. The Vibrating String Controversy , (consultado el 13 de noviembre de 2012). Soy. J. Phys., 1987, v55, n1, p33-37.
- ^ Para obtener una colección especial de los 9 artículos innovadores de los tres autores, consulte Primera aparición de la ecuación de onda: D'Alembert, Leonhard Euler, Daniel Bernoulli. - la controversia sobre las cuerdas vibrantes (consultado el 13 de noviembre de 2012). Herman HJ Lynge e hijo.
- ↑ Para conocer las contribuciones de De Lagrange a la ecuación de ondas acústicas, se pueden consultar Acústica: Introducción a sus principios físicos y aplicaciones Allan D. Pierce, Soc acústica de América, 1989; página 18 (consultado el 9 de diciembre de 2012).
- ^ a b c Speiser, David. Descubriendo los principios de la mecánica 1600–1800 , pág. 191 (Basilea: Birkhäuser, 2008).
- ^ Tipler, Paul y Mosca, Gene. Física para científicos e ingenieros, Volumen 1: Mecánica, oscilaciones y ondas; Termodinámica , págs. 470–471 (Macmillan, 2004).
- ^ Eric W. Weisstein . "Solución de d'Alembert" . MathWorld . Consultado el 21 de enero de 2009 .
- ↑ D'Alembert (1747) "Recherches sur la courbe que forme une corde tenduë mise en vibration" (Investigaciones sobre la curva que forma una cuerda tensa [cuando] se pone en vibración), Histoire de l'académie royale des sciences et belles lettres de Berlín , vol. 3, páginas 214–219.
- Véase también: D'Alembert (1747) "Suite des recherches sur la courbe que forme une corde tenduë mise en vibration" (Más investigaciones sobre la curva que forma una cuerda tensa [cuando] se pone en vibración), Histoire de l'académie royale des sciences et belles lettres de Berlin , vol. 3, páginas 220–249.
- Véase también: D'Alembert (1750) "Addition au mémoire sur la courbe que forme une corde tenduë mise en vibration", Histoire de l'académie royale des sciences et belles lettres de Berlin , vol. 6, páginas 355–360.
- ^ http://math.arizona.edu/~kglasner/math456/linearwave.pdf .
- ^ Bschorr, Oskar; Raida, Hans-Joachim (marzo de 2020). "Ecuación de onda unidireccional derivada del teorema de impedancia" . Acústica . 2 (1): 164-170. doi : 10.3390 / acoustics2010012 .
- ^ V Guruprasad (2015), "Evidencia de observación para los modos de onda viajera con cambios proporcionales de distancia de rodamiento", EPL , 110 (5): 54001, arXiv : 1507.08222 , Bibcode : 2015EL .... 11054001G , doi : 10.1209 / 0295-5075 / 110 / 54001 , S2CID 42285652
- ^ Jackson, John David. Electrodinámica clásica (3ª ed.). Wiley. pag. 425. ISBN 978-0-471-30932-1.
- ^ El estado inicial de "Investigación por métodos numéricos" se establece con splines cuadráticos de lasiguiente manera:
- por
- por
- por
Referencias
- MF Atiyah, R. Bott, L. Garding, " Lagunas para operadores diferenciales hiperbólicos con coeficientes constantes I ", Acta Math. , 124 (1,970), 109-189.
- MF Atiyah, R. Bott y L. Garding, " Lagunas para operadores diferenciales hiperbólicos con coeficientes II constantes ", Acta Math. , 131 (1973), 145-206.
- R. Courant, D. Hilbert, Métodos de física matemática, vol II . Interscience (Wiley) Nueva York, 1962.
- L. Evans, "Ecuaciones diferenciales parciales". Providencia de la Sociedad Americana de Matemáticas, 1998.
- " Ecuaciones de onda lineal ", EqWorld: El mundo de las ecuaciones matemáticas.
- " Ecuaciones de onda no lineales ", EqWorld: El mundo de las ecuaciones matemáticas.
- William C. Lane, " MISN-0-201 La ecuación de onda y sus soluciones ", Proyecto PHYSNET .
enlaces externos
- Ecuaciones de onda no lineales de Stephen Wolfram y Rob Knapp, Explorador de ecuaciones de onda no lineal de Wolfram Demonstrations Project .
- Los aspectos matemáticos de las ecuaciones de onda se discuten en Dispersive PDE Wiki .
- Graham W Griffiths y William E. Schiesser (2009). Ondas lineales y no lineales . Scholarpedia , 4 (7): 4308. doi: 10.4249 / scholarpedia.4308