Aproximación de Boussinesq (ondas de agua)

Una aproximación válida para ondas débilmente no lineales y bastante largas.

Este artículo trata sobre la aproximación de Boussinesq para ondas en una superficie de fluido en movimiento libre. Para otros usos, consulte Aproximación de Boussinesq .

Simulación de olas periódicas sobre un banco submarino con un modelo tipo Boussinesq. Las olas se propagan sobre un banco submarino de forma elíptica en una playa plana. Este ejemplo combina varios efectos de las olas y el agua poco profunda , incluida la refracción , la difracción , la falta de linealidad y la falta de linealidad débil .

En dinámica de fluidos , la aproximación de Boussinesq para ondas de agua es una aproximación válida para ondas débilmente no lineales y bastante largas . La aproximación lleva el nombre de Joseph Boussinesq , quien las derivó por primera vez en respuesta a la observación de John Scott Russell de la onda de traducción (también conocida como onda solitaria o solitón ). El artículo de 1872 de Boussinesq presenta las ecuaciones ahora conocidas como ecuaciones de Boussinesq . ^[1]

La aproximación de Boussinesq para ondas de agua tiene en cuenta la estructura vertical de la velocidad de flujo horizontal y vertical . Esto da como resultado ecuaciones diferenciales parciales no lineales , llamadas ecuaciones de tipo Boussinesq , que incorporan dispersión de frecuencia (como opuesto a las ecuaciones de aguas poco profundas , que no son de dispersión de frecuencia). En la ingeniería costera , las ecuaciones de tipo Boussinesq se utilizan con frecuencia en modelos informáticos para la simulación de olas de agua en mares y puertos poco profundos .

Si bien la aproximación de Boussinesq es aplicable a ondas bastante largas, es decir, cuando la longitud de onda es grande en comparación con la profundidad del agua, la expansión de Stokes es más apropiada para ondas cortas (cuando la longitud de onda es del mismo orden que la profundidad del agua, o más corta ).

Aproximación de Boussinesq

Ondas periódicas en la aproximación de Boussinesq, mostradas en una sección transversal vertical en la dirección de propagación de la onda . Observe los valles planos y las crestas afiladas , debido a la no linealidad de la onda. Este caso (dibujado a escala ) muestra una ola con una longitud de onda igual a 39,1  m , la altura de la ola es de 1,8 m ( es decir, la diferencia entre la elevación de la cresta y la vaguada) y la profundidad media del agua es de 5 m, mientras que la aceleración gravitacional es de 9,81 m / s ² .

La idea esencial en la aproximación de Boussinesq es la eliminación de la coordenada vertical de las ecuaciones de flujo, conservando algunas de las influencias de la estructura vertical del flujo bajo las olas del agua . Esto es útil porque las ondas se propagan en el plano horizontal y tienen un comportamiento diferente (no ondulatorio) en la dirección vertical. A menudo, como en el caso de Boussinesq, el interés se centra principalmente en la propagación de ondas.

Esta eliminación de la coordenada vertical fue realizada por primera vez por Joseph Boussinesq en 1871, para construir una solución aproximada para la onda solitaria (u onda de traslación ). Posteriormente, en 1872, Boussinesq derivó las ecuaciones conocidas hoy en día como ecuaciones de Boussinesq.

Los pasos en la aproximación de Boussinesq son:

• una expansión de Taylor está hecha de la velocidad del flujo horizontal y vertical (o potencial de velocidad ) alrededor de una cierta elevación ,
• esta expansión de Taylor se trunca a un número finito de términos,
• se utilizan la conservación de masa (ver ecuación de continuidad ) para un flujo incompresible y la condición de rizo cero para un flujo irrotacional , para reemplazar las derivadas parciales verticales de cantidades en la expansión de Taylor con derivadas parciales horizontales .

Posteriormente, se aplica la aproximación de Boussinesq al resto de las ecuaciones de flujo, con el fin de eliminar la dependencia de la coordenada vertical. Como resultado, las ecuaciones diferenciales parciales resultantes están en términos de funciones de las coordenadas horizontales (y tiempo ).

Como un ejemplo, considere flujo potencial sobre un lecho horizontal en el ( x, z avión), con x la horizontal y z de la vertical de coordenadas . El lecho está ubicado en z = - h , donde h es la profundidad media del agua. Se hace una expansión de Taylor del potencial de velocidad φ (x, z, t) alrededor del nivel del lecho z = - h : ^[2]

${\displaystyle {\begin{aligned}\varphi \,=\,&\varphi _{b}\,+\,(z+h)\,\left[{\frac {\partial \varphi }{\partial z}}\right]_{z=-h}\,+\,{\frac {1}{2}}\,(z+h)^{2}\,\left[{\frac {\partial ^{2}\varphi }{\partial z^{2}}}\right]_{z=-h}\,\\&+\,{\frac {1}{6}}\,(z+h)^{3}\,\left[{\frac {\partial ^{3}\varphi }{\partial z^{3}}}\right]_{z=-h}\,+\,{\frac {1}{24}}\,(z+h)^{4}\,\left[{\frac {\partial ^{4}\varphi }{\partial z^{4}}}\right]_{z=-h}\,+\,\cdots ,\end{aligned}}}$

donde φ _b (x, t) es el potencial de velocidad en el lecho. Invocando la ecuación de Laplace para φ , como válida para flujo incompresible , da:

${\displaystyle {\begin{aligned}\varphi \,=\,&\left\{\,\varphi _{b}\,-\,{\frac {1}{2}}\,(z+h)^{2}\,{\frac {\partial ^{2}\varphi _{b}}{\partial x^{2}}}\,+\,{\frac {1}{24}}\,(z+h)^{4}\,{\frac {\partial ^{4}\varphi _{b}}{\partial x^{4}}}\,+\,\cdots \,\right\}\,\\&+\,\left\{\,(z+h)\,\left[{\frac {\partial \varphi }{\partial z}}\right]_{z=-h}\,-\,{\frac {1}{6}}\,(z+h)^{3}\,{\frac {\partial ^{2}}{\partial x^{2}}}\left[{\frac {\partial \varphi }{\partial z}}\right]_{z=-h}\,+\,\cdots \,\right\}\\=\,&\left\{\,\varphi _{b}\,-\,{\frac {1}{2}}\,(z+h)^{2}\,{\frac {\partial ^{2}\varphi _{b}}{\partial x^{2}}}\,+\,{\frac {1}{24}}\,(z+h)^{4}\,{\frac {\partial ^{4}\varphi _{b}}{\partial x^{4}}}\,+\,\cdots \,\right\},\end{aligned}}}$

dado que la velocidad vertical ∂ φ / ∂ z es cero en el lecho horizontal - impermeable - z = - h . Esta serie se puede truncar posteriormente a un número finito de términos.

Ecuaciones originales de Boussinesq

Derivación

Para ondas de agua en un fluido incompresible y flujo de irritación en el plano ( x , z ), las condiciones de contorno en la elevación de la superficie libre z = η ( x , t ) son: ^[3]

${\displaystyle {\begin{aligned}{\frac {\partial \eta }{\partial t}}\,&+\,u\,{\frac {\partial \eta }{\partial x}}\,-\,w\,=\,0\\{\frac {\partial \varphi }{\partial t}}\,&+\,{\frac {1}{2}}\,\left(u^{2}+w^{2}\right)\,+\,g\,\eta \,=\,0,\end{aligned}}}$

dónde:

u es el componente de velocidad de flujo horizontal : u = ∂ φ / ∂ x ,

w es el componente de la velocidad de flujo vertical : w = ∂ φ / ∂ z ,

g es la aceleración por gravedad .

Ahora, la aproximación de Boussinesq para el potencial de velocidad φ , como se indicó anteriormente, se aplica en estas condiciones de contorno . Además, en las ecuaciones resultantes sólo se retienen los términos lineal y cuadrático con respecto a η y u _b (con u _b = ∂ φ _b / ∂ x la velocidad horizontal en el lecho z = - h ). Se supone que los términos cúbicos y de orden superior son insignificantes. Entonces, se obtienen las siguientes ecuaciones diferenciales parciales :

conjunto A - Boussinesq (1872), ecuación (25)

${\displaystyle {\begin{aligned}{\frac {\partial \eta }{\partial t}}\,&+\,{\frac {\partial }{\partial x}}\,\left[\left(h+\eta \right)\,u_{b}\right]\,=\,{\frac {1}{6}}\,h^{3}\,{\frac {\partial ^{3}u_{b}}{\partial x^{3}}},\\{\frac {\partial u_{b}}{\partial t}}\,&+\,u_{b}\,{\frac {\partial u_{b}}{\partial x}}\,+\,g\,{\frac {\partial \eta }{\partial x}}\,=\,{\frac {1}{2}}\,h^{2}\,{\frac {\partial ^{3}u_{b}}{\partial t\,\partial x^{2}}}.\end{aligned}}}$

Este conjunto de ecuaciones se ha derivado para un lecho horizontal plano, es decir , la profundidad media h es una constante independiente de la posición x . Cuando los lados derechos de las ecuaciones anteriores se establecen en cero, se reducen a las ecuaciones de aguas poco profundas .

Con algunas aproximaciones adicionales, pero en el mismo orden de precisión, el conjunto A anterior se puede reducir a una única ecuación diferencial parcial para la elevación de la superficie libre η :

conjunto B - Boussinesq (1872), ecuación (26)

${\displaystyle {\frac {\partial ^{2}\eta }{\partial t^{2}}}\,-\,gh\,{\frac {\partial ^{2}\eta }{\partial x^{2}}}\,-\,gh\,{\frac {\partial ^{2}}{\partial x^{2}}}\left({\frac {3}{2}}\,{\frac {\eta ^{2}}{h}}\,+\,{\frac {1}{3}}\,h^{2}\,{\frac {\partial ^{2}\eta }{\partial x^{2}}}\right)\,=\,0.}$

A partir de los términos entre paréntesis, la importancia de la no linealidad de la ecuación se puede expresar en términos del número de Ursell . En cantidades adimensionales , utilizando la profundidad del agua hy la aceleración gravitacional g para la no dimensionalización, esta ecuación dice, después de la normalización : ^[4]

${\displaystyle {\frac {\partial ^{2}\psi }{\partial \tau ^{2}}}\,-\,{\frac {\partial ^{2}\psi }{\partial \xi ^{2}}}\,-\,{\frac {\partial ^{2}}{\partial \xi ^{2}}}\left(\,3\,\psi ^{2}\,+\,{\frac {\partial ^{2}\psi }{\partial \xi ^{2}}}\,\right)\,=\,0,}$

con:

${\displaystyle \psi \,=\,{\frac {1}{2}}\,{\frac {\eta }{h}}}$	: la elevación de la superficie adimensional,
${\displaystyle \tau \,=\,{\sqrt {3}}\,t\,{\sqrt {\frac {g}{h}}}}$	: el tiempo adimensional, y
${\displaystyle \xi \,=\,{\sqrt {3}}\,{\frac {x}{h}}}$	: la posición horizontal adimensional.

Velocidad de fase lineal al cuadrado c ² / ( gh ) en función del número de onda relativo kh .
A = Boussinesq (1872), ecuación (25),
B = Boussinesq (1872), ecuación (26),
C = teoría de onda lineal completa, ver dispersión (ondas de agua)

Dispersión de frecuencia lineal

Las ondas de agua de diferentes longitudes de onda viajan con diferentes velocidades de fase , un fenómeno conocido como dispersión de frecuencia . Para el caso de amplitud de onda infinitesimal , la terminología es dispersión de frecuencia lineal . Las características de dispersión de frecuencia de una ecuación de tipo Boussinesq se pueden utilizar para determinar el rango de longitudes de onda, para lo cual es una aproximación válida .

Las características de dispersión de frecuencia lineal para el conjunto A de ecuaciones anterior son: ^[5]

${\displaystyle c^{2}\,=\;gh\,{\frac {1\,+\,{\frac {1}{6}}\,k^{2}h^{2}}{1\,+\,{\frac {1}{2}}\,k^{2}h^{2}}},}$

con:

• c la velocidad de fase ,
• k el número de onda ( k = 2π / λ , con λ la longitud de onda ).

El error relativo en la velocidad de fase c para el conjunto A , en comparación con la teoría lineal para ondas de agua , es menos del 4% para un número de onda relativo kh <½ π . Entonces, en aplicaciones de ingeniería , el conjunto A es válido para longitudes de onda λ mayores que 4 veces la profundidad del agua h .

Las características de dispersión de frecuencia lineal de la ecuación B son: ^[5]

${\displaystyle c^{2}\,=\,gh\,\left(1\,-\,{\frac {1}{3}}\,k^{2}h^{2}\right).}$

El error relativo en la velocidad de fase para la ecuación B es menos del 4% para kh <2π / 7 , equivalente a longitudes de onda λ mayores que 7 veces la profundidad del agua h , llamadas ondas bastante largas . ^[6]

Para ondas cortas con k ² h ² > 3, la ecuación B pierde sentido físicamente, porque ya no existen soluciones con valores reales de la velocidad de fase . El conjunto original de dos ecuaciones diferenciales parciales (Boussinesq, 1872, ecuación 25, ver el conjunto A anterior) no tiene este inconveniente.

Las ecuaciones de aguas poco profundas tienen un error relativo en la velocidad de fase inferior al 4% para longitudes de onda λ superiores a 13 veces la profundidad del agua h .

Ecuaciones y extensiones de tipo Boussinesq

Existe una abrumadora cantidad de modelos matemáticos a los que se hace referencia como ecuaciones de Boussinesq. Esto puede llevar fácilmente a confusión, ya que a menudo se las conoce como ecuaciones de Boussinesq, mientras que de hecho se considera una variante de las mismas. Por tanto, es más apropiado llamarlas ecuaciones de tipo Boussinesq . Estrictamente hablando, las ecuaciones de Boussinesq son el conjunto B mencionado anteriormente , ya que se utiliza en el análisis en el resto de su artículo de 1872.

Algunas direcciones en las que se han extendido las ecuaciones de Boussinesq son:

• batimetría variable ,
• dispersión de frecuencia mejorada ,
• comportamiento no lineal mejorado ,
• haciendo una expansión de Taylor alrededor de diferentes elevaciones verticales ,
• dividiendo el dominio de los fluidos en capas y aplicando la aproximación de Boussinesq en cada capa por separado,
• inclusión de rompimiento de olas ,
• inclusión de tensión superficial ,
• extensión a ondas internas en una interfaz entre dominios fluidos de diferente densidad de masa ,
• derivación de un principio variacional .

Aproximaciones adicionales para la propagación de ondas unidireccionales

Si bien las ecuaciones de Boussinesq permiten que las ondas viajen simultáneamente en direcciones opuestas, a menudo es ventajoso considerar solo las ondas que viajan en una dirección. Bajo pequeños supuestos adicionales, las ecuaciones de Boussinesq se reducen a:

• la ecuación de Korteweg-de Vries para la propagación de ondas en una dimensión horizontal ,
• la ecuación de Kadomtsev-Petviashvili para la propagación de ondas (casi unidireccional) en dos dimensiones horizontales ,
• la ecuación de Schrödinger no lineal ( ecuación NLS) para la amplitud de valores complejos de las ondas de banda estrecha ( ondas de modulación lenta ).

Además de las soluciones de ondas solitarias, la ecuación de Korteweg-de Vries también tiene soluciones periódicas y exactas, llamadas ondas cnoidales . Estas son soluciones aproximadas de la ecuación de Boussinesq.

Modelos numéricos

Una simulación con un modelo de olas tipo Boussinesq de olas cercanas a la costa que viajan hacia la entrada del puerto. La simulación es con el módulo BOUSS-2D de SMS .

Más rápido que la simulación en tiempo real con el módulo Boussinesq de Celeris, que muestra el rompimiento y la refracción de las olas cerca de la playa. El modelo proporciona un entorno interactivo.

Para la simulación del movimiento de las olas cerca de costas y puertos, existen modelos numéricos, tanto comerciales como académicos, que emplean ecuaciones de tipo Boussinesq. Algunos ejemplos comerciales son los módulos de ondas tipo Boussinesq en MIKE 21 y SMS . Algunos de los modelos gratuitos de Boussinesq son Celeris, ^[7] COULWAVE, ^[8] y FUNWAVE. ^[9] La mayoría de los modelos numéricos emplean diferencias finitas , de volúmenes finitos o elementos finitos técnicas para la discretización de las ecuaciones del modelo. Las revisiones científicas y las intercomparaciones de varias ecuaciones de tipo Boussinesq, su aproximación numérica y su rendimiento son, por ejemplo, Kirby (2003), Dingemans (1997 , Parte 2, Capítulo 5) y Hamm, Madsen & Peregrine (1993) .

Notas

1. Este artículo (Boussinesq, 1872) comienza con: "Tous les ingénieurs connaissent les belles expériences de J. Scott Russell et M. Basin sur la production et la propagation des ondes solitaires" ( "Todos los ingenieros conocen los hermosos experimentos de J. Scott Russell y M. Basin sobre la generación y propagación de ondas solitarias " ).
2. Dingemans (1997), p. 477.
3. Dingemans (1997), p. 475.
4. Johnson (1997), p. 219
5. Dingemans (1997), p. 521.
6. Dingemans (1997), p. 473 y 516.
7. "Celeria.org - Modelo de onda Celeris Boussinesq" . Celeria.org - Modelo de onda Celeris Boussinesq .
8. "ISEC - Modelos" . isec.nacse.org .
9. "James T. Kirby, programa Funwave" . www1.udel.edu .

Referencias

• Boussinesq, J. (1871). "Théorie de l'intumescence liquide, applelée onde solitaire ou de translation, se propageant dans un canal rectangulaire" . Comptes Rendus de l'Académie des Sciences . 72 : 755–759.
• Boussinesq, J. (1872). "Théorie des ondes et des remous qui se propagent le long d'un canal rectangulaire horizontal, en communiquant au liquide contenu dans ce canal des vitesses sensiblement pareilles de la surface au fond" . Journal de Mathématiques Pures et Appliquées . Deuxième Série. 17 : 55-108.
• Dingemans, MW (1997). Propagación de ondas sobre fondos irregulares . Serie avanzada sobre ingeniería oceánica 13 . World Scientific, Singapur. ISBN 978-981-02-0427-3. Archivado desde el original el 8 de febrero de 2012 . Consultado el 21 de enero de 2008 . Consulte la Parte 2, Capítulo 5 .
• Hamm, L .; Madsen, PA; Peregrine, DH (1993). "Transformación de olas en la zona cercana a la costa: una revisión". Ingeniería Costera . 21 (1-3): 5-39. doi : 10.1016 / 0378-3839 (93) 90044-9 .
• Johnson, RS (1997). Una introducción moderna a la teoría matemática de las ondas de agua . Textos de Cambridge en Matemática Aplicada. 19 . Prensa de la Universidad de Cambridge. ISBN 0-521-59832-X.
• Kirby, JT (2003). "Modelos y aplicaciones de Boussinesq a la propagación de olas cercanas a la costa, procesos de zonas de surf y corrientes inducidas por olas". En Lakhan, VC (ed.). Avances en la modelización costera . Serie de oceanografía de Elsevier. 67 . Elsevier. págs. 1-41. ISBN 0-444-51149-0.
• Peregrine, DH (1967). "Largas olas en una playa". Revista de Mecánica de Fluidos . 27 (4): 815–827. Código bibliográfico : 1967JFM .... 27..815P . doi : 10.1017 / S0022112067002605 .
• Peregrine, DH (1972). "Ecuaciones para ondas de agua y las aproximaciones detrás de ellas". En Meyer, RE (ed.). Olas en las playas y transporte de sedimentos resultante . Prensa académica. págs. 95-122. ISBN 0-12-493250-9.